Reduced Mass

Definition

Two bodies of mass m₁ and m₂ interact through a force that depends only on their separation — gravity between a star and a planet, or the Coulomb attraction between a proton and an electron. Naively that is a coupled six-coordinate problem. The reduced mass is the trick that shrinks it to a one-body problem.

μ = m₁·m₂ / (m₁ + m₂)     equivalently     1/μ = 1/m₁ + 1/m₂

Here:

• μ — reduced mass, in kilograms (kg) or, for atoms/molecules, in unified atomic mass units (u) or electron masses.
• m₁, m₂ — the two body masses (kg).

The reciprocal form 1/μ = 1/m₁ + 1/m₂ is worth memorizing: reduced masses add like resistors in parallel or capacitors in series. Two reciprocals summed give a reciprocal larger than either, so μ is always smaller than the lighter of the two masses.

Why the two-body problem collapses to one body

Write Newton's second law for each body under a mutual force F(r) directed along the separation vector r = r₁ − r₂:

m₁ r̈₁ = +F(r) r̂
m₂ r̈₂ = −F(r) r̂

Now change variables to the center-of-mass coordinate R = (m₁r₁ + m₂r₂)/(m₁+m₂) and the relative coordinate r = r₁ − r₂. Adding the two equations shows the center of mass feels no net internal force:

(m₁ + m₂) R̈ = 0     →     R moves at constant velocity (or stays fixed)

Subtracting them (after dividing each by its mass) gives the payoff:

r̈ = r̈₁ − r̈₂ = F(r)/m₁ + F(r)/m₂ = F(r)·(1/m₁ + 1/m₂)
⇒   μ r̈ = F(r)     with     1/μ = 1/m₁ + 1/m₂

That last line is Newton's second law for a single particle of mass μ at position r, moving under exactly the same force law. The messy coupled system splits cleanly into (1) a boring free particle at the center of mass and (2) one fictitious particle of mass μ carrying all the interesting dynamics. Every orbit, every scattering angle, every bound-state energy of the real two-body system is read straight off this equivalent one-body problem.

Step by step: setting up the equivalent one-body problem

1. Identify the two masses m₁ and m₂ and the interaction that depends only on |r₁ − r₂|.
2. Separate coordinates into R (center of mass) and r (relative). The kinetic energy separates cleanly: T = ½(m₁+m₂)Ṙ² + ½μṙ².
3. Discard the center-of-mass motion. It is force-free — a straight line at constant velocity — so it carries no dynamics of interest. Work in the CM frame where R is fixed.
4. Solve the one-body problem μr̈ = F(r): a single particle of mass μ in a central potential V(r). This yields Kepler ellipses (gravity), Rutherford hyperbolas (Coulomb scattering), or the hydrogen energy ladder (quantum Coulomb).
5. Map back to the real bodies. Once r(t) is known, recover each body: r₁ = R + (m₂/(m₁+m₂))r and r₂ = R − (m₁/(m₁+m₂))r. Both bodies orbit the common center of mass; the heavier one traces the smaller ellipse.

The reduced mass is always smaller than either mass

This is the single most useful sanity check. From 1/μ = 1/m₁ + 1/m₂, both terms on the right are positive, so 1/μ exceeds each of them individually — hence μ < m₁ and μ < m₂ simultaneously. Two limiting cases anchor the intuition:

Mass regime	Reduced mass μ	Physical meaning
Equal masses (m₁ = m₂ = m)	μ = m/2	Both bodies matter equally; μ is exactly half a single mass
One much heavier (m₁ ≫ m₂)	μ ≈ m₂ · (1 − m₂/m₁)	μ collapses onto the lighter mass; heavy body ≈ fixed
Proton + electron	μ ≈ 0.99946 mₑ	Nucleus nearly fixed; 0.054% correction
Sun + Earth	μ ≈ 0.999997 M_Earth	Sun essentially at the CM; ppm-level correction
Equal-mass binary stars	μ = M_star/2	Both stars orbit the midpoint; the correction is huge

Worked example: correcting the hydrogen spectrum

The Bohr–Schrödinger energy levels of a hydrogen-like atom are, with the finite-nuclear-mass correction built in through μ:

E_n = − (μ · e⁴) / (8 ε₀² h² n²)   =   − (μ/mₑ) · 13.6057 eV / n²

where e = 1.602 × 10⁻¹⁹ C is the elementary charge, ε₀ = 8.854 × 10⁻¹² F/m is the vacuum permittivity, h = 6.626 × 10⁻³⁴ J·s is Planck's constant, and n = 1, 2, 3, … is the principal quantum number. The bracketed 13.6057 eV is the ionization energy computed with the infinite-mass electron mass mₑ.

The proton is m_p = 1836.15 mₑ, so the reduced mass of the proton–electron system is

μ_H = m_p mₑ /(m_p + mₑ) = mₑ /(1 + mₑ/m_p) = mₑ · (1836.15/1837.15) = 0.99946 mₑ

Every level, and therefore every spectral line, is scaled down by this 0.99946 factor. In terms of the Rydberg constant that means the physically observed value for hydrogen is

R_H = R∞ · μ_H/mₑ = 10 973 731.6 m⁻¹ × 0.99946 = 10 967 758.3 m⁻¹

Without the reduced-mass correction the predicted Balmer and Lyman lines would sit at wavelengths off by roughly one part in 2000 — a discrepancy far larger than spectroscopic precision. The correction is not optional; it is the difference between a theory that matches the lab and one that does not.

History & the deuterium isotope shift

In 1931 Harold Urey reasoned that if a heavy isotope of hydrogen existed, its heavier nucleus would give a reduced mass closer to the bare electron mass, so every spectral line would shift to slightly shorter wavelength. The deuteron is roughly twice the proton mass (m_D ≈ 3670 mₑ), giving

μ_D/μ_H = [mₑ/(1+mₑ/m_D)] / [mₑ/(1+mₑ/m_p)] ≈ 1.000272

Since energies scale with μ, the deuterium Balmer-alpha line sits at a slightly higher energy — a wavelength shorter than hydrogen's Balmer-α (656.279 nm in air, 656.470 nm in vacuum) by about 0.179 nm. Urey photographed exactly this faint companion line in 1931 and won the 1934 Nobel Prize in Chemistry for discovering deuterium. The entire discovery hinged on the reduced-mass shift.

System	Nucleus mass	μ (in mₑ)	Balmer-α wavelength
Infinite-mass (ideal)	∞	1.000000	656.112 nm
Hydrogen (¹H)	1836.15 mₑ	0.999456	656.470 nm
Deuterium (²H)	3670.5 mₑ	0.999728	656.291 nm
Positronium (e⁺e⁻)	1 mₑ	0.500000	~1312 nm (½ scaling)

Positronium is the extreme case: an electron bound to its antiparticle, both with mass mₑ, gives μ = mₑ/2 exactly — halving every hydrogen energy and doubling every wavelength. It is the cleanest possible demonstration that reduced mass, not the electron mass, sets the scale.

Orbital periods and Kepler's third law

For two gravitating bodies the equivalent one-body problem has potential V(r) = −Gm₁m₂/r acting on the reduced mass μ. Solving it recovers Kepler's laws for the relative orbit, and the exact third law reads:

T² = 4π² a³ / [G (m₁ + m₂)]

where T is the orbital period (s), a is the semi-major axis of the relative orbit (m), and G = 6.674 × 10⁻¹¹ N·m²/kg². Notice the subtlety: the total mass (m₁+m₂) sets the period, while the reduced mass μ appears in the orbit's energy E = −Gm₁m₂/(2a) and angular momentum L = μ√(GM a(1−e²)). Treating the Sun as fixed uses M_Sun instead of M_Sun+m_planet; for Earth that is a 3-part-per-million error, but for a binary star of two equal masses it changes the period by a full factor of √2. Reduced mass and total mass play distinct, non-interchangeable roles.

Molecular vibrations and rotations

A diatomic molecule is two atoms on a chemical bond that behaves, near equilibrium, like a spring of stiffness k. Its vibration is a textbook one-body oscillator once you use the reduced mass of the two atoms:

ω = √(k/μ)     with     μ = m_A m_B / (m_A + m_B)

For a rigid rotor the same μ enters the moment of inertia I = μr₀², setting rotational line spacing. This is why isotope substitution is diagnostic: replacing hydrogen (1 u) with deuterium (2 u) in a C–H bond roughly doubles that atom's mass, raises μ, and drops the stretch frequency by a factor near √2 — the kinetic isotope effect that pervades infrared and Raman spectroscopy.

Common misconceptions

• "Reduced mass is the average of the two masses." No — it is the harmonic-style combination m₁m₂/(m₁+m₂), which is always smaller than either mass, never between them. For equal masses it equals m/2, half of each, not m.
• "Total mass and reduced mass are interchangeable in Kepler's law." No — the period depends on the total mass (m₁+m₂); the reduced mass μ governs the energy and angular momentum of the relative orbit. Swapping them silently is a classic error.
• "For hydrogen you can always just use the electron mass." It's a good 0.05% approximation, but the reduced-mass correction is essential for matching precision spectroscopy and is the whole basis of the hydrogen–deuterium isotope shift.
• "Reduced mass only matters in quantum mechanics." It is a purely classical bookkeeping device — it appears identically in Newtonian orbits, scattering cross-sections, and coupled oscillators. Quantum mechanics simply inherits it.
• "The center of mass is where the lighter body orbits around." Both bodies orbit the common center of mass; the heavier body traces the smaller ellipse. The reduced-mass particle at r describes their separation, not either body directly.

JavaScript — reduced-mass calculations

// Reduced mass of two bodies (same units in, same units out)
function reducedMass(m1, m2) {
  return (m1 * m2) / (m1 + m2);
}

// --- Atoms (masses in electron masses) ---
const m_e = 1;              // electron mass (reference)
const m_p = 1836.15267;     // proton  / electron mass ratio
const m_d = 3670.48296;     // deuteron / electron mass ratio

const mu_H = reducedMass(m_p, m_e);   // ~0.99946  mₑ
const mu_D = reducedMass(m_d, m_e);   // ~0.99973  mₑ
console.log(`μ_H = ${mu_H.toFixed(6)} mₑ`);      // 0.999456
console.log(`μ_D = ${mu_D.toFixed(6)} mₑ`);      // 0.999728

// Rydberg constant scales linearly with reduced mass
const R_inf = 10973731.568; // m⁻¹ (infinite-mass Rydberg)
const R_H = R_inf * mu_H;   // 10 967 758 m⁻¹
const R_D = R_inf * mu_D;   // 10 970 742 m⁻¹
console.log(`R_H = ${R_H.toFixed(1)} m⁻¹`);

// Hydrogen level energy: E_n = -(μ/mₑ)·13.6057 eV / n²
function levelEnergy(mu_over_me, n) {
  return -(mu_over_me) * 13.6057 / (n * n);
}
console.log(`E_1(H) = ${levelEnergy(mu_H, 1).toFixed(4)} eV`);  // -13.5983

// Balmer-alpha (n=3→2) wavelength via 1/λ = R(1/4 - 1/9)
function balmerAlpha(R) {
  const invLambda = R * (1/4 - 1/9);
  return 1e9 / invLambda; // nm
}
console.log(`H  Hα = ${balmerAlpha(R_H).toFixed(3)} nm`);  // 656.470 (vacuum)
console.log(`D  Hα = ${balmerAlpha(R_D).toFixed(3)} nm`);  // 656.291 -> ~0.18 nm isotope shift

// --- Diatomic vibration: ω = sqrt(k/μ) ---
const u = 1.66053907e-27;   // kg per atomic mass unit
function vibFreqHz(k, mA_u, mB_u) {
  const mu = reducedMass(mA_u, mB_u) * u;   // kg
  return Math.sqrt(k / mu) / (2 * Math.PI); // Hz
}
// CO bond: k ≈ 1860 N/m, masses 12 u and 16 u
console.log(`CO stretch ≈ ${(vibFreqHz(1860, 12, 16)/1e12).toFixed(1)} THz`); // ~64 THz

// --- Kepler period: total mass sets T, not μ ---
const G = 6.674e-11;
function orbitalPeriod(a, m1, m2) {
  return 2 * Math.PI * Math.sqrt(a**3 / (G * (m1 + m2)));
}

Where reduced mass shows up

• Celestial mechanics. Binary stars, exoplanet radial-velocity wobble, and the Earth–Moon barycenter all require μ and the total mass, not a fixed-center approximation.
• Atomic spectroscopy. Rydberg constant, hydrogen/deuterium/tritium isotope shifts, muonic and exotic atoms (a muon is 207 mₑ, so μ shifts dramatically).
• Molecular physics & chemistry. Vibrational and rotational spectra, kinetic isotope effects, zero-point energy of bonds.
• Scattering theory. Rutherford and quantum scattering cross-sections are computed in the CM frame with the reduced mass.
• Gravitational-wave astronomy. The "chirp mass" of a merging binary combines total and reduced mass; μ controls the orbital dynamics that shape the waveform.