| Bogoliubov–de Gennes equations | |
|---|---|
| Name | Bogoliubov–de Gennes equations |
| Caption | Schematic representation of quasiparticle spectra in a superconductor |
| Field | Quantum physics |
| Introduced | 1950s |
| Derived from | BCS theory, Bogoliubov transformation |
| Applications | Superconductivity, Superfluidity, Topological superconductivity |
Bogoliubov–de Gennes equations
The Bogoliubov–de Gennes equations are a set of coupled mean-field eigenvalue equations that describe quasiparticle excitations in inhomogeneous superconductors and superfluids. They generalize the BCS description to spatially varying order parameters and are central for describing phenomena such as vortex cores, Andreev bound states, and edge modes. Their solutions provide access to local density of states, pairing correlations, and topological properties relevant to experiments in condensed matter and cold atoms.
The Bogoliubov–de Gennes (BdG) framework arises in the study of fermionic pairing where a spatially varying pairing potential couples particle and hole degrees of freedom. It combines concepts from the Bogoliubov transformation for superconductors and single-particle quantum mechanics to form a matrix Hamiltonian acting on Nambu spinors. BdG equations are widely used in contexts ranging from conventional BCS theory superconductors to unconventional pairing in heavy-fermion compounds, cuprates, and iron-based superconductors, and in engineered systems such as semiconductor–superconductor heterostructures studied for Majorana physics.
In second-quantized notation the BdG equations follow from a mean-field Hamiltonian expressed in terms of Nambu spinors Ψ = (ψ↑, ψ†↓)ᵀ. The BdG Hamiltonian is a 2×2 (or 4×4 with spin) matrix H_BdG = ĥ(r) - μ, Δ(r)],[Δ*(r), -ĥ*(r) + μ, where ĥ is the single-particle operator (kinetic energy plus potential), μ the chemical potential, and Δ(r) the pairing potential (order parameter). Solving H_BdG u_n = E_n u_n yields quasiparticle amplitudes and energies E_n. Observables such as the local density of states and pair amplitude are constructed from the eigenfunctions through sums over quasiparticle occupations determined by the Fermi–Dirac statistics.
The BdG equations are obtained by applying a mean-field decoupling of the two-body interaction in the BCS Hamiltonian and performing a Bogoliubov transformation that diagonalizes the quadratic Hamiltonian. Starting from a microscopic model (e.g., attractive Hubbard model or electron–phonon coupling models used by Bardeen, Cooper and Schrieffer), a self-consistency condition ties Δ(r) to the anomalous expectation value ⟨ψ↓(r)ψ↑(r)⟩. The derivation invokes concepts from Green's functions, Gor'kov equations, and mean-field theory, and requires careful treatment of ultraviolet regularization in continuum models, often related to renormalization of the pairing interaction as in cold-atom contexts.
Homogeneous solutions recover the BCS quasiparticle dispersion with an energy gap Δ. Inhomogeneous solutions reveal spatially localized states: vortex core spectra in type-II superconductors harbor low-energy Caroli–de Gennes–Matricon states; surfaces or interfaces host Andreev bound states and gapless surface bands in unconventional pairing symmetries (e.g., d-wave). Impurity-induced Yu–Shiba–Rusinov states arise when magnetic impurities break time-reversal symmetry. BdG calculations also predict discrete subgap spectra in quantum dots proximitized by superconductors and bound Majorana zero modes at defects or ends of one-dimensional systems under suitable conditions.
BdG methods are applied to model tunneling spectroscopy (e.g., scanning tunneling microscopy) via local density of states, Josephson junctions and proximity effects in hybrid devices, and the structure of vortex lattices in materials studied at institutions such as Bell Labs, IBM Research, and university research groups at Stanford University and Harvard University. In ultracold gases, BdG-like approaches describe fermionic pairing across the BEC–BCS crossover in experiments at MIT, JILA, and the École Normale Supérieure. They are also foundational for proposals to realize topological superconductivity in devices combining spin–orbit coupling, Zeeman fields, and s-wave pairing, as in semiconductor nanowire platforms explored by groups at Microsoft Station Q and multiple university laboratories.
Numerical solution of BdG equations employs strategies such as diagonalization on finite lattices, kernel polynomial methods, real-space self-consistent iteration, and Bogoliubov–de Gennes tight-binding implementations. Efficient techniques leverage sparse linear algebra, Chebyshev expansion, and parallel computing on clusters or GPUs; software packages and codes are developed within condensed matter groups and national labs like Los Alamos National Laboratory. Boundary conditions, discretization of continuum models, and treatment of large systems for vortex lattices or disordered samples require careful convergence checks and regularization, especially when computing temperature-dependent properties.
The BdG formalism naturally exhibits particle–hole symmetry, placing superconducting systems into symmetry classes of the tenfold way classification; this enables characterization of topological superconductors via invariants (Chern numbers, Z2 indices). Solutions of BdG models predict localized zero-energy Majorana modes at edges, vortices, or domain walls in systems such as the Kitaev chain and proximitized nanowires with strong spin–orbit coupling and Zeeman splitting. These theoretical predictions have motivated experimental searches for Majorana bound states in hybrid devices involving InSb and InAs nanowires coupled to conventional superconductors and are central to proposals for topological quantum computation pursued by academic and industrial groups.