/usr/share/doc/cp2k/examples/H2O-qmmm-gauss-10.inp is in cp2k 2.4.0-2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 | &FORCE_EVAL
METHOD QMMM
&DFT
BASIS_SET_FILE_NAME ../GTH_BASIS_SETS
POTENTIAL_FILE_NAME ../POTENTIAL
&MGRID
COMMENSURATE
CUTOFF 20
&END MGRID
&QS
&END QS
&SCF
SCF_GUESS atomic
&END SCF
&XC
&XC_FUNCTIONAL pade
&END XC_FUNCTIONAL
&END XC
&END DFT
&MM
&FORCEFIELD
&BEND
ATOMS H O H
K 0.
THETA0 1.8
&END BEND
&BOND
ATOMS O H
K 0.
R0 1.8
&END BOND
&CHARGE
ATOM O
CHARGE -0.8476
&END CHARGE
&CHARGE
ATOM H
CHARGE 0.4238
&END CHARGE
&NONBONDED
&LENNARD-JONES
atoms O O
EPSILON 78.198
SIGMA 3.166
RCUT 11.4
&END LENNARD-JONES
&LENNARD-JONES
atoms O H
EPSILON 0.0
SIGMA 3.6705
RCUT 11.4
&END LENNARD-JONES
&LENNARD-JONES
atoms H H
EPSILON 0.0
SIGMA 3.30523
RCUT 11.4
&END LENNARD-JONES
&END NONBONDED
&END FORCEFIELD
&POISSON
&EWALD
EWALD_TYPE ewald
ALPHA .44
GMAX 21
&END EWALD
&END POISSON
&END MM
&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
#
# QM_KINDS
#
#
# MM_KINDS
#
#
&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
&END FORCE_EVAL
&GLOBAL
PROJECT H2O-qmmm-gauss-10
PRINT_LEVEL MEDIUM
RUN_TYPE MD
&END GLOBAL
&MOTION
&MD
ENSEMBLE NVE
STEPS 1
TIMESTEP 0.5
TEMPERATURE 298
&END MD
&END MOTION
|