<div dir="ltr"><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Hello All,<br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">My SCF calculation was fine and it converged in 23 iterations. But mu nscf calculations (band structure), is showing the following error. Please have a look.<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Here is my input file,<br><br> &CONTROL<br>                 calculation = 'nscf' ,<br>                restart_mode = 'from_scratch' ,<br>                      outdir = '/home/manu/espresso-5.1/bin/Ga2O3_bands/' ,<br>                  pseudo_dir = '/home/manu/espresso-5.1/pseudo/' ,<br>                      prefix = 'gallium' ,<br>                     disk_io = 'default' ,<br>                   verbosity = 'default' ,<br> /<br> &SYSTEM<br>                       ibrav = 13,<br>                           A = 12.208 ,<br>                           B = 3.031 ,<br>                           C = 5.751 ,<br>                       cosAB = 0 ,<br>                       cosAC = -0.23 ,<br>                       cosBC = 0 ,<br>                         nat = 20,<br>                        ntyp = 2,<br>                     ecutwfc = 70 ,<br>                     ecutrho = 800 ,<br>                        nbnd = 90,<br>                 occupations = 'smearing' ,<br>                     degauss = 0.001 ,<br>                    smearing = 'gaussian' ,<br> /<br> &ELECTRONS<br>            electron_maxstep = 200,<br>                    conv_thr = 5.D-12 ,<br>                 startingpot = 'file' ,<br>                 startingwfc = 'atomic' ,<br>                 mixing_mode = 'TF' ,<br>                 mixing_beta = 0.5 ,<br>                 mixing_ndim = 10,<br>             diagonalization = 'david' ,<br> /<br>ATOMIC_SPECIES<br>   Ga   69.72300  Ga.pbe-n-van.UPF <br>    O   15.99400  O.pbe-van_ak.UPF <br>ATOMIC_POSITIONS crystal <br>   Ga      0.090500000    0.000000000    0.794600000    1  1  1 <br>   Ga      0.909500000    0.000000000    0.205400000    1  1  1 <br>   Ga      0.590500000    0.500000000    0.794600000    1  1  1 <br>   Ga      0.409500000    0.500000000    0.205400000    1  1  1 <br>   Ga      0.158600000    0.500000000    0.314000000    1  1  1 <br>   Ga      0.841400000    0.500000000    0.686000000    1  1  1 <br>   Ga      0.658600000    0.000000000    0.314000000    1  1  1 <br>   Ga      0.341400000    0.000000000    0.686000000    1  1  1 <br>    O      0.164500000    0.000000000    0.109800000    1  1  1 <br>    O      0.835500000    0.000000000    0.890200000    1  1  1 <br>    O      0.664500000    0.500000000    0.109800000    1  1  1 <br>    O      0.335500000    0.500000000    0.890200000    1  1  1 <br>    O      0.173300000    0.000000000    0.563200000    1  1  1 <br>    O      0.826700000    0.000000000    0.436800000    1  1  1 <br>    O      0.673300000    0.500000000    0.563200000    1  1  1 <br>    O      0.326700000    0.500000000    0.436800000    1  1  1 <br>    O      0.996000000    0.500000000    0.256600000    1  1  1 <br>    O      0.004000000    0.500000000    0.743400000    1  1  1 <br>    O      0.496000000    0.000000000    0.256600000    1  1  1 <br>    O      0.504000000    0.000000000    0.747340000    1  1  1 <br>K_POINTS crystal <br>6 <br>   0.000000000    0.500000000    0.000000000      2.000000000 <br>   0.500000000    0.500000000    0.000000000      2.000000000 <br>   0.000000000    0.000000000    0.000000000      1.000000000 <br>   0.000000000    0.000000000    0.500000000      1.000000000 <br>   0.500000000    0.500000000    0.500000000      2.000000000 <br>   0.000000000    0.500000000    0.500000000      2.000000000 <br><br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Here is my output file,<br><br><br>     Program PWSCF v.5.1 starts on 26Nov2014 at 13:36:59 <br><br>     This program is part of the open-source Quantum ESPRESSO suite<br>     for quantum simulation of materials; please cite<br>         "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);<br>          URL <a href="http://www.quantum-espresso.org">http://www.quantum-espresso.org</a>", <br>     in publications or presentations arising from this work. More details at<br>     <a href="http://www.quantum-espresso.org/quote">http://www.quantum-espresso.org/quote</a><br><br>     Parallel version (MPI), running on     1 processors<br>     Waiting for input...<br>     Reading input from standard input<br><br>     Current dimensions of program PWSCF are:<br>     Max number of different atomic species (ntypx) = 10<br>     Max number of k-points (npk) =  40000<br>     Max angular momentum in pseudopotentials (lmaxx) =  3<br><br>     Atomic positions and unit cell read from directory:<br>     /home/manu/espresso-5.1/bin/Ga2O3_bands/gallium.save/<br><br><br>     Subspace diagonalization in iterative solution of the eigenvalue problem:<br>     a serial algorithm will be used<br><br><br>     G-vector sticks info<br>     --------------------<br>     sticks:   dense  smooth     PW     G-vecs:    dense   smooth      PW<br>     Sum        4645    1629    477               274229    56803    8851<br><br><br><br>     bravais-lattice index     =           13<br>     lattice parameter (alat)  =      23.0698  a.u.<br>     unit-cell volume          =     718.0264 (a.u.)^3<br>     number of atoms/cell      =           20<br>     number of atomic types    =            2<br>     number of electrons       =       176.00<br>     number of Kohn-Sham states=           90<br>     kinetic-energy cutoff     =      70.0000  Ry<br>     charge density cutoff     =     800.0000  Ry<br>     Exchange-correlation      =  SLA  PW   PBE  PBE ( 1  4  3  4 0)<br><br>     celldm(1)=  23.069777  celldm(2)=   0.248280  celldm(3)=   0.471085<br>     celldm(4)=   0.000000  celldm(5)=   0.000000  celldm(6)=   0.000000<br><br>     crystal axes: (cart. coord. in units of alat)<br>               a(1) = (   0.500000   0.000000  -0.235542 )  <br>               a(2) = (   0.000000   0.248280   0.000000 )  <br>               a(3) = (   0.500000   0.000000   0.235542 )  <br><br>     reciprocal axes: (cart. coord. in units 2 pi/alat)<br>               b(1) = (  1.000000  0.000000 -2.122761 )  <br>               b(2) = ( -0.000000  4.027714  0.000000 )  <br>               b(3) = (  1.000000 -0.000000  2.122761 )  <br><br><br>     PseudoPot. # 1 for Ga read from file:<br>     /home/manu/espresso-5.1/pseudo/Ga.pbe-n-van.UPF<br>     MD5 check sum: c9a755f92899fce861e122f249cbad98<br>     Pseudo is Ultrasoft + core correction, Zval = 13.0<br>     Generated by new atomic code, or converted to UPF format<br>     Using radial grid of  871 points,  6 beta functions with: <br>                l(1) =   0<br>                l(2) =   0<br>                l(3) =   1<br>                l(4) =   1<br>                l(5) =   2<br>                l(6) =   2<br>     Q(r) pseudized with  8 coefficients,  rinner =    1.000   1.000   1.000<br>                                                       1.000   1.000<br><br>     PseudoPot. # 2 for O  read from file:<br>     /home/manu/espresso-5.1/pseudo/O.pbe-van_ak.UPF<br>     MD5 check sum: b3bd5b423a8218f475a9b6107bb6a37c<br>     Pseudo is Ultrasoft, Zval =  6.0<br>     Generated by new atomic code, or converted to UPF format<br>     Using radial grid of  737 points,  4 beta functions with: <br>                l(1) =   0<br>                l(2) =   0<br>                l(3) =   1<br>                l(4) =   1<br>     Q(r) pseudized with  8 coefficients,  rinner =    0.800   0.800   0.800<br><br><br>     atomic species   valence    mass     pseudopotential<br>        Ga            13.00    69.72300     Ga( 1.00)<br>        O              6.00    15.99400     O ( 1.00)<br><br>      2 Sym. Ops. (no inversion) found<br><br><br><br>   Cartesian axes<br><br>     site n.     atom                  positions (alat units)<br>         1           Ga  tau(   1) = (   0.4425500   0.0000000   0.1658453  )<br>         2           Ga  tau(   2) = (   0.5574500   0.0000000  -0.1658453  )<br>         3           Ga  tau(   3) = (   0.6925500   0.1241399   0.0480742  )<br>         4           Ga  tau(   4) = (   0.3074500   0.1241399  -0.0480742  )<br>         5           Ga  tau(   5) = (   0.2363000   0.1241399   0.0366033  )<br>         6           Ga  tau(   6) = (   0.7637000   0.1241399  -0.0366033  )<br>         7           Ga  tau(   7) = (   0.4863000   0.0000000  -0.0811679  )<br>         8           Ga  tau(   8) = (   0.5137000   0.0000000   0.0811679  )<br>         9           O   tau(   9) = (   0.1371500   0.0000000  -0.0128842  )<br>        10           O   tau(  10) = (   0.8628500   0.0000000   0.0128842  )<br>        11           O   tau(  11) = (   0.3871500   0.1241399  -0.1306553  )<br>        12           O   tau(  12) = (   0.6128500   0.1241399   0.1306553  )<br>        13           O   tau(  13) = (   0.3682500   0.0000000   0.0918379  )<br>        14           O   tau(  14) = (   0.6317500   0.0000000  -0.0918379  )<br>        15           O   tau(  15) = (   0.6182500   0.1241399  -0.0259332  )<br>        16           O   tau(  16) = (   0.3817500   0.1241399   0.0259332  )<br>        17           O   tau(  17) = (   0.6263000   0.1241399  -0.1741600  )<br>        18           O   tau(  18) = (   0.3737000   0.1241399   0.1741600  )<br>        19           O   tau(  19) = (   0.3763000   0.0000000  -0.0563888  )<br>        20           O   tau(  20) = (   0.6256700   0.0000000   0.0573169  )<br><br>     number of k points=     9  gaussian smearing, width (Ry)=  0.0010<br>                       cart. coord. in units 2pi/alat<br>        k(    1) = (   0.0000000   2.0138568   0.0000000), wk =   0.4000000<br>        k(    2) = (   0.5000000   2.0138568  -1.0613806), wk =   0.2000000<br>        k(    3) = (   0.0000000   0.0000000   0.0000000), wk =   0.2000000<br>        k(    4) = (   0.5000000   0.0000000   1.0613806), wk =   0.1000000<br>        k(    5) = (   1.0000000   2.0138568   0.0000000), wk =   0.4000000<br>        k(    6) = (   0.5000000   2.0138568   1.0613806), wk =   0.2000000<br>        k(    7) = (   0.5000000  -2.0138568   1.0613806), wk =   0.2000000<br>        k(    8) = (   0.5000000   0.0000000  -1.0613806), wk =   0.1000000<br>        k(    9) = (   0.5000000  -2.0138568  -1.0613806), wk =   0.2000000<br><br>     Dense  grid:   274229 G-vectors     FFT dimensions: ( 120,  54, 120)<br><br>     Smooth grid:    56803 G-vectors     FFT dimensions: (  72,  32,  72)<br><br>     Largest allocated arrays     est. size (Mb)     dimensions<br>        Kohn-Sham Wavefunctions         9.77 Mb     (    7116,   90)<br>        NL pseudopotentials            26.06 Mb     (    7116,  240)<br>        Each V/rho on FFT grid         11.87 Mb     (  777600)<br>        Each G-vector array             2.09 Mb     (  274229)<br>        G-vector shells                 0.28 Mb     (   36072)<br>     Largest temporary arrays     est. size (Mb)     dimensions<br>        Auxiliary wavefunctions        39.09 Mb     (    7116,  360)<br>        Each subspace H/S matrix        1.98 Mb     (     360,  360)<br>        Each <psi_i|beta_j> matrix      0.33 Mb     (     240,   90)<br><br>     The potential is recalculated from file :<br>     /home/manu/espresso-5.1/bin/Ga2O3_bands/gallium.save/charge-density.dat<br><br>     Starting wfc are  120 atomic wfcs<br><br>     Band Structure Calculation<br>     Davidson diagonalization with overlap<br><br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>     Error in routine cdiaghg (497):<br>     S matrix not positive definite<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br><br>     stopping ...<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Thanks,<br><br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif">Manu<br></div><div class="gmail_default" style="font-family:arial,helvetica,sans-serif"><br></div></div>