<div dir="ltr">Hi guys,<div><br></div><div>I need to include inner-shell electrons of Mg into my GPW calculations. That is the total number of electrons on magnesium must be 10, not just 2 valence 3s electrons.</div><div><br></div><div>This can be readily achieved using by using the following pseudopotential</div><div><br></div><div>







<p class="p1"><span class="s1">&KIND Mg</span></p>
<p class="p1"><span class="s1">  BASIS_SET DZVP-GTH</span></p>
<p class="p1"><span class="s1">  POTENTIAL GTH-BLYP </span># 2s2 2p6 3s2 configuration</p><p class="p1"><span class="s1">&END</span></p><p class="p1"><span class="s1"><br></span></p><p class="p1"><span class="s1">instead of typical</span></p><p class="p1"><span class="s1"><br></span></p><p class="p1"><span class="s1">&KIND Mg</span></p><p class="p1"><span class="s1">  BASIS_SET DZVP-GTH</span></p><p class="p1"><span class="s1">  POTENTIAL GTH-BLYP-q2 # 3s2 </span></p><p class="p1"><span class="s1"></span></p><p class="p1"><span class="s1">&END</span></p></div><div><span class="s1"><br></span></div><div>I figured (by trial and error) that including inner electrons requires significantly higher CUTOFF setting. The attached plot shows energy dependence on the cell size of MgH2.</div><div><br></div><div>A. 2-electron calculations are fine with 1000 Ry cutoff</div><div>b. 10-electron calculations produce awful results unless the cutoff is above 8000 Ry.</div><div><br></div><div>QUESTION: Is it really necessary to use such a high cutoff to get smooth curves? Or are there any other settings one can tweak to get accurate results with lower (i.e. cheaper) cutoff?</div><div><br></div><div>Thanks!</div><div>Rustam</div><div><br></div></div>