Dear All<div><br></div><div>I require the average potential energy per water molecule for thermochemical calculations. However, I find that the average PE differs significantly (by ~1-8 kJ/mol) for different system sizes. Based on classical simulations, this is quite unexpected. There is the possibility that the simulations have not equilibrated yet, but the PEs seem to have reasonably plateaued (see attached figure). The length of the trajectories are also already similar to those reported in (older) papers.   </div><div><br></div><div>The DFT settings used for each system are the same, and are summarized by the input lines appended to the end of this email. In all cases the temperature set to 330 K. For reference, I've also attached restart files for each system. </div><div><br></div><div>Any insight as to why these finite size effects are observed would me much appreciated. Could there be a "DFT reason" that I am not aware of?</div><div><br></div><div>Best Regards,</div><div>James</div><div>________________________________________</div><div><br></div><div><div>     &MGRID</div><div>       CUTOFF     3.2000000000000000E+02</div><div>       REL_CUTOFF     4.0000000000000000E+01</div><div>     &END MGRID</div><div>     &XC</div><div>       DENSITY_CUTOFF     1.0000000000000000E-10</div><div>       GRADIENT_CUTOFF     1.0000000000000000E-10</div><div>       TAU_CUTOFF     1.0000000000000000E-10</div><div>       &XC_GRID</div><div>         XC_SMOOTH_RHO  NN10</div><div>         XC_DERIV  SPLINE2_SMOOTH</div><div>       &END XC_GRID</div><div>       &XC_FUNCTIONAL  NO_SHORTCUT</div><div>         &PBE  T</div><div>         &END PBE</div><div>       &END XC_FUNCTIONAL</div><div>       &VDW_POTENTIAL</div><div>         POTENTIAL_TYPE  PAIR_POTENTIAL</div><div>         &PAIR_POTENTIAL</div><div>           R_CUTOFF     1.6000000000000000E+01</div><div>           TYPE  DFTD3(BJ)</div><div>           PARAMETER_FILE_NAME ../input/dftd3.dat</div><div>           REFERENCE_FUNCTIONAL PBE</div><div>         &END PAIR_POTENTIAL</div><div>       &END VDW_POTENTIAL</div></div><div><br></div><div><br></div><div><br></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion on the web visit <a href="https://groups.google.com/d/msgid/cp2k/1f3b60dd-d06f-4b10-810d-1bc9cde4d9e1n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/1f3b60dd-d06f-4b10-810d-1bc9cde4d9e1n%40googlegroups.com</a>.<br />