Derivation of renormalized Hartree-Fock-Bogoliubov and quantum Boltzmann equations in an interacting Bose gas

Technical Summary: Derivation of Renormalized Hartree-Fock-Bogoliubov and Quantum Boltzmann Equations in an Interacting Bose Gas

Problem Statement
The paper addresses the rigorous derivation of effective kinetic equations for an interacting Bose gas, specifically focusing on the emergence of the Quantum Boltzmann Equation (QBE) from the many-body Schrödinger dynamics in the presence of a Bose-Einstein Condensate (BEC). While previous work by the authors [37] established the derivation of QBEs near a BEC, it relied on specific assumptions regarding the propagation of factorization (quasifreeness) and yielded error bounds that limited the time of validity to $t \sim (\log N / \log \log N)^2$. The central challenge is to characterize the leading-order fluctuation dynamics more precisely, separating the fast oscillating Hartree-Fock-Bogoliubov (HFB) dynamics from the slow, dissipative Boltzmann dynamics, and to improve the error estimates to extend the validity of the derived equations.

Methodology
The authors employ a second-quantization framework on a 3-torus $\Lambda$ with $N$ particles. The analysis centers on the decomposition of the full many-body state into a BEC component and thermal fluctuations. The methodology involves the following key steps:

1. Transformation and Fluctuation Dynamics: The authors introduce a unitary transformation to the "relative" evolution, defined by $U_{\text{fluc}}(t) = e^{iS_t} U_{\text{Bog}}^\dagger(t) T^\dagger[k_t] W^\dagger[\sqrt{N|\Lambda|}\phi_t] e^{-itH_N} W[\sqrt{N|\Lambda|}\phi_0] T[k_0]$. Here, $W$ is the Weyl operator (shifting the condensate), $T$ is the Bogoliubov rotation (handling pair correlations), and $U_{\text{Bog}}$ is the Bogoliubov propagation (handling the dispersion).
2. Renormalization Strategy: The core innovation is a recursive renormalization procedure. By expanding the total density and pair correlations in powers of $N^{-1/2}$, the authors identify terms in the Duhamel expansion of the fluctuation dynamics that correspond to HFB contributions.
   • First Order: Terms of order $N^{-1/2}$ in the BEC wave function evolution are eliminated by choosing the Weyl shift $\phi_t$ to satisfy a specific equation.
   • Second Order: Terms of order $N^{-1}$ in the pair-correlation evolution are eliminated by choosing the Bogoliubov parameters $(\gamma_t, \sigma_t)$ and the dispersion $\Omega_t$ to satisfy renormalized HFB equations.
   • Higher Orders: This process is iterated to absorb contributions into the definition of the renormalized fields $(\phi, \gamma, \sigma, \Omega)$.
3. Separation of Scales: After renormalization, the remaining terms in the fluctuation dynamics are shown to be purely of the Quantum Boltzmann type (cubic collision terms). The authors demonstrate that the "error" terms identified in their previous work [37] are actually contributions to the renormalized HFB equations.
4. A Priori Estimates: To control the remainder terms, the authors establish global well-posedness for the renormalized HFB system. They utilize symplectic descriptions and energy/mass conservation laws to derive uniform bounds on the HFB fields, which are then used to bound the tail of the perturbation expansion.

Key Contributions and Results

• Renormalized HFB Equations: The paper derives a set of renormalized HFB equations (Equation 2.16) that include corrections from the initial thermal state and the BEC density. These equations govern the leading-order dynamics of the condensate wave function and the pair correlations.
• Derivation of Quantum Boltzmann Equations: It is shown that the subleading dynamics of the thermal fluctuations and the BEC wave function are governed by cubic Quantum Boltzmann collision operators ($Q_3$). Specifically:
  • The thermal density $f_t$ evolves according to a cubic Boltzmann equation with a collision kernel dependent on the renormalized HFB fields.
  • The BEC wave function $\Phi_t$ and the pair-absorption rate $g_t$ evolve according to equations driven by these same collision terms, with error bounds of order $O(N^{-3/2})$ and $O(N^{-2})$ respectively.
• Improved Error Bounds and Time of Validity:
  • The error bounds in the main theorem (Theorem 2.8) are sharp in powers of $N$.
  • Crucially, the time of validity for the derived equations is extended from $t \sim (\log N / \log \log N)^2$ (as in [37]) to $t \sim (\log N)^2$. This improvement is achieved by correctly identifying and absorbing previously unaccounted terms into the renormalized HFB dynamics.
• Global Well-Posedness: The authors prove the global well-posedness of the renormalized HFB system (Proposition 2.7) under the assumption of a non-negative interaction potential, ensuring the stability of the leading-order dynamics required for the derivation.

Significance
The paper claims to provide a more complete and rigorous characterization of the fluctuation dynamics around a BEC. By distinguishing between HFB contributions and Boltzmann collision terms through a systematic renormalization strategy, the authors resolve ambiguities in the error terms of previous derivations. The extension of the validity time to $(\log N)^2$ represents a significant improvement in the mathematical control of the many-body dynamics in the mesoscopic regime. The work confirms the phenomenological paradigm that HFB dynamics (fast oscillations) and QBE dynamics (slow relaxation) evolve on distinct time scales. The authors note that their approach is expected to be extensible to smaller orders in $1/N$, though terms emerging at order $N^{-2}$ may not be purely of HFB or Boltzmann type.

The results are presented as unconditional for sufficiently short times, relying on the propagation of restricted quasifreeness, which is rigorously proved for the derived time window. The work builds upon and corrects the authors' previous derivation [37] and aligns with the approach of Grillakis-Machedon et al. [60] regarding Bogoliubov rotations, while extending it to include the specific quasifree state structure necessary for the emergence of the Boltzmann equation.