<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7638.1">
<TITLE>Calculating Projected density of states</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->
<BR>

<P><FONT SIZE=2>Hi QEs users,<BR>
I have tried to calculate projected density of states first by performing scf calculation, then nscf calculation and then Projected density of states, but I am getting the following error.<BR>
<BR>
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<BR>
&nbsp;&nbsp;&nbsp;&nbsp; task #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<BR>
&nbsp;&nbsp;&nbsp;&nbsp; from projwfc : error #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4324<BR>
&nbsp;&nbsp;&nbsp;&nbsp; reading inputpp namelist<BR>
<BR>
&nbsp;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<BR>
<BR>
my pdos input file is<BR>
<BR>
&amp;inputpp<BR>
&nbsp; prefix = 'unoptimisedsilane',<BR>
&nbsp; outdir = '/sfiwork/dsharma/PDOSSilane/scfn/scfunopt'<BR>
&nbsp;degauss = 0.0,<BR>
&nbsp;Emin = -12.6653,<BR>
&nbsp;Emax = 12,<BR>
ngauss = 0,<BR>
DeltaE = 0.01,<BR>
filpdos = 'pdosout',<BR>
<BR>
&nbsp;and my script file is<BR>
<BR>
#!/bin/bash<BR>
#PBS -N pdos<BR>
#PBS -l nodes=4<BR>
##PBS -1 waltime=24:00:00<BR>
#PBS -q walton<BR>
#PBS -V<BR>
<BR>
## #PATH of the parallel abinit executable<BR>
export pw_path='/usr/local/espresso-4.0.5/bin'<BR>
<BR>
### PATH to the mpirun binary:<BR>
export mpirun_path='/opt/mpi/bin'<BR>
<BR>
export calc_name='pdos'<BR>
<BR>
### go in the working directory:<BR>
cd $PBS_O_WORKDIR<BR>
<BR>
## Count the number of cpu requested<BR>
NP=`wc -l $PBS_NODEFILE`<BR>
echo $NP &gt; $calc_name.txt<BR>
<BR>
$mpirun_path/mpirun -machinefile $PBS_NODEFILE -np $NP $pw_path/projwfc.x -npool1 &lt; pdos &gt; pdosout<BR>
<BR>
My scf output file is<BR>
<BR>
&nbsp;&nbsp; End of self-consistent calculation<BR>
<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k = 0.0000 0.0000 0.0000 (&nbsp; 6983 PWs)&nbsp;&nbsp; bands (ev):<BR>
<BR>
&nbsp;&nbsp; -12.6653&nbsp; -9.5380&nbsp; -7.2586&nbsp; -6.0976&nbsp; -5.0908&nbsp; -3.8305&nbsp; -0.6212&nbsp;&nbsp; 0.5917<BR>
&nbsp;&nbsp;&nbsp;&nbsp; 0.8368&nbsp;&nbsp; 1.6400<BR>
<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k = 0.2500 0.0000 0.0000 (&nbsp; 6959 PWs)&nbsp;&nbsp; bands (ev):<BR>
<BR>
&nbsp;&nbsp; -12.1861&nbsp; -9.7498&nbsp; -6.9945&nbsp; -6.2237&nbsp; -5.5081&nbsp; -5.1444&nbsp;&nbsp; 0.3355&nbsp;&nbsp; 1.0465<BR>
&nbsp;&nbsp;&nbsp;&nbsp; 1.4373&nbsp;&nbsp; 1.9502<BR>
<BR>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; k =-0.5000 0.0000 0.0000 (&nbsp; 6966 PWs)&nbsp;&nbsp; bands (ev):<BR>
<BR>
&nbsp;&nbsp; -10.9343 -10.8970&nbsp; -6.6489&nbsp; -6.3817&nbsp; -5.9690&nbsp; -5.5686&nbsp;&nbsp; 0.3422&nbsp;&nbsp; 0.8657<BR>
&nbsp;&nbsp;&nbsp;&nbsp; 1.1277&nbsp;&nbsp; 1.7350<BR>
<BR>
&nbsp;&nbsp;&nbsp;&nbsp; highest occupied, lowest unoccupied level (ev):&nbsp;&nbsp;&nbsp; -3.8305&nbsp;&nbsp; -0.6212<BR>
<BR>
One more thing I want to ask, I have choosen the Emin and Emax value from the scf calculation, is this correct?<BR>
<BR>
suggestion welcome<BR>
<BR>
Thanks and regards<BR>
<BR>
Dimpy<BR>
</FONT>
</P>

</BODY>
</HTML>