<div dir="ltr">Dear all, <div><br></div><div>              I am trying to find the equilibrium lattice parameters of metallic alpha Uranium from LDA+U calculations using 'vc-relax'. I started with values close to the values obtained from other DFT studies, but it is not converging even after 100 iterations. Also, i am not sure whether i chose the right parameters for spin-orbit coupling (nspin=2, tot_mag=2). Please let me know if i am missing something. </div><div><br></div><div>Thanks</div><div>Ajit</div><div><br></div><div>****************************************************************</div><div><div>U</div><div>Uranium</div><div> &control</div><div>    calculation = 'vc-relax'</div><div>    restart_mode='from_scratch',</div><div>    prefix='UraniumLDAVC',</div><div>    tstress = .true.</div><div>    tprnfor = .true.</div><div>    pseudo_dir = './',</div><div>    outdir='./'</div><div>   etot_conv_thr=1.0E-4</div><div>    forc_conv_thr=1.0E-3 </div><div> /</div><div> &system    </div><div>    ibrav=  9,</div><div>    celldm(1) = 5.325316810814</div><div>    celldm(2) = 2.049741868141169</div><div>    celldm(3) = 1.742041343325397</div><div>    nat=  2, ntyp= 1,</div><div>    nbnd=24 </div><div>    ecutwfc =50.0</div><div>    smearing='mp'</div><div>    degauss=0.02 </div><div>    lda_plus_u=.TRUE.</div><div>    Hubbard_U(1)=3.66107</div><div>    nspin=2</div><div>    tot_magnetization=2 </div><div> /</div><div> &electrons</div><div>    electron_maxstep=100</div><div>    mixing_mode = 'plain'</div><div>    mixing_beta = 0.7 </div><div>    conv_thr =  1.0d-5</div><div> /</div><div>&ions</div><div> ion_dynamics='bfgs'</div><div>/</div><div>&CELL</div><div>  cell_dynamics='bfgs'</div><div>  press=0.01</div><div>/</div><div>ATOMIC_SPECIES</div><div> U 238.051 U.LDA-PW-paw.UPF</div><div><br></div><div>ATOMIC_POSITIONS {alat}</div><div> U 0.00 0.417311046418597 0.871020671662698 </div><div> U 1.00 0.00 0.00 </div><div><br></div><div>K_POINTS {automatic} </div><div>  12 12 12 0 0 0 </div></div><div><br></div></div>