| Bogoliubov–de Gennes equation | |
|---|---|
| Name | Bogoliubov–de Gennes equation |
| Field | Quantum mechanics; Condensed matter physics |
| Introduced | 1950s |
| Authors | Nikolay Bogoliubov; Pierre-Gilles de Gennes |
Bogoliubov–de Gennes equation
The Bogoliubov–de Gennes equation is a set of coupled linear differential equations describing fermionic quasiparticle excitations in spatially inhomogeneous superconducting or superfluidity systems. It generalizes the BCS theory mean-field description to position-dependent order parameters and is central to understanding bound states, vortex cores, and proximity effects in mesoscopic superconductors.
The Bogoliubov–de Gennes (BdG) formalism arises from applying the Bogoliubov transformation to a fermionic mean-field Hamiltonian with pairing terms, producing coupled equations for particle-like and hole-like amplitudes. It is foundational for analyzing phenomena in type-II superconductors, unconventional superconductivity, and fermionic ultracold atoms with pairing. The BdG approach links microscopic models such as the Hubbard model and the Anderson model to experimentally observable quasiparticle properties measured in scanning tunneling microscopy and spectroscopy experiments performed at institutions like Bell Labs and CERN-affiliated condensed matter groups.
In the BdG framework one writes the mean-field Hamiltonian in Nambu spinor form and derives the eigenproblem H_BdG Ψ_n = E_n Ψ_n. For a single-band, spinful superconductor the BdG Hamiltonian typically includes the kinetic term (with Bloch theorem if periodic), a pairing potential Δ(r) and external potentials V(r). The equations couple u_n(r) and v_n(r), the particle and hole components, through Δ(r) and its complex conjugate. When spin-orbit coupling or Zeeman fields are present, the matrix extends to include spin indices, linking the BdG formalism to models studied by groups at MIT and Harvard University investigating topological superconductors. Boundary conditions and self-consistency require that Δ(r) = g Σ_n u_n(r) v*_n(r) f(E_n), where g is the interaction strength and f is the Fermi–Dirac distribution.
Solutions of the BdG equations produce a spectrum of quasiparticle energies E_n, symmetric about zero due to particle-hole redundancy. In homogeneous limits one recovers the BCS dispersion with a gap at the Fermi surface; in inhomogeneous or finite systems discrete bound states appear, such as Caroli–de Gennes–Matricon states in vortex cores. The existence and energy of Andreev bound states at interfaces follows from solving BdG with spatially varying Δ and scattering potentials; these states underpin conductance features in Josephson junctions and Andreev reflection. The spectrum is used to compute observable quantities like local density of states and thermal conductance.
The BdG Hamiltonian possesses a built-in particle-hole symmetry (an antiunitary operator C with C H_BdG C^{-1} = −H_BdG), placing it in the tenfold way classification of topological phases. Additional symmetries—time-reversal, chiral, and crystalline symmetries—determine the topological class and the possibility of protected boundary modes such as Majorana bound states. Studies connecting BdG models to Kitaev chain and p-wave superconductivity have motivated searches for topological superconductivity at interfaces of semiconductor nanowires with strong spin–orbit coupling and proximitized by conventional superconductors (e.g., experiments at Microsoft Station Q collaborations and academic labs).
BdG methods are applied across contexts: conventional s-wave superconductors, unconventional d-wave cuprates, multiband iron-based superconductors, and superfluid phases of Helium-3. They underpin theoretical predictions for Josephson effects, vortex lattice structures, core-level spectroscopy, and proximity-induced pairing in heterostructures combining superconductors and ferromagnets. The formalism is also central to modeling engineered platforms for quantum computation that rely on Majorana modes, as proposed in works by Alexei Kitaev, Roman Lutchyn, and Yuval Oreg.
Solving BdG equations requires numerical linear algebra and self-consistent iteration. Techniques include finite-difference and finite-element discretizations, plane-wave and tight-binding bases, and diagonalization via Lanczos or sparse eigensolvers. For large-scale systems, methods such as the kernel polynomial method, Chebyshev expansion, and recursive Green's function approaches are employed; these are implemented in codes developed at institutions like Argonne National Laboratory and community packages inspired by Kwant and Wannier90. Parallel computing and GPU acceleration are often necessary to treat realistic three-dimensional heterostructures and disorder.
Experimental signatures of BdG-predicted quasiparticles include conductance peaks from Andreev bound states, zero-bias peaks associated with candidate Majorana modes in nanowires and atomic chains, and spatially resolved density-of-states maps of vortex cores imaged by scanning tunneling microscopy groups at Stanford University and IBM Research. Cold-atom experiments in optical lattices provide tunable realizations of BdG physics with controlled interactions and population imbalance, enabling tests of FFLO-like states. Measurements of thermal and electrical transport, tunneling spectroscopy, and Josephson critical currents serve to validate BdG-based modeling across condensed matter and atomic physics platforms.