<html xmlns:v="urn:schemas-microsoft-com:vml" 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=utf-8">
<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:"Calibri Light";
        panose-1:2 15 3 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0cm;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;
        mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:#0563C1;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:#954F72;
        text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
        {mso-style-name:msonormal;
        mso-margin-top-alt:auto;
        margin-right:0cm;
        mso-margin-bottom-alt:auto;
        margin-left:0cm;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
span.EmailStyle18
        {mso-style-type:personal;
        font-family:"Calibri",sans-serif;
        color:windowtext;}
span.EmailStyle20
        {mso-style-type:personal-reply;
        font-family:"Calibri",sans-serif;
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:72.0pt 72.0pt 72.0pt 72.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="EN-CA" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA">Dear Professor Andreussi,<o:p></o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"><o:p> </o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA">Thank you for your detailed response.<o:p></o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"><o:p> </o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA">I would just like to confirm that the correction for which you are referring is the line that states "the Fermi energy shift due to the parabolic
 pbc-correction is  -0.1819 ev". <o:p></o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"><o:p> </o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA">For the nscf output, however, I cannot seem to find the equivalent correction line, and the Fermi energy is output without any additional lines
 after it. <o:p></o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"><o:p> </o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA">The printed Fermi energy in the scf run file is -3.9197 eV, and with the correction becomes -4.1016 eV. However, in the nscf output the Fermi
 energy is printed as -1.7056 eV.<o:p></o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"><o:p> </o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA">Best,<o:p></o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"><o:p> </o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA">Sankha Mukherjee<o:p></o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"><o:p> </o:p></span></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA">University of Toronto<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span lang="EN-US" style="mso-fareast-language:EN-CA">From:</span></b><span lang="EN-US" style="mso-fareast-language:EN-CA"> users <users-bounces@lists.quantum-espresso.org>
<b>On Behalf Of </b>Andreussi, Oliviero<br>
<b>Sent:</b> March-05-19 8:57 PM<br>
<b>To:</b> Quantum Espresso users Forum <users@lists.quantum-espresso.org><br>
<b>Subject:</b> Re: [QE-users] [EXT] Help with Environ<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Hi Sankha,<span style="mso-fareast-language:EN-CA"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Thanks for reporting your doubts with Environ. Similarly to the potential, the Fermi energy is defined up to an additive constant which is not uniquely defined. If you introduce the dielectric, the Fermi energy changes. If you look at the
 value of the Fermi energy with respect to the value of the potential in the vacuum, I.e. the workfunction, and you compare it to the same quantity in solution (a.k.a. the potential of zero charge) you should expect them to be different. Having said this, the
 way the potential is shifted inside of Environ is probably different from a standard calculation in vacuum and depends on the pbc correction that you use. Please pay attention to the fact that in an Environ calculation there is a correction that needs to be
 added to the Fermi energy and is printed out in the output a few lines after the Fermi energy itself. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I hope this helps, if you have other problems or doubts, please let me know.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Regards,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Oliviero Andreussi<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">University of North Texas<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
On Mar 5, 2019, at 7:39 PM, Sankha Mukherjee <<a href="mailto:sankha.mukherjee@utoronto.ca">sankha.mukherjee@utoronto.ca</a>> wrote:<o:p></o:p></p>
</div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p style="margin:0cm;margin-bottom:.0001pt;background:white"><span style="font-size:12.0pt;font-family:"Calibri Light",sans-serif;color:black">Hello all,</span><o:p></o:p></p>
<p style="margin:0cm;margin-bottom:.0001pt;background:white;font-variant-ligatures: normal;font-variant-caps: normal;orphans: 2;text-align:start;widows: 2;-webkit-text-stroke-width: 0px;text-decoration-style: initial;text-decoration-color: initial;word-spacing:0px">
<span style="font-size:12.0pt;font-family:"Calibri Light",sans-serif;color:black"> </span><o:p></o:p></p>
<p style="margin:0cm;margin-bottom:.0001pt;background:white;font-variant-ligatures: normal;font-variant-caps: normal;orphans: 2;text-align:start;widows: 2;-webkit-text-stroke-width: 0px;text-decoration-style: initial;text-decoration-color: initial;word-spacing:0px">
<span style="font-size:12.0pt;font-family:"Calibri Light",sans-serif;color:black">I am working on catalysis related projects using Environ. When running SCF and subsequently NSCF using Environ I have noticed a discrepancy between the values of the Fermi levels
 when introducing a dielectric constant >1. Considering the vacuum case, these Fermi energies are nearly identical however.</span><o:p></o:p></p>
<p style="margin:0cm;margin-bottom:.0001pt;background:white;font-variant-ligatures: normal;font-variant-caps: normal;orphans: 2;text-align:start;widows: 2;-webkit-text-stroke-width: 0px;text-decoration-style: initial;text-decoration-color: initial;word-spacing:0px">
<span style="font-size:12.0pt;font-family:"Calibri Light",sans-serif;color:black"> </span><o:p></o:p></p>
<p style="margin:0cm;margin-bottom:.0001pt;background:white;font-variant-ligatures: normal;font-variant-caps: normal;orphans: 2;text-align:start;widows: 2;-webkit-text-stroke-width: 0px;text-decoration-style: initial;text-decoration-color: initial;word-spacing:0px">
<span style="font-size:12.0pt;font-family:"Calibri Light",sans-serif;color:black">I was wondering how the permittivity of the medium could be causing this contrast, and if there is one of these values I can trust more than the other.</span><o:p></o:p></p>
<p style="margin:0cm;margin-bottom:.0001pt;background:white;font-variant-ligatures: normal;font-variant-caps: normal;orphans: 2;text-align:start;widows: 2;-webkit-text-stroke-width: 0px;text-decoration-style: initial;text-decoration-color: initial;word-spacing:0px">
<span style="font-size:12.0pt;font-family:"Calibri Light",sans-serif;color:black"> </span><o:p></o:p></p>
<p style="margin:0cm;margin-bottom:.0001pt;background:white;font-variant-ligatures: normal;font-variant-caps: normal;orphans: 2;text-align:start;widows: 2;-webkit-text-stroke-width: 0px;text-decoration-style: initial;text-decoration-color: initial;word-spacing:0px">
<span style="font-size:12.0pt;font-family:"Calibri Light",sans-serif;color:black">Thanks,</span><o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Calibri Light",sans-serif;color:black;mso-fareast-language:EN-CA"> </span><o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Calibri Light",sans-serif;color:black;mso-fareast-language:EN-CA">Sankha Mukherjee</span><o:p></o:p></p>
</div>
</blockquote>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p class="MsoNormal"><span style="mso-fareast-language:EN-CA">_______________________________________________<br>
users mailing list<br>
<a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a><br>
<a href="https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&amp;data=02%7C01%7Coliviero.andreussi%40unt.edu%7C04ab224297344351ea8f08d6a1d498ff%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C636874331832731528&amp;sdata=8%2FhZKvMMmvN3HlyMSyMYnVmhIRCMasv0DNnCsVBdOF0%3D&amp;reserved=0">https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&amp;data=02%7C01%7Coliviero.andreussi%40unt.edu%7C04ab224297344351ea8f08d6a1d498ff%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C636874331832731528&amp;sdata=8%2FhZKvMMmvN3HlyMSyMYnVmhIRCMasv0DNnCsVBdOF0%3D&amp;reserved=0</a><o:p></o:p></span></p>
</div>
</blockquote>
</div>
</body>
</html>