<div dir="ltr"><div><div><div><div><div>dear Paolo,<br><br></div>Tsung-Lung Li explained the atoms will be moved idependently according to the forces acting on them.<br></div>And you said that is correct. It seems to me that somehow contradicts to say that this still preserves the symmetry.<br></div>Please clarify if I'm wrong about this. Because I understand right to keep the symmetry the atoms should be moved not only according to their forces, but also preserving their coordinates, which suggests that their should move isotropically, right?<br><br></div>Thank you very much,<br><br></div>Ricardo.<br></div><div class="gmail_extra"><br><div class="gmail_quote">On 5 October 2017 at 14:09, Paolo Giannozzi <span dir="ltr"><<a href="mailto:p.giannozzi@gmail.com" target="_blank">p.giannozzi@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><span class="">On Thu, Oct 5, 2017 at 4:03 PM, Tsung-Lung Li <span dir="ltr"><<a href="mailto:quantum@mail.ncyu.edu.tw" target="_blank">quantum@mail.ncyu.edu.tw</a>></span> wrote:<br><br></span><div class="gmail_extra"><div class="gmail_quote"><span class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
It seems to me QE does the following.<br>
<br>
1. Initially, QE uses space_group and the 7 inequivalent atoms in<br>
    ATOMIC_POSITIONS to produce 24 atoms in the crystal unit cell.<br></blockquote><div><br></div></span><div>correct<br> <br></div><span class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
2. In each subsequent geometric step of the 'relax' calculation,<br>
    the 24 atoms are moved independently according to the forces<br>
    acting on each atom.<br></blockquote><div><br></div></span><div>correct<br> <br></div><span class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Therefore, at the end of the 'relax' calculation, the space_group=15 symmetry is not always preserved.<br></blockquote><div><br></div></span><div>not correct. Structural optimization based on forces does not break the symmetry, with occasional unfortunate exceptions due to numerical noise<br> <br></div><span class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Can QE only move the 7 inequivalent atoms and use space_group symmetry to construct the other 17 atoms for next geometric step?<br></blockquote><div><br></div></span><div>It could be done in principle, but it would require some extensive changes to the structural optimization algorithm, with dubious advantages in return.<br></div><div> <br></div><div>Paolo<br><br></div><div><div class="h5"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
This way, the space_group symmetry will be preserved at the end of the simulation.<br>
<br>
Thank you in advance.<br>
<br>
Best Regards,<br>
Tsung-Lung Li<br>
<br>
------------------------------<wbr>------------------------------<wbr>------<br>
&control<br>
   calculation = 'relax'<br>
   restart_mode = 'from_scratch'<br>
/<br>
&system<br>
   space_group = 15<br>
   celldm(1) = 11.86181088389489 ! bohr<br>
   celldm(2) = 2.41851202803887<br>
   celldm(3) = 0.903616377250279<br>
   celldm(5) = -0.408489773780522<br>
   ntyp = 4<br>
   nat = 7<br>
/<br>
<br>
ATOMIC_SPECIES<br>
    Ca   40.0780000000  Ca_pbe_v1.uspp.F.UPF<br>
     S   32.0650000000  S_pbe_v1.2.uspp.F.UPF<br>
     O   15.9994000000  O.pbe-n-kjpaw_psl.0.1.UPF<br>
     H    1.0079000000  H.pbe-rrkjus_psl.0.1.UPF<br>
<br>
ATOMIC_POSITIONS crystal_sg<br>
   Ca  0.00000   0.17050   0.25000<br>
   S   0.00000   0.32727   0.75000<br>
   O   0.08319   0.27218   0.59103<br>
   O   0.19997   0.38195   0.91298<br>
   O  -0.20823   0.06826  -0.07831<br>
   H  -0.25800   0.08700  -0.23400<br>
   H  -0.24400   0.02000  -0.07700<br>
<br>
______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org" target="_blank">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listi<wbr>nfo/pw_forum</a></blockquote></div></div></div><span class="HOEnZb"><font color="#888888"><br clear="all"><br>-- <br><div class="m_-4362391197076070643gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>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></div></div></div></div></div>
</font></span></div></div>
<br>______________________________<wbr>_________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/<wbr>listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div><div>Ricardo Afonso<br></div>Student of Magnetism and Superconductivity Group<br></div><div>Federal University of Sao Carlos<br></div></div></div></div></div></div>
</div>