pub fn proton_transfer_evb(s: f64, k: f64, coupling: f64) -> f64
1D Marcus/EVB adiabatic potential for proton transfer s in [-1, 1], symmetric TS exactly at s=0