<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"CMU Bright Roman";
panose-1:2 0 6 3 0 0 0 0 0 0;}
@font-face
{font-family:"Times New Roman \(Body CS\)";
panose-1:2 11 6 4 2 2 2 2 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:12.0pt;
font-family:"Calibri",sans-serif;
mso-ligatures:standardcontextual;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"CMU Bright Roman";
color:windowtext;
font-weight:normal;
font-style:normal;
text-decoration:none none;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:12.0pt;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman"">Dear all,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman""><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman"">Is it possible to obtain the (pseudo-orbital approximation to the) all-electron density from a norm-conserving pseudopotential calculation? I thought at first that rho_core ought
to provide the rest of the approximated density, but (rho%of_r + rho_core) doesn’t seem to integrate to the total number of electrons in the small case I tried, Si with 2x2x2 k-points: rho%of_r integrates to 4, but rho%of_r + rho_core to somewhat less than
4.75. <o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman""><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman"">Is rho_core only for nonlinear core correction? Must one use projector-augmented waves or ultrasoft PPs to get explicit core-electron properties?<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman""><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman"">Of course, I do not expect the core density to be very accurate! But I do need to access some effect from the core anyway. If it helps, I’m using norm-conserving Vanderbilt optimized
PPs from abinit and PBE, with scalar relativistic correction; my ecutwfc = 75 Ry. I’m also doing LSDA calculations, so the actual quantity I’ve been looking at is 0.5 * (rho_core(:) + rho%of_r(:, 1) +/- rho%of_r(:, 2)).<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman""><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman"">Thanks for your time!<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman""><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman";color:black;mso-ligatures:none">Sincerely yours,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman";color:black;mso-ligatures:none"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman";color:black;mso-ligatures:none">Jacob Williams<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman";color:black;mso-ligatures:none">Ph.D. candidate, Yang group<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"CMU Bright Roman";color:black;mso-ligatures:none">Duke University Department of Chemistry
</span><o:p></o:p></p>
</div>
</body>
</html>