<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">I forgot that a problem may arise if the wavefunctions are not bound: even it the 4p levels are not bound, the code can use them to generate the beta projector since it needs the wavefunctions only up to the core radius. If you want to make projections on it you'll need instead to integrate the wavefunction also beyond the core radius. If the wavefunction is not bound this integral will not make much sense. If this is the case, what you can do is to fix the 4p pseudization energy to a small negative value. This might however complicate or hinder the PP generation procedure...<div><br></div><div><br></div><div>GS<br><div><br></div><div><br><div><div>Il giorno 12/nov/2010, alle ore 09.41, Gabriele Sclauzero ha scritto:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Dear Eduardo,<div><br></div><div> I think that the PP you are using does not include the 4p wavefunctions</div><div><span class="Apple-style-span" style="font-family: Times; "><pre style="word-wrap: break-word; white-space: pre-wrap; "><PP_HEADER>
0 Version Number
Ni Element
US Ultrasoft pseudopotential
T Nonlinear Core Correction
SLA PW PBE PBE PBE Exchange-Correlation functional
10.00000000000 Z valence
-72.87827664100 Total energy
0.0000000 0.0000000 Suggested cutoff for wfc and rho
2 Max angular momentum component
1203 Number of points in mesh
2 6 Number of Wavefunctions, Number of Projectors
Wavefunctions nl l occ
4S 0 1.00
3D 2 9.00
</PP_HEADER></pre><pre style="word-wrap: break-word; white-space: pre-wrap; "><br></pre></span><div>probably because the 4p states have 0 occupation and the ld1.x code stores in the PP file only atomic wavefunctions with occupation > 0. Hence the projections on the 4p atomic wavefunctions cannot be computed because they are not available to the code.</div><div>However, I think that the 4p wavefunctions are still computed by ld1.x when generating the PP, since there are beta projectors on the l=1 channel (grep BETA inside the PP file). You could try to modify the default behavior of ld1.x so that it will include those wavefunctions in the PP. Then the projection code should detect them in the PP data and compute the additional projections.</div><div><br></div><div><br></div><div>HTH</div><div><br></div><div>GS</div><div><br></div><div><div>Il giorno 11/nov/2010, alle ore 21.45, Eduardo Ariel Menendez Proupin ha scritto:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div><div>Hi,</div><div>I am doing a tutorial example of calculation of the projected DOS in Ni. It produces the PDOS projected on S and D states, but not on P. </div><div>There are no bound P valence states in nickel. Is that the cause for not having a P partial DOS? </div>
<div><br></div><div>States with P-symmetry component may be present in an empty band. Is there an option to produce the P-PDOS?</div><div><br></div><div><br></div><div>Support information:</div><div>this is from the help file INPUT_PROJWFC.txt</div>
<div><br></div><div> In the collinear case and the non-collinear, non spin-orbit case</div><div> projected DOS are written to file "filpdos".pdos_atm#N(X)_wfc#M(l),</div><div> where N = atom number , X = atom symbol, M = wfc number, l=s,p,d,f</div>
<div> (one file per atomic wavefunction found in the pseudopotential file)</div></div><div><br></div><div><br></div><div>this is the head of the pseudopotential</div><div><PP_INFO></div><div>Generated using Andrea Dal Corso code (rrkj3) </div>
<div>Author: Andrea Dal Corso Generation date: unknown </div><div>Info: Ni </div><div> 1 The Pseudo was generated with a Scalar-Relativistic Calculation</div>
<div> 1.70000000000E+00 Local Potential cutoff radius</div><div>nl pn l occ Rcut Rcut US E pseu</div><div>4S 1 0 1.00 2.00000000000 2.50000000000 0.00000000000</div>
<div>4S 1 0 0.00 2.00000000000 2.50000000000 0.00000000000</div><div>4P 2 1 0.00 2.40000000000 2.60000000000 0.00000000000</div><div>4P 2 1 0.00 2.40000000000 2.60000000000 0.00000000000</div>
<div>3D 3 2 9.00 1.60000000000 2.50000000000 0.00000000000</div><div>3D 3 2 0.00 1.60000000000 2.50000000000 0.00000000000</div><div></PP_INFO></div><div><br></div><div><br></div>
Thanks<br>-- <br><div><br></div>
<div><br></div>Eduardo Menendez<br>Departamento de Fisica<br>Facultad de Ciencias<br>Universidad de Chile<br>Phone: (56)(2)9787439<br>URL: <a href="http://fisica.ciencias.uchile.cl/~emenendez" target="_blank">http://fisica.ciencias.uchile.cl/~emenendez</a>
<div><br></div>
<div>They did it!</div><br>
_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum">http://www.democritos.it/mailman/listinfo/pw_forum</a><br></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div><span class="Apple-style-span" style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic; "><br class="Apple-interchange-newline">§ Gabriele Sclauzero, EPFL SB ITP CSEA</span></div><div><font class="Apple-style-span" color="#7E7E7E"><i> PH H2 462, Station 3, CH-1015 Lausanne</i></font></div></span>
</div>
<br></div></div>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>http://www.democritos.it/mailman/listinfo/pw_forum<br></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div><span class="Apple-style-span" style="color: rgb(126, 126, 126); font-size: 16px; font-style: italic; "><br class="Apple-interchange-newline">§ Gabriele Sclauzero, EPFL SB ITP CSEA</span></div><div><font class="Apple-style-span" color="#7E7E7E"><i> PH H2 462, Station 3, CH-1015 Lausanne</i></font></div></span>
</div>
<br></div></div></body></html>