Syntax:
angle_style style
Examples:
angle_style harmonic angle_style charmm
Description:
Set the formula LAMMPS will use to compute angle interactions between triplets of atoms. The list of atom triplets is specified in the data or restart file and is read in by a read_data or read_restart command. The coefficients for the formula for each angle type can also be specified in those files or via the angle_coeff command. In all the formulas to follow, theta is the angle defined by the triplet of aotms.
A style of none means angle forces are not computed, even if angles are defined.
The charmm style uses the potential
with an additional Urey_Bradley term based on the distance r between the 1st and 3rd atoms in the angle. K, theta0, Kub, and Rub are coefficients defined for each angle type.
The class2 style uses the potential
where Ea is the angle term, Ebb is a bond-bond term, and Eba is a bond-angle term. Theta0 is the equilibrium angle and r1 and r2 are the equilibrium bond lengths. Kn, M, Nn, theta0, r1, r2 are coefficients defined for each angle type.
The cosine style uses the potential
where K is defined for each angle type.
The harmonic style uses the potential
where theta0 is the equilibrium value of the angle, and K is a prefactor. Note that the usual 1/2 factor is included in K. K and theta0 are coefficients defined for each angle type.
Restrictions:
Angle styles can only be set for atom_styles that allow angles to be defined.
Angle styles are part of the "molecular" package. They are only enabled if LAMMPS was built with that package. See the Making LAMMPS section for more info.
Related commands:
Default:
angle_style none