Entanglement in Quantum Chemistry

Quantum chemistry studies many-electron states. Those states are not ordinary products of labeled particles: electrons are identical fermions, so fixed-$N$ electronic wavefunctions live in an antisymmetric subspace. Entanglement enters in several related but distinct ways:

• antisymmetry and exchange are built into the electronic state space;
• a single Slater determinant is the reference point for independent-particle approximations;
• correlated wavefunctions are superpositions of determinants or occupation configurations;
• orbital and mode entanglement quantify how strongly chosen orbitals are entangled with one another;
• spin entanglement appears between regions, centers, orbitals, or effective local moments;
• bond formation and bond breaking often reveal the limits of single-determinant pictures.

This page is a bridge. The canonical construction of Slater determinants, spin and spatial wavefunctions, and many-particle Hamiltonians lives elsewhere in this volume. Detailed Hartree–Fock theory, configuration interaction, coupled cluster, density-matrix renormalization group, active-space methods, and molecular applications belong in later atomic, molecular, and computational treatments.

The guiding distinction is:

$$\text{fermionic antisymmetry} \quad\neq\quad \text{electron correlation} \quad\neq\quad \text{orbital entanglement}.$$

Confusing these three is the fastest way to overclaim what an entanglement measure says about a molecule.

Electronic Hilbert Space

A one-electron spin-orbital depends on spatial and spin variables,

$$q=(\mathbf r,s),$$

and belongs to a one-particle Hilbert space $\mathcal h$. An $N$-electron wavefunction belongs to the antisymmetric sector

$$\wedge^N\mathcal h.$$

Equivalently, in a chosen spin-orbital basis, one may use fermionic Fock space and occupation-number states. A determinant occupying spin-orbitals $i_1,\ldots,i_N$ is represented by

$$a_{i_1}^\dagger a_{i_2}^\dagger \cdots a_{i_N}^\dagger \lvert0\rangle,$$

with a fixed ordering convention for signs.

The spin-orbital basis is not unique. Canonical molecular orbitals, localized orbitals, atomic orbitals, natural orbitals, and active-space orbitals can all be useful. Entanglement between orbitals is therefore a statement about a chosen mode decomposition, not an invariant property of electron labels.

Antisymmetry Is Not the Same as Correlation

A single Slater determinant is already antisymmetric. For two occupied spin-orbitals $a,b$,

$$\Psi(q_1,q_2) = \frac{1}{\sqrt2} \bigl[ \varphi_a(q_1)\varphi_b(q_2) - \varphi_b(q_1)\varphi_a(q_2) \bigr].$$

This state is not a product in formal particle slots, but those slots are not physical electron subsystems. In the occupation basis defined by $a,b$, the same determinant is simply

$$\lvert1_a,1_b\rangle.$$

It has definite occupations of those two spin-orbitals. The antisymmetric sign is essential physics, but it is not by itself the same as electron correlation in the quantum-chemistry sense.

Quantum chemistry usually reserves “correlation” for effects missing from the best single-determinant mean-field reference for a chosen electronic Hamiltonian. Exchange is already included in Hartree–Fock through antisymmetry. Dynamical and static correlation require more than one determinant, more flexible orbitals, or other correlated ansatz structure. Exchange and Correlation develops the atomic energy, pair-hole, and density-functional accounting behind this distinction.

Electronic Hamiltonian

In a spin-orbital basis, a standard nonrelativistic electronic Hamiltonian after the Born–Oppenheimer separation has the second-quantized form

$$H = \sum_{pq} h_{pq} a_p^\dagger a_q + \frac12 \sum_{pqrs} g_{pqrs} a_p^\dagger a_q^\dagger a_s a_r.$$

The one-electron integrals contain kinetic energy and electron-nuclear attraction. In atomic units, schematically,

$$h_{pq} = \int dq\, \varphi_p^*(q) \left( -\frac12\nabla^2 +v_{\mathrm{nuc}}(\mathbf r) \right) \varphi_q(q).$$

The two-electron integrals contain Coulomb repulsion:

$$g_{pqrs} = \int dq_1\,dq_2\, \varphi_p^*(q_1) \varphi_q^*(q_2) \frac{1}{r_{12}} \varphi_r(q_1) \varphi_s(q_2).$$

Details of integral symmetries, antisymmetrized integrals, basis-set convergence, relativistic corrections, nuclear motion, and effective core approximations belong to quantum-chemistry treatments. For this volume, the key point is structural: the electronic Hamiltonian acts on fermionic modes and couples their occupations through two-body interactions.

Hartree–Fock Reference

Hartree–Fock approximates the electronic ground state by one optimized Slater determinant,

$$\lvert\Phi_0\rangle = a_{i_1}^\dagger \cdots a_{i_N}^\dagger \lvert0\rangle.$$

Within that determinant, exchange symmetry is exact. The approximation lies in describing the many-electron state by independent occupied spin-orbitals moving in a self-consistent mean field.

For a fixed Hamiltonian and basis, the correlation energy is often defined as

$$E_{\mathrm{corr}} = E_0-E_{\mathrm{HF}},$$

where $E_0$ is the exact ground-state energy in that same model space. Variationally,

$$E_{\mathrm{corr}}\le0.$$

This number is not itself an entanglement measure. It is an energy difference. It can be related to correlation and multiconfigurational structure, but the relationship depends on the Hamiltonian, basis, state, and reference determinant.

Configuration Interaction Preview

A general correlated wavefunction in a finite spin-orbital basis can be expanded as a linear combination of determinants:

$$\lvert\Psi\rangle = \sum_I C_I\lvert\Phi_I\rangle.$$

Relative to a reference determinant $\lvert\Phi_0\rangle$, singly and doubly excited determinants have the schematic form

$$\lvert\Phi_i^a\rangle = a_a^\dagger a_i \lvert\Phi_0\rangle,$$

and

$$\lvert\Phi_{ij}^{ab}\rangle = a_a^\dagger a_b^\dagger a_j a_i \lvert\Phi_0\rangle.$$

Here $i,j$ usually label occupied spin-orbitals in the reference and $a,b$ virtual spin-orbitals. Full configuration interaction in a finite basis includes all determinants consistent with particle number and spin-orbital space. Truncated CI, coupled-cluster theory, perturbation theory, selected CI, and tensor-network methods choose different ways to approximate the same enormous correlated state space.

From the entanglement viewpoint, the important point is that determinant superpositions create nontrivial reduced states for orbitals, groups of orbitals, spins, or spatial regions. A one-determinant state may have zero occupation-mode entanglement in its own spin-orbital basis; a correlated superposition generally need not.

Orbital Entanglement

Orbital entanglement treats chosen orbitals or groups of orbitals as subsystems. For a spatial orbital $p$ with spin-up and spin-down modes, the local occupation basis is

$$\lvert0\rangle, \qquad \lvert\uparrow\rangle, \qquad \lvert\downarrow\rangle, \qquad \lvert\uparrow\downarrow\rangle.$$

Given a many-electron state $\rho$, the one-orbital reduced state is

$$\rho_p = \operatorname{Tr}_{\bar p}\rho,$$

and the one-orbital entropy is

$$s_p = -\operatorname{Tr}(\rho_p\log\rho_p).$$

For two orbitals $p,q$, define

$$\rho_{pq} = \operatorname{Tr}_{\overline{pq}}\rho, \qquad s_{pq} = -\operatorname{Tr}(\rho_{pq}\log\rho_{pq}).$$

The orbital mutual information is

$$I(p:q) = s_p+s_q-s_{pq}.$$

These quantities are useful diagnostics. They help identify strongly entangled active orbitals, choose orbital orderings in density-matrix renormalization group calculations, and distinguish nearly single-reference regions from multireference regions. They are also basis dependent. Localized orbitals and canonical orbitals can give different entanglement patterns for the same physical state.

Spin Entanglement

Spin entanglement in chemistry is rarely about labeled electrons. It is usually about spins associated with orbitals, sites, fragments, centers, or spatial regions. The two-electron singlet,

$$\lvert S\rangle = \frac{1}{\sqrt2} \left( \lvert\uparrow\downarrow\rangle - \lvert\downarrow\uparrow\rangle \right),$$

is maximally entangled across two distinguishable spin degrees of freedom when those spins are assigned to two physical modes or regions. For electrons, the complete state must also satisfy fermionic antisymmetry. Thus spin entanglement must be discussed together with the spatial or orbital part that identifies the parties.

This distinction matters in molecules. A closed-shell determinant can contain spin-singlet pairing due to the exchange-symmetry structure, while an effective model of two separated radical centers can contain operational spin entanglement between centers. The symbols may look similar, but the subsystem split is different.

Open-shell molecules add another layer: a single determinant may not be an eigenstate of total spin. Spin-adapted configuration state functions build proper spin symmetry from determinant combinations. That technology belongs to quantum chemistry; this page only records why spin entanglement and spin adaptation must be handled with exchange symmetry in view.

Bonding and Correlation

The hydrogen molecule is the standard warning against making single-determinant intuition too absolute. In a simple Heitler–London picture with localized orbitals $a$ and $b$, the covalent spatial factor has the schematic symmetric form

$$\psi_+(\mathbf r_1,\mathbf r_2) \propto a(\mathbf r_1)b(\mathbf r_2) + b(\mathbf r_1)a(\mathbf r_2),$$

paired with a spin singlet to make an antisymmetric two-electron state. Molecular-orbital pictures instead form bonding and antibonding orbitals and often start from a closed-shell determinant.

Near equilibrium, a single-reference description may be a good starting point. During bond stretching or bond breaking, however, near-degenerate configurations become important. A correlated state may require a superposition of determinants to avoid unphysical ionic contributions and to recover the correct separated-atom behavior.

This is a central chemistry use of entanglement diagnostics: large orbital entropies and strong orbital mutual informations often signal static correlation or multireference character. They do not replace chemical analysis, but they help expose where a one-determinant picture is structurally strained.

Active Spaces and Tensor Networks

In practical electronic-structure calculations, one often chooses an active space: a subset of orbitals treated with high-level correlation while the rest are kept inactive, weakly correlated, or perturbative. Entanglement diagnostics can help choose that active set.

If an orbital has large $s_p$, it is strongly entangled with the rest of the molecule in the chosen orbital basis. If $I(p:q)$ is large, orbitals $p$ and $q$ share strong total correlation. These facts can guide active-space construction, orbital localization, and ordering in tensor-network algorithms.

Matrix product state methods represent the many-electron coefficient tensor as a chain of orbital tensors. As in many-body physics, the required bond dimension is controlled by entanglement across orbital cuts:

$$S_A\le\log\chi.$$

The same warning applies here as in the many-body bridge: a tensor network is efficient when the chosen ordering and orbital basis keep the relevant entanglement manageable.

What This Page Owns

This page owns the conceptual bridge from composite-system entanglement to electronic-structure language:

• why electrons require antisymmetric state spaces;
• why a single determinant is not automatically a correlated wavefunction;
• how determinant expansions create orbital and spin reduced states;
• how one-orbital entropy and orbital mutual information are defined;
• why bond breaking and multireference character are entanglement-rich situations;
• why orbital entanglement is basis dependent.

It does not own:

• derivations of Hartree–Fock equations;
• Slater–Condon rules;
• full configuration-interaction algorithms;
• coupled-cluster theory;
• density-functional theory;
• complete active-space methods;
• quantum-chemistry DMRG implementations;
• chemical bonding theory beyond simple previews.

Those pages should link here for subsystem and entanglement language, but their canonical development belongs in electronic-structure and computational volumes.

Common Mistakes

• Calling every Slater determinant entangled because it is antisymmetric in particle slots.
• Confusing exchange with electron correlation.
• Treating correlation energy as an entanglement measure.
• Quoting orbital entanglement without stating the orbital basis and partition.
• Treating spin-up and spin-down electrons as different particle species.
• Ignoring spin adaptation when using determinant expansions for open-shell states.
• Assuming a low orbital entropy in one basis means the molecule has no correlation in every representation.
• Treating entanglement diagnostics as substitutes for convergence tests, basis checks, and chemical interpretation.

Cross-Links

• Slater Determinants
• Exchange and Correlation
• Spin and Spatial Wavefunctions
• Valence Bond Theory connects localized spin-coupled structures and resonance weights to basis-dependent correlation diagnostics.
• Identical-Particle Entanglement Cautions
• Occupation-Number Basis
• Mode Occupations
• Fermionic Fock Space
• One-Body Operators
• Two-Body Operators
• Two-Body Operators in Many-Body Models
• Many-Particle Hamiltonians
• Entanglement Entropy
• Mutual Information
• Entanglement in Many-Body Physics
• Formula Sheet

References

• J. C. Slater, “The Theory of Complex Spectra”, Physical Review 34, 1293-1322, 1929, doi:10.1103/PhysRev.34.1293.
• P.-O. Loewdin, “Quantum theory of many-particle systems. I. Physical interpretations by means of density matrices, natural spin-orbitals, and convergence problems in the method of configurational interaction”, Physical Review 97, 1474-1489, 1955, doi:10.1103/PhysRev.97.1474.
• A. Szabo and N. S. Ostlund, Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory, Dover, 1996.
• T. Helgaker, P. Jorgensen, and J. Olsen, Molecular Electronic-Structure Theory, Wiley, 2000.
• R. J. Bartlett and M. Musial, “Coupled-cluster theory in quantum chemistry”, Reviews of Modern Physics 79, 291-352, 2007, doi:10.1103/RevModPhys.79.291.
• J. Rissler, R. M. Noack, and S. R. White, “Measuring orbital interaction using quantum information theory”, Chemical Physics 323, 519-531, 2006, doi:10.1016/j.chemphys.2005.10.018.
• O. Legeza and J. Solyom, “Optimizing the density-matrix renormalization group method using quantum information entropy”, Physical Review B 68, 195116, 2003, doi:10.1103/PhysRevB.68.195116.
• G. K.-L. Chan and S. Sharma, “The density matrix renormalization group in quantum chemistry”, Annual Review of Physical Chemistry 62, 465-481, 2011, doi:10.1146/annurev-physchem-032210-103338.
• K. Boguslawski, P. Tecmer, O. Legeza, and M. Reiher, “Entanglement measures for single- and multireference correlation effects”, Journal of Physical Chemistry Letters 3, 3129-3135, 2012, doi:10.1021/jz301319v.

Exercises

1. Single determinant in its own spin-orbital basis. Consider a determinant occupying spin-orbitals $1$ and $2$ in a four-spin-orbital basis:

$$\lvert\Phi\rangle = \lvert1100\rangle.$$

What is the entropy of spin-orbital mode $1$?

Solution

Mode $1$ has definite occupation $n_1=1$. Its reduced state is the pure one-mode state

$$\rho_1 = \lvert1\rangle\langle1\rvert.$$

Therefore

$$S(\rho_1)=0.$$

The determinant is antisymmetric as a fermionic state, but in this chosen occupation basis each spin-orbital has definite occupation.

2. A two-configuration superposition. Let

$$\lvert\Psi\rangle = \frac{ \lvert1100\rangle + \lvert0011\rangle }{\sqrt2}.$$

Treat modes $1,2$ as subsystem $A$ and modes $3,4$ as subsystem $B$. Compute $S_A$ in bits.

Solution

The state is already in Schmidt form across $A\vert B$:

$$\lvert\Psi\rangle = \frac{ \lvert11\rangle_A\lvert00\rangle_B + \lvert00\rangle_A\lvert11\rangle_B }{\sqrt2}.$$

Thus

$$\rho_A = \frac12\lvert11\rangle\langle11\rvert + \frac12\lvert00\rangle\langle00\rvert.$$

The nonzero eigenvalues are $1/2$ and $1/2$, so

$$S_A = 1 \quad \text{bit}.$$

3. Exchange versus correlation. Explain why a single Slater determinant can include exchange but still miss electron correlation.

Solution

The determinant is antisymmetric, so it exactly includes the exchange constraint required for fermions. However, it restricts the many-electron state to one antisymmetrized product of spin-orbitals. Correlated electronic motion generally requires superpositions of determinants, changed occupation patterns, or more flexible ansatz structure beyond one mean-field determinant.

4. One-orbital entropy. A spatial orbital has reduced density matrix diagonal in the basis $\lvert0\rangle,\lvert\uparrow\rangle,\lvert\downarrow\rangle,\lvert\uparrow\downarrow\rangle$ with probabilities $1/2,0,0,1/2$. Compute its entropy in bits.

Solution

Only two eigenvalues are nonzero, both equal to $1/2$. Therefore

$$s_p = -2 \left( \frac12\log_2\frac12 \right) = 1.$$

5. Why basis matters. Give one reason orbital entanglement is basis dependent.

Solution

Orbitals define the mode subsystems. A unitary rotation among spin-orbitals changes which occupation modes are called local subsystems, so the reduced density matrices of individual orbitals or orbital groups can change. Therefore orbital entropies must always be reported together with the chosen orbital basis.