<div dir="ltr"><div>Dear all,</div><div><br></div><div>I am (again...) fighting with the dipole correction. I have a system consisting of Ag and MgO centered in the cell around z=0.5</div><div><br></div><div>ATOMIC_POSITIONS (crystal)<br>Ag       0.000000000   0.000000000   0.323271580<br>Ag       0.000000000   0.500000000   0.385512833<br>....<br>Mg      -0.000000000   0.000000000   0.721978968<br>O        0.500000000   0.000000000   0.723372273<br>O       -0.000000000   0.500000000   0.723371161<br>Mg       0.500000000   0.500000000   0.721975961<br></div><div><br></div><div>Then I apply tefield and dipfield with</div><div><br></div><div>  eamp        = 0.00<br>  edir        = 3<br>  emaxpos     = 0.95<br>  eopreg      = 0.10<br></div><div><br></div><div>which means my system has no discontinuity anywhere in the slab and the field is "physical" between 0.95+0.1-1=0.05 and 0.95 of the z direction.</div><div><br></div><div>The SCF cycle converges but the pot_num=11 never becomes horizontal, rather it has a very small slope. I tried to increase the dipole length to see if that makes a difference but apparently it does not help either. With above input the potential amplitude at scf convergence is <br></div><div><br></div><div>     Adding external electric field<br><br>     Computed dipole along edir(3) : <br>        Elec. dipole          0.4894 Ry au,          1.2440 Debye<br>        Ion. dipole           0.4902 Ry au,          1.2460 Debye<br>        Dipole                0.2274 Ry au,          0.5779 Debye<br>        Dipole field          0.0008 Ry au, <br> <br>        Potential amp.       -0.0855 Ry<br>        Total length         55.5874 bohr<br><br></div><div><br></div><div>Is there a scenario where the dipole correction could fail if the potential is "almost flat" in the vacuum region? If I run the same with VASP I get a small "jump" in the potential making it flat in the vacuum region, so I was wondering why this does not happen (probably different approach in PWSCF)?</div><div><br></div><div>Input and resulting potential is attached - any help is welcome! <br></div><div><br></div><div>Thanks in advance and with best wishes,</div><div>Chris <br></div><div>-- <br><div class="gmail_signature"><div dir="ltr">Postdoctoral Researcher<br>Center for Quantum Nanoscience, Institute for Basic Science<br>Ewha Womans University, Seoul, South Korea<blockquote type="cite" style="font-size:12.8px"><div dir="ltr"><div><div dir="ltr"></div></div></div></blockquote></div></div>
</div></div>