Метод: scf/cc-pvtz
Квадратичный fit: E(r) = a*r^2 + b*r + c

a = 0.51672280 Hartree/Å^2
b = -1.58676880 Hartree/Å
c = -77.93642438 Hartree

r_eq (из fit) = 1.53542 Å
E0 (из fit)   = -79.15459940 Hartree

Если использовать запись:
E = 1/2 * k * (r - r0)^2

то:
k = 1.03344560 Hartree/Å^2
k = 648.4969 kcal/(mol·Å^2)

Если использовать запись:
E = k * (r - r0)^2

то:
k = 0.51672280 Hartree/Å^2
k = 324.2485 kcal/(mol·Å^2)