<div dir="ltr">Hi all dear QE users<div>I want to calculate the vacancy migration in alpha-Fe. I defined 2by2by2 supercell and used NEB but i could not achieve any convergence. searching in archives did not help. at the end of this message is my input and a part of the output. I tested various changes like k-point, ds,  num_of_images ... but nothing happened.  Please guide me if possible.</div><div>------------------------------------------------------</div><div>BEGIN<br>BEGIN_PATH_INPUT<br>&PATH<br>  restart_mode      = 'from_scratch'<br>  string_method     = 'neb',<br>  nstep_path        = 50,<br>  ds                = 1.D0,<br>  opt_scheme        = "broyden",<br>  num_of_images     = 12,<br>  CI_scheme         = "no-CI",<br>  path_thr          = 0.1,<br>/<br>END_PATH_INPUT<br>BEGIN_ENGINE_INPUT<br>&CONTROL<br>  prefix         = ""<br>  pseudo_dir     = "",<br>/<br>&SYSTEM<br>    a                         =  5.69623e+00<br>    degauss                   =  1.00000e-02<br>    ecutrho                   =  2.25000e+02<br>    ecutwfc                   =  2.50000e+01<br>    ibrav                     = 1<br>    nat                       = 15<br>    nspin                     = 2<br>    ntyp                      = 1<br>    occupations               = "smearing"<br>    smearing                  = "gaussian"<br>    starting_magnetization(1) =  2.00000e-01<br>/<br>&ELECTRONS<br>  mixing_beta = 0.01<br>  conv_thr =  1.0d-8<br>  electron_maxstep=200<br>/<br>&IONS<br>/<br>ATOMIC_SPECIES<br>Fe     55.84500  Fe.pbe-nd-rrkjus.UPF<br>BEGIN_POSITIONS<br>FIRST_IMAGE<br>ATOMIC_POSITIONS (angstrom)<br>Fe      0.000000   0.000000   0.000000<br>Fe      1.424058   1.424058   1.424058<br>Fe      0.000000   0.000000   2.848116<br>Fe      1.424058   1.424058   4.272174<br>Fe      0.000000   2.848116   0.000000<br>Fe      1.424058   4.272174   1.424058<br>Fe      0.000000   2.848116   2.848116<br>Fe      1.424058   4.272174   4.272174<br>Fe      2.848116   0.000000   0.000000<br>Fe      4.272174   1.424058   1.424058<br>Fe      2.848116   0.000000   2.848116<br>Fe      4.272174   1.424058   4.272174<br>Fe      2.848116   2.848116   0.000000<br>Fe      4.272174   4.272174   1.424058<br>Fe      2.848116   2.848116   2.848116<br>LAST_IMAGE<br>ATOMIC_POSITIONS (angstrom)<br>Fe      0.000000   0.000000   0.000000<br>Fe      1.424058   1.424058   1.424058<br>Fe      0.000000   0.000000   2.848116<br>Fe      1.424058   1.424058   4.272174<br>Fe      0.000000   2.848116   0.000000<br>Fe      1.424058   4.272174   1.424058<br>Fe      0.000000   2.848116   2.848116<br>Fe      1.424058   4.272174   4.272174<br>Fe      2.848116   0.000000   0.000000<br>Fe      4.272174   1.424058   1.424058<br>Fe      2.848116   0.000000   2.848116<br>Fe      4.272174   1.424058   4.272174<br>Fe      2.848116   2.848116   0.000000<br>Fe      4.272174   4.272174   1.424058<br>Fe      4.272174   4.272174   4.272174<br>END_POSITIONS<br>K_POINTS (automatic)<br> 5 5 5 0 0 0<br>END_ENGINE_INPUT<br>END<br></div><div>-----------------------------------------------------------------------</div><div>  tcpu =      0.4    self-consistency for image   1<br><br>     WARNING :  scf convergence NOT achieved on image   1<br><br>     cleaning-up extrapolation files<br></div><div>  -----------------------------------------------------------------------  <br></div><div>Best regards</div><div>S.M.Zamzamian</div><div>Iran, Tehran</div><div>Sharif University of Technology</div></div>