1 |
<OpenMD version=1> |
2 |
<MetaData> |
3 |
|
4 |
molecule { |
5 |
name = "TIP4P_FQ"; |
6 |
|
7 |
atom[0]{ |
8 |
type = "O_TIP4P_FQ"; |
9 |
position( 0.0, 0.0, -0.06556 ); |
10 |
} |
11 |
atom[1]{ |
12 |
type = "H_TIP4P_FQ"; |
13 |
position( 0.0, 0.75695, 0.52032 ); |
14 |
} |
15 |
atom[2]{ |
16 |
type = "H_TIP4P_FQ"; |
17 |
position( 0.0, -0.75695, 0.52032 ); |
18 |
} |
19 |
atom[3]{ |
20 |
type = "M_TIP4P_FQ"; |
21 |
position( 0.0, 0.0, 0.08444 ); |
22 |
} |
23 |
|
24 |
rigidBody[0]{ |
25 |
members(0,1,2,3); |
26 |
} |
27 |
} |
28 |
component{ |
29 |
type = "TIP4P_FQ"; |
30 |
nMol = 2; |
31 |
} |
32 |
|
33 |
|
34 |
ensemble = NVE; |
35 |
forceField = "FlucQ"; |
36 |
forceFieldFileName = "FQ.frc"; |
37 |
cutoffMethod = "shifted_force"; |
38 |
dielectric = 80.0; |
39 |
cutoffRadius = 9.0; |
40 |
switchingRadius = 7.7; |
41 |
outputFluctuatingCharges = true; |
42 |
runTime = 5; |
43 |
sampleTime = 1; |
44 |
statusTime = 1; |
45 |
|
46 |
targetTemp = 300; |
47 |
targetPressure = 1.0; |
48 |
|
49 |
tauThermostat = 1e3; |
50 |
tauBarostat = 1e4; |
51 |
|
52 |
dt = 1.0; |
53 |
useInitialTime = "false"; |
54 |
useInitialExtendedSystemState = "false"; |
55 |
|
56 |
//tempSet = "true"; |
57 |
//thermalTime = 10; |
58 |
</MetaData> |
59 |
<Snapshot> |
60 |
<FrameData> |
61 |
Time: 10000 |
62 |
Hmat: {{ 19.713, 0, 0 }, { 0, 19.713, 0 }, { 0, 0, 19.713 }} |
63 |
Thermostat: 0 , 0 |
64 |
Barostat: {{ 0, 0, 0 }, { 0, 0, 0 }, { 0, 0, 0 }} |
65 |
</FrameData> |
66 |
<StuntDoubles> |
67 |
0 pvqjcw 23.95523961 -2.701653649 -15.70063919 2.409292e-03 9.908152e-06 -1.574237e-03 2.074455e-01 5.294415e-01 5.318552e-01 6.275253e-01 -6.757743e-04 -5.046488e-03 1.001059e-02 0.000000e+00 0.000000e+00 |
68 |
1 pvqjcw -18.95523961 2.701653649 15.70063919 -2.409292e-03 -9.908152e-06 1.574237e-03 1.748160e-03 5.106060e-01 -8.393383e-01 -1.865201e-01 5.628458e-03 -8.785288e-03 -1.926888e-02 0.000000e+00 0.000000e+00 |
69 |
</StuntDoubles> |
70 |
</Snapshot> |
71 |
</OpenMD> |