<div dir="ltr"><div>Hi  Mohnish,<br><br><br></div>Thank you very much for your reply. After I check my input to:<br><br><div>&control<br>    prefix='Zr',<br>    pseudo_dir = '/home/zhang/Software/espresso_PP/',<br>

    outdir='/share/scratch/tmpZr/eCut30SM0.01/',<br> /<br> &system    <br>    ibrav=  1, celldm(1) =50, nat=  1, ntyp= 1,<br>    occupations='smearing', smearing='marzari-vanderbilt', degauss=0.01,<br>

    ecutwfc =30, <br>    nspin=2,  starting_magnetization(1)=0.0<br> /<br> &electrons<br> /<br>ATOMIC_SPECIES<br> Zr  91.224 Zr.pbe-mt_fhi.UPF<br>ATOMIC_POSITIONS  <br>  Zr  0.000000000000 0.000000000000 0.000000000000<br>

K_POINTS automatic<br>  18  18  18 0 0 0<br><br><br></div><div>The program just hung there. The following is the last part of output:<br><br><br><br>     Pseudo is Norm-conserving, Zval =  4.0<br>     Generated using FHI98PP, converted with fhi2upf.x v.5.0.2<br>

     Using radial grid of  539 points,  3 beta functions with: <br>                l(1) =   0<br>                l(2) =   2<br>                l(3) =   3<br><br>     atomic species   valence    mass     pseudopotential<br>

        Zr             4.00    91.22400     Zr( 1.00)<br><br>     Starting magnetic structure <br>     atomic species   magnetization<br>        Zr           0.000<br><br>     48 Sym. Ops., with inversion, found<br><br><br>

<br>   Cartesian axes<br><br>     site n.     atom                  positions (alat units)<br>         1           Zr  tau(   1) = (   0.0000000   0.0000000   0.0000000  )<br><br>     number of k points=   440  Marzari-Vanderbilt smearing, width (Ry)=  0.0100<br>

<br>     Number of k-points >= 100: set verbosity='high' to print them.<br><br>     Dense  grid:  2774807 G-vectors     FFT dimensions: ( 180, 180, 180)<br><br>     Largest allocated arrays     est. size (Mb)     dimensions<br>

        Kohn-Sham Wavefunctions         7.94 Mb     (   86754,    6)<br>        NL pseudopotentials            17.21 Mb     (   86754,   13)<br>        Each V/rho on FFT grid         44.49 Mb     ( 1458000,   2)<br>        Each G-vector array             5.29 Mb     (  693702)<br>

        G-vector shells                 0.05 Mb     (    6335)<br>     Largest temporary arrays     est. size (Mb)     dimensions<br>        Auxiliary wavefunctions        31.77 Mb     (   86754,   24)<br>        Each subspace H/S matrix        0.01 Mb     (      24,   24)<br>

        Each <psi_i|beta_j> matrix      0.00 Mb     (      13,    6)<br>        Arrays for rho mixing         177.98 Mb     ( 1458000,    8)<br><br>     Initial potential from superposition of free atoms<br>     Check: negative starting charge=(component1):   -0.000223<br>

     Check: negative starting charge=(component2):   -0.000223<br><br>     starting charge    3.99791, renormalised to    4.00000<br><br>     negative rho (up, down):  2.231E-04 2.231E-04<br>     Starting wfc are   16 randomized atomic wfcs<br>

<br></div><div>Could you please tell me if there is any mistake in my input file? BTW, I can't get your reply from my email address. Is it possible to get the copy of reply when anyone replies to the forum? <br><br></div>

<div>Best,<br><br>Dongsheng<br></div><div><br><br></div><div><b></b></div></div>