"LAMMPS WWW Site"_lws - "LAMMPS Documentation"_ld - "LAMMPS Commands"_lc :c

:link(lws,http://lammps.sandia.gov)
:link(ld,Manual.html)
:link(lc,Commands_all.html)

:line

create_bonds command :h3

[Syntax:]

create_bonds style args ... keyword value ... :pre

style = {many} or {single/bond} or {single/angle} or {single/dihedral} :ule,l
  {many} args = group-ID group2-ID btype rmin rmax
    group-ID = ID of first group
    group2-ID = ID of second group, bonds will be between atoms in the 2 groups
    btype = bond type of created bonds
    rmin = minimum distance between pair of atoms to bond together
    rmax = maximum distance between pair of atoms to bond together
  {single/bond} args = btype batom1 batom2
    btype = bond type of new bond
    batom1,batom2 = atom IDs for two atoms in bond
  {single/angle} args = atype aatom1 aatom2 aatom3
    atype = bond type of new angle
    aatom1,aatom2,aatom3 = atom IDs for three atoms in angle
  {single/dihedral} args = dtype datom1 datom2 datom3 datom4
    dtype = bond type of new dihedral
    datom1,datom2,datom3,datom4 = atom IDs for four atoms in dihedral :pre
zero or more keyword/value pairs may be appended :l
keyword = {special} :l
  {special} value = {yes} or {no} :pre
:ule

[Examples:]

create_bonds many all all 1 1.0 1.2
create_bonds many surf solvent 3 2.0 2.4
create_bonds single/bond 1 1 2
create_bonds single/angle 5 52 98 107 special no :pre

[Description:]

Create bonds between pairs of atoms that meet a specified distance
criteria.  Or create a single bond, angle, or dihedral between 2, 3,
or 4 specified atoms.

The new bond (angle, dihedral) interactions will then be computed
during a simulation by the bond (angle, dihedral) potential defined by
the "bond_style"_bond_style.html, "bond_coeff"_bond_coeff.html,
"angle_style"_angle_style.html, "angle_coeff"_angle_coeff.html,
"dihedral_style"_dihedral_style.html,
"dihedral_coeff"_dihedral_coeff.html commands.

The {many} style is useful for adding bonds to a system, e.g. between
nearest neighbors in a lattice of atoms, without having to enumerate
all the bonds in the data file read by the "read_data"_read_data.html
command.

The {single} styles are useful for adding bonds, angles, dihedrals
to a system incrementally, then continuing a simulation.

Note that this command does not auto-create any angle or dihedral
interactions when a bond is added.  Nor does it auto-create any bonds
when an angle or dihedral is added.  Or auto-create any angles when a
dihedral is added.  Thus the flexibility of this command is limited.
It can be used several times to create different types of bond at
different distances.  But it cannot typically auto-create all the
bonds or angles or dihedral that would normally be defined in a data
file for a complex system of molecules.

NOTE: If the system has no bonds (angles, dihedrals) to begin with, or
if more bonds per atom are being added than currently exist, then you
must insure that the number of bond types and the maximum number of
bonds per atom are set to large enough values.  And similarly for
angles and dihedrals.  Otherwise an error may occur when too many
bonds (angles, dihedrals) are added to an atom.  If the
"read_data"_read_data.html command is used to define the system, these
parameters can be set via the "bond types" and "extra bond per atom"
fields in the header section of the data file.  If the
"create_box"_create_box.html command is used to define the system,
these 2 parameters can be set via its optional "bond/types" and
"extra/bond/per/atom" arguments.  And similarly for angles and
dihedrals.  See the doc pages for these 2 commands for details.

:line

The {many} style will create bonds between pairs of atoms I,J where I
is in one of the two specified groups, and J is in the other.  The two
groups can be the same, e.g. group "all".  The created bonds will be
of bond type {btype}, where {btype} must be a value between 1 and the
number of bond types defined.

For a bond to be created, an I,J pair of atoms must be a distance D
apart such that {rmin} <= D <= {rmax}.

The following settings must have been made in an input script before
this style is used:

special_bonds weight for 1-2 interactions must be 0.0
a "pair_style"_pair_style.html must be defined
no "kspace_style"_kspace_style.html defined
minimum "pair_style"_pair_style.html cutoff + "neighbor"_neighbor.html skin >= {rmax} :ul

These settings are required so that a neighbor list can be created to
search for nearby atoms.  Pairs of atoms that are already bonded
cannot appear in the neighbor list, to avoid creation of duplicate
bonds.  The neighbor list for all atom type pairs must also extend to
a distance that encompasses the {rmax} for new bonds to create.

An additional requirement for this style is that your system must be
ready to perform a simulation.  This means, for example, that all
"pair_style"_pair_style.html coefficients be set via the
"pair_coeff"_pair_coeff.html command.  A "bond_style"_bond_style.html
command and all bond coefficients must also be set, even if no bonds
exist before this command is invoked.  This is because the building of
neighbor list requires initialization and setup of a simulation,
similar to what a "run"_run.html command would require.

Note that you can change any of these settings after this command
executes, e.g. if you wish to use long-range Coulombic interactions
via the "kspace_style"_kspace_style.html command for your subsequent
simulation.

:line

The {single/bond} style creates a single bond of type {btype} between
two atoms with IDs {batom1} and {batom2}.  {Btype} must be a value
between 1 and the number of bond types defined.

The {single/angle} style creates a single angle of type {atype}
between three atoms with IDs {aatom1}, {aatom2}, and {aatom3}.  The
ordering of the atoms is the same as in the {Angles} section of a data
file read by the "read_data"_read_data.html command.  I.e. the 3 atoms are
ordered linearly within the angle; the central atom is {aatom2}.
{Atype} must be a value between 1 and the number of angle types
defined.

The {single/dihedral} style creates a single dihedral of type {btype}
between two atoms with IDs {batom1} and {batom2}.  The ordering of the
atoms is the same as in the {Dihedrals} section of a data file read by
the "read_data"_read_data.html command.  I.e. the 4 atoms are ordered
linearly within the dihedral.  {Dtype} must be a value between 1 and
the number of dihedral types defined.

:line

The keyword {special} controls whether an internal list of special
bonds is created after one or more bonds, or a single angle or
dihedral is added to the system.

The default value is {yes}.  A value of {no} cannot be used
with the {many} style.

This is an expensive operation since the bond topology for the system
must be walked to find all 1-2, 1-3, 1-4 interactions to store in an
internal list, which is used when pairwise interactions are weighted;
see the "special_bonds"_special_bonds.html command for details.

Thus if you are adding a few bonds or a large list of angles all at
the same time, by using this command repeatedly, it is more efficient
to only trigger the internal list to be created once, after the last
bond (or angle, or dihedral) is added:

create_bonds single/bond 5 52 98 special no
create_bonds single/bond  5 73 74 special no
...
create_bonds single/bond 5 17 386 special no
create_bonds single/bond 4 112 183 special yes :pre

Note that you MUST insure the internal list is re-built after the last
bond (angle, dihedral) is added, before performing a simulation.
Otherwise pairwise interactions will not be properly excluded or
weighted.  LAMMPS does NOT check that you have done this correctly.

:line

[Restrictions:]

This command cannot be used with molecular systems defined using
molecule template files via the "molecule"_molecule.html and
"atom_style template"_atom_style.html commands.

[Related commands:]

"create_atoms"_create_atoms.html, "delete_bonds"_delete_bonds.html

[Default:]

The keyword default is special = yes.
