class AngleHarmonic : public feasst::BondThreeBody

U(angle) = k_energy_per_radian_sq*(angle - equilibrium_degrees)^2 with parameters given in Angle Properties.

The usual 1/2 factor is not included, but can be incorporated into the k parameter manually by the user input to the particle file.

Public Functions

double energy(const double radians, const Bond &angle) const

Return the energy of interaction given an angle in radians.