Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F81967896
input_commands.html
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Sun, Sep 8, 21:21
Size
82 KB
Mime Type
text/html
Expires
Tue, Sep 10, 21:21 (2 d)
Engine
blob
Format
Raw Data
Handle
20414246
Attached To
rLAMMPS lammps
input_commands.html
View Options
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META
NAME=
"Generator"
CONTENT=
"Cosmo Create 1.0.3"
>
</HEAD>
<BODY>
<H2>
LAMMPS Input Commands
</H2>
<P>
<A
HREF=
"README.html"
>
Return
</A>
to top-level of LAMMPS documentation.
</P>
<P>
This page contains a complete list of valid LAMMPS commands which are
read-in from an input script. It will be easiest to understand if you
read it while looking at sample input scripts in the examples
directory.
</P>
<P>
The script of input commands is read by LAMMPS, one line at a time.
Each command causes LAMMPS to take some action. Usually it simply
causes some internal variable(s) to be set. Or it may cause a data file
to be read in or a simulation to be run. Note that most commands have
default settings, which means you only need use a particular command if
you do not want the default setting.
</P>
<P>
Each LAMMPS input script contains exactly one
"
read data
"
(or
"
read restart
"
) command which defines the problem to be
simulated. All other commands can be split into three categories: (a)
commands that (if used) must appear before the
"
read data
"
command because they define settings needed to correctly read-in the
problem and allocate memory for it, (b) commands that must appear after
the
"
read data
"
command because they act on the specified
problem, and (c) commands that can appear either before or after the
"
read data
"
command. Commands in category (c) are used before
the
"
read data
"
command if a default setting needs to be
changed before the problem description is read-in. They can be used
after the
"
read data
"
command if the user wishes to change a
setting before the next
"
run
"
or
"
minimize
"
command
is used. Other than these restrictions, commands can generally appear
in any order in the input script, although some commands require others
to have been previously specified.
</P>
<P>
Each LAMMPS input script also contains one or more
"
run
"
or
"
minimize
"
commands. These trigger an actual dynamics or
minimization computation to be done. Following a run, new commands from
categories (b) and (c) can be used to change various settings, and
additional
"
run
"
commands can then be used to continue the
previous simulation. LAMMPS continues to read successive lines from the
input script until the end-of-file is reached, which causes LAMMPS to
terminate.
</P>
<P>
This page gives examples of each command, some of which can be
specified in multiple styles. Typically the commands take one or more
parameters. The keyword for each command should begin in the leftmost
column and all characters in the command and its parameters should be
in lower-case (except the word NULL or characters in filenames).
Parameters can be separated by arbitrary numbers of spaces and/or tabs,
so long as the command fits on one line. The remainder of the line
after the last parameter is ignored.
</P>
<P>
The next section outlines the structure of a LAMMPS input script. The
final section gives a detailed description of the commands in
alphabetic order, each with its associated parameters and default
settings.
</P>
<UL>
<LI>
<A
HREF=
"#_cch3_951156975"
>
Structure of a LAMMPS input script
</A>
<LI>
<A
HREF=
"#_cch3_931277455"
>
Alphabetic Listing of Commands
</A>
</UL>
<HR>
<H3>
<A
NAME=
"_cch3_951156975"
>
Structure of a LAMMPS input script
</A></H3>
<P>
Any line starting with a # is a
<A
HREF=
"#_cch3_931276588"
>
comment
</A>
.
Comments can appear anywhere in the input script.
</P>
<P>
(1)
<A
HREF=
"#_cch3_930960479"
>
Initialization settings
</A>
(must appear
before
"
read data
"
or
"
read restart
"
)
</P>
<P>
(2)
<A
HREF=
"#_cch3_951435622"
>
Optional Settings
</A>
(can appear before
and/or after
"
read data
"
or
"
read restart
"
)
</P>
<P>
(3)
<A
HREF=
"#_cch3_951435906"
>
Read in a Problem
</A>
via a
"
<A
HREF=
"#_cch3_931277059"
>
read data
</A>
"
or
"
<A
HREF=
"#_cch3_931277070"
>
read restart
</A>
"
command
</P>
<P>
(4)
<A
HREF=
"#_cch3_951435622"
>
Optional Settings
</A>
(same as (2))
</P>
<P>
(5) Problem Settings (must appear after
"
read data
"
or
"
read restart
"
)
</P>
<UL>
<LI>
<A
HREF=
"#_cch3_930960510"
>
Velocity Creation
</A>
<LI>
<A
HREF=
"#_cch3_951435663"
>
Force Field Parameters
</A>
<LI>
<A
HREF=
"#_cch3_930960516"
>
Constraints
</A>
<LI>
<A
HREF=
"#_cch3_930960490"
>
Ensemble Control
</A>
<LI>
<A
HREF=
"#_cch3_930960485"
>
Output Control
</A>
<LI>
<A
HREF=
"#_cch3_951435640"
>
Integrator Settings
</A>
<LI>
<A
HREF=
"#_cch3_951435646"
>
Minimizer Settings
</A>
</UL>
<P>
(6)
<A
HREF=
"#_cch3_951435914"
>
Perform a Simulation
</A>
via a
"
<A
HREF=
"#_cch3_931277194"
>
run
</A>
"
or
"
<A
HREF=
"#_cch3_931277212"
>
minimize
</A>
"
command
</P>
<P>
Repeat (4), (5), and (6) as desired ...
</P>
<HR>
<H3>
<A
NAME=
"_cch3_930960479"
>
Initialization Settings
</A></H3>
<P>
(if used, must appear before
"
read data
"
or
"
read
restart
"
command)
</P>
<PRE>
<A
HREF=
"#_cch3_931276596"
>
units
</A>
real
<A
HREF=
"#_cch3_951437269"
>
extra memory
</A>
2.0 1.5 2.0 2.5
<A
HREF=
"#_cch3_931276604"
>
dimension
</A>
3
<A
HREF=
"#_cch3_931276624"
>
processor grid
</A>
10 10 10
<A
HREF=
"#_cch3_999182956"
>
periodicity
</A>
0 0 0
<A
HREF=
"#_cch3_999182965"
>
slab volume
</A>
3.0
<A
HREF=
"#_cch3_931276632"
>
newton flag
</A>
3
<A
HREF=
"#_cch3_931276687"
>
true flag
</A>
0
<A
HREF=
"#_cch3_951437278"
>
maximum cutoff
</A>
10.0
<A
HREF=
"#_cch3_931276900"
>
mixing style
</A>
geometric
<A
HREF=
"#_cch3_951437286"
>
restart version
</A>
5
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_951435622"
>
Optional Settings
</A></H3>
<P>
(if used, can appear before and/or after
"
read data
"
or
"
read restart
"
command)
</P>
<PRE>
<A
HREF=
"#_cch3_931276654"
>
neighbor
</A>
2.0 1 1 10 1
<A
HREF=
"#_cch3_931276833"
>
nonbond style
</A>
none
<A
HREF=
"#_cch3_931276833"
>
nonbond style
</A>
lj/cutoff 10.0 0
<A
HREF=
"#_cch3_931276833"
>
nonbond style
</A>
lj/smooth 8.0 10.0
<A
HREF=
"#_cch3_931276833"
>
nonbond style
</A>
lj/shift 10.0 0
<A
HREF=
"#_cch3_931276833"
>
nonbond style
</A>
soft 2.5
<A
HREF=
"#_cch3_931276833"
>
nonbond style
</A>
class2/cutoff 10.0 0
<A
HREF=
"#_cch3_931276833"
>
nonbond style
</A>
lj/charmm 15.0 15.1
<A
HREF=
"#_cch3_931276910"
>
coulomb style
</A>
none
<A
HREF=
"#_cch3_931276910"
>
coulomb style
</A>
cutoff 10.0
<A
HREF=
"#_cch3_931276910"
>
coulomb style
</A>
smooth 8.0 10.0
<A
HREF=
"#_cch3_931276910"
>
coulomb style
</A>
ewald 10.0 1.0E-4
<A
HREF=
"#_cch3_931276910"
>
coulomb style
</A>
pppm 10.0 1.0E-4
<A
HREF=
"#_cch3_931276910"
>
coulomb style
</A>
charmm/switch 15.0 15.1
<A
HREF=
"#_cch3_931276910"
>
coulomb style
</A>
debye 10.0 0.5
<A
HREF=
"#_cch3_931276958"
>
bond style
</A>
none
<A
HREF=
"#_cch3_931276958"
>
bond style
</A>
harmonic
<A
HREF=
"#_cch3_931276958"
>
bond style
</A>
fene/standard
<A
HREF=
"#_cch3_931276958"
>
bond style
</A>
fene/shift
<A
HREF=
"#_cch3_931276958"
>
bond style
</A>
nonlinear
<A
HREF=
"#_cch3_931276958"
>
bond style
</A>
class2
<A
HREF=
"#_cch3_931277007"
>
angle style
</A>
none
<A
HREF=
"#_cch3_931277007"
>
angle style
</A>
harmonic
<A
HREF=
"#_cch3_931277007"
>
angle style
</A>
class2
<A
HREF=
"#_cch3_931277007"
>
angle style
</A>
charmm
<A
HREF=
"#_cch3_931277007"
>
angle style
</A>
cosine
<A
HREF=
"#_cch3_931277020"
>
dihedral style
</A>
none
<A
HREF=
"#_cch3_931277020"
>
dihedral style
</A>
harmonic
<A
HREF=
"#_cch3_931277020"
>
dihedral style
</A>
mutliharmonic
<A
HREF=
"#_cch3_931277020"
>
dihedral style
</A>
class2
<A
HREF=
"#_cch3_931277020"
>
dihedral style
</A>
charmm
<A
HREF=
"#_cch3_931277042"
>
improper style
</A>
none
<A
HREF=
"#_cch3_931277042"
>
improper style
</A>
harmonic
<A
HREF=
"#_cch3_931277042"
>
improper style
</A>
cvff
<A
HREF=
"#_cch3_931277042"
>
improper style
</A>
class2
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_951435906"
>
Read in a Problem
</A></H3>
<PRE>
<A
HREF=
"#_cch3_931277059"
>
read data
</A>
data.lj
<A
HREF=
"#_cch3_931277070"
>
read restart
</A>
restart.100000
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_930960510"
>
Velocity Creation
</A></H3>
<P>
(if used, must appear after
"
read data
"
or
"
read
restart
"
command)
</P>
<PRE>
<A
HREF=
"#_cch3_931277080"
>
create group
</A>
types 1 3
<A
HREF=
"#_cch3_931277080"
>
create group
</A>
molecules 200 300
<A
HREF=
"#_cch3_931277080"
>
create group
</A>
region 0.0 1.0 0.0 1.0 INF 1.0
<A
HREF=
"#_cch3_931277080"
>
create group
</A>
remainder
<A
HREF=
"#_cch3_931299999"
>
rotation zero
</A>
1
<A
HREF=
"#_cch3_931277097"
>
create temp
</A>
uniform 300.0 12345678
<A
HREF=
"#_cch3_931277097"
>
create temp
</A>
gaussian 300.0 12345678
<A
HREF=
"#_cch3_931277097"
>
create temp
</A>
velocity 0.0 0.0 0.0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_951435663"
>
Force Field Parameters
</A></H3>
<P>
(if used, must appear after
"
read data
"
or
"
read
restart
"
command)
</P>
<PRE>
<A
HREF=
"#_cch3_931276848"
>
nonbond coeff
</A>
1 2 1.0 3.45 10.0 (nonbond style lj/cutoff)
<A
HREF=
"#_cch3_931276848"
>
nonbond coeff
</A>
1 2 1.0 3.45 8.0 10.0 (nonbond style lj/smooth)
<A
HREF=
"#_cch3_931276848"
>
nonbond coeff
</A>
1 2 1.0 3.45 2.0 10.0 (nonbond style lj/shift)
<A
HREF=
"#_cch3_931276848"
>
nonbond coeff
</A>
1 2 1.0 30.0 2.5 (nonbond style soft)
<A
HREF=
"#_cch3_931276848"
>
nonbond coeff
</A>
1 2 1.0 3.45 10.0 (nonbond style class2/cutoff)
<A
HREF=
"#_cch3_931276848"
>
nonbond coeff
</A>
1 2 1.0 3.45 1.0 3.45 (nonbond style lj/charmm)
<A
HREF=
"#_cch3_931276666"
>
special bonds
</A>
amber
<A
HREF=
"#_cch3_931276666"
>
special bonds
</A>
0.0 0.0 0.5
<A
HREF=
"#_cch3_931276941"
>
pppm mesh
</A>
32 32 64
<A
HREF=
"#_cch3_931276947"
>
pppm order
</A>
5
<A
HREF=
"#_cch3_931276953"
>
dielectric
</A>
1.0
<A
HREF=
"#_cch3_931276970"
>
bond coeff
</A>
1 100.0 3.45 (bond style harmonic)
<A
HREF=
"#_cch3_931276970"
>
bond coeff
</A>
1 30.0 1.5 1.0 1.0 (bond style fene/standard )
<A
HREF=
"#_cch3_931276970"
>
bond coeff
</A>
1 30.0 1.5 1.0 1.0 0.2 (bond style fene/shift)
<A
HREF=
"#_cch3_931276970"
>
bond coeff
</A>
1 28.0 0.748308 0.166667 (bond style nonlinear)
<A
HREF=
"#_cch3_999724447"
>
angle coeff
</A>
1 30.0 108.0 (angle style harmonic)
<A
HREF=
"#_cch3_999724447"
>
angle coeff
</A>
1 30.0 108.0 30.0 2.5 (angle style charmm)
<A
HREF=
"#_cch3_999724447"
>
angle coeff
</A>
1 30.0 (angle style cosine)
<A
HREF=
"#_cch3_999724456"
>
dihedral coeff
</A>
1 10.0 1 3 (dihedral style harmonic)
<A
HREF=
"#_cch3_999724456"
>
dihedral coeff
</A>
1 2.0 2.0 2.0 2.0 2.0 (dihedral style multiharmonic)
<A
HREF=
"#_cch3_999724456"
>
dihedral coeff
</A>
1 2.0 5 180.0 0.5 (dihedral style charmm)
<A
HREF=
"#_cch3_999724473"
>
improper coeff
</A>
1 20.0 0.0 (improper style harmonic)
<A
HREF=
"#_cch3_999724473"
>
improper coeff
</A>
1 20.0 10.0 (improper style cvff)
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_930960516"
>
Constraints
</A></H3>
<P>
(if used, must appear after
"
read data
"
or
"
read
restart
"
command)
</P>
<PRE>
<A
HREF=
"#_cch3_931277114"
>
fix style
</A>
none
<A
HREF=
"#_cch3_931277114"
>
fix style
</A>
1 setforce 0.0 NULL 0.0
<A
HREF=
"#_cch3_931277114"
>
fix style
</A>
1 addforce 1.0 0.0 0.0
<A
HREF=
"#_cch3_931277114"
>
fix style
</A>
1 aveforce 1.0 0.0 0.0
<A
HREF=
"#_cch3_931277114"
>
fix style
</A>
1 rescale 300.0 300.0 100 20.0 0.5
<A
HREF=
"#_cch3_931277114"
>
fix style
</A>
1 hoover/drag 50.0 50.0 0.001
<A
HREF=
"#_cch3_931277114"
>
fix style
</A>
1 langevin 50.0 50.0 0.01 12345 1 1 1
<A
HREF=
"#_cch3_931277114"
>
fix style
</A>
1 springforce 10.0 NULL NULL 1.0
<A
HREF=
"#_cch3_931277114"
>
fix style
</A>
1 dragforce 10.0 -5.0 NULL 2.0 1.0
<A
HREF=
"#_cch3_931277114"
>
fix style
</A>
1 shake 3 0.001 100
<A
HREF=
"#_cch3_931277163"
>
assign fix
</A>
1 atom 200
<A
HREF=
"#_cch3_931277163"
>
assign fix
</A>
1 molecule 50
<A
HREF=
"#_cch3_931277163"
>
assign fix
</A>
1 type 2
<A
HREF=
"#_cch3_931277163"
>
assign fix
</A>
1 region 0.0 1.0 INF INF 0.0 1.0
<A
HREF=
"#_cch3_931277163"
>
assign fix
</A>
1 bondtype 4
<A
HREF=
"#_cch3_931277163"
>
assign fix
</A>
1 angletype 18 10
<A
HREF=
"#_cch3_931277163"
>
assign fix
</A>
1 remainder
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_930960490"
>
Ensemble Control
</A></H3>
<P>
(if used, must appear after
"
read data
"
or
"
read
restart
"
command)
</P>
<PRE>
<A
HREF=
"#_cch3_931276742"
>
temp control
</A>
none
<A
HREF=
"#_cch3_931276742"
>
temp control
</A>
rescale 300.0 300.0 100 20.0 0.5
<A
HREF=
"#_cch3_931276742"
>
temp control
</A>
replace 300.0 300.0 50 12345678
<A
HREF=
"#_cch3_931276742"
>
temp control
</A>
langevin 50.0 50.0 0.01 123456
<A
HREF=
"#_cch3_931276742"
>
temp control
</A>
nose/hoover 300.0 300.0 0.01
<A
HREF=
"#_cch3_931276784"
>
press control
</A>
none
<A
HREF=
"#_cch3_931276784"
>
press control
</A>
nose/hoover xyz 0.0 0.0 0.001
<A
HREF=
"#_cch3_931276784"
>
press control
</A>
nose/hoover xz 0.0 10.0 5.0 5.0 0.0 10.0 0.001
<A
HREF=
"#_cch3_931276784"
>
press control
</A>
nose/hoover yz NULL NULL 5.0 5.0 0.0 10.0 0.001
<A
HREF=
"#_cch3_931276784"
>
press control
</A>
nose/hoover aniso 0.0 0.0 0.0 0.0 1.0 10.0 0.001
<A
HREF=
"#_cch3_931276784"
>
press control
</A>
nose/hoover aniso 0.0 0.0 0.0 0.0 NULL NULL 0.001
<A
HREF=
"#_cch3_999724492"
>
volume control
</A>
none
<A
HREF=
"#_cch3_999724492"
>
volume control
</A>
linear x 0.0 10.0
<A
HREF=
"#_cch3_999724492"
>
volume control
</A>
linear y -1.0 15.0
<A
HREF=
"#_cch3_999724492"
>
volume control
</A>
linear z -10.0 10.0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_930960485"
>
Output Control
</A></H3>
<P>
(if used, must appear after
"
read data
"
or
"
read
restart
"
command)
</P>
<PRE>
<A
HREF=
"#_cch3_931276675"
>
thermo flag
</A>
50
<A
HREF=
"#_cch3_931276681"
>
thermo style
</A>
0
<A
HREF=
"#_cch3_931276696"
>
dump atoms
</A>
100 filename
<A
HREF=
"#_cch3_931276703"
>
dump velocities
</A>
100 filename
<A
HREF=
"#_cch3_931276712"
>
dump forces
</A>
100 filename
<A
HREF=
"#_cch3_931276719"
>
restart
</A>
1000 1 filename
<A
HREF=
"#_cch3_931276719"
>
restart
</A>
1000 2 file1 file2
<A
HREF=
"#_cch3_931276727"
>
diagnostic
</A>
diffusion 100 filename 3 1.0 -1.0 2.5
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_951435640"
>
Integrator Settings
</A></H3>
<P>
(if used, must appear after
"
read data
"
or
"
read
restart
"
command)
</P>
<PRE>
<A
HREF=
"#_cch3_931276638"
>
timestep
</A>
1.0
<A
HREF=
"#_cch3_931276645"
>
respa
</A>
2 2 4
<A
HREF=
"#_cch3_931277185"
>
reset timestep
</A>
0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_951435646"
>
Minimizer Settings
</A></H3>
<P>
(if used, must appear after
"
read data
"
or
"
read
restart
"
command)
</P>
<PRE>
<A
HREF=
"#_cch3_931277200"
>
min style
</A>
hftn
<A
HREF=
"#_cch3_1001972012"
>
min flag
</A>
10
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_951435914"
>
Perform a Simulation
</A></H3>
<PRE>
<A
HREF=
"#_cch3_931277194"
>
run
</A>
10000
<A
HREF=
"#_cch3_931277212"
>
minimize
</A>
0.0001 9999 50000
</PRE>
<HR>
<HR>
<H3>
<A
NAME=
"_cch3_931277455"
>
Alphabetic Listing of Commands:
</A></H3>
<HR>
<H3>
<A
NAME=
"_cch3_999724447"
>
angle coeff
</A></H3>
<UL>
<LI>
1st parameter = angle type #
<LI>
2nd-Nth parameters = coeffs 1 to N-1
</UL>
<PRE>
coeffs: harmonic
(1) K (energy units)
(2) theta (degrees)
class2
currently not enabled for
"
angle coeff
"
command
must be specified in data file (see
"
read data
"
command)
charmm
(1) K (energy units)
(2) theta (degrees)
(3) K_UB (energy/distance^2)
(4) r_UB (distance)
cosine
(1) K (energy units)
define (or override) angle coefficients for an individual angle type
use appropriate number of coeffs for a particular style
see force_fields.html for meaning of coefficients for each style
these coefficients can also be set in data file
by a
"
Angle Coeffs
"
entry, the most recently defined
coefficients are used
cannot use this command before a
"
read data
"
or
"
read restart
"
is performed,
because memory is not yet allocated for the necessary arrays
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931277007"
>
angle style
</A></H3>
<UL>
<LI>
none = compute no angles
<LI>
harmonic = harmonic angles (class 1)
<LI>
class2 = class 2 angles (and associated cross terms)
<LI>
charmm = harmonic + Urey-Bradley
<LI>
cosine = (1 + cos(theta))
</UL>
<PRE>
define style of angle interactions to use for all 3-body terms
must be used before the
"
read data
"
command (if not using the
default) to tell the program how to read the
"
Angle Coeffs
"
entry
in the data file
can be used after the
"
read data
"
command to change the style to none
coefficients for all angle types must be defined in the data (or restart)
file by a
"
Angle Coeffs
"
entry or by
"
angle coeff
"
commands before a run is performed
Default = harmonic
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931277163"
>
assign fix
</A></H3>
<UL>
<LI>
1st parameter = constraint #
<LI>
2nd parameter = style of group of atoms or bondtype
<LI>
3rd-Nth parameters = coeffs 1 to N-2
</UL>
<PRE>
styles:
</PRE>
<UL>
<LI>
atom = single atom
<LI>
molecule = all atoms in a particular molecule
<LI>
type = single atom type
<LI>
region = geometric region of atoms
<LI>
bondtype = bondtype (only for assigning to fix style SHAKE)
<LI>
angletype = angletype (only for assigning to fix style SHAKE)
<LI>
remainder = rest of unconstrained atoms
</UL>
<PRE>
coeffs: atom
(1) global atom #
molecule
(1) molecule #
type
(1) atom type
region
(1) lower x bound of region
(2) upper x bound of region
(3) lower y bound of region
(4) upper y bound of region
(5) lower z bound of region
(6) upper z bound of region
bondtype
(1) bond type
angletype
(1) angle type
(2) bond type used within that angle
remainder
no other parameters required
assign a group of atoms or a bond type to a particular constraint
use appropriate number of coeffs for a particular style
the constraint itself must first be defined by a
"
fix style
"
command
multiple groups of atoms or bond types can be assigned to the same constraint
the bondtype option can only be assigned to a "fix style" of "shake",
multiple bondtypes can be SHAKEn, so long as the size of clusters of
atoms does not exceed the limit described in the "fix style" command
the angletype option can only be assigned to a "fix style" of "shake",
only a single angletype can be SHAKEn, it is designed to be used
in conjunction with "fix style bondtype" to make clusters of size 3
entirely rigid (e.g. water)
the angletype option enables an additional check when SHAKE constraints
are computed: if a cluster is of size 3 and both bonds in
the cluster are of a bondtype specified by the 2nd paramter of
angletype, then the cluster is SHAKEn with an additional angle
constraint that makes it rigid, using the equilibrium angle appropriate
to the specified angletype
IMPORTANT NOTE: the angletype option has one additional affect, namely
that no angle forces for any angle of type angletype are computed
(since it is assumed those angles will be frozen by being SHAKEn), thus
it will likely cause unintended behavior if the bonds in some atom pairs
within angles of type angletype do not have the appropriate bondtype,
since they will not be SHAKEn but neither will the angle force by computed
for style region, a coeff of INF means + or - infinity (all the way
to the boundary)
an atom can be assigned to multiple constraints, the contraints will be
applied in the reverse order they are assigned to that atom
(e.g. each timestep, the last fix assigned to an atom will be applied
to it first, then the next-to-last applied second, etc)
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276970"
>
bond coeff
</A></H3>
<UL>
<LI>
1st parameter = bond type #
<LI>
2nd-Nth parameters = coeffs 1 to N-1
</UL>
<PRE>
coeffs: harmonic
(1) K (energy units)
(2) r0 (distance units)
fene/standard
(1) k for FENE portion (energy/distance^2 units)
(2) r0 for FENE portion (distance units)
(3) epsilon for LJ portion (energy units)
(4) sigma for LJ portion (distance units)
fene/shift
(1) k for FENE (energy/distance^2 units)
(2) r0 for FENE after shift is performed (distance units)
(3) epsilon for LJ (energy units)
(4) sigma for LJ after shift is performed (distance units)
(5) delta shift distance (distance units)
nonlinear
(1) epsilon (energy units)
(2) r0 (distance units)
(3) lamda (distance units)
class2
currently not enabled for
"
bond coeff
"
command
must be specified in data file (see
"
read data
"
command)
define (or override) bond coefficients for an individual bond type
use appropriate number of coeffs for a particular style
see force_fields.html for meaning of coefficients for each style
these coefficients can also be set in data file
by a
"
Bond Coeffs
"
entry, the most recently defined
coefficients are used
cannot use this command before a
"
read data
"
or
"
read restart
"
is performed,
because memory is not yet allocated for the necessary arrays
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276958"
>
bond style
</A></H3>
<UL>
<LI>
none = compute no bonds
<LI>
harmonic = harmonic springs
<LI>
fene/standard = attractive logarithmic term, repulsive LJ
<LI>
fene/shift = same as fene/standard with shift of bond distance
<LI>
nonlinear = non-linear finite-extension spring (van Swol)
<LI>
class2 = class 2 bonds
</UL>
<PRE>
define style of bond interactions to use between all bonded atoms
must be used before the
"
read data
"
command (if not using the
default) to tell the program how to read the
"
Bond
Coeffs
"
entry in the data file (if one exists)
can be used after the
"
read data
"
command to change the style,
in this case
"
bond coeff
"
commands must also be used to set new
coefficients for each bond type (unless the new style is
"
none
"
)
coefficients for all bond types must be defined in the data (or restart)
file by a
"
Bond Coeffs
"
entry or by
"
bond coeff
"
commands before a run is performed
Default = harmonic
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276588"
>
comments
</A></H3>
<PRE>
blank lines are ignored
lines starting with a # are echoed into the log file
for commands, everything on a line after the last parameter is ignored
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276910"
>
coulomb style
</A></H3>
<UL>
<LI>
1st parameter = style of pairwise Coulomb interactions
<LI>
2nd-Nth parameters = coeffs 1 to N-1
</UL>
<PRE>
styles:
</PRE>
<UL>
<LI>
none = no Coulomb interactions are computed
<LI>
cutoff = use a simple cutoff
<LI>
smooth = use a switch region that goes smoothly to zero
<LI>
ewald = use Ewald summations for long-range effects
<LI>
pppm = use particle-mesh Ewald for long-range effects
<LI>
charmm/switch = use the charmm switch to go smoothly to zero
<LI>
debye = add a Debye/Huckel screening exponential
</UL>
<PRE>
coeffs: none
no other parameters required
cutoff
(1) cutoff distance (distance units)
smooth
(1) inner cutoff (distance units)
(2) outer cutoff (distance units)
ewald
(1) cutoff distance for near-field portion (distance units)
(2) accuracy criterion
pppm
(1) cutoff distance for near-field portion (distance units)
(2) accuracy criterion
charmm/switch
(1) inner cutoff (distance units)
(2) outer cutoff (distance units)
debye
(1) cutoff distance (distance units)
(2) kappa (inverse distance units)
use appropriate number of coeffs for a particular style
normally this command should be used before
"
read data
"
or
"
read restart
"
(if simulating a charged system) to tell LAMMPS how big a force cutoff
is being used, the
"
maximum cutoff
"
command can also serve this
purpose
restart files do not store
"
coulomb style
"
choice or cutoff, so
this should be specified in the input script when running from a restart
file
this command can also be used after
"
read data
"
or
"
read restart
"
to
change the style of Coulomb interactions or the cutoff
if simulated system has no charges, should set
"
coulomb style none
"
to
prevent LAMMPS from doing useless nonbond work, LAMMPS will set
this for you and issue a warning
cutoff distance can be smaller or larger than simulation box dimensions
accuracy criterion means
"
one part in value
"
- e.g. 1.0E-4
Ewald and PPPM accuracy criterion are used in conjunction with cutoff
to partition work between short-range and long-range routines
accuracy criterion effectively determines how many k-space vectors are used
to approximate the energy and forces
for PPPM, accuracy criterion determines mesh spacing (see
"
particle mesh
"
command)
3-d periodic boundary conditions are normally used in conjunction with
Ewald and PPPM, see
"
slab volume
"
command for 2-d Ewald/PPPM
cannot use any Coulomb styles other than none with nonbond style = lj/shift or
nonbond style = soft
Coulomb style = smooth should be used with nonbond style = lj/smooth,
and both should use same inner and outer cutoffs
nonbond style = lj/charmm should be used with coulomb style = charmm/switch
for smooth and charmm/switch styles, outer cutoff must be
>
inner cutoff
for smooth and charmm/switch styles, atom pairs less than the inner cutoff
distance use usual Coulomb, pairs between inner and outer are smoothed,
and the potential goes to 0.0 at the outer cutoff
for smooth style, force is continuously differentiable everywhere
for debye style, an exp(-kappa*r) screening is added to the Coulombic
interaction
Default = cutoff 10.0 for real units
cutoff 2.5 for lj units
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931277080"
>
create group
</A></H3>
<UL>
<LI>
1st parameter = style of group of atoms
<LI>
2nd-Nth parameters = coeffs 1 to N-1
</UL>
<PRE>
styles:
</PRE>
<UL>
<LI>
types = range of atom types
<LI>
molecules = range of molecule IDs
<LI>
region = geometric region of atoms
<LI>
remainder = rest of uninitialized atoms
</UL>
<PRE>
coeffs: types
(1) lowest atom type
(2) highest atom type
molecules
(1) lowest molecule ID
(2) highest molecule ID
region
(1) lower x bound of region
(2) upper x bound of region
(3) lower y bound of region
(4) upper y bound of region
(5) lower z bound of region
(6) upper z bound of region
remainder
no other parameters required
used with
"
create temp
"
commmand to initialize velocities of atoms
by default, the
"
create temp
"
command initializes the velocities of all atoms,
this command limits the initialization to a group of atoms
this command is only in force for the next
"
create temp
"
command, any
subsequent
"
create temp
"
command is applied to all atoms (unless the
"
create group
"
command is used again)
for style types, only atoms with a type such that lo-type
<
= type
<
= hi-type
will be initialized by "create temp"
for style types, lo-type can equal hi-type if just want to specify one type
for style molecules, only atoms belonging to molecules with an ID # such
that lo-ID
<
= type
<
= hi-ID will be initialized by "create temp"
for style molecules, lo-ID can equal hi-ID if just want to specify one molecule
for style region, only atoms within the specified spatial region
will be initialized by "create temp"
for style region, a coeff of INF means + or - infinity (all the way
to the boundary)
for style remainder, only previously uninitialized atoms
will be initialized by "create temp"
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931277097"
>
create temp
</A></H3>
<UL>
<LI>
1st parameter = style of temperature creation
<LI>
2nd-Nth parameters = coeffs 1 to N-1
</UL>
<PRE>
styles:
</PRE>
<UL>
<LI>
uniform = uniform distribution of velocities
<LI>
gaussian = gaussian distribution of velocities
<LI>
velocity = assign specific initial velocity to each atom
</UL>
<PRE>
coeffs: uniform
(1) target T (temperature units)
(2) random # seed (0
< seed
<=
8
digits
)
gaussian
(
1
)
target
T
(
temperature
units
)
(
2
)
random
#
seed
(
0
<
seed
<=
8
digits
)
velocity
(
1
)
x
velocity
component
(
velocity
units
)
(
2
)
y
velocity
component
(
velocity
units
)
(
3
)
z
velocity
component
(
velocity
units
)
initialize
velocities
of
atoms
to
a
specified
temperature
use
appropriate
number
of
coeffs
for
a
particular
style
cannot
be
done
before
a
data
or
restart
file
is
read
by
default
,
velocities
are
created
for
all
atoms
-
this
can
be
overridden
by
first
using
a
"
create
group
"
command
for
uniform
and
Gaussian
styles
velocities
are
created
in
processor-independent
fashion
-
is
slower
but
gives
the
same
initial
state
independent
of
#
of
processors
for
uniform
and
Gaussian
styles
the
momentum
of
the
initialized
atoms
is
also
zeroed
,
but
only
if
all
atoms
are
being
initialized
for
uniform
and
Gaussian
styles
,
RN
are
generated
with
Park
/
Miller
RNG
for
velocity
style
in
2-d
simulations
,
still
specify
z
velocity
component
,
even
though
it
is
ignored
</
PRE
>
<HR>
<H3>
<A
NAME=
"_cch3_931276727"
>
diagnostic
</A></H3>
<UL>
<LI>
1st parameter = nametag of a user routine added to diagnostic.f file
<LI>
2nd parameter = call this user routine every this # of timesteps
<LI>
3rd parameter = file name for this routine's diagnostic output
<LI>
4th parameter = # of remaining parameters (0 to 5)
<LI>
5th-9th parameters = optional parameters to pass to user routine
</UL>
<PRE>
call a user-defined diagnostic routine every this many timesteps
this command can be used multiple times to call different routines
at different frequencies, that use different parameters, and that
send output to different files
value of 0 for 2nd parameter means never call this particular routine
this command causes any previous file associated with this user routine
to be closed
new filename can exist, will be overwritten
if the file name specified is
"
none
"
, then no file is opened
each routine that is added to diagnostic.f and enabled with a
"
diagnostic
"
command will be called at the beginning and end of
each
"
run
"
and every so many timesteps during the run
see *** comments in diagnostic.f for changes that must be made in
that file to enable user diagnostics, LAMMPS must then be re-compiled
and re-linked
see the diagnostic.f file for further information on how to create
routines that operate on internal LAMMPS data, do their own file output,
perform different operations (e.g. setup and clean-up) depending
on when they are called, etc
the optional 5th-9th parameters are stored as internal LAMMPS variables
which can be accessed by the diagnostic routine
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276953"
>
dielectric
</A></H3>
<PRE>
set dielectric constant to this value
Default = 1.0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_999724456"
>
dihedral coeff
</A></H3>
<UL>
<LI>
1st parameter = dihedral type #
<LI>
2nd-Nth parameters = coeffs 1 to N-1
</UL>
<PRE>
coeffs: harmonic
(1) K (energy units)
(2) d (+1 or -1)
(3) n (1,2,3,4,6)
multiharmonic
(1) A_1 (energy units)
(2) A_2 (energy units)
(3) A_3 (energy units)
(4) A_4 (energy units)
(5) A_5 (energy units)
class2
currently not enabled for
"
dihedral coeff
"
command
must be specified in data file (see
"
read data
"
command)
charmm
(1) K (energy units)
(2) n (1,2,3,4,6)
(3) d (0 or 180 degrees) (converted to radians within LAMMPS)
(4) weighting factor to turn on/off 1-4 neighbor nonbond interactions
define (or override) dihedral coefficients for an individual dihedral type
use appropriate number of coeffs for a particular style
see force_fields.html for meaning of coefficients for each style
these coefficients can also be set in data file
by a
"
Dihedral Coeffs
"
entry, the most recently defined
coefficients are used
cannot use this command before a
"
read data
"
or
"
read restart
"
is performed,
because memory is not yet allocated for the necessary arrays
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931277020"
>
dihedral style
</A></H3>
<UL>
<LI>
none = compute no dihedrals
<LI>
harmonic = simple harmonic dihedrals (class 1)
<LI>
multiharmonic = multiple simple harmonic dihedrals (class 1)
<LI>
class2 = class 2 dihedrals (and associated cross terms)
<LI>
charmm= simple harmonic dihedrals + charmm 1-4 interactions
</UL>
<PRE>
define style of dihedral interactions to use for all 4-body terms
must be used before the
"
read data
"
command (if not using the
default) to tell the program how to read the
"
Dihedral
Coeffs
"
entry in the data file
can be used after the
"
read data
"
command to change the style to none
coefficients for all dihedral types must be defined in the data (or restart)
file by a
"
Dihedral Coeffs
"
entry or by
"
dihedral coeff
"
commands before a run is performed
Default = harmonic
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276604"
>
dimension
</A></H3>
<UL>
<LI>
specify 3 for 3-d or 2 for 2-d run
</UL>
<PRE>
for a 2-d run, assumes all z-coords are set to 0.0 in
"
read data
"
or
"
read restart
"
files and program creates no z velocities
this command sets the processor grid to default values for 2-d or 3-d
so must be used before
"
processor grid
"
command
must be set before data or restart file is read
Default = 3
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276696"
>
dump atoms
</A></H3>
<UL>
<LI>
1st parameter = # of timesteps
<LI>
2nd parameter = file name
</UL>
<PRE>
dump all atom positions to a file every this many timesteps
(every this many iteration when the minimizer is invoked)
when rRESPA is enabled, this is steps of outermost loop (longest timesteps)
positions are also dumped at the start and end of each run
when dumped during minimization, all dumps will have the same timestamp
since the timestep does not change during minimization
value of 0 means never dump
any previous file is closed
new filename can exist, will be overwritten
atom positions in dump file are in
"
box
"
units (0.0 to 1.0) in each dimension
IMPORTANT NOTE: due to the way periodic boundary conditions are enforced
(only when neighbor lists are rebuilt), atom coords appearing in the dump
file can be slightly outside the specified box
Default = 0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276712"
>
dump forces
</A></H3>
<UL>
<LI>
1st parameter = # of timesteps
<LI>
2nd parameter = file name
</UL>
<PRE>
dump all atom forces to a file every this many timesteps
(every this many iteration when the minimizer is invoked)
when rRESPA is enabled, this is steps of outermost loop (longest timesteps)
forces are also dumped at the start and end of each run
when dumped during minimization, all dumps will have the same timestamp
since the timestep does not change during minimization
any previous file is closed
new filename can exist, will be overwritten
value of 0 means never dump
Default = 0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276703"
>
dump velocities
</A></H3>
<UL>
<LI>
1st parameter = # of timesteps
<LI>
2nd parameter = file name
</UL>
<PRE>
dump all atom velocities to a file every this many timesteps
when rRESPA is enabled, this is steps of outermost loop (longest timesteps)
velocities are also dumped at the start and end of every run
any previous file is closed
new filename can exist, will be overwritten
value of 0 means never dump
Default = 0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_951437269"
>
extra memory
</A></H3>
<UL>
<LI>
1st parameter = extra_own = padding factor on allocation of owned atom
arrays
<LI>
2nd parameter = extra_ghost = padding factor on allocation of ghost
atom arrays
<LI>
3rd parameter = extra_neigh = padding factor on allocation of neighbor
lists
<LI>
4th parameter = extra_buf = padding factor on allocation of
communication buffers
</UL>
<PRE>
factors that affect how much extra memory is allocated when a problem is setup
factor of 1.0 means no padding (use exactly what LAMMPS estimates is
needed), factor of 2.0 means 2x longer arrays, etc
typically don't need to change default settings unless LAMMPS tells you
to
"
boost
"
some factor at run-time
final section of log file lists optimal settings for these parameters,
i.e. the job could have been run with those
"
extra memory
"
settings
and would have used minimal memory
must be set before data or restart file is read
Default = 1.5 for all 4 parameters
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931277114"
>
fix style
</A></H3>
<UL>
<LI>
1st parameter = constraint # (except for none)
<LI>
2nd parameter = style of that constraint
<LI>
3rd-Nth parameters = coeffs 1 to N-2
</UL>
<PRE>
styles:
</PRE>
<UL>
<LI>
none = erase all constraints and all atom and bond assignments
<LI>
setforce = set force on each atom in group
<LI>
addforce = add a force to each atom in group
<LI>
aveforce = apply an external force to group of atoms such that every
atom is accelerated the same
<LI>
rescale = thermostat a group of atoms by rescaling their velocities
<LI>
hoover/drag = thermostat a group of atoms by the Hoover method
<LI>
langevin = thermostat a group of atoms by the Langevin method
<LI>
springforce = apply a spring force to each atom in group
<LI>
dragforce = drag each atom in group to a specified position
<LI>
shake = apply bond length constraints to certain bonds, enabling longer
timesteps
</UL>
<PRE>
coeffs: none
no other parameters required (use
"
none
"
as 1st parameter)
setforce
(1) x component of set force (in force units)
(2) y component of set force (in force units)
(3) z component of set force (in force units)
addforce
(1) x component of added force (in force units)
(2) y component of added force (in force units)
(3) z component of added force (in force units)
aveforce
(1) x comp of added average force per atom (in force units)
(2) y comp of added average force per atom (in force units)
(3) z comp of added average force per atom (in force units)
rescale
(1) desired T at beginning of run
(2) desired T at end of run
(3) check for rescaling every this many timesteps
(4) T window outside of which velocities will be rescaled
(5) fractional amount (0.0 to 1.0) of rescaling to perform
hoover/drag
(1) desired T at beginning of run
(2) desired T at end of run
(3) damping constant for drag (roughly inverse time units)
langevin
(1) desired T at beginning of run
(2) desired T at end of run
(3) Langevin damping parameter (inverse time units)
(4) random seed to use for white noise (0
< seed
<=
8
digits
)
(
5
)
0
/
1 =
off/on
x
dimension
(
6
)
0
/
1 =
off/on
y
dimension
(
7
)
0
/
1 =
off/on
z
dimension
springforce
(
1
)
x
position
of
spring
origin
(
2
)
y
position
(
3
)
z
position
(
4
)
force
constant
k
(
so
that
k
*
distance =
force
units
)
dragforce
(
1
)
x
position
to
drag
atom
towards
(
2
)
y
position
(
3
)
z
position
(
4
)
force
magnitude
f
(
in
force
units
)
(
5
)
delta
outside
of
which
to
apply
force
(
in
distance
units
)
shake
(
1
)
max
#
of
SHAKE
iterations
within
each
atom
cluster
(
2
)
SHAKE
tolerance
(
accuracy
of
one
part
in
tolerance
)
(
3
)
print
bond
statistics
every
this
many
steps
(
0 =
never)
define
a
constraint
cannot
skip
a
constraint
number
,
all
must
be
used
before
a
run
is
performed
use
appropriate
number
of
coeffs
for
a
particular
style
which
atoms
or
bonds
the
constraint
will
affect
is
set
by
the
"
assign
fix
"
command
all
of
the
constraints
(
except
for
rescale
)
are
applied
every
timestep
all
specified
temperatures
are
in
temperature
units
for
style
setforce
,
a
coeff
of
NULL
means
do
not
alter
that
force
component
for
style
aveforce
,
average
force
on
the
group
of
fixed
atoms
is
computed
,
then
new
average
force
is
added
in
and
actual
force
on
each
atom
is
set
to
new
total
value
-
>
has effect of applying same force to entire group
of atoms
thermostatting constraints (rescale, hoover/drag, langevin) cannot be used in
conjuction with global
"
temp control
"
, since they conflict and will
cause atom velocities to be reset twice
thermostatting constraints (rescale, hoover/drag, langevin) cannot be used
when performing a minimization
if multiple Langevin constraints are specified the Marsaglia RNG will
only use the last RNG seed specified for initialization
meaning of rescale and Langevin thermostatting coefficients is same as in
"
temp control
"
command
for rescale style, it can be used as a coarse temperature rescaler,
for example
"
rescale 200.0 300.0 100 10.0 1.0
"
will ramp the temperature
up during the simulation, resetting it to the target temperatue as needed
for rescale style, it can be used to create an instantaneous
drag force that slowly rescales the temperature without oscillation,
for example
"
rescale 300.0 300.0 1 0.0 0.0001
"
will force (or keep)
the temperature to be 300.0, the time frame over which this occurs
will become longer as the last parameter is made smaller
for hoover/drag style, the drag force accumulates over time so some
oscillation in temperature can occur, for example
"
rescale 300.0 300.0 1 0.0 0.0001
"
will force (or keep)
the temperature to be 300.0, the time frame over which the oscillations
occur will become longer as the last parameter is made smaller
style springforce is designed to be applied to an entire group of atoms
en masse (e.g. an umbrella force on an entire molecule)
for springforce style, the center of mass r0 of the group of atoms is computed,
then a restoring force = -k*(r-r0)*mass/masstotal is applied to each
atom in the group where mass = mass of the atom and masstotal = mass of
all the atoms in the group - thus
"
k
"
should represent the total
force on the group of atoms (not per atom)
for springforce style, a xyz position of NULL means do not include that
dimension in the distance or force computation
for dragforce style, apply a drag force of magnitude f to each atom in the
group in the direction (r-r0) where r0 = (x,y,z) - do not apply the force if
the atom is within a distance delta of r0
for dragforce style, a xyz position of NULL means do not include that
dimension in the distance or force computation
for shake style, certain bonds in the system are constrained every timestep
to be at their equilibrium length, this is done by applying a SHAKE-like
constraint to the forces on the atoms so that their position at the next
timestep will preserve the atom separations
for shake style, only atoms in small clusters can be constrained -
e.g. water molecules, CH3 groups, but not the C backbone of a
long polymer chain - a cluster is defined as a central atom
connected to others in the cluster by constrained bonds connected
together by constrained bonds - the max size of such a cluster is
4 atoms to enable easier parallelization
for shake style, the max iteration count need not be large (e.g. 3) since
iterations are only done within a cluster and converge quickly
see the
"
minimize
"
command for what constraints are allowed for use
with the minimizer
see the
"
respa
"
command for how constraints are applied when rRESPA
timestepping is enabled
Default = none
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_999724473"
>
improper coeff
</A></H3>
<UL>
<LI>
1st parameter = improper type #
<LI>
2nd-Nth parameters = coeffs 1 to N-1
</UL>
<PRE>
coeffs: harmonic
(1) K (energy units)
(2) chi (degrees)
cvff
(1) K (energy units)
(2) d (+1 or -1)
(3) n (0,1,2,3,4,6)
class2
currently not enabled for
"
improper coeff
"
command
must be specified in data file (see
"
read data
"
command)
define (or override) improper coefficients for an individual improper type
use appropriate number of coeffs for a particular style
see force_fields.html for meaning of coefficients for each style
these coefficients can also be set in data file
by a
"
Improper Coeffs
"
entry, the most recently defined
coefficients are used
cannot use this command before a
"
read data
"
or
"
read restart
"
is performed,
because memory is not yet allocated for the necessary arrays
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931277042"
>
improper style
</A></H3>
<UL>
<LI>
none = compute no impropers
<LI>
harmonic = harmonic impropers
<LI>
cvff = cvff improper (class 1 variant)
<LI>
class2 = class 2 Wilson out-of-plane
</UL>
<PRE>
define style of improper interactions to use for all trigonal centers
in class2 case, dictates that angle-angle terms be included for all
trigonal and tetrahedral centers
angle for harmonic is improper torsion, angle for class2 is Wilson out-of-plane
must be used before the
"
read data
"
command (if not using the
default) to tell the program how to read the
"
Improper
Coeffs
"
entry in the data file
can be used after the
"
read data
"
command to change the style to none
coefficients for all improper types must be defined in the data (or restart)
file by a
"
Improper Coeffs
"
entry or by
"
improper coeff
"
commands before a run is performed
Default = harmonic
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_951437278"
>
maximum cutoff
</A></H3>
<PRE>
specifies the longest force cutoff that will be used in any runs
this value is used by LAMMPS to accurately allocate memory
for neighbor arrays
if the value is inaccurate (e.g. the command is not used), it is not an
error, but LAMMPS may allocate insufficient memory for neighbor lists
this command is not typically needed if the
"
nonbond style
"
and
"
coulomb style
"
commands are used before the
"
read data
"
or
"
read restart
"
command, since
they specify the appropriate cutoffs
an exception to this is if a short cutoff is used initially,
but a longer cutoff will be used for a subsequent run (in the same
input script), in this case the
"
maximum cutoff
"
command should be
used to insure enough memory is allocated for the later run
note that a restart file contains nonbond cutoffs (so it is not necessary
to use a
"
nonbond style
"
command before
"
read restart
"
), but LAMMPS
still needs to know what the maximum cutoff will be before the
restart file is read
must be set before data or restart file is read
Default = cutoffs for nonbond and Coulomb styles
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_1001972012"
>
min flag
</A></H3>
<PRE>
write out minimization info every this many iterations
value of 0 means never write
Default = 1
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931277200"
>
min style
</A></H3>
<UL>
<LI>
hftn = Hessian-free truncated Newton method
</UL>
<PRE>
choose minimization algorithm to use when
"
minimize
"
command is performed
currently, the hftn style is the only option available
Default = hftn
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931277212"
>
minimize
</A></H3>
<UL>
<LI>
1st parameter = stopping tolerance (in force units)
<LI>
2nd parameter = max iterations of minimizer
<LI>
3rd parameter = max number of force or energy evaluations
</UL>
<PRE>
perform an energy minimization of the atomic coordinates of the system
uses algorithm selected with
"
min style
"
command
minimize commands can be interspersed with
"
run
"
commands
to alternate between dynamics and relaxation of the system
minimization stops if any of 3 criteria are met:
(1) largest force component
< stopping
tolerance
(
2
)
#
of
iterations
>
max iterations
(3) # of force and energy evaluations
>
max evaluations
output from the minimizer is specified by the
"
dump atoms
"
,
"
dump forces
"
,
and
"
restart
"
commands
when using constraints with the minimizer, fixes are
applied when atoms move except for the following
fixes associated with temperature control are not allowed
(rescale, hoover/drag, langevin)
the minimizer does not invoke the
"
fix style shake
"
contraints on
bond lengths
the minimizer does not invoke pressure control or volume control settings
for good convergence, should specify use of smooth nonbond force fields
that have continuous second derivatives, e.g. set
"
coulomb style
"
to
"
smooth
"
or
"
pppm
"
, set
"
nonbond style
"
to
"
lj/smooth
"
or
use a long cutoff
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276900"
>
mixing style
</A></H3>
<UL>
<LI>
1st parameter = style of mixing used to generate i-j nonbond
interactions
</UL>
<PRE>
styles:
</PRE>
<UL>
<LI>
geometric = sqrt(i*j) for both epsilon and sigma
<LI>
arithmetic = sqrt(i*j) for epsilon, (i+j)/2 for sigma
<LI>
sixthpower = see force_fields file for details
</UL>
<PRE>
determine the kind of mixing rule that is applied to generate nonbond
coefficients for interactions between type i and type j atoms
mixing rules are applied only when nonbond coeffs are input in a
"
read data
"
file
for nonbond style "soft", only epsilons (prefactor A) are input - they are
always mixed geometrically, regardless of mixing style setting
for nonbond style "lj/charmm", mixing style is always arithmetic,
regardless of mixing style setting
must be set before data file is read
Default = geometric for all nonbond styles except
arithmetic for nonbond style lj/charmm
sixthpower for nonbond style class2/cutoff
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276654"
>
neighbor
</A></H3>
<UL>
<LI>
1st parameter = skin distance in distance units
<LI>
2nd parameter = neighboring style: 0 = N^2, 1 = binning
<LI>
3rd parameter = build neighbor list every this many steps (see next
param)
<LI>
4th parameter = delay building until after this many steps since last
build
<LI>
5th parameter = build criteria: 0 = always build, 1 = only build if
some atom has moved 1/2 or more of the skin thickness
</UL>
<PRE>
factors that affect how and when neighbor lists are constructed
the binning style is almost always faster than the N^2 style
skin must be large enough that all atoms needed for bond
interactions are also acquired by interprocessor communication
last parameter incurs extra checking and communication to test against
skin thickness, but may mean neighbor list is created less often
when rRESPA is run, the 3rd and 4th parameters refer to the
nonbond (short-range) timestepping
normally this command should be used before the data or restart file is read,
since the skin distance is used to estimate memory needed for
neighbor lists
this command can also be used after the
"
read data
"
or
"
read restart
"
command
to change the style of neighbor list construction, but if the
skin distance is changed it can cause LAMMPS to run out of neighbor
list memory, the
"
maximum cutoff
"
command can be used to avoid this
problem
Default = 2.0 1 1 10 1 for real units
0.3 1 1 10 1 for lj units
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276632"
>
newton flag
</A></H3>
<PRE>
turn off or on Newton's 3rd law for bond and non-bond force computation
</PRE>
<UL>
<LI>
value = 0 = no Newton's 3rd law for either
<LI>
value = 1 = Newton's 3rd law only for bonded computations
<LI>
value = 2 = Newton's 3rd law only for non-bonded computations
<LI>
value = 3 = Newton's 3rd law for both bonded and non-bonded
computations
</UL>
<PRE>
no Newton's 3rd law means more force computation and less communication
yes Newton's 3rd law means less force computation and more communication
which choice is faster is problem dependent on N, # of processors,
and cutoff length(s)
expect for round-off errors, setting this flag should not affect answers,
only run time
must be set before data or restart file is read
Default = 3
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276848"
>
nonbond coeff
</A></H3>
<UL>
<LI>
1st parameter = 1st atom type
<LI>
2nd parameter = 2nd atom type
<LI>
3rd-Nth parameters = coeffs 1 to N-2
</UL>
<PRE>
coeffs: lj/cutoff
(1) epsilon (energy units)
(2) sigma (distance units)
(3) cutoff (distance units)
lj/smooth
(1) epsilon (energy units)
(2) sigma (distance units)
(3) inner cutoff (distance units)
(4) outer cutoff (distance units)
lj/shift
(1) epsilon (energy units)
(2) sigma (distance units)
(3) delta shift distance (distance units)
(4) cutoff (distance units)
soft
(1) prefactor A at start of run (energy units)
(2) prefactor A at end of run (energy units)
(3) cutoff (distance units)
class2/cutoff
(1) epsilon (energy units)
(2) sigma (distance units)
(3) cutoff (distance units)
lj/charmm
(1) epsilon (energy units)
(2) sigma (distance units)
(3) epsilon for 1-4 interactions (energy units)
(4) sigma for 1-4 interactions (distance units)
define (or override) nonbond coefficients for an individual atom type pair
use appropriate number of coeffs for a particular style
1st atom type must be
<
= 2nd atom type
all cutoffs are in global units, not local sigma units
(e.g. in reduced units a setting of "lj/cutoff 1.0 1.2 2.5" means a
cutoff of 2.5, not 1.2*2.5)
turn off a particular type pair interaction by setting the
cutoff to 0.0 (both cutoffs to zero for lj/smooth option)
for soft style, prefactor A is ramped from starting value to
ending value during run
these coefficients (except the cutoffs) can also be set in data file
by a "Nonbond Coeffs" entry and associated mixing rules, the cutoffs can
be set (globally) via the "nonbond style" command, the most
recently defined coefficients/cutoffs are used
cannot use this command before a "read data" or "read restart" is performed,
because memory is not yet allocated for the necessary arrays
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276833"
>
nonbond style
</A></H3>
<UL>
<LI>
1st parameter = style of pairwise nonbond interactions (other than
Coulombic)
<LI>
2nd-Nth parameters = coeffs 1 to N-1
</UL>
<PRE>
styles:
</PRE>
<UL>
<LI>
none = no nonbond interactions are computed
<LI>
lj/cutoff = LJ with a cutoff
<LI>
lj/smooth = LJ with a switched region that goes smoothly to zero
<LI>
lj/shift = same as lj/cutoff with shift of interparticle distance
<LI>
soft = cosine potential with time-varying prefactor
<LI>
class2/cutoff
<LI>
lj/charmm = LJ with charmm switched region that goes smoothly to zero
</UL>
<PRE>
coeffs: none
no other parameters required
lj/cutoff
(1) cutoff (distance units)
(2) offset flag (0 or 1)
lj/smooth
(1) inner cutoff (distance units)
(2) outer cutoff (distance units)
lj/shift
(1) cutoff (distance units)
(2) offset flag (0 or 1)
soft
(1) cutoff (distance units)
class2/cutoff
(1) cutoff (distance units)
(2) offset flag (0 or 1)
lj/charmm
(1) inner cutoff (distance units)
(2) outer cutoff (distance units)
define style of pairwise nonbond interactions to use between all atom types
use appropriate number of coeffs for a particular style
this is separate from charge interactions (see
"
coulomb style
"
command)
normally this command should be used before
"
read data
"
to tell LAMMPS how big a force cutoff is being used, the
"
maximum cutoff
"
command can also serve this purpose
when running from a restart file, the restart file contains the nonbond
style and nonbond cutoffs (but not the offset flag), so it is often
not necessary to use a
"
nonbond style
"
command before
"
read restart
"
,
however LAMMPS still needs to know what the maximum cutoff will be
before the restart file is read, see
"
maximum cutoff
"
command
for more details
this command can also be used after
"
read data
"
or
"
read restart
"
to
change the style of nonbond interactions and/or the cutoff
cutoff distance can be smaller or larger than simulation box dimensions
nonbond style determines how many nonbond coefficients the program expects to
find in a
"
Nonbond Coeffs
"
entry in the data file or when using the
"
nonbond coeff
"
command, thus the style must be set (if not using default)
before using the
"
read data
"
command (if the data file contains a
"
Nonbond Coeffs
"
entry) or a
"
nonbond coeff
"
command
coefficients for all atom type pairs must be defined in data (or restart)
file by a
"
Nonbond Coeffs
"
entry or by
"
nonbond coeffs
"
commands before
a run is performed
this command sets the cutoff(s) for all type pair interactions, thus
overriding any previous settings by a
"
nonbond coeff
"
command or
that were read in from a data or restart file
for lj/cutoff, lj/shift, class2/cutoff styles,
offset flag only affects printout of thermodynamic energy
(not forces or dynamics), determines whether offset energy
is added in to LJ potential to make value at cutoff = 0.0,
flag = 0 -
>
do not add in offset energy,
flag = 1 -
>
add in offset energy
for lj/smooth and lj/charmm styles, outer cutoff must be
>
inner cutoff
for lj/smooth and lj/charmm styles, atom pairs less than the inner cutoff
distance use straight LJ, pairs between inner and outer use a smoothed LJ,
and the potential goes to 0.0 at the outer cutoff
for lj/smooth and lj/charmm styles, energy and forces are continuous at inner
cutoff and go smoothly to zero at outer cutoff
for lj/shift and soft styles, must set
"
coulomb style
"
to
"
none
"
for lj/charmm style, must set
"
coulomb style
"
to
"
charmm/switch
"
,
"
pppm
"
,
or
"
ewald
"
for lj/shift style, delta shift distances for each atom pair are set by
"
Nonbond Coeffs
"
entry in data file or by
"
nonbond coeffs
"
command
for soft style, values of the prefactor
"
A
"
, which is ramped from one
value to another during the run, are set by
"
Nonbond Coeffs
"
entry
in data file or by
"
nonbond coeffs
"
command
Default = lj/cutoff 10.0 0 for real units
lj/cutoff 2.5 0 for lj units
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_999182956"
>
periodicity
</A></H3>
<UL>
<LI>
1st parameter = periodic BC in x direction (0) yes, (1) no
<LI>
2nd parameter = periodic BC in y direction (0) yes, (1) no
<LI>
3rd parameter = periodic BC in z direction (0) yes, (1) no
</UL>
<PRE>
turn on/off periodicity in any of three dimensions
used in inter-particle distance computation and when particles move
to map (or not map) them back into periodic box
for a 2-d run (see
"
dimension
"
command), 3rd parameter must be
specified, but doesn't matter if it is 0 or 1
must be set before data or restart file is read
Default = 0 0 0 (periodic in all dimensions)
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276941"
>
pppm mesh
</A></H3>
<UL>
<LI>
1st parameter = # of mesh points in x direction
<LI>
2nd parameter = # of mesh points in y direction
<LI>
3rd parameter = # of mesh points in z direction
</UL>
<PRE>
specify the mesh size used by
"
coulomb style pppm
"
mesh dimensions that are power-of-two are fastest for FFTs, but any sizes
can be used that are supported by native machine libraries
this command is optional - if not used, a default
mesh size will be chosen to satisfy accuracy criterion - if used, the
specifed mesh size will override the default
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276947"
>
pppm order
</A></H3>
<PRE>
specify the order of the interpolation function that is used by
"
coulomb
style pppm
"
to map particle charge to the particle mesh
order is roughly equivalent to how many mesh points a point charge
overlaps onto in each dimension
Default = 5
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276784"
>
press control
</A></H3>
<UL>
<LI>
1st parameter = style of pressure control
<LI>
2nd parameter = pressure coupling
<LI>
3rd-9th parameters = coeffs 1 to 7
</UL>
<PRE>
styles:
</PRE>
<UL>
<LI>
none = no control (constant volume)
<LI>
nose/hoover = Nose-Hoover constant P
</UL>
<PRE>
coupling:
</PRE>
<UL>
<LI>
xyz = couple all 3 dimensions together (isotropic)
<LI>
xy or yz or xz = couple 2 dimensions together, other is independent
<LI>
aniso = all 3 dimensions are independent (anisotropic)
</UL>
<PRE>
coeffs: none
no other parameters required
nose/hoover xyz
(1) desired P at beginning of run
(2) desired P at end of run
nose/hoover xy or yz or xz or aniso
(1) desired Px at beginning of run (or NULL, see below)
(2) desired Px at end of run
(3) desired Py at beginning of run
(4) desired Py at end of run
(5) desired Pz at beginning of run
(6) desired Pz at end of run
(7) frequency constant for volume adjust (inverse time units)
enable constant pressure simulations
all specified pressures are in pressure units
any dimension being varied by pressure control must be periodic
for xyz coupling, all 3 dimensions expand/contract together uniformly
using total scalar pressure as the driving force
for xy/yz/xz coupling, the 2 specified dimensions expand/contract together
uniformly using pressure components averaged over those 2 dimensions
as the driving force, the non-specified dimension will expand/contract
independently using its pressure component as the driving force
for anisotropic, all 3 dimensions expand/contract independently using
individual pressure components as the 3 driving forces
in all cases, the simulation box stays rectilinear (not Parinello-Rahman)
for dimensions coupled together, their specified P values should be the same
a non-coupled dimension (e.g. dimension z for xy option or any dimension
for aniso option) can have 2 NULL values as specified pressures,
which means apply no pressure control in that dimension (constant volume)
target pressure at intermediate points during a run is a ramped value
between the beginning and ending pressure(s)
for nose/hoover style, frequency constant is like an inverse
"
piston
"
mass which determines how rapidly the pressure fluctuates in response to a
restoring force, large frequency -
>
small mass -
>
rapid fluctations
for nose/hoover style, units of frequency/damping constant are
inverse time, so a value of 0.001 means relax in a timespan on the
order of 1000 fmsec (real units) or 1000 tau (LJ units)
IMPORTANT NOTE: the computation of P in LAMMPS does not include
a long-range Van der Waals correction, this introduces a known
error when performing constant P simulations since the correction
factor changes as the box size varies
Default = none
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276624"
>
processor grid
</A></H3>
<UL>
<LI>
1st parameter = # of processors in x dimension
<LI>
2nd parameter = # of processors in y dimension
<LI>
3rd parameter = # of processors in z dimension
</UL>
<PRE>
specify 3-d grid of processors to map to physical simulation domain
for 2-d problem, specify N by M by 1 grid
program will choose these values to best map processor grid to physical
simulation box, only use this command if wish to override program choice
product of 3 parameters must equal total # of processors
must be set before data or restart file is read
Default = none
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931277059"
>
read data
</A></H3>
<PRE>
read the initial atom positions and bond info from the specified file
the format for the data file is specified in the file data_format
if a
"
Velocities
"
entry is not in data file, all atom velocities
are set to 0.0
if a
"
Coeffs
"
entry is in data file, the appropriate
"
style
"
command
command must be used first (unless default setting is used) to tell
LAMMPS how many coefficients to expect
a
"
Nonbond Coeffs
"
entry only contains one set of coefficients for each
atom type, after being read-in mixing rules are applied to
compute the cross-type coefficients, see the
"
mixing style
"
command
and data_format file for more information
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931277070"
>
read restart
</A></H3>
<PRE>
read atom and force-field information from specified file
allows continuation of a previous run
file is binary to enable exact restarts
do not have to restart on same # of processors, but can only do exact
restarts on same # of processors due to roundoff
when restart file is read, warnings are issued if certain parameters
in the restart file do not match current settings (e.g. newton flag,
dimension, periodicity, units) - this usually indicates an error
the restart file stores the
"
nonbond style
"
and many-body styles and
coefficients and cutoffs, so these do not have to be re-specified in the
input script, unless you want to change them
the restart file does not store
"
coulomb style
"
choice or cutoff, so
this should be re-specified in the input script
the restart file stores the constraint assignments for each atom generated
by using the "assign fix" command, it does NOT store the constraint
parameters themselves, so they must be re-specified with "fix style"
commands after the restart file is read - one exception to this is that
SHAKE constraints (bondtype or angletype) are not stored with the
atoms, so they must be re-specified when performing a restart with both
the "fix style" and "assign fix" commands
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931277185"
>
reset timestep
</A></H3>
<PRE>
explicitly reset the timestep to this value
the
"
read data
"
and
"
read restart
"
commands set the timestep to zero
and the file value respectively, so this should be done after those commands
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276645"
>
respa
</A></H3>
<UL>
<LI>
1st parameter = compute bond forces this many times for every one
3/4-body force call
<LI>
2nd parameter = compute 3/4-body forces this many times for every one
nonbond (short-range) force call
<LI>
3rd parameter = compute nonbond (short-range) forces this many times
for every one long-range force call
</UL>
<PRE>
factors that affect sub-cycling of force calculations within rRESPA hierarchy
bonded intramolecular forces are calculated every innermost sub-timestep
bonded 3- and 4-body forces are computed every 1st parameter sub-timesteps
short-range nonbond pairwise forces (LJ, Coulombic) are computed every
(2nd parameter * 1st parameter) sub-timesteps
long-range (Ewald, PPPM) forces are computed every
(3rd parameter * 2nd parameter * 1st parameter) sub-timesteps
the timestepping for all 3 inner loops (bond, 3/4-body, nonbond) is performed
as sub-cycling within the long-range timestepping loop
the fastest (innermost) timestep size is set by the
"
timestep
"
command
when running rRESPA, all input commands that specify numbers of timesteps
(e.g. run, thermo flag, restart, etc) refer to the outermost loop
of long-range timestepping
the only exception to this rule is the
"
neighbor
"
command, where the timestep
parameters refer to the short-range (nonbond) timestepping
when using constraints (via the
"
fix style
"
and
"
assign fix
"
commands)
with rRESPA, the setforce and aveforce constraints are applied at every
level of the hierarchy (whenever forces are computed), the other
constraints are applied only at the short-range (nonbond) level
when using
"
temp control langevin
"
with rRESPA, thermostatting is applied
at the short-range (nonbond) level
rRESPA cannot be used with
"
fix style shake
"
setting all 3 parameters to 1 turns off rRESPA
Default = 1 1 1 (no rRESPA)
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276719"
>
restart
</A></H3>
<UL>
<LI>
1st parameter = # of timesteps
<LI>
2nd parameter = 1 or 2 = naming convention for restart files
<LI>
3rd (and 4th) parameters = file name(s)
</UL>
<PRE>
create a restart file every this many timesteps
value of 0 means never create one
if the style is 1, restart information will be written to files
named filename.timestep and no 4th parameter is needed
if the style is 2, restart information will be written alternately to files
given by the 3rd and 4th parameters, so only 2 restart files ever exist
when the minimizer is invoked this command means create a restart file
at the end of the minimization with the filename filename.timestep.min
a restart file stores atom and force-field information in binary form
allows program to restart from where it left off (see
"
read restart
"
commmand)
Default = 0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_951437286"
>
restart version
</A></H3>
<PRE>
tell LAMMPS that a restart file from an older version of LAMMPS will be read-in
via a
"
read restart
"
command
this command is necessary because older restart files have a different format
valid settings are 2001 (LAMMPS 2001), 2000 (LAMMPS 2000),
6 (LAMMPS 99) or 5 (LAMMPS 5.0)
restart files from earlier versions of LAMMPS are not readable without
some source code modifications
restart files are always written out in the current-version format
regardless of this setting
this must be set before the
"
read restart
"
command is executed
Default = current version of code = 2001
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931299999"
>
rotation zero
</A></H3>
<PRE>
zero out angular momentum when creating velocities for a group of atoms
value of 0 means don't zero out, value of 1 means zero it
Default = 0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931277194"
>
run
</A></H3>
<PRE>
run or continue dynamics for specified # of timesteps
when rRESPA is enabled, this is steps of outermost loop (longest timesteps)
must have performed
"
read data
"
or
"
read restart
"
command first
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_999182965"
>
slab volume
</A></H3>
<PRE>
invoke 2-d slab Ewald/PPPM and set extended slab volume via this ratio
2-d slab Ewald/PPPM can be used for a system that is periodic in x-y,
but not in z
this ratio dampens inter-slab interactions in the z dimension
by providing empty volume between slabs and removing
dipole inter-slab interactions
ratio value is the size of the extended dimension in z divided by
the actual dimension in z
recommended ratio value is 3.0: larger is inefficient, smaller
risks unwanted inter-slab interactions
when 2-d slab Ewald/PPPM is used, z-direction periodicity must be
turned off - e.g. periodicity 0 0 1
when 2-d slab Ewald/PPPM is used, user must prevent particle migration
beyond initial z-bounds, typically by providing walls
2-d slab Ewald/PPPM can only be used only with electrostatically
neutral systems
2-d slab Ewald/PPPM can only be used (for the moment) with constant
volume simulations (no pressure control) - the pressure computation
(printed as thermodynamic data) does not include any slab correction
factor or a volume correction for the extended z direction
must be set before data or restart file is read
Default = none (normal 3-D Ewald/PPPM)
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276666"
>
special bonds
</A></H3>
<UL>
<LI>
charmm (0.0 0.0 0.0)
<LI>
amber (0.0 0.0 0.5/0.8333)
<LI>
1st parameter = nonbond weight applied to 1-2 neighbors
<LI>
2nd parameter = nonbond weight applied to 1-3 neighbors
<LI>
3rd parameter = nonbond weight applied to 1-4 neighbors
</UL>
<PRE>
weighting factors to turn on/off nonbond interactions of atom pairs that
are
"
close
"
in the molecular topology
1-2 neighbors are a pair of atoms connected by a bond
1-3 neighbors are a pair of atoms 2 hops away, etc.
weight values are from 0.0 to 1.0 and are used to multiply the
energy and force interaction (both Coulombic and LJ) between the 2 atoms
weight of 0.0 means no interaction
weight of 1.0 means full interaction
can either specify a single keyword (charmm, amber) or can give
3 numeric values
using the charmm keyword means use the CHARMM force field
settings of 0.0 0.0 0.0, requiring that pair-specific 1-4 interactions
be read in individually (see
"
dihedral style charmm
"
command)
using the amber keyword means use the AMBER force field
settings of 0.0 0.0 N, where N = 0.5 for Van der Waals 1-4 interactions
and 1.0/1.2 for Coulombic 1-4 interactions
Default = CHARMM force field values of 0.0 0.0 0.0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276742"
>
temp control
</A></H3>
<UL>
<LI>
1st parameter = style of temperature control
<LI>
2nd-Nth parameters = coeffs 1 to N-1
</UL>
<PRE>
styles:
</PRE>
<UL>
<LI>
none = no control
<LI>
rescale = instantaneous rescaling
<LI>
replace = Gaussian replacement
<LI>
langevin = Langevin white noise
<LI>
nose/hoover = Nose-Hoover constant T
</UL>
<PRE>
coeffs: none
no other parameters required
rescale
(1) desired T at beginning of run
(2) desired T at end of run
(3) check for rescaling every this many timesteps
(4) T window outside of which velocities will be rescaled
(5) fractional amount (0.0 to 1.0) of rescaling to perform
replace
(1) desired T at beginning of run
(2) desired T at end of run
(3) do Gaussian replacement every this many timesteps
(4) random # seed to use for replacement (0
< seed
<=
8
digits
)
langevin
(
1
)
desired
T
at
beginning
of
run
(
2
)
desired
T
at
end
of
run
(
3
)
Langevin
damping
parameter
(
inverse
time
units
)
(
4
)
random
seed
to
use
for
white
noise
(
0
<
seed
<=
8
digits
)
nose
/
hoover
(
1
)
desired
T
at
beginning
of
run
(
2
)
desired
T
at
end
of
run
(
3
)
frequency
constant
for
friction
force
(
inverse
time
units
)
enable
constant
temperature
simulations
use
appropriate
number
of
coeffs
for
a
particular
style
all
specified
temperatures
are
in
temperature
units
target
temperature
at
intermediate
points
during
run
is
a
ramped
value
between
the
beginning
and
ending
temperatures
for
rescale
style
,
temperature
is
controlled
by
explicitly
rescaling
velocities
towards
the
target
temperature
for
rescale
style
,
rescaling
is
only
done
if
current
temperature
is
beyond
the
target
temperature
plus
or
minus
the
window
value
for
rescale
style
,
the
amount
of
rescaling
is
contfolled
by
the
fractional
amount
(
0
.
0
to
1
.
0
),
e
.
g
.
a
value
of
0
.
5
means
set
the
velocities
to
halfway
between
the
current
and
target
temperature
for
rescale
style
,
it
can
be
used
as
a
coarse
temperature
rescaler
,
for
example
"
rescale
200
.
0
300
.
0
100
10
.
0
1
.
0
"
will
ramp
the
temperature
up
during
the
simulation
,
resetting
it
to
the
target
temperatue
as
needed
for
rescale
style
,
it
can
be
used
to
create
an
instantaneous
drag
force
that
slowly
rescales
the
temperature
without
oscillation
,
for
example
"
rescale
300
.
0
300
.
0
1
0
.
0
0
.
0001
"
will
force
(
or
keep
)
the
temperature
to
be
300
.
0
,
the
time
frame
over
which
this
occurs
will
become
longer
as
the
last
parameter
is
made
smaller
for
replace
style
,
Gaussian
RNs
from
the
Marsaglia
RNG
are
used
for
langevin
style
,
uniform
RNs
from
the
Marsaglia
RNG
are
used
for
replace
and
langevin
styles
,
the
seed
is
used
to
initialize
the
Marsaglia
RNG
,
on
successive
runs
the
RNG
will
just
continue
on
for
replace
and
langevin
styles
,
generated
RNs
depend
on
#
of
processors
so
will
not
get
same
answers
independent
of
#
of
processors
for
replace
and
langevin
styles
,
RNG
states
are
not
saved
in
restart
file
,
so
cannot
do
an
exact
restart
for
langevin
style
,
damping
parameter
means
small
value
-
>
less damping
for nose/hoover style, frequency constant is like an inverse
"
piston
"
mass which determines how rapidly the temperature
fluctuates in response to a restoring force, large frequency -
>
small mass -
>
rapid fluctations
for nose/hoover style, cannot use a end-of-run T of 0.0, must be finite
for langevin and nose/hoover styles, units of frequency/damping constant are
inverse time, so a value of 0.01 means relax in a timespan on the
order of 100 fmsec (real units) or 100 tau (LJ units)
Default = none
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276675"
>
thermo flag
</A></H3>
<PRE>
print thermodynamic info to screen and log file every this many timesteps
value of 0 means never print
Default = 0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276681"
>
thermo style
</A></H3>
<PRE>
determines format of thermodynamic output to screen and log file
</PRE>
<UL>
<LI>
style = 0 -
>
standard output - about 5 lines per entry
<LI>
style = 1 -
>
reduced output - 1 line per entry
<LI>
style = 2 -
>
output with class 2 terms - about 8 lines per entry
</UL>
<PRE>
Default = 0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276638"
>
timestep
</A></H3>
<PRE>
timestep size for MD run (time units)
when rRESPA is enabled, the timestep size is for the innermost (bond) loop
Default = 1.0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276687"
>
true flag
</A></H3>
<PRE>
read atom positions (see
"
read data
"
command) and dump atom positions
(see
"
dump atoms
"
command) in one of 2 formats
</PRE>
<UL>
<LI>
flag = 0 -
>
read/dump only atom positions (remapped to periodic box)
<LI>
flag = 1 -
>
dump atom positions plus integer box counts
<LI>
flag = 2 -
>
read atom positions plus integer box counts
<LI>
flag = 3 -
>
read/dump atom positions plus integer box counts
</UL>
<PRE>
for each dimension, box count of
"
n
"
means add that many box lengths
to get
"
true
"
un-remapped position,
"
n
"
can be positive, negative, or zero
must be set before data or restart file is read
Default = 0
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_931276596"
>
units
</A></H3>
<UL>
<LI>
real or lj
</UL>
<PRE>
set units to one of two options for all subsequent input parameters
option real = conventional units:
</PRE>
<UL>
<LI>
distance = Angstroms
<LI>
time = femtoseconds
<LI>
mass = grams/mole
<LI>
temperature = degrees K
<LI>
pressure = atmospheres
<LI>
energy = Kcal/mole
<LI>
velocity = Angstroms/femtosecond
<LI>
force = grams/mole * Angstroms/femtosecond^2
</UL>
<PRE>
option lj = LJ reduced units:
</PRE>
<UL>
<LI>
distance = sigmas
<LI>
time = reduced LJ tau
<LI>
temperature = reduced LJ temp
<LI>
pressure = reduced LJ pressure
<LI>
energy = epsilons
<LI>
velocity = sigmas/tau
<LI>
force = reduced LJ force (sigmas/tau^2)
</UL>
<PRE>
for LJ units, LAMMPS sets global epsilon,sigma,mass all equal to 1.0
subsequent input numbers in data and command file must be in these units
output numbers to screen and log and dump files will be in these units
this command (if it appears) must be the first command (aside from
comments) in the input script
must be set before data or restart file is read
Default = real
</PRE>
<HR>
<H3>
<A
NAME=
"_cch3_999724492"
>
volume control
</A></H3>
<UL>
<LI>
1st parameter = style of volume control
<LI>
2nd parameter = dimension to control (x,y,z)
<LI>
3rd-4th parameters = lo/hi simulation box boundaries in this dimension
</UL>
<PRE>
styles:
</PRE>
<UL>
<LI>
none = no variation in any dimension (constant volume)
<LI>
linear = uniform expansion or contraction
</UL>
<PRE>
enable volume changes (density changes) during a simulation
specified box boundaries are in distance units
each dimension is controlled separately
dimensions not specified by a
"
volume control
"
command can be left
alone (constant volume or nonperiodic) or controlled by
a
"
press control
"
command
any dimension being varied by volume control must be periodic
the lo/hi values are the desired global simulation box boundaries at
the end of the simulation run
at each timestep, the box is expanded/contracted uniformly from its initial
lo/hi values to the specified ending lo/hi values
initial lo/hi values are specified in the data or restart file or
inherited from the end of the previous run
at each timestep, all atom coordinates are also scaled to the new box
Default = none
</PRE>
</BODY>
</HTML>
Event Timeline
Log In to Comment