[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