<html>
<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;}
p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph
{mso-style-priority:34;
margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:36.0pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
/* List Definitions */
@list l0
{mso-list-id:1903326374;
mso-list-type:hybrid;
mso-list-template-ids:-562012854 536870929 536870937 536870939 536870927 536870937 536870939 536870927 536870937 536870939;}
@list l0:level1
{mso-level-text:"%1\)";
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level2
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level3
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
text-indent:-9.0pt;}
@list l0:level4
{mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level5
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level6
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
text-indent:-9.0pt;}
@list l0:level7
{mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level8
{mso-level-number-format:alpha-lower;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;}
@list l0:level9
{mso-level-number-format:roman-lower;
mso-level-tab-stop:none;
mso-level-number-position:right;
text-indent:-9.0pt;}
ol
{margin-bottom:0cm;}
ul
{margin-bottom:0cm;}
--></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-BE" link="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">Dear QE users,<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">I’m currently performing convergence tests for my system by looking at the energies, forces and stresses for different values of ecutwfc, ecutrho and k-points to get the right optimized parameters for my calculations.
However, during my research I will be mainly focused on properties of the electron density as I’ll be performing QTAIM, ELF and NCI calculations, so I would like to converge my parameters with respect to the electron density distribution of my system but I’m
unsure about the correct procedure in this situation. I have a couple of things I’m uncertain about:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<ol style="margin-top:0cm" start="1" type="1">
<li class="MsoListParagraph" style="margin-left:0cm;mso-list:l0 level1 lfo1"><span lang="EN-US">I use the pp.x code to create a .cube file containing the electron density of my system, which uses the values of nr1,nr2,nr3 to construct the .cube file. The values
of nr1,nr2,nr3 are themselves dependent on the values of ecutwfc and ecutrho, so by varying these values , I would obtain different grids in consecutive calculations, which makes it hard to compare the results. Would manually assigning values to nr1,nr2,nr3
during the SCF calculations be a valid option to ensure the grids of the obtained cube files are the same? Or could this hinder the convergence as the FFT-grids are now sub-optimal in some cases?<o:p></o:p></span></li><li class="MsoListParagraph" style="margin-left:0cm;mso-list:l0 level1 lfo1"><span lang="EN-US">Once the .cube files for the different calculations are generated, how would one choose which values to compare to get an idea of the convergence? Would it be sufficient
to choose a given point inside the .cube file and see how this value converges with the given parameters, or would you have to (in some way) look at all the values inside the .cube file to assess convergence? Perhaps you could do this by calculating the deviation
in the electron density at each point in the .cube file w.r.t. some reference calculation with very high values for ecutwfc, ecutrho and k-points and then taking the average of the deviation at all points inside the .cube file to get an actual number to compare
between calculations.<o:p></o:p></span></li></ol>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">If anyone has some experience with convergence testing w.r.t the electron density I would be very happy to get some help on the subject!<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">Thanks in advance and happy holidays,<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>
<p class="MsoNormal"><span lang="EN-US">Belgium<o:p></o:p></span></p>
</div>
</body>
</html>