<div dir="ltr"><div>I think there was a problem exactly with that case. It has been corrected recently:<br><br>
* ESM energy and forces for 'bc2' case and nonzero esm_efield were not<br>
correct (r13727). Also: problem with restart in NEB with ESM fixed<br><br></div>Paolo</div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Sep 14, 2017 at 2:39 PM, Michele Re Fiorentin <span dir="ltr"><<a href="mailto:michele.refiorentin@iit.it" target="_blank">michele.refiorentin@iit.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div><div><div><div>Dear all,<br><br></div>I tried to compare the results of the ESM example <span class="m_5648432870061838249gmail-final-path">Al001_bc2_efield.out (Al slab within a capacitor) with those you get when using tefield.<br></span></div><span class="m_5648432870061838249gmail-final-path">I can reproduce the reference values reported in the ESM example, getting a total energy of -49.2553 Ry.<br></span></div><span class="m_5648432870061838249gmail-final-path">Then, I setup a simulation using tefield ad dipfield (and no ESM of course) with the same amount of vacuum around the Al slab as in the ESM output and eamp set to half esm_efield (accounting for Ha/a.u instead of Ry/a.u.).<br></span></div><span class="m_5648432870061838249gmail-final-path">I get a final total energy of -49.5166, which is rather far from the ESM case. I attach here my two input files.<br></span></div><span class="m_5648432870061838249gmail-final-path">What am I missing?<br></span></div><span class="m_5648432870061838249gmail-final-path">Thank you very much for your cooperation.<br></span></div><span class="m_5648432870061838249gmail-final-path">Best,<br><br></span></div><span class="m_5648432870061838249gmail-final-path">Michele<br><br><br></span></div><span class="m_5648432870061838249gmail-final-path">***** ESM INPUT FILE *****<br>&control<br> calculation='scf',<br> restart_mode='from_scratch',<br> prefix='Al_esm_efieldM01',<br> outdir='./Al_esm_efieldM01/',<br> tprnfor = .TRUE.,<br> /<br> &system<br> ibrav = 0,<br> nat= 4, ntyp= 1,<br> ecutwfc = 20.0,<br> nosym=.TRUE.<br> occupations='smearing', <br> smearing='mp', <br> degauss=0.03<br> assume_isolated='esm', <br> esm_bc='bc2', <br> esm_efield= 0.00192148511256006<br> /<br> &electrons<br> mixing_beta = 0.3<br> /<br> &ions<br> /<br>ATOMIC_SPECIES<br> Al 26.981538 Al.pbe-n-van.UPF<br>CELL_PARAMETERS bohr<br> 10.82227686 0.00000000 0.00000000<br> 0.00000000 10.82227686 0.00000000<br>0.00000000 0.00000000 22.67672253<br>ATOMIC_POSITIONS bohr<br>Al 0.00000000 0.00000000 0<br>Al 5.41113843 0.00000000 0<br>Al 0.00000000 5.41113843 0<br>Al 5.41113843 5.41113843 0<br>K_POINTS automatic<br>6 6 1 1 1 0<br></span></div><div><br><span class="m_5648432870061838249gmail-final-path"></span></div><div><span class="m_5648432870061838249gmail-final-path"><br></span></div><div><span class="m_5648432870061838249gmail-final-path"><br></span></div><span class="m_5648432870061838249gmail-final-path">**** tefield INPUT FILE ****<br><br>&control<br> calculation='scf',<br> restart_mode='from_scratch',<br> prefix='Al_efieldM01',<br> outdir='./Al_efieldM01/',<br> tprnfor = .true.,<br> tefield = .true.,<br> dipfield = .true.<br> /<br> &system<br> ibrav = 6,<br> celldm(1) = 10.82227686,<br> celldm(3) = 4.<a href="tel:(907)%20489-1972" value="+19074891972" target="_blank">19074891972</a>,<br> nat= 4, ntyp= 1,<br> ecutwfc = 20.0,<br> nosym=.true.<br> occupations='smearing', <br> smearing='mp', <br> degauss=0.03<br> edir = 3,<br> emaxpos = 0.9,<br> eopreg = 0.1,<br> eamp = -0.00097234533<br> /<br> &electrons<br> mixing_beta = 0.3<br> /<br> &ions<br> /<br>ATOMIC_SPECIES<br> Al 26.981538 Al.pbe-n-van.UPF<br>ATOMIC_POSITIONS bohr<br>Al 0.00000000 0.00000000 22.67672253<br>Al 5.41113843 0.00000000 22.67672253<br>Al 0.00000000 5.41113843 22.67672253<br>Al 5.41113843 5.41113843 22.67672253<br>K_POINTS automatic<br>6 6 1 1 1 0<br><br><br></span><div><div><div><div><div><div><div><div><div><div><div>
<div><br clear="all"><div><div><div class="m_5648432870061838249gmail_signature"><div dir="ltr"><font face="arial, helvetica, sans-serif" size="2">--<br>Michele Re Fiorentin, PhD<br> <br>Istituto Italiano di Tecnologia (IIT)<br>Center for Sustainable Future Technologies @ PoliTO<br>Corso Trento, 21<br>10129 Torino (Italy)<br>Tel.: <a href="tel:+39%20011%20509%201904" value="+390115091904" target="_blank">+39 011/5091904</a><br>email: <a href="mailto:michele.refiorentin@iit.it" target="_blank">michele.refiorentin@iit.it</a></font></div></div></div>
<div dir="ltr"><br></div></div></div></div></div></div></div></div></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><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div>