<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=iso-8859-1">
<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;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
span.EmailStyle18
{mso-style-type:personal-compose;
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-GB" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">Hi Giovanni,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Thank you for your answer! I do have some follow-up questions if you don’t mind.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">When converging the clean slab, which electronic properties do you propose to look at for convergence? Because up until now, I have only done convergence tests by looking at the total energies of my system… I come from
molecular calculations, so I’m still a bit unfamiliar with things like the density of states and band structures. I’m mainly concerned about adsorption energies and analyzing the bonds inside the system by looking at the electron density and its derivatives
(i.e. theories like QTAIM, ELF, NCI).<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">You also say that I should increase my unit cell size in the direction perpendicular to the slab, but can you just change just change your unit cell without consequences? I have been told that all your calculations have
to be done in the same unit cell to be able to compare the values of properties such as the adsorption energy… Wouldn’t changing the unit cell also change the plane wave basis set (i.e. larger unit cell = denser reciprocal lattice = more plane waves for a
similar cutoff)?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">For the dipole correction, I currently use the assume_isolated = ‘2D’ flag enabled, while freezing the atoms in the bottom row of my slab to their bulk positions to simulate the bulk.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Kind regards,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Léon Luntadila Lufungula<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Structural Chemistry group<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">University of Antwerp<o:p></o:p></span></p>
</div>
</body>
</html>