<div dir="ltr">Hi Rustam, <div><br></div><div>the Mg basis set has functions with rather large exponent, this requires high cutoff with GPW. </div><div>The alternative is to use GAPW, with the same PP and BS, then the hard part of the density is not collocated on the regular grids and a lower cutoff can be used. </div><div>However, GAPW has some overheads with respect to GPW. </div><div>Possibly another solution are the NLCC pseudo potentials, but you should check and test what is available. </div><div><br></div><div>cheers</div><div>Marcella<br><br>On Friday, June 24, 2016 at 7:47:04 PM UTC+2, Rustam wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><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><span>&KIND Mg</span></p>
<p><span>  BASIS_SET DZVP-GTH</span></p>
<p><span>  POTENTIAL GTH-BLYP </span># 2s2 2p6 3s2 configuration</p><p><span>&END</span></p><p><span><br></span></p><p><span>instead of typical</span></p><p><span><br></span></p><p><span>&KIND Mg</span></p><p><span>  BASIS_SET DZVP-GTH</span></p><p><span>  POTENTIAL GTH-BLYP-q2 # 3s2 </span></p><p><span></span></p><p><span>&END</span></p></div><div><span><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></blockquote></div></div>