<div class="gmail_quote">
<div><font size="4"><font color="#000000" face="verdana,sans-serif">Dear all,</font></font></div>
<div><font size="4" face="Verdana"></font><font color="#000000"> </font></div>
<div><font color="#000000" size="4" face="Verdana">I meet a big problem when doing ESM calculations that all the calculations couldn't get converged.</font></div>
<div><font color="#000000" size="4" face="Verdana">The input is posted below and the efforts I made are,</font></div>
<div><font size="4" face="Verdana"></font><font color="#000000"> </font></div>
<div><font color="#000000" size="4" face="Verdana">1, using different degauss</font></div>
<div><font color="#000000" size="4" face="Verdana">2, using different ecutwfc</font></div>
<div><font color="#000000" size="4" face="Verdana">3, using different esm_bc</font></div>
<div><font color="#000000" size="4" face="Verdana">4, using different number of cores</font></div>
<div><font color="#000000" size="4" face="Verdana">5, using different tot_charge</font></div>
<div><font color="#000000" size="4" face="Verdana">6, searching on the pw_forum</font></div>
<div><font size="4" face="Verdana"></font><font color="#000000"> </font></div>
<div><font color="#000000" size="4" face="Verdana">It can be converged when tot_charge=0. But otherwise, it can't and just run out of all the steps. The strange thing is the total energies in each step vary quite differently. </font></div>


<div><font color="#000000" size="4" face="Verdana">For example, </font></div>
<div><font color="#000000" size="4" face="Verdana">at step #47, the total energy is -7776.16619281 Ry,</font></div>
<div><font color="#000000" size="4" face="Verdana">at step #52, the total energy is -5612.51355642 Ry, </font></div>
<div><font color="#000000" size="4" face="Verdana">but step #53, the total energy is -8825.13517388 Ry.</font></div>
<div><font color="#000000" size="4" face="Verdana">just doesn't have trend to converge.</font></div>
<div><font color="#000000" size="4" face="Verdana"></font> </div>
<div><font size="4" face="Verdana">For esm_bc, pbc and bc1 were tested, but none of them worked.</font></div>
<div><font color="#000000"> </font></div>
<div><font color="#000000" size="4" face="Verdana">Hope you can help me out of this issue.</font></div>
<div><font color="#000000" size="4" face="Verdana">Thanks a lot.</font></div>
<div><font size="4" face="Verdana"></font><font color="#000000"> </font></div>
<div><font color="#000000" size="4" face="Verdana">Best,</font></div>
<div><font color="#000000" size="4" face="Verdana">Xin</font></div>
<div><font color="#000000" size="4" face="Verdana">LBNL, CA</font></div>
<div><font size="4" face="Verdana"></font><font color="#000000"> </font></div>
<div><font color="#000000">##############INPUT, </font><font color="#000000"><a href="http://pw.in###################">pw.in###################</a></font></div>
<div><font color="#000000">&control<br>    calculation='scf'<br>    prefix='HOPG'<br>    pseudo_dir='/'<br>    outdir='./'<br>    wf_collect=.TRUE.<br>    disk_io='low'<br>    restart_mode='from_scratch'<br>

    <br>/ <br>&system<br>    ibrav=8 , a=4.2384 , b=2.447 , c=32.0875<br>    nat=16, ntyp=1, tot_charge=-2<br>    nbnd=96, occupations='smearing', degauss=0.0019<br>    ecutwfc=35, ecutrho=300 <br>    assume_isolated = 'esm', esm_bc='bc1'    <br>

/<br>&electrons<br>    electron_maxstep=300<br>    conv_thr=1.0d-7<br>    mixing_mode='local-TF'<br>    mixing_beta=0.1<br>    <br>/ <br>ATOMIC_SPECIES<br>C 12.00 C.pbe-van_bm.UPF<br>ATOMIC_POSITIONS angstrom<br>

C      0.000000000    1.223500013   11.521874905 <br>C      1.412785889    1.223500013   11.521874905 <br>C      2.119199991    0.000000000   11.521874905 <br>C      3.531985754    0.000000000   11.521874905 <br>C      0.000000000    1.223500013   14.834999979 <br>

C      0.706414102    0.000000000   14.834999979 <br>C      2.119199991    0.000000000   14.834999979 <br>C      2.825614220    1.223500013   14.834999979 <br>C      0.000000000    1.223500013   18.148125053 <br>C      1.412785889    1.223500013   18.148125053 <br>

C      2.119199991    0.000000000   18.148125053 <br>C      3.531985754    0.000000000   18.148125053 <br>C      0.000000000    1.223500013   21.461249470 <br>C      0.706414102    0.000000000   21.461249470 <br>C      2.119199991    0.000000000   21.461249470 <br>

C      2.825614220    1.223500013   21.461249470 <br>K_POINTS automatic<br> 17 31 1 0 0 0<br>##############################################</font></div></div><br>