Skip to content

What Are Strong Correlations?

A quantum material is strongly correlated when interactions and collective constraints reorganize the states relevant to an observable so thoroughly that an independent-particle description is not a controlled starting point at that scale. Typical consequences include suppressed charge fluctuations, local moments, narrow coherence scales, strongly renormalized or absent quasiparticles, large transfers of spectral weight, and low-energy degrees of freedom that are not bare electrons.

This is an operational definition, not a universal inequality. A material may be strongly correlated for low-energy transport yet admit a useful high-energy band description. A heavy Fermi liquid can have sharp quasiparticles and still be strongly correlated because their mass and coherence scale arise collectively. Conversely, a broad spectral line can come from disorder, phonons, instrumental resolution, or an unresolved band manifold rather than strong electronic correlation.

We connect interaction scales to local moments, quasiparticle weight and measurable response. Mott Insulators examines interaction-driven charge incompressibility; Hubbard Physics in Materials explains how a crystal becomes an effective lattice model.

Required background. Interacting Many-Body Systems Overview supplies dimensionless interaction scales and the distinction between weak and strong coupling. Hubbard Model supplies the hopping, local-repulsion, and atomic-limit benchmark. Lifetime and Spectral Weight supplies pole widths and coherent weight for testing quasiparticle validity, while Metals, Insulators, and Semiconductors distinguishes band, interaction, and localization mechanisms.

The phrase “strongly correlated” is used in three related but non-equivalent senses.

Meaning Physical question
Competing energy scales Are interactions large compared with kinetic energy or a coherence scale?
Failure of an approximation Does an independent-particle or weak-coupling calculation reproduce the relevant response?
Emergent degrees of freedom Do local moments, heavy particles or fractional excitations describe the low-energy states?

None is sufficient by itself. A large bare Coulomb energy can be screened into a conventional metal. A weak microscopic interaction can become decisive near a nested Fermi surface, in one dimension, or at a critical point. A local-density approximation can miss a gap for reasons that do not prove a particular Mott mechanism. Strong correlation is therefore a statement about scales, observables, and the failure of specified approximations.

The distinction between strong interaction and strong correlation is especially important. Interaction strength is a Hamiltonian parameter. Correlation is structure in a state or response that cannot be reduced to independent occupations in the chosen description. A symmetry-broken mean-field state may capture a large interaction with one determinant, while a modest interaction in a frustrated or low-dimensional system can generate highly nonclassical correlations.

A downfolded electronic Hamiltonian often has the schematic form

H=∑ij,ab,σtijabciaσ†cjbσ+∑i,aUania↑nia↓+12∑i≠ja,bVijabnianjb+HHund+Hso+He−ph.H ={} \sum_{ij,ab,\sigma} t_{ij}^{ab} c_{ia\sigma}^\dagger c_{jb\sigma} + \sum_{i,a} U_a n_{ia\uparrow}n_{ia\downarrow} + \frac{1}{2} \sum_{\substack{i\ne j\\a,b}} V_{ij}^{ab}n_{ia}n_{jb} + H_{\mathrm{Hund}} + H_{\mathrm{so}} + H_{\mathrm{e-ph}}.

The hopping matrix includes the chosen one-body orbital structure and sets bare bandwidths and velocities. Interorbital local repulsions, not written separately here, are also needed in a general multiorbital model. The local repulsion UU, intersite interaction VV, Hund coupling, spin–orbit coupling, crystal-field splittings, and electron–phonon terms compete with those kinetic scales. Temperature, disorder, pressure, dimensionality, filling, and frustration decide which competition is visible.

For a one-band lattice model,

λU=UW\lambda_U = \frac{U}{W}

is a useful first coordinate. It is not a phase criterion. Both UU and WW depend on the chosen orbital basis and downfolding window. The critical ratio for a metal–insulator transition depends on lattice geometry, density of states, magnetic order, temperature, degeneracy, and whether longer-range interactions are retained.

Other systems require different coordinates:

rs∼Coulomb energyFermi energyr_s \sim \frac{\text{Coulomb energy}}{\text{Fermi energy}}

organizes dilute continuum electrons, while a Hund metal may be controlled by JHJ_H, orbital filling, and a suppressed coherence scale even before UU reaches a single-band Mott threshold (Georges, de’ Medici and Mravlje, 2013). A Kondo material can be strongly correlated because an exponentially small screening scale emerges from moderate microscopic exchange. Flat bands amplify interactions by reducing WW. One-dimensional fermions can lose the Landau quasiparticle pole through infrared collective behavior without a universal large bare coupling.

Independent-particle band theory assumes that the many-electron problem can be organized by occupied one-electron states, perhaps after replacing bare electrons by weakly dressed quasiparticles. It fails in several distinct ways:

  • charge configurations become locally constrained rather than approximately independent;
  • several atomic multiplets or orbital sectors remain dynamically active;
  • the self-energy varies strongly with frequency or momentum;
  • quasiparticle weight becomes small or the pole broadens into a continuum;
  • spectral weight moves across energies comparable to UU, not merely near the Fermi level;
  • doping changes the available low-energy Hilbert space rather than rigidly shifting a chemical potential;
  • competing orders, frustration, or critical fluctuations prevent a single stable reference state.

Failure is observable dependent. A density-functional calculation may predict the lattice and broad orbital character well while failing for a charge gap or magnetic entropy. A renormalized band may fit a narrow low-temperature window while missing the incoherent response above the coherence scale.

For one band, write

GR(k,ω)=1ℏω+μ−εk−ΣR(k,ω).G^R(\mathbf k,\omega) = \frac{1}{ \hbar\omega+\mu -\varepsilon_{\mathbf k} -\Sigma^R(\mathbf k,\omega) }.

The pole energy EkE_{\mathbf k} is measured relative to the chemical potential and solves Ek+μ−εk−Re⁡ΣR(k,Ek/ℏ)=0E_{\mathbf k}+\mu-\varepsilon_{\mathbf k}-\operatorname{Re}\Sigma^R(\mathbf k,E_{\mathbf k}/\hbar)=0. When the imaginary part is small and varies slowly across the pole, the residue is

Zk=[1−∂Re⁡ΣR(k,ω)∂(ℏω)]ω=Ek/ℏ−1.Z_{\mathbf k} = \left[ 1 - \frac{\partial \operatorname{Re}\Sigma^R(\mathbf k,\omega) }{ \partial(\hbar\omega) } \right]_{\omega=E_{\mathbf k}/\hbar}^{-1}.

The pole half-width in energy units is approximately

Γk⋆=−ZkIm⁡ΣR(k,Ek/ℏ).\Gamma_{\mathbf k}^{\star} = - Z_{\mathbf k} \operatorname{Im} \Sigma^R( \mathbf k,E_{\mathbf k}/\hbar ).

A useful quasiparticle requires more than Zk>0Z_{\mathbf k}>0: its width must be small compared with the energy scale on which it is used. Small ZZ means that a bare-electron insertion overlaps only weakly with the coherent excitation; the missing weight resides in incoherent many-body states. It does not imply a broad peak. Nor does ZZ alone determine the effective mass: momentum dependence of the self-energy also changes the band velocity. Heavy-fermion metals provide a central example of small-ZZ, large-mass quasiparticles that become sharply coherent only below a low temperature.

The converse warning is equally important. A broad peak or mass enhancement does not uniquely identify electron–electron correlations. Electron–phonon coupling, disorder, multiband hybridization, surface reconstruction, and finite experimental resolution must be tested.

The half-filled repulsive Hubbard model gives the cleanest local diagnostic. Define the double occupancy

D=⟨ni↑ni↓⟩D = \left\langle n_{i\uparrow}n_{i\downarrow} \right\rangle

and the local spin-polarization fluctuation

mloc2=⟨(ni↑−ni↓)2⟩=⟨ni⟩−2D.\begin{aligned} m_{\mathrm{loc}}^2 &= \left\langle \left( n_{i\uparrow}-n_{i\downarrow} \right)^2 \right\rangle \\ &= \langle n_i\rangle - 2D. \end{aligned}

The quantity mloc2m_{\mathrm{loc}}^2 is dimensionless: it equals 4⟨(Siz)2⟩/ℏ24\langle(S_i^z)^2\rangle/\hbar^2, not ⟨Si2⟩\langle\mathbf S_i^2\rangle. It measures an instantaneous local fluctuation; a persistent local moment also has a characteristic correlation time.

Repulsion suppresses DD. Near one electron per site, this increases the local-moment amplitude even before the moments develop long-range order. At large U/tU/t, charge-changing states cost energy of order UU, while the low-energy sector is approximately restricted to one electron per site.

Virtual hopping still acts inside that constrained sector. In the simplest half-filled one-band limit it produces an antiferromagnetic exchange scale

J∼4t2U.J \sim \frac{4t^2}{U}.

This scale separation is a model example of emergence: high-energy charge fluctuations generate low-energy spin dynamics. The controlled projection and its corrections are derived in Effective Hamiltonians in Many-Body Systems, while Exchange Interactions places the result among direct, superexchange, carrier-mediated, and itinerant mechanisms.

Local moments do not by themselves prove a Mott insulator. Hund metals, dilute magnetic alloys, charge-transfer materials, and bad metals can support substantial fluctuating moments. Conversely, a nonmagnetic singlet can be strongly correlated. The evidence must connect filling, charge response, spectra, and spin dynamics.

The displayed Green function has units of inverse energy. For angular frequency ω\omega, define the frequency spectral density

A(k,ω)=−ℏπIm⁡GR(k,ω).A(\mathbf k,\omega) =-\frac{\hbar}{\pi}\operatorname{Im}G^R(\mathbf k,\omega).

Equivalently, A(k,ω)=ℏAE(k,E)A(\mathbf k,\omega)=\hbar A_E(\mathbf k,E) at E=ℏωE=\hbar\omega, where AE=−Im⁡GR/πA_E=-\operatorname{Im}G^R/\pi is normalized using dEdE. For a normalized single-particle orbital, the exact frequency sum rule is

∫−∞∞dω A(k,ω)=1.\int_{-\infty}^{\infty} d\omega\, A(\mathbf k,\omega) = 1.

Interactions do not remove this weight; they redistribute it. A weakly dressed band concentrates most of the weight near a quasiparticle pole. A correlated metal may show a narrow coherent feature plus broad removal and addition continua. Approaching a Mott regime can suppress low-energy weight and build lower and upper Hubbard sectors. Doping can transfer weight between these sectors over an energy range far larger than the chemical-potential shift.

Two normalized spectra with coherent weights 0.8 and 0.25; reducing the central peak area increases the areas of the broad satellites.

The normalized Lorentzian model in Exercise 4, with η=0.15ω0\eta=0.15\omega_0, Γ=0.45ω0\Gamma=0.45\omega_0 and satellite centers ±2.5ω0\pm2.5\omega_0. Both spectra have total area one over all frequencies. The central component has area ZZ, not height ZZ; its height also depends on its width. These illustrative curves are not spectra calculated for a specific material.

Optical data provide a complementary partial sum,

W(Ωc)=∫0Ωcdω σ1(ω).W(\Omega_c) = \int_0^{\Omega_c} d\omega\, \sigma_1(\omega).

The cutoff Ωc\Omega_c must be stated. A lost Drude weight may reappear in a mid-infrared feature or in higher interband structure. Matrix elements, multiple orbitals and temperature-dependent lattice effects complicate the comparison. A changing partial integral is not a violation of the full optical sum rule. Broad frequency coverage helps distinguish redistribution from an apparent loss outside the measured window (Basov and collaborators, 2011).

This redistribution distinguishes a correlated evolution from a simple rigid-band picture, but it is not unique to a single microscopic model. Coupling to phonons produces satellites; disorder broadens and localizes; density-wave order reconstructs bands. The strongest diagnosis combines one-particle spectra with charge, spin, optical, and thermodynamic information.

An emergent degree of freedom is a variable that efficiently describes the low-energy sector but is not a freely propagating constituent of the microscopic Hamiltonian. Strong correlations often create such variables by freezing some fluctuations and retaining others.

Localized charge can leave spin or orbital moments; a coherent Kondo lattice supports heavy quasiparticles; paired fermions have collective phase fluctuations and Bogoliubov excitations. These are different reductions of the microscopic dynamics. Quantum Spin Liquids and Fractionalization examine more radical possibilities in which excitations belong to distinct deconfined sectors.

Fractionalization is a stronger claim than “the electron is dressed.” An electron-like Landau quasiparticle retains the electron’s conserved quantum-number sector with renormalized overlap. Fractionalized excitations separate those quantum numbers or introduce new topological sectors, and they require evidence for deconfinement rather than a suggestive continuum alone.

Emergence also occurs without exotic fractionalization. Screening converts a bare Coulomb interaction into collective dielectric response. Local moments interact through exchange. A Fermi surface organizes low-temperature thermodynamics even though it is a many-electron object. Anderson’s maxim that “more is different” is useful precisely when it is read as a demand for a correct effective description, not as permission to detach phenomenology from microscopic constraints.

Several complementary measurements are needed because different mechanisms can produce similar anomalies.

Determine composition, valence, crystal structure, dimensionality, carrier density, and orbital content. Compare several electronic-structure approximations where possible. State whether the baseline predicts a metal, band insulator, semimetal, or ordered state and which observable it fails to reproduce.

Kohn–Sham eigenvalues are not exact excitation energies, and a failed semilocal functional is not itself proof of Mott physics. Structural errors, magnetic order, spin–orbit coupling, nonlocal exchange, and an incomplete orbital window can change the baseline.

Photoemission and tunneling constrain occupied or local spectral weight; inverse probes access addition spectra; optical spectroscopy follows current-carrying weight; resonant x-ray and neutron methods probe orbital, charge, and spin sectors. A narrow energy window can mistake a pseudogap, hybridization gap, density-wave gap, or surface feature for a universal correlation gap.

Specific heat, compressibility, magnetic susceptibility, entropy, charge stiffness, and transport test different parts of the effective theory. Large linear specific heat can indicate a heavy mass, but phonons and nuclear terms must be removed. Bad-metal resistivity can signal loss of coherent transport, but saturation, inhomogeneity, and multiband conduction matter. Linear-in-temperature resistivity is not a standalone definition of strong correlation or “Planckian” dynamics.

Observation What else can produce it?
Large effective mass A flat bare band or electron–phonon dressing
Broad or missing spectral peak Disorder, matrix elements or finite resolution
Insulation at partial band filling Symmetry breaking, structural change or localization
Mid-infrared optical weight Interband transitions or polaron absorption
Fluctuating local moments Dilute impurities or mixed valence
A scattering continuum Ordinary multiparticle excitations or decay

The material diagnosis should survive tuning. Pressure often increases hopping and bandwidth; doping changes filling and screening; magnetic field changes spin and orbital scales; strain changes geometry and crystal fields. A model becomes more convincing when the same parameters explain several trends without being independently refitted for every observable.

Example: local moments in a two-site singlet

Section titled “Example: local moments in a two-site singlet”

Two sites already illustrate how charge fluctuations and spin dynamics separate. Take two electrons with hopping t>0t>0 and repulsion U≥0U\geq0:

H=−t∑σ(c1σ†c2σ+c2σ†c1σ)+U∑i=12ni↑ni↓.H=-t\sum_\sigma(c_{1\sigma}^\dagger c_{2\sigma}+c_{2\sigma}^\dagger c_{1\sigma}) +U\sum_{i=1}^2n_{i\uparrow}n_{i\downarrow}.

The singlet ground-state energy is Es=(U−U2+16t2)/2E_s=(U-\sqrt{U^2+16t^2})/2 and the triplet energy is zero. The Hubbard dimer derives the exact spectrum. For our purpose, use the Hellmann–Feynman derivative: ∂UEs\partial_U E_s is the total double occupancy, so the occupancy per site is

D=12∂Es∂U=14(1−UU2+16t2).D=\frac12\frac{\partial E_s}{\partial U} =\frac14\left(1-\frac{U}{\sqrt{U^2+16t^2}}\right).

At U=0U=0, D=1/4D=1/4. At U=8tU=8t, D≃0.0264D\simeq0.0264 and mloc2=1−2D≃0.947m_{\mathrm{loc}}^2=1-2D\simeq0.947. The large local fluctuation coexists with a total-spin singlet and ⟨Siz⟩=0\langle S_i^z\rangle=0: local moments do not require a nonzero ordered moment.

The singlet–triplet gap, with its large-U/tU/t expansion, is

Δst=U2+16t2−U2=4t2U+O(t4/U3).\Delta_{st}=\frac{\sqrt{U^2+16t^2}-U}{2} =\frac{4t^2}{U}+O(t^4/U^3).

At U=8tU=8t, its exact value 0.472t0.472t is close to the leading exchange estimate 0.500t0.500t. Increasing UU suppresses charge fluctuations but also lowers the spin scale. A two-site system does not establish a bulk Mott phase; it isolates the local mechanism before lattice geometry and collective order enter.

For NN spinful orbitals, the full Fock space has dimension

dim⁡H=4N.\dim\mathcal H = 4^N.

The kinetic term is simplest in momentum space, local interactions in real space, and neither basis diagonalizes the full problem. Fermionic signs, frustration, long-range entanglement, low coherence scales, and competing nearly degenerate phases make both analysis and numerics difficult. Real-frequency experiments add an inverse problem because many finite-temperature methods work in imaginary time.

Different methods control different approximations. Exact diagonalization controls a finite cluster; Monte Carlo can control statistical and size errors in favorable sign regimes; tensor networks control entanglement truncation; dynamical mean-field theory treats local temporal fluctuations while approximating spatial correlations. Computational Quantum Matter explains how to choose and compare them. LeBlanc and collaborators (2015) show why independent methods and size extrapolation remain valuable even for the simple square-lattice Hubbard model.

At half filling, evaluate mloc2=⟨ni⟩−2Dm_{\mathrm{loc}}^2=\langle n_i\rangle-2D for an uncorrelated paramagnet with D=1/4D=1/4 and for an ideal singly occupied limit with D=0D=0. Interpret the change.

Solution

At half filling, ⟨ni⟩=1\langle n_i\rangle=1. For independent up and down occupations,

mloc2=1−2(14)=12.m_{\mathrm{loc}}^2 = 1 - 2\left( \frac{1}{4} \right) = \frac{1}{2}.

In the singly occupied limit,

mloc2=1.m_{\mathrm{loc}}^2 = 1.

Repulsion suppresses empty–double fluctuations and makes a spinful singly occupied state more probable. This establishes local-moment amplitude, not long-range magnetic order.

Suppose near the Fermi level

ΣR(ω)=Σ0+(1−Z0−1)ℏω−iγ,\Sigma^R(\omega) = \Sigma_0 + \left( 1-Z_0^{-1} \right) \hbar\omega - i\gamma,

with real Σ0\Sigma_0, 0<Z0≤10<Z_0\le1, and γ>0\gamma>0. Find the quasiparticle residue and the renormalized width.

Solution

The derivative of the real self-energy is

∂Re⁡ΣR∂(ℏω)=1−Z0−1,\frac{\partial \operatorname{Re}\Sigma^R }{ \partial(\hbar\omega) } = 1-Z_0^{-1},

so

Z=[1−(1−Z0−1)]−1=Z0.Z = \left[ 1- \left( 1-Z_0^{-1} \right) \right]^{-1} = Z_0.

The pole denominator can be multiplied by Z0Z_0, showing that the renormalized half-width is

Γ⋆=Z0γ.\Gamma^\star = Z_0\gamma.

A small residue does not automatically make the pole broad. For a low-temperature excitation away from the Fermi surface, compare Γ⋆\Gamma^\star with its excitation energy; near the Fermi surface, temperature and experimental resolution also determine what can be resolved. A finite constant γ\gamma here describes a local spectral window, not the asymptotic clean Fermi-liquid scattering law.

A one-band fit gives U=4.0 eVU=4.0\,\mathrm{eV} and W=2.0 eVW=2.0\,\mathrm{eV} at ambient pressure. Under pressure, assume U=3.8 eVU=3.8\,\mathrm{eV} and W=3.0 eVW=3.0\,\mathrm{eV}. Compute U/WU/W and explain what can and cannot be concluded.

Solution

The ratios are

UW∣P=0=2.0,UW∣P≈1.27.\left. \frac{U}{W} \right|_{P=0} = 2.0, \qquad \left. \frac{U}{W} \right|_{P} \approx 1.27.

Pressure moves this model toward greater itinerancy. The numbers do not by themselves locate a phase boundary. One still needs filling, orbital degeneracy, magnetic and structural order, temperature, longer-range interactions, and the downfolding convention. Spectral, transport, and thermodynamic evolution must test the proposed trend.

Model a momentum-resolved spectrum as

A(ω)=ZLη(ω)+1−Z2[LΓ(ω−Ω)+LΓ(ω+Ω)],A(\omega) ={} ZL_\eta(\omega) + \frac{1-Z}{2} \left[ L_\Gamma(\omega-\Omega) + L_\Gamma(\omega+\Omega) \right],

where La(x)=a/[π(x2+a2)]L_a(x)=a/[\pi(x^2+a^2)] integrates to one. Verify the sum rule and identify the coherent and incoherent weights.

Solution

Integrating term by term gives

∫dω A(ω)=Z+1−Z2(1+1)=1.\int d\omega\,A(\omega) = Z + \frac{1-Z}{2} \left( 1+1 \right) = 1.

The central coherent feature carries weight ZZ. The two incoherent satellites together carry 1−Z1-Z. Correlations redistribute weight without violating the single-orbital sum rule. In a real multiorbital experiment, matrix elements and a finite energy window complicate this direct partition.

A half-filled material becomes antiferromagnetic and insulating at the same temperature. List measurements that would help distinguish a correlation-driven local-moment insulator from a weak-coupling Slater insulator.

Solution

Useful tests include whether a charge gap or pseudogap persists above the ordering temperature, whether local moments and substantial magnetic entropy exist above it, whether optical weight moves over an interaction-scale window, whether the gap follows the magnetic order parameter, and whether realistic ordered band theory reproduces both charge and spin scales. Pressure or frustration that suppresses order is particularly informative. The two mechanisms can interpolate or coexist, so the goal is a quantitative hierarchy of scales rather than a forced binary label.