<div dir="ltr"><div dir="ltr">Dear Lorenzo, Paolo, <br></div><div dir="ltr"><br></div><div>I have successfully sent PHonon code after uncomenting the line you suggested, thank you very much for your help (it will take some time to get the final result, though, but I am confident that it will work). And for your answer below, as well. The garnet structure is quite unique as it is at the same time very big and very symmetric. So I guess this type of configurations has little chance to happen in other circumstances. I guess I will compute this structure "the old way" and re-comment the line afterwards, to prevent pathological cases. <br></div><div><br></div><div>Thanks again, <br></div><div><br></div><div>Merlin <br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Le jeu. 7 mai 2020 à 10:14, Lorenzo Paulatto <<a href="mailto:paulatz@gmail.com">paulatz@gmail.com</a>> a écrit :<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">I forgot to answer to the list before.<br>
I confirme that using the "old" random matrix it "works" in the sense <br>
that this specific case works. It will fail in other cases though.<br>
<br>
cheers<br>
<br>
On 5/7/20 9:31 AM, Lorenzo Paulatto wrote:<br>
> Yes,<br>
> the phonon code checks if modes are degenerate by symmetrizing a random <br>
> matrix and checking if the eigenvalues are the same. For a sufficiently <br>
> large systems there are going to be accidental degeneracies.<br>
> <br>
> Can you try to uncomment<br>
> !!#define __UNIFORM_DISTRIB<br>
> in PHonon/PH/random_matrix.f90 and see if it fixes, or at least changes, <br>
> the error?<br>
> <br>
> cheers<br>
> <br>
> <br>
> <br>
> <br>
> On 5/7/20 9:25 AM, Paolo Giannozzi wrote:<br>
>> This phonon error shows up quite frequently, especially in large cells:<br>
>> <a href="https://www.mail-archive.com/search?a=1&l=users%40lists.quantum-espresso.org&haswords=set_irr_sym_new&x=14&y=14&from=&subject=&datewithin=1d&date=&notwords=&o=newest" rel="noreferrer" target="_blank">https://www.mail-archive.com/search?a=1&l=users%40lists.quantum-espresso.org&haswords=set_irr_sym_new&x=14&y=14&from=&subject=&datewithin=1d&date=&notwords=&o=newest</a> <br>
>><br>
>> While in most cases it is due to "quasi-symmetric" atomic positions <br>
>> and cells, in this case:<br>
>> <a href="https://www.mail-archive.com/users@lists.quantum-espresso.org/msg37662.html" rel="noreferrer" target="_blank">https://www.mail-archive.com/users@lists.quantum-espresso.org/msg37662.html</a> <br>
>><br>
>> the reasons is no so obvious. Any idea?<br>
>><br>
>> Paolo<br>
>> -- <br>
>> Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>
>> Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
>> Phone +39-0432-558216, fax +39-0432-558222<br>
>><br>
>><br>
>> _______________________________________________<br>
>> developers mailing list<br>
>> <a href="mailto:developers@lists.quantum-espresso.org" target="_blank">developers@lists.quantum-espresso.org</a><br>
>> <a href="https://lists.quantum-espresso.org/mailman/listinfo/developers" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/developers</a><br>
>><br>
> <br>
<br>
-- <br>
Lorenzo Paulatto - Paris<br>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr">Merlin Méheut<br>adresse labo:<br>GET - OMP   - Université Paul Sabatier <br>14 avenue Edouard Belin<br>31400 Toulouse<br>FRANCE <br>tel: (+33) 5 61 33 26 17<br><br></div></div></div>