<div dir="ltr">Dear all,<div><br></div><div>I'm facing a bad convergence issue for a very simple system consisting of a WS2 bilayer.</div><div><br></div><div>While with NO vdW interaction turned on everything works fine, as I use the input_dft variable the optimization task fails. In practice, the pw.x tries to optimize the bilayer for some steps, when the Total force decreases, then the Total force starts increasing, until at some point (=geometry), the electron minimization does not converge within electron_maxstep = 300 steps.</div><div><br></div><div>Here, a typical input file I use:</div><div><br></div><div>&CONTROL<br>    calculation   = 'relax'<br>    title         = 'WS1-2L'<br>    restart_mode  = 'from_scratch'<br>    pseudo_dir   = '/m100_work/IscrB_TWITTER/pseudo/'<br>    outdir        = './tmp/'<br>    prefix        = 'WS2-2L'<br>    verbosity     = 'high'<br>    tprnfor       = .true.<br>    tstress       = .true.<br>    !etot_conv_thr = 1.D-5<br>    forc_conv_thr = 1.D-6<br>/<br>&SYSTEM<br>    ibrav             = 4<br>    a                 = __A<br>    c                 = 24.450<br>    nat               = 6<br>    ntyp              = 2<br>    ecutwfc           = 50<br>    ecutrho           = 500<br>    occupations       = 'smearing'<br>    smearing          = 'mv'<br>    degauss           = 0.00735<br>    input_dft         = 'vdw-df2-b86r'<br>/<br>&ELECTRONS<br>    conv_thr         = 1.d-8<br>    electron_maxstep = 300<br>    mixing_beta      = 0.15<br>/<br>&IONS<br>/<br>ATOMIC_SPECIES<br>        W    183.84       W.pbe-spn-kjpaw_psl.1.0.0.UPF<br>        S     32.065      S.pbe-n-kjpaw_psl.1.0.0.UPF<br>ATOMIC_POSITIONS { crystal }<br>S    0.666666667   0.333333333   0.080882686    <br>W    0.333333333   0.666666667   0.145220672    <br>S    0.666666667   0.333333333   0.209558657    <br>S    0.333333333   0.666666667   0.371324029    <br>W    0.666666667   0.333333333   0.435662015    <br>S    0.333333333   0.666666667   0.500000000    <br>K_POINTS { automatic }<br>12  12  1    0  0  0<br><br><br>I tried to decrease the beta (down to 0.01), to increase degauss (even though the system is semiconducting and this should not affect the convergence, since the energy gap is always much larger than degauss), to increase the k-point grid, to increase or decrease the vacuum, to increase ecutrho.</div><div><br></div><div>I suspect some stupid error I'm doing and that I cannot see, since most of the above listed recipies has worked for much more difficult systems that a simple WS2 bilayer, unless WS2 requires some special care that I ignore, or unless there is some issue with the W and/or S pseudopotential.</div><div>I also tried to run the run with no vdw and then that with vdw as a restart of the previous one.</div><div><br></div><div>WS2 bilayer seems for example to have been successfully studied with several (and even more recent) vdW functionals for example in this paper: <a href="https://pubs.acs.org/doi/pdf/10.1021/acs.jctc.0c00471">https://pubs.acs.org/doi/pdf/10.1021/acs.jctc.0c00471</a></div><div><br></div><div>I would appreciate any hint/help.</div><div><br></div><div>Thank you.</div><div>Giovanni</div><div><br clear="all"><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><span style="color:rgb(0,0,0)">-- <br><br>Giovanni Cantele, PhD<br>CNR-SPIN<br>c/o Dipartimento di Fisica<br>Universita' di Napoli "Federico II"<br>Complesso Universitario M. S. Angelo - Ed. 6<br>Via Cintia, I-80126, Napoli, Italy<br><a href="mailto:giovanni.cantele@spin.cnr.it" style="color:rgb(17,85,204)" target="_blank">e-mail: giovanni.cantele@spin.cnr.it</a><br>Phone: +39 081 676910<br>Skype contact: giocan74<br><br>ResearcherID: <a href="http://www.researcherid.com/rid/A-1951-2009" style="color:rgb(17,85,204)" target="_blank">http://www.researcherid.com/rid/A-1951-2009</a><br>Web page: </span><a href="https://sites.google.com/view/giovanni-cantele/home" style="color:rgb(17,85,204)" target="_blank">https://sites.google.com/view/giovanni-cantele/home</a><br></div></div></div></div></div>