---
title: Two-Mode Squeezed LG Beams
url: https://www.emergentmind.com/topics/two-mode-squeezed-laguerre-gaussian-beams
type: topic
---

# Two-Mode Squeezed LG Beams

Searching arXiv for the cited papers and closely related work on two-mode squeezed Laguerre-Gaussian beams.
Two-mode squeezed Laguerre–Gaussian beams are structured optical fields in which continuous-variable squeezing is carried by transverse modes with Laguerre–Gaussian (LG) spatial structure, so that quadrature correlations, orbital-angular-momentum content, and radial-mode structure become coupled within a single paraxial framework. In the literature represented here, the concept appears in three closely connected forms: as vortex-bearing output states obtained by rotating two-mode photon-number squeezed inputs in quadrature space [2104.01967]; as an experimentally observed and theoretically modeled multimode LG decomposition of atom-generated squeezed light that can be projected into effective two-mode LG pairs [1507.02247]; and as a group-theoretic beam-engineering construction in which continuous-variable two-mode squeezed states are mapped directly onto scalar paraxial beams expanded in an LG basis [2509.12426]. Across these formulations, the central theme is that squeezing redistributes phase-space weight and induces correlations among transverse degrees of freedom while preserving the paraxial modal structure characteristic of LG beams.

## 1. Conceptual definition and modal setting

In paraxial optics, LG modes form a complete orthonormal basis for cylindrically symmetric beams. A standard expression is
\[
\mathrm{LG}_{p}^{l}(r,\phi,z) = \frac{C_{pl}}{w(z)}\left(\frac{\sqrt{2}\,r}{w(z)}\right)^{|l|} L_{p}^{|l|}\!\left(\frac{2r^2}{w^2(z)}\right)\, \exp\!\left(-\frac{r^2}{w^2(z)}\right)\, \exp\!\left(i l\,\phi\right)\, \exp\!\left[i k\frac{r^2}{2R(z)} - i\,(2p+|l|+1)\,\zeta(z)\right],
\]
where \(C_{pl}=\sqrt{\frac{2\,p!}{\pi\,(p+|l|)!}}\), \(w(z)=w_0\sqrt{1+(z/z_R)^2}\), \(R(z)=z\,[1+(z_R/z)^2]\), \(\zeta(z)=\arctan(z/z_R)\), and \(z_R=\pi w_0^2/\lambda\) [1507.02247]. In quantum-optical notation, the positive-frequency field operator may be decomposed as
\[
\hat{E}^{(+)}(\mathbf{r},z) = \sum_{p,l} u_{p,l}(\mathbf{r},z)\, \hat{a}_{p,l}(z),
\]
with \(u_{p,l}\equiv \mathrm{LG}_{p}^{l}\) and bosonic mode operators \(\hat a_{p,l}\) [1507.02247].

The “two-mode” descriptor is used in two technically distinct but related senses. In spontaneous parametric down-conversion, it refers to two correlated bosonic modes \(a\) and \(b\), whose squeezing is generated by
\[
S_2(\xi) = \exp\big(\xi \, \hat a \hat b - \xi^* \, \hat a^\dagger \hat b^\dagger\big), \qquad \xi = r e^{i\theta},
\]
producing the two-mode squeezed vacuum
\[
|\Psi_{2\mathrm{MS}}\rangle = \sqrt{1-\lambda^2} \sum_{n=0}^{\infty} \lambda^n \, |n\rangle_a \, |n\rangle_b, \qquad \lambda = \tanh r
\]
[2104.01967]. In the spatial beam-engineering framework, the two modes are the two transverse degrees of freedom \(x\) and \(y\), represented by canonical quadratures \((q_x,p_x)\) and \((q_y,p_y)\), with two-mode squeezing implemented through an \(\mathrm{su}(1,1)\) algebraic action [2509.12426].

This dual usage is not contradictory. It indicates that two-mode squeezing may be represented either in Hilbert-space mode labels or in coupled transverse canonical variables. A plausible implication is that the structured-light interpretation and the nonclassical-state interpretation are different projections of the same underlying covariance structure.

## 2. Two-mode squeezing formalisms and quadrature correlations

For the ideal two-mode squeezed vacuum, the covariance matrix in canonical ordering \((X_a,P_a,X_b,P_b)\) is
\[
V = \frac{1}{2}\begin{pmatrix}
\cosh(2r)\, \mathbb{1}_2 & \sinh(2r)\, \sigma_z \\
\sinh(2r)\, \sigma_z & \cosh(2r)\, \mathbb{1}_2
\end{pmatrix}, \qquad \sigma_z=\mathrm{diag}(1,-1),
\]
which yields the EPR-like correlations
\[
\mathrm{Var}(\hat X_a - \hat X_b) = \frac{e^{-2r}}{2}, \qquad \mathrm{Var}(\hat P_a + \hat P_b) = \frac{e^{-2r}}{2}
\]
[2104.01967]. The associated Wigner function is Gaussian,
\[
W(\boldsymbol{\xi}) = \frac{1}{\pi^2 \sqrt{\det V}} \exp\!\big(-\boldsymbol{\xi}^T V^{-1}\boldsymbol{\xi}\big), \qquad \boldsymbol{\xi}=(X_a,P_a,X_b,P_b)^T
\]
[2104.01967].

In the beam-theoretic formulation, the \(\mathrm{su}(1,1)\) generators for the two transverse degrees of freedom are
\[
K_0 = \frac{1}{4}(p_x^2 + q_x^2 + p_y^2 + q_y^2), \qquad
K_1 = \frac{1}{2}(q_x q_y - p_x p_y), \qquad
K_2 = -\frac{1}{2}(q_x p_y + q_y p_x),
\]
with commutators
\[
[K_0, K_\pm] = \pm K_\pm, \qquad [K_+, K_-] = -2 K_0
\]
[2509.12426]. Under the squeezing action generated by \(K_2\),
\[
e^{i \zeta_2 K_2} q_x e^{-i \zeta_2 K_2} = q_x \cosh(\zeta_2/2) - q_y \sinh(\zeta_2/2),
\]
\[
e^{i \zeta_2 K_2} q_y e^{-i \zeta_2 K_2} = q_y \cosh(\zeta_2/2) - q_x \sinh(\zeta_2/2),
\]
\[
e^{i \zeta_2 K_2} p_x e^{-i \zeta_2 K_2} = p_x \cosh(\zeta_2/2) + p_y \sinh(\zeta_2/2),
\]
\[
e^{i \zeta_2 K_2} p_y e^{-i \zeta_2 K_2} = p_y \cosh(\zeta_2/2) + p_x \sinh(\zeta_2/2)
\]
[2509.12426].

These exact operator relations show that two-mode squeezing does not merely rescale a single quadrature pair. It mixes the transverse degrees of freedom hyperbolically, thereby creating joint spatial and spectral correlations. In the LG basis, those correlations are represented as redistribution among radial orders at fixed azimuthal index \(\ell\) [2509.12426].

## 3. Generation routes

Two main physical routes are described in the cited work.

The first route begins from two-mode photon-number squeezed states produced by spontaneous parametric down-conversion. The working finite-\(N\) input used in the vortex-generation scheme is
\[
|\psi\rangle = \frac{D}{\cosh r} \sum_{j=0}^{N/2} (\tanh r)^j \, |j,j\rangle,
\]
with normalization factor \(D\) and squeezing parameter \(r\) [2104.01967]. This state exhibits sub-Poissonian number statistics and a Gaussian quadrature-space profile without vortices [2104.01967]. The quadrature-space wavefunction is built from Hermite–Gaussian functions,
\[
\psi_{n_x,n_y}(x,y)= \sqrt{\frac{1}{\pi\,2^{n_x+n_y}\,n_x!\,n_y!}}\,H_{n_x}(x)\,H_{n_y}(y)\,e^{-(x^2+y^2)/2},
\]
so that
\[
\psi(x,y)=\langle x,y|\psi\rangle = \frac{D}{\cosh r}\sum_{j=0}^{N/2}(\tanh r)^j\,H_j(x)\,H_j(y)\,e^{-(x^2+y^2)/2}
\]
[2104.01967].

The second route is polarization self-rotation in dense rubidium vapor. In that experiment, a squeezed vacuum field is generated via polarization self-rotation interaction of an ensemble of Rb atoms with a strong near-resonant linearly polarized laser field. The orthogonal polarization after the cell contains a squeezed vacuum state measured by balanced homodyne detection using the pump as the local oscillator, with typical results of \(1.9\pm0.2\) dB of squeezing and \(11\) dB of anti-squeezing at \(1\) MHz [1507.02247]. The experiment found that the spatial structure is multimode rather than single-mode and can be decomposed into LG modes [1507.02247].

The PSR theory derives an inhomogeneous paraxial propagation equation,
\[
\left(\frac{\partial}{\partial z} - \frac{i}{2k}\nabla_\perp^2\right)\tilde{E}(\mathbf{r},z) = \frac{i k}{2\epsilon_0}\,\frac{N}{V}\sum_{m,n} \mu_{mn}\,\tilde{\rho}_{mn}(\mathbf{r},z),
\]
and, in a weak-scattering approximation, obtains a superposition of radial LG modes,
\[
\Omega_{x,y}(\mathbf{r},z)= \Omega_{x,y}^{\mathrm{homo}}(\mathbf{r},z)+\sum_{p}\,u_{0,p}(\mathbf{r},z)\,\int_{\mathrm{cell}}dz'\,c_{0,p}(z')
\]
[1507.02247]. Selection rules favor \(l=0\) because the pump is cylindrically symmetric and the source term carries no \(\phi\) dependence [1507.02247].

A more abstract but exact construction is given by the scalar-beam mapping of continuous-variable states. There, two-mode squeezed LG beams are defined directly by the expansion
\[
\Psi_{p,\ell,\zeta}(\rho, \theta, z) = \sum_{j=0}^\infty c_j(p, |\ell|, \alpha, \beta) \Psi_{j,\ell}(\rho, \theta, z), \qquad \zeta=\alpha e^{i\beta},
\]
with coefficients
\[
c_j(p, \ell, \alpha, \beta) = (-1)^j
\sqrt{ C(|\ell| + p - 1, p) C(|\ell| + j - 1, j) }
\, e^{i \beta (p - j)} (\mathrm{sech}\,\alpha)^{|\ell| + 1} (\tanh \alpha)^{p + j}
\, {}_2F_1( -j, -p; |\ell| + 1; -\mathrm{csch}^2 \alpha )
\]
[2509.12426].

Taken together, these routes show three experimentally or formally relevant regimes: direct quantum-state generation in SPDC, atom-generated multimode squeezing in PSR, and mode-engineered optical analogs based on \(\mathrm{su}(1,1)\) beam synthesis.

## 4. Hermite–Gaussian to Laguerre–Gaussian conversion and vortex formation

A defining mechanism in the SPDC-based construction is a quadrature-space “\(\phi\) converter,” implemented by
\[
\hat C(\phi) = \exp\!\big(i2\phi\, [\hat a^\dagger \hat b + \hat b^\dagger \hat a]\big)
\]
[2104.01967]. In the Heisenberg picture,
\[
\hat a \to \hat a\cos\phi - i\,\hat b\sin\phi, \qquad
\hat b \to \hat b\cos\phi - i\,\hat a\sin\phi,
\]
so that at \(\phi=\pi/4\),
\[
\hat a \to \frac{\hat a - i\hat b}{\sqrt{2}}, \qquad
\hat b \to \frac{\hat b - i\hat a}{\sqrt{2}}
\]
[2104.01967]. In quadrature space this corresponds to the rotation
\[
\begin{pmatrix} X' \\ Y' \end{pmatrix}
=
\begin{pmatrix}
\cos\phi & -\sin\phi \\
\sin\phi & \cos\phi
\end{pmatrix}
\begin{pmatrix} X \\ Y \end{pmatrix}
\]
[2104.01967].

The paper invokes the standard astigmatic HG\(\to\)LG conversion result that Hermite–Gaussian mode products are converted into single Laguerre–Gaussian modes by imparting a \(\pi/2\) phase shift in the mode profile, and states that choosing \(\phi=\pi/4\) implements exactly that phase condition [2104.01967]. After rotation, the state becomes a coherent superposition of Fock states \(|j-(l-k),\,j+(l-k)\rangle\) with \(\phi\)-dependent coefficients [2104.01967]. In quadrature space, the relevant Hermite–Gaussian pairs \(H_{j-m}(x)H_{j+m}(y)\) map to LG terms \(L_{j-m}^{\,2m}(r^2)\) under the Gaussian envelope, yielding
\[
\psi_{\mathrm{LG}}(x,y) = \frac{1}{\cosh r}\sum_{j=0}^{N/2}\sum_{m=0}^{j} A^{r,N}_{j}\, L_{j-m}^{\,2m}(x^2+y^2)\, e^{-(x^2+y^2)/2}
\]
[2104.01967].

The main numerical finding is that, for sufficiently small squeezing parameter \(r\), the transformed quadrature-space field develops \(N\) phase singularities for an input state truncated at total photon number \(N\) [2104.01967]. In the reported numerics, the dominant component at small \(r\) corresponds to \(m=N/2\), hence \(\ell=N\), so the observed number of vortices equals the input photon number \(N\) [2104.01967]. The article’s own interpretation is that reduced uncertainty in photon number in Fock space is mapped to reduced uncertainty in vortex number in quadrature space [2104.01967].

This mapping should not be confused with a generic statement that any squeezed LG beam must contain \(N\) vortices. The claim is specific to the truncated input family, the \(\phi=\pi/4\) converter, and the small-\(r\) regime treated numerically in that work [2104.01967]. By contrast, the general \(\mathrm{su}(1,1)\) two-mode squeezed LG construction preserves \(\ell\) while redistributing radial order \(j\) [2509.12426].

## 5. Phase-space structure, Wigner analysis, and inseparability

The cited literature treats two-mode squeezed LG beams as phase-space objects as much as spatial modes. In the scalar-beam mapping, the optical Wigner distribution is defined as
\[
W(x, y, k_x, k_y) = \frac{1}{(2\pi)^2} \int\!\!\int du\, dv\, e^{-i (k_x u + k_y v)}
\Psi^*(x - u/2, y - v/2)\Psi(x + u/2, y + v/2)
\]
[2509.12426]. Squeezing deforms this distribution hyperbolically along coupled quadratures while preserving area, and for LG beams at fixed \(\ell\) yields tilted or elliptical contours in spatial and spectral marginals, signaling intermodal correlations analogous to continuous-variable entanglement [2509.12426].

For non-Gaussian states, the same work states that Wigner marginals can exhibit negativity due to phase-space interference from nodal structure; in that framework, such negativity acts as a witness of classical continuous-variable inseparability rather than quantum nonclassicality [2509.12426]. This distinction is important because the paper explicitly constructs classical analogs of continuous-variable quantum states, not quantum light sources [2509.12426].

A Duan–Simon-type criterion is also presented. Defining
\[
\hat u = q_x + q_y, \qquad \hat v = p_x - p_y,
\]
the separability bound for Gaussian states is
\[
\langle(\Delta \hat u)^2\rangle + \langle(\Delta \hat v)^2\rangle \ge 2
\]
[2509.12426]. The paper states that for the Gaussian case \((p,\ell)=(0,0)\), any \(\alpha>0\) yields a violation of the bound in the optical analog [2509.12426].

The PSR work formulates the same issue in a quantum-optical multimode setting. Its effective Hamiltonian can include cross-mode bilinear couplings,
\[
\hat{H}_{\mathrm{eff}} \sim i\hbar \sum_{(p,l),(p',l')} \big( \kappa_{(p,l),(p',l')}\,\hat{a}_{p,l}\hat{a}_{p',l'}\,e^{-i\phi} - \kappa^*_{(p,l),(p',l')}\,\hat{a}^{\dagger}_{p,l}\hat{a}^{\dagger}_{p',l'}\,e^{i\phi} \big),
\]
leading to a two-mode squeezer
\[
\hat{S}_2(r,\phi)=\exp\!\big[r\,e^{i\phi}\,\hat{a}_{m}\hat{a}_{n}-r\,e^{-i\phi}\,\hat{a}^{\dagger}_{m}\hat{a}^{\dagger}_{n}\big]
\]
with joint quadrature variances
\[
V(\hat{X}_\pm)=e^{-2r},\qquad V(\hat{P}_\mp)=e^{-2r}
\]
for appropriate combinations [1507.02247]. The same paper notes that the reported experiment demonstrates multimode squeezing and LG content, but does not directly extract two-mode entanglement [1507.02247]. This is a useful boundary condition on interpretation: multimode LG squeezing and the availability of a two-mode LG description do not, by themselves, constitute a direct entanglement measurement.

## 6. Beam parameters, propagation, and numerical-aperture constraints

For the squeezed Gaussian seed \((p,\ell)=(0,0)\), the \(\mathrm{su}(1,1)\) mapping yields a closed-form propagated field,
\[
\Psi_{0,0,\zeta}(\rho, \theta, z) = \frac{1}{\sqrt{\pi}} \frac{\sqrt{2}}{w(z; \alpha, \beta)}
\exp\!\left[ - \frac{i k \rho^2}{2 R(z; \alpha, \beta)} \right]
\exp\!\left[- i \phi(z; \alpha, \beta) \right]
\exp\!\left[ - \frac{\rho^2}{w^2(z; \alpha, \beta)} \right]
\]
with exact modified parameters
\[
w(z; \alpha, \beta) = \left| 1 - e^{-i 2 \phi(z)} e^{i \beta} \coth \alpha \right| \sinh \alpha \cdot w(z),
\]
\[
R(z; \alpha, \beta) =
\frac{ \left| 1 - e^{-i 2 \phi(z)} e^{i \beta} \coth \alpha \right|^2 w^2(z) \cdot R(z)}
{ \left| 1 - e^{-i 2 \phi(z)} e^{i \beta} \coth \alpha \right|^2 w^2(z) + 4 \sin(2 \phi(z) - \beta) \coth \alpha },
\]
\[
\phi(z; \alpha, \beta) = \phi(z) + \arctan\!\left(
\frac{\sin(2 \phi(z) - \beta)\coth \alpha}{\tanh \alpha - \cos(2 \phi(z) - \beta)}
\right)
\]
[2509.12426]. The minimum waist at focus occurs for \(\beta=0\) [2509.12426].

At the same time, the work stresses that finite numerical aperture limits the highest accessible radial order and therefore the achievable squeezing. The relevant cutoff relations are
\[
k_{\rho,\max}^{(A)} = \frac{2\pi}{\lambda_0\, NA^{(A)}}, \qquad
NA^{(G)} \approx \frac{\lambda_0}{\pi w_0}
\]
[2509.12426]. The operational bound is expressed through truncation of the LG expansion at \(j\le j_{\max}\), with \(\alpha_{\max}\) set by the requirement
\[
\sum_{j=0}^{j_{\max}} |c_j(p, |\ell|, \alpha_{\max}, \beta)|^2 \approx 1
\]
[2509.12426]. The paper states that stronger squeezing demands higher numerical aperture because squeezing pushes energy into higher-order LG modes with larger divergence and higher spatial-frequency content [2509.12426]. Numerical propagation of aperture-truncated beams confirms that diffraction remains the ultimate limit [2509.12426].

The PSR study presents a related but experimentally grounded propagation constraint: different LG radial orders acquire different Gouy-phase shifts near the focus, so their squeezing angles become mode dependent [1507.02247]. This is directly relevant to two-mode LG beams because joint-quadrature measurements require those modal squeezing axes to be phase aligned. The reported data showed strong asymmetries and sharp noise changes when scanning an iris along \(z\) near the focus, behavior consistent with Gouy-phase evolution [1507.02247].

## 7. Detection, applications, limitations, and comparison of regimes

In the SPDC-based vortex scheme, the resulting LG vortex states can be characterized by measuring the OAM spectrum via mode projectors or phase-flattening and coupling into single-mode fibers, by interferometric diagnostics of azimuthal phase and topological charge such as fork holograms, by direct camera imaging of intensity and phase to count vortices and verify azimuthal phase winding, and by homodyne tomography in quadrature space [2104.01967]. In the PSR setting, homodyne detection with a shaped or masked local oscillator performs mode-selective projection onto desired transverse profiles, with photocurrent variance given as a weighted sum of modal quadrature variances through overlap integrals \(O_p\) [1507.02247].

The following comparison organizes the roles played by LG structure in the three cited formulations.

| Regime | Principal LG feature | Key technical point |
|---|---|---|
| SPDC plus \(\phi\) converter [2104.01967] | Vortex formation in quadrature space | Small-\(r\) numerics give \(N\) vortices for input truncated at photon number \(N\) |
| PSR in Rb vapor [1507.02247] | Multimode LG decomposition dominated by \(l=0\) | Spatial masks and LO shaping reveal hidden LG squeezing and enable effective two-mode mode pairs |
| \(\mathrm{su}(1,1)\) beam mapping [2509.12426] | Exact two-mode squeezed LG expansion | Squeezing redistributes radial order at fixed \(\ell\), with finite NA limiting achievable squeezing |

Applications highlighted across the sources include quantum metrology and sensing, quantum cryptography, OAM- and radial-mode–based communications, imaging, and interfaces to atomic memories [2104.01967; 1507.02247; 2509.12426]. The SPDC paper argues that the scheme “has the potential to improve measurement sensitivity beyond the Standard Quantum Limit (SQL \(\propto \sqrt{N}\)),” but it does not provide an explicit Fisher-information or Cramér–Rao analysis and does not claim Heisenberg scaling [2104.01967]. The PSR paper emphasizes that spatial multimode structure is not merely a complication but can be harnessed as a resource by selective projection onto LG components [1507.02247]. The beam-mapping paper frames the corresponding correlations as a classical analog of continuous-variable entanglement and explicitly states that diffraction constraints remain in force [2509.12426].

Several limitations recur. Optical loss and detector inefficiency wash out number correlations and reduce effective squeezing; imperfect HG\(\to\)LG conversion and misalignment lead to mixed \(\ell\)-content and degraded vortex visibility; decoherence and phase noise smear azimuthal phase and obscure singularities [2104.01967]. In PSR, near-resonant absorption and technical noise limit observable squeezing, while mode-dependent Gouy-phase rotations misalign squeezing angles and degrade Gaussian-LO measurements [1507.02247]. In beam-engineered analogs, finite numerical aperture truncates the high-order modal content required for strong squeezing and enforces the diffraction limit [2509.12426].

A common misconception is that LG squeezing necessarily implies single-mode OAM squeezing in a pure \((p,\ell)\) state. The cited work instead shows three more nuanced possibilities: squeezing may be encoded in correlated number-state inputs and revealed as vortex structure after quadrature rotation [2104.01967]; it may be distributed across multiple LG radial orders and hidden by an unshaped local oscillator [1507.02247]; or it may be represented as an exact radial-order superposition at fixed \(\ell\) in a classical analog of continuous-variable two-mode squeezing [2509.12426]. This suggests that “two-mode squeezed Laguerre–Gaussian beams” are best understood not as a single canonical object but as a family of structured squeezed states and beam constructions linked by common transverse-mode algebra, quadrature correlations, and phase-space geometry.

Source: https://www.emergentmind.com/topics/two-mode-squeezed-laguerre-gaussian-beams