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>