class LennardJones : public feasst::ModelTwoBody

The Lennard-Jones potential is given by

\( U_{LJ} = 4\epsilon [ (\sigma/r)^{2\alpha} - (\sigma/r)^\alpha ] \),

where \(\alpha=6\) is assumed by this class for optimization.

Subclassed by feasst::LennardJonesAlpha

Public Functions

double hard_sphere_threshold() const

Return the threshold for hard sphere interaction.

void serialize(std::ostream &ostr) const

Output a serialized version of the existing model.


  • hard_sphere_threshold: when r < threshold*sigma, return NEAR_INFINITY (default: 0.2).