<div dir="ltr">I am run MD simulation of WO2(OH)2(H2O) with tetrahedral configuration. I have run it up to 90 ps. This is my input file：<div><div>&FORCE_EVAL</div><div>  METHOD Quickstep</div><div>  &DFT</div><div>    CHARGE -1</div><div>    MULTIPLICITY 1</div><div>    BASIS_SET_FILE_NAME GTH_BASIS_SETS</div><div> ####   BASIS_SET_FILE_NAME BASIS_MOLOPT</div><div>    POTENTIAL_FILE_NAME POTENTIAL</div><div><br></div><div>    &MGRID</div><div>      CUTOFF 280</div><div>      REL_CUTOFF 30</div><div>    &END MGRID</div><div>    &QS</div><div>      EPS_DEFAULT 1.0E-12</div><div>      WF_INTERPOLATION ASPC</div><div>      EXTRAPOLATION_ORDER 3 </div><div>    &END QS</div><div>    &SCF</div><div>      MAX_SCF 150</div><div>      SCF_GUESS ATOMIC</div><div>      &OT ON</div><div>        PRECONDITIONER FULL_KINETIC</div><div>        MINIMIZER DIIS</div><div>      &END OT</div><div>      &OUTER_SCF</div><div>        MAX_SCF 50</div><div>        EPS_SCF 1.0E-5</div><div>      &END OUTER_SCF</div><div> </div><div>      &PRINT</div><div>        &RESTART OFF</div><div>        &END</div><div>      &END</div><div><br></div><div>    # SCF_GUESS        RESTART</div><div>    # EPS_SCF      1.0E-7</div><div>    &END SCF</div><div>    &XC</div><div>      &XC_FUNCTIONAL BLYP</div><div>      &END XC_FUNCTIONAL</div><div>      &XC_GRID</div><div>        XC_DERIV SPLINE2_SMOOTH</div><div>        XC_SMOOTH_RHO NN10</div><div>      &END XC_GRID</div><div>    &END XC</div><div>  &END DFT</div><div><br></div><div>  &SUBSYS</div><div>    &CELL</div><div>      PERIODIC XYZ</div><div>      ABC 9.86 9.86 9.86</div><div>    &END CELL</div><div>    &TOPOLOGY</div><div>      COORD_FILE_NAME   WO2-2OH-2H2O-box-dH-H2O.xyz</div><div>      COORD_FILE_FORMAT xyz</div><div>    &END TOPOLOGY</div><div>   &KIND H</div><div>      BASIS_SET DZVP-GTH</div><div>      POTENTIAL GTH-BLYP-q1</div><div>    &END KIND</div><div>    &KIND O</div><div>      BASIS_SET DZVP-GTH</div><div>      POTENTIAL GTH-BLYP-q6</div><div>    &END KIND</div><div>    &KIND W</div><div>      BASIS_SET DZVP-GTH</div><div>      POTENTIAL GTH-BLYP-q14</div><div>    &END KIND</div><div> ###    &COLVAR</div><div> ###    &DISTANCE</div><div> ###      ATOMS 7 8</div><div> ###    &END DISTANCE</div><div> ###   &END COLVAR</div><div> ###   &COLVAR</div><div> ###    &ANGLE</div><div> ###      ATOMS 8 7 9</div><div> ###    &END ANGLE</div><div> ###   &END COLVAR</div><div> ###   &COLVAR</div><div> ###    &ANGLE</div><div> ###      ATOMS 1 7 8</div><div> ###    &END ANGLE</div><div> ###   &END COLVAR</div><div>  &END SUBSYS</div><div>&END FORCE_EVAL</div><div><br></div><div>&GLOBAL</div><div>  PROJECT WO2-2OH-2H2O-dH-H2O_MD</div><div>  RUN_TYPE MD</div><div>  PRINT_LEVEL LOW</div><div><br></div><div>  &TIMINGS</div><div>     THRESHOLD 0.000001</div><div>  &END</div><div>&END GLOBAL</div><div><br></div><div>&MOTION</div><div>  &MD</div><div>    ENSEMBLE NVT</div><div>    STEPS 100000</div><div>    TIMESTEP 0.5</div><div>    TEMPERATURE 300.0    </div><div>  ###  TEMP_TOL 50.0</div><div>    &THERMOSTAT</div><div>     TYPE NOSE</div><div>     &NOSE</div><div>      TIMECON 100</div><div>     &END NOSE</div><div>    &END THERMOSTAT</div><div>  &END MD</div><div> ## &CONSTRAINT</div><div> ##  &COLLECTIVE</div><div> ##   &RESTRAINT</div><div> ##     K [hartree] 0.5</div><div> ##   &END RESTRAINT</div><div> ##    COLVAR 1</div><div> ##    INTERMOLECULAR</div><div> ##    TARGET [angstrom] 1.89</div><div> ##  &END COLLECTIVE</div><div> ##  &COLLECTIVE</div><div> ##   &RESTRAINT</div><div> ##     K [hartree] 0.05</div><div> ##   &END RESTRAINT         </div><div> ##    COLVAR 2       </div><div> ##    INTERMOLECULAR </div><div> ##    TARGET [rad] 1.87    </div><div> ##  &END COLLECTIVE  </div><div> ##  &COLLECTIVE </div><div> ##   &RESTRAINT</div><div> ##     K [hartree] 0.05</div><div> ##   &END RESTRAINT     </div><div> ##    COLVAR 3       </div><div> ##    INTERMOLECULAR </div><div> ##    TARGET [rad] 2.04    </div><div> ##  &END COLLECTIVE </div><div>####   &LAGRANGE_MULTIPLIERS ON</div><div>####     COMMON_ITERATION_LEVELS  3 </div><div>####   &END   </div><div>##  &END CONSTRAINT</div><div>##  &PRINT                   </div><div>##   &TRAJECTORY             </div><div>##    &EACH                  </div><div>##      MD 100               </div><div>##    &END                   </div><div>##   &END                    </div><div>##   &VELOCITIES OFF         </div><div>##   &END                    </div><div>##   &FORCES OFF             </div><div>##   &END                    </div><div>##   &RESTART_HISTORY OFF    </div><div>##   &END                    </div><div>##   &RESTART                </div><div>##    &EACH                  </div><div>##      MD 10                </div><div>##    &END                   </div><div>##   &END                    </div><div>##  &END                     </div><div>&END MOTION</div><div><br></div><div>&EXT_RESTART</div><div>  RESTART_FILE_NAME WO2-2OH-2H2O-dH-H2O_MD-1.restart</div><div>&END EXT_RESTART</div><div><br></div><br>在 2015年8月10日星期一 UTC+8下午4:59:33，Samuel Andermatt写道：<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Please provide the information listed here: <a href="http://www.cp2k.org/howto:gethelp" target="_blank" rel="nofollow" onmousedown="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwww.cp2k.org%2Fhowto%3Agethelp\46sa\75D\46sntz\0751\46usg\75AFQjCNENrmojKUmuYSeSZtdcx9PT7pSCWg';return true;" onclick="this.href='http://www.google.com/url?q\75http%3A%2F%2Fwww.cp2k.org%2Fhowto%3Agethelp\46sa\75D\46sntz\0751\46usg\75AFQjCNENrmojKUmuYSeSZtdcx9PT7pSCWg';return true;">http://www.cp2k.org/howto:<wbr>gethelp</a><br></div></blockquote></div></div>