<div dir="ltr">Dear QE users,<div>  I am trying to perform vc-relax calculation with dft+u for Ce doped LaPO4 (Quantum Espresso Version: q-e 5.1.1) and I am getting unknown error after first scf is completed and program tries to calculate stress and forces. I am writing the input and part of output file. </div><div><br></div><div>Can someone suggest me what is wrong with my input file and help me correcting it?</div><div><br></div><div><br></div><div>Thanks and Regards</div><div><br></div><div>Input file: </div><div><br></div><div><div>&control</div><div>    calculation = 'vc-relax'</div><div>    restart_mode='from_scratch',</div><div>    prefix='Ce-LaPO4',</div><div>    pseudo_dir ='./',</div><div>    outdir='./'</div><div>    etot_conv_thr = 1.0D-4 ,</div><div>    forc_conv_thr = 1.0D-5 ,</div><div>    nstep = 200 ,</div><div>    tstress = .true. ,</div><div>    tprnfor = .true. ,</div><div> /</div><div> &SYSTEM</div><div>                       ibrav= 0,</div><div>                   celldm(1) = 13.0,</div><div>                         nat = 24,</div><div>                        ntyp = 4,</div><div>                     ecutwfc = 100,</div><div>                     ecutrho = 500,</div><div>                 occupations = 'smearing' ,</div><div>                     degauss = 0.05 ,</div><div>                    smearing = 'methfessel-paxton' ,</div><div>                  lda_plus_u = .true. ,</div><div>                Hubbard_U(2) = 0.5 ,</div><div><br></div><div> /</div><div> &ELECTRONS</div><div>                    conv_thr = 1D-8 ,</div><div>                 mixing_mode = 'TF' ,</div><div>                 mixing_beta = 0.7 ,</div><div> /</div><div> &IONS</div><div>                ion_dynamics = 'bfgs' ,</div><div> /</div><div> &CELL</div><div>               cell_dynamics = 'bfgs' ,</div><div> /</div><div>CELL_PARAMETERS alat</div><div>   1.008440563   0.000150030  -0.000707679</div><div>   0.000156146   1.037667024  -0.000743249</div><div>  -0.226206279  -0.000695433   0.921551306</div><div>ATOMIC_SPECIES</div><div>   La  138.90000  La.pbe-n-bpaw.UPF</div><div>   Ce  140.11600  Ce.GGA-PBE-paw-v1.0.UPF</div><div>    P   30.97300  P.pbe-n-kjpaw_psl.0.1.UPF</div><div>    O   15.99900  O.pbe-n-kjpaw_psl.0.1.UPF</div></div><div><div>ATOMIC_POSITIONS crystal</div><div>La       0.284838828   0.158253303   0.098300462</div><div>La       0.714536203   0.841484885   0.901173763</div><div>La       0.214814257   0.658205334   0.402230013</div><div>Ce       0.786626475   0.342226342   0.596739916</div><div>P        0.303273015   0.161632905   0.610568451</div><div>P        0.696786814   0.839904425   0.389731049</div><div>P        0.196338454   0.660568442   0.889571945</div><div>P        0.803629132   0.337924352   0.110095440</div><div>O        0.249488041   0.006113849   0.440535214</div><div>O        0.749904274   0.996627918   0.558841456</div><div>O        0.250658325   0.505658401   0.059644684</div><div>O        0.749910831   0.492042692   0.938777940</div><div>O        0.381694123   0.328988749   0.498608502</div><div>O        0.620480629   0.672521537   0.503747032</div><div>O        0.119023170   0.828565349   0.001936043</div><div>O        0.880537528   0.169145778   0.998817825</div><div>O        0.467991298   0.102102220   0.805713568</div><div>O        0.530191034   0.897120766   0.195217822</div><div>O        0.029728485   0.601755089   0.695280125</div><div>O        0.970086008   0.397279910   0.304172329</div><div>O        0.126638385   0.213230905   0.706276196</div><div>O        0.872495834   0.788789541   0.292888819</div><div>O        0.371878660   0.711047196   0.792338643</div><div>O        0.628450197   0.288810079   0.208792769</div><div><br></div><div><br></div><div>K_POINTS automatic</div><div>  8 8 8 0 0 0</div></div><div><br></div><div><br></div><div>Output file: </div><div><br></div><div><br></div><div><br></div><div><div>     the Fermi energy is    10.7359 ev</div><div><br></div><div>!    total energy              =   -2748.50526154 Ry</div><div>     Harris-Foulkes estimate   =   -2748.50526154 Ry</div><div>     estimated scf accuracy    <          3.9E-09 Ry</div><div><br></div><div>     total all-electron energy =     -7004.709326 Ry</div><div><br></div><div>     The total energy is the sum of the following terms:</div><div><br></div><div>     one-electron contribution =    -170.14060729 Ry</div><div>     hartree contribution      =     196.91626823 Ry</div><div>     xc contribution           =    -181.65417944 Ry</div><div>     ewald contribution        =    -672.75077535 Ry</div><div>     Hubbard energy            =       0.02176371 Ry</div><div>     one-center paw contrib.   =   -1920.88383324 Ry</div><div>     smearing contrib. (-TS)   =      -0.01389816 Ry</div><div><br></div><div>     convergence has been achieved in  23 iterations</div><div><br></div><div>     Forces acting on atoms (Ry/au):</div><div><br></div><div>     atom    1 type  1   force =     0.00009711   -0.00000294    0.00005862</div><div>     atom    2 type  1   force =    -0.00010873   -0.00007113    0.00007847</div><div>     atom    3 type  1   force =    -0.00004895    0.00005462    0.00000052</div><div>     atom    4 type  2   force =    -0.00032522   -0.00010967    0.00002218</div><div>     atom    5 type  3   force =    -0.00013072    0.00001515   -0.00007378</div><div>     atom    6 type  3   force =     0.00002085   -0.00000456    0.00002304</div><div>     atom    7 type  3   force =    -0.00006498   -0.00007657   -0.00004950</div><div>     atom    8 type  3   force =    -0.00007293    0.00022264    0.00009364</div><div>     atom    9 type  4   force =    -0.00004581    0.00000134    0.00000356</div><div>     atom   10 type  4   force =     0.00002917   -0.00059913   -0.00000776</div><div>     atom   11 type  4   force =     0.00001827    0.00001823    0.00004762</div><div>     atom   12 type  4   force =     0.00003066    0.00009104    0.00013979</div><div>     atom   13 type  4   force =    -0.00000072   -0.00008827    0.00004357</div><div>     atom   14 type  4   force =    -0.00010515    0.00026136   -0.00018367</div><div>     atom   15 type  4   force =    -0.00006743    0.00007268    0.00001610</div><div>     atom   16 type  4   force =     0.00005974    0.00008529   -0.00013874</div><div>     atom   17 type  4   force =    -0.00004507    0.00003833   -0.00004261</div><div>     atom   18 type  4   force =    -0.00003720   -0.00001576   -0.00003803</div><div>     atom   19 type  4   force =     0.00035582    0.00032602    0.00016180</div><div>     atom   20 type  4   force =     0.00008210   -0.00001638   -0.00012772</div><div>     atom   21 type  4   force =     0.00036060   -0.00015974    0.00002964</div><div>     atom   22 type  4   force =     0.00007277   -0.00005677    0.00010798</div><div>     atom   23 type  4   force =    -0.00001852   -0.00002585    0.00009441</div><div>     atom   24 type  4   force =    -0.00005567    0.00004007   -0.00025913</div><div><br></div><div>     Total force =     0.001168     Total SCF correction =     0.000122</div><div>     SCF correction compared to forces is large: reduce conv_thr to get better values</div><div><br></div><div><br></div><div>     entering subroutine stress ...</div><div><br></div><div><br></div><div> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%<wbr>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<wbr>%%%%%%%%%%%%%%%%%%%</div><div>     Error in routine gen_at_dj (1):</div><div>     unexpected error</div><div> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%<wbr>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<wbr>%%%%%%%%%%%%%%%%%%%</div><div><br></div><div>     stopping ...</div></div><div><br></div></div>