<div dir="ltr"><div>Hello everyone,</div><div><br></div><div>I am trying to optimize this structure . I have run successful vc-relax optimization with low criteria (etot_conv_thr = 1.0D-5, forc_conv_thr = 1.0D-4) ,but when I increase the criteria bfgs fails. This structure exists in literature with buckling that I can't seem to get here.Due to bfgs failure in vc-relax,  I had decided to run a relax followed by a vc-relax. </div><div>In my current relax calculation the optimization is stuck at bfgs 7 with scf step 138. My question is if I want to successfully relax this structure (I would be doing phonon afterwards ) what should I do here ? Should I stop the current run and rerun ? Please let me know the best solution.</div><div><br></div>Input:<div>&CONTROL<br>    calculation   = "relax"<br>    restart_mode  = "from_scratch"<br>    outdir="/tmp/"<br>    prefix="p"<br>    pseudo_dir="./"<br>    verbosity="high"<br>    etot_conv_thr = 1.0D-6<br>    forc_conv_thr = 1.0D-4<br>    nstep           = 1000<br>    tprnfor       = .true.<br>    disk_io       = "none"<br>/<br><br>&SYSTEM<br>  ibrav           = 0<br>  celldm(1) = 1<br>  nat         = 8<br>  ntyp          = 2<br>  ecutwfc   = 210<br>/<br><br>&ELECTRONS<br>    electron_maxstep     = 600<br>    mixing_beta        = 0.2<br>    diagonalization    = "cg"<br>    diago_cg_maxiter        = 500<br>    conv_thr           =  1.0d-5<br>/<br><br>&IONS<br>    ion_dynamics = "bfgs"<br>/<br><br>ATOMIC_SPECIES<br> B 00d0 B.pbe-n-rrkjus_psl.1.0.0.UPF<br> A 00d0 A.pbe-dnl-kjpaw_psl.1.0.0.UPF<br><br>CELL_PARAMETERS (alat=  1.00000000)<br>  15.342766198   0.000000000   0.000000000<br>   0.000000000  15.342766198   0.000000000<br>   0.000000000   0.000000000  37.794630380<br><br>ATOMIC_POSITIONS (crystal)<br>A               0.6461706388        0.1462855988        0.5000000000<br>B               0.6495390386        0.8505502843        0.5000000000<br>B               0.3504609614        0.1494497157        0.5000000000<br>A               0.3538293612        0.8537144012        0.5000000000<br>A               0.1462855988        0.3538293612        0.5000000000<br>B               0.8505502843        0.3504609614        0.5000000000<br>B               0.1494497157        0.6495390386        0.5000000000<br>A               0.8537144012        0.6461706388        0.5000000000<br><br>K_POINTS {automatic}<br>  14 14 3 0 0 0<br><br><br><div><br></div><div><br></div><div>Current situation of output:</div><div><br></div><div>     Total force =     0.000970     Total SCF correction =     0.002258<br>     SCF correction compared to forces is large: reduce conv_thr to get better values<br>     Energy error            =      7.3E-08 Ry<br>     Gradient error          =      3.9E-04 Ry/Bohr<br><br>     number of scf cycles    = 138<br>     number of bfgs steps    =   7<br><br>     energy             old  =   -1193.3574486014 Ry<br>     energy             new  =   -1193.3574485288 Ry<br><br>     CASE: energy            _new > energy            _old<br><br>     new trust radius        =       0.0000504644 bohr<br><br>     trust_radius < trust_radius_min<br><br>     resetting bfgs history<br><br>     new conv_thr            =       0.0000001000 Ry<br></div></div><div><br></div><div><br></div><div><br></div><div>Kind Regards,</div><div>Ashley Cooper</div><div>Research Scholar</div><div>Assam University Silchar</div><div>Assam, India</div><div><br></div></div>