My metadynamics is not running properly <br />Please check my input file<br /><br />&GLOBAL<br />  PROJECT 1a<br />  RUN_TYPE MD<br />  PRINT_LEVEL MEDIUM<br />&END GLOBAL<br /><br />#&EXT_RESTART<br />#RESTART_FILE_NAME 1a-1.restart<br />#&END EXT_RESTART<br /><br />&MOTION<br />  &MD<br />    ENSEMBLE NVT # ENSEMBLE NVE<br />    STEPS 200000<br />    TIMESTEP 1.0<br />    TEMPERATURE 295.0<br />    &THERMOSTAT<br />      &NOSE<br />        LENGTH 3 # 3 is the default<br />        MTS 2 # 2 is the default<br />        TIMECON 300.0<br />        YOSHIDA 3 # 3 is the default<br />      &END NOSE<br />    &END THERMOSTAT<br />  &END MD<br /> <br />  &FREE_ENERGY<br />    &METADYN<br />      DO_HILLS<br />      NT_HILLS 100<br />      WW 3.0e-3<br /><br /><br />      &METAVAR<br />        SCALE 0.3<br />        COLVAR 1<br /><br />        &WALL<br />          POSITION [angstrom] 4.0<br />          TYPE QUADRATIC<br />          &QUADRATIC<br />            DIRECTION WALL_PLUS<br />            K [kcalmol] 150<br />          &END QUADRATIC<br /><span style="white-space: pre;">        </span>&END WALL<br />      &END METAVAR<br /><br />      &PRINT<br />        &COLVAR<br />           COMMON_ITERATION_LEVELS 3<br />           &EACH<br />             MD 1<br />           &END<br />        &END<br />      &END<br />    &END METADYN<br />  &END<br /><br /><br />  &PRINT<br />    &RESTART_HISTORY<br />      &EACH<br />        MD 100<br />      &END<br />    &END RESTART_HISTORY<br />    &VELOCITIES<br />      &EACH<br />        MD 1<br />      &END<br />    &END<br />  # Normal restart file<br />  &RESTART<br />    &EACH<br />      MD 1<br />    &END<br />  &END RESTART<br />  &END PRINT<br />&END MOTION<br /><br />&FORCE_EVAL<br />  METHOD Quickstep              ! Electronic structure method (DFT,...)<br />  &DFT<br />    CHARGE = 0<br />    MULTIPLICITY = 2<br />    LSD<br />    BASIS_SET_FILE_NAME  BASIS_SET.txt<br />    POTENTIAL_FILE_NAME  POTENTIAL.txt<br />    # if WFN restart is required<br />    WFN_RESTART_FILE_NAME 1a-RESTART.wfn<br />    &MGRID<br />     #CUTOFF 1000<br />      CUTOFF 300<br />      NGRIDS 4<br />      REL_CUTOFF 60<br />    &END MGRID<br />    &QS<br />      METHOD PM6<br />      EPS_DEFAULT 1.0E-10 # 1.0E-10 is the default<br />    &END QS<br />   &SCF<br />     SCF_GUESS RESTART<br />     EPS_SCF 1.0E-5 # 1.0E-5 is the default<br />     MAX_SCF 600 # 50 is the default<br />     &OT<br />     MINIMIZER CG<br />     PRECONDITIONER FULL_ALL<br />     LINESEARCH 3PNT<br />     &END<br />     &OUTER_SCF<br />       EPS_SCF 1.0E-5 # 1.0E-5 is the default<br />       MAX_SCF 20 # Keep > 100 when OT is on<br />     &END<br />   &END SCF  <br />   &XC<br />     &XC_FUNCTIONAL BLYP<br />     &END XC_FUNCTIONAL<br />     &vdW_POTENTIAL<br />       DISPERSION_FUNCTIONAL PAIR_POTENTIAL # POTENTIAL_TYPE is alias to DISPERSION_FUNCTIONAL<br />       &PAIR_POTENTIAL<br />         TYPE DFTD3<br />         PARAMETER_FILE_NAME dftd3.dat<br />         REFERENCE_FUNCTIONAL PBE<br />       &END PAIR_POTENTIAL<br />    &END vdW_POTENTIAL<br />   &END XC<br />   &POISSON                    ! Solver requested for non periodic calculations<br />      PERIODIC XYZ          ! Type of solver<br />   &END POISSON<br />  &END DFT  <br />  &SUBSYS<br />    &CELL<br /> <span style="white-space: pre;">     </span>A     20.0000000     0.00000000     0.00000000<br /> <span style="white-space: pre;"> </span>B     0.00000000     20.0000000     0.00000000<br /> <span style="white-space: pre;"> </span>C     0.00000000     0.00000000     20.00000000<br /><br />      PERIODIC XYZ              ! Non periodic calculations. That's why the POISSON section is needed <br />    &END CELL<br />    &TOPOLOGY                    ! Section used to center the atomic coordinates in the given box. Useful for big molecules<br />      COORD_FILE_FORMAT xyz<br />      COORD_FILE_NAME  1a.xyz<br />    &END<br />    &KIND H<br />      ELEMENT H<br />      BASIS_SET DZVP-GTH-BLYP<br />      POTENTIAL GTH-BLYP-q1<br />    &END KIND    <br />    &KIND O<br />      ELEMENT O<br />      BASIS_SET DZVP-GTH-BLYP<br />      POTENTIAL GTH-BLYP-q6<br />    &END KIND    <br />    &KIND C<br />      ELEMENT C<br />      BASIS_SET DZVP-GTH-BLYP<br />      POTENTIAL GTH-BLYP-q4<br />    &END KIND   <br />      &KIND Cl<br />      ELEMENT Cl<br />      BASIS_SET DZVP-GTH-BLYP<br />      POTENTIAL GTH-BLYP-q7<br />    &END KIND  <br /><br />     &COLVAR<br />      &DISTANCE<br />        ATOMS 2 5<br />      &END DISTANCE<br />     &END COLVAR<br />  &END SUBSYS<br />&END FORCE_EVAL<br /><br /><br />

<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/d6bda0e3-9068-4e04-9802-2bd1ce3e5aa5n%40googlegroups.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/d6bda0e3-9068-4e04-9802-2bd1ce3e5aa5n%40googlegroups.com</a>.<br />