<div dir="ltr"><div>Dear QE community,</div><div><br></div><div>I am running a 2x2x2 supercell of Cu2O but  show this after some run<br></div><div>"Error in routine  cegterg (5014):       cannot allocate psi" after.</div><div><br></div><div>Can anybody put me through on what to do? My input and output files are as below #************************************************************************************<br>#*                  Generated by cif2cell 1.2.10 2019-05-21 19:56                   *<br>#*  T. Bjorkman, Comp. Phys. Commun. 182, 1183-1186 (2011). Please cite generously. *<br>#*                                                                                  *<br>#*                Data obtained from COD. Reference number : 1000063                *<br>#*                            Cu2 O   (Copper(I) oxide)                             *<br>#*            Neuburger, M C, Zeitschrift fuer Physik 67, 845-850 (1930)            *<br>#************************************************************************************<br>&CONTROL<br> calculation='scf'<br> outdir='.',<br> prefix='basic'<br> pseudo_dir='.'<br> verbosity='high'<br> tprnfor=.true.,<br> tstress=.true.,<br>/<br>&SYSTEM<br>  ibrav = 0<br>  A =    4.25200<br>  nat = 47<br>  ntyp = 2<br>  ecutwfc=90,<br>  ecutrho=450,<br>  input_dft='pbe',<br>  degauss=0.005d0,<br>  occupations='smearing'<br>  smearing='mv',<br>/<br> &ELECTRONS<br>   conv_thr=1d-08,<br>   mixing_beta=0.7d0,<br>/<br>CELL_PARAMETERS {alat}<br>  2.000000000000000   0.000000000000000   0.000000000000000 <br>  0.000000000000000   2.000000000000000   0.000000000000000 <br>  0.000000000000000   0.000000000000000   2.000000000000000 <br>ATOMIC_SPECIES<br>  O    15.99900   O.pbe-n-kjpaw_psl.0.1.UPF <br>  Cu   63.54600  Cu.pbe-dn-kjpaw_psl.0.2.UPF<br>ATOMIC_POSITIONS {crystal}<br>Cu   0.125000000000000   0.125000000000000   0.125000000000000 <br>Cu   0.375000000000000   0.375000000000000   0.125000000000000 <br>Cu   0.375000000000000   0.125000000000000   0.375000000000000 <br>Cu   0.125000000000000   0.375000000000000   0.375000000000000 <br>Cu   0.625000000000000   0.625000000000000   0.125000000000000 <br>Cu   0.625000000000000   0.125000000000000   0.125000000000000 <br>Cu   0.125000000000000   0.625000000000000   0.625000000000000 <br>Cu   0.125000000000000   0.125000000000000   0.625000000000000 <br>Cu   0.625000000000000   0.625000000000000   0.625000000000000 <br>Cu   0.625000000000000   0.125000000000000   0.625000000000000 <br>Cu   0.125000000000000   0.625000000000000   0.125000000000000 <br>Cu   0.875000000000000   0.875000000000000   0.125000000000000 <br>Cu   0.875000000000000   0.375000000000000   0.125000000000000 <br>Cu   0.375000000000000   0.875000000000000   0.625000000000000 <br>Cu   0.375000000000000   0.375000000000000   0.625000000000000 <br>Cu   0.875000000000000   0.875000000000000   0.625000000000000 <br>Cu   0.875000000000000   0.375000000000000   0.625000000000000 <br>Cu   0.375000000000000   0.875000000000000   0.125000000000000 <br>Cu   0.875000000000000   0.625000000000000   0.375000000000000 <br>Cu   0.875000000000000   0.125000000000000   0.375000000000000 <br>Cu   0.375000000000000   0.625000000000000   0.875000000000000 <br>Cu   0.375000000000000   0.125000000000000   0.875000000000000 <br>Cu   0.875000000000000   0.625000000000000   0.875000000000000 <br>Cu   0.875000000000000   0.125000000000000   0.875000000000000 <br>Cu   0.375000000000000   0.625000000000000   0.375000000000000 <br>Cu   0.625000000000000   0.875000000000000   0.375000000000000 <br>Cu   0.625000000000000   0.375000000000000   0.375000000000000 <br>Cu   0.125000000000000   0.875000000000000   0.875000000000000 <br>Cu   0.125000000000000   0.375000000000000   0.875000000000000 <br>Cu   0.625000000000000   0.875000000000000   0.875000000000000 <br>Cu   0.625000000000000   0.375000000000000   0.875000000000000 <br>Cu   0.125000000000000   0.875000000000000   0.375000000000000 <br> O   0.000000000000000   0.000000000000000   0.000000000000000 <br> O   0.250000000000000   0.250000000000000   0.250000000000000 <br> O   0.500000000000000   0.500000000000000   0.000000000000000 <br> O   0.500000000000000   0.000000000000000   0.000000000000000 <br> O   0.000000000000000   0.500000000000000   0.500000000000000 <br> O   0.000000000000000   0.000000000000000   0.500000000000000 <br>  <br> O   0.500000000000000   0.000000000000000   0.500000000000000 <br> O   0.000000000000000   0.500000000000000   0.000000000000000 <br> O   0.750000000000000   0.750000000000000   0.250000000000000 <br> O   0.750000000000000   0.250000000000000   0.250000000000000 <br> O   0.250000000000000   0.750000000000000   0.750000000000000 <br> O   0.250000000000000   0.250000000000000   0.750000000000000 <br> O   0.750000000000000   0.750000000000000   0.750000000000000 <br> O   0.750000000000000   0.250000000000000   0.750000000000000 <br> O   0.250000000000000   0.750000000000000   0.250000000000000 <br>K_POINTS  {automatic}<br>3 3 3 0 0 0</div><div><br></div><div><br></div><div>out put</div><div>Program PWSCF v.5.1 starts on 21May2019 at 21:38: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>     Reading input from <a href="http://modify_cuprous-conventional-222reg.in">modify_cuprous-conventional-222reg.in</a><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>     IMPORTANT: XC functional enforced from input :<br>     Exchange-correlation      = PBE ( 1  4  3  4 0)<br>     Any further DFT definition will be discarded<br>     Please, verify this is what you really want<br><br>               file O.pbe-n-kjpaw_psl.0.1.UPF: wavefunction(s)  2P renormalized<br>               file Cu.pbe-dn-kjpaw_psl.0.2.UPF: wavefunction(s)  4S 3D renormalized<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        9265    7409   1941               668829   478785   64229<br><br><br><br>     bravais-lattice index     =            0<br>     lattice parameter (alat)  =       8.0351  a.u.<br>     unit-cell volume          =    4150.1745 (a.u.)^3<br>     number of atoms/cell      =           47<br>     number of atomic types    =            2<br>     number of electrons       =       442.00<br>     number of Kohn-Sham states=          265<br>     kinetic-energy cutoff     =      90.0000  Ry<br>     charge density cutoff     =     450.0000  Ry<br>     convergence threshold     =      1.0E-08<br>     mixing beta               =       0.7000<br>     number of iterations used =            8  plain     mixing<br>     Exchange-correlation      = PBE ( 1  4  3  4 0)<br><br>     celldm(1)=   8.035116  celldm(2)=   0.000000  celldm(3)=   0.000000<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) = (   2.000000   0.000000   0.000000 )  <br>               a(2) = (   0.000000   2.000000   0.000000 )  <br>               a(3) = (   0.000000   0.000000   2.000000 )  <br><br>     reciprocal axes: (cart. coord. in units 2 pi/alat)<br>               b(1) = (  0.500000  0.000000  0.000000 )  <br>               b(2) = (  0.000000  0.500000  0.000000 )  <br>               b(3) = (  0.000000  0.000000  0.500000 )  <br><br><br>     PseudoPot. # 1 for  O read from file:<br>     ./O.pbe-n-kjpaw_psl.0.1.UPF<br>     MD5 check sum: 833c495a5d0f9b75bf046115549bffff<br>     Pseudo is Projector augmented-wave + core cor, Zval =  6.0<br>     Generated using "atomic" code by A. Dal Corso  v.6.3<br>     Shape of augmentation charge: BESSEL<br>     Using radial grid of 1095 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 0 coefficients <br><br><br>     PseudoPot. # 2 for Cu read from file:<br>     ./Cu.pbe-dn-kjpaw_psl.0.2.UPF<br>     MD5 check sum: 315c06c7b579c7c336aac2006b8584c5<br>     Pseudo is Projector augmented-wave + core cor, Zval = 11.0<br>     Generated using "atomic" code by A. Dal Corso  v.5.0.2 svn rev. 9415<br>     Shape of augmentation charge: BESSEL<br>     Using radial grid of 1199 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 0 coefficients <br><br><br>     atomic species   valence    mass     pseudopotential<br>        O              6.00    15.99900      O( 1.00)<br>        Cu            11.00    63.54600     Cu( 1.00)<br><br>     24 Sym. Ops. (no inversion) found<br><br><br>                                    s                        frac. trans.<br><br>      isym =  1     identity                                     <br><br> cryst.   s( 1) = (     1          0          0      )<br>                  (     0          1          0      )<br>                  (     0          0          1      )<br><br> cart.    s( 1) = (  1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000  1.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000  1.0000000 )<br><br><br>      isym =  2     180 deg rotation - cart. axis [0,0,1]        <br><br> cryst.   s( 2) = (    -1          0          0      )<br>                  (     0         -1          0      )<br>                  (     0          0          1      )<br><br> cart.    s( 2) = ( -1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000 -1.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000  1.0000000 )<br><br><br>      isym =  3     180 deg rotation - cart. axis [0,1,0]        <br><br> cryst.   s( 3) = (    -1          0          0      )<br>                  (     0          1          0      )<br>                  (     0          0         -1      )<br><br> cart.    s( 3) = ( -1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000  1.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000 -1.0000000 )<br><br><br>      isym =  4     180 deg rotation - cart. axis [1,0,0]        <br><br> cryst.   s( 4) = (     1          0          0      )<br>                  (     0         -1          0      )<br>                  (     0          0         -1      )<br><br> cart.    s( 4) = (  1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000 -1.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000 -1.0000000 )<br><br><br>      isym =  5     120 deg rotation - cart. axis [-1,-1,-1]     <br><br> cryst.   s( 5) = (     0          0          1      )<br>                  (     1          0          0      )<br>                  (     0          1          0      )<br><br> cart.    s( 5) = (  0.0000000  1.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000  1.0000000 )<br>                  (  1.0000000  0.0000000  0.0000000 )<br><br><br>      isym =  6     120 deg rotation - cart. axis [-1,1,1]       <br><br> cryst.   s( 6) = (     0          0         -1      )<br>                  (    -1          0          0      )<br>                  (     0          1          0      )<br><br> cart.    s( 6) = (  0.0000000 -1.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000  1.0000000 )<br>                  ( -1.0000000  0.0000000  0.0000000 )<br><br><br>      isym =  7     120 deg rotation - cart. axis [1,1,-1]       <br><br> cryst.   s( 7) = (     0          0         -1      )<br>                  (     1          0          0      )<br>                  (     0         -1          0      )<br><br> cart.    s( 7) = (  0.0000000  1.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000 -1.0000000 )<br>                  ( -1.0000000  0.0000000  0.0000000 )<br><br><br>      isym =  8     120 deg rotation - cart. axis [1,-1,1]       <br><br> cryst.   s( 8) = (     0          0          1      )<br>                  (    -1          0          0      )<br>                  (     0         -1          0      )<br><br> cart.    s( 8) = (  0.0000000 -1.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000 -1.0000000 )<br>                  (  1.0000000  0.0000000  0.0000000 )<br><br><br>      isym =  9     120 deg rotation - cart. axis [1,1,1]        <br><br> cryst.   s( 9) = (     0          1          0      )<br>                  (     0          0          1      )<br>                  (     1          0          0      )<br><br> cart.    s( 9) = (  0.0000000  0.0000000  1.0000000 )<br>                  (  1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000  1.0000000  0.0000000 )<br><br><br>      isym = 10     120 deg rotation - cart. axis [-1,1,-1]      <br><br> cryst.   s(10) = (     0         -1          0      )<br>                  (     0          0         -1      )<br>                  (     1          0          0      )<br><br> cart.    s(10) = (  0.0000000  0.0000000  1.0000000 )<br>                  ( -1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000 -1.0000000  0.0000000 )<br><br><br>      isym = 11     120 deg rotation - cart. axis [1,-1,-1]      <br><br> cryst.   s(11) = (     0         -1          0      )<br>                  (     0          0          1      )<br>                  (    -1          0          0      )<br><br> cart.    s(11) = (  0.0000000  0.0000000 -1.0000000 )<br>                  ( -1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000  1.0000000  0.0000000 )<br><br><br>      isym = 12     120 deg rotation - cart. axis [-1,-1,1]      <br><br> cryst.   s(12) = (     0          1          0      )<br>                  (     0          0         -1      )<br>                  (    -1          0          0      )<br><br> cart.    s(12) = (  0.0000000  0.0000000 -1.0000000 )<br>                  (  1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000 -1.0000000  0.0000000 )<br><br><br>      isym = 13     inv. 180 deg rotation - cart. axis [1,1,0]   <br><br> cryst.   s(13) = (     0         -1          0      )<br>                  (    -1          0          0      )<br>                  (     0          0          1      )<br><br> cart.    s(13) = (  0.0000000 -1.0000000  0.0000000 )<br>                  ( -1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000  1.0000000 )<br><br><br>      isym = 14     inv. 180 deg rotation - cart. axis [1,-1,0]  <br><br> cryst.   s(14) = (     0          1          0      )<br>                  (     1          0          0      )<br>                  (     0          0          1      )<br><br> cart.    s(14) = (  0.0000000  1.0000000  0.0000000 )<br>                  (  1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000  1.0000000 )<br><br><br>      isym = 15     inv.  90 deg rotation - cart. axis [0,0,-1]  <br><br> cryst.   s(15) = (     0          1          0      )<br>                  (    -1          0          0      )<br>                  (     0          0         -1      )<br><br> cart.    s(15) = (  0.0000000 -1.0000000  0.0000000 )<br>                  (  1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000 -1.0000000 )<br><br><br>      isym = 16     inv.  90 deg rotation - cart. axis [0,0,1]   <br><br> cryst.   s(16) = (     0         -1          0      )<br>                  (     1          0          0      )<br>                  (     0          0         -1      )<br><br> cart.    s(16) = (  0.0000000  1.0000000  0.0000000 )<br>                  ( -1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000 -1.0000000 )<br><br><br>      isym = 17     inv. 180 deg rotation - cart. axis [1,0,1]   <br><br> cryst.   s(17) = (     0          0         -1      )<br>                  (     0          1          0      )<br>                  (    -1          0          0      )<br><br> cart.    s(17) = (  0.0000000  0.0000000 -1.0000000 )<br>                  (  0.0000000  1.0000000  0.0000000 )<br>                  ( -1.0000000  0.0000000  0.0000000 )<br><br><br>      isym = 18     inv. 180 deg rotation - cart. axis [-1,0,1]  <br><br> cryst.   s(18) = (     0          0          1      )<br>                  (     0          1          0      )<br>                  (     1          0          0      )<br><br> cart.    s(18) = (  0.0000000  0.0000000  1.0000000 )<br>                  (  0.0000000  1.0000000  0.0000000 )<br>                  (  1.0000000  0.0000000  0.0000000 )<br><br><br>      isym = 19     inv.  90 deg rotation - cart. axis [0,1,0]   <br><br> cryst.   s(19) = (     0          0          1      )<br>                  (     0         -1          0      )<br>                  (    -1          0          0      )<br><br> cart.    s(19) = (  0.0000000  0.0000000 -1.0000000 )<br>                  (  0.0000000 -1.0000000  0.0000000 )<br>                  (  1.0000000  0.0000000  0.0000000 )<br><br><br>      isym = 20     inv.  90 deg rotation - cart. axis [0,-1,0]  <br><br> cryst.   s(20) = (     0          0         -1      )<br>                  (     0         -1          0      )<br>                  (     1          0          0      )<br><br> cart.    s(20) = (  0.0000000  0.0000000  1.0000000 )<br>                  (  0.0000000 -1.0000000  0.0000000 )<br>                  ( -1.0000000  0.0000000  0.0000000 )<br><br><br>      isym = 21     inv. 180 deg rotation - cart. axis [0,1,1]   <br><br> cryst.   s(21) = (     1          0          0      )<br>                  (     0          0         -1      )<br>                  (     0         -1          0      )<br><br> cart.    s(21) = (  1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000 -1.0000000 )<br>                  (  0.0000000 -1.0000000  0.0000000 )<br><br><br>      isym = 22     inv. 180 deg rotation - cart. axis [0,1,-1]  <br><br> cryst.   s(22) = (     1          0          0      )<br>                  (     0          0          1      )<br>                  (     0          1          0      )<br><br> cart.    s(22) = (  1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000  1.0000000 )<br>                  (  0.0000000  1.0000000  0.0000000 )<br><br><br>      isym = 23     inv.  90 deg rotation - cart. axis [-1,0,0]  <br><br> cryst.   s(23) = (    -1          0          0      )<br>                  (     0          0          1      )<br>                  (     0         -1          0      )<br><br> cart.    s(23) = ( -1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000 -1.0000000 )<br>                  (  0.0000000  1.0000000  0.0000000 )<br><br><br>      isym = 24     inv.  90 deg rotation - cart. axis [1,0,0]   <br><br> cryst.   s(24) = (    -1          0          0      )<br>                  (     0          0         -1      )<br>                  (     0          1          0      )<br><br> cart.    s(24) = ( -1.0000000  0.0000000  0.0000000 )<br>                  (  0.0000000  0.0000000  1.0000000 )<br>                  (  0.0000000 -1.0000000  0.0000000 )<br><br><br>     point group T_d (-43m) <br>     there are  5 classes<br>     the character table:<br><br>       E     8C3   3C2   6S4   6s_d <br>A_1    1.00  1.00  1.00  1.00  1.00<br>A_2    1.00  1.00  1.00 -1.00 -1.00<br>E      2.00 -1.00  2.00  0.00  0.00<br>T_1    3.00  0.00 -1.00  1.00 -1.00<br>T_2    3.00  0.00 -1.00 -1.00  1.00<br><br>     the symmetry operations in each class:<br>     E        1<br>     3C2      2    3    4<br>     8C3      5    7    8    6   12    9   10   11<br>     6s_d    13   14   17   18   21   22<br>     6S4     15   16   20   19   23   24<br><br>   Cartesian axes<br><br>     site n.     atom                  positions (alat units)<br>         1           Cu  tau(   1) = (   0.2500000   0.2500000   0.2500000  )<br>         2           Cu  tau(   2) = (   0.7500000   0.7500000   0.2500000  )<br>         3           Cu  tau(   3) = (   0.7500000   0.2500000   0.7500000  )<br>         4           Cu  tau(   4) = (   0.2500000   0.7500000   0.7500000  )<br>         5           Cu  tau(   5) = (   1.2500000   1.2500000   0.2500000  )<br>         6           Cu  tau(   6) = (   1.2500000   0.2500000   0.2500000  )<br>         7           Cu  tau(   7) = (   0.2500000   1.2500000   1.2500000  )<br>         8           Cu  tau(   8) = (   0.2500000   0.2500000   1.2500000  )<br>         9           Cu  tau(   9) = (   1.2500000   1.2500000   1.2500000  )<br>        10           Cu  tau(  10) = (   1.2500000   0.2500000   1.2500000  )<br>        11           Cu  tau(  11) = (   0.2500000   1.2500000   0.2500000  )<br>        12           Cu  tau(  12) = (   1.7500000   1.7500000   0.2500000  )<br>        13           Cu  tau(  13) = (   1.7500000   0.7500000   0.2500000  )<br>        14           Cu  tau(  14) = (   0.7500000   1.7500000   1.2500000  )<br>        15           Cu  tau(  15) = (   0.7500000   0.7500000   1.2500000  )<br>        16           Cu  tau(  16) = (   1.7500000   1.7500000   1.2500000  )<br>        17           Cu  tau(  17) = (   1.7500000   0.7500000   1.2500000  )<br>        18           Cu  tau(  18) = (   0.7500000   1.7500000   0.2500000  )<br>        19           Cu  tau(  19) = (   1.7500000   1.2500000   0.7500000  )<br>        20           Cu  tau(  20) = (   1.7500000   0.2500000   0.7500000  )<br>        21           Cu  tau(  21) = (   0.7500000   1.2500000   1.7500000  )<br>        22           Cu  tau(  22) = (   0.7500000   0.2500000   1.7500000  )<br>        23           Cu  tau(  23) = (   1.7500000   1.2500000   1.7500000  )<br>        24           Cu  tau(  24) = (   1.7500000   0.2500000   1.7500000  )<br>        25           Cu  tau(  25) = (   0.7500000   1.2500000   0.7500000  )<br>        26           Cu  tau(  26) = (   1.2500000   1.7500000   0.7500000  )<br>        27           Cu  tau(  27) = (   1.2500000   0.7500000   0.7500000  )<br>        28           Cu  tau(  28) = (   0.2500000   1.7500000   1.7500000  )<br>        29           Cu  tau(  29) = (   0.2500000   0.7500000   1.7500000  )<br>        30           Cu  tau(  30) = (   1.2500000   1.7500000   1.7500000  )<br>        31           Cu  tau(  31) = (   1.2500000   0.7500000   1.7500000  )<br>        32           Cu  tau(  32) = (   0.2500000   1.7500000   0.7500000  )<br>        33           O   tau(  33) = (   0.0000000   0.0000000   0.0000000  )<br>        34           O   tau(  34) = (   0.5000000   0.5000000   0.5000000  )<br>        35           O   tau(  35) = (   1.0000000   1.0000000   0.0000000  )<br>        36           O   tau(  36) = (   1.0000000   0.0000000   0.0000000  )<br>        37           O   tau(  37) = (   0.0000000   1.0000000   1.0000000  )<br>        38           O   tau(  38) = (   0.0000000   0.0000000   1.0000000  )<br>        39           O   tau(  39) = (   1.0000000   0.0000000   1.0000000  )<br>        40           O   tau(  40) = (   0.0000000   1.0000000   0.0000000  )<br>        41           O   tau(  41) = (   1.5000000   1.5000000   0.5000000  )<br>        42           O   tau(  42) = (   1.5000000   0.5000000   0.5000000  )<br>        43           O   tau(  43) = (   0.5000000   1.5000000   1.5000000  )<br>        44           O   tau(  44) = (   0.5000000   0.5000000   1.5000000  )<br>        45           O   tau(  45) = (   1.5000000   1.5000000   1.5000000  )<br>        46           O   tau(  46) = (   1.5000000   0.5000000   1.5000000  )<br>        47           O   tau(  47) = (   0.5000000   1.5000000   0.5000000  )<br><br>   Crystallographic axes<br><br>     site n.     atom                  positions (cryst. coord.)<br>         1           Cu  tau(   1) = (  0.1250000  0.1250000  0.1250000  )<br>         2           Cu  tau(   2) = (  0.3750000  0.3750000  0.1250000  )<br>         3           Cu  tau(   3) = (  0.3750000  0.1250000  0.3750000  )<br>         4           Cu  tau(   4) = (  0.1250000  0.3750000  0.3750000  )<br>         5           Cu  tau(   5) = (  0.6250000  0.6250000  0.1250000  )<br>         6           Cu  tau(   6) = (  0.6250000  0.1250000  0.1250000  )<br>         7           Cu  tau(   7) = (  0.1250000  0.6250000  0.6250000  )<br>         8           Cu  tau(   8) = (  0.1250000  0.1250000  0.6250000  )<br>         9           Cu  tau(   9) = (  0.6250000  0.6250000  0.6250000  )<br>        10           Cu  tau(  10) = (  0.6250000  0.1250000  0.6250000  )<br>        11           Cu  tau(  11) = (  0.1250000  0.6250000  0.1250000  )<br>        12           Cu  tau(  12) = (  0.8750000  0.8750000  0.1250000  )<br>        13           Cu  tau(  13) = (  0.8750000  0.3750000  0.1250000  )<br>        14           Cu  tau(  14) = (  0.3750000  0.8750000  0.6250000  )<br>        15           Cu  tau(  15) = (  0.3750000  0.3750000  0.6250000  )<br>        16           Cu  tau(  16) = (  0.8750000  0.8750000  0.6250000  )<br>        17           Cu  tau(  17) = (  0.8750000  0.3750000  0.6250000  )<br>        18           Cu  tau(  18) = (  0.3750000  0.8750000  0.1250000  )<br>        19           Cu  tau(  19) = (  0.8750000  0.6250000  0.3750000  )<br>        20           Cu  tau(  20) = (  0.8750000  0.1250000  0.3750000  )<br>        21           Cu  tau(  21) = (  0.3750000  0.6250000  0.8750000  )<br>        22           Cu  tau(  22) = (  0.3750000  0.1250000  0.8750000  )<br>        23           Cu  tau(  23) = (  0.8750000  0.6250000  0.8750000  )<br>        24           Cu  tau(  24) = (  0.8750000  0.1250000  0.8750000  )<br>        25           Cu  tau(  25) = (  0.3750000  0.6250000  0.3750000  )<br>        26           Cu  tau(  26) = (  0.6250000  0.8750000  0.3750000  )<br>        27           Cu  tau(  27) = (  0.6250000  0.3750000  0.3750000  )<br>        28           Cu  tau(  28) = (  0.1250000  0.8750000  0.8750000  )<br>        29           Cu  tau(  29) = (  0.1250000  0.3750000  0.8750000  )<br>        30           Cu  tau(  30) = (  0.6250000  0.8750000  0.8750000  )<br>        31           Cu  tau(  31) = (  0.6250000  0.3750000  0.8750000  )<br>        32           Cu  tau(  32) = (  0.1250000  0.8750000  0.3750000  )<br>        33           O   tau(  33) = (  0.0000000  0.0000000  0.0000000  )<br>        34           O   tau(  34) = (  0.2500000  0.2500000  0.2500000  )<br>        35           O   tau(  35) = (  0.5000000  0.5000000  0.0000000  )<br>        36           O   tau(  36) = (  0.5000000  0.0000000  0.0000000  )<br>        37           O   tau(  37) = (  0.0000000  0.5000000  0.5000000  )<br>        38           O   tau(  38) = (  0.0000000  0.0000000  0.5000000  )<br>        39           O   tau(  39) = (  0.5000000  0.0000000  0.5000000  )<br>        40           O   tau(  40) = (  0.0000000  0.5000000  0.0000000  )<br>        41           O   tau(  41) = (  0.7500000  0.7500000  0.2500000  )<br>        42           O   tau(  42) = (  0.7500000  0.2500000  0.2500000  )<br>        43           O   tau(  43) = (  0.2500000  0.7500000  0.7500000  )<br>        44           O   tau(  44) = (  0.2500000  0.2500000  0.7500000  )<br>        45           O   tau(  45) = (  0.7500000  0.7500000  0.7500000  )<br>        46           O   tau(  46) = (  0.7500000  0.2500000  0.7500000  )<br>        47           O   tau(  47) = (  0.2500000  0.7500000  0.2500000  )<br><br>     number of k points=     4  Marzari-Vanderbilt smearing, width (Ry)=  0.0050<br>                       cart. coord. in units 2pi/alat<br>        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   0.0740741<br>        k(    2) = (   0.0000000   0.0000000   0.1666667), wk =   0.4444444<br>        k(    3) = (   0.0000000   0.1666667   0.1666667), wk =   0.8888889<br>        k(    4) = (   0.1666667   0.1666667   0.1666667), wk =   0.5925926<br><br>                       cryst. coord.<br>        k(    1) = (   0.0000000   0.0000000   0.0000000), wk =   0.0740741<br>        k(    2) = (   0.0000000   0.0000000   0.3333333), wk =   0.4444444<br>        k(    3) = (   0.0000000   0.3333333   0.3333333), wk =   0.8888889<br>        k(    4) = (   0.3333333   0.3333333   0.3333333), wk =   0.5925926<br><br>     Dense  grid:   668829 G-vectors     FFT dimensions: ( 120, 120, 120)<br><br>     Smooth grid:   478785 G-vectors     FFT dimensions: ( 100, 100, 100)<br><br>     Largest allocated arrays     est. size (Mb)     dimensions<br>        Kohn-Sham Wavefunctions       242.01 Mb     (   59850,  265)<br>        NL pseudopotentials           635.61 Mb     (   59850,  696)<br>        Each V/rho on FFT grid         26.37 Mb     ( 1728000)<br>        Each G-vector array             5.10 Mb     (  668829)<br>        G-vector shells                 0.02 Mb     (    2455)<br>     Largest temporary arrays     est. size (Mb)     dimensions<br>        Auxiliary wavefunctions       968.03 Mb     (   59850, 1060)<br>        Each subspace H/S matrix       17.14 Mb     (    1060, 1060)<br>        Each <psi_i|beta_j> matrix      2.81 Mb     (     696,  265)<br>        Arrays for rho mixing         210.94 Mb     ( 1728000,    8)<br><br>     Check: negative/imaginary core charge=   -0.000025    0.000000<br><br>     Initial potential from superposition of free atoms<br><br>     starting charge  441.99268, renormalised to  442.00000<br>     Starting wfc are  348 randomized atomic wfcs<br>     Checking if some PAW data can be deallocated... <br><br>     total cpu time spent up to now is     1375.2 secs<br><br>     Self-consistent Calculation<br><br>     iteration #  1     ecut=    90.00 Ry     beta=0.70<br>     Davidson diagonalization with overlap<br><br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>     Error in routine  cegterg (5014):<br>      cannot allocate psi<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br><br>     stopping ...<br>--------------------------------------------------------------------------<br>MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD <br>with errorcode 1.<br><br>NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.<br>You may or may not see output from other processes, depending on<br>exactly when Open MPI kills them.<br>--------------------------------------------------------------------------</div><div>thanks</div><div>Ismail Ibrahim Idowu</div><div>Ahmadu Bello University Zaria, Nigeria<br></div></div>