<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I fear your enterprise is destined for failure.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Performing periodic calculations on charged materials does not converge. They can be made to converge by adding a jellium neutralising background, but I do not think that it has been proven that that is an accurate approach. It's worse for slabs, where the
results depend explicitly on the amount of vacuum between images.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I do not have a viable alternative suggestion.</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Ciao</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Terry</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
UNSW Canberra<br>
</div>
<div>
<div id="appendonsend"></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt; color:rgb(0,0,0)">
<br>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> users <users-bounces@lists.quantum-espresso.org> on behalf of Zhenyao Fang <zhenyaof@sas.upenn.edu><br>
<b>Sent:</b> Tuesday, 17 August 2021 9:30 AM<br>
<b>To:</b> users@lists.quantum-espresso.org <users@lists.quantum-espresso.org><br>
<b>Subject:</b> [QE-users] Total energy of charged slabs</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt">
<div class="PlainText">Dear pwscf users,<br>
<br>
I am performing some calculations on negatively charged slabs (with one additional electron), and I’m mainly focusing on the energy difference between charged system and neutral system. With ions fixed, I found that the charged system has a higher energy than
the neutral system, instead of being lower. <br>
<br>
The energy difference dE = E(q=-1) - E(q=0) is around 3.944 eV. The LUMO position of the neutral slab is 3.462 eV, and the vacuum level for the neutral slab is 7.654 eV. Besides, by inspecting the charge difference plots, the additional charge is localized
inside the slab. I understand that there could be problems about jellium background, so I used the package “sxdefectalign2d” to correct the interactions due to periodic images, and the total energy of the charged slab remains the same as the size of vacuum
increases. <br>
<br>
I’m quite confused by the fact that the total energy would increase with one more electron. Since the additional charge is mainly inside the slab, it should lower the energy and be stabilised inside the slab. Therefore, I was wondering how I can reconcile these
controversies. Besides, another related question is how we define the zero energy in QE? In other words, does LUMO larger or smaller than zero possibly imply the stability of the system?<br>
<br>
I am looking forward to your replies. Thanks in advance.<br>
<br>
Best,<br>
Zhenyao Fang<br>
<br>
<br>
<br>
<br>
<br>
_______________________________________________<br>
Quantum ESPRESSO is supported by MaX (<a href="http://www.max-centre.eu">www.max-centre.eu</a>)<br>
users mailing list users@lists.quantum-espresso.org<br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/users">https://lists.quantum-espresso.org/mailman/listinfo/users</a></div>
</span></font></div>
</div>
</body>
</html>