[CP2K-user] [CP2K:20987] METADYNAMICS run to move atoms apart

Vallabh Vasudevan vallav at umich.edu
Wed Dec 18 22:04:32 UTC 2024


Hi,

I am setting up a metadynamics simulation to move a Hydrogen away from a 
specific Oxygen. I am planning on setting up a distance based collective 
variable. The DISTANCE function under COLVAR does not give a final value to 
move towards. How do I do this?

&MOTION
&MD
        @INCLUDE            ${intg_file}

         STEPS              1000
         TIMESTEP           0.5
         TEMPERATURE        300

&END MD
&FREE_ENERGY
&METADYN
      DO_HILLS .TRUE. ! unbiased MTD run
      NT_HILLS 60 ! 
      WW 0.5E-3 ! height of the gaussian to spawn - default = 0.1 hartree

 

      &METAVAR ! for each COLVAR, a METAVAR section has to be set
        COLVAR 1
        SCALE  1 ! change only when do_hills true
&END METAVAR
&PRINT
&COLVAR
           COMMON_ITERATION_LEVELS 3
&EACH
             MD 1
&END EACH
&END COLVAR

 

        &HILLS
           COMMON_ITERATION_LEVELS 3
&EACH
             MD 1
&END EACH
&END HILLS
&END PRINT
&END METADYN
&END FREE_ENERGY
&PRINT
        @INCLUDE            ${mdio_file}
&END
&END MOTION

 

 

&FORCE_EVAL
   METHOD QS
   STRESS_TENSOR ANALYTICAL
   @INCLUDE                 ${qmps_file}

&SUBSYS
      @INCLUDE              ${qmbs_file}
&CELL
         @INCLUDE           ${cell_file}
         CELL_FILE_FORMAT   cp2k
         PERIODIC XYZ
&END

&TOPOLOGY
        COORD_FILE_NAME     ${crds_file}
        COORD_FILE_FORMAT   xyz
&END TOPOLOGY

 

      &COLVAR
&DISTANCE
          ATOMS 0 15 ! Select the middle O and H to pull
          AXIS XYZ 
&END DISTANCE
&END COLVAR
&END SUBSYS

 

&END FORCE_EVAL
has context menu

-- 
You received this message because you are subscribed to the Google Groups "cp2k" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cp2k+unsubscribe at googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/cp2k/8115c8cc-cdb9-492b-a8d8-eb38c4eded67n%40googlegroups.com.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20241218/15a7868c/attachment.htm>


More information about the CP2K-user mailing list