<div dir="ltr"><div><div><div>In serial execution, and for the last QE version, I cannot reproduce your problem (that error message should NEVER appear). Note however that<br></div>- 40 Ry for Ni HGH is way too low, and<br></div>- There are no atomic wavefunctions in the HGH pseudopotential files, so you cannot compute the PDOS<br><br></div>Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Sep 30, 2015 at 2:19 PM, DHIRENDRA VAIDYA <span dir="ltr"><<a href="mailto:dhirendra22121987@gmail.com" target="_blank">dhirendra22121987@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div>Hi,<br><br></div>I am trying to plot density of state for Nickel Silicide.<br><br></div>I am getting this error<br><br>     Error in routine  descla_init (2):<br>      wrong value for computed nrcx<br><br></div>What is it?<br><br></div>My script is belo<br><br>fname=NiSi2_Si<br>OutDir=NiSi2_Si<br>PREFIX=NiSi2_Si<br><br><br>cat > $<a href="http://fname.scf.in" target="_blank">fname.scf.in</a> <<EOF<br>&control<br> calculation = 'scf',<br> verbosity = 'low',<br> restart_mode = 'from_scratch',<br> prefix = '$PREFIX',<br> tstress = .TRUE.,<br> tprnfor = .TRUE.,<br> pseudo_dir = '$PSEUDO_DIR/',<br> outdir = '$OutDir'<br>/ <br><br>&system<br> ibrav=2,<br> celldm(1) = 10.2381,<br> nat=3, ntyp=2,<br> ecutwfc=40,<br> nbnd=20<br> occupations='smearing', degauss=0.01<br>/<br><br>&electrons<br> diagonalization='david',<br> mixing_mode = 'plain',<br> mixing_beta = 0.7,<br> conv_thr =  1.0d-6,<br>/ <br>&ions<br>/<br>&cell<br>cell_factor=2.5<br>/<br>ATOMIC_SPECIES<br> Ni 58.693400 Ni.pz-hgh.UPF <br> Si 28.085500 Si.pz-hgh.UPF<br><br>ATOMIC_POSITIONS crystal<br> Ni 0.000000 0.000000 0.000000<br> Si 0.250000 0.250000 0.250000<br> Si 0.750000 0.750000 0.750000<br><br><br>K_POINTS automatic<br>15 15 15 0 0 0<br><br><br><br>EOF<br> $PW_command < $<a href="http://fname.scf.in" target="_blank">fname.scf.in</a> > $fname.scf.out<br><br><br>mkdir pdos<br>cat > <a href="http://proj.in" target="_blank">proj.in</a> <<EOF<br>&projwfc<br>prefix = '$PREFIX',<br>outdir = '$OutDir'<br>ngauss=0<br>degauss=0.00735<br>filpdos='./pdos/dos.txt'<br>filproj='./pdos/proj.txt'<br>DeltaE=0.001,<br>Emin=-15,<br>Emax=15,<br>/<br>EOF<br><br>$proj_command < <a href="http://proj.in" target="_blank">proj.in</a> > proj.out<span class="HOEnZb"><font color="#888888"><br clear="all"><div><div><div><div><div><br>-- <br><div><div dir="ltr">--<br>Dhirendra<br></div></div>
</div></div></div></div></div></font></span></div>
<br>_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span><span><font color="#888888">Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a></font></span></span></div></div></div></div>
</div>