Bogoliubov-Huang-Meng Approximation

• Bogoliubov-Huang-Meng approximation is a weak-coupling, weak-disorder framework for Bose-condensed systems that combines a quadratic Bogoliubov expansion with a perturbative disorder treatment.
• It employs a canonical Bogoliubov transformation to diagonalize the interacting Hamiltonian and integrates a Huang-Meng shift to capture disorder-induced condensate deformation and depletion.
• The approximation is valid under small depletion and weak disorder, making it a key perturbative tool for studying equilibrium properties and dynamics in disordered Bose gases.

The Bogoliubov-Huang-Meng approximation is a weak-coupling, weak-disorder framework for interacting Bose-condensed systems in random media. In the formulation used for disordered quadrupolar condensates, it assumes a macroscopically occupied uniform condensate, expands the Hamiltonian to quadratic order in fluctuations, diagonalizes the clean interacting part by a Bogoliubov transformation, and incorporates weak static disorder through a disorder-dependent shift of the quasiparticle operators (Boudjemaa, 5 Aug 2025). In the broader dirty-boson literature, the Huang-Meng component denotes the perturbative treatment of disorder fluctuations on top of the Bogoliubov quasiparticle background; it is therefore a disorder extension of standard Bogoliubov theory rather than a synonym for clean-gas Lee-Huang-Yang or Hartree-Fock-Bogoliubov approximations (Nagler et al., 2019).

1. Definition, scope, and historical placement

In the usage documented for disordered Bose gases, the approximation is introduced in the regime of weak interactions, weak disorder, nearly complete condensation, and small fluctuations around a uniform condensate. Its structural content is the combination of a Bogoliubov quadratic truncation with a Huang-Meng treatment of disorder as a weak source that deforms the condensate and generates an additional disorder-induced depletion channel (Boudjemaa, 5 Aug 2025).

This scope matters because several neighboring literatures use similar language while addressing different objects. A rigorous paper on interacting bosons in a homogeneous random medium proves thermodynamic exactness of a generalized Bogoliubov $c$-number substitution for the pressure, but explicitly does not derive the usual Huang-Meng weak-disorder depletion formulas or the disorder-corrected excitation spectrum; its focus is the exactness of the substitution step itself, especially when generalized Bose-Einstein condensation is spread over an infinitesimal low-energy band (Jaeck et al., 2010). Conversely, some dilute-gas papers use “Huang” in the Lee-Huang-Yang sense and contain no disorder at all, so they should not be conflated with Huang-Meng dirty-boson theory (Brietzke et al., 2019).

Within this landscape, the Bogoliubov-Huang-Meng approximation is best understood as a perturbative dirty-boson scheme: it keeps the standard condensate-plus-fluctuation architecture of Bogoliubov theory, but supplements it by a weak-disorder sector that produces condensate deformation and disorder depletion. The approximation is therefore neither a purely spectral Bogoliubov theory nor a general theory of strong disorder, localization, or Bose-glass physics.

2. Quadratic construction around the condensate

The standard starting point is a second-quantized Bose Hamiltonian containing kinetic energy, a random external potential, and two-body interactions. In the quadrupolar extension, the Hamiltonian is written as

$$\hat H =\sum_{\bf k} E_k\hat a^\dagger_{\bf k}\hat a_{\bf k} +\frac{1}{ {\cal V}}\sum_{\bf k,\bf p} U ({\bf k-\bf p}) \hat a^\dagger_{\bf k} \hat a_{\bf p} +\frac{1}{2 {\cal V}}\sum_{\bf k,\bf p,\bf q} V ({\bf q})\, \hat a^\dagger_{\bf k+\bf q} \hat a^\dagger_{\bf p-\bf q}\hat a_{\bf p}\hat a_{\bf k},$$

with $E_k=\hbar^2k^2/(2m)$, zero disorder average $\langle U(\mathbf r)\rangle=0$, and disorder correlator $$\langle U(\mathbf r)U(\mathbf r')\rangle=R(\mathbf r-\mathbf r')$$ (Boudjemaa, 5 Aug 2025).

The Bogoliubov step is the assumption that almost all particles occupy the zero-momentum state, so that

$\hat a_0=\hat a_0^\dagger=\sqrt{N}.$

Keeping only terms quadratic in the nonzero-momentum operators gives

$$\hat H = \frac{1} {2} N n V (0)+n U_0 +\frac{\sqrt{ N} }{ {\cal V} }\sum_{\mathbf k \neq 0} U_{\bf k}\, \left( \hat a^\dagger_{\bf k} +\hat a_{-\bf k}\right) + \sum_{\mathbf k \neq 0} \left [ E_k+ V ({\bf k}) n \right]\hat a^\dagger_{\bf k}\hat a_{\bf k} + \frac{1}{2} n\sum_{\mathbf k \neq 0} V ({\bf k}) \left(\hat a^\dagger_{\bf k} \hat a^\dagger_{-\bf k}+ \hat a_{\bf k} \hat a_{-\bf k} \right).$$

The terms $$\hat a^\dagger_{\bf k}\hat a^\dagger_{-\bf k}+\hat a_{\bf k}\hat a_{-\bf k}$$ are the anomalous pairing sector, while the linear source $$\propto U_{\bf k}(\hat a^\dagger_{\bf k}+\hat a_{-\bf k})$$ is the Huang-Meng disorder contribution (Boudjemaa, 5 Aug 2025).

Diagonalization proceeds by a canonical transformation of the form

$$\hat a_{k}= u_{k} \hat b_{k}- v_{k} \hat b_{-k}^\dagger-\beta_{\bf k},$$

where $$\hat H =\sum_{\bf k} E_k\hat a^\dagger_{\bf k}\hat a_{\bf k} +\frac{1}{ {\cal V}}\sum_{\bf k,\bf p} U ({\bf k-\bf p}) \hat a^\dagger_{\bf k} \hat a_{\bf p} +\frac{1}{2 {\cal V}}\sum_{\bf k,\bf p,\bf q} V ({\bf q})\, \hat a^\dagger_{\bf k+\bf q} \hat a^\dagger_{\bf p-\bf q}\hat a_{\bf p}\hat a_{\bf k},$$0 are quasiparticle operators and $$\hat H =\sum_{\bf k} E_k\hat a^\dagger_{\bf k}\hat a_{\bf k} +\frac{1}{ {\cal V}}\sum_{\bf k,\bf p} U ({\bf k-\bf p}) \hat a^\dagger_{\bf k} \hat a_{\bf p} +\frac{1}{2 {\cal V}}\sum_{\bf k,\bf p,\bf q} V ({\bf q})\, \hat a^\dagger_{\bf k+\bf q} \hat a^\dagger_{\bf p-\bf q}\hat a_{\bf p}\hat a_{\bf k},$$1 is a disorder-induced shift. After this step,

$$\hat H =\sum_{\bf k} E_k\hat a^\dagger_{\bf k}\hat a_{\bf k} +\frac{1}{ {\cal V}}\sum_{\bf k,\bf p} U ({\bf k-\bf p}) \hat a^\dagger_{\bf k} \hat a_{\bf p} +\frac{1}{2 {\cal V}}\sum_{\bf k,\bf p,\bf q} V ({\bf q})\, \hat a^\dagger_{\bf k+\bf q} \hat a^\dagger_{\bf p-\bf q}\hat a_{\bf p}\hat a_{\bf k},$$2

In this sense, the Bogoliubov-Huang-Meng approximation is literally a Bogoliubov diagonalization of the interacting condensate plus a Huang-Meng disorder translation (Boudjemaa, 5 Aug 2025).

3. Disorder-induced deformation, depletion, and equilibrium observables

A defining equilibrium observable is the disorder-induced contribution to the noncondensed density,

$$\hat H =\sum_{\bf k} E_k\hat a^\dagger_{\bf k}\hat a_{\bf k} +\frac{1}{ {\cal V}}\sum_{\bf k,\bf p} U ({\bf k-\bf p}) \hat a^\dagger_{\bf k} \hat a_{\bf p} +\frac{1}{2 {\cal V}}\sum_{\bf k,\bf p,\bf q} V ({\bf q})\, \hat a^\dagger_{\bf k+\bf q} \hat a^\dagger_{\bf p-\bf q}\hat a_{\bf p}\hat a_{\bf k},$$3

This quantity is identified as the “density of the disorder-averaged condensate” or “condensate deformation,” and it is distinct from interaction-induced quantum depletion (Boudjemaa, 5 Aug 2025). The total depletion is decomposed as

$$\hat H =\sum_{\bf k} E_k\hat a^\dagger_{\bf k}\hat a_{\bf k} +\frac{1}{ {\cal V}}\sum_{\bf k,\bf p} U ({\bf k-\bf p}) \hat a^\dagger_{\bf k} \hat a_{\bf p} +\frac{1}{2 {\cal V}}\sum_{\bf k,\bf p,\bf q} V ({\bf q})\, \hat a^\dagger_{\bf k+\bf q} \hat a^\dagger_{\bf p-\bf q}\hat a_{\bf p}\hat a_{\bf k},$$4

with $$\hat H =\sum_{\bf k} E_k\hat a^\dagger_{\bf k}\hat a_{\bf k} +\frac{1}{ {\cal V}}\sum_{\bf k,\bf p} U ({\bf k-\bf p}) \hat a^\dagger_{\bf k} \hat a_{\bf p} +\frac{1}{2 {\cal V}}\sum_{\bf k,\bf p,\bf q} V ({\bf q})\, \hat a^\dagger_{\bf k+\bf q} \hat a^\dagger_{\bf p-\bf q}\hat a_{\bf p}\hat a_{\bf k},$$5 the zero-temperature interaction depletion, $$\hat H =\sum_{\bf k} E_k\hat a^\dagger_{\bf k}\hat a_{\bf k} +\frac{1}{ {\cal V}}\sum_{\bf k,\bf p} U ({\bf k-\bf p}) \hat a^\dagger_{\bf k} \hat a_{\bf p} +\frac{1}{2 {\cal V}}\sum_{\bf k,\bf p,\bf q} V ({\bf q})\, \hat a^\dagger_{\bf k+\bf q} \hat a^\dagger_{\bf p-\bf q}\hat a_{\bf p}\hat a_{\bf k},$$6 the thermal depletion, and $$\hat H =\sum_{\bf k} E_k\hat a^\dagger_{\bf k}\hat a_{\bf k} +\frac{1}{ {\cal V}}\sum_{\bf k,\bf p} U ({\bf k-\bf p}) \hat a^\dagger_{\bf k} \hat a_{\bf p} +\frac{1}{2 {\cal V}}\sum_{\bf k,\bf p,\bf q} V ({\bf q})\, \hat a^\dagger_{\bf k+\bf q} \hat a^\dagger_{\bf p-\bf q}\hat a_{\bf p}\hat a_{\bf k},$$7 the disorder-induced contribution (Boudjemaa, 5 Aug 2025).

For three-dimensional isotropic speckle disorder with

$$\hat H =\sum_{\bf k} E_k\hat a^\dagger_{\bf k}\hat a_{\bf k} +\frac{1}{ {\cal V}}\sum_{\bf k,\bf p} U ({\bf k-\bf p}) \hat a^\dagger_{\bf k} \hat a_{\bf p} +\frac{1}{2 {\cal V}}\sum_{\bf k,\bf p,\bf q} V ({\bf q})\, \hat a^\dagger_{\bf k+\bf q} \hat a^\dagger_{\bf p-\bf q}\hat a_{\bf p}\hat a_{\bf k},$$8

the disorder deformation in the quadrupolar generalization becomes

$$\hat H =\sum_{\bf k} E_k\hat a^\dagger_{\bf k}\hat a_{\bf k} +\frac{1}{ {\cal V}}\sum_{\bf k,\bf p} U ({\bf k-\bf p}) \hat a^\dagger_{\bf k} \hat a_{\bf p} +\frac{1}{2 {\cal V}}\sum_{\bf k,\bf p,\bf q} V ({\bf q})\, \hat a^\dagger_{\bf k+\bf q} \hat a^\dagger_{\bf p-\bf q}\hat a_{\bf p}\hat a_{\bf k},$$9

where $E_k=\hbar^2k^2/(2m)$0 is the disorder correlation length, $E_k=\hbar^2k^2/(2m)$1, and $E_k=\hbar^2k^2/(2m)$2 contains the angular quadrupolar correction (Boudjemaa, 5 Aug 2025). In the contact limit,

$E_k=\hbar^2k^2/(2m)$3

one recovers $E_k=\hbar^2k^2/(2m)$4, so the quadrupolar theory reduces to the standard Huang-Meng result (Boudjemaa, 5 Aug 2025).

The same framework yields a disorder correction to the ground-state energy,

$E_k=\hbar^2k^2/(2m)$5

For speckle disorder, the paper writes

$E_k=\hbar^2k^2/(2m)$6

so the disorder contribution lowers the ground-state energy (Boudjemaa, 5 Aug 2025).

The approximation is explicitly accompanied by a validity criterion. At $E_k=\hbar^2k^2/(2m)$7, it requires

$E_k=\hbar^2k^2/(2m)$8

or equivalently

$E_k=\hbar^2k^2/(2m)$9

These inequalities encode, respectively, the small-depletion requirement and the weak-disorder requirement (Boudjemaa, 5 Aug 2025).

4. Extensions: anisotropic interactions, nonequilibrium quenches, and trapped clouds

A recent extension formulates the approximation for a homogeneous three-dimensional quadrupolar Bose gas. In momentum space the interaction is

$\langle U(\mathbf r)\rangle=0$0

and the Bogoliubov spectrum becomes

$\langle U(\mathbf r)\rangle=0$1

The resulting spectrum is anisotropic at higher momentum because of the $\langle U(\mathbf r)\rangle=0$2 term, but at small $\langle U(\mathbf r)\rangle=0$3 the sound velocity remains isotropic: $\langle U(\mathbf r)\rangle=0$4 This contrasts with dipolar condensates, for which the low-$\langle U(\mathbf r)\rangle=0$5 sound velocity is anisotropic (Boudjemaa, 5 Aug 2025).

The same paper extends the construction to quench dynamics under a time-dependent Bogoliubov-Huang-Meng approximation. After a sudden interaction change, one writes

$\langle U(\mathbf r)\rangle=0$6

and follows the time evolution of the interaction-induced depletion

$\langle U(\mathbf r)\rangle=0$7

and the disorder deformation

$\langle U(\mathbf r)\rangle=0$8

Within that quadratic theory, both quantities are enhanced in the asymptotic steady state relative to equilibrium, while $\langle U(\mathbf r)\rangle=0$9 exhibits damped oscillations whose amplitudes depend strongly on the disorder correlation length and the relative quadrupolar interaction (Boudjemaa, 5 Aug 2025).

In trapped dirty-boson experiments, the perturbative Huang-Meng sector has also been tested indirectly through static cloud geometry. For a harmonically trapped condensate of $$\langle U(\mathbf r)U(\mathbf r')\rangle=R(\mathbf r-\mathbf r')$$0 molecules in laser speckle, quantitative agreement with the perturbative Huang-Meng approach is reported for small disorder strengths, roughly

$$\langle U(\mathbf r)U(\mathbf r')\rangle=R(\mathbf r-\mathbf r')$$1

where the geometric mean of the measured transverse widths is reproduced well (Nagler et al., 2019). At stronger disorder, a nonperturbative disorder treatment becomes necessary, and even then the use of a local-density approximation yields a constant theoretical aspect ratio while the measured cloud aspect ratio decreases monotonously with increasing disorder strength. The discrepancy is attributed to failure of the local-density approximation in the strong-disorder regime (Nagler et al., 2019).

5. Mathematical status and rigorous underpinnings

The Bogoliubov-Huang-Meng approximation is not itself a fully rigorous disorder theory, but two nearby mathematical results clarify its foundations. First, for an interacting Bose gas in a homogeneous random medium, the generalized Bogoliubov $$\langle U(\mathbf r)U(\mathbf r')\rangle=R(\mathbf r-\mathbf r')$$2-number substitution over the low-energy set

$$\langle U(\mathbf r)U(\mathbf r')\rangle=R(\mathbf r-\mathbf r')$$3

is thermodynamically exact for the pressure after variational optimization and the double limit $$\langle U(\mathbf r)U(\mathbf r')\rangle=R(\mathbf r-\mathbf r')$$4, $$\langle U(\mathbf r)U(\mathbf r')\rangle=R(\mathbf r-\mathbf r')$$5: $$\langle U(\mathbf r)U(\mathbf r')\rangle=R(\mathbf r-\mathbf r')$$6 This result is motivated by generalized, especially type III, Bose-Einstein condensation in random media and provides a rigorous underpinning for the condensate-substitution step, though it does not derive Huang-Meng depletion formulas or disorder-corrected spectra (Jaeck et al., 2010).

Second, the dynamical validity of the Bogoliubov approximation without disorder has been established for $$\langle U(\mathbf r)U(\mathbf r')\rangle=R(\mathbf r-\mathbf r')$$7-boson Schrödinger evolution with interaction scaling $$\langle U(\mathbf r)U(\mathbf r')\rangle=R(\mathbf r-\mathbf r')$$8, $$\langle U(\mathbf r)U(\mathbf r')\rangle=R(\mathbf r-\mathbf r')$$9. After factoring out the condensate, the exact fluctuation vector satisfies

$\hat a_0=\hat a_0^\dagger=\sqrt{N}.$0

while the approximating quasifree fluctuation vector satisfies

$\hat a_0=\hat a_0^\dagger=\sqrt{N}.$1

with a quadratic Bogoliubov Hamiltonian containing both $\hat a_0=\hat a_0^\dagger=\sqrt{N}.$2 and pairing terms. The exact many-body state is then approximated in norm by a condensate plus Bogoliubov fluctuations,

$\hat a_0=\hat a_0^\dagger=\sqrt{N}.$3

This result is disorder-free, but it rigorously justifies the condensate-plus-quasi-free-fluctuations structure on which Huang-Meng-type approximations also rely (Nam et al., 2016).

Taken together, these two strands suggest a layered picture: the generalized $\hat a_0=\hat a_0^\dagger=\sqrt{N}.$4-number substitution in random media can be exact for thermodynamic pressure, while the quadratic Bogoliubov fluctuation theory can be dynamically valid in clean systems over an explicit scaling regime. The full disorder-sensitive Bogoliubov-Huang-Meng approximation lies between these rigorous endpoints.

6. Relation to neighboring approximations, misconceptions, and limits of validity

A recurrent source of confusion is the word “Huang.” In clean dilute-gas theory it often refers to the Huang in Lee-Huang-Yang, as in rigorous results showing that the Bogoliubov approximation captures the correct second-order correction to the ground-state energy of the dilute Bose gas. Such works address the clean-gas LHY correction and contain no disorder, no random external potential, and no Huang-Meng depletion sector (Brietzke et al., 2019). Similarly, self-consistent clean-system derivations of modified Gross-Pitaevskii equations with depletion and anomalous density are built on the same condensate-plus-fluctuation logic, but their potential $\hat a_0=\hat a_0^\dagger=\sqrt{N}.$5 is a generic external trap rather than a quenched random field (Salasnich, 2018).

Within its own domain, the approximation has clear limits. The disordered quadrupolar formulation states them explicitly: weak disorder only, small fluctuations and small depletion, low temperature, a homogeneous condensate background treated perturbatively, and a quadratic theory that yields prethermal steady states after a quench rather than true long-time thermalization (Boudjemaa, 5 Aug 2025). The trapped dirty-boson study adds an independent limitation: even when disorder is treated beyond simple perturbation theory, a local-density approximation can fail for strong disorder and anisotropy-sensitive observables such as the cloud aspect ratio (Nagler et al., 2019).

A further misconception is to interpret the exactness of $\hat a_0=\hat a_0^\dagger=\sqrt{N}.$6-number substitution in random media as a full derivation of the Huang-Meng approximation. The rigorous random-media result concerns thermodynamic exactness of a generalized substitution for the pressure and the relation of the maximizing $\hat a_0=\hat a_0^\dagger=\sqrt{N}.$7-numbers to generalized condensation, not the calculation of disorder-induced depletion, disorder-corrected excitation spectra, or weak-disorder quasiparticle occupations (Jaeck et al., 2010). Conversely, a clean Bogoliubov or HFB derivation with anomalous density and Lee-Huang-Yang terms can illuminate the fluctuation backbone of the theory, but it remains silent about the Huang-Meng disorder component unless a quenched random potential is added explicitly (Salasnich, 2018).

In its most precise current usage, then, the Bogoliubov-Huang-Meng approximation is a perturbative dirty-boson scheme in which the condensate is treated at mean-field level, the excitations at quadratic Bogoliubov level, and the random potential through a lowest-order disorder translation. Its chief outputs are the quasiparticle spectrum, the interaction depletion, the disorder-induced condensate deformation $\hat a_0=\hat a_0^\dagger=\sqrt{N}.$8, and related thermodynamic or dynamical observables. Its chief restrictions are the weak-disorder, small-depletion, and low-temperature assumptions under which the quadratic closure remains reliable.