<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">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">fname.scf.in</a> > $fname.scf.out<br><br><br>mkdir pdos<br>cat > <a href="http://proj.in">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">proj.in</a> > proj.out<br clear="all"><div><div><div><div><div><br>-- <br><div class="gmail_signature"><div dir="ltr">--<br>Dhirendra<br></div></div>
</div></div></div></div></div></div>