<div dir="ltr"><div class="gmail_quote"><div dir="ltr">Dear all,<br></div><div dir="ltr"><div>I'm experiencing a problem with total force minimization when the vdw_corr = 'TS' is used. Without the vdw_corr option, force minimization completes smoothly, as shown below:</div><div><br></div><div><div>     Total force =     0.144455     Total SCF correction =     0.002938</div><div>     Total force =     0.096620     Total SCF correction =     0.004258</div><div>     Total force =     0.023833     Total SCF correction =     0.000694</div><div>     Total force =     0.013572     Total SCF correction =     0.000196</div><div>     Total force =     0.013176     Total SCF correction =     0.000775</div><div>     Total force =     0.013520     Total SCF correction =     0.001316</div><div>     Total force =     0.013322     Total SCF correction =     0.000670</div><div>     Total force =     0.013164     Total SCF correction =     0.000897</div><div>     Total force =     0.013634     Total SCF correction =     0.000800</div><div>     Total force =     0.016822     Total SCF correction =     0.000639</div><div>     Total force =     0.015558     Total SCF correction =     0.000354</div><div>     Total force =     0.012359     Total SCF correction =     0.000858</div><div>     Total force =     0.009014     Total SCF correction =     0.000322</div><div>     Total force =     0.007140     Total SCF correction =     0.000338</div><div>     Total force =     0.006073     Total SCF correction =     0.000556</div><div>     Total force =     0.005047     Total SCF correction =     0.000465</div><div>     Total force =     0.004639     Total SCF correction =     0.000050</div><div>     Total force =     0.003113     Total SCF correction =     0.000382</div><div>     Total force =     0.002225     Total SCF correction =     0.000268</div><div>     Total force =     0.001698     Total SCF correction =     0.000156</div><div>     Total force =     0.001505     Total SCF correction =     0.000372</div><div>     Total force =     0.001422     Total SCF correction =     0.000244</div><div>     Total force =     0.001442     Total SCF correction =     0.000420</div><div>However, after implementing the VdW correction, it seems that the forces cannot be minimized further even after 200 ionic steps:</div><div><br></div><div><div>     Total force =     0.146620     Total SCF correction =     0.003488</div><div>     Total force =     0.108842     Total SCF correction =     0.004843</div><div>     Total force =     0.069519     Total SCF correction =     0.000923</div><div>     Total force =     0.016339     Total SCF correction =     0.002655</div><div>     Total force =     0.015353     Total SCF correction =     0.001041</div><div>     Total force =     0.016188     Total SCF correction =     0.000361</div><div>     Total force =     0.016350     Total SCF correction =     0.000832</div><div>     Total force =     0.015878     Total SCF correction =     0.001654</div><div>     Total force =     0.016273     Total SCF correction =     0.000560</div><div>     Total force =     0.016349     Total SCF correction =     0.000406</div><div>     Total force =     0.015879     Total SCF correction =     0.001403</div><div>     .</div><div>     .</div><div>     .</div><div>     .</div><div>     Total force =     0.015866     Total SCF correction =     0.001214</div><div>     Total force =     0.016279     Total SCF correction =     0.000263</div><div>     Total force =     0.016351     Total SCF correction =     0.000394</div><div>     Total force =     0.015879     Total SCF correction =     0.000394</div><div>     Total force =     0.016277     Total SCF correction =     0.001323</div><div>     Total force =     0.016347     Total SCF correction =     0.000194</div><div>     Total force =     0.015878     Total SCF correction =     0.000606</div></div><div>I was wondering if anyone can make a comment on this? I'm using QE 6.2 on a cluster. Input copied below:</div><div><div> &CONTROL</div><div>                       title = 'Li0.5NiO2 NCPP FM VdW ' ,</div><div>                 calculation = 'vc-relax' ,</div><div>                 max_seconds = 86300 ,</div><div>              !    wf_collect = .true. ,</div><div>                      outdir = '/' ,</div><div>                      wfcdir = '/' ,</div><div>                  pseudo_dir = '/scratch/05246/pseudo/' ,</div><div><span style="white-space:pre-wrap"> </span>!<span style="white-space:pre-wrap">       </span>restart_mode = 'restart' ,</div><div>                      prefix = 'LNO' ,</div><div>                     disk_io = 'low' ,</div><div>                   verbosity = 'high' ,</div><div><span style="white-space:pre-wrap">                       </span>nstep = 200 ,</div><div> /</div><div> &SYSTEM</div><div>                       ibrav = 0,</div><div>                   celldm(1) = 1.889726,</div><div>                         nat = 21,<br></div><div>                        ntyp = 3,</div><div>                     ecutwfc = 80 ,<br></div><div>                     ecutrho = 320 ,</div><div>                occupations = 'smearing' ,</div><div>                     degauss = 0.001 ,</div><div>                    smearing = 'gaussian' ,</div><div>                       nspin = 2 ,</div><div>   starting_magnetization(1) = 0.25,</div><div><span style="white-space:pre-wrap"> </span>             vdw_corr = 'TS' ,<br></div><div> /</div><div> &ELECTRONS</div><div>            electron_maxstep = 500,</div><div>                    conv_thr = 1.0d-6 ,</div><div>/</div><div> &IONS</div><div>/</div><div> &CELL</div><div>/</div><div>CELL_PARAMETERS {alat}</div><div>  5.790000   0.000000   0.000000</div><div> -1.447500   2.507144   0.000000</div><div>  0.000000   0.000000  14.221000</div><div>ATOMIC_SPECIES </div><div>   Ni   58.69000   Ni_ONCV_PBE-1.0.upf  </div><div>    O    15.99000  O_ONCV_PBE-1.0.upf </div><div>   Li    6.94000   Li_ONCV_PBE-1.0.upf </div><div>ATOMIC_POSITIONS crystal </div><div>O       0.000000   0.000000   0.232000</div><div>O       0.000000   0.000000   0.768000</div><div>O       0.333333   0.333333   0.565333</div><div>O       0.333333   0.333333   0.101333</div><div>O       0.166667   0.666667   0.898667</div><div>O       0.166667   0.666667   0.434667</div><div>Ni     0.000000   0.000000   0.500000</div><div>Ni     0.333333   0.333333   0.833333</div><div>Ni     0.166667   0.666667   0.166667</div><div>Li      0.166667   0.666667   0.666667</div><div>O       0.500000   0.000000   0.232000</div><div>O       0.500000   0.000000   0.768000</div><div>O       0.833333   0.333333   0.565333</div><div>O       0.833333   0.333333   0.101333</div><div>O       0.666667   0.666667   0.898667</div><div>O       0.666667   0.666667   0.434667</div><div>Ni     0.500000   0.000000   0.500000</div><div>Ni     0.833333   0.333333   0.833333</div><div>Ni     0.666667   0.666667   0.166667</div><div>Li      0.500000   0.000000   0.000000</div><div>Li      0.833333   0.333333   0.333333</div><div>K_POINTS automatic </div><div>4 8 2   0 0 0 </div></div><div><br></div><div>Thank you</div><div>P.S. I benchmarked the wf cut-off energy vs binary oxides for each metal and 80 Ry gave values comparable to 320 Ry.</div></div><div><br></div></div></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div><div><div><font size="1"><span style="color:rgb(153,153,153)"><b>Hooman Yaghoobnejad<br></b></span></font></div><font size="1"><span style="color:rgb(153,153,153)"><b>PhD, Department of Chemistry<br></b></span></font></div><font size="1"><span style="color:rgb(153,153,153)"><b>Missouri University of Science and Technology<br></b></span></font></div><font size="1"><span style="color:rgb(153,153,153)"><b>Rolla, MO 65409<br></b></span></font></div><font size="1"><span style="color:rgb(153,153,153)"><b>USA</b></span></font><br></div></div></div></div></div>