LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Commands

angle_style command

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:

angle_coeff

Default:

angle_style none