<div dir="ltr">Dear Lorenzo,<div><br><div>Thank you for your email and suggestion.</div><div><br></div><div>The calculation flow for XAS with Xspexctra as follows:</div><div>Firstly, prepare the GIPAW Pseudopotential.</div><div>Secondly, extract the core wavefunction</div><div>Thirdly, prepare input file and run SCF</div><div>Finally, prepare input file and run Xspectra.</div><div><br></div><div>Considering the calculation flow, the problem (or unclear part) is to get the core wavefunction. </div><div>As far as I understand, there are two ways to extract core wave functions:</div><div><br></div><div>One is to use "upf2plotcore.sh" script. (upf2plotcore.sh < Ga.PBE.PP.UPF > Ga.wfc),  which is your suggestion. Am I correct? </div><div>If yes, </div><div>upf2plotcore.sh < Ga.PBE.PP.UPF(without core hole) > Ga.wfc<br></div><div>upf2plotcore.sh < Ga.PBE.PP.UPF(with core hole) > h_Ga.wfc.<br></div><div>which one is correct for spectroscopy.</div><div><br></div><div>The other is to extract from all electron (AE) wavefunction which also includes core wave functions. So, can we use the second way to extract core wavefunction? </div><div>If yes, why the output file by ld1.x does not include 1S core wfc information using my input.</div><div>Your thoughts are welcomed.</div><div><br></div><div>Best regards,</div><div><br></div><div>Mehmet Emin Kilic, PhD</div><div>-----------------------------</div><div><pre style="white-space:pre-wrap;color:rgb(0,0,0)">Hello,
If you want to do spectroscopy with the xspectra code, it does not work
like that. You do not need an all-electron pseudopotential, but a pseudo
generated with a 1s1 core.

-- 
Lorenzo Paulatto

On Fri, 31 May 2019, 07:35 emin klc, <<a href="https://lists.quantum-espresso.org/mailman/listinfo/users" target="_blank">klcmemin at gmail.com</a>> wrote:

><i> Dear QE experts,
</i>><i>
</i>><i> I have tried to generate a Ga Pseudopotentials for theoretical
</i>><i> spectroscopy by ld1.x.
</i>><i> Firstly, I need to have all electron (AE) wave functions including 1s 2s
</i>><i> 2p 3s 3p core states and valence states.
</i>><i> I'm having some troubles since the output file (ld1.wfc)  only gives 2s 2p
</i>><i> 3s 3p as a core states except for 1s, but which is important for K edge
</i>><i> spectroscopy.
</i>><i> I include below my input file for generating all electron wave functions.
</i>><i>
</i>><i> Could you say what is wrong with my input?
</i>><i> -----------------------
</i>><i> &input
</i>><i>       title = 'Ga',
</i>><i>       iswitch = 1,
</i>><i>       rel = 1,
</i>><i>       zed  = 31.0,
</i>><i>       config = '[Ar] 4s2 4p1 3d10.0',
</i>><i>       dft = 'PBE'
</i>><i> ------------------------
</i>><i>
</i>><i> Mehmet Emin Kilic, PhD
</i>></pre></div><div class="gmail-yj6qo"></div><br class="gmail-Apple-interchange-newline"></div></div>