I would like to repeat calculations of the static dielectric constant that someone did years ago with a different code, for use in other calculations.  My basic setup is below.<div><br></div><div>I'm not an expert and it's a challenge.  One previous problem was that I did not realize I was not calculating "traditional" dipole moments that are compatible with the fluctuation formula used in exercises like the one at the link in red below.  I now include "PERIODIC FALSE" (please see my setup at the bottom of this message) and the problem seems to have been fixed, but I mention it here in case it is not.</div><div><font color="#ff0000">https://www.cp2k.org/exercises:2018_uzh_acpc2:mol_sol</font></div><div><br></div><div><div>The issue now is that I have realized I should be using Ewald
summation
in my periodic boundary conditions.  I thought I was (that it was the default).  In my setup below that I have been using, I have:</div><div><br></div><div><p><font color="#0000ff">    &POISSON</font></p><p><font color="#0000ff">       PERIODIC XYZ </font></p><p><font color="#0000ff">    &END</font></p></div><div><br></div><div>However, I came across an example input file in a CP2K exercise online and saw that he included explicit language related to Ewald summation, so it appears I have not been using it so far.  The example uses:</div><div><br></div><div>





<p><font color="#ff0000">    &POISSON</font></p>
<p><font color="#ff0000">      &EWALD</font></p>
<p><font color="#ff0000">        EWALD_TYPE ewald</font></p>
<p><font color="#ff0000">        ALPHA .45</font></p>
<p><font color="#ff0000">        GMAX  19</font></p>
<p><font color="#ff0000">      &END EWALD</font></p>
<p><font color="#ff0000">    &END POISSON</font></p><p><b>Could someone please point me in the right direction?  </b>Should I replace my simple "PERIODIC XYZ" with the example directly above (in red)?  How can I figure out reasonable values for ALPHA and GMAX if so?  Any advice would be very helpful.  Thank you.<br></p><p><br></p><p>My setup (below).</p></div><div>





<p><font color="#0000ff" face="Arial">&GLOBAL</font></p>
<p><font color="#0000ff" face="Arial">  PROJECT X</font></p>
<p><font color="#0000ff" face="Arial">  RUN_TYPE MD</font></p>
<p><font color="#0000ff" face="Arial">&END GLOBAL</font></p>
<p><font color="#0000ff" face="Arial">&FORCE_EVAL</font></p>
<p><font color="#0000ff" face="Arial">  STRESS_TENSOR  ANALYTICAL</font></p>
<p><font color="#0000ff" face="Arial">  METHOD  Quickstep</font></p>
<p><font color="#0000ff" face="Arial">  &DFT</font></p>
<p><font color="#0000ff" face="Arial">    &PRINT</font></p>
<p><font color="#0000ff" face="Arial">      &MOMENTS</font></p>
<p><font color="#0000ff" face="Arial">      FILENAME = moments</font></p>
<p><font color="#0000ff" face="Arial">      ADD_LAST NUMERIC</font></p>
<p><font color="#0000ff" face="Arial">      PERIODIC FALSE</font></p>
<p><font color="#0000ff" face="Arial">      &EACH</font></p>
<p><font color="#0000ff" face="Arial">        MD 1</font></p>
<p><font color="#0000ff" face="Arial">      &END</font></p>
<p><font color="#0000ff" face="Arial">      &END MOMENTS</font></p>
<p><font color="#0000ff" face="Arial">    &END PRINT</font></p>
<p><font color="#0000ff" face="Arial">    BASIS_SET_FILE_NAME BASIS_MOLOPT</font></p>
<p><font color="#0000ff" face="Arial">    POTENTIAL_FILE_NAME POTENTIAL_SET</font></p>
<p><font color="#0000ff" face="Arial">    CHARGE 0</font></p>
<p><font color="#0000ff" face="Arial">    MULTIPLICITY 1</font></p>
<p><font color="#0000ff" face="Arial">    &MGRID</font></p>
<p><font color="#0000ff" face="Arial">       CUTOFF [Ry] 600</font></p>
<p><font color="#0000ff" face="Arial">       REL_CUTOFF [Ry] 60</font></p>
<p><font color="#0000ff" face="Arial">    &END</font></p>
<p><font color="#0000ff" face="Arial">    &QS</font></p>
<p><font color="#0000ff" face="Arial">       METHOD GPW</font></p>
<p><font color="#0000ff" face="Arial">       EPS_DEFAULT 1.0E-10</font></p>
<p><font color="#0000ff" face="Arial">       EXTRAPOLATION ASPC</font></p>
<p><font color="#0000ff" face="Arial">       EXTRAPOLATION_ORDER 4</font></p>
<p><font color="#0000ff" face="Arial">    &END QS</font></p>
<p><font color="#0000ff" face="Arial">    &POISSON</font></p>
<p><font color="#0000ff" face="Arial">       PERIODIC XYZ </font></p>
<p><font color="#0000ff" face="Arial">    &END</font></p>
<p><font color="#0000ff" face="Arial">    &SCF</font></p>
<p><font color="#0000ff" face="Arial">      SCF_GUESS ATOMIC </font></p>
<p><font color="#0000ff" face="Arial">      MAX_SCF 30</font></p>
<p><font color="#0000ff" face="Arial">      EPS_SCF 1.0E-6 </font></p>
<p><font color="#0000ff" face="Arial">      &OT</font></p>
<p><font color="#0000ff" face="Arial">        PRECONDITIONER FULL_SINGLE_INVERSE</font></p>
<p><font color="#0000ff" face="Arial">        MINIMIZER DIIS</font></p>
<p><font color="#0000ff" face="Arial">      &END OT</font></p>
<p><font color="#0000ff" face="Arial">      &OUTER_SCF </font></p>
<p><font color="#0000ff" face="Arial">        MAX_SCF 10</font></p>
<p><font color="#0000ff" face="Arial">        EPS_SCF 1.0E-6 </font></p>
<p><font color="#0000ff" face="Arial">      &END</font></p>
<p><font color="#0000ff" face="Arial">      &PRINT</font></p>
<p><font color="#0000ff" face="Arial">        &RESTART OFF</font></p>
<p><font color="#0000ff" face="Arial">        &END</font></p>
<p><font color="#0000ff" face="Arial">      &END</font></p>
<p><font color="#0000ff" face="Arial">    &END SCF</font></p>
<p><font color="#0000ff" face="Arial">    &XC</font></p>
<p><font color="#0000ff" face="Arial">      &XC_FUNCTIONAL</font></p>
<p><font color="#0000ff" face="Arial">         &PBE</font></p>
<p><font color="#0000ff" face="Arial">         &END</font></p>
<p><font color="#0000ff" face="Arial">      &END XC_FUNCTIONAL</font></p>
<p><font color="#0000ff" face="Arial">      &XC_GRID</font></p>
<p><font color="#0000ff" face="Arial">            XC_DERIV  SPLINE2_SMOOTH</font></p>
<p><font color="#0000ff" face="Arial">            XC_SMOOTH_RHO NN50</font></p>
<p><font color="#0000ff" face="Arial">      &END XC_GRID</font></p>
<p><font color="#0000ff" face="Arial">      &VDW_POTENTIAL</font></p>
<p><font color="#0000ff" face="Arial">         POTENTIAL_TYPE PAIR_POTENTIAL</font></p>
<p><font color="#0000ff" face="Arial">         &PAIR_POTENTIAL</font></p>
<p><font color="#0000ff" face="Arial">            PARAMETER_FILE_NAME dftd3.dat</font></p>
<p><font color="#0000ff" face="Arial">            TYPE DFTD3</font></p>
<p><font color="#0000ff" face="Arial">            REFERENCE_FUNCTIONAL PBE</font></p>
<p><font color="#0000ff" face="Arial">            R_CUTOFF [angstrom] 16</font></p>
<p><font color="#0000ff" face="Arial">         &END</font></p>
<p><font color="#0000ff" face="Arial">      &END VDW_POTENTIAL</font></p>
<p><font color="#0000ff" face="Arial">    &END XC</font></p>
<p><font color="#0000ff" face="Arial">  &END DFT</font></p>
<p><font color="#0000ff" face="Arial">  &SUBSYS</font></p>
<p><font color="#0000ff" face="Arial">    &CELL </font></p>
<p><font color="#0000ff" face="Arial">      ABC X X X</font></p>
<p><font color="#0000ff" face="Arial">    &END CELL</font></p>
<p><font color="#0000ff" face="Arial">    &COORD</font></p>
<p><font color="#0000ff" face="Arial">      XXX</font></p>
<p><font color="#0000ff" face="Arial">    &END COORD</font></p>
<p><font color="#0000ff" face="Arial">    &TOPOLOGY</font></p>
<p><font color="#0000ff" face="Arial">    &END</font></p>
<p><font color="#0000ff" face="Arial">    &KIND X</font></p>
<p><font color="#0000ff" face="Arial">      BASIS_SET DZVP-MOLOPT-SR-GTH</font></p>
<p><font color="#0000ff" face="Arial">      POTENTIAL GTH-PBE-q9</font></p>
<p><font color="#0000ff" face="Arial">    &END KIND</font></p>
<p><font color="#0000ff" face="Arial">  &END SUBSYS</font></p>
<p><font color="#0000ff" face="Arial">&END FORCE_EVAL</font></p>
<p><font color="#0000ff" face="Arial">&MOTION</font></p>
<p><font color="#0000ff" face="Arial">  &MD</font></p>
<p><font color="#0000ff" face="Arial">     XXX</font></p>
<p><font color="#0000ff" face="Arial">  &END MD</font></p>
<p><font color="#0000ff" face="Arial">  &PRINT</font></p>
<p><font color="#0000ff" face="Arial">    …</font></p>
<p><font color="#0000ff" face="Arial">  &END PRINT</font></p>
<p><font color="#0000ff" face="Arial">&END MOTION</font></p></div></div>