<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><br class="webkit-block-placeholder"></div><div><br class="webkit-block-placeholder"></div>Hi Axel, <div><br class="webkit-block-placeholder"></div><div>Playing with the GAPW parameters requires some insight in the formalism.</div><div>I will try to summarize briefly the role of the most relevant GAPW parameters.</div><div>You are right they are grouped in DFT%QS and in KIND</div><div>In QS section</div><div>- <span class="Apple-style-span" style="font-family: Times; font-size: 16px; -webkit-border-horizontal-spacing: 2px; -webkit-border-vertical-spacing: 2px; ">EPSFIT: </span>this parameter is used to determine which primitive Gaussians are considered hard and which one soft. With the default values exponents >= 4 are considered to be hard. </div><div>- <span class="Apple-style-span" style="font-family: Times; font-size: 16px; -webkit-border-horizontal-spacing: 2px; -webkit-border-vertical-spacing: 2px; ">EPSISO, EPS_SVD and MAX_RAD_LOCAL <span class="Apple-style-span" style="font-family: Helvetica; font-size: 12px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; ">are parameters used in the construction of the GAPW projectors. Playing with them might help in difficult cases that turn out to have instabilities. Otherwise they can be safely left to default values</span></span></div><div>- in GAPW local compensation densities are needed in order to compute correctly the hartree contributions that involve the local densities. In the original scheme we have a hard compensation density, expanded in Gaussians with fixed exponent <span class="Apple-style-span" style="font-family: Times; font-size: 16px; -webkit-border-horizontal-spacing: 2px; -webkit-border-vertical-spacing: 2px; ">ALPHA0_HARD, <span class="Apple-style-span" style="font-family: Helvetica; font-size: 12px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; ">and a soft compensation density that has the same multipole expansion as the hard one and it is use also on the fft. Having two different compensation densities introduces some additional terms in the expression of the total energy. These terms  are rather expensive to be computed because they involve three center neighbor lists. By default  ALPHA0_HARD is set to 10</span></span></div><div>and ALPHA0_SOFT=ALPHA0_HARD, so that all the complex 3-center terms disappear.  Since we noticed very small differences by having or not having these terms, by default we set equal hard and soft exponents ? IF FULL_GAPW is set to true, all the terms are computed. If the exponents are not given from input, they are automatically set to meaningful values. </div><div>- LMAXN0 and LMAXN1 are maximum angular momenta used in the expansion of compensation charge and local charges. The maximum possible value  is 2 times the maximum angular momentum present in the basis set. Having LMAXN0 or LMAXN1 smaller than this can improve the performance, but might have an effect on the accuracy.</div><div>-<span class="Apple-style-span" style="font-family: Times; font-size: 16px; -webkit-border-horizontal-spacing: 2px; -webkit-border-vertical-spacing: 2px; ">QUADRATURE <span class="Apple-style-span" style="font-family: Helvetica; font-size: 12px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; "> this variable defines the algorithm used to assign the points on  the local radial grid</span></span></div><div><br class="webkit-block-placeholder"></div><div><br class="webkit-block-placeholder"></div><div>in KIND</div><div>HARD_EXP_RADIUS is the radius of the region around the atom where the local densities are defined. The  spherical regions centered on neighboring atoms do not have to overlap, therefore the radius cannot be too extended..</div><div><br class="webkit-block-placeholder"></div><div><span class="Apple-style-span" style="font-family: Times; font-size: 16px; -webkit-border-horizontal-spacing: 2px; -webkit-border-vertical-spacing: 2px; ">LEBEDEV_GRID RADIAL_GRID </span>are parameters that determine se size of the local grids, The optimal choice  depends on the system and on the maximum angular momentum that need to be collocated. For not too complex systems  one can use about 100 points on the radial grids, and 50 on the angular.</div><div><br class="webkit-block-placeholder"></div><div><br class="webkit-block-placeholder"></div><div>Other useful info:</div><div>- GAPW does not need a large PW cut-off, since the hard density is computed on the local atomic grids.</div><div>- Comparing GAPW and Gaussian calculations for small molecules you should find differences of a few microhartree. Due to the definition of the local densities, GAPW improves systematically by improving the basis set</div><div>- if you use all electron and  OT, convergence is significantly improved by setting the preconditioner to FULL_ALL</div><div><br class="webkit-block-placeholder"></div><div>I am available for further clarifications on the method and its  implementation.</div><div><br class="webkit-block-placeholder"></div><div>ciao</div><div> Marcella</div><div><br class="webkit-block-placeholder"></div><div><br><div><div>On 10 Mar 2008, at 18:35, Axel wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><br>hi!<br><br>i'd appreciate it, if somebody could shed some light on how to<br>justify<br>parameter choices for GAPW calculations. the input reference mostly<br>contains statements of the kind XXX sets the value for 'xxx', which<br>are<br>probably of most use to the person who knows the implemention very<br>well<br>(which i obviously don't :-) ).<br><br>from looking through the regtests and the input reference it seems to<br>me<br>that ALHPA0_HARD and ALHPA0_SOFT, LMAXN0, LMAXN1 in %FORCE_EVAL%DFT<br>%QS<br>seem to be the most important global options, and that LEBEDEV_GRID,<br>RADIAL_GRID, and HARD_EXP_RADIUS in %FORCE_EVAL%SUBSYS%KIND seem to be<br>the most<br>important per atom type parameters.<br><br>are there any do's and don'ts associated with those parameters, e.g.<br>what needs to<br>be changed when going from a bulk water system to a compound with<br>"heavy" elements?<br>which of those parameters are related to higher angular momentum and<br>thus need to<br>be increased or decreased? can the HARD_EXP_RADIUS and its sibling<br>RHO0_EXP_RADIUS<br>be understood in a similar way as the ultrasoft cutoff radius for<br>USPPs?<br><br>thanks in advance,<br>     axel.<br><br><br><br></blockquote></div><br></div></body></html>