pub fn lennard_jones_potential(r: f64, epsilon: f64, sigma: f64) -> f64
Lennard-Jones pair potential energy between two atoms.
U(r) = 4ε[(σ/r)^12 - (σ/r)^6]