Bogoliubov Excitations in BECs

Updated 17 January 2026

Bogoliubov excitations are fundamental quantum fluctuations emerging from linearizing the mean-field solution of an interacting Bose–Einstein condensate.
They are obtained by a quadratic expansion of the Hamiltonian leading to coupled particle-hole modes diagonalized via paraunitary transformations, revealing key dispersion and stability features.
Their study is crucial for understanding superfluidity, topological edge modes, and decay processes in equilibrium and non-equilibrium quantum systems.

Bogoliubov excitations are the fundamental quantum fluctuations (quasiparticles) arising on top of a Bose–Einstein condensate (BEC) or, more generally, the mean-field ground state of an interacting bosonic system. They are obtained by linearizing the dynamics around the condensate solution, yielding an effective quadratic Hamiltonian whose elementary excitations mix particle and hole amplitudes. The properties, stability, topological character, and decay dynamics of these modes underpin the understanding of condensed matter, cold atom, and quantum optical systems, and are central to phenomena such as superfluidity, Anderson localization, and topological bosonic phases.

1. Mathematical Structure of Bogoliubov Excitations

Starting from a weakly interacting Bose gas with a macroscopically occupied condensate, the field operator is split as $\hat{a}_{m,s}\to\sqrt{N_0/M}\,\zeta_{m,s}+\hat{a}_{m,s}$ . Expanding the Hamiltonian to quadratic order in the fluctuation operators gives, in momentum space, a Bogoliubov Hamiltonian of the form

$H^{(B)} = \frac12\sum_k (\hat{a}_k^\dagger,\, \hat{a}_{-k})\,H_{\rm Bog}(k)\,(\hat{a}_k,\,\hat{a}_{-k}^\dagger)^T,$

where $H_{\rm Bog}(k)$ is a $2\mathcal{N}\times2\mathcal{N}$ Hermitian matrix, with $\mathcal{N}$ the number of internal degrees of freedom (spin, sublattice, component).

Because $H_{\rm Bog}(k)$ couples particle and hole sectors, it is diagonalized by a paraunitary (symplectic) transformation $T(k)$ satisfying

$T(k)^\dagger \sigma_z T(k) = \sigma_z, \quad \sigma_z = \mathrm{diag}(+1_\mathcal{N},-1_\mathcal{N}),$

yielding transformed quasiparticle operators $(b_k,b_{-k}^\dagger)^T$ , and a diagonalized spectrum comprising positive and negative eigenvalues $\{E_j(k),-E_j(-k)\}$ . The associated eigenvectors have definite bi-orthogonality in the $\sigma_z$ -weighted inner product (Engelhardt et al., 2015).

In spatially inhomogeneous systems, notably disordered or trapped BECs, the diagonalization generalizes to the Bogoliubov–de Gennes (BdG) system, where eigenmodes $(u,v)$ satisfy

$H_{\rm BdG}\begin{pmatrix} u \ v \end{pmatrix} = E \sigma_z \begin{pmatrix} u \ v \end{pmatrix},$

with $H_{\rm BdG}$ a generalized Hermitian matrix (or operator).

2. Dispersion and Physical Regimes

Homogeneous Bose Gas

For a translation-invariant BEC with $g>0$ , the celebrated Bogoliubov dispersion reads

$\omega_k = \sqrt{\epsilon_k(\epsilon_k+2gn)}\,,\quad \epsilon_k = \frac{\hbar^2k^2}{2m},$

with amplitudes $u_k^2 = \frac{\epsilon_k+gn}{2\omega_k}+1/2$ , $v_k^2 = \frac{\epsilon_k+gn}{2\omega_k}-1/2$ (Frérot et al., 2023). At low $k$ , $\omega_k\sim c_s k$ with sound velocity $c_s=\sqrt{gn/m}$ ; for $k\xi\gg1$ , the spectrum is particle-like.

Inhomogeneous and Multi-Component Systems

In systems with trapping, disorder, or multiple components, the excitation energies are given by solving the BdG eigenproblem for the relevant potential or interaction structure (Xie et al., 15 Jun 2025, Gaul et al., 2011). The bi-orthogonality and spectrum symmetry $\omega \leftrightarrow -\omega^*$ persist but require careful handling of zero modes (Goldstone, dipole, etc.) (Xie et al., 15 Jun 2025).

In driven-dissipative or non-equilibrium systems (e.g., polariton condensates), the spectrum can become complex, and the modes may acquire finite lifetimes due to coupling to environment or reservoir degrees of freedom (Pieczarka et al., 2021, Frérot et al., 2023).

3. Topological and Symmetry-Protected Excitations

Symplectic Polarization and Inversion Symmetry

For inversion-symmetric systems, the Bogoliubov bands inherit a nontrivial topology characterized by a quantized symplectic polarization $P_s$ . The Berry connection for the positive-energy bands is

$A(k) = i \sum_{\lambda\leq\lambda_{\rm max}} \mathrm{Tr} [\Gamma_\lambda \sigma_z t_k^\dagger \sigma_z \partial_k t_k],$

with the integrated polarization

$P_s = \frac{1}{2\pi} \int_{\rm BZ} dk\,A(k),$

which quantizes to $P_s=0,\,\frac12$ (mod 1) in the presence of inversion symmetry. This quantization can be expressed by the eigenvalues of the inversion operator at time-reversal invariant momenta: $P_s = \frac{i}{2\pi} \ln \prod_{\lambda\leq\lambda_{\rm max}} \eta_\lambda(0)\eta_\lambda(\pi),$ where $\eta_\lambda=±1$ are inversion eigenvalues (Engelhardt et al., 2015).

Bulk–Boundary Correspondence

A nontrivial $P_s$ guarantees edge-localized Bogoliubov modes in systems with open boundaries. The presence of a gap and topological invariant ensures midgap collective excitations, robust to disorder preserving the symmetry. Experimentally, such edge states can be probed via spectroscopic measurements such as Bragg scattering.

Higher-Order Topological Bogoliubov Modes

Bosonic topological excitations extend beyond conventional edge states. For instance, a trivial $s$ -wave superfluid, when subjected to a mirror-symmetric onsite potential, can realize gapped excitation bands supporting symmetry-protected corner modes localized at zero-dimensional regions. The protection is afforded by mirror or sublattice symmetry, enforcing a nontrivial topological index in excitation (not ground state) bands (Tu et al., 2023, Guo et al., 2024).

Squeezing-Induced Topological Gaps

Pairing (squeezing) terms in the BdG Hamiltonian can induce effective spin-orbit or Zeeman-type couplings, opening topological gaps in bosonic excitation spectra. The associated topological invariants (Chern numbers, $\mathbb{Z}_2$ indices) are directly computable from the block-diagonalized particle (number-conserving) part of the transformed Hamiltonian, with the squeezed representation maintaining the full spectrum and topology (Wan et al., 2020).

4. Decay, Interactions, and Non-Equilibrium Effects

Decay Processes in One-Dimensional Gases

In 1D, the leading decay channel for Bogoliubov excitations is via disintegration into three lower-energy quasiparticles. The decay rate for phonons scales as $Q^7$ , is negligible at low energy, but becomes prominent with integrability-breaking three-body or finite-range two-body interactions. True integrable models (Lieb–Liniger with delta interactions) exhibit exact absence of decay for all momenta, underscoring the impact of quantum integrability (Ristivojevic et al., 2016).

Driven and Dissipative Systems

In driven-dissipative condensates (photon or polariton fluids), coupling to reservoirs (thermal phonons, vacuum fluctuations) alters both the excitation spectrum and their generation rates. For instance, in microcavity polariton systems, the steady-state occupation of Bogoliubov modes can be dominated by thermal lattice phonons well above quantum-vacuum contributions, with implications for the observation of quantum correlations and analog Hawking radiation (Frérot et al., 2023).

In paraxial fluids of light, interference between counter-propagating Bogoliubov excitations is directly observed, with the linear, non-saturating phase shift providing a "smoking gun" for superfluidity and collective behavior (Fontaine et al., 2020).

5. Numerical and Variational Approaches

Bi-Orthogonal Spectral Methods

The BdG system is non-Hermitian but possesses a bi-orthogonal structure in the Krein inner product, enabling the development of efficient, spectrally accurate numerical solvers leveraging the fast Fourier transform and structure-preserving orthogonalization. Such approaches allow high-precision calculation of the full excitation spectrum and modes in multi-component and high-dimensional trapped condensates (Xie et al., 15 Jun 2025).

Time-Dependent Variational Principle

Alternatively, variational approaches based on time-dependent Gaussian (or extended Gaussian) ansätze map the Gross–Pitaevskii equation to a dynamical system for the width and deformation parameters. The eigenvalues of the linearized Jacobian at the fixed point coincide with the Bogoliubov excitation frequencies, reproducing all low-lying modes, including angular rotons in dipolar condensates, and providing a computationally tractable route for high-dimensional or complex systems (Kreibich et al., 2012, Kreibich et al., 2012, Kreibich et al., 2012).

Quantum-Defect and Rydberg Analysis of Excitation Spectra

In self-trapped condensates with long-range (e.g., monopolar $1/r$) interactions, the Bogoliubov spectrum exhibits a Rydberg-like structure with quantum defects encoding deviation from pure Coulombic asymptotics. This provides a physical and computational framework for understanding high-lying excitations and their deviation from simple models (Kreibich et al., 2012).

6. Disorder and Quasiperiodic Dynamics

Bogoliubov excitations in disordered or quasiperiodic systems display rich localization and delocalization physics. In the presence of spatially correlated potential disorder, the excitation spectrum and sound velocity acquire non-universal shifts, and low-energy modes can localize with disorder-dependent mean-free paths and localization lengths. In quasiperiodically driven systems (e.g., the kicked rotor), Bogoliubov quasiparticles undergo an Anderson-like localization–diffusion transition matching the universality class (critical exponent $\nu \approx 1.6$ ) of non-interacting systems for weak interactions, even as the interacting condensate remains dynamically stable over experimental timescales (Gaul et al., 2011, Vermersch et al., 2014).

7. Experimental Signatures and Physical Implications

Table: Bogoliubov Excitations—Key Physical Regimes and Features

System/Regime	Distinctive Bogoliubov Properties	References
Homogeneous BEC	Linear sound mode, $\omega_k \sim c_s k$	(Frérot et al., 2023)
Inversion-symmetric	Quantized symplectic polarization, edge modes	(Engelhardt et al., 2015)
Driven-dissipative	Complex energies, thermal/vacuum excitation rates	(Pieczarka et al., 2021, Frérot et al., 2023)
Topological (AIII)	Winding number in 1D, symmetry-protected edge modes	(Guo et al., 2024)
Weyl/Dirac systems	Nodal excitations, bosonic arcs, surface states	(Wu et al., 2016)
Disordered/Quasiperiodic	Anderson localization, identical critical exponents	(Gaul et al., 2011, Vermersch et al., 2014)
Dipolar/Long-range	Angular rotons, Rydberg-like spectra, quantum defects	(Kreibich et al., 2012, Kreibich et al., 2012)

Physical detection of Bogoliubov excitations employs Bragg spectroscopy, time-of-flight measurements, imaging of density-density correlations, and the observation of midgap or edge-localized modes. In complex platforms—e.g., quantum fluids of light, atomic Bose gases, or magnonic systems—Bogoliubov excitations provide essential diagnostics for superfluidity, topological properties, and the character of quantum fluctuations across regimes.

References:

(Engelhardt et al., 2015, Frérot et al., 2023, Gaul et al., 2011, Guo et al., 2024, Wu et al., 2016, Wan et al., 2020, Pieczarka et al., 2021, Xie et al., 15 Jun 2025, Kreibich et al., 2012, Kreibich et al., 2012, Vermersch et al., 2014, Fontaine et al., 2020, Ristivojevic et al., 2016)