<div dir="ltr"><div>Dear QE users and experts,</div><div>I was doing band structure computation of a solid solution CdS and CdTe compounds using VCA and supercell respectively with the above mentioned hybrid functional. The computation went well with the steps 0.25 and 0.50 for both VCA and supercell approaches. However, at 0.75 step of both VCA and supercell, the computation aborted with no error message. <br></div><div>I will appreciate your useful inputs in this regard.</div><div>My input file and some portion of the output file are pasted here for your comments.</div><div>Thank you.</div><div><br></div><div>&CONTROL<br>  calculation = 'scf'<br>  restart_mode = 'from_scratch'<br>  etot_conv_thr =   8.0000000000d-05<br>  forc_conv_thr =   1.0000000000d-04<br>  outdir = './CdS.out/'<br>  prefix = 'CdS'<br>  pseudo_dir = './pseudo/'<br>  tprnfor = .true.<br>  tstress = .true.<br>  verbosity = 'high'<br>/<br>&SYSTEM<br>  degauss =   7.3498618000d-03<br>  ecutrho =   2.4000000000d+02<br>  ecutwfc =   6.0000000000d+01<br>  ibrav = 0<br>  nat = 8<br>  ntyp = 3<br>  nbnd = 40<br>  occupations = 'smearing'<br>  smearing = 'cold'<br>  input_dft='pbe0', nqx1=3, nqx2=3, nqx3=3, <br>  exxdiv_treatment='gygi-baldereschi'<br>  x_gamma_extrapolation = .true.<br>/<br>&ELECTRONS<br>  conv_thr =   1.6000000000d-09<br>  electron_maxstep = 80<br>  mixing_beta =   4.0000000000d-01<br>/<br>&IONS<br>  ion_dynamics = 'bfgs'<br>/<br>ATOMIC_SPECIES<br>Cd     112.41 Cd.pbe-n-nc.UPF <br>S      32.066 S.pbe-n-nc.UPF<br>Te     127.6 Te.pbe-n-nc.UPF<br>ATOMIC_POSITIONS (crystal)<br>Cd      -0.014782148   0.014782148  -0.014782148<br>Cd      -0.014782148   0.485217852   0.514782148<br>Cd       0.514782148   0.014782148   0.514782148<br>Cd       0.514782148   0.485217852  -0.014782148<br>S        0.750000000   0.250000000   0.750000000<br>Te       0.750000000   0.750000000   0.250000000<br>Te       0.250000000   0.250000000   0.250000000<br>Te       0.250000000   0.750000000   0.750000000<br>K_POINTS automatic<br>6 6 6 0 0 0<br>CELL_PARAMETERS angstrom<br>      6.45800             0.0000000000       0.0000000000<br>      0.0000000000       6.45800            0.0000000000<br>      0.0000000000       0.0000000000       6.45800</div><div><br></div><div><br></div><div>    Program PWSCF v.6.4.1 starts on  5Aug2020 at 14:55:44 <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>         "P. Giannozzi et al., J. Phys.:Condens. Matter 29 465901 (2017);<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 & OpenMP), running on      40 processor cores<br>     Number of MPI processes:                40<br>     Threads/MPI process:                     1<br><br>     MPI processes distributed on     1 nodes<br>     R & G space division:  proc/nbgrp/npool/nimage =      40<br>     Reading input from <a href="http://CdTe0.75S.scf.in">CdTe0.75S.scf.in</a><br>Warning: card &IONS ignored<br>Warning: card   ION_DYNAMICS = 'BFGS' ignored<br>Warning: card / ignored<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>     Message from routine read_upf::<br>     Pseudo file Cd.pbe-n-nc.UPF has been fixed on the fly.<br>     IMPORTANT: XC functional enforced from input :<br>     Exchange-correlation      = PBE0 ( 6  4  8  4 0 0)<br>     EXX-fraction              =        0.25<br>     Any further DFT definition will be discarded<br>     Please, verify this is what you really want<br><br>     Message from routine setup :<br>     BEWARE: nonlinear core correction is not consistent with hybrid XC<br><br>     Subspace diagonalization in iterative solution of the eigenvalue problem:<br>     a serial algorithm will be used<br><br>     Estimated static dynamical RAM per process >     771.51 MB<br><br>     Estimated max dynamical RAM per process >     774.46 MB<br><br>     Estimated total dynamical RAM >      30.25 GB<br><br>     Check: negative core charge=   -0.000018<br><br>     Initial potential from superposition of free atoms<br><br>     starting charge   71.99836, renormalised to   72.00000<br>     Starting wfcs are   52 randomized atomic wfcs<br><br>     total cpu time spent up to now is        3.9 secs<br><br>     Self-consistent Calculation<br><br>     iteration #  1     ecut=    60.00 Ry     beta= 0.40<br>     Davidson diagonalization with overlap<br>     ethr =  1.00E-02,  avg # of iterations =  2.5<br><br>     Threshold (ethr) on eigenvalues was too large:<br>     Diagonalizing with lowered threshold<br><br>     Davidson diagonalization with overlap<br>     ethr =  4.52E-04,  avg # of iterations =  2.0<br><br>     total cpu time spent up to now is       14.0 secs<br><br>     total energy              =    -536.55746986 Ry<br>     Harris-Foulkes estimate   =    -536.72101146 Ry<br>     estimated scf accuracy    <       0.32597706 Ry<br><br>     iteration #  2     ecut=    60.00 Ry     beta= 0.40<br>     Davidson diagonalization with overlap<br>     ethr =  4.53E-04,  avg # of iterations =  2.1<br><br>     total cpu time spent up to now is       19.7 secs<br><br>     total energy              =    -536.61905325 Ry<br>     Harris-Foulkes estimate   =    -536.62228319 Ry<br>     estimated scf accuracy    <       0.02002484 Ry<br><br>     iteration #  3     ecut=    60.00 Ry     beta= 0.40<br>     Davidson diagonalization with overlap<br>     ethr =  2.78E-05,  avg # of iterations =  4.0<br><br>     total cpu time spent up to now is       25.9 secs<br><br>     total energy              =    -536.62107999 Ry<br>     Harris-Foulkes estimate   =    -536.62149111 Ry<br>     estimated scf accuracy    <       0.00571317 Ry<br><br>     iteration #  4     ecut=    60.00 Ry     beta= 0.40<br>     Davidson diagonalization with overlap<br>     ethr =  7.93E-06,  avg # of iterations =  2.4<br><br>     total cpu time spent up to now is       31.0 secs<br><br>     total energy              =    -536.62144049 Ry<br>     Harris-Foulkes estimate   =    -536.62142198 Ry<br>     estimated scf accuracy    <       0.00134111 Ry<br><br>     iteration #  5     ecut=    60.00 Ry     beta= 0.40<br>     Davidson diagonalization with overlap<br>     ethr =  1.86E-06,  avg # of iterations =  2.5<br><br>     total cpu time spent up to now is       36.0 secs<br><br>     total energy              =    -536.62150194 Ry<br>     Harris-Foulkes estimate   =    -536.62150114 Ry<br>     estimated scf accuracy    <       0.00001534 Ry<br><br>     iteration #  6     ecut=    60.00 Ry     beta= 0.40<br>     Davidson diagonalization with overlap<br>     ethr =  2.13E-08,  avg # of iterations =  3.1<br><br>     total cpu time spent up to now is       43.2 secs<br><br>     total energy              =    -536.62150954 Ry<br>     Harris-Foulkes estimate   =    -536.62151057 Ry<br>     estimated scf accuracy    <       0.00000476 Ry<br><br>     iteration #  7     ecut=    60.00 Ry     beta= 0.40<br>     Davidson diagonalization with overlap<br>     ethr =  6.61E-09,  avg # of iterations =  3.0<br><br>     total cpu time spent up to now is       49.0 secs<br><br>     total energy              =    -536.62151076 Ry<br>     Harris-Foulkes estimate   =    -536.62151119 Ry<br>     estimated scf accuracy    <       0.00000130 Ry<br><br>     iteration #  8     ecut=    60.00 Ry     beta= 0.40<br>     Davidson diagonalization with overlap<br>     ethr =  1.80E-09,  avg # of iterations =  2.3<br><br>     total cpu time spent up to now is       54.6 secs<br><br>     total energy              =    -536.62151109 Ry<br>     Harris-Foulkes estimate   =    -536.62151107 Ry<br>     estimated scf accuracy    <       0.00000007 Ry<br><br>     iteration #  9     ecut=    60.00 Ry     beta= 0.40<br>     Davidson diagonalization with overlap<br>     ethr =  9.83E-11,  avg # of iterations =  3.0<br><br>     total cpu time spent up to now is       60.7 secs<br><br>     total energy              =    -536.62151111 Ry<br>     Harris-Foulkes estimate   =    -536.62151111 Ry<br>     estimated scf accuracy    <          7.3E-09 Ry<br><br>     iteration # 10     ecut=    60.00 Ry     beta= 0.40<br>     Davidson diagonalization with overlap<br>     ethr =  1.01E-11,  avg # of iterations =  2.0<br><br>     total cpu time spent up to now is       66.0 secs<br><br>     End of self-consistent calculation<br><br><br>     occupation numbers <br>     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br>     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br>     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br>     1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000   1.0000<br>     1.0000   1.0000   1.0001   1.0001   0.0000   0.0000   0.0000   0.0000<br><br>     the Fermi energy is     4.8584 ev<br><br>!    total energy              =    -536.62151111 Ry<br>     Harris-Foulkes estimate   =    -536.62151111 Ry<br>     estimated scf accuracy    <          9.4E-10 Ry<br><br>     convergence has been achieved in  10 iterations<br><br>     Using ACE for calculation of exact exchange<br><br>     EXX grid:   114089 G-vectors     FFT dimensions: (  64,  64,  64)</div></div>