<div dir="ltr"><div dir="ltr"><div dir="ltr">With regards<div>Maybe my question is very repetitive but i could found nothing in the forum to help. I am trying to run a vc-relax for cementite structure using GGA-USPP. The problem is that i cannot achieve any convergence even after 1000 iteration. The total energy fluctuates in the order of ~0.1! By using other PP, i simply achieved to the convergence, but i need to use these pp ones. I would appreciate any suggestions in this regard.</div><div>This is my input:</div><div><div>&CONTROL</div><div>    calculation   = "vc-relax"</div><div>    forc_conv_thr =  1.00000e-03</div><div>    max_seconds   =  8.64000e+07</div><div>    nstep         = 1000</div><div>    pseudo_dir    = "..../bin"</div><div>/</div><div><br></div><div>&SYSTEM</div><div>    a                         =  5.03018e+00</div><div>    angle1(1)                 =  0.00000e+00</div><div>    angle2(1)                 =  0.00000e+00</div><div>    b                         =  6.73931e+00</div><div>    c                         =  4.49084e+00</div><div>    degauss                   =  1.00000e-02</div><div>    ecutrho                   =  6.37306e+02</div><div>    ecutwfc                   =  7.08118e+01</div><div>    ibrav                     = 8</div><div>    lda_plus_u                = .FALSE.</div><div>    nat                       = 16</div><div>    nspin                     = 2</div><div>    ntyp                      = 2</div><div>    occupations               = "smearing"</div><div>    smearing                  = "gaussian"</div><div>    starting_magnetization(1) =  4.00000e-01</div><div>    starting_magnetization(2) =  0.00000e+00</div><div>/</div><div><br></div><div>&ELECTRONS</div><div>    conv_thr         =  1.00000e-06</div><div>    electron_maxstep = 2000</div><div>    mixing_beta      =  4.00000e-01</div><div>    startingpot      = "atomic"</div><div>    startingwfc      = "atomic+random"</div><div>/</div><div><br></div><div>&IONS</div><div>    ion_dynamics = "bfgs"</div><div>/</div><div><br></div><div>&CELL</div><div>    cell_dofree    = "xyz"</div><div>    cell_dynamics  = "bfgs"</div><div>    press          =  0.00000e+00</div><div>    press_conv_thr =  5.00000e-01</div><div>/</div><div><br></div><div>K_POINTS {automatic}</div><div> 2  2  3  0 0 0</div><div><br></div><div>ATOMIC_SPECIES</div><div>Fe     55.84500  Fe.pbe-spn-rrkjus_psl.1.0.0.UPF</div><div>C      12.01070  C.pbe-n-rrkjus_psl.1.0.0.UPF</div><div><br></div><div>ATOMIC_POSITIONS {angstrom}</div><div>Fe      0.892203   0.460793   1.491414</div><div>Fe      4.137977   6.278514   2.999431</div><div>Fe      1.622887   6.278514   3.736836</div><div>Fe      3.407293   0.460793   0.754008</div><div>Fe      3.407293   2.908860   0.754008</div><div>Fe      1.622887   3.830447   3.736836</div><div>Fe      4.137977   3.830447   2.999431</div><div>Fe      0.892203   2.908860   1.491414</div><div>Fe      0.184411   1.684827   3.759376</div><div>Fe      4.845769   5.054480   0.731469</div><div>Fe      2.330679   5.054480   1.513954</div><div>Fe      2.699502   1.684827   2.976891</div><div>C       0.618742   5.054480   2.519058</div><div>C       4.411438   1.684827   1.971786</div><div>C       1.896348   1.684827   0.273636</div><div>C       3.133833   5.054480   4.217209</div></div><div><br></div><div><div style="color:rgb(80,0,80)">Thanks</div><div style="color:rgb(80,0,80)"><br></div><div style="color:rgb(80,0,80)">Seyed Mehrdad Zamzamian</div><div style="color:rgb(80,0,80)">Sharif University of Technology, Tehran, Iran</div><div style="color:rgb(80,0,80)">Energy engineering department</div><div style="color:rgb(80,0,80)">E-mail: <a href="mailto:mehrdad.zamzamian@gmail.com" target="_blank">mehrdad.zamzamian@gmail.com</a></div></div></div></div></div>