<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.EmailStyle19
{mso-style-type:personal;
font-family:"Calibri",sans-serif;
color:windowtext;}
span.EmailStyle21
{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">I realized I did not include that additional flag environ_restart for the nscf runs. Thank you.<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 levels are closer in value now, but the ‘pbc correction’ is still only output for the scf run and not the nscf run.<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">scf gives a Fermi level of -3.9197 eV which needs to be corrected by -0.1819 eV, and nscf gives only a Fermi level of -3.9185 eV without the
correction. Please advise, which value should we use.<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<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 10:20 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">I am sorry I misunderstood your previous message. The discrepancy that worries you is between nscf and scf, not between dielectric and vacuum, correct? My guess is that you are not running the nscf with Environ. You need to specify environ_restart
= true in the Environ input file, and maybe you should also set environ_thr to a large value, otherwise PW will not pass through the Environ calculation. Let me know if this solves the problem or if you had already setup the calculation like this. As for the
printout of the Fermi correction for the parabolic correction, I am not sure why is not in the output, it may be due to the problem above or a different issue, but it should be the same as in the scf calculation. <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Best,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Oliviero<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
On Mar 5, 2019, at 9:08 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 class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA">Dear Professor Andreussi,</span><o:p></o:p></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"> </span><o:p></o:p></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.</span><o:p></o:p></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"> </span><o:p></o:p></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". </span><o:p></o:p></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"> </span><o:p></o:p></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. </span><o:p></o:p></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"> </span><o:p></o:p></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.</span><o:p></o:p></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"> </span><o:p></o:p></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA">Best,</span><o:p></o:p></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"> </span><o:p></o:p></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA">Sankha Mukherjee</span><o:p></o:p></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA"> </span><o:p></o:p></p>
<p class="MsoNormal" style="background:white"><span style="font-size:12.0pt;color:black;mso-fareast-language:EN-CA">University of Toronto</span><o:p></o:p></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 <<a href="mailto:users-bounces@lists.quantum-espresso.org">users-bounces@lists.quantum-espresso.org</a>>
<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 <<a href="mailto:users@lists.quantum-espresso.org">users@lists.quantum-espresso.org</a>><br>
<b>Subject:</b> Re: [QE-users] [EXT] Help with Environ</span><o:p></o:p></p>
</div>
</div>
<p class="MsoNormal"> <o:p></o:p></p>
<div>
<p class="MsoNormal">Hi Sankha,<o:p></o:p></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&data=02%7C01%7Coliviero.andreussi%40unt.edu%7Ce451276c3b014173e18c08d6a1e0ef71%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C636874384818282402&sdata=YCdrBpofHsJ5ek5mJbqb5%2BqwcT%2FH1KpIFhcxW55Swgk%3D&reserved=0">https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&data=02%7C01%7Coliviero.andreussi%40unt.edu%7C04ab224297344351ea8f08d6a1d498ff%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C636874331832731528&sdata=8%2FhZKvMMmvN3HlyMSyMYnVmhIRCMasv0DNnCsVBdOF0%3D&reserved=0</a></span><o:p></o:p></p>
</div>
</blockquote>
</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&data=02%7C01%7Coliviero.andreussi%40unt.edu%7Ce451276c3b014173e18c08d6a1e0ef71%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C636874384818302410&sdata=EE%2FGhlij9%2F7DRuL2FYNO36VJdp0nKYP53FkLLTxdAjk%3D&reserved=0">https://nam04.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.quantum-espresso.org%2Fmailman%2Flistinfo%2Fusers&data=02%7C01%7Coliviero.andreussi%40unt.edu%7Ce451276c3b014173e18c08d6a1e0ef71%7C70de199207c6480fa318a1afcba03983%7C0%7C0%7C636874384818302410&sdata=EE%2FGhlij9%2F7DRuL2FYNO36VJdp0nKYP53FkLLTxdAjk%3D&reserved=0</a><o:p></o:p></span></p>
</div>
</blockquote>
</div>
</body>
</html>