V(r) = D_e(1−e^(−a(r−r_e)))² — models diatomic bond with anharmonicity, dissociation limit, and unequal level spacing. Compare quantum Morse levels vs harmonic approximation.