Constraint
rupa
shyam... at googlemail.com
Thu Feb 25 15:59:49 UTC 2010
Hi all,
I have some problem with constrained MD of water-dimer with O--O
distance as constraint, which is required to extent up to 5 Å.
However, the required criteria was not satisfied in the simulation,
the extension of O--O distance was beyond 5 Å.
Any suggestion in this regard is appreciated.
Thanks
Rupa
P.S. Input file is provided below for reference.
&FORCE_EVAL
METHOD Quickstep
&DFT
BASIS_SET_FILE_NAME GTH_BASIS_SETS
POTENTIAL_FILE_NAME POTENTIAL
&MGRID
CUTOFF 280
NGRIDS 1
REL_CUTOFF 50
&END MGRID
&QS
EPS_DEFAULT 1.0E-12
EPS_GVG 1.0E-6
EPS_PGF_ORB 1.0E-6
&END QS
&SCF
EPS_SCF 1.0E-6
MAX_SCF 350
SCF_GUESS atomic
&END SCF
&XC
&XC_FUNCTIONAL BLYP
&END XC_FUNCTIONAL
&VDW_POTENTIAL
POTENTIAL_TYPE PAIR_POTENTIAL
&PAIR_POTENTIAL
&END PAIR_POTENTIAL
&END VDW_POTENTIAL
&END XC
&POISSON
PERIODIC NONE
POISSON_SOLVER MT
&END POISSON
&END DFT
&SUBSYS
&CELL
ABC 20.000 20.000 20.000
PERIODIC NONE
&END CELL
&COORD
H 0.828857203 1.510765153 1.850758135
O 0.223273347 1.612802255 1.087529046
H 0.746091252 2.082274953 0.419190658
O 1.531077052 1.395421022 3.663417772
H 1.345748937 0.482975975 3.940795058
H 0.806283419 1.920091495 4.042933787
&END COORD
&KIND O
BASIS_SET TZVP-GTH
POTENTIAL GTH-BLYP-q6
&END KIND
&KIND H
BASIS_SET TZVP-GTH
POTENTIAL GTH-BLYP-q1
&END KIND
&COLVAR
&DISTANCE
ATOMS 2 4
&END DISTANCE
&END COLVAR
&END SUBSYS
&END FORCE_EVAL
&GLOBAL
PROJECT md
RUN_TYPE MD
PRINT_LEVEL LOW
&END GLOBAL
&MOTION
&MD
ENSEMBLE NVT
STEPS 4000
TIMESTEP 1.0
TEMPERATURE 300.0
&THERMOSTAT
&NOSE
LENGTH 3
YOSHIDA 3
TIMECON 100.
MTS 2
&END NOSE
&END THERMOSTAT
&END MD
&CONSTRAINT
SHAKE_TOLERANCE 1.0E-08
&COLLECTIVE
COLVAR 1
INTERMOLECULAR
TARGET 5
TARGET_GROWTH 0.0001
&END COLLECTIVE
&END CONSTRAINT
&END MOTION
More information about the CP2K-user
mailing list