Two-Mode Squeezed Thermal States

• Two-mode squeezed thermal states are bipartite continuous-variable Gaussian states formed by applying a two-mode squeezing transformation to independent thermal states, combining quantum correlations with thermal noise.
• They are central to continuous-variable quantum protocols such as teleportation, metrology, and optomechanical reservoir engineering, with properties characterized by covariance matrices and phase-space methods.
• Non-Gaussian modifications, including photon addition and subtraction, provide routes to tune entanglement and fidelity, offering experimental strategies to enhance quantum communication and state reconstruction.

Searching arXiv for the specified paper and closely related work on two-mode squeezed thermal states to ground the article in cited literature. arxiv_search query="Two-Mode Squeezed Thermal States (Hu et al., 2012) photon-added two-mode squeezed thermal state teleportation covariance matrix" max_results=10 Two-mode squeezed thermal states (TMSTSs), also denoted TMST or two-mode squeezed thermal states in the literature, are bipartite continuous-variable Gaussian states obtained by applying a two-mode squeezing transformation to a product of single-mode thermal states. In the symmetric case, both input modes carry the same mean thermal population, while the squeezing parameter controls the intermode correlations; the pure limit recovers the two-mode squeezed vacuum, and the limit $r=0$ gives two independent thermal modes (Hu et al., 2012). Because they incorporate both EPR-type correlations and input thermal noise, TMSTSs serve as a realistic resource model for continuous-variable teleportation, Gaussian-state metrology, optomechanical reservoir engineering, and spectral-mode quantum state reconstruction (Kumar, 2024).

1. Definition, notation, and limiting cases

The general two-mode squeezed thermal state can be written as

$$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$

where the input single-mode thermal state is

$$\hat\rho_T(n)=\frac{1}{n+1}\exp\!\Bigl[-\ln\!\Bigl(\frac{n+1}{n}\Bigr)\hat a^\dagger\hat a\Bigr],\qquad n\ge 0,$$

and $r,\phi$ are the squeezing strength and phase, respectively (Marian et al., 2016). In the symmetric case used extensively in teleportation and non-Gaussian-state studies,

$$\rho_{\rm TMSTS}(r,\bar n)=S(r)\,\rho_{th,1}\otimes\rho_{th,2}\,S^\dagger(r),$$

with

$$S(r)=\exp\!\bigl[r\,(a^\dagger b^\dagger-a\,b)\bigr], \qquad \rho_{th} =\sum_{k=0}^\infty \frac{\bar n^k}{(\bar n+1)^{k+1}}\,|k\rangle\langle k|.$$

The literature also uses the phase-convention

$$S_2(r)=\exp\!\bigl[r\,(a_1a_2-a_1^\dagger a_2^\dagger)\bigr]$$

for the same class of states (Kumar, 2024).

Several limits organize the subject. When $n_{1}=n_{2}=0$, the state reduces to the two-mode squeezed vacuum (Agasti, 2024). When $r=0$, it reduces to the input product of thermal states (Agasti, 2024). In the special symmetric thermal case $n_{1}=n_{2}=\bar n$, the state is zero-mean and fully specified by a $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$0 covariance matrix, which makes it part of the standard two-mode Gaussian-state manifold studied in metrology and information geometry (Marian et al., 2016).

A recurrent notational distinction is between the broad class of squeezed thermal states with possibly unequal thermal inputs $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$1 and the symmetric TMSTS/TMST family with $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$2. This distinction matters for covariance asymmetry, entanglement thresholds, and the angle of maximally squeezed hybrid quadratures in filtered settings (Agasti, 2024).

2. Gaussian phase-space structure

In quadrature form, a zero-mean two-mode squeezed thermal state has covariance matrix

$$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$3

with

$$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$4

$$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$5

$$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$6

For $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$7, this becomes the standard-form matrix with opposite signs in the $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$8-quadrature correlations (Marian et al., 2016).

For the symmetric TMST parameterization used in teleportation studies, defining $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$9, the covariance matrix in the ordering $$\hat\rho_T(n)=\frac{1}{n+1}\exp\!\Bigl[-\ln\!\Bigl(\frac{n+1}{n}\Bigr)\hat a^\dagger\hat a\Bigr],\qquad n\ge 0,$$0 is

$$\hat\rho_T(n)=\frac{1}{n+1}\exp\!\Bigl[-\ln\!\Bigl(\frac{n+1}{n}\Bigr)\hat a^\dagger\hat a\Bigr],\qquad n\ge 0,$$1

and the Wigner function is

$$\hat\rho_T(n)=\frac{1}{n+1}\exp\!\Bigl[-\ln\!\Bigl(\frac{n+1}{n}\Bigr)\hat a^\dagger\hat a\Bigr],\qquad n\ge 0,$$2

The corresponding characteristic function is

$$\hat\rho_T(n)=\frac{1}{n+1}\exp\!\Bigl[-\ln\!\Bigl(\frac{n+1}{n}\Bigr)\hat a^\dagger\hat a\Bigr],\qquad n\ge 0,$$3

These formulas make explicit that the Gaussian TMST is fully encoded by second moments (Kumar, 2024).

An alternative derivation starts from the Glauber-Sudarshan representation of the thermal input,

$$\hat\rho_T(n)=\frac{1}{n+1}\exp\!\Bigl[-\ln\!\Bigl(\frac{n+1}{n}\Bigr)\hat a^\dagger\hat a\Bigr],\qquad n\ge 0,$$4

which leads, after two-mode squeezing, to a two-mode $$\hat\rho_T(n)=\frac{1}{n+1}\exp\!\Bigl[-\ln\!\Bigl(\frac{n+1}{n}\Bigr)\hat a^\dagger\hat a\Bigr],\qquad n\ge 0,$$5-function

$$\hat\rho_T(n)=\frac{1}{n+1}\exp\!\Bigl[-\ln\!\Bigl(\frac{n+1}{n}\Bigr)\hat a^\dagger\hat a\Bigr],\qquad n\ge 0,$$6

This form is especially useful in normal-order calculations and in deriving closed expressions for non-Gaussian deformations of the state (Hu et al., 2012).

The Gaussian structure also fixes collective-quadrature squeezing. For symmetric mechanical realizations, the collective quadratures

$$\hat\rho_T(n)=\frac{1}{n+1}\exp\!\Bigl[-\ln\!\Bigl(\frac{n+1}{n}\Bigr)\hat a^\dagger\hat a\Bigr],\qquad n\ge 0,$$7

have stationary variances

$$\hat\rho_T(n)=\frac{1}{n+1}\exp\!\Bigl[-\ln\!\Bigl(\frac{n+1}{n}\Bigr)\hat a^\dagger\hat a\Bigr],\qquad n\ge 0,$$8

showing the familiar coexistence of squeezed and anti-squeezed directions (Pontin et al., 2015).

3. Non-Gaussian deformations: photon addition, subtraction, and catalysis

A central non-Gaussian extension is the photon-added TMSTS,

$$\hat\rho_T(n)=\frac{1}{n+1}\exp\!\Bigl[-\ln\!\Bigl(\frac{n+1}{n}\Bigr)\hat a^\dagger\hat a\Bigr],\qquad n\ge 0,$$9

where

$r,\phi$0

is the normalization factor (Hu et al., 2012). Defining

$r,\phi$1

and, for $r,\phi$2,

$r,\phi$3

one obtains the closed form

$r,\phi$4

so the normalization is expressed through a Jacobi polynomial (Hu et al., 2012).

The Wigner function of the photon-added state factorizes as

$r,\phi$5

where $r,\phi$6 is the Gaussian background of the underlying TMSTS and $r,\phi$7 is a non-Gaussian multiplier written in terms of two-variable Hermite polynomials. For single-photon addition $r,\phi$8,

$r,\phi$9

which gives an analytic certificate of Wigner negativity (Hu et al., 2012). This sharply separates the Gaussian TMSTS, whose Wigner function is Gaussian, from its photon-added descendants, which are non-Gaussian and can be Wigner-negative.

A broader experimentally realistic scheme applies three operations to a TMST resource: photon subtraction, photon addition, and photon catalysis. Each TMST mode $$\rho_{\rm TMSTS}(r,\bar n)=S(r)\,\rho_{th,1}\otimes\rho_{th,2}\,S^\dagger(r),$$0 is mixed with an ancillary Fock state $$\rho_{\rm TMSTS}(r,\bar n)=S(r)\,\rho_{th,1}\otimes\rho_{th,2}\,S^\dagger(r),$$1 at a beam splitter of transmissivity $$\rho_{\rm TMSTS}(r,\bar n)=S(r)\,\rho_{th,1}\otimes\rho_{th,2}\,S^\dagger(r),$$2, and one conditions on detecting $$\rho_{\rm TMSTS}(r,\bar n)=S(r)\,\rho_{th,1}\otimes\rho_{th,2}\,S^\dagger(r),$$3 photons in the auxiliary output. In the idealized limit $$\rho_{\rm TMSTS}(r,\bar n)=S(r)\,\rho_{th,1}\otimes\rho_{th,2}\,S^\dagger(r),$$4, photon subtraction corresponds to

$$\rho_{\rm TMSTS}(r,\bar n)=S(r)\,\rho_{th,1}\otimes\rho_{th,2}\,S^\dagger(r),$$5

photon addition to

$$\rho_{\rm TMSTS}(r,\bar n)=S(r)\,\rho_{th,1}\otimes\rho_{th,2}\,S^\dagger(r),$$6

and photon catalysis tends to an operator proportional to $$\rho_{\rm TMSTS}(r,\bar n)=S(r)\,\rho_{th,1}\otimes\rho_{th,2}\,S^\dagger(r),$$7 (Kumar, 2024).

These operations are not interchangeable. Photon addition and subtraction both increase non-Gaussianity, but their entanglement and teleportation effects differ quantitatively; photon catalysis introduces a transmissivity-dependent mixture of identity and photon-number-dependent damping rather than a simple ladder-operator action (Kumar, 2024).

4. Entanglement, EPR squeezing, and non-locality

For Gaussian TMSTSs, entanglement is naturally formulated through partial transpose at the covariance-matrix level. In the symmetric unfiltered case,

$$\rho_{\rm TMSTS}(r,\bar n)=S(r)\,\rho_{th,1}\otimes\rho_{th,2}\,S^\dagger(r),$$8

where $$\rho_{\rm TMSTS}(r,\bar n)=S(r)\,\rho_{th,1}\otimes\rho_{th,2}\,S^\dagger(r),$$9 is the smallest symplectic eigenvalue of the partially transposed covariance matrix, and the logarithmic negativity is

$$S(r)=\exp\!\bigl[r\,(a^\dagger b^\dagger-a\,b)\bigr], \qquad \rho_{th} =\sum_{k=0}^\infty \frac{\bar n^k}{(\bar n+1)^{k+1}}\,|k\rangle\langle k|.$$0

Thus $$S(r)=\exp\!\bigl[r\,(a^\dagger b^\dagger-a\,b)\bigr], \qquad \rho_{th} =\sum_{k=0}^\infty \frac{\bar n^k}{(\bar n+1)^{k+1}}\,|k\rangle\langle k|.$$1 is necessary and sufficient for two-mode Gaussian entanglement in this setting (Agasti, 2024).

Equivalent criteria appear in different normalizations. In the Gaussian optomechanical treatment, the Duan condition reads

$$S(r)=\exp\!\bigl[r\,(a^\dagger b^\dagger-a\,b)\bigr], \qquad \rho_{th} =\sum_{k=0}^\infty \frac{\bar n^k}{(\bar n+1)^{k+1}}\,|k\rangle\langle k|.$$2

while in spectral homodyne reconstruction the Simon PPT test is implemented through the smallest partially transposed symplectic eigenvalue $$S(r)=\exp\!\bigl[r\,(a^\dagger b^\dagger-a\,b)\bigr], \qquad \rho_{th} =\sum_{k=0}^\infty \frac{\bar n^k}{(\bar n+1)^{k+1}}\,|k\rangle\langle k|.$$3, with entanglement present iff $$S(r)=\exp\!\bigl[r\,(a^\dagger b^\dagger-a\,b)\bigr], \qquad \rho_{th} =\sum_{k=0}^\infty \frac{\bar n^k}{(\bar n+1)^{k+1}}\,|k\rangle\langle k|.$$4 in that quadrature convention (Woolley et al., 2014). In the electromechanical parametrically driven setting, the analytic entanglement condition is

$$S(r)=\exp\!\bigl[r\,(a^\dagger b^\dagger-a\,b)\bigr], \qquad \rho_{th} =\sum_{k=0}^\infty \frac{\bar n^k}{(\bar n+1)^{k+1}}\,|k\rangle\langle k|.$$5

and for equal thermal populations $$S(r)=\exp\!\bigl[r\,(a^\dagger b^\dagger-a\,b)\bigr], \qquad \rho_{th} =\sum_{k=0}^\infty \frac{\bar n^k}{(\bar n+1)^{k+1}}\,|k\rangle\langle k|.$$6 it is reported as

$$S(r)=\exp\!\bigl[r\,(a^\dagger b^\dagger-a\,b)\bigr], \qquad \rho_{th} =\sum_{k=0}^\infty \frac{\bar n^k}{(\bar n+1)^{k+1}}\,|k\rangle\langle k|.$$7

(Mahboob et al., 2014).

For non-Gaussian photon-added TMSTSs, one convenient sufficient inseparability test is the Shchukin-Vogel quantity

$$S(r)=\exp\!\bigl[r\,(a^\dagger b^\dagger-a\,b)\bigr], \qquad \rho_{th} =\sum_{k=0}^\infty \frac{\bar n^k}{(\bar n+1)^{k+1}}\,|k\rangle\langle k|.$$8

with inseparability certified when $$S(r)=\exp\!\bigl[r\,(a^\dagger b^\dagger-a\,b)\bigr], \qquad \rho_{th} =\sum_{k=0}^\infty \frac{\bar n^k}{(\bar n+1)^{k+1}}\,|k\rangle\langle k|.$$9. Numerical plots show that for realistic $$S_2(r)=\exp\!\bigl[r\,(a_1a_2-a_1^\dagger a_2^\dagger)\bigr]$$0, this threshold is reached at lower squeezing $$S_2(r)=\exp\!\bigl[r\,(a_1a_2-a_1^\dagger a_2^\dagger)\bigr]$$1 for photon addition than for photon subtraction; for $$S_2(r)=\exp\!\bigl[r\,(a_1a_2-a_1^\dagger a_2^\dagger)\bigr]$$2, the quoted thresholds are $$S_2(r)=\exp\!\bigl[r\,(a_1a_2-a_1^\dagger a_2^\dagger)\bigr]$$3 for single-photon addition and $$S_2(r)=\exp\!\bigl[r\,(a_1a_2-a_1^\dagger a_2^\dagger)\bigr]$$4 for photon subtraction (Hu et al., 2012).

Two points qualify the common expectation that “more squeezing always helps.” First, Bell-CHSH non-locality is stricter than entanglement. Using displaced parity operators, the phase-space CHSH parameter satisfies $$S_2(r)=\exp\!\bigl[r\,(a_1a_2-a_1^\dagger a_2^\dagger)\bigr]$$5 under local realism, and for a pure two-mode squeezed vacuum one has

$$S_2(r)=\exp\!\bigl[r\,(a_1a_2-a_1^\dagger a_2^\dagger)\bigr]$$6

For mixed or filtered states, the Bell-violating region lies strictly inside the entangled region (Agasti, 2024). Second, with non-identical filters the filtered logarithmic negativity becomes bell-shaped in $$S_2(r)=\exp\!\bigl[r\,(a_1a_2-a_1^\dagger a_2^\dagger)\bigr]$$7: it rises, reaches a maximum near

$$S_2(r)=\exp\!\bigl[r\,(a_1a_2-a_1^\dagger a_2^\dagger)\bigr]$$8

and vanishes at

$$S_2(r)=\exp\!\bigl[r\,(a_1a_2-a_1^\dagger a_2^\dagger)\bigr]$$9

where $n_{1}=n_{2}=0$0 is the filter-overlap factor. This means that arbitrarily large input squeezing can be counterproductive once spectral overlap is reduced (Agasti, 2024).

5. Teleportation resource and statistical geometry

In the Braunstein-Kimble protocol, a TMST or its non-Gaussian variants may be used as the shared entangled resource. For a coherent-state input, the Gaussian benchmark fidelity is

$n_{1}=n_{2}=0$1

which makes the degradation from thermal noise explicit (Kumar, 2024). In the photon-added analysis, the coherent-state $n_{1}=n_{2}=0$2 fidelity is written as $n_{1}=n_{2}=0$3, with special cases

$n_{1}=n_{2}=0$4

$n_{1}=n_{2}=0$5

The numerical comparison shows that symmetric photon addition $n_{1}=n_{2}=0$6 raises fidelity for $n_{1}=n_{2}=0$7 above a modest threshold, whereas non-symmetric addition $n_{1}=n_{2}=0$8 typically degrades fidelity (Hu et al., 2012).

When success probability is incorporated, the comparison among non-Gaussian operations changes. Defining the fidelity enhancement $n_{1}=n_{2}=0$9, the success probability $r=0$0, and the figure of merit

$r=0$1

numerical optimization for teleporting a coherent state with $r=0$2 yields the following quoted optima: single-photon catalysis is optimal at low squeezing with $r=0$3, $r=0$4, $r=0$5, $r=0$6, $r=0$7, and $r=0$8; single-photon subtraction is optimal at intermediate squeezing with $r=0$9, $n_{1}=n_{2}=\bar n$0, $n_{1}=n_{2}=\bar n$1, $n_{1}=n_{2}=\bar n$2, $n_{1}=n_{2}=\bar n$3, and $n_{1}=n_{2}=\bar n$4 (Kumar, 2024).

Beyond communication, TMSTSs define a four-parameter Gaussian-state manifold with parameters $n_{1}=n_{2}=\bar n$5. Marian and Marian showed that the quantum Fisher information metric on this manifold is diagonal. Writing $n_{1}=n_{2}=\bar n$6,

$n_{1}=n_{2}=\bar n$7

with

$n_{1}=n_{2}=\bar n$8

$n_{1}=n_{2}=\bar n$9

The diagonal form implies orthogonality of the parameters in the local metric sense, and the quantum Cramér-Rao bound becomes

$$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$00

for $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$01 independent copies (Marian et al., 2016).

6. Preparation, reconstruction, filtering, and dynamical settings

TMSTSs appear both as directly generated optical states and as effective steady states in engineered dissipative systems. In optical experiments, symmetric two-mode squeezed thermal states of spectral sideband modes $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$02 can be reconstructed with a single homodyne detector and active cavity stabilization. The measured covariance matrix is assembled in the symmetric/antisymmetric spectral-mode basis and then transformed back to the original sideband basis without iterative maximum-likelihood reconstruction, because Gaussianity reduces the problem to moment inversion (Cialdi et al., 2015). The same experiment reports successful tests on states ranging from uncorrelated coherent states to entangled states, with reconstructed examples including a pure squeezed state with $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$03 and $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$04 bit (Cialdi et al., 2015).

In optomechanics, TMSTS-like states arise as targets of reservoir engineering. A three-mode setup with two mechanical oscillators coupled to a single cavity mode admits an effective linearized Hamiltonian

$$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$05

which is rewritten in terms of Bogoliubov operators

$$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$06

with $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$07. Cooling the sum mode $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$08 together with coherent coupling from $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$09 yields a highly pure two-mode squeezed mechanical steady state; in the resolved-sideband regime, cooperativities $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$10 are quoted as sufficient for $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$11 dB of two-mode squeezing, $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$12, purity $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$13, and teleportation fidelity $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$14 (Woolley et al., 2014).

Closely related dissipation-driven proposals show that finite-temperature steady states can still be written formally as

$$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$15

and that ground-state precooling is not necessary for two-mode squeezing. In the large-cooperativity limit, the maximal tolerable thermal occupation $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$16 can become very large even at modest $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$17; one worked example quotes $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$18 for $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$19, $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$20, $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$21, and $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$22 (Tan et al., 2013).

Parametric modulation provides a complementary route. In a cavity opto-mechanical system, modulation of the optical spring at $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$23 generates the two-mode squeezing Hamiltonian

$$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$24

and pulsed excitation can transiently surpass the stationary $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$25 dB limit. The reported experiment reaches a minimal normalized spectral density $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$26 at $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$27, corresponding to $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$28 dB, and in the pulsed regime the minimal variance falls transiently to $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$29 at $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$30 ms (Pontin et al., 2015). In an electromechanical resonator implementing the analogous parametric-down-conversion Hamiltonian

$$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$31

the measured suppression of joint quadrature noise reaches $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$32 dB below the bare thermal level at $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$33 K, with correlation coefficients approaching unity as the pump is increased (Mahboob et al., 2014).

Filtered and open-system variants further refine the picture. For filtered two-mode squeezed mixed states, entanglement and non-locality are maximized when the two filters are identical, and the filter overlap factor $$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$34 rescales the intermode correlations (Agasti, 2024). Under thermalization dynamics in a Markovian bath,

$$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$35

so the smallest partially transposed symplectic eigenvalue evolves as

$$\hat\rho_{\rm ST}(n_{1},n_{2},r,\phi) =\hat S_{12}(r,\phi)\, \bigl[\hat\rho_{T}(n_{1})\otimes\hat\rho_{T}(n_{2})\bigr]\, \hat S_{12}^\dagger(r,\phi),$$36

The reported conclusions are that increasing the thermal population of the environment enhances the rate of dissipation, stronger interaction slows dissipation in a normalized dimensionless time scale, and identical filters keep entanglement and non-locality at their peak (Agasti, 2024).

Taken together, these results place TMSTSs at the intersection of Gaussian-state theory, non-Gaussian state engineering, CV teleportation, and engineered dissipation. Their defining feature is not merely the coexistence of squeezing and thermal noise, but the fact that this coexistence remains analytically tractable across covariance-matrix methods, entanglement criteria, homodyne tomography, and several experimentally relevant non-Gaussian and filtered generalizations.