<div dir="ltr">There is one more question...<div><br><div><br></div><div><b>6) No Electrostatics</b><br><div>If I want a system with no electrostatics at all (e.g. a dummy molecules), is this the right way to do it?</div><div><br></div><div><div><font size="1"> &POISSON</font></div><div><font size="1"> PERIODIC none</font></div><div><font size="1"> &EWALD </font></div><div><font size="1"> EWALD_TYPE none </font></div><div><font size="1"> &END EWALD </font></div><div><font size="1"> &END POISSON</font></div></div><div><br></div><div><br></div><div>Of course I would also be very happy if someone would only answer one or two of the above questions. <br></div><div><br></div><div>Many thanks,</div><div>Jadzia</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br>On Wednesday, September 27, 2017 at 3:29:47 PM UTC+2, Jadzia wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Hi everyone, <div><br></div><div>I am trying to simulate a system with an isolated molecule, therefore I need a non-periodic setup. I have read the previous posts about the topic, but there are a number of questions open. They are both for the wavelet and the MT Poisson solver. </div><div><br></div><div><b>1) System size</b></div><div>My simulations will be quite long (both MM and QM) and due to the movement of the molecule, to my understanding it needs to be avoided that the molecule leaves the box. Will it be fine if I set the box size to 1000 Angstrom for security reasons or is this not recommended? </div><div><br></div><div><b>2) Periodic cells</b></div><div>This question is closely related to question 1). If I set the Poisson solver to non-periodic, could one set the CELL to periodic, and thus solve the problem that the molecule leaves the box?</div><div><br></div><div><b>3) Ewald methods</b></div><div>Also closely related to question 1). If I set the box size to 1000, even with Ewald as Poisson solver the system should be effectively non-periodic due to the large distance, right? Would this be a plausible alternative to the MT or wavelet solvers?</div><div><br></div><div><b>4) Mixed force evaluations</b></div><div>In the case that I have a mixed force evaluation setup and use the Wavelet solver (see the config file below), do I need to include the &CENTER_COORDINATES section also in the mixed force eval section, or only in the real force evaluation section?</div><div><br></div><div><b>5) MT solver for periodic systems</b></div><div>Is the MT solver only for nonperiodic systems in CP2K, or does it also work for periodic systems? In the publication it is at least written it works for periodic systems and even with ewald/spme. </div><div><br></div><div>Looking forward to hearing from you.</div><div><br></div><div>Best wishes,</div><div>Jadzia</div><font size="1"><br><br></font><div dir="ltr"><font size="1">&GLOBAL</font></div><div dir="ltr"><font size="1"> PROJECT system</font></div><div dir="ltr"><font size="1"> RUN_TYPE md</font></div><div dir="ltr"><font size="1"> PRINT_LEVEL medium</font></div><div dir="ltr"><font size="1">&END global</font></div><div dir="ltr"><font size="1"><br></font></div><div dir="ltr"><font size="1">&MOTION</font></div><div dir="ltr"><font size="1"> &MD</font></div><div dir="ltr"><font size="1"> ENSEMBLE langevin</font></div><div dir="ltr"><font size="1"> &LANGEVIN </font></div><div dir="ltr"><font size="1"> GAMMA [fs^-1] 1.0E-3 # "Langevin" gamma (adapt for your system) </font></div><div dir="ltr"><font size="1"> &END LANGEVIN</font></div><div dir="ltr"><font size="1"> STEPS 5000000</font></div><div dir="ltr"><font size="1"> TIMESTEP 0.5</font></div><div dir="ltr"><font size="1"> TEMPERATURE 300.0</font></div><div dir="ltr"><font size="1"> &END MD </font></div><div dir="ltr"><font size="1">&END MOTION</font></div><div dir="ltr"><font size="1"><br></font></div><div dir="ltr"><font size="1">&MULTIPLE_FORCE_EVALS</font></div><div dir="ltr"><font size="1"> FORCE_EVAL_ORDER 2 3</font></div><div dir="ltr"><font size="1"> MULTIPLE_SUBSYS true</font></div><div dir="ltr"><font size="1">&END</font></div><div dir="ltr"><font size="1"><br></font></div><div dir="ltr"><font size="1">&FORCE_EVAL</font></div><div dir="ltr"><font size="1"> STRESS_TENSOR analytical</font></div><div dir="ltr"><font size="1"> METHOD mixed</font></div><div dir="ltr"><font size="1"> &MIXED</font></div><div dir="ltr"><font size="1"> MIXING_TYPE genmix</font></div><div dir="ltr"><font size="1"> &GENERIC</font></div><div dir="ltr"><font size="1"> VARIABLES a b</font></div><div dir="ltr"><font size="1"> MIXING_FUNCTION a+b</font></div><div dir="ltr"><font size="1"> &END</font></div><div dir="ltr"><font size="1"> @include cp2k.in.mapping</font></div><div dir="ltr"><font size="1"> &END</font></div><div dir="ltr"><font size="1"> &SUBSYS</font></div><div dir="ltr"><font size="1"> &TOPOLOGY</font></div><div dir="ltr"><font size="1"> &CENTER_COORDINATES </font></div><div dir="ltr"><font size="1"> &END </font></div><div dir="ltr"><font size="1"> CONNECTIVITY off</font></div><div dir="ltr"><font size="1"> COORDINATE pdb</font></div><div dir="ltr"><font size="1"> COORD_FILE_NAME system.pdb</font></div><div dir="ltr"><font size="1"> &END TOPOLOGY</font></div><div dir="ltr"><font size="1"> &CELL</font></div><div dir="ltr"><font size="1"> PERIODIC NONE</font></div><div dir="ltr"><font size="1"> ABC 20 20 20 </font></div><div dir="ltr"><font size="1"> &END CELL</font></div><div dir="ltr"><font size="1"> &END SUBSYS</font></div><div dir="ltr"><font size="1">&END FORCE_EVAL</font></div><div dir="ltr"><font size="1"><br></font></div><div dir="ltr"><font size="1">&FORCE_EVAL</font></div><div dir="ltr"><font size="1"> STRESS_TENSOR analytical</font></div><div dir="ltr"><font size="1"> METHOD fist</font></div><div dir="ltr"><font size="1"> &MM</font></div><div dir="ltr"><font size="1"> &FORCEFIELD</font></div><div dir="ltr"><font size="1"> PARM_FILE_NAME system1.prm</font></div><div dir="ltr"><font size="1"> PARMTYPE chm</font></div><div dir="ltr"><font size="1"> &SPLINE</font></div><div dir="ltr"><font size="1"> EMAX_SPLINE 100000</font></div><div dir="ltr"><font size="1"> RCUT_NB 12.</font></div><div dir="ltr"><font size="1"> &END</font></div><div dir="ltr"><font size="1"> &END FORCEFIELD </font></div><div dir="ltr"><font size="1"> &POISSON</font></div><div dir="ltr"><font size="1"> POISSON_SOLVER WAVELET</font></div><div dir="ltr"><font size="1"> PERIODIC NONE</font></div><div dir="ltr"><font size="1"> &EWALD </font></div><div dir="ltr"><font size="1"> EWALD_TYPE NONE </font></div><div dir="ltr"><font size="1"> &END EWALD </font></div><div dir="ltr"><font size="1"> &END POISSON</font></div><div dir="ltr"><font size="1"> &END MM</font></div><div dir="ltr"><font size="1"> &SUBSYS</font></div><div dir="ltr"><font size="1"> &CELL</font></div><div dir="ltr"><font size="1"> PERIODIC NONE</font></div><div dir="ltr"><font size="1"> ABC 20 20 20</font></div><div dir="ltr"><font size="1"> &END CELL</font></div><div dir="ltr"><font size="1"> &TOPOLOGY</font></div><div dir="ltr"><font size="1"> &CENTER_COORDINATES </font></div><div dir="ltr"><font size="1"> &END </font></div><div dir="ltr"><font size="1"> COORDINATE pdb</font></div><div dir="ltr"><font size="1"> COORD_FILE_NAME system1.pdb</font></div><div dir="ltr"><font size="1"> CONN_FILE_FORMAT psf</font></div><div dir="ltr"><font size="1"> CONN_FILE_NAME system1.psf</font></div><div dir="ltr"><font size="1"> &END TOPOLOGY</font></div><div dir="ltr"><font size="1"> &END SUBSYS</font></div><div dir="ltr"><font size="1">&END FORCE_EVAL</font></div><div dir="ltr"><font size="1"><br></font></div><div dir="ltr"><font size="1">&FORCE_EVAL</font></div><div dir="ltr"><font size="1"> STRESS_TENSOR analytical</font></div><div dir="ltr"><font size="1"> METHOD fist</font></div><div dir="ltr"><font size="1"> &MM</font></div><div dir="ltr"><font size="1"> &FORCEFIELD</font></div><div dir="ltr"><font size="1"> PARM_FILE_NAME system2.prm</font></div><div dir="ltr"><font size="1"> PARMTYPE chm</font></div><div dir="ltr"><font size="1"> &SPLINE</font></div><div dir="ltr"><font size="1"> EMAX_SPLINE 100000</font></div><div dir="ltr"><font size="1"> RCUT_NB 12.</font></div><div dir="ltr"><font size="1"> &END</font></div><div dir="ltr"><font size="1"> &END FORCEFIELD </font></div><div dir="ltr"><font size="1"> &POISSON</font></div><div dir="ltr"><font size="1"> POISSON_SOLVER WAVELET</font></div><div dir="ltr"><font size="1"> PERIODIC NONE</font></div><div dir="ltr"><font size="1"> &EWALD </font></div><div dir="ltr"><font size="1"> EWALD_TYPE NONE </font></div><div dir="ltr"><font size="1"> &END EWALD </font></div><div dir="ltr"><font size="1"> &END POISSON</font></div><div dir="ltr"><font size="1"> &END MM</font></div><div dir="ltr"><font size="1"> &SUBSYS</font></div><div dir="ltr"><font size="1"> &CELL</font></div><div dir="ltr"><font size="1"> PERIODIC NONE</font></div><div dir="ltr"><font size="1"> ABC 20 20 20</font></div><div dir="ltr"><font size="1"> &END CELL</font></div><div dir="ltr"><font size="1"> &TOPOLOGY</font></div><div dir="ltr"><font size="1"> &CENTER_COORDINATES </font></div><div dir="ltr"><font size="1"> &END </font></div><div dir="ltr"><font size="1"> COORDINATE pdb</font></div><div dir="ltr"><font size="1"> COORD_FILE_NAME system2.pdb</font></div><div dir="ltr"><font size="1"> CONN_FILE_FORMAT psf</font></div><div dir="ltr"><font size="1"> CONN_FILE_NAME system2.psf</font></div><div dir="ltr"><font size="1"> &END TOPOLOGY</font></div><div dir="ltr"><font size="1"> &END SUBSYS</font></div><div dir="ltr"><font size="1">&END FORCE_EVAL</font></div></div></blockquote></div></div></div></div>