pub fn sn2_model_1d(x: f64) -> f64
1D SN2 model potential: V(x) = -2exp(-2(x+1.5)^2) - 3exp(-2(x-1.5)^2) + 1.5*exp(-8x^2)