<div dir="ltr">Maybe the following note in file Doc/release-notes is relevant for your problem?<div>---</div><div>...</div><div>Fixed in 6.8 version:</div><div>...<br>   * The rho => 0 limit of spin-polarized BEEF XC energy was not correct,<br></div>    leading to funny total energy numbers and problems in structural optimization<br>    (fixed by Gabriel S. Gusmão, Georgia Tech)<div>---</div><div><br></div><div>Paolo</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Nov 18, 2021 at 3:06 PM Tom Demeyere <<a href="mailto:T.Demeyere@soton.ac.uk">T.Demeyere@soton.ac.uk</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear users,<br>
<br>
I am trying to calculate binding energies using the BEEF-vdW and the RISM formalism (QE-6.7), for this aim I need to compute different molecules (H2O, O2, …) in H2O with HCl 1M.<br>
<br>
O2 shows large energy fluctuations when I ask for:<br>
<br>
- input_dft: “BEEF-vdW”<br>
- nspin = 2<br>
<br>
(The calculation converges when asking for nspin = 0, or changing the functional to PBE.)<br>
<br>
Typical output:<br>
<br>
…<br>
<br>
iteration #  8 <br>
...<br>
     total energy              =    -180.99829260 Ry<br>
     estimated scf accuracy    <       0.00029814 Ry<br>
…<br>
iteration #  9 <br>
...<br>
     total energy              =    -151.60027495 Ry<br>
     estimated scf accuracy    <       0.00021978 Ry<br>
…<br>
iteration # 10<br>
     total energy              =    -139.03397732 Ry<br>
     estimated scf accuracy    <       0.00010331 Ry<br>
…<br>
<br>
The calculation ends at some point with: “history already reset at previous step: stopping”, the final energy given is wrong.<br>
<br>
Full input:<br>
<br>
&CONTROL<br>
  calculation = 'relax'<br>
  etot_conv_thr = 1.0000000000d-05<br>
  forc_conv_thr = 5.0000000000d-04<br>
  pseudo_dir = '/work/e89/e89/td5g20/PseudoPotentials/'<br>
  tprnfor = .true.<br>
  verbosity = 'high'<br>
  trism = .true.<br>
/<br>
&SYSTEM<br>
  ecutwfc = 60, ecutrho = 480<br>
  occupations = ‘smearing’, smearing = ‘fd’, degauss = 0.01<br>
  ibrav = 1, A = 20<br>
  ntyp = 1, nat = 2<br>
  nosym = .true.<br>
  input_dft = "BEEF-vdW"<br>
/<br>
&ELECTRONS<br>
  conv_thr = 1.00D-12<br>
  mixing_mode = ‘local-TF’,  mixing_beta = 0.2<br>
/<br>
&IONS<br>
  ion_dynamics = "bfgs"<br>
/<br>
&RISM<br>
  nsolv = 3  <br>
  closure = 'kh'  <br>
  tempv = 298.15<br>
  ecutsolv = 160.0<br>
  starting1d = 'zero'  <br>
  rism1d_conv_thr = 1.0e-8  <br>
  rism1d_maxstep = 100000  <br>
  mdiis1d_size = 20  <br>
  mdiis1d_step = 0.1<br>
  starting3d = 'zero'<br>
  rism3d_maxstep = 10000  <br>
  rism3d_conv_thr = 1.0e-6  <br>
  rism3d_conv_level = 0.5  <br>
  mdiis3d_size = 20  <br>
  mdiis3d_step = 0.8<br>
  solute_lj(1) = 'none'<br>
  solute_epsilon(1) = 0.1554<br>
  solute_sigma(1) = 3.1660<br>
/<br>
ATOMIC_SPECIES<br>
O   15.999 O.pbe-n-kjpaw_psl.1.0.0.UPF<br>
ATOMIC_POSITIONS angstrom<br>
O             1.6891044673        1.7570200245        0.6442147904<br>
O             2.3377455327        0.7138399755        0.6439152096<br>
K_POINTS gamma<br>
SOLVENTS mol/L<br>
H2O  -1.0  H2O.spc.MOL<br>
H3O+  1.0  H3O+.aq.MOL<br>
Cl-   1.0  Cl-.aq.MOL<br>
<br>
Best regards,<br>
Tom Demeyere<br>
<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu" rel="noreferrer" target="_blank">www.max-centre.eu</a>)<br>
users mailing list <a href="mailto:users@lists.quantum-espresso.org" target="_blank">users@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/users</a></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 206, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>