Dear all,<div><br></div><div>I'm trying to do finite electric field calculation for a thin slab with hexagonal symmetry. I consider a saw-like electric field (as large as 0.0005 au~0.03 eV/A) exerted along c-axis of my hexagonal unit cell. My input settings (partly) are as follows:</div>
<div><br></div><div><div> &CONTROL</div><div>                       title = Ta_slab ,</div><div>                 calculation = 'scf' ,</div><div>                restart_mode = 'from_scratch' ,</div><div>
                      outdir = './tmp/' ,</div><div>                  pseudo_dir = './' ,</div><div>                      prefix = 'tao' ,</div><div>                     tstress = .true. ,</div><div>
                     tprnfor = .true. ,</div><div>                     tefield = .true. ,</div><div>                    dipfield = .true. ,</div><div><br></div><div> /</div><div> &SYSTEM</div><div>                       ibrav = 4,</div>
<div>                   celldm(1) = 6.90415879,</div><div>                   celldm(3) = 21.57221207,</div><div>                         nat = 18,</div><div>                        ntyp = 3,</div><div>                     ecutwfc = 50 ,</div>
<div>                    noncolin = .true. ,</div><div>                    lspinorb = .true. ,</div><div>                        nbnd = 152,</div><div>                        edir = 3 ,</div><div>                        eamp = 0.0005834 ,</div>
<div>                     emaxpos = 0.75 ,</div><div>                      eopreg = 0.4  ,</div><div> /</div><div> &ELECTRONS</div><div>                    conv_thr = 1.0d-6 ,</div></div><div>.</div><div>.</div><div>.</div>
<div><br></div><div>Using these settings, the calculation appears to be going well for the first few iterations with rather sharp convergence in total energy, but then it becomes slower and slower and eventually ends up with a dE much larger than  conv_thr. Without electric field, there is no such a problem. Any idea how to improve the SCF convergence while tefield is on?</div>
<div><br></div><div>Any help is appreciated,</div><div>Sincerely</div><div>-- <br>Soroush Pakseresht<div>Advanced Science Institute (RIKEN)</div><div>Saitama-Japan</div><br>
</div>