Download ReaxFF User Manual

Example 2.8: ffield-input file.
Force field identifier
Reactive MD-force field: Hydrocarbon parameters
39
! Number of general parameters
50.0000 !Overcoordination parameter
9.8407 !Overcoordination parameter
21.2839 !Valency angle conjugation parameter
3.0000 !Triple bond stabilisation parameter
6.5000 !Triple bond stabilisation parameter
1.0000 !Not used
0.9782 !Undercoordination parameter
1.0250 !Triple bond stabilisation parameter
General
6.3452 !Undercoordination parameter
11.6274 !Undercoordination parameter
0.0000 !Triple bond stabilization energy
0.0000 !Lower Taper-radius
10.0000 !Upper Taper-radius
2.8793 !Not used
33.8667 !Valency undercoordination
88.6186 !Valency angle/lone pair parameter
1.0563 !Valency angle
2.0384 !Valency angle parameter
6.1431 !Not used
7.5203 !Double bond/angle parameter
0.3989 !Double bond/angle parameter: overcoord
3.9954 !Double bond/angle parameter: overcoord
-2.4837 !Not used
4.7120 !Torsion/BO parameter
10.0000 !Torsion overcoordination
2.3170 !Torsion overcoordination
-1.2635 !Conjugation 0 (not used)
2.1645 !Conjugation
1.4553 !vdWaals shielding
0.1000 !Cutoff for bond order (*100)
2.8921 !Valency angle conjugation parameter
7.1783 !Overcoordination parameter
1.4473 !Overcoordination parameter
3.1353 !Valency/lone pair parameter
0.5000 !Not used
Parameter identifiers
20.0000 !Not used
5.0000 !Molecular energy (not used)
0.0000 !Molecular energy (not used)
1.6052 !Valency angle conjugation parameter
2
cov.r; valency; a.m;
Rvdw;
Evdw;
gammaEEM; cov.r2; #el
alfa; gammavdW; valency; Eunder; n.u.;
chiEEM;
etaEEM; n.u.
cov r3;Elp;
Heat inc.;13BO1; 13BO2; 13BO3;
n.u.;
n.u.
ov/un; val1;
n.u.;
val3; vval4; n.u.;
n.u.;
n.u.
C
1.3826
4.0000 12.0000
2.0195
0.0763
0.8712
1.2360
4.0000
10.6359
1.9232
4.0000 40.5154
0.0000
5.7254
6.9235
0.0000
Atom
1.1663
0.0000 200.0498
6.1551 28.6991 12.1086
0.0000
0.0000
-14.1953
3.5288
0.0000
6.2998
2.9663
0.0000
0.0000
0.0000
H
0.6510
1.0000
1.0080
1.7693
0.0244
0.7625 -0.1000
1.0000
10.0482
5.2587
1.0000
0.0000
0.0000
3.8196
9.8832
1.0000
-0.1000
0.0000 65.0500
3.7647
2.7644
1.0000
0.0000
0.0000
-13.3669
3.6915
0.0000
6.2998
2.8793
0.0000
0.0000
0.0000
3
Edis1;
Edis2;
Edis3;
pbe1;
pbo5;
13corr; pbo6;
kov
pbe2;
pbo3;
pbo4;
n.u.;
pbo1;
pbo2;
ovcorr;
n.u.
1 1 152.0140 104.0507 72.1693
0.2447 -0.7132
1.0000 23.5135
0.3545
0.1152 -0.2069
9.2317
1.0000 -0.1042
5.9159
1.0000
Bond
1 2 174.2967
0.0000
0.0000 -0.5193
0.0000
1.0000
6.0000
0.4401
18.9231
1.0000
0.0000
1.0000 -0.0099
8.2733
0.0000
2 2 177.8312
0.0000
0.0000 -0.3029
0.0000
1.0000
6.0000
0.6891
10.6518
1.0000
0.0000
1.0000 -0.0191
5.4288
0.0000
1
Evdw;
Rvdw;
alfa;
cov.r;
cov.r2;
cov.r3
Off-diagonal
1 2
0.0404
1.8583 10.3804
1.0376 -1.0000 -1.0000
3
Thetao; ka;
kb;
pconj;
pv2;
kpenal; pv3
1 1 1 70.2140 14.0458
2.0508
0.0000
0.0000 35.9933
1.0400
1 1 2 71.6289 18.4967
8.4619
0.0000
0.0000
0.0000
1.0400
Angle
2 1 2 72.7374 18.0638
2.9517
0.0000
0.2000
0.0000
1.0400
4
V1;
V2;
V3;
V2(BO); vconj;
n.u;
n.u.
1 1 1 1
0.0000 28.8256
0.1796 -4.6957 -1.3130
0.0000
0.0000
1 1 1 2
0.0000 32.8083
0.4536 -4.6087 -1.7236
0.0000
0.0000
Torsion
2 1 1 2
0.0000 36.7455
0.3087 -4.7435 -0.7311
0.0000
0.0000
0 1 2 0
0.0000 00.0000
0.1000 -4.7435
0.0000
0.0000
0.0000
1
Rhb;
Dehb;
vhb1;
vhb2
H-bond
1 2 1
2.0347
0.0000
4.9076
4.2357
}
}
General parameters. The first section of the force field contains the general
parameters, which affect all interactions regardless of atom type. The first line of this
16