Dear Professor Andreussi, Thank you for your detailed response.
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". 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. 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. Best, Sankha Mukherjee University of Toronto From: users <[email protected]> On Behalf Of Andreussi, Oliviero Sent: March-05-19 8:57 PM To: Quantum Espresso users Forum <[email protected]> Subject: Re: [QE-users] [EXT] Help with Environ Hi Sankha, 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. I hope this helps, if you have other problems or doubts, please let me know. Regards, Oliviero Andreussi University of North Texas On Mar 5, 2019, at 7:39 PM, Sankha Mukherjee <[email protected]<mailto:[email protected]>> wrote: Hello all, 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. 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. Thanks, Sankha Mukherjee _______________________________________________ users mailing list [email protected]<mailto:[email protected]> 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
_______________________________________________ users mailing list [email protected] https://lists.quantum-espresso.org/mailman/listinfo/users
