Diatomic bond stretch ($q=r-r_e,\omega =\sqrt{k/\mu }$)

Let $r$ be the instantaneous bond length, $r_e>0$ its equilibrium value, and $q=r-r_e$ the stretching displacement. Assume one fixed electronic quantum state, a smooth potential energy $V(r)$ with a nondegenerate local minimum at $r_e$, and omit molecular rotation initially. The Taylor series has no linear term because $V^{\prime }(r_e)=0$:Keeping only the quadratic term is the harmonic approximation near a potential minimum. The bond force constant $k$ has units $\mathrm{N}{\mathrm{m}}^{-1}$. It measures local stiffness; it is not the bond-dissociation energy, which concerns the depth and distant shape of the entire well.

The two-particle relative coordinate explains the reduced mass. For stretching along a fixed axis in a stationary centre of mass frame, let $u_1,u_2$ be the nuclear displacements from equilibrium. Then $u_2-u_1=q$ and $m_1{u}_1+m_2{u}_2=0$, soHere $T$ is the internal kinetic energy, and $m_1,m_2>0$ are the two nuclear masses. The lighter partner moves farther; the reduced mass is smaller than either individual mass. The restoring force $-dV/dq\simeq -kq$ gives simple harmonic motion:The angular frequency is $\omega$, the ordinary frequency is $\nu$, and ${\tilde{\nu }}_{\mathrm{harm}}$ is the corresponding spectroscopic wavenumber; $c$ is the speed of light. Use $c$ in $\mathrm{cm}{\mathrm{s}}^{-1}$ to obtain ${\tilde{\nu }}_{\mathrm{harm}}$ in $\mathrm{c}{\mathrm{m}}^{-1}$. The symbol $k$ here is a bond force constant, not an angular wavenumber.

For an original dimensionless example, choose $m_1=m_2=2$ and $k=4$ in consistent model units. Then $\mu =1$ and $\omega =2$. Releasing at $q(0)=0.1$ with $\dot{q}(0)=0$ gives $q(t)=0.1\cos (2t)$ and $u_1=-q/2$, $u_2=q/2$. Each atomic nucleus moves only half the stretching amplitude, while the conserved internal mechanical energy above $V_0$ is $\frac{1}{2}k(0.1{)}^2=0.02$. Fixing one partner instead would give $\omega =\sqrt{2}$, demonstrating why a one-body calculation changes the model.

The bond force constant is the curvature of the effective potential energy at the equilibrium bond length. For a diatomic bond stretch, it plays the role of a spring constant. Its definition is local: two wells with the same curvature can have different bond-dissociation energies.

Quantizing the stretching coordinate gives a quantum harmonic oscillator with reduced mass $\mu$:Here $v$ is the vibrational quantum number and $\hslash =h/(2\pi )$ the reduced Planck constant. Consecutive energy eigenvalues have gap $\hslash \omega$; the ground state has zero-point energy $\hslash \omega /2$ above the well minimum. The offset $V_0$ cancels from every transition energy.

A short proof uses creation and annihilation operators: with $\ell =\sqrt{\hslash /(\mu \omega )}$, write $q=\ell (a+a^{†})/\sqrt{2}$. The Quantum Hamiltonian becomes $\hat{H}-V_0=\hslash \omega (a^{†}a+1/2)$. Nonnegativity of $a^{†}a$ gives the lowest level; repeated application of $a^{†}$ raises its energy by $\hslash \omega$ each time. This also gives $⟨q{⟩}_v=0$ andThe ground state therefore has a nonzero spread of bond lengths, even though its mean stretching displacement vanishes. A stationary energy eigenstate has a time-independent position probability density; the classical oscillation curve is not a trajectory of an atomic nucleus in that state.

The parabolic model formally extends $q$ to the whole real line, whereas an actual bond length satisfies $r>0$. It is useful for low levels concentrated near $r_e$, with typical stretching displacements much smaller than $r_e$. It fails near breakup of the chemical bond and when neglected anharmonicity becomes substantial.

Keep the same bond force constant $k$ and approximate the nuclear masses of ${}^1\mathrm{H}$, ${}^2\mathrm{H}$, and ${}^{35}\mathrm{Cl}$ by $1u$, $2u$, and $35u$, respectively. Here $u$ is the unified atomic mass unit, and $\mathrm{D}$ denotes deuterium. These are chosen mass number approximations, not measured atomic masses or natural-abundance averages. With the same chlorine isotope in both molecules,This is the isotope scaling of diatomic vibration prediction. Doubling the light partner's mass nearly doubles the reduced mass, lowering the frequency by about $28.3\%$. The factor differs from $1/\sqrt{2}$ because the chlorine partner also moves.

For the same small initial stretch $A$ and zero velocity, both models start with the same internal mechanical energy $k{A}^2/2$, but the DCl oscillation period is longer by $\sqrt{72/37}\approx 1.394972$. In the centre of mass frame, $u_{\mathrm{H}}=-35q/36$ and $u_{\mathrm{Cl}}=q/36$ for HCl; $u_{\mathrm{D}}=-35q/37$ and $u_{\mathrm{Cl}}=2q/37$ for DCl. These ratios explain which partner moves and are separate from the change in frequency.

The common-$k$ assumption belongs to the leading Born-Oppenheimer approximation: fixed nuclear positions and electric charges determine the electronic potential energy, while nuclear masses govern motion on that curve. Quantitative comparison requires actual atomic masses, departures from that approximation, and separation of anharmonicity and molecular rotation from the measured band. Even with common $k$, the quantum harmonic oscillator predicts a different zero-point energy for each isotope.^[ref]