<div dir="ltr">Dear Sir,<div>      I am trying to do position optimization for BaMgF4 on 24 cores. I have performed scf calculation successfully using quantum espresso v 5.1.1. I am getting the following error after I run the code. Input file is given in the last. Please help me running this correctly.</div><div><br></div><div><div>     bravais-lattice index     =            9</div><div>     lattice parameter (alat)  =       7.8500  a.u.</div><div>     unit-cell volume          =    1227.8445 (a.u.)^3</div><div>     number of atoms/cell      =           24</div><div>     number of atomic types    =            3</div><div>     number of electrons       =       160.00</div><div>     number of Kohn-Sham states=           96</div><div>     kinetic-energy cutoff     =     100.0000  Ry</div><div>     charge density cutoff     =    1000.0000  Ry</div><div>     convergence threshold     =      1.0E-08</div><div>     mixing beta               =       0.7000</div><div>     number of iterations used =            8  plain     mixing</div><div>     Exchange-correlation      =  SLA  PW   PBX  PBC ( 1  4  3  4 0 0)</div><div>     nstep                     =           50</div><div><br></div><div><br></div><div>     celldm(1)=   7.850000  celldm(2)=   3.550000  celldm(3)=   1.430000</div><div>     celldm(4)=   0.000000  celldm(5)=   0.000000  celldm(6)=   0.000000</div><div><br></div><div>     crystal axes: (cart. coord. in units of alat)</div><div>               a(1) = (   0.500000   1.775000   0.000000 )</div><div>               a(2) = (  -0.500000   1.775000   0.000000 )</div><div>               a(3) = (   0.000000   0.000000   1.430000 )</div><div><br></div><div>     reciprocal axes: (cart. coord. in units 2 pi/alat)</div><div>               b(1) = (  1.000000  0.281690  0.000000 )</div><div>               b(2) = ( -1.000000  0.281690  0.000000 )</div><div>               b(3) = (  0.000000  0.000000  0.699301 )</div><div><br></div><div><br></div><div>     PseudoPot. # 1 for Ba read from file:</div><div>     /home2/yogesh/RRCAT/BaMgF4/relax/Ba.pbe-spn-kjpaw_psl.0.2.3.upf</div><div>     MD5 check sum: d0937a9c57b96273d14250f3d6e8b30f</div><div>     Pseudo is Projector augmented-wave + core cor, Zval = 10.0</div><div>     Generated using "atomic" code by A. Dal Corso  v.5.0.99 svn rev. 10869</div><div>     Shape of augmentation charge: PSQ</div><div>     Using radial grid of 1251 points,  6 beta functions with:</div><div>                l(1) =   0</div><div>                l(2) =   0</div><div>                l(3) =   1</div><div>                l(4) =   1</div><div>                l(5) =   2</div><div>                l(6) =   2</div><div>     Q(r) pseudized with 0 coefficients</div><div><br></div><div><br></div><div>     PseudoPot. # 2 for Mg read from file:</div><div>     /home2/yogesh/RRCAT/BaMgF4/relax/Mg.pbe-n-kjpaw_psl.0.3.0.upf</div><div>     MD5 check sum: 9eaefceb6cf3903db1b93e4e5c6adf9a</div><div>     Pseudo is Projector augmented-wave + core cor, Zval =  2.0</div><div>     Generated using "atomic" code by A. Dal Corso  v.5.0.99 svn rev. 10869</div><div>     Shape of augmentation charge: PSQ</div><div>     Using radial grid of 1129 points,  4 beta functions with:</div><div>                l(1) =   0</div><div>                l(2) =   0</div><div>                l(3) =   1</div><div>                l(4) =   1</div><div>     Q(r) pseudized with 0 coefficients</div><div><br></div><div><br></div><div>     PseudoPot. # 3 for  F read from file:</div><div>     /home2/yogesh/RRCAT/BaMgF4/relax/F.pbe-n-kjpaw_psl.0.1.upf</div><div>     MD5 check sum: 0bca836995e818da6ac81be02c66119a</div><div>     Pseudo is Projector augmented-wave + core cor, Zval =  7.0</div><div>     Generated using "atomic" code by A. Dal Corso  v.5.0.99 svn rev. 10869</div><div>     Shape of augmentation charge: PSQ</div><div>     Using radial grid of 1105 points,  4 beta functions with:</div><div>                l(1) =   0</div><div>                l(2) =   0</div><div>                l(3) =   1</div><div>                l(4) =   1</div><div>     Q(r) pseudized with 0 coefficients</div><div><br></div><div><br></div><div>     atomic species   valence    mass     pseudopotential</div><div>        Ba            10.00   137.33000     Ba( 1.00)</div><div>        Mg             2.00    24.30500     Mg( 1.00)</div><div>        F              7.00    18.99800      F( 1.00)</div><div><br></div><div>      4 Sym. Ops. (no inversion) found ( 2 have fractional translation)</div><div><br></div><div><br></div><div><br></div><div>   Cartesian axes</div><div><br></div><div>     site n.     atom                  positions (alat units)</div><div>         1           Ba  tau(   1) = (  -0.0200000   3.0884290   0.6141993  )</div><div>         2           Ba  tau(   2) = (   0.0200000   0.4615710   1.3291993  )</div><div>         3           Ba  tau(   3) = (  -0.0200000   0.4615710   1.3291993  )</div><div>         4           Ba  tau(   4) = (   0.0200000   3.0884290   0.6141993  )</div><div>         5           Mg  tau(   5) = (   0.0200000   1.5841520   1.4276548  )</div><div>         6           Mg  tau(   6) = (  -0.0200000   1.9658480   0.7126548  )</div><div>         7           Mg  tau(   7) = (   0.0200000   1.9658480   0.7126548  )</div><div>         8           Mg  tau(   8) = (  -0.0200000   1.5841520   1.4276548  )</div><div>         9           F   tau(   9) = (  -0.0200000   1.2709000   1.0210200  )</div><div>        10           F   tau(  10) = (   0.0200000   2.2791000   0.3060200  )</div><div>        11           F   tau(  11) = (  -0.0200000   2.2791000   0.3060200  )</div><div>        12           F   tau(  12) = (   0.0200000   1.2709000   1.0210200  )</div><div>        13           F   tau(  13) = (   0.0200000   1.1147000   0.2445300  )</div><div>        14           F   tau(  14) = (  -0.0200000   2.4353000   0.9595300  )</div><div>        15           F   tau(  15) = (   0.0200000   2.4353000   0.9595300  )</div><div>        16           F   tau(  16) = (  -0.0200000   1.1147000   0.2445300  )</div><div>        17           F   tau(  17) = (  -0.0200000   1.7998500   1.1969100  )</div><div>        18           F   tau(  18) = (   0.0200000   1.7501500   0.4819100  )</div><div>        19           F   tau(  19) = (  -0.0200000   1.7501500   0.4819100  )</div><div>        20           F   tau(  20) = (   0.0200000   1.7998500   1.1969100  )</div><div>        21           F   tau(  21) = (   0.0300000   3.3441000   0.0357500  )</div><div>        22           F   tau(  22) = (  -0.0300000   0.2059000   0.7507500  )</div><div>        23           F   tau(  23) = (   0.0300000   0.2059000   0.7507500  )</div><div>        24           F   tau(  24) = (  -0.0300000   3.3441000   0.0357500  )</div><div><br></div><div>     number of k points=   258  Methfessel-Paxton smearing, width (Ry)=  0.0500</div><div><br></div><div>     Number of k-points >= 100: set verbosity='high' to print them.</div><div><br></div><div>     Dense  grid:   655723 G-vectors     FFT dimensions: ( 150, 150, 120)</div><div><br></div><div>     Smooth grid:   165831 G-vectors     FFT dimensions: (  96,  96,  72)</div><div><br></div><div>     Largest allocated arrays     est. size (Mb)     dimensions</div><div>        Kohn-Sham Wavefunctions         1.32 Mb     (     899,   96)</div><div>        NL pseudopotentials             3.18 Mb     (     899,  232)</div><div>        Each V/rho on FFT grid          1.72 Mb     (  112500)</div><div>        Each G-vector array             0.21 Mb     (   27323)</div><div>        G-vector shells                 0.10 Mb     (   13366)</div><div>     Largest temporary arrays     est. size (Mb)     dimensions</div><div>        Auxiliary wavefunctions         5.27 Mb     (     899,  384)</div><div>        Each subspace H/S matrix        2.25 Mb     (     384,  384)</div><div>        Each <psi_i|beta_j> matrix      0.34 Mb     (     232,   96)</div><div>        Arrays for rho mixing          13.73 Mb     (  112500,    8)</div><div><br></div><div>     Check: negative/imaginary core charge=   -0.000002    0.000000</div><div><br></div><div>     Initial potential from superposition of free atoms</div><div><br></div><div>     starting charge  159.92828, renormalised to  160.00000</div><div>     Starting wfc are  112 randomized atomic wfcs</div><div>     Checking if some PAW data can be deallocated...</div><div><br></div><div>     total cpu time spent up to now is      141.6 secs</div><div><br></div><div>     Self-consistent Calculation</div><div><br></div><div>     iteration #  1     ecut=   100.00 Ry     beta=0.70</div><div>     Davidson diagonalization with overlap</div><div><br></div><div> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div><div>     Error in routine cdiaghg (337):</div><div>     S matrix not positive definite</div><div> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div><div><br></div><div>     stopping ...</div></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>Input file :</div><div><br></div><div><div>&control</div><div>    calculation = 'relax',</div><div>    prefix='BaMgF4',</div><div>    pseudo_dir ='/home2/yogesh/RRCAT/BaMgF4/relax/',</div><div>    outdir='/home2/yogesh/RRCAT/BaMgF4/relax/'</div><div>    wfcdir='/home2/yogesh/RRCAT/BaMgF4/relax/'</div><div>    tstress = .true. ,</div><div>    tprnfor = .true. ,</div><div> /</div><div> &system</div><div>    ibrav=  9, celldm(1) = 7.85, celldm(2) = 3.55 ,  celldm(3) = 1.43 , nat=  6, ntyp= 3,</div><div>    ecutwfc = 100.0,</div><div>    ecutrho = 1000.0 ,</div><div>    occupations = 'smearing' ,</div><div>    degauss = 0.05 ,</div><div>    smearing = 'methfessel-paxton' ,</div><div>    space_group = 36</div><div> /</div><div> &electrons</div><div>    mixing_mode = 'plain'</div><div>    mixing_beta = 0.7</div><div>    conv_thr =  1.0d-8</div><div> /</div><div> &IONS</div><div>                ion_dynamics = 'bfgs' ,</div><div> /</div><div>ATOMIC_SPECIES</div><div>    Ba  137.33 Ba.pbe-spn-kjpaw_psl.0.2.3.upf</div><div>    Mg  24.305 Mg.pbe-n-kjpaw_psl.0.3.0.upf</div><div>    F   18.998 F.pbe-n-kjpaw_psl.0.1.upf</div><div>ATOMIC_POSITIONS crystal_sg</div><div>    Ba 0.48000    0.36998    0.42951</div><div>    Mg 0.02000    0.44624    -0.00164</div><div>    F  -0.0200    0.35800    0.71400</div><div>    F  0.02000    0.31400    0.17100</div><div>    F  -0.0200    0.50700    0.83700</div><div>    F  0.53000    0.44200    0.02500</div><div>K_POINTS automatic</div><div> 14 6 10 0 0 0</div></div><div><br></div><div><br></div><div><br></div><div>Thanks and Regards</div><div>Yogesh Kumar</div></div>