File: pair_reax_c.html

package info (click to toggle)
lammps 0~20120615.gite442279-1
  • links: PTS, VCS
  • area: main
  • in suites: wheezy
  • size: 128,448 kB
  • sloc: cpp: 321,874; fortran: 15,187; ansic: 11,007; python: 7,889; perl: 2,915; sh: 2,088; makefile: 924; f90: 374; objc: 238; lisp: 169; csh: 16; tcl: 6
file content (310 lines) | stat: -rw-r--r-- 13,461 bytes parent folder | download
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
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
<HTML>
<CENTER><A HREF = "http://lammps.sandia.gov">LAMMPS WWW Site</A> - <A HREF = "Manual.html">LAMMPS Documentation</A> - <A HREF = "Section_commands.html#comm">LAMMPS Commands</A> 
</CENTER>






<HR>

<H3>pair_style reax/c command 
</H3>
<P><B>Syntax:</B>
</P>
<PRE>pair_style reax/c cfile keyword value 
</PRE>
<UL><LI>cfile = NULL or name of a control file 

<LI>zero or more keyword/value pairs may be appended 

<PRE>keyword = <I>checkqeq</I> or <I>lgvdw</I>
  <I>checkqeq</I> value = <I>yes</I> or <I>no</I> = whether or not to require qeq/reax fix 
  <I>lgvdw</I> value = <I>yes</I> or <I>no</I> = whether or not to use a low gradient vdW correction 
</PRE>

</UL>
<P><B>Examples:</B>
</P>
<PRE>pair_style reax/c NULL
pair_style reax/c controlfile checkqeq no
pair_style reax/c NULL lgvdw yes
pair_coeff * * ffield.reax 1 2 2 3 
</PRE>
<P><B>Description:</B>
</P>
<P>Style <I>reax/c</I> computes the ReaxFF potential of van Duin, Goddard and
co-workers.  ReaxFF uses distance-dependent bond-order functions to
represent the contributions of chemical bonding to the potential
energy. There is more than one version of ReaxFF.  The version
implemented in LAMMPS uses the functional forms documented in the
supplemental information of the following paper: <A HREF = "#Chenoweth_2008">(Chenoweth et al.,
2008)</A>.  The version integrated into LAMMPS matches
the most up-to-date version of ReaxFF as of summer 2010.  For more
technical details about the pair reax/c implementation of ReaxFF, see
the "(Aktulga)" paper.
</P>
<P>The <I>reax/c</I> style differs from the <A HREF = "pair_reax.html">pair_style reax</A>
command in the lo-level implementation details.  The <I>reax</I> style is a
Fortran library, linked to LAMMPS.  The <I>reax/c</I> style was initially
implemented as stand-alone C code and is now integrated into LAMMPS as
a package.
</P>
<P>LAMMPS provides several different versions of ffield.reax in its
potentials dir, each called potentials/ffield.reax.label.  These are
documented in potentials/README.reax.  The default ffield.reax
contains parameterizations for the following elements: C, H, O, N, S.
</P>
<P>The format of these files is identical to that used originally by van
Duin.  We have tested the accuracy of <I>pair_style reax/c</I> potential
against the original ReaxFF code for the systems mentioned above.  You
can use other ffield files for specific chemical systems that may be
available elsewhere (but note that their accuracy may not have been
tested).
</P>
<P>IMPORTANT NOTE: We do not distribute a wide variety of ReaxFF force
field files with LAMMPS.  Adri van Duin's group at PSU is the central
repository for this kind of data as they are continuously deriving and
updating parameterizations for different classes of materials.  You
can visit their WWW site at
<A HREF = "http://www.engr.psu.edu/adri">http://www.engr.psu.edu/adri</A>, register
as a "new user", and then submit a request to their group describing
material(s) you are interested in modeling with ReaxFF.  They can tell
you what is currently available or what it would take to create a
suitable ReaxFF parameterization.
</P>
<P>The <I>cfile</I> setting can be specified as NULL, in which case default
settings are used.  Or a control file can be specified which contains
cutoff values for the ReaxFF potential in addition to some performance
and output controls.  Each line in the control specifies the value for
a control variable.  The format of the control file is described
below.  
</P>
<P>Two examples using <I>pair_style reax/c</I> are provided in the examples/reax
sub-directory, along with corresponding examples for 
<A HREF = "pair_reax.html">pair_style reax</A>.
</P>
<P>Use of this pair style requires that a charge be defined for every
atom.  See the <A HREF = "atom_style.html">atom_style</A> and
<A HREF = "read_data.html">read_data</A> commands for details on how to specify
charges.
</P>
<P>The ReaxFF parameter files provided were created using a charge
equilibration (QEq) model for handling the electrostatic interactions.
Therefore, by default, LAMMPS requires that the <A HREF = "fix_qeq_reax.html">fix
qeq/reax</A> command be used with <I>pair_style reax/c</I>
when simulating a ReaxFF model, to equilibrate charge each timestep.
Using the keyword <I>checkqeq</I> with the value <I>no</I> 
turns off the check for <I>fix qeq/reax</I>,
allowing a simulation to be run without charge equilibration.
In this case, the static charges you
assign to each atom will be used for computing the electrostatic
interactions in the system.
See the <A HREF = "fix_qeq_reax.html">fix qeq/reax</A> command for details.
</P>
<P>Using the optional keyword <I>lgvdw</I> with the value <I>yes</I> turns on 
the low-gradient correction of the ReaxFF/C for long-range 
London Dispersion, as described in the <A HREF = "#Liu">(Liu)</A> paper. Force field 
file <I>ffield.reax.lg</I> is designed for this correction, and is trained
for several energetic materials (see <A HREF = "#Liu">Liu</A> ). 
When using lg-correction, 
recommended value for parameter <I>thb</I> is 0.01, which can be set in 
the control file.  Note: Force field files 
are different for the original or lg corrected pair styles, using
wrong ffield file generates an error message.
</P>
<P>The thermo variable <I>evdwl</I> stores the sum of all the ReaxFF potential
energy contributions, with the exception of the Coulombic and charge
equilibration contributions which are stored in the thermo variable
<I>ecoul</I>.  The output of these quantities is controlled by the
<A HREF = "thermo.html">thermo</A> command.
</P>
<P>This pair style tallies a breakdown of the total ReaxFF potential
energy into sub-categories, which can be accessed via the <A HREF = "compute_pair.html">compute
pair</A> command as a vector of values of length 14.
The 14 values correspond to the following sub-categories (the variable
names in italics match those used in the original FORTRAN ReaxFF code):
</P>
<OL><LI><I>eb</I> = bond energy
<LI><I>ea</I> = atom energy
<LI><I>elp</I> = lone-pair energy
<LI><I>emol</I> = molecule energy (always 0.0)
<LI><I>ev</I> = valence angle energy
<LI><I>epen</I> = double-bond valence angle penalty
<LI><I>ecoa</I> = valence angle conjugation energy
<LI><I>ehb</I> = hydrogen bond energy
<LI><I>et</I> = torsion energy
<LI><I>eco</I> = conjugation energy 
<LI><I>ew</I> = van der Waals energy
<LI><I>ep</I> = Coulomb energy
<LI><I>efi</I> = electric field energy (always 0.0)
<LI><I>eqeq</I> = charge equilibration energy 
</OL>
<P>To print these quantities to the log file (with descriptive column
headings) the following commands could be included in an input script:
</P>
<PRE>compute reax all pair reax/c
variable eb  	 equal c_reax[1]
variable ea  	 equal c_reax[2] 
...
variable eqeq 	 equal c_reax[14]
thermo_style custom step temp epair v_eb v_ea ... v_eqeq 
</PRE>
<P>Only a single pair_coeff command is used with the <I>reax/c</I> style which
specifies a ReaxFF potential file with parameters for all needed
elements.  These are mapped to LAMMPS atom types by specifying N
additional arguments after the filename in the pair_coeff command,
where N is the number of LAMMPS atom types:
</P>
<UL><LI>filename
<LI>N indices = mapping of ReaxFF elements to atom types 
</UL>
<P>The filename is the ReaxFF potential file.  Unlike for the <I>reax</I>
pair style, any filename can be used.
</P>
<P>In the ReaxFF potential file, near the top, after the general
parameters, is the atomic parameters section that contains element
names, each with a couple dozen numeric parameters.  If there are M
elements specified in the <I>ffield</I> file, think of these as numbered 1
to M. Each of the N indices you specify for the N atom types of LAMMPS
atoms must be an integer from 1 to M.  Atoms with LAMMPS type 1 will
be mapped to whatever element you specify as the first index value,
etc.  If a mapping value is specified as NULL, the mapping is not
performed.  This can be used when the <I>reax/c</I> style is used as part
of the <I>hybrid</I> pair style.  The NULL values are placeholders for atom
types that will be used with other potentials.
</P>
<P>IMPORTANT NOTE: Currently the reax/c pair style cannot be used as part
of the <I>hybrid</I> pair style.  Some additional work still need to be
done to enable this.
</P>
<P>As an example, say your LAMMPS simulation has 4 atom types and the
elements are ordered as C, H, O, N in the <I>ffield</I> file.  If you want
the LAMMPS atom type 1 and 2 to be C, type 3 to be N, and type 4 to be
H, you would use the following pair_coeff command:
</P>
<PRE>pair_coeff * * ffield.reax 1 1 4 2 
</PRE>
<HR>

<P>The format of a line in the control file is as follows:
</P>
<PRE>variable_name value 
</PRE>
<P>and it may be followed by an "!" character and a trailing comment.
</P>
<P>If the value of a control variable is not specified, then default
values are used.  What follows is the list of variables along with a
brief description of their use and default values.
</P>
<P>simulation_name: Output files produced by <I>pair_style reax/c</I> carry
this name + extensions specific to their contents.  Partial energies
are reported with a ".pot" extension, while the trajectory file has
".trj" extension.
</P>
<P>tabulate_long_range: To improve performance, long range interactions
can optionally be tabulated (0 means no tabulation). Value of this
variable denotes the size of the long range interaction table.  The
range from 0 to long range cutoff (defined in the <I>ffield</I> file) is
divided into <I>tabulate_long_range</I> points.  Then at the start of
simulation, we fill in the entries of the long range interaction table
by computing the energies and forces resulting from van der Waals and
Coulomb interactions between every possible atom type pairs present in
the input system.  During the simulation we consult to the long range
interaction table to estimate the energy and forces between a pair of
atoms. Linear interpolation is used for estimation. (default value =
0)
</P>
<P>energy_update_freq: Denotes the frequency (in number of steps) of
writes into the partial energies file. (default value = 0)
</P>
<P>nbrhood_cutoff: Denotes the near neighbors cutoff (in Angstroms)
regarding the bonded interactions. (default value = 5.0)
</P>
<P>hbond_cutoff: Denotes the cutoff distance (in Angstroms) for hydrogen
bond interactions.(default value = 7.5. Value of 0.0 turns off hydrogen 
		bonds)
</P>
<P>bond_graph_cutoff: is the threshold used in determining what is a
physical bond, what is not. Bonds and angles reported in the
trajectory file rely on this cutoff. (default value = 0.3)
</P>
<P>thb_cutoff: cutoff value for the strength of bonds to be considered in
three body interactions. (default value = 0.001)
</P>
<P>thb_cutoff_sq: cutoff value for the strength of bond order products 
to be considered in three body interactions. (default value = 0.00001)
</P>
<P>write_freq: Frequency of writes into the trajectory file. (default
value = 0)
</P>
<P>traj_title: Title of the trajectory - not the name of the trajectory
file.
</P>
<P>atom_info: 1 means print only atomic positions + charge (default = 0)
</P>
<P>atom_forces: 1 adds net forces to atom lines in the trajectory file
(default = 0)
</P>
<P>atom_velocities: 1 adds atomic velocities to atoms line (default = 0)
</P>
<P>bond_info: 1 prints bonds in the trajectory file (default = 0)
</P>
<P>angle_info: 1 prints angles in the trajectory file (default = 0)
</P>
<HR>

<P><B>Mixing, shift, table, tail correction, restart, rRESPA info</B>:
</P>
<P>This pair style does not support the <A HREF = "pair_modify.html">pair_modify</A>
mix, shift, table, and tail options.
</P>
<P>This pair style does not write its information to <A HREF = "restart.html">binary restart
files</A>, since it is stored in potential files.  Thus, you
need to re-specify the pair_style and pair_coeff commands in an input
script that reads a restart file.
</P>
<P>This pair style can only be used via the <I>pair</I> keyword of the
<A HREF = "run_style.html">run_style respa</A> command.  It does not support the
<I>inner</I>, <I>middle</I>, <I>outer</I> keywords.
</P>
<P><B>Restrictions:</B>
</P>
<P>This pair style is part of the USER-REAXC package.  It is only enabled
if LAMMPS was built with that package.  See the <A HREF = "Section_start.html#start_3">Making
LAMMPS</A> section for more info.
</P>
<P>The ReaxFF potential files provided with LAMMPS in the potentials
directory are parameterized for real <A HREF = "units.html">units</A>.  You can use
the ReaxFF potential with any LAMMPS units, but you would need to
create your own potential file with coefficients listed in the
appropriate units if your simulation doesn't use "real" units.
</P>
<P><B>Related commands:</B>
</P>
<P><A HREF = "pair_coeff.html">pair_coeff</A>, <A HREF = "fix_qeq_reax.html">fix_qeq_reax</A>,
<A HREF = "pair_reax.html">pair_style reax</A>
</P>
<P><B>Default:</B>
</P>
<P>The keyword defaults are checkqeq = yes, lgvdw = no.
</P>
<HR>

<A NAME = "Chenoweth_2008"></A>

<P><B>(Chenoweth_2008)</B> Chenoweth, van Duin and Goddard, 
Journal of Physical Chemistry A, 112, 1040-1053 (2008).
</P>
<A NAME = "Aktulga"></A>

<P><B>(Aktulga)</B> Aktulga, Fogarty, Pandit, Grama, Parallel Computing, to
appear (2011).
</P>
<A NAME = "Liu"></A>

<P><B>(Liu)</B> L. Liu, Y. Liu, S. V. Zybin, H. Sun and W. A. Goddard, Journal
of Physical Chemistry A, 115, 11016-11022 (2011).
</P>
</HTML>