Hello all:<div>I'm a newbie from the classical MD world. Recently, I've been attempting to run QS dynamics on the rock salt structure of LiF at 300K to get my feet wet. My simulations are on the 2x2x2 supercell of the cubic crystal structure. The simulation runs and the thermodynamic properties (temp, energy etc) are all converged after an initial equilibration, however the system melts, transforming from the rock salt structure into an amorphous blend after only 1ps or so of dynamics. Since this structure has been simulated before with QS, I must be doing something wrong. Can anybody see anything obviously wrong with the following input?</div><div><br></div><div>Best,</div><div>Mike</div><div><br></div><div>============================================================================================</div><div><div>@SET MD_SAVE_FREQ  1</div><div>@SET LA 8.2</div><div>@SET LB 8.2</div><div>@SET LC 8.2</div><div>@SET RTEMP 298</div><div>@SET PREFIX LiF.2x2x2</div><div>@SET NSTEPS 10000</div><div><br></div><div>&GLOBAL</div><div> PRINT_LEVEL medium</div><div> PROJECT ${PREFIX}.${RTEMP}</div><div> RUN_TYPE MD</div><div> WALLTIME 36000000</div><div>&END GLOBAL</div><div><br></div><div>&MOTION</div><div> &CELL_OPT</div><div>  EXTERNAL_PRESSURE 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0</div><div>  KEEP_ANGLES TRUE</div><div>  OPTIMIZER BFGS</div><div>  TYPE MD</div><div>  &PRINT</div><div>   &CELL</div><div>    FILENAME ./${PREFIX}.${RTEMP}K</div><div>    &EACH</div><div>     CELL_OPT ${MD_SAVE_FREQ}</div></div><div><div>    &END EACH</div><div>   &END CELL</div><div>  &END PRINT</div><div> &END CELL_OPT</div><div><br></div><div> &PRINT</div><div>  &VELOCITIES</div><div>   FILENAME ./${PREFIX}.${RTEMP}K</div><div>   FORMAT XYZ</div><div>   &EACH</div><div>     MD ${MD_SAVE_FREQ}</div><div>   &END EACH</div><div>  &END VELOCITIES</div><div>  &TRAJECTORY</div><div>   FILENAME ./${PREFIX}.${RTEMP}K</div><div>   FORMAT XYZ</div><div>   &EACH</div><div>     MD ${MD_SAVE_FREQ}</div><div>   &END EACH</div><div>  &END TRAJECTORY</div><div> &END PRINT</div></div><div><div> &MD</div><div>  ENSEMBLE NVT</div><div>  STEPS ${NSTEPS}</div><div>  TIMESTEP 0.5</div><div>  TEMPERATURE ${RTEMP}</div><div>  &LANGEVIN</div><div>   GAMMA 0.0001</div><div>  &END LANGEVIN</div><div>  &THERMOSTAT</div><div>   REGION MOLECULE</div><div>   TYPE NOSE</div><div>   &NOSE</div><div>    LENGTH 3</div><div>    YOSHIDA 3</div><div>    TIMECON 100</div><div>   &END NOSE</div><div>  &END THERMOSTAT</div><div> &END MD</div><div>&END MOTION</div><div><br></div><div>&FORCE_EVAL</div><div> METHOD Quickstep</div><div> &DFT</div><div>  CHARGE 0</div><div>  MULTIPLICITY 1</div></div><div><div>  WFN_RESTART_FILE_NAME ${PREFIX}.${RTEMP}.md.wfn.restart</div><div>  BASIS_SET_FILE_NAME /global/home/users/mike/codes/cp2k/tests/QS/BASIS_SET</div><div>  POTENTIAL_FILE_NAME /global/home/users/mike/codes/cp2k/tests/QS/GTH_POTENTIALS</div><div>  &MGRID</div><div>   CUTOFF [Ry] 320</div><div>  &END MGRID</div><div>  &PRINT</div><div>   &MULLIKEN</div><div>    FILENAME ./${PREFIX}.${RTEMP}</div><div>    &EACH</div><div>     MD 1</div><div>    &END EACH</div><div>   &END MULLIKEN</div><div>  &END PRINT</div><div><br></div><div>  &QS</div><div>   EPS_DEFAULT 1.0E-12</div><div>   EXTRAPOLATION ASPC</div><div>   EXTRAPOLATION_ORDER 3</div><div>   MAP_CONSISTENT</div><div>   METHOD GPW</div><div>  &END QS</div><div>  &SCF</div><div>   EPS_SCF 1.0E-5</div><div>   MAX_SCF 50</div><div>   SCF_GUESS ATOMIC</div></div><div><div>   &DIAGONALIZATION</div><div>    ALGORITHM DAVIDSON</div><div>   &END DIAGONALIZATION</div><div><br></div><div>   &OUTER_SCF</div><div>    EPS_SCF 1.0E-5</div><div>    MAX_SCF 20</div><div>   &END OUTER_SCF</div><div><br></div><div>   #&OT</div><div>   # MINIMIZER DIIS</div><div>   # PRECONDITIONER FULL_ALL</div><div>   #&END OT</div><div>  &END SCF</div><div>  &XC</div><div>   &XC_FUNCTIONAL PADE</div><div>   &END XC_FUNCTIONAL</div></div><div><div>  &END XC</div><div><br></div><div>  &POISSON</div><div>   PERIODIC XYZ</div><div>   POISSON_SOLVER PERIODIC</div><div>  &END POISSON</div><div> &END DFT</div><div> &SUBSYS</div><div>  &CELL</div><div>   ABC ${LA} ${LB} ${LC}</div><div>   PERIODIC XYZ</div><div>  &END CELL</div><div><br></div><div>  &COORD</div><div>   @INCLUDE ${PREFIX}.cp2k.xyz</div><div>   SCALED .TRUE.</div><div>  &END COORD</div><div><br></div><div>  &KIND Li</div><div>   BASIS_SET DZV-GTH-PADE</div><div>   POTENTIAL GTH-PADE-q3</div><div>   #BASIS_SET DZVP-GTH-PADE</div><div>   #POTENTIAL GTH-PADE-q3</div><div>  &END KIND</div></div><div><div>  &KIND F</div><div>   BASIS_SET DZVP-GTH-PADE</div><div>   POTENTIAL GTH-PADE-q7</div><div>  &END KIND</div><div> &END SUBSYS</div><div>&END FORCE_EVAL</div></div><div><br></div>