---
title: Bogoliubov-Huang-Meng Approximation
url: https://www.emergentmind.com/topics/bogoliubov-huang-meng-approximation
type: topic
---

# Bogoliubov-Huang-Meng Approximation

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 [2508.03791]. 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 [1911.02626].

## 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 [2508.03791].

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 [1010.0425]. 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 [1901.00537].

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')\) [2508.03791].

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 [2508.03791].

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 b_{\bf k}\) are quasiparticle operators and \(\beta_{\bf k}\) is a disorder-induced shift. After this step,
\[
\hat H = E+\sum_{\mathbf k\neq0}\varepsilon_k\, \hat b_{\bf k}^\dagger \hat b_{\bf k}.
\]
In this sense, the Bogoliubov-Huang-Meng approximation is literally a Bogoliubov diagonalization of the interacting condensate plus a Huang-Meng disorder translation [2508.03791].

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

A defining equilibrium observable is the disorder-induced contribution to the noncondensed density,
\[
n_R={ \cal V}^{-1}\sum_{\mathbf k \neq 0} \langle |\beta_{\bf k}|^2 \rangle = n \int \frac{ d^3k}{ (2\pi)^3} \frac{ E_k^2}{ \varepsilon_{k}^4} R_{\bf k}.
\]
This quantity is identified as the “density of the disorder-averaged condensate” or “condensate deformation,” and it is distinct from interaction-induced quantum depletion [2508.03791]. The total depletion is decomposed as
\[
\tilde n=\tilde n_0+\tilde n_T+n_R,
\]
with \(\tilde n_0\) the zero-temperature interaction depletion, \(\tilde n_T\) the thermal depletion, and \(n_R\) the disorder-induced contribution [2508.03791].

For three-dimensional isotropic speckle disorder with
\[
R(\mathbf r)=R\,\mathrm{sinc}^2(\mathbf r/\sigma),
\qquad
R_{\mathbf k}= \frac{\pi ^2 R}{ k \sigma} \Theta\left(2- k\sigma\right),
\]
the disorder deformation in the quadrupolar generalization becomes
\[
n_R = n_{\text{HM}\, f_R(\sigma/\xi,\beta(\theta)),
\qquad
n_{\text{HM}= \frac{R \, m^2 \sqrt{n/a}} {8\pi^{3/2} \hbar^4},
\]
where \(\sigma\) is the disorder correlation length, \(\xi=\hbar/\sqrt{mng}\), and \(\beta(\theta)\) contains the angular quadrupolar correction [2508.03791]. In the contact limit,
\[
\epsilon_{\mathrm{qq}}=0,\qquad \beta(\theta)=\frac14,
\]
one recovers \(n_R\to n_{\mathrm{HM}}\), so the quadrupolar theory reduces to the standard Huang-Meng result [2508.03791].

The same framework yields a disorder correction to the ground-state energy,
\[
E=E_0+E_{\mathrm{LHY}+E_R,
\qquad
\frac{E_0}{\cal V}=\frac{n^2g}{2},
\qquad
\frac{E_R}{\cal V}=- n \int \frac{d^3 k}{(2\pi)^3}  R_{\bf k} \frac{ E_k}{\varepsilon_{k}^2}.
\]
For speckle disorder, the paper writes
\[
\frac{E_R}{N}=-  \frac{R m \sqrt{na/\pi}}{\hbar^2}  h_R( \beta(\theta)),
\]
so the disorder contribution lowers the ground-state energy [2508.03791].

The approximation is explicitly accompanied by a validity criterion. At \(T=0\), it requires
\[
\frac{f_0(\beta(\theta))}{8\pi^2\xi^3}\ll 1,
\qquad
n_{\mathrm{HM}} f_R(\sigma/\xi,\beta(\theta))\ll 1,
\]
or equivalently
\[
R<R_c= \frac{8\pi^{3/2}\hbar^4} {m^2\sqrt{n/a}\,f_R(\sigma/\xi,\beta(\theta))}.
\]
These inequalities encode, respectively, the small-depletion requirement and the weak-disorder requirement [2508.03791].

## 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
\[
V(\mathbf k)= g \bigg[ 1+\frac{4 \pi}{105} \epsilon_{\text{qq}}\, Y_4^0 (\theta)\, k^2 \bigg],
\]
and the Bogoliubov spectrum becomes
\[
\varepsilon_{k}= \sqrt{E_k^2 + 2E_k n g \left[1+ \frac{4 \pi \epsilon_{\text{qq}}}{105} Y_4^0 (\theta_{\bf k}) k^2 \right]}.
\]
The resulting spectrum is anisotropic at higher momentum because of the \(Y_4^0(\theta_{\bf k})\) term, but at small \(k\) the sound velocity remains isotropic:
\[
\varepsilon_k=\hbar c_s k,
\qquad
c_s=\sqrt{\frac{gn}{m}}.
\]
This contrasts with dipolar condensates, for which the low-\(k\) sound velocity is anisotropic [2508.03791].

The same paper extends the construction to quench dynamics under a time-dependent Bogoliubov-Huang-Meng approximation. After a sudden interaction change, one writes
\[
\hat a_{\bf k} (t)= u_k(t) \hat b_{\bf k}-v_k(t) \hat b^\dagger_{-\bf k}-\beta_{\bf k} (t),
\]
and follows the time evolution of the interaction-induced depletion
\[
\tilde n_0(t)=\frac1{\cal V}\sum_{\mathbf k\neq0}|v_k(t)|^2
\]
and the disorder deformation
\[
n_R(t)=\frac1{\cal V}\sum_{\mathbf k\neq0}\langle |\beta_{\bf k}(t)|^2\rangle.
\]
Within that quadratic theory, both quantities are enhanced in the asymptotic steady state relative to equilibrium, while \(n_R(t)\) exhibits damped oscillations whose amplitudes depend strongly on the disorder correlation length and the relative quadrupolar interaction [2508.03791].

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 \(^{6}\mathrm{Li}_2\) molecules in laser speckle, quantitative agreement with the perturbative Huang-Meng approach is reported for small disorder strengths, roughly
\[
\bar u/k_B \lesssim 30\,\mathrm{nK},
\]
where the geometric mean of the measured transverse widths is reproduced well [1911.02626]. 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 [1911.02626].

## 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 \(c\)-number substitution over the low-energy set
\[
I_\delta=\left\{k\in\Lambda_l^*:\frac{k^2}{2}\le \delta\right\}
\]
is thermodynamically exact for the pressure after variational optimization and the double limit \(l\to\infty\), \(\delta\downarrow0\):
\[
\text{a.s.}\quad \lim_{l\to\infty} p_l(\beta,\mu) = \lim_{\delta\downarrow 0}\,\liminf_{l\to\infty}\, \max_{\{c_k\}} p_{l,\delta}^{\mathrm{low}}(\mu,\{c_k\}) = \lim_{\delta\downarrow 0}\,\limsup_{l\to\infty}\, \max_{\{c_k\}} p_{l,\delta}^{\mathrm{low}}(\mu,\{c_k\}).
\]
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 [1010.0425].

Second, the dynamical validity of the Bogoliubov approximation without disorder has been established for \(N\)-boson Schrödinger evolution with interaction scaling \(N^{3\beta-1}w(N^\beta(x-y))\), \(0\le\beta<1/2\). After factoring out the condensate, the exact fluctuation vector satisfies
\[
i\partial_t \Phi_N(t)=\widetilde H_N(t)\Phi_N(t),
\]
while the approximating quasifree fluctuation vector satisfies
\[
i\partial_t\Phi(t)=\mathbb H(t)\Phi(t),
\]
with a quadratic Bogoliubov Hamiltonian containing both \(d\Gamma(h(t))\) and pairing terms. The exact many-body state is then approximated in norm by a condensate plus Bogoliubov fluctuations,
\[
\left\| \Psi_N(t)- \sum_{n=0}^N u(t)^{\otimes (N-n)}\otimes_s \psi_n(t) \right\|_{\mathfrak H^N}^2 \le C_\varepsilon (1+t)^{1+\varepsilon} N^{(2\beta-1+\varepsilon)/2}.
\]
This result is disorder-free, but it rigorously justifies the condensate-plus-quasi-free-fluctuations structure on which Huang-Meng-type approximations also rely [1604.05240].

Taken together, these two strands suggest a layered picture: the generalized \(c\)-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 [1901.00537]. 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 \(U(\mathbf r)\) is a generic external trap rather than a quenched random field [1810.02636].

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 [2508.03791]. 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 [1911.02626].

A further misconception is to interpret the exactness of \(c\)-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 \(c\)-numbers to generalized condensation, not the calculation of disorder-induced depletion, disorder-corrected excitation spectra, or weak-disorder quasiparticle occupations [1010.0425]. 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 [1810.02636].

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 \(n_R\), 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.

Source: https://www.emergentmind.com/topics/bogoliubov-huang-meng-approximation