## Abstract

Four-wave mixing near resonance in an atomic vapor can produce relative intensity squeezed light suitable for precision measurements beyond the shot-noise limit. We develop an analytic distributed gain/loss model to describe the competition of mixing and absorption through the non-linear medium. Using a novel matrix calculus, we present closed-form expressions for the degree of relative intensity squeezing produced by this system. We use these theoretical results to analyze experimentally measured squeezing from a ^{85}Rb vapor and demonstrate the analytic model’s utility as an experimental diagnostic.

© 2011 Optical Society of America

## 1. Introduction

Relative intensity squeezing by four-wave mixing in an atomic vapor is emerging as a promising technique for performing high-precision measurements beyond the shot-noise limit. First demonstrated by McCormick et al. [1], the technique uses atomic coherences to produce quantum correlated “twin beams”, enabling the shot-noise of one beam to be measured and subtracted from the other to obtain a low-noise differential measurement; for example of a weakly absorbing sample. This scheme was recently shown to reduce the relative intensity noise by 9.2±0.5dB below the shot-noise limit [2], and noise reduction has been observed in both the low Fourier frequency [3] and multi-mode imaging [4, 5] domains.

Furthermore, as one of the twin beams is near-resonant with the atoms, this squeezing technique has promising applications in quantum information processing [6, 7]. However, absorption near resonance degrades the quantum correlations. Both mixing gain and absorption losses occur simultaneously as the beams propagate through the vapor, and are therefore competing processes.

Earlier theoretical investigations of this system have applied numerical methods [3] and the Heisenberg-Langevin formalism [8] to predict the resulting degree of squeezing. The numerical model demonstrated excellent agreement with experimental results, but it can be difficult to gain insight into the competing processes from numerical calculations. The Heisenberg-Langevin model provided a microscopic description of a specific four-wave mixing configuration in a cold atomic gas, which accurately predicted the resulting gain profiles. However, calculation of the predicted squeezing required complex matrix integrals and no comparison to experimentally measured squeezing was presented.

In this work, we present a very general approach for determining the squeezing produced by a four-wave mixing system, and develop a matrix-based analysis method to include arbitrarily many injected vacuum modes. Considering special cases, simple closed-form expressions are easily obtained. Finally, we present experimentally measured squeezing from four-wave mixing in a rubidium-85 vapor, and demonstrate how the model can be used as a diagnostic tool to determine the limiting technical factors.

## 2. Relative intensity squeezing

The “double-Λ” four-wave mixing scheme introduced by McCormick et al [1] uses a high-intensity “pump” beam to drive a cycle of four off-resonant transitions in a hot rubidium vapor, causing the emission of correlated “probe” and “conjugate” photons (Fig. 1A). The probe transition is stimulated by a seed laser incident at an angle *θ* to the pump, resulting in the spontaneous emission of the conjugate on the opposite side of the pump beam (Fig. 1B). The beam powers are measured individually and subtracted to obtain the relative intensity noise as measured on a spectrum analyzer (S.A.).

Labelling the Fock-space annihilation operators of the probe, conjugate and pump by *â*, *b̂* and *ĉ* respectively and the interaction strength by *ξ*, the interaction picture Hamiltonian is

*ψ*〉 and the substitution

_{c}*ĉ*→

*ψ*can be made:

_{c}*τ*is

*â*and

*b̂*, where

*s*is the “squeezing parameter” [9]. The four-wave mixing system therefore produces a two-mode squeezed state, reducing amplitude difference noise at the expense of increasing phase difference noise [5].

The phase of *s* results in a rotation of the (arbitrary) measurement quadratures, so *s* may be taken as real and positive. The probe and conjugate modes *â* and *b̂* are then transformed as

*N̂*

_{0}〉 ≫ 1, the number operators after squeezing become

*G*≡ cosh

^{2}

*s*is the increase in probe intensity, termed the “mixing gain”.

The relative intensity operator *N̂ _{a}* –

*N̂*is unchanged by

_{b}*Ŝ*, so Var (

*N̂*–

_{a}*N̂*) = Var (

_{b}*N̂*)

_{0}. Hence the beams have been amplified without increasing the relative intensity noise; they are relative intensity squeezed. The noise figure of the process (or “degree of squeezing”) is the ratio of the measured noise to the corresponding shot-noise level for equal optical power. Assuming an initially shot-noise limited probe, the noise figure is

## 3. Optical losses after squeezing

We initially consider only losses that occur after mixing, such as from imperfect optical transmission or detection efficiency. These losses are modelled by a beamsplitter with an empty port [10] whose output state is a combination of the input and vacuum modes, contributing “vacuum fluctuations” to the transmitted beam [11]. Denoting the vacuum modes introduced by losses on the probe and conjugate by the annihilation operators *x̂* and *ŷ* respectively, the standard beam-splitter input-output relations [12] give

*η*and

_{a}*η*are the fractions of the probe and conjugate intensities transmitted. The relative intensity noise can then be expressed in terms of the individual beam variances and covariance to give

_{b}*η*≠

_{a}*η*) result in detection of amplified noise instead of squeezing.

_{b}## 4. Optical losses during squeezing: Interleaved gain/loss model

The four-wave mixing process consists of Raman transitions between the hyperfine ground states (Fig. 1A), which are most efficient when the intermediate virtual level is tuned close to resonance. However, this also increases direct absorption from the Doppler broadened transition, increasing losses and reducing correlations. To analyze this trade-off, we develop a model for the effect of competing mixing and absorption on relative intensity squeezing.

Following the approach of the numerical model presented in Ref. [3], the beam trajectories through the medium are divided into *N* discrete interleaved stages of gain and loss (Fig. 2). Distributed models of this type were first proposed by Loudon [13], and applied by Caves and Crouch [14] to model distributed squeezing losses in a single-mode parametric amplifier. We present a fully analytical model of the relative squeezing produced by this system in the continuum limit, including losses on both beams.

Each stage comprises ideal squeezing (by parameter *s*) followed by loss (represented by transmission coefficients *t _{a}* on the probe and

*t*on the conjugate). Consecutive stages of the model are related by combining Eq. (2) with Eq. (4) to give

_{b}*S*in the absence of losses, and transmissions

*T*and

_{a}*T*in the absence of squeezing. These are related to the incremental coefficients above by

_{b}*s*=

*S/N*, ${t}_{a}={T}_{a}^{1/2N}$ and ${t}_{b}={T}_{b}^{1/2N}$.

This transformation can be written in matrix form and applied recursively to express the output beam operators *â _{N}* and
${\widehat{b}}_{N}^{\u2020}$ as a sum of the incident beam operators (

*â*

_{0}and ${\widehat{b}}_{0}^{\u2020}$) and the injected vacuum operators (

*x̂*and ${\widehat{y}}_{i}^{\u2020}$) as

_{i}*ẑ*consisting of the probe annihilation and conjugate creation operators

_{i}*ẑ*

_{1}=

*â*

_{0}and ${\widehat{z}}_{2}={\widehat{b}}_{0}^{\u2020}$, followed by the injected vacuum operators

*ẑ*

_{2i+1}=

*x̂*and ${\widehat{z}}_{2i+2}={\widehat{y}}_{i}^{\u2020}$ for 1 ≤

_{i}*i*≤

*N*. Expanding Eq. (7) in terms of this operator and a set of coefficients

*α*and

_{i}*β*gives

_{i}*N̂*

_{0}≫ 1), this variance simplifies to

*N̂*

_{0}) = 〈

*N̂*

_{0}〉 is

*α*,

_{i}*β*coefficients in terms of the model parameters

_{i}*T*,

_{a}*T*and

_{b}*S*, and hence obtain an ab-initio expression for the degree of squeezing.

Equating the coefficients of Eq. (7) and Eq. (9) leads to

*N*vacuum modes

*x̂*contribute a term to the variance in Eq. (10):

_{i}*i*= 1, 2 contributions to the variance in Eq. (10), and summing over the vacuum contributions gives

*N*→ ∞. To obtain a closed form expression for the sum, the infinitesimal parameters are expanded as a power series in 1/

*N*. Expanding the elements of

**A**in Eq. (8) gives

*N*using Eq. (13) gives

*N*→ ∞, the neglected

*O*(1/

*N*

^{2}) terms vanish and

**A**

*→ exp(*

^{N}**A**

_{0}). Hence the sum

**X**converges and obeys

**X**in terms of the model parameters

*T*,

_{a}*T*and

_{b}*S*, and can be solved algebraically.

The sum **X** contains all terms in the variance of Eq. (10) except *i* = 1,2 which correspond to the probe and conjugate coefficients. The probe contribution is

*α*

_{1}

*α*

_{2}–

*β*

_{1}

*β*

_{2})

^{2}has diagonal ⌈0 1⌋. Computing the full variance sum in Eq. (10) yields

Evaluating all the contributions to Eq. (10), we obtain an analytic expression for the degree of relative intensity squeezing produced by this system. This expression is algebraic in *T _{a}*,

*T*and

_{b}*S*, but runs to a dozen typeset lines. However, special cases are readily derived and provide physical insight not readily accessible from numerical models.

In the experimentally studied case (Refs. [1–7]), detection efficiencies are carefully balanced (*η _{a}* =

*η*≡

_{b}*η*) and the far-detuned conjugate experiences negligible absorption (

*T*= 1). The corresponding degree of squeezing is

_{b}*χ*= (log

*T*– log

_{a}*T*)/4

_{b}*ξ*. The three terms describe the shot-noise limit, correlations from four-wave mixing, and injected vacuum noise.

Figure 3 shows the noise figure as a function of the probe transmission *T _{a}* and intrinsic mixing gain

*G*= cosh

^{2}

*S*(in excellent agreement with the numerical model of Ref. [3]). Note the counter-intuitive result that the strongest squeezing is obtained with imperfect transmission (

*T*< 1). This is because the shot noise carried by the incident probe beam is also amplified by the mixing process, and a small “optimal” level of probe loss decreases this contribution to the measured noise power before injected vacuum noise dominates. This optimal level is easily obtained by minimizing Eq. (16).

_{a}The model can also be applied to other four-wave mixing systems. For example, interchanging the probe and conjugate wavelengths produces the “reverse” configuration (Fig. 4A). This system is interesting as the weakly-coupled conjugate transition is brought closer to resonance, producing much higher intrinsic mixing gain for the same beam powers. In this case, the probe experiences negligible absorption (*T _{a}* = 1) compared to the conjugate (

*T*< 1), and the predicted squeezing is

_{b}Finally, it is worth noting for consistency that in the limit of both *T _{a}*,

*T*→ 1, the post-mixing optical-loss result of Eq. (5) is obtained.

_{b}## 5. Experimental diagnostic

The analytic model derived above provides a simple yet powerful tool for optimizing the degree of squeezing obtained experimentally. It not only predicts which parameters will provide optimal results, but can also be used as a diagnostic tool to determine which factors are limiting the experimentally measured degree of squeezing.

To demonstrate this, we constructed a four-wave mixing apparatus as described in Ref. [1]. A 400mW pump beam intersects a 100 *μ*W probe beam with 1/*e*^{2} beam waists of 630 *μ*m and 375 *μ*m respectively at an angle of 0.3° within a pure ^{85}Rb vapor cell of internal length 7mm heated to 130°C (Fig. 1B). The probe was generated by an AOM with fixed detuning 3040MHz below the pump, which was scanned across the Doppler broadened *D*_{1} resonance at 795nm. The relative intensity between probe and conjugate was measured with a balanced photodetector (Thorlabs PDB150A), refitted with high-efficiency photodiodes (Hamamatsu S3883, net efficiency 95%). The overall detection efficiency of the system was *η* = 85 ± 1%. The relative intensity noise was measured with a Rhode & Schwarz FSP7 spectrum analyzer at an analysis frequency of 1MHz with 30kHz resolution bandwidth.

The effective probe and conjugate gains (*G _{a}* and

*G*) were measured as a function of pump beam detuning (Fig. 5A) and used to simultaneously solve Eq. (15) for the intrinsic gain

_{b}*G*and probe transmission

*T*(Fig. 5B) via the coefficients

_{a}*α*

_{1}and

*α*

_{2}. Note that the gain resonance extends well into the Doppler-broadened absorption resonance for detunings below 600MHz from the line-centre, demonstrating strong competition between the processes.

The model predicts significantly stronger squeezing should be possible for these parameters than measured experimentally (Fig. 5C). As our incident probe beam was measured to be shot-noise limited, this demonstrates that our noise measurement was limited by technical factors and not by insufficient intrinsic gain or excessive absorption.

One cause for our discrepancy is the technical difficulty in eliminating the bright pump beam after the vapor cell. Detection of the pump introduces uncorrelated fluctuations and increases the measured noise level [3]. This is evident around zero detuning, where the probe is fully absorbed (Fig. 5A) but the measured noise is far above the standard quantum limit (Fig. 5D). Upon blocking the pump beam, the relative intensity noise was measured at the shot-noise limit. Subtracting this background level implies that −4.1dB of squeezing would be obtained by eliminating cross-beam detection.

The model assumes that gain occurs uniformly throughout the vapor cell, which requires that the beams be overlapped over the entire region. If the beams are not properly overlapped, the mixing strength decreases and gain becomes spatially varying. As losses are unchanged, the relative intensity noise increases as a result. This is likely the cause of the remaining discrepancy, and improved squeezing could be achieved by manipulating beam alignment and waists to ensure proper overlap throughout the cell.

## 6. Conclusions

We have presented a method for analytically calculating the degree of squeezing produced by a four-wave mixing system in the presence of absorption. Our model included the contributions from arbitrarily many injected vacuum modes that were subsequently squeezed by the system, producing an ab-initio quantum mechanical description of the introduced losses. Our general result is not reliant on implementation details and can be applied to analyze any four-wave mixing scheme, while the matrix methods techniques we developed can be applied to model other systems both within quantum optics and more generally.

We considered two special cases, corresponding to the experimentally studied system (introduced in Ref. [1]) and the “reverse” configuration, with probe and conjugate wavelengths interchanged. We presented closed-form expressions for the relative intensity noise in these cases, and demonstrated that a small level of probe loss was desirable to suppress amplification of the initial shot-noise. The reverse configuration was shown to produce the same squeezing at equal gain for ideal transmission, but was significantly more sensitive to losses.

The model was applied to analyze experimentally measured squeezing and determine the intrinsic mixing gain and transmission factors of the four-wave mixing resonance. Comparing the expected squeezing to measured results provided insight into the factors limiting our measurement and hence where to direct effort in optimizing the many free parameters of the system.

While it should be noted that this model considers the propagation of a Gaussian beam mode only, an arbitrary beam can be analyzed as a product state of orthogonal spatial modes, with each mode independently squeezed [4]. The model established in this paper may therefore be applied to each pair of spatial modes and the resulting noise powers summed to obtain the overall relative intensity squeezing for a multi-mode beam. Such multi-mode squeezing has been experimentally demonstrated [5], with promising applications in quantum imaging.

## References and links

**1. **C. F. McCormick, V. Boyer, E. Arimondo, and P. D. Lett, “Strong relative intensity squeezing by four-wave mixing in rubidium vapor,” Opt. Lett. **32**, 178–180 (2007). [CrossRef]

**2. **Q. Glorieux, L. Guidoni, S. Guibal, J.-P. Likforman, and T. Coudreau, “Strong quantum correlations in four wave mixing in ^{85}Rb vapor,” (SPIE, 2010), vol. 7727 of *Proc. SPIE*, p. 772703.

**3. **C. F. McCormick, A. M. Marino, V. Boyer, and P. D. Lett, “Strong low-frequency quantum correlations from a four-wave-mixing amplifier,” Phys. Rev. A **78**, 043816 (2008). [CrossRef]

**4. **V. Boyer, A. M. Marino, and P. D. Lett, “Generation of spatially broadband twin beams for quantum imaging,” Phys. Rev. Lett. **100**, 143601 (2008). [CrossRef] [PubMed]

**5. **V. Boyer, A. M. Marino, R. C. Pooser, and P. D. Lett, “Entangled images from four-wave mixing,” Science **321**, 544–547 (2008). [CrossRef] [PubMed]

**6. **V. Boyer, C. F. McCormick, E. Arimondo, and P. D. Lett, “Ultraslow propagation of matched pulses by four-wave mixing in an atomic vapor,” Phys. Rev. Lett. **99**, 143601 (2007). [CrossRef] [PubMed]

**7. **R. C. Pooser, A. M. Marino, V. Boyer, K. M. Jones, and P. D. Lett, “Low-noise amplification of a continuous-variable quantum state,” Phys. Rev. Lett. **103**, 010501 (2009). [CrossRef] [PubMed]

**8. **Q. Glorieux, R. Dubessy, S. Guibal, L. Guidoni, J.-P. Likforman, T. Coudreau, and E. Arimondo, “Double-Λ microscopic model for entangled light generation by four-wave mixing,” Phys. Rev. A **82**, 033819 (2010). [CrossRef]

**9. **C. C. Gerry and P. L. Knight, *Introductory Quantum Optics*, (Cambridge University Press, 2005).

**10. **H. A. Bachor and T. C. Ralph, *A Guide to Experiments in Quantum Optics*, (Wiley-VCH Verlag, 2004). [CrossRef]

**11. **C. M. Caves, “Quantum-mechanical radiation-pressure fluctuations in an interferometer,” Phys. Rev. Lett. **45**, 75–79 (1980). [CrossRef]

**12. **R. Loudon, *The Quantum Theory of Light* (Oxford University Press, 1983), 2nd ed.

**13. **R. Loudon, “Theory of noise accumulation in linear optical-amplifier chains,” IEEE J. Quantum Electron. **21**, 766–773 (1985). [CrossRef]

**14. **C. M. Caves and D. D. Crouch, “Quantum wideband traveling-wave analysis of a degenerate parametric amplifier,” J. Opt. Soc. Am. B **4**, 1535–1545 (1987). [CrossRef]