<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:\5B8B\4F53;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:\7B49\7EBF;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"\@\5B8B\4F53";
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"\@\7B49\7EBF";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
text-align:justify;
text-justify:inter-ideograph;
font-size:10.5pt;
font-family:\7B49\7EBF;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
pre
{mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:\5B8B\4F53;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:\7B49\7EBF;
color:windowtext;}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:\5B8B\4F53;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:\7B49\7EBF;}
/* Page Definitions */
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 90.0pt 72.0pt 90.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="ZH-CN" link="#0563C1" vlink="#954F72" style="word-wrap:break-word;text-justify-trim:punctuation">
<div class="WordSection1">
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:#201F1E;background:white">Dear QE developers,<o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:#201F1E;background:white"><o:p> </o:p></span></pre>
<p class="MsoNormal" style="line-height:120%"><span lang="EN-US" style="font-size:12.0pt;line-height:120%;font-family:"Times New Roman",serif;color:#201F1E;background:white">I am a postdoc researcher working at the Brookhaven National Lab under the supervision
of Dr. Deyu Lu. My current research project involves benchmark calculations of Ti K-edge XANES using the XSpectra code (QE 6.5 and 6.6), which was compiled on our local cluster using intel/psxe suite, mvapich2 and intel mkl library. This project is in collaboration
with Dr. John Vinson from NIST and Exciting developers, including Dr. Christian Vorwerk from U. Chicago, Benedikt Maurer and Fabian Peschel from
</span><span lang="EN-US"><a href="https://scholar.google.com/citations?view_op=view_org&hl=en&org=15466593082447993016"><span style="font-size:12.0pt;line-height:120%;font-family:"Times New Roman",serif;color:#201F1E;background:white">Humboldt-Universität
zu Berlin</span></a></span><span lang="EN-US" style="font-size:12.0pt;line-height:120%;font-family:"Times New Roman",serif;color:#201F1E;background:white">.<o:p></o:p></span></p>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:#201F1E;background:white">I noticed the default of the XSpectra parameter </span><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black;background:#D9D9D9">r_paw(1:...)</span><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"> was set to <span style="background:#D9D9D9">1.5*rc</span>, where rc is the cutoff radius for the pseudo-potential. XSpectra performs a pretreatment of the pseudopotential in two steps in subroutine <span style="background:#D9D9D9">init_gipaw_1.f90</span>: First, it normalizes the all electron partial wave functions. Second, it calculates the projectors using the normalized partial waves according to the dual orthonormality. This pretreatment is done with the default paw cutoff radius at 1.5 times the cutoff stored in the paw pseudopotential file (<span style="background:#D9D9D9">r_paw(1:...) = 1.5*rc</span>).<o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black">To understand the effects of the cutoff radius scaling, we did several tests for Ti K-edge XANES. We considered two Ti core hole PAW pseudopotentials (see attachment). The first one is generated in our group that has two 4p projectors, with cutoff radius at 1.75 and 1.79 Bohr and reference energy at 7.06 and 0.78 Ryd, respectively. This potential was used previously for the XANES simulations published in Nano Letters 19 (6), 3457-3463, 2019. In the second potential, we only kept one 4p projector with reference energy at 0.78 Ryd and cutoff radius at 1.79 Bohr.<o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black">To understand the effects of normalization of the all-electron partial wave and projector generation, we consider two cases. In the first case referred to as <b>original</b>, we modified the source code (init_gipaw_1.f90 and paw_gipaw.f90; attached) to disable the partial wave normalization and projector generation and take the partial waves and projectors as those in the pseudopotential file. In the second case referred to as <b>regenerated</b>, we follow the standard XSpectra pretreatment of the partial wave and projector. In the regenerated mode, we denote r_paw = 1.0*rc as <b>hard cutoff</b> and r_paw = 1.5*rc as <b>soft cutoff</b> by specifying r_paw in the input file using standard XSpectra code.<o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p> </o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p> </o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p> </o:p></span></pre>
<pre style="text-align:center;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif"><img border="0" width="449" height="407" style="width:4.677in;height:4.2343in" id="Picture_x0020_4" src="cid:image001.jpg@01D828B3.381077D0" alt="Chart, line chart
Description automatically generated"></span><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black">Figure 1 Comparison of projectors (a and c) and the partial waves (b and d) stored in the pseudopotential file with one 4p projector (labelled as original) with those regenerated in Xspectra (labelled as regenerated). (a) and (b): hard cutoff; (c) and (d): soft cutoff.<o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p> </o:p></span></pre>
<pre style="text-align:center;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif"><img border="0" width="449" height="404" style="width:4.677in;height:4.2031in" id="Picture_x0020_3" src="cid:image002.jpg@01D828B3.381077D0" alt="Diagram
Description automatically generated"></span><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black">Figure 2 Comparison of projectors (a and c) and the partial waves (b and d) stored in the pseudopotential file for the Ti core hole pseudopotential with two projectors.<o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p> </o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black">Next we compare partial waves and projectors obtained for anatase TiO<sub>2</sub> (mp-390) in Figure 1 (with one 4p projector) and Figure 2 (with two 4p projectors). We note that the results shouldn’t depend on the system we choose, as it only concerns the PAW pseudopotential. We can see that the shape of partial waves has a trivial dependence on r_paw and the PSP pretreatment, and the differences only come from a scaling factor. However, the shape of the projector shows much stronger dependence on these two factors, due to the enforced orthogonality condition applied onto matrix S in subroutine init_gipaw_1.f90. As shown in Fig. 1a with hard cutoff, regenerated projector has a sharp drop to 0 at rc, while the tail in the original projector is smooth. The same trend is also found in Fig. 2a. With the soft cutoff in Figs. 1c and 2c, one can see that the generated projectors are extended to 1.5 rc. While the original projectors are mostly out-of-phase in Fig. 2c, the regenerated projectors become mostly in-phase below 2 Bohr.<o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black">These subtle changes in partial waves and projectors can have a non-trivial impact on calculated XANES spectra. We performed three calculations for anatase TiO<sub>2</sub> (mp-390): i) original, ii) regenerated hard cutoff and iii) regenerated and soft cutoff. The calculated XANES spectra are shown in Figure 3. <o:p></o:p></span></pre>
<pre style="text-align:center;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif"><img border="0" width="554" height="229" style="width:5.7656in;height:2.3802in" id="Picture_x0020_2" src="cid:image003.jpg@01D828B3.381077D0" alt="Chart, line chart, histogram
Description automatically generated"></span><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black">Figure 3 The Ti K-edge XANES spectra of anatase TiO<sub>2</sub> (mp-390) calculated using (a) one and (b) two 4p projectors.<o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p> </o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black">We consider the reference spectrum as those calculated from the original mode, in which rc is always the same as in the PSP file and the projector from the PSP file goes to zero at 1.0*rc. However, after the pretreatment, we can see differences in intensities compared with the original spectrum. In Fig 3a, both regenerated hard cutoff and soft cutoff will lead to noticeable differences compared to the reference spectra throughout the full energy range. When using two 4p projectors, however only the spectra calculated soft cutoff deviates from the original spectrum, while the spectra calculated with hard cutoff is identical compared with the reference as shown in Fig 3b. The reason why the PSP pretreatment behaves differently with 1 or 2 4p projectors has not been fully understood.<o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black">In addition, we compared the XANES spectra for anatase TiO<sub>2 </sub>(mp-390), rutile TiO<sub>2 </sub>(mp-2657) and brookite TiO<sub>2</sub> (mp-1840) obtained with XSpectra and with OCEAN and EXCITING, which are two codes widely used to calculate the XANES spectra by solving BSE. In this test, we only used the pseudopotential with two 4p projectors. Since XSpectra outputs the cross-section while OCEAN and EXCITING output the imaginary part of the dielectric constant, the XSpectra data are rescaled in order to compare the results directly with the other two codes as shown in Fig 4. Overall, OCEAN and EXCITING match very well as expected since they are both BSE codes. It is also clear that XANES spectra obtained from XSpectra with hard cutoff agree better than with soft cutoff, especially in terms of peak intensities. <o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p> </o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US"><img border="0" width="554" height="424" style="width:5.7656in;height:4.4114in" id="Picture_x0020_1" src="cid:image004.jpg@01D828B3.381077D0" alt="Chart
Description automatically generated"></span><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black">Figure 4 Comparison of XANES spectra obtained with XSpectra, OCEAN and EXCITING. (a) Anatase TiO<sub>2</sub>, (b) Rutile TiO<sub>2 </sub>and (c) Brookite TiO<sub>2.<o:p></o:p></sub></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><sub><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p> </o:p></span></sub></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black">Based on our tests results, we found that the XSpectra results are more consistent without the reconstruction of the all-electron partial waves and the projectors in the core-hole potential. Therefore, we believe that the default r_paw=1.5*rc together with PSP pretreatment are likely introduce some artifact in the XANES spectra such that the main edge intensity is weakened, as compared to the original spectra from XSpectra and reference spectra from other codes. We appreciate your feedback on this matter.<o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p> </o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black">Thanks,<o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black">Fanchen<o:p></o:p></span></pre>
<pre style="text-align:justify;text-justify:inter-ideograph;line-height:120%"><span lang="EN-US" style="font-family:"Times New Roman",serif;color:black"><o:p> </o:p></span></pre>
<p class="MsoNormal" style="line-height:120%"><span lang="EN-US" style="font-size:12.0pt;line-height:120%;font-family:"Times New Roman",serif"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
</div>
</body>
</html>