<div dir="ltr">Hi,<div>Thank you for the tips! I am not using EXTRAPOLATION or PRECONDITIONER currently. If I include them now will they affect the accuracy significantly? I am in the middle of long MD trajectories of TiO2 crystal with water and KCl ions. Also running the water box with ions. My current setup is below (same for both systems). The TiO2 system uses 1-80 scf steps for each MD step. Would you recommend these changes for that system also? Best regards. <div><br></div><div><br></div><div><div> &GLOBAL</div><div> PRINT_LEVEL MEDIUM</div><div> PROJECT_NAME brookite</div><div> RUN_TYPE MD</div><div> &END GLOBAL</div><div> &MOTION</div><div> &MD</div><div> ENSEMBLE NVT</div><div> STEPS 250</div><div> TIMESTEP 4.9999999999999989E-01</div><div> STEP_START_VAL 3751</div><div> TIME_START_VAL 4.3234999999995443E+03</div><div> ECONS_START_VAL -1.5206763646667368E+04</div><div> TEMPERATURE 3.0000000000000000E+02</div><div> &THERMOSTAT</div><div> TYPE NOSE</div><div> REGION GLOBAL</div><div> &NOSE</div><div> LENGTH 3</div><div> TIMECON 1.9999999999999993E+01</div><div> &COORD</div><div> 4.2157834584794962E-01 -3.8977051434234506E+00 3.0959934319656378E+02</div><div> &END COORD</div><div> &VELOCITY</div><div> 5.7927878733942134E-05 7.8527005854180164E-04 5.0048257140974906E-04</div><div> &END VELOCITY</div><div> &MASS</div><div> 1.4749962622867215E+06 6.4949196930282767E+02 6.4949196930282767E+02</div><div> &END MASS</div><div> &FORCE</div><div> 4.2152207720701857E-08 6.1579061284564603E-06 -8.4610128141038230E-07</div><div> &END FORCE</div><div> &END NOSE</div><div> &END THERMOSTAT</div><div> &AVERAGES T</div><div> &RESTART_AVERAGES</div><div> ITIMES_START 1</div><div> AVECPU 9.1875550803485453E+01</div><div> AVEHUGONIOT 0.0000000000000000E+00</div><div> AVETEMP_BARO 0.0000000000000000E+00</div><div> AVEPOT -1.5208725197223592E+04</div><div> AVEKIN 1.0826866556597368E+00</div><div> AVETEMP 3.0108762437707566E+02</div><div> AVEKIN_QM 0.0000000000000000E+00</div><div> AVETEMP_QM 0.0000000000000000E+00</div><div> AVEVOL 5.5116376409126919E+04</div><div> AVECELL_A 2.7174261790895553E+01</div><div> AVECELL_B 6.9659084710430690E+01</div><div> AVECELL_C 2.9116900255502028E+01</div><div> AVEALPHA 9.0000000000000071E+01</div><div> AVEBETA 9.0000000000000071E+01</div><div> AVEGAMMA 9.0000000000000071E+01</div><div> AVE_ECONS 1.6039009887889247E+02</div><div> AVE_PRESS 0.0000000000000000E+00</div><div> AVE_PXX 0.0000000000000000E+00</div><div> &END RESTART_AVERAGES</div><div> &END AVERAGES</div><div> &END MD</div><div> &CONSTRAINT</div><div> &FIXED_ATOMS</div><div> COMPONENTS_TO_FIX XYZ</div><div> LIST 148 172 196 150 174 198 64 88 112 \</div><div> 259 270 281 221 234 247 149 173 197 289 \</div><div> 313 337 260 271 282 225 238 251 155 179 \</div><div> 203 163 187 211 164 188 212 290 314 338 \</div><div> 71 95 119 266 277 288 265 276 287</div><div> &END FIXED_ATOMS</div><div> &END CONSTRAINT</div><div> &PRINT</div><div> &TRAJECTORY ON</div><div> ADD_LAST NUMERIC</div><div> FILENAME trajectory</div><div> &END TRAJECTORY</div><div> &RESTART SILENT</div><div> ADD_LAST NUMERIC</div><div> &EACH</div><div> MD 1</div><div> &END EACH</div><div> &END RESTART</div><div> &END PRINT</div><div> &END MOTION</div><div> &FORCE_EVAL</div><div> METHOD QS</div><div> &DFT</div><div> BASIS_SET_FILE_NAME BASIS_MOLOPT</div><div> POTENTIAL_FILE_NAME GTH_POTENTIALS</div><div> &SCF</div><div> MAX_SCF 300</div><div> EPS_SCF 1.0000000000000001E-05</div><div> SCF_GUESS RESTART</div><div> &OT T</div><div> PRECONDITIONER NONE</div><div> &END OT</div><div> &PRINT</div><div> &RESTART ON</div><div> ADD_LAST NUMERIC</div><div> BACKUP_COPIES 1</div><div> &EACH</div><div> MD 1</div><div> &END EACH</div><div> &END RESTART</div><div> &END PRINT</div><div> &END SCF</div><div> &QS</div><div> EPS_DEFAULT 1.0000000000000000E-10</div><div> METHOD GPW</div><div> &END QS</div><div> &MGRID</div><div> NGRIDS 4</div><div> CUTOFF 4.0000000000000000E+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_FUNCTIONAL NO_SHORTCUT</div><div> &PBE T</div><div> PARAMETRIZATION ORIG</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.5000000000000005E+01</div><div> TYPE DFTD3</div><div> PARAMETER_FILE_NAME dftd3.dat</div><div> REFERENCE_FUNCTIONAL PBE</div><div> CALCULATE_C9_TERM T</div><div> REFERENCE_C9_TERM T</div><div> &END PAIR_POTENTIAL</div><div> &END VDW_POTENTIAL</div><div> &END XC</div><div> &POISSON</div><div> PERIODIC XYZ</div><div> &END POISSON</div><div> &PRINT</div><div> &E_DENSITY_CUBE SILENT</div><div> ADD_LAST NUMERIC</div><div> STRIDE 1 1 1</div><div> &EACH</div><div> MD 99999999</div><div> &END EACH</div><div> &END E_DENSITY_CUBE</div><div> &PDOS SILENT</div><div> ADD_LAST NUMERIC</div><div> FILENAME dosfile</div><div> LOG_PRINT_KEY T</div><div> COMPONENTS F</div><div> NLUMO -1</div><div> &EACH</div><div> MD 99999999</div><div> &END EACH</div><div> &END PDOS</div><div> &END PRINT</div><div> &END DFT</div></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br><br>On Thursday, August 18, 2016 at 3:05:57 PM UTC+3, jgh wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Hi
<br>
<br>The PBE-D3 setup looks fine.
<br>
<br>Yes, 10-20 SCF iterations during an MD are on the very high end for
<br>such a well behaved system. I would try the following setup (in case you
<br>don't use it already)
<br>
<br>&QS
<br> ....
<br> EXTRAPOLATION ASPC
<br> EXTRAPOLATION_ORDER 4
<br>&END QS
<br>
<br>&OT
<br> PRECONDITIONER FULL_SINGLE_INVERSE
<br> MINIMIZER DIIS
<br>&END OT
<br>
<br>regards
<br>
<br>Juerg
<br>------------------------------<wbr>------------------------------<wbr>--
<br>Juerg Hutter <wbr> Phone : ++41 44 635 4491
<br>Institut für Chemie C FAX : ++41 44 635 6838
<br>Universität Zürich E-<wbr>mail: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="mMn6c0BwCgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">hut...@chem.uzh.ch</a>
<br>Winterthurerstrasse 190
<br>CH-8057 Zürich, Switzerland
<br>------------------------------<wbr>------------------------------<wbr>---
<br>
<br>-----<a href="javascript:" target="_blank" gdf-obfuscated-mailto="mMn6c0BwCgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a> wrote: -----To: cp2k <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="mMn6c0BwCgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>>
<br>From: Simiam Ghan
<br>Sent by: <a href="javascript:" target="_blank" gdf-obfuscated-mailto="mMn6c0BwCgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>
<br>Date: 08/18/2016 01:08AM
<br>Subject: Re: [CP2K:8075] SCF, MD run-time verses atomic species.
<br>
<br>Hello Juerg,Thank you for your reply. I added REFERENCE_C9_TERM .TRUE. to the vdw section as you suggested and found that my MD simulation became twice (!) as fast. That is, time per MD step went from ~60 sec to ~30 sec for box of water with 64 molecules and K,Cl ions. I made a quick comparison of Energy_Force results with and without this flag and indeed forces and energies do not change significantly. This is great news. Why is this not a default? What's the catch? Also, would you say this is now a correct way to declare a PBE-D3 setup?
<br>
<br>&VDW_POTENTIAL POTENTIAL_TYPE PAIR_POTENTIAL &PAIR_POTENTIAL R_CUTOFF 1.5000000000000005E+01 TYPE DFTD3 PARAMETER_FILE_NAME dftd3.dat REFERENCE_FUNCTIONAL PBE CALCULATE_C9_TERM T REFERENCE_C9_TERM T &END PAIR_POTENTIAL &END VDW_POTENTIAL
<br>Finally, I mentioned earlier that 30 scf steps were needed to converge the same box of water with ions. I was then using EPS_SCF = 1E-06. Reducing this to EPS_SCF = 1E-05 brought the number down to 10-20 scf necessary. Does that still sound high? Is the definition of EPS_SCF documented somewhere, as it is apparently not just Energy convergence.
<br>Cheers,Simiam Ghan
<br>
<br>
<br>
<br>On Friday, July 1, 2016 at 3:10:52 PM UTC+3, jgh wrote:Hi
<br>
<br>
<br>
<br>no this shouldn't be, but without more information I will have to guess.
<br>
<br>You could also have a look at the timings at the end of the output to
<br>
<br>see if some routines got slower or if all parts of the run were affected.
<br>
<br>
<br>
<br>Two things to consider:
<br>
<br>
<br>
<br>1) Use REFERENCE_C9_TERM TRUE in order to reduce the time for vdW in MD.
<br>
<br>2) 30 SCF iterations in MD for such a simple system is pointing to a problem
<br>
<br> with your setup.
<br>
<br>
<br>
<br>regards
<br>
<br>
<br>
<br>Juerg
<br>
<br>------------------------------<wbr>------------------------------<wbr>--
<br>
<br>Juerg Hutter <wbr> Phone : ++41 44 635 4491
<br>
<br>Institut für Chemie C FAX : ++41 44 635 6838
<br>
<br>Universität Zürich E-<wbr>mail: <a>hut...@chem.uzh.ch</a>
<br>
<br>Winterthurerstrasse 190
<br>
<br>CH-8057 Zürich, Switzerland
<br>
<br>------------------------------<wbr>------------------------------<wbr>---
<br>
<br>
<br>
<br>-----<a>cp...@googlegroups.com</a> wrote: -----To: cp2k <<a>cp...@googlegroups.com</a>>
<br>
<br>From: Simiam Ghan
<br>
<br>Sent by: <a>cp...@googlegroups.com</a>
<br>
<br>Date: 07/01/2016 01:38PM
<br>
<br>Subject: [CP2K:7883] SCF, MD run-time verses atomic species.
<br>
<br>
<br>
<br>Dear all,I am running NVT MD with Quickstep on a box of water with 64 molecules. If I replace a water molecule with a KCl ion pair, i observe that the MD and SCF step times more than double on my setup. (MD from ~20 sec to ~46 sec). SCF iterations are converging (except the very first MD step) in around 30 steps in both cases but each SCF now takes over twice as long as before. Is there an explanation of why such a 'small' change in system could double the run time? Is KCl really so heavy to calculate compared to H2O? The number of electrons in the system increases by 8, from 512 to 520. Below my KCl input file.
<br>
<br>Greetings,Simiam
<br>
<br>
<br>
<br>
<br>
<br>&GLOBAL
<br>
<br> PROJECT H2O_KCl RUN_TYPE MD PRINT_LEVEL MEDIUM
<br>
<br>&END GLOBAL
<br>
<br>&FORCE_EVAL
<br>
<br> METHOD Quickstep ! GPW method.
<br>
<br> &SUBSYS ! A subsystem: coordinates, topology, molecules and cell.
<br>
<br> &CELL ! Supercell setup. ABC [angstrom] 12.414 12.414 12.414 ! Using 64 H2O molecules, we thus get a density of 1g/cm^3. PERIODIC XYZ ! Use PBC in all dimensions. &END CELL
<br>
<br> &COORD UNIT angstromH -0.567712 -0.469646 -0.645913H 0.626116 -0.687796 0.308193O 0 0 0(...)K 2.1035 2.1035 4.2735Cl 4.1035 4.1035 2.6035 ###H 3.73881 3.10388 5.46104
<br>
<br>H 3.65742 2.89924 6.989O 3.1035 3.1035 6.207(...) &END COORD
<br>
<br> &KIND O BASIS_SET DZVP-MOLOPT-GTH-q6 POTENTIAL GTH-PBE-q6 &END KIND &KIND H BASIS_SET DZVP-MOLOPT-GTH-q1 POTENTIAL GTH-PBE-q1 &END KIND &KIND K BASIS_SET DZVP-MOLOPT-SR-GTH-q9 POTENTIAL GTH-PBE-q9 &END KIND &KIND Cl BASIS_SET DZVP-MOLOPT-GTH-q7 POTENTIAL GTH-PBE-q7 &END KIND
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br> &END SUBSYS
<br>
<br> &DFT
<br>
<br> BASIS_SET_FILE_NAME BASIS_MOLOPT POTENTIAL_FILE_NAME GTH_POTENTIALS ! SPIN_POLARIZED ! Do spin-polarized calculation
<br>
<br> &POISSON PERIODIC XYZ &END POISSON
<br>
<br> &QS METHOD GPW EPS_DEFAULT 1.0E-10 ! Set various epsilons for QS to values that will lead ! to energy correct up to 1e-10. &END QS
<br>
<br> &MGRID CUTOFF 400 ! This is Ecut of eq. 39 in VandeVondele (2005), i.e., plane-wave cutoff ! that determines size of finest grid (see caption of Fig. 1). Cutoffs for ! the subsequent, coarser grid levels are given by eq. 39. NGRIDS 4 ! This is N of eq. 39 in VandeVondele (2005), i.e., number of grids used. REL_CUTOFF 40 ! This controls the grid level onto which Gaussians will be mapped. &END MGRID
<br>
<br> &XC
<br>
<br> &XC_FUNCTIONAL &PBE PARAMETRIZATION ORIG &END PBE &END XC_FUNCTIONAL
<br>
<br> &VDW_POTENTIAL
<br>
<br> POTENTIAL_TYPE PAIR_POTENTIAL
<br>
<br> &PAIR_POTENTIAL TYPE DFTD3 REFERENCE_FUNCTIONAL PBE CALCULATE_C9_TERM .TRUE. PARAMETER_FILE_NAME dftd3.dat R_CUTOFF 15.0 &END PAIR_POTENTIAL
<br>
<br> &END VDW_POTENTIAL
<br>
<br> &END XC
<br>
<br> &SCF
<br>
<br> SCF_GUESS RESTART ! Use data from previous run as initial guess for wavefunction. EPS_SCF 1.0E-6 ! Threshold for converged total energy. MAX_SCF 300 ! Maximum number of SCF iterations performed.
<br>
<br> &OT PRECONDITIONER NONE ! This should be stable with respect to the "Cholesky errors" &END OT
<br>
<br> &PRINT &RESTART ON
<br>
<br> BACKUP_COPIES 1
<br>
<br> &EACH MD 1 &END EACH
<br>
<br> ADD_LAST NUMERIC
<br>
<br> &END RESTART
<br>
<br> &END PRINT
<br>
<br> &END SCF
<br>
<br> &PRINT &E_DENSITY_CUBE
<br>
<br> STRIDE 1 1 1
<br>
<br> &EACH MD 99999999 &END EACH
<br>
<br> ADD_LAST NUMERIC &END E_DENSITY_CUBE
<br>
<br> &PDOS COMPONENTS .FALSE. NLUMO = -1 FILENAME dosfile LOG_PRINT_KEY TRUE
<br>
<br> &EACH MD 99999999 &END EACH
<br>
<br> ADD_LAST NUMERIC &END PDOS
<br>
<br> &END PRINT
<br>
<br> &END DFT
<br>
<br>&END FORCE_EVAL
<br>
<br>&MOTION
<br>
<br> &MD ENSEMBLE NVT STEPS 10000 TEMPERATURE 300.0 ! K TIMESTEP 0.5 ! fs &THERMOSTAT
<br>
<br> REGION GLOBAL TYPE NOSE
<br>
<br> &NOSE LENGTH 3 ! Length of Nose-Hoover chain TIMECON 20.0 ! Period of typical vibrational motion in system in fs &END NOSE &END THERMOSTAT &END MD
<br>
<br> &PRINT
<br>
<br> &RESTART
<br>
<br> &EACH MD 1 &END EACH
<br>
<br> ADD_LAST NUMERIC
<br>
<br> &END RESTART
<br>
<br> &TRAJECTORY ON ADD_LAST NUMERIC FILENAME trajectory &END TRAJECTORY
<br>
<br> &END PRINT
<br>
<br>&END MOTION
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>--
<br>
<br>
<br>
<br>You received this message because you are subscribed to the Google Groups "cp2k" group.
<br>
<br>
<br>
<br>To unsubscribe from this group and stop receiving emails from it, send an email to <a>cp2k+...@googlegroups.com</a>.
<br>
<br>
<br>
<br>To post to this group, send email to <a>cp...@googlegroups.com</a>.
<br>
<br>
<br>
<br>Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>.
<br>
<br>
<br>
<br>For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>--
<br>
<br>You received this message because you are subscribed to the Google Groups "cp2k" group.
<br>
<br>To unsubscribe from this group and stop receiving emails from it, send an email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="mMn6c0BwCgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+uns...@googlegroups.<wbr>com</a>.
<br>
<br>To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="mMn6c0BwCgAJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.
<br>
<br>Visit this group at <a href="https://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/group/cp2k';return true;" onclick="this.href='https://groups.google.com/group/cp2k';return true;">https://groups.google.com/<wbr>group/cp2k</a>.
<br>
<br>For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.
<br>
<br>
<br></blockquote></div></div></div>