How to set the positions of QMMM cell explicitly in QMMM calculation
Geng Sun
sungen... at gmail.com
Wed Jan 13 09:15:22 UTC 2016
Hello everyone,
I have a question about how to set the QMMM calculation,
Under the example folder from the cp2k distributions, one of the
input file looks like below : (only the QMMM and SUBSYS are shown)
In the input file, we have two CELL sections in QMMM and SUBSYS. It
seems that the cell in QMMM is encompassed by cell in SUBSYS
Here is my question, what's the relative positions of cell/QMMM in
cell/SUBSYS ?
is cell/QMMM movable in MD simulations ? what's the meaning of the
keywors INITIAL_TRANSLATION_VECTOR
<https://manual.cp2k.org/cp2k-3_0-branch/CP2K_INPUT/FORCE_EVAL/QMMM.html#INITIAL_TRANSLATION_VECTOR>
?
Thanks
Geng Sun
&QMMM
MM_POTENTIAL_FILE_NAME ../MM_POTENTIAL
&CELL
ABC 6.0 6.0 6.0
&END CELL
ECOUPL GAUSS
NOCOMPATIBILITY
&INTERPOLATOR
EPS_R 1.0e-14
EPS_X 1.0e-14
MAXITER 100
&END INTERPOLATOR
&MM_KIND H
RADIUS 0.44
&END MM_KIND
&MM_KIND O
RADIUS 0.78
&END MM_KIND
&QM_KIND H
MM_INDEX 2
MM_INDEX 3
&END QM_KIND
&QM_KIND O
MM_INDEX 1
&END QM_KIND
&END QMMM
&SUBSYS
&CELL
ABC 24.955 24.955 24.955
&END CELL
&COORD
O 0.000000 0.000000 0.000000 H2O1
H 0.000000 0.000000 1.000000 H2O1
H 0.942809 0.000000 -0.333333 H2O1
O -1.617979 -0.948062 -2.341650 H2O2
H -2.529195 -1.296822 -2.122437 H2O2
H -1.534288 -0.833088 -3.331486 H2O2
O -1.447990 2.117783 1.555094 H2O3
H -1.501128 2.645178 2.403050 H2O3
H -2.090603 1.352766 1.597519 H2O3
&END COORD
&KIND H
BASIS_SET SZV-GTH
POTENTIAL GTH-PADE-q1
&END KIND
&KIND O
BASIS_SET SZV-GTH
POTENTIAL GTH-PADE-q6
&END KIND
&TOPOLOGY
&END TOPOLOGY
&END SUBSYS
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20160113/13e9d454/attachment.htm>
More information about the CP2K-user
mailing list