Thermal Bogoliubov Transformation

• Thermal Bogoliubov Transformation is a finite-temperature method that generalizes the canonical transformation by incorporating thermal occupation and mixing operators in a doubled Hilbert space.
• It employs a two-mode squeezing operator to construct a thermal vacuum state, facilitating the diagonalization of quasiparticle operators and accurate statistical descriptions at nonzero temperatures.
• This transformation is pivotal in thermo-field dynamics, mean-field theories, and real-time quantum field formulations, underpinning calculations of pairing correlations and collective excitations.

The thermal Bogoliubov transformation is a universal structure underpinning finite-temperature quantum many-body theory, central to formalisms such as thermo-field dynamics (TFD), resonating Hartree-Bogoliubov (Res-HB) theory, and the real-time/imaginary-time formulations of quantum field theory. It generalizes the canonical Bogoliubov transformation by encoding thermal occupation, enabling pure-state representations and diagonal quasinormal bases at nonzero temperature. The transformation operates at the operator level, in many cases as a two-mode squeezing transformation in a doubled (real+tilde) Hilbert space, and re-expresses the statistical density matrix as a vacuum state with respect to temperature-dependent quasiparticle operators.

1. Formal Structure of the Thermal Bogoliubov Transformation

The thermal Bogoliubov transformation acts on a doubled Hilbert space $\mathcal{H}\otimes\widetilde{\mathcal{H}}$. For each original mode (bosonic or fermionic), a “tilde” partner is introduced. In TFD, the generator is a two-mode squeezing operator

$$B(\theta) = \exp\left[\theta\,(\hat{a}\,\tilde{a} - \hat{a}^\dagger \tilde{a}^\dagger)\right]$$

where $\theta$ is a real, temperature-dependent parameter determined by requiring that the physical reduced density matrix matches the thermal (Gibbs) state. For a bosonic mode,

$$\theta(\beta) = \mathrm{arctanh}\left(e^{-\beta\hbar\omega/2}\right), \quad \langle \hat{a}^\dagger \hat{a} \rangle = \sinh^2\theta = (e^{\beta\hbar\omega} - 1)^{-1}$$

Analogous formulae exist for fermions, with sign and distribution changes reflecting Fermi-Dirac statistics. The Bogoliubov transformation mixes annihilation and creation operators from both original and tilde spaces: $$B(\theta)\,\hat{a}\,B^\dagger(\theta) = \hat{a}\,\cosh\theta + \tilde{a}^\dagger\,\sinh\theta$$ The corresponding “thermal vacuum” $|0(\beta)\rangle$ in the doubled space is annihilated by all temperature-dependent quasiparticle operators and yields thermal averages of observables (Błasiak et al., 27 May 2025, Vdovin et al., 2010, Hou et al., 2018).

2. Implementation in Thermo-Field Dynamics (TFD) and Hilbert Space Doubling

TFD is formulated by doubling the Hilbert space and introducing both real (physical) and tilde (auxiliary) operators, subject to commutation/anticommutation and tilde conjugation rules. The thermal Hamiltonian

$\mathcal{H} = H - \widetilde{H}$

drives real-time evolution in the doubled space. The thermal Bogoliubov transformation generates new temperature-dependent modes

$$\begin{aligned} a_k(T) &= u_k(T)\,a_k - v_k(T)\,\tilde{a}^\dagger_k \ \tilde{a}_k(T) &= u_k(T)\,\tilde{a}_k + v_k(T)\,a^\dagger_k \end{aligned}$$

where $u_k(T)^2 + v_k(T)^2 = 1$ and $v_k(T)^2 = 1/(1 + e^{E_k/T})$ for fermions. The unique thermal vacuum $$B(\theta) = \exp\left[\theta\,(\hat{a}\,\tilde{a} - \hat{a}^\dagger \tilde{a}^\dagger)\right]$$0 is annihilated by all $$B(\theta) = \exp\left[\theta\,(\hat{a}\,\tilde{a} - \hat{a}^\dagger \tilde{a}^\dagger)\right]$$1 and $$B(\theta) = \exp\left[\theta\,(\hat{a}\,\tilde{a} - \hat{a}^\dagger \tilde{a}^\dagger)\right]$$2 and forms the starting point for formulating perturbation theory, reduced density matrices, and observable distributions at finite $$B(\theta) = \exp\left[\theta\,(\hat{a}\,\tilde{a} - \hat{a}^\dagger \tilde{a}^\dagger)\right]$$3 (Vdovin et al., 2010, Błasiak et al., 27 May 2025, Hou et al., 2018).

3. Thermal Bogoliubov Transformation in Mean-Field and Many-Body Theory

In mean-field approaches and Hartree-Bogoliubov (HB) theory, the thermal Bogoliubov transformation enables a consistent description of statistical occupations across the quasiparticle spectrum. The Res-HB theory extends the standard mean-field construction by working in a subspace spanned by superpositions of non-orthogonal HB wavefunctions, each itself obtained by a temperature-dependent Bogoliubov transformation that diagonalizes the thermal Fock-Bogoliubov operator $$B(\theta) = \exp\left[\theta\,(\hat{a}\,\tilde{a} - \hat{a}^\dagger \tilde{a}^\dagger)\right]$$4 (Nishiyama et al., 2013).

The thermal HB density matrix in each resonating sector takes the form

$$B(\theta) = \exp\left[\theta\,(\hat{a}\,\tilde{a} - \hat{a}^\dagger \tilde{a}^\dagger)\right]$$5

where $$B(\theta) = \exp\left[\theta\,(\hat{a}\,\tilde{a} - \hat{a}^\dagger \tilde{a}^\dagger)\right]$$6 is constructed from the variational derivatives of the energy functional with respect to the normal and anomalous thermal densities. A Bogoliubov rotation $$B(\theta) = \exp\left[\theta\,(\hat{a}\,\tilde{a} - \hat{a}^\dagger \tilde{a}^\dagger)\right]$$7 simultaneously diagonalizes $$B(\theta) = \exp\left[\theta\,(\hat{a}\,\tilde{a} - \hat{a}^\dagger \tilde{a}^\dagger)\right]$$8 and $$B(\theta) = \exp\left[\theta\,(\hat{a}\,\tilde{a} - \hat{a}^\dagger \tilde{a}^\dagger)\right]$$9, yielding temperature-dependent quasiparticle energies $\theta$0 and occupation numbers $\theta$1 (Nishiyama et al., 2013). This structure generalizes the zero-temperature prescription, with temperature entering both the spectrum and the population statistics.

4. Applications in Real-Time Formalism and Quantum Field Theory

In the real-time (Schwinger-Keldysh) formalism, the thermal Bogoliubov transformation acts as a $\theta$2 matrix mixing on thermal doublet fields. For the closed-time-path ($\theta$3), the finite-temperature propagator $\theta$4 is generated by

$\theta$5

where $\theta$6 is explicitly given via the Bose-Einstein or Fermi-Dirac distributions and depends only on $\theta$7. In this case, the temperature dependence of all Green’s functions can be factorized via a scalar thermal operator $\theta$8: $\theta$9 with $$\theta(\beta) = \mathrm{arctanh}\left(e^{-\beta\hbar\omega/2}\right), \quad \langle \hat{a}^\dagger \hat{a} \rangle = \sinh^2\theta = (e^{\beta\hbar\omega} - 1)^{-1}$$0. For general real-time contours $$\theta(\beta) = \mathrm{arctanh}\left(e^{-\beta\hbar\omega/2}\right), \quad \langle \hat{a}^\dagger \hat{a} \rangle = \sinh^2\theta = (e^{\beta\hbar\omega} - 1)^{-1}$$1, the structure becomes more complex, requiring an extended “generalized zero-temperature” propagator before the scalar operator applies (Das et al., 2018).

5. Reduced Density Matrices, Wigner Distributions, and Generalizations

Within TFD and the iBT (inverse Bogoliubov transformation) variant, the thermal state is treated as either a transformed vacuum or by moving the transformation into the propagators. The reduced one-particle density matrix (1-RDM) and Wigner phase-space distributions are extracted by tracing out tilde components and back-transforming via the Bogoliubov operator: $$\theta(\beta) = \mathrm{arctanh}\left(e^{-\beta\hbar\omega/2}\right), \quad \langle \hat{a}^\dagger \hat{a} \rangle = \sinh^2\theta = (e^{\beta\hbar\omega} - 1)^{-1}$$2 For harmonic oscillators, analytic results recover thermal Gaussian distributions, while in anharmonic or multi-mode settings, approximate schemes (correlation neglect or moment expansion) are necessary (Błasiak et al., 27 May 2025). These techniques provide a route to finite-temperature quantum statistics within a pure-state wavefunction framework.

6. Physical Interpretation and Impact

The thermal Bogoliubov transformation establishes a direct connection between operator algebraic properties and finite-temperature statistical mechanics. In mean-field and TFD approaches, it enables the explicit construction of thermal vacua, squeezing/coherent state structures, and systematic perturbation theory at nonzero $$\theta(\beta) = \mathrm{arctanh}\left(e^{-\beta\hbar\omega/2}\right), \quad \langle \hat{a}^\dagger \hat{a} \rangle = \sinh^2\theta = (e^{\beta\hbar\omega} - 1)^{-1}$$3. It underlies the calculation of thermodynamic quantities, pairing correlations, and the emergence of phenomena such as the pseudogap in fermionic superfluids (Hou et al., 2018), and is vital for studying collective excitations and damping in hot nuclei beyond standard random-phase approximations (Vdovin et al., 2010). In the resonating mean-field context, the transformation supports a rigorous and generalizable mean-field description capable of treating both normal and anomalous thermal correlations (Nishiyama et al., 2013).

7. Computational Strategies and Practical Approximations

Exact evaluation of thermal partition functions in resonating spaces typically leads to computationally prohibitive continued fractions, as each level introduces increasingly complex string operators in the projected subspace. Approximations such as diagonal dominance—whereby the full Hamiltonian is replaced by a sum over quadratic forms defined by the thermal FB operators for each resonating configuration—allow efficient calculation of properties while capturing the essential physics of low-lying configurations (Nishiyama et al., 2013). In high-dimensional TFD applications, neglecting non-trivial mode correlations or reconstructing densities via Hermite-moment expansions provides scalability in complex quantum systems (Błasiak et al., 27 May 2025).

---

Key References: (Nishiyama et al., 2013, Das et al., 2018, Hou et al., 2018, Błasiak et al., 27 May 2025, Vdovin et al., 2010)