<div dir="ltr">The error I have got is:<div>  <div>Message from routine bfgs:<br>     history already reset at previous step: exiting<br>     Energy error            =      3.4E-09 Ry<br>     Gradient error          =      4.0E-05 Ry/Bohr<br>     Cell gradient error     =      2.8E-03 kbar<br><br>     bfgs failed after  25 scf cycles and  23 bfgs steps, convergence not achieved<br>     (criteria: energy <  1.0E-06 Ry, force <  1.0E-05 Ry/Bohr, cell <  5.0E-03 kbar)<br><br>     End of BFGS Geometry Optimization<br><br>     Final enthalpy           =    -

**  Ry<br><br>     File ./vc/g.bfgs deleted, as requested<br>Begin final coordinates<br>     new unit-cell volume =    

**  a.u.^3 (   **Ang^3 )<br>     density =   

**   g/cm^3<br><br>CELL_PARAMETERS (alat= )<br><br><br>ATOMIC_POSITIONS (crystal)<br><br>End final coordinates<br>---------<br><br>I get this message for a spin orbit coupling calculations, I had used all relativistic pseudos. Input file description is as follows:</div><div><br></div><div><br></div><div>&CONTROL<br>    calculation   = "vc-relax"<br>    etot_conv_thr =  1.00e-06  <br>    forc_conv_thr =  1.00e-05  <br>    nstep         =  100       <br>    outdir        = "./vc"<br>    prefix        = "g"<br>    pseudo_dir    = "./"<br>    verbosity     = "high"<br>/<br><br>&SYSTEM<br>    celldm(1)   =   <br>    degauss     =  0.000735 <br>    ecutrho     =  440<br>    ecutwfc     =  110  <br>    ibrav       =  0<br>    nat         =  <br>    ntyp        =  2<br>    occupations = "smearing"<br>    smearing    = "gaussian"<br>    starting_magnetization(1) =  0.00000<br>    starting_magnetization(2) =  1.00000 <br>    noncolin = .true.               <br>    lspinorb = .true.         <br>/<br><br>&ELECTRONS<br>    conv_thr         =  1.00000e-08<br>    electron_maxstep =  200<br>    mixing_beta      =  7.00000e-01<br>/<br><br>&IONS<br>    ion_dynamics = "bfgs"<br>/<br><br>&CELL<br>    cell_dofree    = "2Dxy"<br>    cell_dynamics  = "bfgs"<br>    press          =  0.00000e+00<br>    press_conv_thr =  5.00000e-03<br>/<br><br>K_POINTS {automatic}<br> 8  8  1  0 0 0     <br><br>CELL_PARAMETERS {alat}<br><br><br>ATOMIC_SPECIES<br>A      *******      A.upf   !using norm-conserving and full relativistic with pbe <br>B      

*******     B.upf   !using norm-conserving and full relativistic with pbe <br><br>ATOMIC_POSITIONS {crystal}<br>  <br><br><div>What should I do next for optimization to be converged?</div><div>Also, please provide some guidelines for applying the magnetization value. How much should it be used?</div><div>If I have two types of atoms and put zero for starting_magnetization for both, how is it taken? If magnetization is given for one type of atom then how is it calculated? And if startine_magnetization is given for both, then what does the code consider?</div><div><br></div><div><br></div><div>Thank you.</div><div>Ashley Cooper</div><div>Phd<br><br></div></div></div></div>