Born–Oppenheimer Approximation¶
The Born–Oppenheimer approximation is one of the most important approximations in quantum chemistry. It forms the foundation of nearly all electronic structure methods implemented in Gaussian, including Hartree–Fock (HF), Density Functional Theory (DFT), Møller–Plesset perturbation theory (MP2), Coupled Cluster (CC), and many multireference methods.
The central idea behind the approximation is remarkably simple: nuclei are much heavier than electrons and therefore move much more slowly. Because of this large difference in mass, the electronic motion can be treated independently of the nuclear motion.
This approximation transforms an otherwise intractable many-particle quantum mechanical problem into one that can be solved efficiently on modern computers.
The Complete Molecular Schrödinger Equation¶
A molecule consists of
- Electrons
- Atomic nuclei
Both particles obey quantum mechanics and are described by a single molecular wavefunction,
where
- \(\mathbf r\) represents all electronic coordinates.
- \(\mathbf R\) represents all nuclear coordinates.
The complete molecular Schrödinger equation is
where
- \(\hat H\) is the molecular Hamiltonian,
- \(E\) is the total molecular energy.
In principle, solving this equation would provide the exact quantum mechanical description of the molecule.
Unfortunately, this equation is impossible to solve exactly for all but the simplest systems.
Why Is the Full Problem Difficult?¶
The molecular Hamiltonian contains interactions between
- every electron,
- every nucleus,
- every electron–electron pair,
- every electron–nucleus pair,
- every nucleus–nucleus pair.
As the number of atoms increases, the number of interactions grows rapidly, making an exact solution computationally impossible.
For example,
- Water contains 10 electrons.
- Benzene contains 42 electrons.
- Proteins may contain tens of thousands of electrons.
Directly solving the complete molecular Schrödinger equation for such systems is beyond the capabilities of even the largest supercomputers.
The Key Observation¶
The Born–Oppenheimer approximation exploits the enormous mass difference between electrons and nuclei.
| Particle | Approximate Mass |
|---|---|
| Electron | 1 |
| Proton | 1836 |
| Carbon nucleus | ≈ 22,000 |
| Oxygen nucleus | ≈ 29,000 |
Because nuclei are thousands of times heavier than electrons,
- electrons move extremely rapidly,
- nuclei move comparatively slowly.
The electrons therefore adjust almost instantaneously whenever the nuclei change position.
Freezing the Nuclei¶
Instead of solving for electrons and nuclei simultaneously, the Born–Oppenheimer approximation assumes that the nuclei remain fixed while solving the electronic problem.
Mathematically,
Nuclear Positions
│
▼
Treat as Constant
│
▼
Solve Electronic Schrödinger Equation
│
▼
Obtain Electronic Energy
For every fixed arrangement of nuclei, Gaussian computes the corresponding electronic energy.
Separation of the Wavefunction¶
The molecular wavefunction is approximated as the product of two independent wavefunctions,
where
- \(\psi_e\) describes the electronic motion,
- \(\chi_n\) describes the nuclear motion.
Notice that the electronic wavefunction still depends on the nuclear positions as parameters, but the nuclei are treated as fixed during the electronic calculation.
This separation dramatically simplifies the problem.
Electronic Schrödinger Equation¶
After applying the Born–Oppenheimer approximation, Gaussian solves only the electronic Schrödinger equation,
where
- \(\hat H_e\) is the electronic Hamiltonian,
- \(E_e\) is the electronic energy.
This is the equation solved in every Single Point Energy calculation.
Potential Energy Surface¶
Repeating the electronic calculation for many different nuclear geometries produces the Potential Energy Surface (PES).
The collection of these energies forms the multidimensional potential energy surface.
Different Gaussian calculations use this surface in different ways:
| Calculation | Use of the Potential Energy Surface |
|---|---|
| Single Point | One point on the surface |
| Optimization | Search for a minimum |
| Frequency | Evaluate the curvature near a stationary point |
| IRC | Follow the reaction path |
| Scan | Map selected regions of the surface |
Thus, the Born–Oppenheimer approximation provides the framework for all of these calculations.
Advantages of the Approximation¶
The Born–Oppenheimer approximation
- greatly simplifies the molecular Schrödinger equation,
- reduces computational cost,
- allows independent treatment of electronic and nuclear motion,
- enables efficient quantum chemical calculations for molecules containing hundreds of atoms.
Without this approximation, routine electronic structure calculations would be computationally impractical.
Limitations¶
Although highly successful, the Born–Oppenheimer approximation is not exact.
It becomes less accurate when
- electronic and nuclear motions become strongly coupled,
- electronic states are nearly degenerate,
- conical intersections are present,
- non-adiabatic transitions occur,
- ultrafast photochemical processes are studied.
In these situations, more advanced methods that explicitly couple electronic and nuclear motion are required.
Relevance to Gaussian¶
Every time a Single Point Energy calculation is performed in Gaussian,
- The nuclear coordinates are read from the input file.
- The nuclei are assumed to remain fixed.
- The electronic Schrödinger equation is solved.
- The electronic energy is computed.
Because the nuclei do not move during a Single Point Energy calculation, the Born–Oppenheimer approximation is particularly well suited for this type of calculation.
Summary¶
The Born–Oppenheimer approximation is the fundamental approximation underlying modern quantum chemistry. By exploiting the large mass difference between electrons and nuclei, it separates the molecular problem into electronic and nuclear components. Gaussian uses this approximation to solve the electronic Schrödinger equation for fixed nuclear geometries, making calculations such as Single Point Energy, Geometry Optimization, Frequency Analysis, and Reaction Path calculations computationally feasible.
Next Section¶
Having separated the electronic and nuclear motions using the Born–Oppenheimer approximation, the next step is to examine the Electronic Schrödinger Equation, which describes the behavior of electrons in the electrostatic field created by the fixed nuclei.