Born–Oppenheimer in Molecules

The Born–Oppenheimer hierarchy turns the coupled quantum mechanics of electrons and nuclei into the working language of molecular structure. At each nuclear geometry, one solves an electronic eigenproblem. Its eigenvalues define potential-energy surfaces; quantum nuclear motion on one or several surfaces then produces vibration, rotation, tunneling, dissociation, and rovibronic spectra.

The approximation is not the claim that nuclei are classical or permanently fixed. “Clamped nuclei” names only the auxiliary electronic problem. A molecular calculation is complete only after nuclear motion, exchange symmetry, and the required corrections have been restored.

This page owns the molecular application:

• how a surface minimum becomes an equilibrium geometry;
• how a surface curvature becomes vibrational frequencies;
• how isotope masses enter nuclear motion;
• how rotations, vibrations, and electronic transitions arise;
• which diagonal and nonadiabatic corrections matter for spectroscopy;
• how to diagnose when one surface is inadequate.

Born–Oppenheimer Approximation as Scale Separation owns the general electronic-channel derivation, gauge covariance, gap estimates, and mathematical approximation logic. Molecular Hamiltonian owns the all-particle Coulomb operator and exact center-of-mass separation. Potential Energy Surfaces owns the geometry of minima and saddles, reaction coordinates, asymptotic channels, surface representations, and uncertainty.

Nonadiabatic Coupling owns molecular avoided-crossing dynamics, wavepacket branching, surface-hopping approximations, vibronic models, and photochemical applications.

Conical Intersections owns exact degeneracy conditions, branching-plane geometry, intersection seams, MECIs, and their molecular Berry-phase consequences.

Starting from the Internal Hamiltonian

After overall translation has been removed, write the field-free molecular Hamiltonian schematically as

$$H_{\mathrm{int}} = T_N+H_e(R).$$

Here:

• $R$ denotes a complete set of internal nuclear coordinates;
• $T_N$ is the corresponding nuclear kinetic operator;
• $H_e(R)$ acts on electronic coordinates while $R$ labels the nuclear geometry.

For transparent formulas, it is common to write

$$T_N = -\sum_A \frac{\hbar^2}{2M_A} \nabla_A^2.$$

In a rigorous internal-coordinate treatment, $T_N$ may also contain reduced masses, mass-polarization terms, rotational terms, and coordinate metrics. Those terms do not alter the scale-separation idea, but they matter in precision work.

Surface convention

Throughout this page, $H_e(R)$ includes internuclear repulsion:

$$\begin{aligned} H_e(R) ={}& -\frac12\sum_i\nabla_i^2 -\sum_{i,A}\frac{Z_A}{r_{iA}} \\ &+ \sum_{i<j}\frac{1}{r_{ij}} + \sum_{A<B}\frac{Z_AZ_B}{R_{AB}} \end{aligned}$$

in atomic units. Its eigenvalue $E_a(R)$ is therefore already a potential-energy surface. If another source excludes $V_{NN}$ from the electronic Hamiltonian, translate by

$$U_a(R) = \mathcal E_a(R)+V_{NN}(R).$$

The physics is identical when the convention is used consistently.

Fixed-Geometry Electronic States

At every admissible geometry, solve

$$H_e(R) \lvert\phi_a(R)\rangle = E_a(R) \lvert\phi_a(R)\rangle.$$

The electronic inner product holds $R$ fixed:

$$\langle\phi_a(R) \vert \phi_b(R)\rangle_e = \delta_{ab}.$$

Each state must have the correct electron antisymmetry and can be classified by spin and by the symmetries of the chosen nuclear geometry. The geometry is a parameter of this electronic calculation, not a measured classical configuration of the exact free molecule.

Geometry space has identifications

Two laboratory configurations related by an overall translation or rotation represent the same internal shape. Configurations related by exchange of identical nuclei may also represent the same physical arrangement, subject to the required permutation symmetry of the total state.

Consequently, a potential-energy surface is naturally a function on nuclear shape space, not on an arbitrary list of labeled laboratory coordinates. Cartesian coordinates remain convenient computationally, but zero modes and permutation-related copies must be handled correctly.

Electronic labels can change with geometry

Point-group irreducible representations classify states at geometries with that symmetry. Along a distortion, the point group may decrease and labels may correlate or mix. Energy ordering alone is not a reliable state label near avoided crossings. State tracking should use symmetry, overlaps, transition properties, or a carefully chosen diabatic representation.

Exact Channel Expansion and the One-Surface Step

A complete fixed-geometry electronic basis gives the exact expansion

$$\Psi(r,R) = \sum_a \chi_a(R)\phi_a(r;R).$$

The expansion itself is not the Born–Oppenheimer approximation. Nuclear derivatives act on both factors, producing diagonal and off-diagonal derivative couplings.

The simplest molecular ansatz retains one electronic channel:

$$\Psi(r,R) \approx \chi_a(R)\phi_a(r;R).$$

At leading order, the nuclear wavefunction obeys

$$\left[ T_N+E_a(R) \right] \chi_{a\nu}(R) = E_{a\nu} \chi_{a\nu}(R).$$

The index $\nu$ collectively labels vibration, rotation, tunneling, and other nuclear quantum numbers on surface $a$.

Where approximation enters

The leading one-surface equation neglects:

• coupling to other electronic surfaces;
• the geometry dependence of the electronic state inside $T_N$;
• finite-mass corrections beyond the chosen internal kinetic operator;
• relativistic, radiative, and nuclear-structure contributions unless added separately.

An accurate electronic energy at one geometry does not validate these omissions. The relevant object is the electronic eigenspace over the entire nuclear region sampled by $\chi_{a\nu}(R)$.

Potential-Energy Surfaces as Molecular Potentials

A surface $E_a(R)$ supplies the effective potential for nuclear motion in electronic state $a$.

Minima and equilibrium structures

An equilibrium geometry $R_e$ on a smooth surface satisfies

$$\left. \nabla_R E_a(R) \right|_{R_e} = 0.$$

For a stable minimum, the Hessian restricted to internal distortions is positive:

$$\left. \delta R^{\mathsf T} \nabla_R^2E_a \delta R \right|_{R_e} >0$$

for every nonzero vibrational displacement $\delta R$.

The subscript $e$ means equilibrium, not electron. $R_e$ is a parameter of the effective surface. It is not generally equal to a vibrationally averaged bond length inferred from a spectrum.

Gradients and forces

For an exact, normalized, nondegenerate electronic eigenstate,

$$\frac{\partial E_a}{\partial R_\lambda} = \left\langle \phi_a \left| \frac{\partial H_e}{\partial R_\lambda} \right| \phi_a \right\rangle_e.$$

This is the Hellmann–Feynman relation. In a finite geometry-dependent basis, additional Pulay terms can appear because the basis itself changes with $R$. A numerically converged energy does not automatically imply a comparably converged gradient.

Saddles and asymptotes

A first-order saddle has one unstable internal Hessian direction and often organizes a reaction barrier. At large fragment separation, the surface must approach the correct separated-fragment channel. Surface construction, stationary points, reaction paths, and global representations deserve their own treatment; here they serve as inputs to nuclear quantum mechanics.

Small Oscillations of a Diatomic

For a diatomic molecule on an isolated electronic surface, let $R$ be the internuclear distance and $\mu_N$ the nuclear reduced mass:

H₂⁺ Ion supplies the simplest explicit example: one solves the one-electron two-center equation at each $R$, adds proton–proton repulsion, and then treats the resulting curve as the potential for nuclear vibration and rotation.

$$\mu_N = \frac{M_AM_B}{M_A+M_B}.$$

Near a stable minimum,

$$E_a(R) \approx E_a(R_e) + \frac12 k_a(R-R_e)^2,$$

where

$$k_a = \left. \frac{d^2E_a}{dR^2} \right|_{R_e}.$$

The leading harmonic frequency is

$$\omega_e = \sqrt{\frac{k_a}{\mu_N}}.$$

This relation separates electronic and isotopic ingredients:

• $k_a$ comes primarily from the electronic surface;
• $\mu_N$ comes from the isotopes;
• the zero-point width and level spacing depend on both.

The harmonic levels are

$$E_v \approx E_a(R_e) + \hbar\omega_e \left( v+\frac12 \right).$$

Anharmonicity, rotation, diagonal corrections, and surface coupling shift these levels. The exact oscillator solution is developed in Quantum Harmonic Oscillator.

Vibrations of Diatomics owns the molecular radial problem beyond this reduction, including term-value conventions, isotope scaling, Morse-model limits, infrared selection rules, and measured level progressions.

Isotope shift

If the electronic surface is held fixed while $\mu_N$ changes,

$$\omega_e\propto\mu_N^{-1/2}.$$

This leading isotope law is powerful but incomplete. High-resolution isotope shifts also probe adiabatic corrections, nonadiabatic effective masses, nuclear radii, and hyperfine interactions.

Diatomic Rotation and Vibration

For an isolated electronic $\Sigma$ state with neglected spin and electronic–rotational coupling, the radial nuclear equation is

$$\begin{aligned} \Bigg[ &- \frac{\hbar^2}{2\mu_N} \frac{d^2}{dR^2} + \frac{\hbar^2J(J+1)} {2\mu_NR^2} \\ &+ E_a(R) \Bigg] u_{vJ}(R) = E_{vJ} u_{vJ}(R). \end{aligned}$$

Here $u_{vJ}(R)$ is the reduced radial nuclear wavefunction. The centrifugal term couples rotation to vibration because $R$ is not fixed.

Near $R_e$, the rigid-rotor constant is

$$B_e = \frac{\hbar^2} {2\mu_NR_e^2}.$$

The leading product formula

$$E_{vJ} \approx E_v+B_eJ(J+1)$$

neglects centrifugal distortion and the variation of the rotational constant across the vibrational wavefunction. A better state-dependent constant samples

$$B_v \propto \left\langle \frac{1}{R^2} \right\rangle_v.$$

For electronic states with nonzero orbital angular momentum about the molecular axis, open shells, or strong spin couplings, the simple $J(J+1)$ term must be replaced by the appropriate coupled angular-momentum Hamiltonian. Rigid Rotor owns the canonical model, while Molecular Rotation Applications develops rotational tensor selection rules.

Polyatomic Normal Modes

For a nonlinear molecule with $N_N$ nuclei, expand the surface around $R_e$ in Cartesian displacements $q_{A\alpha}$:

$$E_a(R) \approx E_a(R_e) + \frac12 \mathbf q^{\mathsf T} K\mathbf q.$$

Define the mass-weighted Hessian

$$F_{A\alpha,B\beta} = \frac{ K_{A\alpha,B\beta} }{ \sqrt{M_AM_B} }.$$

After removing overall translation and rotation, solve

$$\sum_{B\beta} F_{A\alpha,B\beta} L_{B\beta,k} = \omega_k^2 L_{A\alpha,k}.$$

The nonzero eigenvalues give harmonic frequencies $\omega_k$. There are:

$$3N_N-6$$

vibrational modes for a nonlinear equilibrium geometry and

$$3N_N-5$$

for a linear one.

Normal modes are local coordinates near one minimum. Large-amplitude torsion, inversion, proton transfer, and dissociation may require curvilinear coordinates or a global surface. Imaginary harmonic frequency at a stationary point means a negative Hessian eigenvalue, not an exponentially growing quantum energy.

Normal Modes of Polyatomics develops the full mass-matrix formulation, translation and rotation projection, point-group labels, infrared and Raman activity, and computational validation.

Molecular Structure Is an Effective Statement

The surface minimum $R_e$ is the most common theoretical definition of equilibrium structure. Several experimentally reported structures answer different questions:

Symbol or description	Meaning
$r_e$	Minimum of an effective electronic potential surface
Ground-state average	Expectation or effective average over the $v=0$ nuclear distribution
Isotopic substitution structure	Geometry inferred from changes in rotational constants
Effective spectroscopic structure	Parameters fitted within a chosen rovibrational Hamiltonian

These need not coincide. Zero-point motion is isotope dependent, and fitted parameters can absorb higher-order corrections.

An exact field-free molecular eigenstate also has definite total angular momentum and is not pinned to one laboratory orientation. A body-fixed equilibrium geometry is an effective internal description that becomes observable through rotationally invariant quantities, aligned preparations, environmental localization, or conditional measurements.

Electronic Excitation and Vibronic Bands

Different electronic states have different surfaces and nuclear eigenfunctions:

$$\left[ T_N+E_a(R) \right] \chi_{a\nu}(R) = E_{a\nu}\chi_{a\nu}(R).$$

For a transition $a,\nu\to b,\nu'$, the Born–Oppenheimer transition amplitude contains

$$\mathcal M_{b\nu',a\nu} = \int dR\, \chi_{b\nu'}^*(R) \boldsymbol\mu_{ba}(R) \chi_{a\nu}(R),$$

where

$$\boldsymbol\mu_{ba}(R) = \langle \phi_b(R) \vert \widehat{\boldsymbol\mu} \vert \phi_a(R) \rangle_e.$$

If the electronic transition dipole varies slowly over the nuclear wavefunctions, the Condon approximation replaces it by a representative constant. The vibrational intensity is then governed by the Franck–Condon overlap

$$S_{\nu'\nu} = \left| \int dR\, \chi_{b\nu'}^*(R) \chi_{a\nu}(R) \right|^2.$$

The common vertical-transition picture reflects that the light–matter interaction acts on a time scale short compared with substantial nuclear displacement. It does not mean nuclear kinetic energy is zero before or after the transition.

Diagonal Adiabatic Correction

Even when one electronic channel is retained, $T_N$ differentiates the geometry-dependent electronic state. Define the diagonal connection

$$\mathbf d_{aa}^{A} = \langle \phi_a \vert \nabla_A\phi_a \rangle_e.$$

A gauge-invariant scalar correction is obtained from the part of each electronic derivative orthogonal to the retained state. Define

$$\begin{aligned} Q_a &= 1-\lvert\phi_a\rangle \langle\phi_a\rvert, \\ \lvert\nabla_A^\perp\phi_a\rangle &= Q_a \lvert\nabla_A\phi_a\rangle. \end{aligned}$$

Then

$$\Phi_a(R) = \sum_A \frac{\hbar^2}{2M_A} \left\| \nabla_A^\perp\phi_a \right\|_e^2.$$

Locally, for a nondegenerate real electronic state with

$$\mathbf d_{aa}^{A}=0,$$

this reduces to the familiar diagonal Born–Oppenheimer correction:

$$E_{\mathrm{DBOC}}(R) = \sum_A \frac{\hbar^2}{2M_A} \langle \nabla_A\phi_a \vert \nabla_A\phi_a \rangle_e.$$

The correction is geometry and isotope dependent. It modifies equilibrium structures, barriers, and vibrational levels. It is often small on a chemical-energy scale but relevant at spectroscopic accuracy.

Translational consistency

Electronic derivatives with respect to laboratory nuclear coordinates can contain apparent center-of-mass contributions. A rigorous implementation must either separate translation first or use a formulation known to reproduce the same internal correction. The simple appearance of the DBOC formula does not license mixing laboratory and body-fixed derivatives.

Mass conventions

The exact all-particle Hamiltonian contains bare nuclear masses. Leading Born–Oppenheimer nuclear motion formally starts from those masses. Practical high-accuracy rovibrational Hamiltonians may use atomic masses, geometry-dependent effective masses, or explicit nonadiabatic mass tensors to represent electronic inertia.

These choices are not interchangeable bookkeeping tricks. The masses must be matched to the correction terms and to the dissociation threshold convention.

Berry Connection in One Surface

Define

$$\boldsymbol{\mathcal A}_a^A(R) = i \langle \phi_a(R) \vert \nabla_A\phi_a(R) \rangle_e.$$

The gauge-covariant one-surface nuclear Hamiltonian can be written schematically as

$$\begin{aligned} H_{N,a} ={}& \sum_A \frac{ \left[ -i\hbar\nabla_A -\hbar\boldsymbol{\mathcal A}_a^A \right]^2 }{2M_A} \\ &+ E_a(R) + \Phi_a(R). \end{aligned}$$

Under

$$\phi_a \longrightarrow e^{i\theta(R)}\phi_a,$$

the nuclear factor transforms as

$$\chi_a \longrightarrow e^{-i\theta(R)}\chi_a,$$

so the total molecular wavefunction is unchanged.

For a simply connected region with an isolated real electronic state, one can often choose a local gauge with $\boldsymbol{\mathcal A}=0$. Around a conical intersection, a globally single-valued real gauge may not exist. The nuclear wavefunction then acquires a geometric phase. Born–Oppenheimer Berry Phase owns that topology.

Off-Diagonal Nonadiabatic Coupling

Between different electronic channels,

$$\mathbf d_{ab}^{A}(R) = \langle \phi_a(R) \vert \nabla_A\phi_b(R) \rangle_e.$$

For nondegenerate states and an exact differentiable eigenproblem,

$$\mathbf d_{ab}^{A} = \frac{ \left\langle \phi_a \left| \nabla_AH_e \right| \phi_b \right\rangle_e }{ E_b-E_a }, \qquad a\ne b.$$

This relation exposes the central molecular diagnostic:

$$\begin{gathered} \text{small electronic gap} \\ \Downarrow \\ \text{potentially large channel coupling} \end{gathered}$$

Symmetry can force the numerator to vanish, so the gap alone is not sufficient. Conversely, a modest coupling integrated over a long trajectory or amplified by resonance can still produce substantial population transfer.

What nonadiabatic corrections do

Beyond a scalar DBOC, nonadiabatic theory can generate:

• transitions between electronic surfaces;
• geometry-dependent vibrational mass corrections;
• rotational mass corrections;
• level shifts and predissociation widths;
• intensity borrowing;
• effective vector potentials and geometric phases.

No single number called “the Born–Oppenheimer error” describes all these effects.

Approximation Hierarchy

Level	Molecular content retained
Clamped-nuclei electronic problem	$E_a(R)$ and $\phi_a(r;R)$ at specified geometries
Leading one-surface Born–Oppenheimer	Quantum nuclear motion on one bare surface
Adiabatic one-surface	Adds diagonal scalar and, when needed, geometric connection terms
Nonadiabatic perturbation theory	Adds effective masses, scalar corrections, and controlled channel effects
Coupled-surface dynamics	Retains several electronic channels and their derivative couplings
All-particle non-BO calculation	Solves electrons and nuclei together after exact translation removal

The appropriate level depends on the observable. A surface may predict an equilibrium geometry well while missing a weak transition intensity, an isotope shift, or a predissociation lifetime.

Where One Surface Works Well

The one-surface picture is usually strongest when:

• the nuclear wavefunction is localized near a smooth minimum;
• the electronic state is nondegenerate and separated by a robust gap;
• nuclear kinetic energies are small compared with relevant electronic gaps;
• the target observable is not unusually sensitive to a small coupling;
• dissociation and reaction paths do not encounter near-degeneracies.

This regime underlies much of ground-state molecular structure, thermochemistry, rotational spectroscopy, and vibrational spectroscopy.

Precision hydrogen as a layered success

Molecular hydrogen illustrates both the success and the limits of the hierarchy. A highly accurate Born–Oppenheimer surface provides the leading energy. Adiabatic, nonadiabatic, relativistic, radiative, and finite-size corrections are then added with a controlled ledger. Agreement with precision spectra tests the entire layered calculation, not the bare surface alone.

Water and polyatomic spectroscopy

High-accuracy water spectra likewise require more than an electronic minimum: a global surface, variational nuclear motion, an accurate dipole-moment surface, and small adiabatic, nonadiabatic, relativistic, and radiative corrections. The Born–Oppenheimer organization remains indispensable even when empirical refinement is used to reach the final spectroscopic accuracy.

Where One Surface Fails

Conical intersections

Two same-spin electronic surfaces can become degenerate along a conical-intersection seam in polyatomic nuclear space. Derivative couplings become singular in an adiabatic basis, electronic and nuclear motion exchange energy efficiently, and a geometric phase can alter interference.

Avoided crossings

States of the same symmetry generally repel along a one-dimensional coordinate. If the nuclear wavepacket traverses a narrow avoided crossing, population transfer can be appreciable even for heavy nuclei.

Dissociation and rearrangement

At large separation, electronic gaps can close or several fragment arrangements can compete. A single state label followed from the equilibrium geometry may not correlate smoothly with the physically relevant asymptotic channel.

Light-particle transfer

Hydrogen transfer, proton-coupled electron transfer, and systems containing muons or positrons can weaken a naive heavy–light hierarchy. Tunneling may also sample geometries far from the minimum where gaps are smaller.

Excited-state photochemistry

Photoexcited wavepackets often reach intersections on femtosecond time scales. Internal conversion and intersystem crossing then require several electronic states and, for the latter, spin–orbit coupling beyond the spin-free Hamiltonian.

Near-threshold and resonance observables

A tiny energy correction can strongly shift a weakly bound level, scattering length, or predissociation width. “Small compared with an electronic energy” is not the right accuracy test for a near-threshold observable.

Choosing a Molecular Treatment

Use a one-surface calculation only after answering:

1. Which electronic state or isolated electronic subspace is being followed?
2. Over what nuclear region does the target wavefunction have support?
3. What is the minimum relevant electronic gap in that region?
4. Are derivative couplings symmetry allowed?
5. Does the observable depend on energy, phase, intensity, or lifetime?
6. Is the state close to a dissociation threshold or resonance?
7. Are isotope effects or spectroscopic precision required?
8. Does a global smooth electronic gauge exist?

When the answers are unfavorable, options include:

• retain several adiabatic surfaces;
• transform to a diabatic or quasidiabatic basis;
• use nonadiabatic perturbation theory;
• introduce an effective mass tensor;
• solve a reduced all-particle problem;
• use exact factorization as an interpretive or computational framework.

Nuclear Permutation Symmetry

Identical nuclei remain identical after electronic separation. Suppose $P_{AB}$ exchanges two identical nuclei. The electronic factor can transform nontrivially:

$$P_{AB}\phi_a(r;R) = \lambda_a \phi_a(r;R),$$

with $\lambda_a=\pm1$ in a one-dimensional symmetry sector. The nuclear spatial and spin factors must compensate so that the total state has the bosonic or fermionic symmetry required by the isotope.

This bookkeeping produces ortho and para species and restricts allowed rotational levels. Treating nuclei as permanently labeled classical centers erases those restrictions and can give the wrong statistical weights or selection rules.

From Surfaces to Observables

A complete Born–Oppenheimer-based prediction is a pipeline:

$$\begin{gathered} \text{electronic energies} \\ \text{and wavefunctions} \\ \downarrow \\ \text{surfaces, gradients,} \\ \text{couplings, and moments} \\ \downarrow \\ \text{quantum nuclear states} \\ \text{or dynamics} \\ \downarrow \\ \text{frequencies, intensities,} \\ \text{rates, and thresholds} \end{gathered}$$

Different observables need different electronic data:

• equilibrium structures require accurate gradients;
• harmonic frequencies require Hessians;
• vibrational spectra require a surface over the sampled region;
• intensities require dipole or higher-multipole surfaces;
• nonadiabatic dynamics requires several surfaces and coupling vectors;
• dissociation requires correct asymptotes;
• isotope shifts require a mass-consistent correction model.

An accurate energy table alone is not a molecular dynamics model.

Common Misconceptions

“The nuclei do not move”

They do not move inside the auxiliary electronic eigenproblem. Nuclear motion returns as a quantum problem on one or several electronic surfaces.

“The molecular wavefunction is exactly a product”

A single product is an approximation. A complete sum over electronic channels is exact, and nuclear derivatives couple its terms.

“A small mass ratio guarantees validity”

Only together with electronic gaps, coupling matrix elements, nuclear velocities, and the target time or energy scale.

“The potential-energy surface is directly observable”

The surface is an effective theoretical object. Spectra, scattering probabilities, structures inferred under stated conventions, and transition rates are observables.

“The equilibrium bond length is the measured bond length”

$r_e$ is a surface minimum. Ground-state averages and spectroscopically fitted structures include nuclear motion and can be isotope dependent.

“Adding the DBOC makes the result exact”

The DBOC is one diagonal correction. Off-diagonal coupling, effective masses, relativity, QED, and nuclear effects can remain.

“Adiabatic means the same thing everywhere”

In molecular work, “Born–Oppenheimer,” “adiabatic,” and “Born–Huang” are used with varying conventions. State explicitly whether the bare surface, diagonal correction, geometric connection, or coupled channels are retained.

“One can sort electronic states by energy at every geometry”

Energy ordering can swap at avoided crossings or become ambiguous at degeneracies. State tracking requires wavefunction information and symmetry.

Computational Workflow

For a molecule-specific Born–Oppenheimer calculation:

1. Specify isotopes, charge, spin multiplicity, and electronic symmetry.
2. Remove overall translation and declare the internal coordinate convention.
3. Choose the electronic structure method and basis appropriate to the target accuracy.
4. Compute energies over the nuclear region sampled by the intended states or dynamics.
5. Check gradients, Hessians, asymptotes, and state continuity.
6. Locate small gaps and evaluate or estimate nonadiabatic couplings.
7. Construct the nuclear kinetic operator with a stated mass convention.
8. Solve nuclear motion variationally or propagate it dynamically.
9. Add DBOC, nonadiabatic, relativistic, radiative, and nuclear terms as required.
10. Compute transition moments or scattering observables with compatible wavefunctions.
11. Compare isotopologues and limiting cases.
12. Report uncertainty separately for the surface, nuclear solver, and omitted physics.

Exercises

Exercise 1: Surface convention

One source defines

$$H_{\mathrm{el}}\phi_a = \mathcal E_a\phi_a$$

without $V_{NN}$. Another includes $V_{NN}$ in $H_e$ and reports $E_a$. Relate their nuclear Hamiltonians.

Solution

Because $V_{NN}(R)$ is a scalar in the electronic problem,

$$E_a(R) = \mathcal E_a(R)+V_{NN}(R).$$

The two nuclear Hamiltonians are

$$T_N+\mathcal E_a(R)+V_{NN}(R)$$

and

$$T_N+E_a(R),$$

which are identical. Adding $V_{NN}$ to $E_a$ a second time would double count nuclear repulsion.

Exercise 2: Diatomic isotope shift

Assume two isotopologues share the same harmonic force constant $k$. Show that their vibrational-frequency ratio is determined by their reduced masses.

Solution

For isotope $s$,

$$\omega_s = \sqrt{\frac{k}{\mu_s}}.$$

Therefore

$$\frac{\omega_2}{\omega_1} = \sqrt{\frac{\mu_1}{\mu_2}}.$$

The heavier isotopologue has the smaller vibrational spacing. Deviations from this leading relation diagnose mass-dependent corrections or a breakdown of the common-surface approximation.

Exercise 3: Rotation–vibration coupling

Explain why replacing $R$ by $R_e$ in the centrifugal term suppresses rotation–vibration coupling.

Solution

The exact one-surface centrifugal contribution is

$$\frac{\hbar^2J(J+1)} {2\mu_NR^2}.$$

Because $R$ is an operator, its expectation value depends on the vibrational state. Replacing it by the constant $R_e$ gives

$$B_eJ(J+1),$$

which is independent of vibration. The rigid-rotor approximation therefore removes the variation of rotational energy across the vibrational wavefunction and misses centrifugal distortion.

Exercise 4: Gauge-invariant diagonal correction

Show that

$$\left\langle \nabla\phi \left| 1-\lvert\phi\rangle\langle\phi\rvert \right| \nabla\phi \right\rangle$$

is invariant under $\phi\to e^{i\theta(R)}\phi$.

Solution

The derivative transforms as

$$\nabla\phi' = e^{i\theta} \left( \nabla\phi +i(\nabla\theta)\phi \right).$$

The projector

$$Q=1-\lvert\phi\rangle\langle\phi\rvert$$

annihilates $\phi$. Hence

$$Q\nabla\phi' = e^{i\theta}Q\nabla\phi.$$

Taking the norm removes the phase, so the projected derivative norm and $\Phi_a(R)$ are gauge invariant.

Exercise 5: Harmonic normal modes

Why must three translational and, for a nonlinear molecule, three rotational zero modes be removed from a Cartesian Hessian before interpreting its eigenvalues as vibrational frequencies?

Solution

The field-free surface depends only on internal shape. A common translation or infinitesimal rotation does not change its energy, so those displacements have zero restoring force. They produce null directions of the exact Cartesian Hessian. After projecting them out, a nonlinear $N_N$-nucleus molecule has

$$3N_N-6$$

internal vibrational directions. Small nonzero “zero-mode frequencies” in a numerical calculation indicate incomplete optimization, finite precision, or an inconsistent projection.

Exercise 6: Vertical excitation

Starting from the Born–Oppenheimer transition amplitude, state the additional assumption needed to reduce intensity to a Franck–Condon factor.

Solution

The amplitude is

$$\mathcal M = \int dR\, \chi_{b\nu'}^*(R) \boldsymbol\mu_{ba}(R) \chi_{a\nu}(R).$$

The Condon approximation assumes

$$\boldsymbol\mu_{ba}(R) \approx \boldsymbol\mu_{ba}(R_0)$$

over the nuclear region carrying appreciable amplitude. The constant electronic factor then leaves the overlap

$$\int dR\, \chi_{b\nu'}^*(R) \chi_{a\nu}(R),$$

whose squared magnitude is the Franck–Condon factor. Rapid geometry dependence of the transition moment produces non-Condon corrections.

Exercise 7: Diagnose an avoided crossing

A wavepacket approaches a region where $E_b-E_a$ decreases by two orders of magnitude while $\langle\phi_a|\nabla_RH_e|\phi_b\rangle$ remains comparable. What happens to the derivative-coupling estimate?

Solution

For nondegenerate states,

$$d_{ab} \sim \frac{ \langle\phi_a| \nabla_RH_e |\phi_b\rangle }{ E_b-E_a }.$$

If the numerator is unchanged and the gap decreases by a factor of $100$, the estimated coupling grows by a factor of $100$. A one-surface treatment becomes less credible. One should inspect the wavepacket velocity and symmetry, then use coupled surfaces or a suitable diabatic representation if transfer is appreciable.

Exercise 8: Identical nuclei

Two identical spin-$1/2$ nuclei are exchanged. If the electronic-vibrational-rotational spatial factor is symmetric, what symmetry must their nuclear-spin factor have?

Solution

Spin-$1/2$ nuclei are fermions, so the complete molecular state must be antisymmetric under their exchange. If the electronic-vibrational-rotational spatial factor contributes $+1$, the nuclear-spin factor must contribute $-1$. It must therefore lie in the antisymmetric nuclear-spin sector. For two spin-$1/2$ particles, that sector is the singlet.

Summary

• Fixed-geometry electronic eigenvalues become potentials for quantum nuclear motion.
• A surface minimum defines $R_e$; measured effective structures also contain zero-point and model-dependent effects.
• Surface curvature and nuclear masses determine leading vibrational frequencies, while moments of inertia determine rotational scales.
• The diagonal correction, Berry connection, effective masses, and off-diagonal couplings arise because electronic states vary with geometry.
• Small electronic gaps, conical intersections, dissociation, and near-threshold observables demand special care.
• Precision predictions are layered calculations combining electronic, nuclear, nonadiabatic, relativistic, radiative, and nuclear contributions.

Connections

• Common Molecular Hamiltonians
• Molecular Hamiltonian
• Molecular Quantum Mechanics
• Nonadiabatic Coupling
• Conical Intersections
• Potential Energy Surfaces
• Electronic Structure Overview
• Born–Oppenheimer Approximation as Scale Separation
• Born–Oppenheimer Berry Phase
• Molecular Physics
• Quantum Harmonic Oscillator
• Rigid Rotor
• Rotations of Molecules
• Vibrations of Diatomics
• Normal Modes of Polyatomics
• Rovibrational Coupling
• Electronic Spectroscopy
• Molecular Rotation Applications
• Quantum Chemistry Roadmap
• Quantum Chemistry References

References

1. M. Born and R. Oppenheimer, “Zur Quantentheorie der Molekeln,” Annalen der Physik 389, 457–484 (1927), doi:10.1002/andp.19273892002.
2. M. Born and K. Huang, Dynamical Theory of Crystal Lattices (Clarendon Press, 1954), Chapters IV–V.
3. H. C. Longuet-Higgins, U. Öpik, M. H. L. Pryce, and R. A. Sack, “Studies of the Jahn–Teller Effect. II. The Dynamical Problem,” Proceedings of the Royal Society A 244, 1–16 (1958), doi:10.1098/rspa.1958.0022.
4. C. A. Mead and D. G. Truhlar, “On the Determination of Born–Oppenheimer Nuclear Motion Wave Functions Including Complications Due to Conical Intersections and Identical Nuclei,” Journal of Chemical Physics 70, 2284–2296 (1979), doi:10.1063/1.437734.
5. N. C. Handy and A. M. Lee, “The Adiabatic Approximation,” Chemical Physics Letters 252, 425–430 (1996), doi:10.1016/0009-2614(96)00171-6.
6. W. Kutzelnigg, “The Adiabatic Approximation I. The Physical Background of the Born–Handy Ansatz,” Molecular Physics 90, 909–916 (1997), doi:10.1080/002689797171904.
7. B. T. Sutcliffe and R. G. Woolley, “Molecular Structure Calculations without Clamping the Nuclei,” Physical Chemistry Chemical Physics 7, 3664–3676 (2005), doi:10.1039/B509723C.
8. G. A. Hagedorn and A. Joye, “A Time-Dependent Born–Oppenheimer Approximation with Exponentially Small Error Estimates,” Communications in Mathematical Physics 223, 583–626 (2001), doi:10.1007/s002200100562.
9. K. Pachucki and J. Komasa, “Nonadiabatic Corrections to Rovibrational Levels of $\mathrm H_2$,” Journal of Chemical Physics 130, 164113 (2009), doi:10.1063/1.3114680.
10. K. Pachucki and J. Komasa, “Accurate Adiabatic Correction in the Hydrogen Molecule,” Journal of Chemical Physics 141, 224103 (2014), doi:10.1063/1.4902981.
11. K. Pachucki and J. Komasa, “Leading Order Nonadiabatic Corrections to Rovibrational Levels of $\mathrm H_2$, $\mathrm D_2$, and $\mathrm T_2$,” Journal of Chemical Physics 143, 034111 (2015), doi:10.1063/1.4927079.
12. J. Komasa, M. Puchalski, P. Czachorowski, G. Łach, and K. Pachucki, “Rovibrational Energy Levels of the Hydrogen Molecule through Nonadiabatic Perturbation Theory,” Physical Review A 100, 032519 (2019), doi:10.1103/PhysRevA.100.032519.
13. J. Tennyson, P. Barletta, M. A. Kostin, O. L. Polyansky, N. F. Zobov, and S. V. Shirin, “Ab Initio Rotation–Vibration Energy Levels of Triatomics to Spectroscopic Accuracy,” Spectrochimica Acta Part A 58, 663–672 (2002), doi:10.1016/S1386-1425(01)00663-1.
14. J. Tennyson, M. A. Kostin, P. Barletta, G. J. Harris, O. L. Polyansky, J. Ramanlal, and N. F. Zobov, “High-Accuracy Ab Initio Rotation–Vibration Transitions for Water,” Science 299, 539–542 (2003), doi:10.1126/science.1079558.
15. O. L. Polyansky, A. A. Kyuberis, L. Lodi, J. Tennyson, S. N. Yurchenko, and N. F. Zobov, “Calculation of Rotation–Vibration Energy Levels of the Water Molecule with Near-Experimental Accuracy Based on an Ab Initio Potential Energy Surface,” Journal of Physical Chemistry A 117, 9633–9643 (2013), doi:10.1021/jp312343z.
16. P. R. Bunker and P. Jensen, Molecular Symmetry and Spectroscopy, 2nd ed. (NRC Research Press, 1998).
17. H. Lefebvre-Brion and R. W. Field, The Spectra and Dynamics of Diatomic Molecules, revised ed. (Elsevier, 2004).
18. T. Helgaker, P. Jørgensen, and J. Olsen, Molecular Electronic-Structure Theory (Wiley, 2000).