<div dir="ltr">Dear Paolo<div><br></div><div>thanks for your reply.<br><div><br></div><div>at the end of the following document it was mentioned that </div><div><br></div><div><pre>" Obsolete variables, ignored:
io_choice
smoothing "</pre></div><div><br></div><div><a href="http://www.quantum-espresso.org/wp-content/uploads/Doc/INPUT_PROJWFC.html">http://www.quantum-espresso.org/wp-content/uploads/Doc/INPUT_PROJWFC.html</a><br><div>
<br></div><div style>in the following you can see my input files:</div><div><br></div><div><b>scf:</b></div><div><div>&control</div><div> calculation='scf',</div><div> restart_mode='restart',</div>
<div> prefix='B',</div><div> pseudo_dir = '../pseudo/',</div><div> outdir='../tmp/',</div><div> /</div><div> &system </div><div> ibrav= 0, celldm(1)=2.75899995, nat= 2, ntyp= 1,</div>
<div> ecutwfc = 40, ecutrho = 400, occupations='smearing', smearing='gaussian', degauss=0.001,</div><div> nspin=2, starting_magnetization(1)=1 </div><div> /</div><div> &electrons</div><div> conv_thr = 1.0e-10</div>
<div> mixing_beta = 0.3 </div><div> /</div><div>CELL_PARAMETERS </div><div>1.50000000<span class="" style="white-space:pre"> </span>0.86600000<span class="" style="white-space:pre"> </span>0.00000000</div><div>1.50000000 -0.86600000<span class="" style="white-space:pre"> </span>0.00000000</div>
<div>0.00000000<span class="" style="white-space:pre"> </span>0.00000000<span class="" style="white-space:pre"> </span>8.00000000</div><div>ATOMIC_SPECIES</div><div>C 12.0107 C.pbe-rrkjus.UPF </div><div>ATOMIC_POSITIONS {alat}</div>
<div>C 0.000000000 0.000000000 0.000000000 </div><div>C 1.997267588 0.000001338 0.000000000</div><div>K_POINTS {automatic}</div><div>4 4 1 0 0 0</div><div><br></div><div><br></div><div style><b>nscf:</b></div>
<div style><div>&control</div><div> calculation='nscf',</div><div> restart_mode='restart',</div><div> prefix='B',</div><div> pseudo_dir = '../pseudo/',</div><div> outdir='../tmp/',</div>
<div> wf_collect=.TRUE.</div><div> /</div><div> &system </div><div> ibrav= 0, celldm(1)=2.75899995, nat= 2, ntyp= 1,</div><div> ecutwfc = 40, ecutrho = 400, occupations='tetrahedra', </div><div> nspin=2, starting_magnetization(1)=1 </div>
<div> /</div><div> &electrons</div><div> conv_thr = 1.0e-10</div><div> mixing_beta = 0.3 </div><div> /</div><div>CELL_PARAMETERS </div><div>1.50000000<span class="" style="white-space:pre"> </span>0.86600000<span class="" style="white-space:pre"> </span>0.00000000</div>
<div>1.50000000 -0.86600000<span class="" style="white-space:pre"> </span>0.00000000</div><div>0.00000000<span class="" style="white-space:pre"> </span>0.00000000<span class="" style="white-space:pre"> </span>8.00000000</div>
<div>ATOMIC_SPECIES</div><div>C 12.0107 C.pbe-rrkjus.UPF </div><div>ATOMIC_POSITIONS {alat}</div><div>C 0.000000000 0.000000000 0.000000000 </div><div>C 1.997267588 0.000001338 0.000000000</div>
<div>K_POINTS {automatic}</div><div>20 20 1 0 0 0</div><div><br></div><div style><b>dos:</b></div><div style><div> &dos</div><div> outdir='../tmp/'</div><div> prefix='B'</div><div> fildos='B.dos',</div>
<div> Emin=-40.0, Emax=20.0, DeltaE=0.1, </div><div> /</div><div><br></div><div style><b>pdos:</b></div><div style><div> &projwfc</div><div> outdir='../tmp/'</div><div> prefix='B'</div><div>
Emin=-40.0, Emax=20.0, DeltaE=0.1</div><div> ngauss=0, degauss=0.001</div><div> /</div><div><br></div><div style>overall the only way that I can make my PDOS graphs as smooth as DOS ones is reducing degauss value to 0.01 ! but in some cases that we study the proper value for degauss is 0.001 and reducing it we may miss something. mightn't it ? </div>
<div style>We are going to compare a large sets of PDOS and DOS graphs and prefer not to compare one graph with lots of singularities with the others that are too smooth. I appreciate your advice in advance. Many Thanks</div>
<div style><br></div><div style><br></div><div style><br></div><div style><br></div><div style>-- <br>Sincerely Yours<br>David G.<br>JCU<br></div><div><br></div></div></div></div><div><br></div><div><pre style="color:rgb(0,0,0)">
On Fri, 2013-04-26 at 21:50 +0200, David Grifith wrote:
><i> when I set smoothing=0.3 as I did in the old versions of q/e the
</i>><i> recent version ignores the smoothing procedure and does nothing to
</i>><i> smooth the PDOS graph.
</i>><i>
</i>><i> I would appreciate if anyone helps me to make PDOS graphs as smooth as
</i>><i> DOS ones.</i></pre><pre style="color:rgb(0,0,0)"><i><br></i></pre><pre style="color:rgb(0,0,0)"><span style="font-family:'Times New Roman';font-size:medium;white-space:normal">on </span><i style="font-family:'Times New Roman';font-size:medium;white-space:normal">Sat Apr 27 16:35:42 CEST 2013, Paolo Giannozzi wrote: </i><span style="font-family:'Times New Roman';font-size:medium;white-space:normal"> </span><i style="font-family:'Times New Roman';font-size:medium;white-space:normal"><br>
</i><i> </i>
Not sure what you define as "smooth". The DOS in a true crystal is far
from smooth: it has singularities and you need the tetrahedron method
to have a good DOS. With gaussian broadening, you may need a relatively
large broadening or else you will get a series of spikes, but then
you miss the singularities. dos.x can calculate DOS using tetrahedra
(this is what it should do according to your input data). projwfc.x
cannot calculate PDOS with tetrahedra, if I remember correctly. More
exactly: it could, provided somebody take the existing code that
calculates the PDOS for phonons (it should be routine dos_gam in
matdyn.f90), adapts it to electrons.
Anyway: without a simple test (with outputs) and a definition of
what are the "old" and "recent" versions, it is hard to say more.
I do not remember any recent change in the DOS-producing codes
that may explain the phenomenology you describe.
Paolo
--
Paolo Giannozzi, Dept. Chemistry&Physics&Environment,
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy
Phone +39-0432-558216, fax +39-0432-558222 </pre><div><div><br></div><br>
</div></div></div></div></div></div>