<div>Dear Mostafa,</div><div><div>It happens to be that my research is focused on molecules in ultrahigh electric fields such as 1Ry or higher. The SCF calculations exhibit positive total energies in oscillations. My problem is how to get those calculations converged with your help. Do you or anyone have ideas to resolve this problem?</div><div><br></div><div>When I use option of <span style="  ; ; ;; ">occupations = 'smearing' ,</span> the calculations stop and complain that "Error in routine iosys(1): Berry Phase/electric fields only for insulators!". Do you have any ideas to crack this problem? What parameters should I specify?</div><div><br></div><div>Thank you very much!</div><div><br></div><div>Jibiao Li</div><div><br></div><div><div style="font-family: "lucida Grande", Verdana, "Microsoft YaHei";"><div>&CONTROL</div><div>                 calculation = 'scf' ,</div><div>                restart_mode = 'from_scratch' ,</div><div>                      outdir = '/home/jibiao/Calc/molecules/N2/' ,</div><div>                  pseudo_dir = '/home/jibiao/Codes/pseudo/' ,</div><div>                      prefix = 'N2' ,</div><div>                     tstress = .true. ,</div><div>                     tprnfor = .true. ,</div><div>                    lelfield = .true. ,</div><div> /</div><div> &SYSTEM</div><div>                       ibrav = 1,</div><div>                   celldm(1) = 15.<span t="7" data="117815022" style="border-bottom: 1px dashed rgb(204, 204, 204); z-index: 1;">117815022</span>,</div><div>                         nat = 2,</div><div>                        ntyp = 1,</div><div>                     ecutwfc = 29 ,</div><div>                     ecutrho = 180 ,</div><div><font color="#0000ff">                 occupations = 'fixed' ,</font></div><div> /</div><div> &ELECTRONS</div><div>            electron_maxstep = 300,</div><div>                 mixing_beta = 0.2D0 ,</div><div>             diagonalization = 'david' ,</div><div>              efield_cart(1) = 1,</div><div>              efield_cart(2) = 0,</div><div>              efield_cart(3) = 0,</div><div> /</div><div>ATOMIC_SPECIES</div><div>    N   14.<span t="7" data="00700" style="border-bottom: 1px dashed rgb(204, 204, 204); z-index: 1;">00700</span>  N.pbe-van_ak.UPF</div><div>ATOMIC_POSITIONS angstrom</div><div>    N      5.<span t="7" data="122061326" style="border-bottom: 1px dashed rgb(204, 204, 204); z-index: 1;">122061326</span>    3.<span t="7" data="999999675" style="border-bottom: 1px dashed rgb(204, 204, 204); z-index: 1;">999999675</span>    4.<span t="7" data="000000376" style="border-bottom: 1px dashed rgb(204, 204, 204); z-index: 1;">000000376</span>    </div><div>    N      4.<span t="7" data="000000000" style="border-bottom: 1px dashed rgb(204, 204, 204); z-index: 1;">000000000</span>    4.<span t="7" data="000000000" style="border-bottom: 1px dashed rgb(204, 204, 204); z-index: 1;">000000000</span>    4.<span t="7" data="000000000" style="border-bottom: 1px dashed rgb(204, 204, 204); z-index: 1;">000000000</span>    0  0  0</div><div>K_POINTS automatic</div><div>  1 1 1   0 0 0</div><div><br></div></div></div><div><br></div><div style="font-size: 12px;font-family: Arial Narrow;padding:2px 0 2px 0;">------------------ Original ------------------</div><div style="font-size: 12px;background:#efefef;padding:8px;"><div><b>From: </b> "Mostafa Youssef";<mostafa.youssef@aucegypt.edu>;</div><div><b>Date: </b> Tue, Aug 8, 2017 11:44 PM</div><div><b>To: </b> "pw_forum"<pw_forum@pwscf.org>; <wbr></div><div></div><div><b>Subject: </b> Re: [Pw_forum] SCF calculations for molecules in electric fieldsare not converged</div></div><div><br></div><div dir="ltr"><div><div>Dear Jibiao Li,<br><br></div><div>A field of 1 Ry a.u. is giant and atypical, so one should not be surprised do not converge.  Recall,  1 in Ry a.u.=36.3609*10^10 V/m.<br><br></div><div>I also noticed that you impose fixed occupations and then apply Gaussian smearing. Those are in contradiction of each other. Apparently the code enforces fixed occupations and ignores the smearing.<br></div><div><br><br></div>Mostafa<br></div>AUC, New Cairo, Egypt<br></div></div>