Basis Sets

An electronic basis set specifies the one-particle functions from which orbitals and many-electron states are constructed. Its flexibility controls what a calculation can represent before a correlation method is chosen. This page concerns basis construction; the general logic of refinement and extrapolation is in basis-set convergence.

Helpful background. Rayleigh–Ritz explains the overlap metric of a nonorthogonal basis.

Localized Gaussian functions and contractions

A Cartesian primitive Gaussian centered at $\mathbf A$ has the form

$$g(\mathbf r)=N(x-A_x)^l(y-A_y)^m(z-A_z)^n e^{-\alpha|\mathbf r-\mathbf A|^2},\qquad \alpha>0.$$

Its angular degree is $L=l+m+n$. A contracted function is a prescribed linear combination of primitives. Distinguish primitive normalization from contraction normalization, and record whether coefficient tables assume normalized primitives. An omitted normalization changes all integrals consistently enough to evade some superficial tests.

The Gaussian product identity,

$$e^{-\alpha|\mathbf r-\mathbf A|^2} e^{-\beta|\mathbf r-\mathbf B|^2} =e^{-\mu|\mathbf A-\mathbf B|^2} e^{-(\alpha+\beta)|\mathbf r-\mathbf P|^2},$$

uses $\mu=\alpha\beta/(\alpha+\beta)$ and $\mathbf P=(\alpha\mathbf A+\beta\mathbf B)/(\alpha+\beta)$. It underlies efficient molecular integrals. Gaussian primitives are computationally convenient but do not individually reproduce a Coulomb cusp or an exponential long-range tail.

Multiple radial functions of the same angular type provide radial flexibility. Polarization functions add higher angular momenta, allowing the density to deform in a bond. Diffuse functions with small exponents describe extended tails, anions, and weak binding. These serve different purposes; increasing the number of tightly localized functions cannot repair a missing diffuse tail.

Representation conventions and conditioning

Cartesian $d$ shells have six components, while pure spherical $d$ shells have five. Basis names alone do not identify this choice, orbital ordering, or phase conventions. Match these details when comparing integral files or density matrices between programs.

Diffuse functions on neighboring centers can become almost linearly dependent. Diagonalize the overlap matrix and inspect its smallest eigenvalues before interpreting unstable energies. Removing directions below a threshold changes the effective basis; report and vary that threshold. A larger nominal basis with more removed functions need not contain the preceding effective trial space.

For a normalized $s$ primitive, $N=(2\alpha/\pi)^{3/4}$. Two normalized equal-exponent $s$ functions separated by $R$ have overlap

$$S(R)=e^{-\alpha R^2/2}.$$

For $\alpha=0.1$ and $R=1$ in consistent length units, $S\approx0.95123$. Their two-function overlap eigenvalues are approximately $1.95123$ and $0.04877$. This simple pair shows why diffuse neighboring functions approach dependence even when each function is correctly normalized.

Molecular and periodic convergence choices

For localized molecular bases, vary radial quality, polarization, and diffuse augmentation according to the observable. Correlation-consistent families provide an organized hierarchy, but asymptotic extrapolation still requires a tested model and matching treatment of core correlation.

Plane waves instead use a kinetic-energy cutoff and a cell. Their orthogonality simplifies overlap handling, while nuclear cusps generally require pseudopotentials or augmentation and independently converged auxiliary grids. A vacuum-cell test is essential for isolated molecules represented periodically.

In an interaction energy, a fragment can borrow basis functions from its partner. This basis-set superposition effect can exaggerate binding. A counterpoise calculation evaluates each isolated fragment in the combined basis, using partner functions without partner nuclei or electrons. It is a diagnostic and correction convention, not proof that all basis errors have vanished.

Exercises

Why does adding polarization functions not necessarily fix the energy of a weakly bound anion?

Solution

Polarization adds angular flexibility. A weakly bound electron also requires radial support far from the nuclei. If all radial exponents remain large, the basis may lack that tail regardless of angular degree. Diffuse augmentation and a test of the binding energy against the neutral threshold are needed.

References

• Dunning, T. H., Jr., “Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen,” Journal of Chemical Physics 90, 1007–1023 (1989), doi:10.1063/1.456153 — systematic correlation-consistent construction.
• Psi4 developers, “Basis Sets” documentation, accessed 13 October 2026 — orbital and auxiliary basis settings and representation conventions.
• Szabo, A., and N. S. Ostlund, Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory, Dover, 1996 — Gaussian bases, overlap, and electronic integral conventions.