<p dir="ltr">If the quantity you are trying to compute is sensitive to differences of the order of 4*10^{-5} A in atomic positions, you have chosen an hopeless task, in my opinion.</p>
<p dir="ltr">Paolo</p>
<div class="gmail_extra"><br><div class="gmail_quote">Il 28/dic/2016 02:18 PM, "Pablo García Risueño" <<a href="mailto:garcia.risueno@gmail.com" target="_blank">garcia.risueno@gmail.com</a>> ha scritto:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Dear professor<br><br></div>Thank you very much for your reply. The differences are important in this case, my final calculated quantities are very sensitive to these optimized positions. Should I change any of the ***conv_thr variables, or other variable, to have the same result for same inputs?<br><br></div>Is there any random number in the algorithm of relax calculations, can we be sure that different outputs with the same input are not due to any bug?<br><br></div>Thank you very much. Best regards.<br></div><div class="gmail_extra"><br><div class="gmail_quote">2016-12-28 14:10 GMT+01:00 Paolo Giannozzi <span dir="ltr"><<a href="mailto:p.giannozzi@gmail.com" target="_blank">p.giannozzi@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><p dir="ltr">The differences you find are very small. Nothing to worry about in my opinion.</p>
<p dir="ltr">By the way: Davidson diagonalization is typically faster than CG; do not specify incompatible options in K_POINTS (gamma or automatic, not both; gamma should be used unless you have a good reason not to)</p><span class="m_6133018495245676239HOEnZb"><font color="#888888">
</font></span><p dir="ltr"><span class="m_6133018495245676239HOEnZb"><font color="#888888">Paolo</font></span></p><div><div class="m_6133018495245676239h5"><br>
Il 28/dic/2016 01:46 PM, "Pablo García Risueño" <<a href="mailto:garcia.risueno@gmail.com" target="_blank">garcia.risueno@gmail.com</a>> ha scritto:<br>
><br>
> Dear Espresso community<br>
><br>
> I have one problem that is important for me; it is somewhat surprising. I run geometry optimization (relax) calculations with pw.x with identical input files, and I obtain rather different final coordinates. The problem does not happen if ecutwfc is 30, but it does appear for cutoffs of 60, 80 or 90. Below one can see the exact input file, and examples of the difference between the final coordinates for both runs (both run with identical input) for given cutoffs. <br>
><br>
> Could anybody give me a clue on the origin of the problem, and how to solve it?<br>
><br>
> Thank you very much. Best regards.<br>
><br>
><br>
><br>
> Input file:<br>
><br>
>  <br>
><br>
> &CONTROL<br>
><br>
>     calculation = 'relax',<br>
><br>
>     restart_mode = 'from_scratch',<br>
><br>
>     prefix='',<br>
><br>
>     outdir = './',<br>
><br>
>     pseudo_dir = '/path_xxx/PP/',<br>
><br>
>     forc_conv_thr = 1.0D-6 ,<br>
><br>
>     etot_conv_thr = 1.0D-8 ,<br>
><br>
>  /<br>
><br>
> &system<br>
><br>
>     ibrav = 0, a=18.0,<br>
><br>
>     nat= 26, ntyp= 2,<br>
><br>
>     ecutwfc = 30d0,<br>
><br>
>     nbnd = 100,<br>
><br>
> /<br>
><br>
> &electrons<br>
><br>
>     conv_thr = 1.0e-9,<br>
><br>
>     mixing_beta = 0.7,<br>
><br>
>     mixing_mode = 'plain',<br>
><br>
>     diagonalization = 'cg'<br>
><br>
> /<br>
><br>
> &IONS<br>
><br>
> /<br>
><br>
>  <br>
><br>
> ATOMIC_SPECIES<br>
><br>
> C   12.0107   C.pz-vbc.UPF<br>
><br>
> H   1.007825035  H.pz-vbc.UPF<br>
><br>
>  <br>
><br>
> ATOMIC_POSITIONS { angstrom }<br>
><br>
> C  8.891700e+00     8.891700e+00     8.891700e+00<br>
><br>
> C  9.783400e+00     9.783400e+00     8.000000e+00<br>
><br>
> C  9.783400e+00     8.000000e+00     9.783400e+00<br>
><br>
> C  8.000000e+00     9.783400e+00     9.783400e+00<br>
><br>
> C  8.891700e+00     1.067510e+01     1.067510e+01<br>
><br>
> C  1.067510e+01     1.067510e+01     8.891700e+00<br>
><br>
> C  1.067510e+01     8.891700e+00     1.067510e+01<br>
><br>
> C  9.783400e+00     1.156680e+01     9.783400e+00<br>
><br>
> C  9.783400e+00     9.783400e+00     1.156680e+01<br>
><br>
> C  1.156680e+01     9.783400e+00     9.783400e+00<br>
><br>
> H  8.391700e+00     8.391700e+00     8.391700e+00<br>
><br>
> H  7.500000e+00     1.028340e+01     9.283400e+00<br>
><br>
> H  7.500000e+00     9.283400e+00     1.028340e+01<br>
><br>
> H  9.283400e+00     1.028340e+01     7.500000e+00<br>
><br>
> H  1.028340e+01     9.283400e+00     7.500000e+00<br>
><br>
> H  9.283400e+00     7.500000e+00     1.028340e+01<br>
><br>
> H  1.028340e+01     7.500000e+00     9.283400e+00<br>
><br>
> H  9.283400e+00     1.206680e+01     9.283400e+00<br>
><br>
> H  1.206680e+01     9.283400e+00     9.283400e+00<br>
><br>
> H  9.283400e+00     9.283400e+00     1.206680e+01<br>
><br>
> H  8.391700e+00     1.117510e+01     1.117510e+01<br>
><br>
> H  1.117510e+01     1.117510e+01     8.391700e+00<br>
><br>
> H  1.117510e+01     8.391700e+00     1.117510e+01<br>
><br>
> H  1.028340e+01     1.206680e+01     1.028340e+01<br>
><br>
> H  1.028340e+01     1.028340e+01     1.206680e+01<br>
><br>
> H  1.206680e+01     1.028340e+01     1.028340e+01<br>
><br>
>  <br>
><br>
>  <br>
><br>
> CELL_PARAMETERS {cubic}<br>
><br>
>  1.00  0.00  0.00<br>
><br>
>  0.00  1.00  0.00<br>
><br>
>  0.00  0.00  1.00<br>
><br>
>  <br>
><br>
> K_POINTS {gamma} {automatic}<br>
><br>
> 1 1 1  0 0 0<br>
><br>
>  <br>
><br>
>  <br>
><br>
> The program is run with<br>
><br>
>  <br>
><br>
> mpirun -np 32 /path_xxx/pw.x<br>
><br>
><br>
><br>
> Final coordinates: The first two rows for ecutwfc=30 for two runs with identical input are:<br>
><br>
> ATOMIC_POSITIONS (angstrom)<br>
> C  8.898432121  8.898432121  8.898432121 <br>
> C  9.783306562  9.783306562  8.019670107<br>
><br>
> ATOMIC_POSITIONS (angstrom)<br>
> C  8.898432121  8.898432121  8.898432121  <br>
> C  9.783306562  9.783306562  8.019670109<br>
><br>
><br>
> Final coordinates: The first two rows for ecutwfc=60 for two runs with identical input are:<br>
><br>
> ATOMIC_POSITIONS (angstrom)<br>
> C  8.904988579  8.904988579  8.904988579 <br>
> C  9.783426086  9.783426086  8.031809869<br>
><br>
> ATOMIC_POSITIONS (angstrom)<br>
> C  8.904962246  8.904962246  8.904962246<br>
> C  9.783425401  9.783425401  8.031847251<br>
><br>
> -- <br>
> --<br>
><br>
> Dr. Pablo García Risueño<br>
><br>
> Institut für Physikalische Chemie, Universität Hamburg, Grindelallee 117, 20146 Hamburg<br>
><br>
> Tel. +49 040 42 83 84 82 7<br>
><br></div></div><span>
> ______________________________<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" target="_blank"> http://pwscf.org/mailman/listi<wbr>nfo/pw_forum</a><br>
</span><p></p>
<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><br></blockquote></div><br><br clear="all"><br>-- <br><div class="m_6133018495245676239gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr">--<br><br>Dr. Pablo García Risueño<br><br>Institut für Physikalische Chemie, Universität Hamburg, Grindelallee 117, 20146 Hamburg<br><br>Tel. +49 040 42 83 84 82 7</div></div></div></div></div></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></div>