BOMD - Total energy increasing
Sandeep Kumar Reddy
kuma... at gmail.com
Fri Jul 1 08:42:40 UTC 2011
Dear all,
I am doing an BOMD run of carboxylic acid in water. After 6ps,
i plotted conserved quantity as a function of step no. which shows total
energy is increasing. Can anybody help
why it is increasing ?
Input file :
&FORCE_EVAL
METHOD Quickstep
&DFT
BASIS_SET_FILE_NAME BASIS_MOLOPT
POTENTIAL_FILE_NAME POTENTIAL
CHARGE = 0
&MGRID
CUTOFF 280
NGRIDS 5
REL_CUTOFF 40
&END MGRID
&QS
METHOD GPW
EPS_DEFAULT 1.0E-10
&END QS
&SCF
SCF_GUESS ATOMIC
EPS_SCF 1.0E-6
MAX_SCF 150
# MIXING 0.3
&OUTER_SCF
EPS_SCF 1.0E-6
MAX_SCF 50
&END
&OT
MINIMIZER CG
PRECONDITIONER FULL_ALL
ENERGY_GAP 0.001
STEPSIZE 0.05
&END
&END SCF
&XC
&XC_FUNCTIONAL BLYP
&END XC_FUNCTIONAL
&vdW_POTENTIAL
DISPERSION_FUNCTIONAL PAIR_POTENTIAL
&PAIR_POTENTIAL
TYPE DFTD2
REFERENCE_FUNCTIONAL BLYP
&END PAIR_POTENTIAL
&END vdW_POTENTIAL
&END XC
&END DFT
&SUBSYS
&CELL
ABC 12.16174 12.16174 12.16174
&END CELL
&COORD
@INCLUDE coor.xyz
&END COORD
&KIND O
BASIS_SET TZV2P-MOLOPT-GTH-q6
POTENTIAL GTH-BLYP-q6
&END KIND
&KIND C
BASIS_SET TZV2P-MOLOPT-GTH-q4
POTENTIAL GTH-BLYP-q4
&END KIND
&KIND H
BASIS_SET TZV2P-MOLOPT-GTH-q1
POTENTIAL GTH-BLYP-q1
MASS 1.00794
&END KIND
&END SUBSYS
&END FORCE_EVAL
&GLOBAL
PROJECT OUT-URACIL-MD
RUN_TYPE MD
PRINT_LEVEL MEDIUM
&END GLOBAL
&MOTION
&MD
ENSEMBLE NVT
STEPS 100000
TIMESTEP 0.5
&THERMOSTAT
&NOSE
LENGTH 3
YOSHIDA 3
TIMECON 100.0
MTS 2
&END NOSE
&END THERMOSTAT
TEMPERATURE 298.15
&END MD
&END MOTION
----
Thanks & Regards,
Sandeep kumar Reddy
Ph. D. student
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20110701/b7cab70f/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Untitled.png
Type: image/png
Size: 8296 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20110701/b7cab70f/attachment.png>
More information about the CP2K-user
mailing list