Squeezed Thermal States in Quantum Systems

Updated 5 February 2026

Squeezed thermal states (STS) are continuous-variable Gaussian quantum states created by applying a squeezing operator to a thermal state, balancing quantum noise and thermal disorder.
STS exhibit reduced quadrature noise and constant entropy, making them pivotal for applications in quantum optics, quantum thermodynamics, and engineered reservoir design.
Experimental protocols using trap-frequency modulation, nonlinear interactions, and machine learning tomography enable precise generation and real-time characterization of STS.

A squeezed thermal state (STS) is a fundamental continuous-variable Gaussian quantum state formed by applying a squeezing operation to a thermal (Gibbs) state of a bosonic mode. STSs generalize both thermal and pure squeezed-vacuum states, capturing essential features of quantum noise suppression and classical thermal disorder. They serve as a central resource in quantum optics, open quantum systems, quantum thermodynamics, optomechanics, and relativistic quantum information.

1. Mathematical Definition and Structure

Let $\rho_{\rm th}$ denote the thermal (Gibbs) state of a harmonic oscillator at frequency $\omega$ and inverse temperature $\beta=(k_BT)^{-1}$ : $\rho_{\rm th} = \frac{e^{-\beta\hbar\omega (a^\dagger a+\frac12)}}{Z}, \qquad Z=\mathrm{Tr}[e^{-\beta\hbar\omega(a^\dagger a+\frac12)}]$ or in the photon-number basis,

$\rho_{\rm th} = \sum_{k=0}^\infty \frac{n_{\rm th}^k}{(n_{\rm th}+1)^{k+1}}\,|k\rangle\langle k|, \qquad n_{\rm th} = \frac{1}{e^{\beta\hbar\omega}-1}.$

The unitary single-mode squeezing operator

$S(r,\theta) = \exp\left[\frac{r}{2}\left(e^{-i\theta} a^2 - e^{i\theta} (a^\dagger)^2\right)\right]$

acts on $\rho_{\rm th}$ to yield the STS: $\rho_{\rm STS} = S(r,\theta)\, \rho_{\rm th}\, S^\dagger(r,\theta).$

STS are zero-displacement Gaussian states described by squeezing amplitude $r \ge 0$ , squeezing phase $\theta$ , and thermal mean occupation $n_{\rm th}$ .

The covariance matrix in the Hermitian quadratures

$X = \frac{1}{\sqrt{2}}(a+a^\dagger), \qquad P = \frac{1}{\sqrt{2}i}(a-a^\dagger)$

takes the form (with $\theta=0$ ): $V = \begin{pmatrix} (n_{\rm th}+\frac12)\,e^{-2r} & 0 \ 0 & (n_{\rm th}+\frac12)\,e^{+2r} \end{pmatrix}.$ The Wigner function is a zero-centered ellipse with widths set by the covariance matrix.

An equivalent Fock-space representation is available via the action of the squeezing operator on the thermal density matrix, with explicit matrix elements involving Hermite polynomials and Legendre polynomials in the context of photon addition/subtraction (Hu et al., 2011, Hu et al., 2011).

2. Physical Properties and Statistical Measures

Quadrature Squeezing:

Squeezing modulates the variances of $X$ and $P$ ,

$\langle X^2\rangle_{\rm STS} = (n_{\rm th}+\frac12)e^{-2r}, \qquad \langle P^2\rangle_{\rm STS} = (n_{\rm th}+\frac12)e^{2r}.$

Genuine squeezing is present if $\langle (\Delta Q)^2\rangle < \frac12$ for some quadrature $Q$ (Sewran et al., 2014).

Purity and Entropy:

The purity is

$\mu = {\rm Tr}[\rho_{\rm STS}^2] = \frac{1}{2n_{\rm th}+1}$

(unchanged by $r$ ). The von Neumann entropy depends only on $n_{\rm th}$ : $S(\rho_{\rm STS}) = (n_{\rm th}+1)\ln(n_{\rm th}+1) - n_{\rm th}\ln n_{\rm th}.$ Thus, squeezing is entropy-preserving.

Photon-Number Statistics:

The Fock basis populations interpolate between a Bose–Einstein (thermal) distribution and an even-only pure squeezed vacuum (Tan et al., 2014, Hughes et al., 2024). Explicit probability expressions involve convolutions of thermal and squeezed vacuum statistics.

Second-Order Coherence:

The equal-time normalized second-order coherence is (Hughes et al., 2024): $g^{(2)}(0) = 2 + \frac{(2n_{\rm th}+1)^2\sinh^2(2r)}{[(2n_{\rm th}+1)\cosh(2r)-1]^2}.$ For large thermal occupation $n_{\rm th}\to\infty$ , $g^{(2)}(0)\to2$ (thermal limit), while for $n_{\rm th}\to0$ , $g^{(2)}(0)\to3$ (squeezed vacuum limit).

3. Generation Protocols and Dynamical Engineering

3.1. Theoretical and Experimental Approaches

Trap-Frequency Modulation: Sudden or engineered changes in the oscillator frequency induce squeezing, transforming a thermal state into an STS. “Shortcuts to adiabaticity” protocols achieve this in finite time by reverse-engineering control waveforms and dissipators, enabling entropy control as well (Dupays et al., 2020).
Nonlinear Interactions: In cavity QED, optomechanics, or trapped-ion systems, parametric Hamiltonians of the form $H_{\rm int}\propto a^2 + a^{\dagger2}$ (two-photon processes) generate squeezing even in the presence of loss and thermal noise; the resulting steady state is an STS (Hughes et al., 2024).
Two-Photon Raman Processes: Two-level ions subject to two-photon driving yield unitary or dissipative generation of STS with independently controlled squeezing strength, phase, and temperature (Dupays et al., 2020).
Levitated Optomechanics: Sudden switching between trap frequencies in levitated nanoparticle systems enables classical squeezing of thermal distributions, and with pre-cooling, the same protocol can reach the quantum STS regime (Rashid et al., 2016).

3.2. Open-System and Non-Hamiltonian Effects

STS naturally emerge as the steady states of quadratic open-system Lindblad master equations with both parametric squeezing and thermalization channels (Hughes et al., 2024). Non-Hamiltonian thermostats such as Nosé–Hoover Chains also effectively generate STS as stationary or transient states in molecular dynamics simulations (Sewran et al., 2014).

The interplay of driving, dissipation, and bath temperature governs the achievable squeezing depth, entropy, and steady-state structure. For example, above a critical pump threshold, the anti-squeezed quadrature diverges, and steady-state squeezing is limited by thermal occupancy (Hughes et al., 2024).

4. Fundamental Limits, Decoherence, and Nonclassicality

Threshold Temperature for Squeezing:

Thermal disorder limits squeezing. In realistic models, there exists a threshold temperature $T_{\rm th}$ : for $T>T_{\rm th}$ , the squeezed quadrature variance cannot fall below the vacuum level $1/2$ (Sewran et al., 2014). For driven oscillator models with bath cutoff $\omega_c=3.93\times10^{13}\,\mathrm{Hz}$ , $T_{\rm th}=311.13$ K, approaching biochemical temperatures.

Decoherence of Non-Gaussian States:

Photon-added STS (PASTS) exhibit negative Wigner function values and sub-Poissonian statistics. Such nonclassicality is more robust to environment-induced decoherence than their photon-subtracted counterparts: the Wigner negativity at the phase-space origin persists longer under loss for PASTS (Hu et al., 2011).

Non-Gaussianity Quantification:

The Hilbert–Schmidt fidelity between photon-added/subtracted STS and the original STS can be calculated in closed form via Legendre polynomials (Hu et al., 2011, Hu et al., 2011), further enabling systematic quantification of distance from Gaussianity and resilience to decoherence channels.

5. Quantum Thermodynamics: Entropy Production and Engineered Reservoirs

STS function as non-equilibrium, phase-sensitive reservoirs within the emerging framework of quantum thermodynamics. When used as thermal baths in information erasure or quantum engines, standard formulations of Landauer's principle must be generalized by introducing an effective Hamiltonian $H_{\rm eff}=S H S^\dagger$ (Xu, 4 Feb 2026). The energy cost of operations—generalized “work” and entropy production—must be calculated with respect to $H_{\rm eff}$ , preserving the validity and non-negativity of the generalized Landauer bound.

Entropy production in quantum channels involving STS reservoirs (e.g., Unruh–DeWitt detectors in relativistic settings) displays explicit spacetime dependence and incorporates the unitary resource provided by squeezing (Xu, 4 Feb 2026).

STS as engineered baths can, depending on their parameters, allow quantum heat engines to transiently exceed the classical Carnot bound under generalized thermodynamic definitions. The framework developed for unitarily transformed (squeezed) reservoirs thus provides analytic control over quantum irreversibility and information–energy conversion in a broad range of applications.

6. Tomographic Characterization and Experimental Realization

Covariance Matrix Tomography via Machine Learning:

Robust reconstruction of the physical covariance matrix of noisy STS from sparse quadrature data is achieved via supervised convolutional neural network architectures, capable of accounting for multimode admixture and providing real-time analysis (Rodrıguez et al., 26 Sep 2025). Such machine-learning-based covariance estimation attains high fidelity (average $\langle F\rangle \approx 0.99$ on simulated data) and is robust to strong admixtures of thermal and squeezed components.

Homodyne and Spectral State Reconstruction:

Experimental schemes using single or multi-channel homodyne detection combined with active cavity stabilization allow full covariance matrix reconstruction of single- and two-mode STS from parametric oscillators (Cialdi et al., 2015). Entanglement and purity are then directly extracted from the reconstructed covariance matrix via Simon or Duan criteria.

Phase-Sensitive Metrology:

STS enable quantum-enhanced phase estimation beyond the standard quantum limit when employed as probes in interferometric settings (Yu et al., 2020, Tan et al., 2014). The precision is set not only by the squeezing parameter $r$ but also by the purity $\mu=1/(2n_{\rm th}+1)$ , with excess thermal noise reducing both achievable phase sensitivity and the width of the quantum advantage interval.

Table: Core STS Properties

Property	Formula/Characterization	Reference
Covariance matrix	$V=\begin{pmatrix}(n_{\rm th}+\tfrac12)e^{-2r}&0\0&(n_{\rm th}+\tfrac12)e^{2r}\end{pmatrix}$	(Hughes et al., 2024)
Purity	$\mu=1/(2n_{\rm th}+1)$	(Tan et al., 2014)
Mean photon number	$\langle a^\dagger a\rangle=(2n_{\rm th}+1)\sinh^2 r+n_{\rm th}$	(Tan et al., 2014)
$g^{(2)}(0)$	$2+\frac{(2n_{\rm th}+1)^2\sinh^2(2r)}{[(2n_{\rm th}+1)\cosh(2r)-1]^2}$	(Hughes et al., 2024)
Threshold temperature $T_{\rm th}$ for squeezing	$T_{\rm th}=311.13$ K (for given $\omega_c$ )	(Sewran et al., 2014)

7. Generalizations, Extensions, and Applications

Multimode and Entangled STS:

Applying multimode (or two-mode) squeezing operators to products of thermal states yields multimode Gaussian STS, fully characterized by covariance matrices with off-diagonal (correlating) entries (Cialdi et al., 2015). Entanglement in such states is governed by the competition between squeezing and thermal noise; partial transpose and Duan inequalities provide necessary and sufficient conditions for Gaussian entanglement.

Photon Engineering:

Photon-addition and subtraction applied to STS yields non-Gaussian states whose photon-number statistics, normalization, and Wigner function can be expressed in closed-form via Legendre polynomials, exploiting the thermofield double formalism (Hu et al., 2011, Hu et al., 2011).

Quantum Cosmology and Relativistic Fields:

STS formalism underpins quantum descriptions of the inflaton in anisotropic Bianchi type-I cosmologies, where both quantum fluctuations (squeezing) and thermal excitations drive particle production and impact the semiclassical Einstein equations (Chand, 2022).

Quantum Sensing and Metrology:

State-of-the-art gravitational wave detectors, quantum optical sensors, and atomic interferometers employ STS or its limiting cases, exploiting the metrological advantage conferred by squeezing, with thermal admixture quantifying the unavoidable environmental decoherence (Yu et al., 2020). The role of purity $\mu$ is critical: only for $\mu\to1$ does the quantum Cramér–Rao bound attain its minimum.

Engineered Reservoirs and Quantum Heat Engines:

STS function as customizable, phase-sensitive quantum reservoirs for studying generalized quantum thermodynamics, entropy production, fluctuation theorems, and information erasure at the quantum limit (Xu, 4 Feb 2026).

Squeezed thermal states serve as a versatile and analytically tractable class of continuous-variable quantum states. Their algebraic structure, entropic properties, and phase-space signatures—combined with experimental accessibility and broad resonance across theoretical and applied domains—establish STS as a key resource for quantum technologies, nonequilibrium thermodynamics, and fundamental studies of decoherence and nonclassicality.