Field-quadrature and photon-number correlations produced by parametric processes

doi:10.1364/OE.18.019792

« Show journal navigation

Field-quadrature and photon-number correlations produced by parametric processes

C. J. McKinstrie, M. Karlsson, and Z. Tong »View Author Affiliations

Optics Express, Vol. 18, Issue 19, pp. 19792-19823 (2010)
http://dx.doi.org/10.1364/OE.18.019792

View Full Text Article

Acrobat PDF (1216 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
1. Introduction
2. General results
3. Applications
4. Discussion
5. Summary
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (26)
Cited By
Figures (15)
Metrics

Abstract

In a previous paper [Opt. Express 13, 4986 (2005)], formulas were derived for the field-quadrature and photon-number variances produced by multiple-mode parametric processes. In this paper, formulas are derived for the quadrature and number correlations. The number formulas are used to analyze the properties of basic devices, such as two-mode amplifiers, attenuators and frequency convertors, and composite systems made from these devices, such as cascaded parametric amplifiers and communication links. Amplifiers generate idlers that are correlated with the amplified signals, or correlate pre-existing pairs of modes, whereas attenuators decorrelate pre-existing modes. Both types of device modify the signal-to-noise ratios (SNRs) of the modes on which they act. Amplifiers decrease or increase the mode SNRs, depending on whether they are operated in phase-insensitive (PI) or phase-sensitive (PS) manners, respectively, whereas attenuators always decrease these SNRs. Two-mode PS links are sequences of transmission fibers (attenuators) followed by two-mode PS amplifiers. Not only do these PS links have noise figures that are 6-dB lower than those of the corresponding PI links, they also produce idlers that are (almost) completely correlated with the signals. By detecting the signals and idlers, one can eliminate the effects of electronic noise in the detectors.

1. Introduction

Parametric devices based on four-wave mixing (FWM) in fibers perform a variety of signal-processing functions that are required by communication and sensing systems [1–4

J. Hansryd, P. A. Andrekson, M. Westland, J. Li, and P. O. Hedekvist, “Fiber-based optical parametric amplifiers and their applications,” IEEE J. Sel. Top. Quantum Electron. 8, 506–520 (2002). [CrossRef]

]. Because they should perform these functions without degrading the signals, their noise properties are important. In previous papers [5

C. J. McKinstrie, S. Radic, and M. G. Raymer, “Quantum noise properties of parametric amplifiers driven by two pump waves,” Opt. Express 12, 5037–5066 (2004). [CrossRef] [PubMed]

C. J. McKinstrie, M. Yu, M. G. Raymer, and S. Radic, “Quantum noise properties of parametric processes,” Opt. Express 13, 4986–5012 (2005). [CrossRef] [PubMed]

], detailed studies were made of the properties of multiple-mode parametric processes. Such processes occur when several two-mode processes occur simultaneously in the same device, or when several two-mode devices are concatenated. Formulas were derived for the field-quadrature variances associated with homodyne detection and the photon-number variances associated with direct detection, for devices of arbitrary complexity. These formulas are sufficient for many applications, which involve measurements of individual modes (signals).

However, there are other applications which require simultaneous measurements of two or more modes. For such applications, the correlations between the mode quadratures and numbers are important. An example is two-mode squeezing, in which neither the signal nor the idler mode is squeezed by itself. Instead, squeezing exists as a correlation between the modes [7

R. Loudon and P. L. Knight, “Squeezed light,” J. Mod. Opt. 34, 709–759 (1987). [CrossRef]

, 8

S. M. Barnett and P. M. Radmore, Methods in Theoretical Quantum Optics (Oxford University Press, 1997).

]. Although the variances of the signal and idler numbers are shot-noise limited, the variance of the number difference is not [9

N. Christensen, R. Leonhardt, and J. D. Harvey, “Noise characteristics of cross-phase modulation instability light,” Opt. Commun. 101, 205–212 (1993). [CrossRef]

, 10

J. E. Sharping, M. Fiorentino, and P. Kumar, “Observation of twin-beam-type quantum correlation in optical fiber,” Opt. Lett. 26, 367–369 (2001). [CrossRef]

This paper is organized as follows: In Sec. 2, formulas are derived for the quadrature and number correlations produced by multiple-mode processes. These formulas complete the results of [5

C. J. McKinstrie, S. Radic, and M. G. Raymer, “Quantum noise properties of parametric amplifiers driven by two pump waves,” Opt. Express 12, 5037–5066 (2004). [CrossRef] [PubMed]

, 6

C. J. McKinstrie, M. Yu, M. G. Raymer, and S. Radic, “Quantum noise properties of parametric processes,” Opt. Express 13, 4986–5012 (2005). [CrossRef] [PubMed]

]. In Sec. 3, they are used to analyze the properties of basic devices, such as two-mode amplifiers, attenuators and frequency-convertors, and composite systems made from these devices, such as cascaded phase-sensitive (PS) amplifiers [11–14

R. Tang, P. Devgan, P. L. Voss, V. S. Grigoryan, and P. Kumar, “In-line frequency-nondegenerate phase-sensitive fiber-optical parametric amplifier,” IEEE Photon. Technol. Lett. 17, 1845–1847 (2005). [CrossRef]

] and communication links [15–18

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

]. Some consequences of these results are discussed in Sec. 4. (It is possible to read Sec. 4 after Sec. 2, and refer to Sec. 3 only as necessary.) Finally, in Sec. 5, the main results of this paper are summarized briefly.

2. General results

Multiple-mode parametric processes are governed by the input-output (IO) equations

b_{i} = Σ_{k} (μ_{ik} a_{k} + ν_{ik} a_{k}^{†}),

(1)

where a_i is an input-mode operator, b_i is an output-mode operator, μ_ik and ν_ik are transfer coefficients, and † is a hermitian conjugate [5

C. J. McKinstrie, S. Radic, and M. G. Raymer, “Quantum noise properties of parametric amplifiers driven by two pump waves,” Opt. Express 12, 5037–5066 (2004). [CrossRef] [PubMed]

, 6

C. J. McKinstrie, M. Yu, M. G. Raymer, and S. Radic, “Quantum noise properties of parametric processes,” Opt. Express 13, 4986–5012 (2005). [CrossRef] [PubMed]

]. The input modes satisfy the boson commutation relations [a_i ,a_j ] = 0 and [a_i ,a ^† _j] = δ_ij , where [ , ] is a commutator and δ_ij is the Kronecker delta function. The output modes satisfy similar commutation relations, which imply that

Σ_{k} (μ_{ik} ν_{jk} - ν_{ik} μ_{jk}) = 0,

(2)

Σ_{k} (μ_{ik} μ_{jk}^{*} - ν_{ik} ν_{jk}^{*}) = δ_{ij} .

(3)

The input quadrature operator

p_{i} = \frac{(a_{i} e^{- i θ_{i}} + a_{i}^{†} e^{i θ_{i}})}{2^{\frac{1}{2}}}

, where θ_i is the phase of a local oscillator (LO), and the input number operator m_i = a ^† _i a_i . (In homodyne detection, a beam splitter is used to combine a signal with a LO, and the difference between the output numbers is proportional to the input quadrature of the signal.) If the inputs are independent coherent states (CS) with amplitudes 〈a_i 〉 α_i , where 〈〉 is an expectation value, the input quadratures

〈 p_{i} 〉 = \frac{(α_{i} e^{- i θ_{i}} + α_{i}^{*} e^{i θ_{i}})}{2^{\frac{1}{2}}}

and the input numbers 〈m_i 〉 = ∣α_i ∣². [If some α_i = 0, those inputs are vacuum states (VS).]

The output operators are defined in the same way as the input operators (a_i →b_i , p_i →q_i and m_i →n_i ). For CS inputs, Eqs. (1) imply that the output amplitudes (first-order moments)

β_{i} = Σ_{k} (μ_{ik} α_{k} + ν_{ik} α_{k}^{*}) .

(4)

In general, the output strengths ∣β_i ∣² depend on the input phases ϕ_k = arg(α_k ). The output quadratures (alternative first-order moments)

〈 q_{i} (θ_{i}) 〉 = \frac{(β_{i} e^{- i θ_{i}} + β_{i}^{*} e^{i θ_{i}})}{2^{\frac{1}{2}}}

(5)

depend on both the input and LO phases.

There are two standard ways to calculate the higher-order output moments (quadrature products, numbers and number products). In the first method, one combines Eqs. (1) and calculates expectation values using the properties of CS (a_i ∣α_i 〉 = α_i ∣α_i 〉). In the second method, one rewrites the mode operators as

a_{i} = α_{i} + v_{i}, b_{i} = β_{i} + w_{i},

(6)

where the auxiliary (noise) operators ν_i and w_i also satisfy Eqs. (1) and the aforementioned commutation relations, and calculates expectation values using the properties of VS (ν_i ∣0〉=0). The second method will be used herein (because it is similar to the semi-classical method, which is familiar to communication engineers). It applies to CS inputs, but can be generalized to other inputs (such as squeezed CS).

2.1. Quadrature fluctuations

The quadrature-deviation operator

δ q_{i} (θ_{i}) = q_{i} (θ_{i}) - 〈 q_{j} (θ_{j}) 〉 = \frac{(w_{i} e^{- i θ_{i}} + w_{i}^{†} e^{i θ_{i}})}{2^{\frac{1}{2}}} .

(7)

It follows from this definition that the output quadrature correlation

〈 δ q_{i} (θ_{i}) δ q_{j} (θ_{j}) 〉 = \frac{〈 (w_{i} e^{- i θ_{i}} + w_{i}^{†} e^{i θ_{i}}) (w_{j} e^{- i θ_{j}} + w_{j}^{†} e^{i θ_{j}}) 〉}{2} .

(8)

By combining Eq. (8) with the noise moments

〈 w_{i} w_{j} 〉 = Σ_{k} μ_{ik} ν_{jk}, 〈 w_{i} w_{j}^{†} 〉 = Σ_{k} μ_{ik} μ_{jk}^{*},

〈 w_{i}^{†} w_{j} 〉 = Σ_{k} ν_{ik}^{*} ν_{jk}, 〈 w_{i}^{†} w_{j}^{†} 〉 = Σ_{k} ν_{ik}^{*} μ_{jk}^{*},

(9)

one finds that

〈 δ q_{i} (θ_{i}) δ q_{j} (θ_{j}) 〉 = \frac{Σ_{k} (μ_{ik} e^{- i θ_{i}} + ν_{ik}^{*} e^{i θ_{i}}) (μ_{jk}^{*} e^{i θ_{j}} + ν_{jk} e^{- i θ_{j}})}{2} .

(10)

When i = j, Eq. (10) reduces to Eq. (40) of [6

C. J. McKinstrie, M. Yu, M. G. Raymer, and S. Radic, “Quantum noise properties of parametric processes,” Opt. Express 13, 4986–5012 (2005). [CrossRef] [PubMed]

], the right side of which is manifestly real. When i ≠ j, the right side of Eq. (10) involves summations of

μ_{ik} μ_{jk}^{*} e^{i (θ_{j} - θ_{i})}

μ_{ik} v_{jk} e^{- i (θ_{i} + θ_{j})}

v_{ik}^{*} μ_{jk}^{*} e^{i (θ_{i} + θ_{j})}

and

v_{ik}^{*} v_{jk} e^{i (θ_{i} - θ_{j})}

. Equation (2) implies that the sum of the second and third terms is real, whereas Eq. (3) implies that the sum of the first and fourth terms is real. Hence, the quadrature-correlation formula predicts real correlations (as it must do) and reduces to the known variance formula in the appropriate limit.

2.2. Number fluctuations

The output number operator

n_{i} = {∣ β_{i} ∣}^{2} + β_{i}^{*} w_{i} + β_{i} w_{i}^{†} + w_{i}^{†} w_{i} .

(11)

The first term on the right side of Eq. (11) is the signal-signal term, the second and third terms are (collectively) the signal-noise term and the fourth term is the noise-noise term. It follows from Eq. (11) that the output number

〈 n_{i} 〉 = {∣ β_{i} ∣}^{2} + 〈 w_{i}^{†} w_{i} 〉,

(12)

where 〈w ^† _i w_i 〉 = ∑_k∣ν_ik ∣² [Eqs. (9)] is the number of noise photons.

The number-deviation operator

δ n_{i} = n_{i} - 〈 n_{i} 〉 = β_{i}^{*} w_{i} + β_{i} w_{i}^{†} + w_{i}^{†} w_{i} - 〈 w_{i}^{†} w_{i} 〉 .

(13)

It follows from this definition, and the fact that the odd-order moments of w ^(†) have zero expectation values, that the output number correlation

〈 δ n_{i} δ n_{j} 〉 = 〈 (β_{i}^{*} w_{i} + β_{i} w_{i}^{†}) (β_{j}^{*} w_{j} + β_{j} w_{j}^{†}) 〉 + 〈 w_{i}^{†} w_{i} w_{j}^{†} w_{j} 〉 - 〈 w_{i}^{†} w_{i} 〉 〈 w_{j}^{†} w_{j} 〉 .

(14)

The first term on the right side of Eq. (14) is the signal-noise term, and the second and third terms are (collectively) the noise-noise term. By comparing Eqs. (8) and (14), one finds that the signal-noise term

{〈 δ n_{i} δ n_{j} 〉}_{sn} = 2 ρ_{i} ρ_{j} 〈 δ q_{i} (ϕ_{i}) δ q_{j} (ϕ_{j}) 〉,

(15)

where ρ_j = ∣β_j ∣ and ϕ_j = arg(β_j ) are the modulus and phase of the output amplitude, respectively. Equation (15) shows that the mode phase in direct (number) detection plays the role of the LO phase in homodyne (quadrature) detection. For many applications, the signal-noise terms are much larger than the noise-noise terms (and are much easier to calculate). By combining Eqs. (9) and (15), one obtains the alternative (explicit) formula

{〈 δ n_{i} δ n_{j} 〉}_{sn} = Σ_{k} (β_{i}^{*} μ_{ik} + β_{i} ν_{ik}^{*}) (β_{j} μ_{jk}^{*} + β_{j}^{*} ν_{jk}) .

(16)

One can facilitate the evaluation of the fourth-order moment 〈w ^† _i w_iw ^† _j w_j 〉 by using the reduced noise operators

{(w_{i}^{†} w_{i})}_{1} = Σ_{k} Σ_{l} (v_{ik}^{*} μ_{il} v_{k} v_{l} + v_{ik}^{*} v_{il} v_{k} v_{l}^{†}),

(17)

{(w_{j}^{†} w_{j})}_{r} = Σ_{k} Σ_{l} (μ_{jk}^{*} v_{jl} v_{k}^{†} v_{l}^{†} + v_{jk}^{*} v_{jl} v_{k} v_{l}^{†}),

(18)

where l denotes an operator that acts (to the left) on the input ket-vector 〈0∣ and r denotes an operator that acts (to the right) on the input bra-vector ∣0〉. Notice that (w ^† _i w_i )_l = (w ^† _i w_i )^† _r (as it must do). By combining Eqs. (17) and (18), using the identities 〈ν_kν_lν ^† _m ν ^† _n〉 = δ_kmδ_ln + δ_knδ_lm and 〈ν_kν ^† _l ν_mν ^† _n〉 = δ_klδ_mn , and collecting terms, one finds that the noise-noise term

{〈 δ n_{i} δ n_{j} 〉}_{nn} = 2 Σ_{k} (μ_{ik} ν_{ik}^{*}) (μ_{jk}^{*} ν_{jk}) + Σ_{k} Σ_{l > k} (μ_{ik} ν_{il}^{*} + μ_{il} ν_{ik}^{*}) (μ_{jk}^{*} ν_{jl} + μ_{jl}^{*} ν_{jk}) .

(19)

It follows from Eqs. (16) and (19) that the output number correlation

〈 δ n_{i} δ n_{j} 〉 = Σ_{k} (β_{i}^{*} μ_{ik} + β_{i} ν_{ik}^{*}) (β_{j} μ_{jk}^{*} + β_{j}^{*} ν_{jk}) + 2 Σ_{k} (μ_{ik} ν_{ik}^{*}) (μ_{jk}^{*} ν_{jk})

+ Σ_{k} Σ_{l > k} (μ_{ik} ν_{il}^{*} + μ_{il} ν_{ik}^{*}) (μ_{jk}^{*} ν_{jl} + μ_{jl}^{*} ν_{jk}) .

(20)

When i = j, Eq. (20) reduces to Eq. (42) of [6

C. J. McKinstrie, M. Yu, M. G. Raymer, and S. Radic, “Quantum noise properties of parametric processes,” Opt. Express 13, 4986–5012 (2005). [CrossRef] [PubMed]

]. When i ≠ j, the right side of Eq. (20) involves three complicated summations. Because the number deviations commute, 〈δn_iδn_j 〉 = 〈δn_jδn_i 〉. It follows from this result, and the fact that interchanging i and j in these summations is equivalent to conjugating them, that the summations are real. (A similar argument could have been made in the context of quadrature correlations.) Hence, the number-correlation formula predicts real correlations (as it must do) and reduces to the known variance formula in the appropriate limit.

By combining the formulas for 〈δn ² _i〉 and 〈δn ² _j〉, which follow from Eq. (20), with the formula for (〈δn_iδn_j 〉 + 〈δn_jδn_i 〉)/2, which depends symmetrically on i and j, one finds that the differential variance

〈 {(δ n_{i} - δ n_{j})}^{2} 〉 = Σ_{k} {∣ β_{i}^{*} μ_{ik} + β_{i} ν_{ik}^{*} - β_{j}^{*} μ_{jk} - β_{j} ν_{jk}^{*} ∣}^{2} + 2 Σ_{k} {∣ μ_{ik} ν_{ik}^{*} - μ_{jk} ν_{jk}^{*} ∣}^{2}

+ Σ_{k} Σ_{l > k} {∣ μ_{ik} ν_{il}^{*} + μ_{il} ν_{ik}^{*} - μ_{jk} ν_{jl}^{*} - μ_{jl} ν_{jk}^{*} ∣}^{2}

(21)

is non-negative (as it must be).

3. Applications

In this section, the consequences of Eq. (16) are determined for basic devices (two-mode amplifiers, attenuators and frequency convertors), and composite systems made from these devices (copiers, cascaded PS amplifiers and PS links). Results are stated for direct detection only, because in the aforementioned applications it is more common than homodyne detection. The results for homodyne detection are similar [Eq. (15)].

For the aforementioned devices (and many others), μ_ik and ν_ik are not nonzero simultaneously, so the equation for the output number variance can be rewritten in the compact form

{〈 δ n_{i}^{2} 〉}_{sn} = {∣ β_{i} ∣}^{2} Σ_{k} {∣ λ_{ik} ∣}^{2},

(22)

where λ_ik = μ_ik if i and k are like (both odd or both even) and λ_ik = ν_ik if i and k are unlike (one odd and the other even). If i and j are like, the output number correlation

{〈 δ n_{i} δ n_{j} 〉}_{sn} = β_{i}^{*} β_{j} Σ_{k} μ_{ik} μ_{jk}^{*} + β_{i} β_{j}^{*} Σ_{l} ν_{il}^{*} ν_{jl},

(23)

where k is like i (or j) and l is unlike. It follows from Eqs. (3) and (23) that

{〈 δ n_{i} δ n_{j} 〉}_{sn} = Re (β_{i}^{*} β_{j} Σ_{k} λ_{ik} λ_{jk}^{*}),

(24)

where k is arbitrary. Conversely, if i and j are unlike,

{〈 δ n_{i} δ n_{j} 〉}_{sn} = β_{i}^{*} β_{j}^{*} Σ_{k} μ_{ik} ν_{jk} + β_{i} β_{j} Σ_{l} ν_{il}^{*} μ_{jl}^{*},

(25)

where k is like i and l is unlike. It follows from Eqs. (2) and (25) that

{〈 δ n_{i} δ n_{j} 〉}_{sn} = Re (β_{i}^{*} β_{j}^{*} Σ_{k} λ_{ik} λ_{jk}),

(26)

where k is arbitrary. Henceforth, the subscript sn will be omitted.

3.1. Two-mode amplifier

A two-mode amplifier (Fig. 1) is governed by the IO equations

b_{1} = μ a_{1} + ν a_{2}^{†},

(27)

b_{2} = ν a_{1}^{†} + μ a_{2},

(28)

where mode 1 is the signal, mode 2 is the idler, and the transfer coefficients μ and ν satisfy the auxiliary equation ∣μ∣² − ∣ν∣² = 1 [7

R. Loudon and P. L. Knight, “Squeezed light,” J. Mod. Opt. 34, 709–759 (1987). [CrossRef]

, 8

S. M. Barnett and P. M. Radmore, Methods in Theoretical Quantum Optics (Oxford University Press, 1997).

]. It is convenient to define the phase-insensitive (PI) gain G = ∣μ∣², in which case ∣ν∣² = G−1.

Fig. 1. In a two-mode amplifier (▷), the signal-mode operator (a ₁) is coupled to the hermitian conjugate of the idler-mode operator (a ^† ₂). Amplifiers are characterized by their transfer coefficients μ and ν.

If both inputs are CS, Eqs. (27) and (28) imply that the output strengths

{∣ β_{1} ∣}^{2} = {∣ μ α_{1} ∣}^{2} + {∣ ν α_{2} ∣}^{2} + 2 ∣ μ ν α_{1} α_{2} ∣ \cos θ,

(29)

{∣ β_{2} ∣}^{2} = {∣ μ α_{2} ∣}^{2} + {∣ ν α_{1} ∣}^{2} + 2 ∣ μ ν α_{1} α_{2} ∣ \cos θ,

(30)

where the phase difference

θ = ϕ_{v} - ϕ_{μ} - ϕ_{α_{1}} - ϕ_{α_{2}}

. The signal and idler (sideband) gains are maximal when θ = 0 and minimal when θ = π. Notice that ∣β ₁∣² − ∣β ₂∣² = ∣α ₁∣² − ∣α ₂∣². This relation is one of the Manley-Rowe-Weiss (MRW) equations [19

J. M. Manley and H. E. Rowe, “Some general properties of nonlinear elements—Part I. General energy relations,” Proc. IRE 44, 904–913 (1956). [CrossRef]

, 20

M. T. Weiss, “Quantum derivation of energy relations analogous to those for nonlinear reactances,” Proc. IRE 45, 1012–1013 (1957).

], and reflects the fact that sideband photons are produced in pairs.

Equations (22) and (26) imply that the output variances and correlation

〈 δ n_{1}^{2} 〉 = ({∣ μ ∣}^{2} + {∣ ν ∣}^{2}) ({∣ μ α_{1} ∣}^{2} + {∣ ν α_{2} ∣}^{2} + 2 ∣ μ ν α_{1} α_{2} ∣ \cos θ),

(31)

〈 δ n_{2}^{2} 〉 = ({∣ μ ∣}^{2} + {∣ ν ∣}^{2}) ({∣ μ α_{2} ∣}^{2} + {∣ ν α_{1} ∣}^{2} + 2 ∣ μ ν α_{1} α_{2} ∣ \cos θ),

(32)

〈 δ n_{1} δ n_{2} 〉 = 2 ∣ μ ν ∣ [∣ μ ν ∣ ({∣ α_{1} ∣}^{2} + {∣ α_{2} ∣}^{2}) + ({∣ μ ∣}^{2} + {∣ ν ∣}^{2}) ∣ α_{1} α_{2} ∣ \cos θ],

(33)

respectively. On the right sides of Eqs. (31)–(33), the first terms stem from the input fluctuations, which combine incoherently, whereas the second terms stem from the input amplitudes, which combine coherently. By combining these equations, one finds that the differential variance

〈 {(δ n_{1} - δ n_{2})}^{2} 〉 = {∣ α_{1} ∣}^{2} + {∣ α_{2} ∣}^{2}

(34)

is a constant, which equals the sum of the variances of the input CS. The number difference is constant (because photons are produced in pairs), so its output variance equals its input variance, which equals the sum of the individual variances (because the inputs are independent).

One can also explain Eq. (34) in terms of superposition modes. By redefining the phases of the input and output modes in Eqs. (27) and (28), one can replace the transfer coefficients by their moduli. The (rephased) sum and difference modes a _± = (a ₁ ± a ₂)/2^1/2 satisfy the IO equations b _± = ∣μ∣a _± ± ∣ν∣a ^† _±. In terms of these modes, the number-difference operator n ₁ − n ₂ = b ^† ₊ b ₋+b ^† ₋ b ₊ = a ^† ₊ a ₋+a ^† ₋ a ₊, which is a constant operator. It is easy to verify that 〈(δn ₁ − δn ₂)²〉 = 〈a ^† ₊ a ₊+a ^† ₋ a ₋〉 = ∣α ₁∣² + ∣α ₂∣². Physically, the number difference is proportional to the product of the sum and difference amplitudes. The sum mode is stretched and the difference mode is squeezed (by the the same amount), so the product of their amplitudes is constant.

For the common case in which the input idler is a VS (α ₂ = 0), the output amplitudes β ₁ = μα ₁ and β ₂ = να ^* ₁. The output variances and correlation

〈 δ n_{1}^{2} 〉 = ({∣ μ ∣}^{2} + {∣ ν ∣}^{2}) {∣ μ α_{1} ∣}^{2},

(35)

〈 δ n_{2}^{2} 〉 = ({∣ μ ∣}^{2} + {∣ ν ∣}^{2}) {∣ ν α_{1} ∣}^{2},

(36)

〈 δ n_{1} δ n_{2} 〉 = 2 {∣ μ ν α_{1} ∣}^{2}

(37)

and the differential variance

〈 {(δ n_{1} - δ n_{2})}^{2} 〉 = {∣ α_{1} ∣}^{2} .

(38)

In the low-gain regime (G−1≪1), the signal variance is of order ∣α ₁∣². The idler variance and correlation are much smaller, because the idler is weak. Conversely, in the high-gain regime (G ≫ 1), the variances and correlation are all of order G ²∣α ₁∣², which shows that the sidebands are strongly correlated. This correlation is also evidenced by the fact that the differential variance is much smaller than the individual variances (or that of a CS with strength G∣α ₁∣²).

For the symmetric case in which the inputs are equal (α ₁ = α ₂ = α), the (common) output strength ∣β∣² = G _θ∣α∣², where G_θ = 2G − 1 + 2[G(G − 1)]^1/2 cos θ is the PS gain. Notice that the maximal PS gain is almost 4 times higher than the PI gain. The (common) output variance and correlation

〈 δ n_{i}^{2} 〉 = ({∣ μ ∣}^{2} + {∣ ν ∣}^{2}) [({∣ μ ∣}^{2} + {∣ ν ∣}^{2}) + 2 ∣ μ ν ∣ \cos θ] {∣ α ∣}^{2},

(39)

〈 δ n_{1} δ n_{2} 〉 = 2 ∣ μ ν ∣ [2 ∣ μ ν ∣ + ({∣ μ ∣}^{2} + {∣ ν ∣}^{2}) \cos θ] {∣ α ∣}^{2} .

(40)

For in-phase input (θ = 0) the output variance and correlation are both approximately 2GG ₀∣α∣², so the sidebands are strongly correlated. Conversely, for out-of-phase input (θ = π), the variance is approximately 2GG_π ∣α∣², but the correlation is approximately −2GG _π∣α∣², so the sidebands are strongly anti-correlated. In the former case, the individual variances are much larger than ∣α∣² and a positive correlation reduces the differential variance to 2∣α∣², whereas in the latter, the individual variances are smaller than ∣α∣² and a negative correlation increases the differential variance to the aforementioned value.

A comparison of the one-and two-input results shows that the output strength is four times larger in the latter configuration than in the former, but the output variance is only two times larger, so the output signal-to-noise ratio is higher in the latter. (This result requires the inputs to be uncorrelated.)

3.2. Two-mode attenuator or frequency convertor

A two-mode attenuator (Fig. 2) is governed by the IO equations

b_{1} = τ a_{1} + ρ a_{3},

(41)

b_{3} = - ρ^{*} a_{1} + τ^{*} a_{3},

(42)

where mode 1 is the signal, mode 3 is the loss mode, and the transfer coefficients τ and ρ satisfy the auxiliary equation ∣τ∣² + ∣ρ∣² = 1 [7

R. Loudon and P. L. Knight, “Squeezed light,” J. Mod. Opt. 34, 709–759 (1987). [CrossRef]

, 8

S. M. Barnett and P. M. Radmore, Methods in Theoretical Quantum Optics (Oxford University Press, 1997).

]. It is convenient to define the transmission T = ∣τ∣², in which case ∣ρ∣² = 1 − T. Equations (41) and (42) also govern a two-mode frequency convertor, in which mode 3 is the idler [5

C. J. McKinstrie, S. Radic, and M. G. Raymer, “Quantum noise properties of parametric amplifiers driven by two pump waves,” Opt. Express 12, 5037–5066 (2004). [CrossRef] [PubMed]

, 6

C. J. McKinstrie, M. Yu, M. G. Raymer, and S. Radic, “Quantum noise properties of parametric processes,” Opt. Express 13, 4986–5012 (2005). [CrossRef] [PubMed]

, 21

C. J. McKinstrie, J. D. Harvey, S. Radic, and M. G. Raymer, “Translation of quantum states by four-wave mixing in fibers,” Opt. Express 13, 9131–9142 (2005). [CrossRef] [PubMed]

Fig. 2. In a two-mode attenuator (◁), the signal-mode operator (a ₁) is coupled to the loss-mode operator (a ₃). Attenuators are characterized by their transfer coefficients τ and ρ.

If both inputs are CS with nonzero amplitudes (frequency exchangers), Eqs. (41) and (42) imply that the output strengths

{∣ β_{1} ∣}^{2} = {∣ τ α_{1} ∣}^{2} + {∣ ρ α_{3} ∣}^{2} + 2 ∣ τ ρ α_{1} α_{3} ∣ \cos θ,

(43)

{∣ β_{3} ∣}^{2} = {∣ τ α_{3} ∣}^{2} + {∣ ρ α_{1} ∣}^{2} - 2 ∣ τ ρ α_{1} α_{3} ∣ \cos θ,

(44)

where the phase difference

θ = ϕ_{ρ} - ϕ_{τ} + ϕ_{α_{3}} - ϕ_{α_{1}}

. The conversion efficiency is minimal when θ = 0 and maximal when θ = π. Notice that ∣β ₁∣² + ∣β ₃∣² = ∣α ₁∣² + ∣α ₃∣². This MRW relation reflects the fact that signal photons are converted to idler photons (or vice versa), but the total number of sideband photons is constant.

Equations (22) and (24) imply that the output variances and correlation

〈 δ n_{i}^{2} 〉 = ({∣ τ ∣}^{2} + {∣ ρ ∣}^{2}) {∣ β_{i} ∣}^{2},

(45)

〈 δ n_{1} δ n_{3} 〉 = 0,

(46)

respectively. On the right side of Eq. (45), the first term stems from the input fluctuations, which combine incoherently (to the value 1), whereas the second term stems from the input amplitudes, which combine coherently. Equations (45) and (46) reflect the well-known fact that a frequency convertor (frequency-changing beam splitter [22

M. G. Raymer, S. J. van Enk, C. J. McKinstrie, and H. J. McGuinness, “Interference of two photons of different color,” Opt. Commun. 283, 747–752 (2010). [CrossRef]

]) converts two uncorrelated input CS to two (different) uncorrelated output CS. By combining these equations, one finds that the differential variance

〈 {({δ n}_{1} - δ n_{3})}^{2} 〉 = {∣ α_{1} ∣}^{2} + {∣ α_{3} ∣}^{2}

(47)

is a constant, which equals the sum of the variances of the input CS. Because the output numbers are uncorrelated, the differential variance equals the total variance. But the total number is constant, so its output variance equals its input variance, which equals the sum of the individual variances (because the inputs are independent). For the common special case in which α ₃ = 0 (attenuators or frequency convertors), ∣β ₁∣² = T∣α ₁∣² and ∣β ₃∣² = (1 − T)∣α ₁∣².

3.3. Amplifier followed by attenuators

For a two-mode amplifier followed by two attenuators in parallel (Fig. 3), the composite IO equations are

b_{1} = (τ_{1} μ) a_{1} + (τ_{1} ν) a_{2}^{†} + ρ_{1} a_{3},

(48)

b_{2} = (τ_{2} ν) a_{1}^{†} + (τ_{2} μ) a_{2} + ρ_{2} a_{4},

(49)

where modes 1 and 2 are the signal and idler, modes 3 and 4 are the loss modes of the signal and idler attenuators, respectively, μ and ν are the transfer coefficients of the amplifier, τ ₁ and ρ ₁ are the transfer coefficients of the signal attenuator, and τ ₂ and ρ ₂ are the transfer coefficients of the idler attenuator. The transfer coefficients satisfy auxiliary equations, which were stated in Secs. 3.1 and 3.2.

Fig. 3. Architecture of a copier. A two-mode amplifier (▷) is followed by two attenuators (◁) in parallel. Mode 1 is the signal, mode 2 is the idler (copied signal), and modes 3 and 4 are the loss modes.

If both inputs are CS, Eqs. (48) and (49) imply that the output strengths

{∣ β_{1} ∣}^{2} = {∣ τ_{1} ∣}^{2} ({∣ μ α_{1} ∣}^{2} + {∣ ν α_{2} ∣}^{2} + 2 ∣ μ ν α_{1} α_{2} ∣ \cos θ),

(50)

{∣ β_{2} ∣}^{2} = {∣ τ_{2} ∣}^{2} ({∣ μ α_{2} ∣}^{2} + {∣ ν α_{1} ∣}^{2} + 2 ∣ μ ν α_{1} α_{2} ∣ \cos θ),

(51)

where the phase difference

θ = ϕ_{ν} - ϕ_{μ} - ϕ_{α_{1}} - ϕ_{α_{2}}

. The sideband gains are maximal when θ = 0 and minimal when θ = π. Notice that the phase difference does not depend on the transmission coefficients and ∣β ₁/τ ₁∣² − ∣β ₂/τ ₂∣² = ∣α ₁∣² − ∣α ₂∣².

Equations (22) and (26) imply that the output variances and correlation

〈 δ n_{1}^{2} 〉 = {∣ τ_{1} ∣}^{2} ({∣ τ_{1} μ ∣}^{2} + {∣ τ_{1} v ∣}^{2} + {∣ ρ_{1} ∣}^{2}) ({∣ μ α_{1} ∣}^{2} {+ ∣ v α_{2} ∣}^{2} + 2 ∣ μ v α_{1} α_{2} ∣ \cos θ),

(52)

〈 δ n_{2}^{2} 〉 = {∣ τ_{2} ∣}^{2} ({∣ τ_{2} μ ∣}^{2} + {∣ τ_{2} v ∣}^{2} + {∣ ρ_{2} ∣}^{2}) ({∣ μ α_{2} ∣}^{2} {+ ∣ v α_{1} ∣}^{2} + 2 ∣ μ v α_{1} α_{2} ∣ \cos θ),

(53)

〈 δ n_{1} δ n_{2} 〉 = 2 {∣ τ_{1} τ_{2} ∣}^{2} ∣ μ v ∣ [∣ μ v ∣ ({∣ α_{1} ∣}^{2} + {∣ α_{2} ∣}^{2}) + ({∣ μ ∣}^{2} + {∣ v ∣}^{2}) ∣ α_{1} α_{2} ∣ \cos θ],

(54)

respectively. The formula for the differential variance is not illuminating.

For the common case in which the idler is a VS (α ₂ = 0) and the attenuators are identical (τ ₁ = τ ₂ = τ), the output strengths ∣β ₁∣² = TG∣α ₁∣² and ∣β ₂∣² = T(G−1)∣α ₁∣², where the PI gain G and transmittance T were defined in Secs. 3.1 and 3.2, respectively. The system is balanced (∣β ₁∣² = ∣α ₁∣²) if TG = 1, in which case loss compensates PI gain. The output variances and correlation

〈 δ n_{1}^{2} 〉 = ({∣ τ μ ∣}^{2} + {∣ τ v ∣}^{2} + {∣ ρ ∣}^{2}) {∣ τ μ α_{1} ∣}^{2},

(55)

〈 δ n_{2}^{2} 〉 = ({∣ τ μ ∣}^{2} + {∣ τ v ∣}^{2} + {∣ ρ ∣}^{2}) {∣ τv α_{1} ∣}^{2},

(56)

〈 δ n_{1} δ n_{2} 〉 = 2 {∣ τ μ τ v α_{1} ∣}^{2} .

(57)

If the reflection coefficient (ρ) were zero, Eqs. (55)–(57) would be just Eqs. (35)–(37), with μ and ν replaced by τμ and τν, respectively. However it is nonzero and the associated terms in Eqs. (55) and (56) represent the effects of loss-mode vacuum fluctuations (which do not contribute to the correlation). By combining Eqs. (55)–(57), one finds that the differential variance

〈 {(δ n_{1} - δ n_{2})}^{2} 〉 = (1 + 2 {∣ ρ v ∣}^{2}) {∣ τ α_{1} ∣}^{2} .

(58)

There is no non-trivial condition under which Eq. (58) reduces to Eq. (38). If the system is balanced, 〈δn ² ₁〉 = (3−2T)∣α ₁∣², 〈δn ² ₂〉 = (3−2T)(1−T)∣α ₁∣² and 〈δn ₁ δn ₂〉 = 2(1−T)∣α ₁∣². In the low-gain (low-loss) regime (1−T ≪ 1), the signal variance is of order ∣α ₁∣2, and the idler variance and correlation are much smaller than the signal variance: Only a weak idler is produced. Conversely, in the high-gain (high-loss) regime (T ≪ 1), both variances are approximately 3, whereas the correlation is approximately 2: The sidebands are both strong, but are not completely correlated. If the system is unbalanced (TG ≪ 1), 〈δn ² _i〉 ≈ 〈n_i 〉 and 〈δn ₁ δn ₂〉 ≪ 〈δn ² _i〉: The signal and idler are (almost) independent CS. This system is the prototype of the copiers that precede realistic PS links (Sec. 3.6). The device described in this section is more flexible than the basic amplifier described in Sec. 3.1: By choosing the values of τ ₁ and τ ₂ judiciously, one can equalize the output strengths of the sidebands and control the degree of correlation between the sidebands.

3.4. Attenuators followed by an amplifier

For two attenuators in parallel followed by a two-mode amplifier (Fig. 4), the composite IO equations are

b_{1} = (μ τ_{1}) a_{1} + (v τ_{2}^{*}) a_{2}^{†} + (μ ρ_{1}) a_{3} + (v ρ_{2}^{*}) a_{4}^{†},

(59)

b_{2} = (v τ_{1}^{*}) a_{1}^{†} + (μ τ_{2}) a_{2} + (v ρ_{1}^{*}) a_{3}^{†} + ({μ ρ}_{2}) a_{4} .

(60)

The symbols in Eqs. (59) and (60) were defined in the previous subsection.

Fig. 4. Architecture of an idealized one-stage link. Two attenuators (◁) in parallel are followed by a two-mode amplifier (▷). Mode 1 is the signal, mode 2 is the idler, and modes 3 and 4 are the loss modes.

If both inputs are CS, Eqs. (59) and (60) imply that the output strengths

{∣ β_{1} ∣}^{2} = {∣ μ τ_{1} α_{1} ∣}^{2} + {∣ v τ_{2} α_{2} ∣}^{2} + 2 ∣ μ v τ_{1} τ_{2} α_{1} α_{2} ∣ \cos θ,

(61)

{∣ β_{2} ∣}^{2} = {∣ μ τ_{2} α_{2} ∣}^{2} + {∣ v τ_{1} α_{1} ∣}^{2} + 2 ∣ μ v τ_{1} τ_{2} α_{1} α_{2} ∣ \cos θ,

(62)

where the phase difference

θ = ϕ_{ν} - ϕ_{μ} - ϕ_{τ_{1}} - ϕ_{τ_{2}} - ϕ_{α_{1}} - ϕ_{α_{2}}

. The sideband gains are maximal when θ = 0 and minimal when θ = π. Notice that ∣β ₁∣² − ∣β ₂∣² = ∣τ ₁ α ₁∣² − ∣τ ₂ α ₂∣².

Equations (22) and (26) imply that the output variances and correlation

〈 δ n_{1}^{2} 〉 = ({∣ μ ∣}^{2} + {∣ v ∣}^{2}) ({∣ μ τ_{1} α_{1} ∣}^{2} + {∣ v τ_{2} α_{2} ∣}^{2} + 2 ∣ μ v τ_{1} τ_{2} α_{1} α_{2} ∣ \cos θ),

(63)

〈 δ n_{2}^{2} 〉 = ({∣ μ ∣}^{2} + {∣ v ∣}^{2}) ({∣ μ τ_{2} α_{2} ∣}^{2} + {∣ v τ_{1} α_{1} ∣}^{2} + 2 ∣ μ v τ_{1} τ_{2} α_{1} α_{2} ∣ \cos θ),

(64)

〈 δ n_{1} δ n_{2} 〉 = 2 ∣ μ v ∣ [∣ μ v ∣ ({∣ τ_{1} α_{1} ∣}^{2} + {∣ τ_{2} α_{2} ∣}^{2}) + ({∣ μ ∣}^{2} + {∣ v ∣}^{2}) ∣ τ_{1} τ_{2} α_{1} α_{2} ∣ \cos θ],

(65)

respectively. By combining Eqs. (63)–(65), one finds that the differential variance

〈 {(δ n_{1} - δ n_{2})}^{2} 〉 = {∣ τ_{1} α_{1} ∣}^{2} + {∣ τ_{2} α_{2} ∣}^{2} .

(66)

Although Eqs. (59) and (60) involve four modes, whereas Eqs. (27) and (28) involve only two, Eqs. (61)–(66) are just Eqs. (29)–(34), with α_j replaced by τ_jα_j . This result reflects the fact that the attenuators convert input CS with amplitudes α_j to output CS with amplitudes τ_jα_j (Sec. 3.2), which are the inputs to the amplifier (Sec. 3.1).

For the symmetric case in which α ₁ = α ₂ = α and τ ₁ = τ ₂ = τ, the (common) output strength ∣β∣² = G_θT∣α∣², where the PS gain G_θ and transmittance T were defined in Secs. 3.1 and 3.2, respectively. The system is balanced (∣β∣² = ∣α∣²) if G ₀ T = 1, in which case in-phase gain compensates loss. This condition is satisfied when G = (L + 2 + 1/L)/4, where L = 1/T is the loss. The (common) output variance and correlation

〈 δ n_{i}^{2} 〉 = ({∣ μ ∣}^{2} + {∣ v ∣}^{2}) [({∣ μ ∣}^{2} + {∣ v ∣}^{2}) + 2 ∣ μ v ∣ \cos θ] {∣ τ α ∣}^{2},

(67)

〈 δ n_{1} δ n_{2} 〉 = 2 ∣ μ v ∣ [2 ∣ μ v ∣ + ({∣ μ ∣}^{2} + {∣ v ∣}^{2}) \cos θ] {∣ τ α ∣}^{2} .

(68)

If the system is balanced, 〈δn_i 〉² = (L + 1/L)∣α∣²/2 and 〈δn ₁ δn ₂〉 = (L − 1/L)∣α∣²/2. In the low-loss regime (L − 1 ≪ 1), the variance is approximately ∣α∣² and the correlation is much smaller than the variance: Only a weak correlation is produced. Conversely, in the high-loss regime (L ≫ 1), the variance and correlation are both approximately L∣α∣²/2: A strong correlation is produced (even though the amplifier only restores the sideband strengths to their input value).

The system described in the preceding paragraph is the prototype of the stages in a PS link (Sec. 3.6). It is instructive to compare the properties of this system with those of a PI link [6

C. J. McKinstrie, M. Yu, M. G. Raymer, and S. Radic, “Quantum noise properties of parametric processes,” Opt. Express 13, 4986–5012 (2005). [CrossRef] [PubMed]

, 15

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

, 18

C. J. McKinstrie, S. Radic, R. M. Jopson, and A. R. Chraplyvy, “Quantum noise limits on optical monitoring with parametric devices,” Opt. Commun. 259, 309–320 (2006). [CrossRef]

], in which the input idler is a VS (α ₂ = 0) and the output idler is discarded. For both types of link, Eq. (63) implies that the output signal variance 〈δn ² ₁〉 = (2G − ₁)〈n ₁〉. However, a balanced PI link requires that G = L, whereas a balanced PS link requires only that G ≈ L/4: The higher efficiency of PS amplification (which requires α ₂ ≠ 0) allows PS links to operate with smaller values of G than PI links, and amplify input fluctuations by smaller amounts.

3.5. Cascaded phase-sensitive amplifier

A cascaded PS amplifier [13

C. Lundström, J. Kakande, P. A. Andrekson, Z. Tong, M. Karlsson, P. Petropoulos, F. Parmigiani, and D. J. Richardson, “Experimental comparison of gain and saturation characteristics of a parametric amplifier in phase-sensitive and phase-insensitive mode,” ECOC 2009, paper 1.1.1.

, 14

J. Kakande, C. Lundström, P. A. Andrekson, Z. Tong, M. Karlsson, P. Petropoulos, F. Parmigiani, and D. J. Richardson, “Detailed characterization of a fiber-optic parametric amplifier in phase-sensitive and phase-insensitive operation,” Opt. Express 18, 4130–4137 (2010). [CrossRef] [PubMed]

] consists of a two-mode amplifier (which amplifies the input signal and generates an idler) followed by an optical processor (which controls the relative phase of the sidebands) and another two-mode amplifier (which provides PS amplification). The optical processor and connecting fibers are modeled as two attenuators in parallel. If the attenuator losses are comparable to the amplifier gain, the first two components produce a pair of sidebands whose amplitudes are comparable to the input signal amplitude. This part of the device is called a copier (Sec. 3.3). If the attenuator losses are equal, the sideband amplitudes differ slightly (because the signal and idler gains differ slightly). However, by choosing the losses judiciously, one can equalize the sideband amplitudes before the PS amplifier.

For a cascaded PS amplifier, the composite IO equations are

b_{1} = (μ τ_{1} μ_{c} + v τ_{2} v_{c}) a_{1} + (μ τ_{1} v_{c} + v τ_{2} μ_{c}) a_{2}^{†} + (μ ρ_{1}) a_{3} + (v ρ_{2}) a_{4}^{†},

(69)

b_{2} = (μ τ_{2} v_{c} + v τ_{1} μ_{c}) a_{1}^{†} + (μ τ_{2} μ_{c} + v τ_{1} v_{c}) a_{2} + ({v ρ}_{1}) a_{3}^{†} + (μ ρ_{2}) a_{4},

(70)

Fig. 5. Architecture of a cascaded phase-sensitive amplifier. A two-mode amplifier (▷) is followed by two attenuators (◁) in parallel and another two-mode amplifier. The signal, idler and loss modes are labeled 1, 2, 3 and 4, respectively.

where modes 1 and 2 are the signal and idler, respectively, modes 3 and 4 are the loss modes of the attenuators, μ_c and ν_c are the transfer coefficients of the first amplifier (copier), τ_j and ρ_j are the transfer coefficients of the attenuators, and μ and ν are the transfer coefficients of the second (PS) amplifier. The dependence of the cascaded PS amplifier on the phases of the input amplitude and transfer coefficients was studied in [23

Z. Tong, A. Bogris, C. Lundström, C. J. McKinstrie, M. Vasilyev, M. Karlsson, and P. A. Andrekson, “Modeling and measurement of the noise figure of a cascaded non-degenerate phase-sensitive parametric amplifier,” Opt. Express 18, 14820–14835 (2010). [CrossRef] [PubMed]

,24

Z. Tong, C. J. McKinstrie, C. Lundström, M. Karlsson, and P. A. Andrekson, “Noise performance of optical fiber transmission links that use non-degenerate cascaded phase-sensative amplifiers,” Opt. Express 18, 15426–15439 (2010). [CrossRef] [PubMed]

]. The results of this paper are based on the simplifying assumption that the amplitudes and coefficients are real, which is appropriate for in-phase (or out-of-phase) amplification.

If the copier is equalized (τ ₁ μ_c = τ ₂ ν_c = σ), the (common) output amplitude β = σ(μ + ν) α. Equations (22) and (26) imply that the output variances and correlation

〈 δ n_{1}^{2} 〉 = [2 σ^{2} {(μ + v)}^{2} + (μ^{2} + v^{2}) - 2 μ^{2} τ_{1}^{2}] β^{2},

(71)

〈 δ n_{2}^{2} 〉 = [2 σ^{2} {(μ + v)}^{2} + (μ^{2} + v^{2}) - 2 v^{2} τ_{1}^{2}] β^{2},

(72)

〈 δ n_{1} δ n_{2} 〉 = [2 σ^{2} {(μ + v)}^{2} + 2 μ v (1 - τ_{1}^{2})] β^{2},

(73)

respectively. The difference between the variances is 2τ ² ₁ β ², which is much smaller than the other contributions to the variances. Hence, both variances are approximately equal to the average variance, which involves the term (μ ² + ν ²)(1 − τ ² ₁). By combining Eqs. (71)–(73), one finds that the differential variance

〈 {(δ n_{1} - δ n_{2})}^{2} 〉 = 2 (1 - τ_{1}^{2}) {(σ α)}^{2} .

(74)

Equation (74) is exact. In these results, μ ² = G is the PI gain and (μ + ν)² = G ₀ is the inphase gain. These parameters are related by the identities G ₀ = 2G − 1 + 2[G(G − 1)]^1/2 and G = (G ₀ + 2 + 1/G ₀)/4.

By comparing Eqs. (71)–(73) with Eqs. (39) and (40), one finds that the σ-terms in the former equations represent the noise penalty (cost) associated with copying. (The τ ₁-terms are smaller than the σ-terms and can be omitted from this discussion.) If the copier is balanced (σ = 1), the amplitude β = (μ + ν)α, and the (common) variance and correlation

〈 δ n_{i}^{2} 〉 \approx [2 G_{0} + \frac{(G_{0} + \frac{1}{G_{0}})}{2}] β^{2},

(75)

〈 δ n_{1} δ n_{2} 〉 \approx [2 G_{0} + \frac{(G_{0} - \frac{1}{G_{0}})}{2}] β^{2} .

(76)

The (common) contribution of the copier to the variance and correlation (2G ₀ β ²) is larger than that of the PS amplifier (≈ G ₀ β ²/2). Although the variance and correlation are increased by the copier, the differential variance is not. For a balanced copier, Eq. (74) predicts that the differential variance is approximately 2α ², in agreement with Eq. (34).

In contrast, if the copier is symmetric (τ ₁ = τ ₂ = τ_c ), the output amplitudes β_i = τ_cγ_iα, where γ ₁ = (μμ_c + νν_c ) and γ ₂ = (μν_c + νμ_c ) are the transfer coefficients associated with concatenated amplifiers, which satisfy the auxiliary equation γ ² ₁ − γ ² ₂ = 1. Equations (22) and (26) imply that the output variances and correlation

〈 δ n_{i}^{2} 〉 = [τ_{c}^{2} (γ_{1}^{2} + γ_{2}^{2}) + ρ_{c}^{2} (μ^{2} + v^{2})] β_{i}^{2},

(77)

〈 δ n_{1} δ n_{2} 〉 = [τ_{c}^{2} (2 γ_{1} γ_{2}) + ρ_{c}^{2} (2 μ v)] β_{1} β_{2} .

(78)

The only difference between the variances is their proportionality to β ² _j. By combining Eqs. (77) and (78), one finds that the differential variance

〈 {({δ n}_{1} - {δ n}_{2})}^{2} 〉 = [{1 + 2 ρ_{c}^{2} (μ γ_{2} - v γ_{1})}^{2}] {(τ_{c} α)}^{2} .

(79)

It is easy to verify that μγ ₂ − νγ ₁ = ν_c , which shows the equivalence of Eqs. (58) and (79), and the similarity between Eqs. (74) and (79). For a balanced, high-gain copier (τ_cμ_c = 1 and μ ² _c ≫ 1) and a high-gain amplifier (μ ² ≫ 1), β ₂ ≈ β ₁ = (μ + ν)α, and the variances and correlation are approximately 5G ₀ β ²/2, in agreement with Eqs. (75) and (76). Just like an equalized copier, a symmetric copier increases the variances and correlation of the output sidebands (relative to those of an amplifier with CS inputs), but does not affect the differential variance significantly: The output sidebands are strongly correlated. If the second amplifier is absent (μ = 1 and ν = 0), γ ₁ = μ_c , γ ₂ = ν_c , and Eqs. (77) and (78) are consistent with Eqs. (55)–(57).

The noise properties of cascaded PS amplifiers were measured experimentally [23

, 24

], and the results are consistent with the preceding analysis (and some straightforward extensions required by the experiment).

3.6. Multiple-stage phase-sensitive link

Communication links are sequences of fibers (attenuators) and amplifiers. Links based on one-mode PS amplifiers were studied in [6

C. J. McKinstrie, M. Yu, M. G. Raymer, and S. Radic, “Quantum noise properties of parametric processes,” Opt. Express 13, 4986–5012 (2005). [CrossRef] [PubMed]

, 15

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

, 16

R. E. Slusher and B. Yurke, “Squeezed light for coherent communications,” J. Lightwave Technol. 8, 466–477 (1990). [CrossRef]

], so in this section only links based on two-mode PS amplifiers (Sec. 3.1) are considered. Each stage in an idealized link consists of two attenuators in parallel, followed by an optical processor and a two-mode amplifier (Sec. 3.4), and both inputs are CS (Fig. 6 without the copier). Hence, one can determine the composite IO equations by iterating Eqs. (59) and (60).

Fig. 6. Architecture of a multiple-stage link. The copier consists of a two-mode amplifier (▷) followed by two attenuators (◁) in parallel, whereas stage r of the link consists of two attenuators in parallel followed by a two-mode amplifier. Mode 1 is the signal, mode 2 is the idler, modes −1 and 0 are the loss modes of the copier, and modes 2r + 1 and 2r + 2 are the loss modes of the stage.

For a link with two identical stages, the IO equations are

b_{1}^{(2)} = τ^{2} (μ^{2} + ν^{2}) a_{1} + τ^{2} (2 μ v) a_{2}^{†} + ρ τ (μ^{2} + v^{2}) a_{3}

+ ρ τ (2 μ v) a_{4}^{†} + (ρ μ) a_{5} + (ρ v) a_{6}^{†},

(80)

b_{2}^{(2)} = τ^{2} (2 μ v) a_{1}^{†} + τ^{2} (μ^{2} + v^{2}) a_{2} + ρ τ (2 μ v) a_{3}^{†}

+ ρ τ (μ^{2} + v^{2}) a_{4} + (ρ v) a_{5}^{†} + (ρ μ) a_{6},

(81)

where modes 1 and 2 are the signal and idler modes, respectively, modes 3 and 4 are the loss modes of stage 1, modes 5 and 6 are the loss modes of stage 2, μ and ν are the (common) transfer coefficients of the amplifiers, τ and ρ are the (common) transfer coefficients of the attenuators, and the transfer coefficients were assumed to be real. This simplification is sufficient to model a link in which the sidebands are in-phase with the transfer coefficients, which is the case of most interest. For a three-stage link,

b_{1}^{(3)} = τ^{3} [μ (μ^{2} + 3 v^{2}) a_{1} + v (3 μ^{2} + v^{2}) a_{2}^{†}] + ρ τ^{2} [μ (μ^{2} + 3 v^{2}) a_{3} + v (3 μ^{2} + v^{2}) a_{4}^{†}]

+ ρ τ [(μ^{2} + v^{2}) a_{5} + (2 μ v) a_{6}^{†}] + ρ (μ a_{7} + {v a}_{8}^{†}),

(82)

b_{2}^{(3)} = τ^{3} [v (3 μ^{2} + v^{2}) a_{1}^{†} + μ (μ^{2} + 3 v^{2}) a_{2}] + ρ τ^{2} [v (3 μ^{2} + v^{2}) a_{3}^{†} + μ (μ^{2} + 3 v^{2}) a_{4}]

+ ρ τ [(2 μ v) a_{5}^{†} + (μ^{2} + v^{2}) a_{6}] + ρ (v a_{7}^{†} + {μ a}_{8}),

(83)

where modes 7 and 8 are the loss modes of stage 3. By continuing this sequence of equations, one finds that for an n-stage PS link, the composite IO equations are

b_{1}^{(n)} = τ^{n} (p_{n} a_{1} + q_{n} a_{2}^{†}) + ρ Σ_{r = 1}^{n} τ^{n - r} (p_{n - r + 1} a_{2 r + 1} + q_{n - r + 1} a_{2 r + 2}^{†}),

(84)

b_{2}^{(n)} = τ^{n} (q_{n} a_{1}^{†} + p_{n} a_{2}) + ρ Σ_{r = 1}^{n} τ^{n - r} (q_{n - r + 1} a_{2 r + 1}^{†} + p_{n - r + 1} a_{2 r + 2}),

(85)

where modes 2r + 1 and 2r + 2 are the loss modes of stage r. The polynomials p_n and q_n are defined by the initial conditions p ₁ = μ and q ₁ = ν, together with the recursion relations p _n+1 = μp_n + νq_n and q _n+1 = μq_n + νp_n . (These polynomials should not be confused with the input and output quadratures of Sec. 2.) It is easy to verify that p_n + q_n = (μ + ν)ⁿ and p_n − q_n = (μ − ν)ⁿ, from which it follows that

p_{n}^{2} + q_{n}^{2} = \frac{[\frac{{(μ + v)}^{2 n} + 1}{{(μ + v)}^{2 n}}]}{2},

(86)

2 p_{n} q_{n} = \frac{[\frac{{(μ + v)}^{2 n} - 1}{{(μ + v)}^{2 n}}]}{2} .

(87)

If the input sidebands are CS with amplitudes α_i , the output amplitudes

β_{1} = τ^{n} (p_{n} α_{1} + q_{n} α_{2}^{*}),

(88)

β_{2} = τ^{n} (q_{n} α_{1}^{*} + p_{n} α_{2}) .

(89)

If the input amplitudes are equal and in-phase with the transfer coefficients (positive), the (common) output amplitude β = τⁿ (μ + ν) ⁿα. Hence, the balanced-link condition [τ(μ + ν) = 1] does not depend on the number of stages in the link.

If the the inputs are equal and in-phase, Eqs. (22) and (26) imply that the (common) output variance and correlation

\frac{〈 {δ n}_{i}^{2} 〉}{β^{2}} = τ^{2 n} (p_{n}^{2} + q_{n}^{2}) + ρ^{2} Σ_{r = 1}^{n} τ^{2 (n - r)} (p_{n - r + 1}^{2} + q_{n - r + 1}^{2}),

(90)

\frac{〈 δ n_{1} δ n_{2} 〉}{β^{2}} = τ^{2 n} (2 p_{n} q_{n}) + ρ^{2} Σ_{r = 1}^{n} τ^{2 (n - r)} ({2 p}_{n - r + 1} q_{n - r + 1}),

(91)

respectively. By combining Eqs. (86) and (87) with Eqs. (90) and (91), one obtains formulas that depend explicitly on the transfer coefficients of the constituent devices. The link is balanced if G ₀ = L, where G ₀ = (μ + ν)² is the in-phase gain and L = 1/τ ² is the loss of each stage. By making these substitutions in Eqs. (90) and (91), and doing the summations, one finds that

〈 δ n_{i}^{2} 〉 = \frac{[1 + \frac{1}{L^{2 n}} + n (L - 1) + \frac{(1 - \frac{1}{L^{2 n}})}{(L + 1)}] {∣ α ∣}^{2}}{2},

= \frac{[1 + n (L - 1) + \frac{(1 + \frac{1}{L^{2 n - 1}})}{(L + 1)}] {∣ α ∣}^{2}}{2},

(92)

〈 δ n_{1} δ n_{2} 〉 = \frac{[1 - \frac{1}{L^{2 n}} + n (L - 1) - \frac{(1 - \frac{1}{L^{2 n}})}{(L + 1)}] {∣ α ∣}^{2}}{2},

= \frac{[1 + n (L - 1) - \frac{(1 + \frac{1}{L^{2 n - 1}})}{(L + 1)}] {∣ α ∣}^{2}}{2} .

(93)

In the first forms of Eqs. (92) and (93), the first pairs of terms stem from the signal and idler fluctuations, which are transmitted through the whole link, whereas the second pairs stem from the loss-mode (vacuum) fluctuations, which are added throughout the link. In the low-loss regime (L − 1 ≪ 1), 〈δn ² _i〉 ≈ ∣α∣² and 〈δn ₁ δn ₂〉 ≈ 0: The outputs are (almost) independent CS, as were the inputs. Conversely, in the high-loss regime (L ≫ 1), 〈δn ² _i〉 ≈ [1 + n(L − 1)]∣α∣²/2 ≈ 〈δn ₁ δn ₂〉: The output variance and correlation are determined primarily by the total loss nL and the outputs are (almost) completely correlated. At the inputs to the second and subsequent stages, the sidebands are strongly correlated. In each stage, the sidebands are diminished and decorrelated by the attenuators, then augmented and recorrelated by the PS amplifiers. (If the sidebands were not decorrelated at the inputs to the amplifiers, the factors of 1/2 would be absent from variance and correlation formulas.) For the special case in which n = 1, 〈δn ² _i〉 = (L + 1/L)∣α∣²/2 and 〈δn ₁ δn ₂〉 = (L − 1/L)∣α∣²/2. These results are consistent with Eqs. (67) and (68).

By combining Eqs. (90) and (91), one finds that the differential variance

\frac{〈 {(δ n_{1} - δ n_{2})}^{2} 〉}{β^{2}} = 2 τ^{2 n} {(p_{n} - q_{n})}^{2} + 2 ρ^{2} Σ_{r = 1}^{n} τ^{2 (n - r)} {(p_{n - r + 1} - q_{n - r + 1})}^{2} .

(94)

For a balanced link,

〈 {(δ n_{1} - δ n_{2})}^{2} 〉 = 2 [\frac{1}{L^{2 n}} + \frac{(1 - \frac{1}{L^{2 n}})}{(L + 1)}] {∣ α ∣}^{2}

= \frac{2 (1 + \frac{1}{L^{2 n - 1}}) {∣ α ∣}^{2}}{(L + 1)} .

(95)

In the first form of Eq. (95), the first term stems from sideband fluctuations and decreases rapidly as L increases (because the sideband fluctuations are strongly correlated), whereas the second term stems from vacuum fluctuations and decreases slowly as L increases (because uncorrelated vacuum fluctuations are added throughout the link). In the low-loss regime, the differential variance (≈ 2α ²) has the value associated with two independent CS. In the high-loss regime, the differential variance (≈ 2α ²/L) does not depend the number of stages in the link. For the special case in which n = 1, Eq. (95) reduces to Eq. (66).

Conventional communication systems are based on single-carrier-frequency signals. Hence, in realistic PS links, the signals must be copied (idlers must be generated) before they are transmitted. For a copier (Sec. 3.3) followed by an n-stage PS link (Fig. 6), the composite IO equations are

b_{1} = τ^{n} [(p_{n} τ_{1} μ_{c} + q_{n} τ_{2} v_{c}) a_{1} + (p_{n} τ_{1} v_{c} + q_{n} τ_{2} μ_{c}) a_{2}^{†} + (p_{n} ρ_{1}) a_{- 1}

+ (q_{n} ρ_{2}) a_{0}^{†}] + ρ Σ_{r = 1}^{n} τ^{n - r} (p_{n - r + 1} a_{2 r + 1} + q_{n - r + 1} a_{2 r + 2}^{†}) .

(96)

b_{2} = τ^{n} [(p_{n} τ_{2} v_{c} + q_{n} τ_{1} μ_{c}) a_{1}^{†} + (p_{n} τ_{2} μ_{c} + q_{n} τ_{1} v_{c}) a_{2} + (q_{n} ρ_{1}) a_{- 1}^{†}

+ (p_{n} ρ_{2}) a_{0}] + ρ Σ_{r = 1}^{n} τ^{n - r} (p_{n - r + 1} a_{2 r + 1}^{†} + q_{n - r + 1} a_{2 r + 2}),

(97)

where modes −1 and 0 are the loss modes of the copier attenuators, μ_c and ν_c are the transfer coefficients of the copier amplifier, τ_j and ρ_j are the transfer coefficients of the copier attenuators, and the output superscripts (n) were omitted. All the other symbols were defined above. The ρ-terms in Eqs. (96) and (97) are identical to those in Eqs. (84) and (85), as are their contributions to the output variances and correlations. (These terms are associated with the loss modes of the link.) Hence, only the first four (copier) terms are retained the following analysis.

For an equalized copier (τ ₁ μ_c = τ ₂ ν_c = σ), the (common) output amplitude β = στⁿ (p_n + q_n )α. If the copier and link are balanced [σ = 1 and (μ + ν)τ = 1, respectively], so also is the combined system, in which case β = α. The contributions of the copier terms to the output variances and correlation are

\frac{〈 δ n_{1}^{2} 〉}{β^{2}} = τ^{2 n} [2 σ^{2} {(p_{n} + q_{n})}^{2} + p_{n}^{2} + q_{n}^{2} - 2 p_{n}^{2} τ_{1}^{2}],

(98)

\frac{〈 δ n_{2}^{2} 〉}{β^{2}} = τ^{2 n} [2 σ^{2} {(p_{n} + q_{n})}^{2} + p_{n}^{2} + q_{n}^{2} - 2 q_{n}^{2} τ_{1}^{2}],

(99)

\frac{〈 δ n_{1} δ n_{2} 〉}{β^{2}} = τ^{2 n} [2 σ^{2} {(p_{n} + q_{n})}^{2} + (2 p_{n} q_{n}) (1 - τ_{1}^{2})],

(100)

respectively. The (normalized) signal and idler variances differ by the amount 2τ ²ⁿ τ ² ₁, which is negligible, so both variances are approximately equal to the average variance, which involves the term (p ² _n + q ² _n)(1 − τ ² ₁). By combining Eqs. (98)–(100), one finds that the copier contributions to the differential variance are

\frac{〈 {(δ n_{1} - δ n_{2})}^{2} 〉}{β^{2}} = 2 τ^{2 n} {(p_{n} - q_{n})}^{2} (1 - τ_{1}^{2}) .

(101)

Equation (101) is exact. For the case in which n = 1 and τ = 1, Eqs. (98)–(101) reduce to Eqs. (71)–(74), respectively.

By comparing Eqs. (90) and (91) with Eqs. (98)–(100), one finds that the copier increases the (normalized) variance and correlation by the amounts 2σ ² τ ²ⁿ(μ + ν)²ⁿ and −τ ²ⁿ[(μ + ν)²ⁿ + 1/(μ + ν)²ⁿ]τ ² ₁/2. For a balanced system, the second amount (≈ −τ ² ₁/2) is negligible and the first amount (2) is much smaller than the other contributions (≈ nL/2). By comparing Eqs. (94) and (101), one finds that the copier has a negligible impact on the differential variance (because it produces correlated sideband photons).

For a symmetric copier (τ ₁ = τ ₂ = τ_c ), the output amplitudes β_i = τⁿτ_cγ_iα, where γ ₁ = (p_nμ_c + q_nν_c ) and γ ₂ = (p_nν_c + q_nμ_c ) are the transfer coefficients associated with concatenated amplifiers, which satisfy the auxiliary equation stated in Sec. 3.5. By using the properties of the constituent transfer coefficients, one finds that

γ_{1}^{2} + γ_{2}^{2} = \frac{[{(μ_{c} + ν_{c})}^{2} {(μ + ν)}^{2 n} + \frac{1}{{(μ_{c} + ν_{c})}^{2} {(μ + ν)}^{2 n}}]}{2},

(102)

2 γ_{1} γ_{2} = \frac{[{(μ_{c} + ν_{c})}^{2} {(μ + ν)}^{2 n} - \frac{1}{{(μ_{c} + ν_{c})}^{2} {(μ + ν)}^{2 n}}]}{2} .

(103)

Equations (102) and (103) are extensions of Eqs. (86) and (87), respectively. If the system is balanced (τⁿτ_cγ ₁ = 1), the amplitudes β ₁ = α and β ₂ = αγ ₂/γ ₁. The contributions of the copier terms to the output variances and correlation are

\frac{〈 δ n_{i}^{2} 〉}{β_{i}^{2}} = τ^{2 n} [τ_{c}^{2} (γ_{1}^{2} + γ_{2}^{2}) + ρ_{c}^{2} (p_{n}^{2} + q_{n}^{2})],

(104)

\frac{〈 δ n_{1} δ n_{2} 〉}{β_{1} β_{2}} = τ^{2 n} [τ_{c}^{2} (2 γ_{1} γ_{2}) + ρ_{c}^{2} (2 p_{n} q_{n})],

(105)

respectively. The only difference between the variances is their dependence on β ² _j. By combining Eqs. (104) and (105), one finds that the copier contributions to the differential variance are

\frac{〈 {(δ n_{1} - δ n_{2})}^{2} 〉}{{(τ^{n} τ_{c} α)}^{2}} = τ^{2 n} [1 + 2 ρ_{c}^{2} {(p_{n} γ_{2} - q_{n} γ_{1})}^{2}],

(106)

where p_nγ ₂ − q_nγ ₁ = ν_c . For the case in which n = 1 and τ = 1, Eqs. (104)–(106) reduce to Eqs. (77)–(79), respectively.

By comparing Eqs. (90) and (91) with Eqs. (104) and (105), one finds that the copier increases the (normalized) variance and correlation by the amounts (τⁿτ_c )²(γ ² ₁ + γ ² ₂ − p ² _n − q ² _n) and (τⁿτ_c )²(2γ ₁ γ ₂ − 2p_nq_n ), respectively. For a balanced system with a high-gain copier (μ ² _c ≫ 1), the (common) copier contribution to the variance and correlation (≈ 2) is much smaller than the other contributions (≈ nL/2). By comparing Eqs. (94) and (106), one finds that the copier makes a negligible contribution to the differential variance. Thus, the performance of a copier and n-stage PS link does not depend on whether the copier is equalized or symmetric and depends only weakly on the whether the copier is present.

For reference, if the link is balanced, Eqs. (102)–(105) imply that

\frac{〈 δ n_{i}^{2} 〉}{β_{i}^{2}} = \frac{[T_{c} (G_{c 0} + \frac{1}{G_{c 0} L^{2 n}}) + (1 - T_{c}) (1 + \frac{1}{L^{2 n}})]}{2},

(107)

\frac{〈 δ n_{1} δ n_{2}^{2} 〉}{β_{1} β_{2}} = \frac{[T_{c} (G_{c 0} - \frac{1}{G_{c 0} L^{2 n}}) + (1 - T_{c}) (1 - \frac{1}{L^{2 n}})]}{2},

(108)

where the in-phase copier gain G _c0 = 2G_c − 1 + 2[G_c (G_c − 1)]^1/2.

4. Discussion

For direct detection, the signal-to-noise ratio of mode i is 〈n_i 〉²/〈δn ² _i〉 and the noise figure associated with mode i is the input ratio of the signal divided by the output ratio of mode i. The correlation coefficient of modes i and j is 〈δn_iδn_j 〉/(〈δn ² _i〉〈δn ² _j〉)^1/2. For a two-mode PI amplifier, which has one CS and one VS input, Eqs. (35)–(37) imply that the sideband noise figures

F_{i} = \frac{(2 G - 1)}{G_{i}},

(109)

where G ₁ = G and G ₂ = G − 1 are the PI signal and idler gains, respectively. The correlation coefficient

C_{12} = \frac{2 {[G (G - 1)]}^{\frac{1}{2}}}{(2 G - 1)} .

(110)

These quantities are plotted as functions of the PI gain parameter G in Fig. 7. As the gain increases, the signal noise figure increases monotonically and the idler noise figure decreases monotonically, to their (common) asymptotic value of 2 (3 dB). This factor of 2 arises because only the signal contributes coherent components to the outputs, whereas the signal and idler both contribute incoherent components (noise). The correlation coefficient tends to 1 rapidly as the gain increases (C ₁₂ ≈ 1 − 1/8G ²), so only a moderate gain is required to produce a strong correlation.

Fig. 7. Properties of a two-mode PI amplifier. (a) Noise figures plotted as functions of gain. The solid and dashed curves represent the signal and idler, respectively. (b) Correlation coefficient plotted as a function of gain.

For a two-mode PS amplifier, which has two CS inputs of equal strength, Eqs. (39) and (40) imply that the (common) noise figure

F = \frac{(2 G - 1)}{G_{θ}},

(111)

where G_θ = 2G − 1 + 2[G(G − 1)]^1/2 cos θ is the PS gain and θ is the phase difference between the pumps and sidebands. G_θ ≥ 1 unless cos θ < 0 and G < 1/(1 − cos² θ). The correlation coefficient

C_{12} = \frac{2 {[G (G - 1)]}^{\frac{1}{2}} C_{θ}}{(2 G - 1) G_{θ}},

(112)

where C_θ = 2[G(G − 1)]^1/2 + (2G − 1)cos θ describes the phase dependence of the correlation. C_θ ≥ 0 unless cos θ < 0 and G < [1 + 1/(1 − cos² θ)^1/2]/2. If the inputs are in-phase with the pumps (θ = 0), C ₀ = 2G − 1 + 2[G(G − 1)]^1/2 = G ₀, whereas if the inputs are out-of-phase (θ = π), C ₀ = −G ₀. The noise figure and correlation coefficient are plotted as functions of the PI gain in Fig. 8. If θ = 0, the noise figure tends rapidly to 1/2 (−3 dB) as the gain increases, because the coherent components of the sidebands increase twice as rapidly as the noise components. (Their gain factors are 4G and 2G, respectively.) This result requires the inputs to be independent (uncorrelated). The correlation coefficient tends rapidly to 1, so (once again) only a moderate gain is required to produce a strong correlation. In contrast, if θ = π, the noise figure increases, because the coherent components decrease while the noise components increase, and the correlation coefficient tends rapidly to −1. The case in which θ = π/2 is intermediate.

Fig. 8. Properties of a two-mode PS amplifier. (a) Noise figures and (b) correlation coefficients plotted as functions of gain. The solid, dot-dashed and dashed curves represent the phase differences 0 (in-phase), π/2 and π (out-of-phase), respectively.

For a two-mode attenuator (frequency convertor), in which the input signal is a CS and the input loss mode (idler) is a VS, Eqs. (45) and (46) imply that the noise figures

F_{i} = \frac{1}{T_{i}},

(113)

where T ₁ = T and T ₃ = 1 − T are the signal and loss-mode (idler) transmissions, respectively. The correlation coefficient

C_{13} = 0,

(114)

for all values of the transmissions. These results reflect the fact that the outputs are also independent CS. The sideband noise figures are plotted as functions of the loss parameter L = 1/T in Fig. 9. F ₁ equals the loss parameter. It exceeds 1 because the attenuator does not decrease the noise component of the output signal as much as the coherent component. (If it did, the properties of a strongly-attenuated signal would violate the Heisenberg uncertainty principle.) In the high-loss regime (L ≫ 1), the coherent component of the output loss mode (idler) is comparable to that of the input signal, so F ₃ ≈ 1. In an attenuator the loss mode is inaccessible, whereas in a frequency convertor the idler is an accessible copy of the signal.

Fig. 9. Noise figures plotted as functions of loss for a two-mode attenuator. The solid and dashed curves represent the signal and loss-mode, respectively.

For a copier (PI amplifier followed by two attenuators in parallel), Eqs. (55)–(57) imply that the noise figures

F_{i} = \frac{[T (2 G - 1) + (1 - T)]}{T G_{i}}

(115)

and the correlation coefficient

C_{12} = \frac{2 T {[G (G - 1)]}^{\frac{1}{2}}}{T (2 G - 1) + (1 - T)} .

(116)

Equations (115) and (116) are based on the simplifying assumption that the attenuators are identical. (One could also use attenuators with slightly different transmissions to equalize the output strengths of the sidebands. However, the performance of the copier does not depend sensitively on whether it is equalized or symmetric.) If the attenuators are absent (T = 1), Eqs. (115) and (116) reduce to Eqs. (109) and (110), respectively. Conversely, if the amplifier is absent (G = 1), the signal version of Eq. (115) reduces to the signal version of Eq. (113). For a balanced copier (TG = 1), F ₁ ≈ 3 − 2T and F ₂ ≈ (3 − 2T)/(1 − T) and C ₁₂ = 2(1 − T)^1/2/(3 − 2T). These quantities are plotted as functions of gain in Fig. 10. As the gain increases, the signal noise figure increases, and the idler noise figure decreases, to their (common) asymptotic value of 3 (4.8 dB). The correlation coefficient starts to increase as the gain increases. However, its growth saturates rapidly as it approaches its asymptotic value of 2/3, because the noise added by the attenuators is uncorrelated. A balanced copier produces output sidebands that have comparable strengths and are partially correlated. For an unbalanced copier (TG ≪ 1), F_i ≈ 1/TG_i and C ₁₂ ≈ 2TG. The output sidebands are much noisier than the input signal (as befits weak nearly-CS), and are only weakly correlated.

Fig. 10. Properties of a balanced copier. (a) Noise figures plotted as functions of gain. The solid and dashed curves represent the signal and idler, respectively. (b) Correlation coefficient plotted as a function of gain.

For the first stage of an idealized PS link (parallel attenuators followed by a two-mode PS amplifier), which has two CS inputs of equal strength, Eqs. (67) and (68) imply that the (common) noise figure

F = \frac{(2 G - 1)}{T G_{θ}},

(117)

where G_θ was defined after Eq. (111). The correlation coefficient

C_{12} = \frac{2 {[G (G - 1)]}^{\frac{1}{2}} C_{θ}}{(2 G - 1) G_{θ}},

(118)

where C_θ was defined after Eq. (112). The right side of Eq. (117) is just the right side of Eq. (111) divided by T, because the attenuators replace CS by diminished CS. Equation (118) is identical to Eq. (112), because the attenuators do not influence the correlation, which is produced by the amplifier. The noise figure and correlation coefficient are plotted as functions of loss in Fig. 11, for a balanced link (TG ₀ = 1). Results are also included for the associated PI link [6

C. J. McKinstrie, M. Yu, M. G. Raymer, and S. Radic, “Quantum noise properties of parametric processes,” Opt. Express 13, 4986–5012 (2005). [CrossRef] [PubMed]

]. As the loss increases, the PS and PI noise figures both increase. However, the PS noise figure increases less rapidly and, for large values of loss, is about 6-dB lower. The PS correlation coefficient also increases less rapidly than its PI counterpart. The PI and PS formulas for these quantities have the same dependences on the gain. Differences between the displayed results exist only because the relations between the gain and loss are different for the two links (G = L and G ≈ L/4, respectively). This example and the preceding one show that the order of amplification and attenuation is important.

Fig. 11. Properties of a balanced one-stage link. (a) Noise figures and (b) correlation coefficients plotted as functions of loss. The solid and dashed curves represent PS and PI links, respectively.

For a cascaded PS amplifier (PI amplifier followed by parallel attenuators and a PS amplifier), Eqs. (77) and (78) imply that the signal noise figure

F_{1} = \frac{[T (2 H_{0} - 1) + (1 - T) (2 G_{2} - 1)]}{T H_{0}},

(119)

where H ₀ = G ₂ G ₁ + (G ₂ − 1)(G ₁ − 1) + 2[G ₂(G ₂ − 1)G ₁(G ₁ − 1)]^1/2 is the in-phase gain of both amplifiers, and G ₁ and G ₂ are the PI gains of the first and second amplifiers, respectively. The formula for the idler noise figure is similar. (In the denominator, H ₀ is replaced by H ₀ − 1.) The correlation coefficient

C_{12} = \frac{2 T {[H_{0} (H_{0} - 1)]}^{\frac{1}{2}} + 2 (1 - T) {[G_{2} (G_{2} - 1)]}^{\frac{1}{2}}}{T (2 H_{0} - 1) + (1 - T) (2 G_{2} - 1)} .

(120)

If the second amplifier is absent (H ₀ = G ₁), Eqs. (119) and (120) reduce to Eqs. (115) and (116), respectively. Conversely, if the first amplifier is absent (H ₀ = G ₂), Eq. (119) reduces to Eq. (109), adjusted for a diminished input signal, and Eq. (120) reduces to Eq. (110).

If the gains are high (G ₁ and G ₂ ≫ 1), H ₀ ≈ 4G ₁ G ₂, the (common) noise figure F ≈ [8G ₂ TG ₁ + 2G ₂(1 − T)]/4G ₂ TG ₁ ≥ 2 and the correlation coefficient C ₁₂ ≈ 1. The outputs have comparable strengths and are strongly correlated, for all values of the gain ratios and transmission. For a low-loss copier (1 − T ≪ 1), F ≈ 2 (3.0 dB), so the cascaded PS amplifier has the same noise figure as a PI amplifier (but different PS gain). For a balanced copier (TG ₁ = 1) with significant loss (T ≪ 1), F ≈ 2.5 (4.0 dB), which is 1.0 dB higher than the noise figure of a PI amplifier, but is 0.8 dB lower than that of the constituent copier. The second (PS) amplifier reduces the sideband noise, even though its inputs (which are the copier outputs) are partially correlated. The sideband noise figures and correlation coefficient are plotted as functions of the PI gain of the second amplifier in Fig. 12. Notice that only moderate values of G ₂ are required to establish the properties of the cascaded PS amplifier:With the exception of the net gain TH ₀, the properties of this device depend only weakly on G ₂. For an unbalanced copier with high loss (TG ₁ ≪ 1), F ≈ 1/2TG ₁, which is much higher than the noise figure of a PI amplifier, but is lower than that of the copier by a factor of 2 (3.0 dB). The second amplifier achieves its maximal noise reduction because its inputs are uncorrelated. The noise figure and correlation coefficient are plotted as functions of loss in Fig. 13. Results are also included for the copier alone. For small values of loss, the noise figures of both devices are approximately 3 dB, which is appropriate for high-gain amplifiers. As the loss increases, the noise figure of the cascaded PS amplifier increases more slowly than that of the copier. This result demonstrates (again) the beneficial effects of PS amplification. For the cascaded PS amplifier, the correlation coefficient is (almost) independent of loss, because its last constituent device is an amplifier. In contrast, for the copier alone, the correlation decreases as the loss increases.

Fig. 12. Properties of a cascaded PS amplifier with a balanced copier. (a) Noise figure and (b) correlation coefficient plotted as functions of the gain G ₂, for the case in which G ₁ = 10 and the transmission T = 0.1.

For an idealized n-stage PS link (sequence of parallel attenuators and two-mode PS amplifiers), which has two CS inputs of equal strength, Eqs. (92) and (93) imply that the (common) noise figure

F \approx \frac{[1 + n (L - 1)]}{2},

(121)

Fig. 13. Properties of a cascaded PS amplifier with high loss. (a) Noise figure and (b) correlation coefficient plotted as functions of loss, for the case in which the gains G ₁ = 10 and G ₂ = 10. The solid and dashed curves represent the PS amplifier and an unbalanced copier, respectively.

where L is the stage loss, and the correlation coefficient

C_{12} \approx 1 - \frac{2}{(L + 1) [1 + n (L - 1)]} .

(122)

Equations (121) and (122) are based on the (realistic) assumptions that the loss is high (L ≫ 1) and the link is balanced (G ₀ = L, where G ₀ is the in-phase stage gain). To be precise, for high losses F ≈ nL/2 and C ₁₂ ≈ 1 − 2/nL ². The extra terms in Eqs. (121) and (122) increase their ranges of accuracy. The noise figure and correlation coefficient are plotted as functions of the stage loss in Fig. 14. Notice that the approximate results [Eqs. (121) and (122)] agree well with the exact results [Eqs. (92) and (93)], even for moderate values of loss. The asymptotic (high-loss) noise figure of an idealized PS link (nL/2) is smaller than that of a PI link (2nL) by a factor of 4 (6.0 dB). For both types of link, the noise figure is 2nG, where G is the PI stage gain. However, a balanced PI link requires that G = L, whereas a balanced PS link requires that G ≈ L/4: The higher efficiency of PS amplification allows PS links to operate with smaller values of G than PI links, and produce less noise. In each stage, the sidebands are diminished and decorrelated by the attenuators, then augmented and recorrelated by the amplifiers. The correlation evolution is nearly periodic, so the correlation coefficient depends only weakly on the number of stages.

Fig. 14. Properties of an idealized 3-stage PS link. (a) Noise figure and (b) correlation coefficient plotted as functions of the stage loss. The solid curves represent the exact results [Eqs. (92) and (93)], whereas the dashed curves represent the approximate results [Eqs. (121) and (122)].

A realistic n-stage PS link (copier followed by a sequence of parallel attenuators and PS amplifiers) has one CS input and one VS input (because the copier generates the idler that is required by the PS link). As stated above, the noise properties of a copier depend only weakly on whether it is equalized or symmetric. For a symmetric copier and a balanced link, Eqs. (92), (93), (107) and (108) imply that the (common) noise figure

F \approx \frac{2 [T_{c} G_{c 0} + (1 - T_{c}) + n (L - 1)]}{T_{c} G_{c 0}},

(123)

where G _c0 = 2G_c − 1 + 2[G_c (G_c − 1)]^1/2 is the in-phase copier gain, G_c is the PI copier gain and T_c is the copier transmission. The correlation coefficient

C_{12} \approx 1 - \frac{2}{(L + 1) [T_{c} G_{c 0} + (1 - T_{c}) + n (L - 1)]} .

(124)

For a balanced high-gain copier (T_cG_c = 1 and G_c ≫ 1, so G _c0 ≈ 4G_c ), F ≈ [5 + n(L − 1)]/2 and C ₁₂ ≈ 1 − 2/(L + 1)[5 + n(L − 1)]. The noise figure and correlation coefficient are plotted as functions of the stage loss in Fig. 15. Once again, the approximate results [Eqs. (123) and (124)] agree well with the exact results [Eqs. (107) and (108)]. For small values of loss the effects of the copier are significant: The noise figure is 3 (rather than 1) and the correlation coefficient is 2/3 (rather than 0). In contrast, for large values of loss the copier effects are minor: The sidebands are still strongly correlated (C ₁₂ ≈ 1 − 2/nL ²) and the copier contribution to the noise figure (2) is much less than the link contribution (nL/2). Thus, one can use a copier to generate the idler that is required by a PS link, without sacrificing the 6-dB advantage of the link (relative to a PI link). This positive result was obtained because the moderate amount of noise produced by the copier is swamped by the large amount of noise produced by the link (loss and gain).

Fig. 15. Properties of a realistic 3-stage PS link with a balanced copier. (a) Noise figure and (b) correlation coefficient plotted as functions of the stage loss. The solid curves represent the exact results [Eqs. (107) and (108)], whereas the dashed curves represent the approximate results [Eqs. (123) and (124)].

A cascaded PS attenuator is a cascaded PS amplifier with more loss than gain. For such a device, it was shown above that F ≈ 1/2TG ₁. If one were to split the total transmission T into two parts, the first T_c associated with copying and the second 1/L_t associated with transmission, and relabel G ₁ as G_c , one would find that F ≈ L_t /2T_cG_c , which is just Eq. (123) with nL replaced by L_t . Furthermore, Eqs. (120) and (124) both imply that C ₁₂ ≈ 1. Thus, the noise properties of a cascaded PS attenuator, which is straightforward to construct, mimic those of a realistic PS link, which is difficult to construct.

A common feature of devices that use two-mode amplifiers is the strong correlation between the output signal and idler. By measuring both sidebands, one can subtract the effects of electrical noise in the detectors [25

L. A. Krivitsky, U. L. Andersen, R. Dong, A. Huck, C. Wittmann, and G. Leuchs, “Electronic noise-free measurements of squeezed light,” Opt. Lett. 33, 2395–2397 (2008). [CrossRef] [PubMed]

] and improve the performances of these devices.

5. Summary

In this paper, formulas were derived for the field-quadrature and photon-number variances and correlations produced by multiple-mode parametric processes. These formulas were used to analyze the properties of basic devices, such as two-mode amplifiers, attenuators and frequency convertors, and composite systems made from these devices, such as cascaded parametric amplifiers and communication links. For these systems (and many others), the general formulas for the variances and correlations [Eqs. (16) and (19)] simplify significantly [Eq. (22) and Eq. (24) or (26)].

Two-mode amplifiers with one coherent-state (CS) input (signal) and one vacuum-state (VS) input are phase insensitive (PI), so the output signal powers do not depend on the input signal phases. These amplifiers generate idlers that are correlated with the amplified signals [Eqs. (35)–(37)]. The noise-figure of a device is the input signal-to-noise ratio divided by the output ratio. PI amplifiers have (high-gain) noise figures of 3 dB [Eq. (109)], because they add amplified VS fluctuations to the signals. In contrast, amplifiers with two CS inputs are phase sensitive (PS). These amplifiers correlate their input modes [Eqs. (39) and (40)] and have (in-phase) noise figures of −3 dB [Eq. (111)], because they combine the input amplitudes coherently, but only combine the CS fluctuations incoherently. (This remarkable performance is only possible if the inputs are uncorrelated.) Two-mode attenuators with one CS input and one VS input produce two uncorrelated CS outputs [Eqs. (45) and (46)]. The noise figure of an attenuator equals its loss factor [Eq. (113)]. Two attenuators acting in parallel on correlated modes (such as those produced by amplifiers) decorrelate the modes.

To operate in a PS manner, a two-mode amplifier requires two nonzero inputs. However, current communication systems are based on one-carrier-frequency signals. A standard way to produce the second input (copy the signal) is to use a two-mode amplifier, which amplifies the signal and generates an idler of comparable strength, and two attenuators, which can equalize the strengths of the output sidebands, reduce them to the level of the input signal and control their degree of correlation (as required). The noise figure of a balanced (zero-net-gain) copier is 4.8 dB [Eq. (115)].

A cascaded PS amplifier is a PI amplifier (which copies the signal as described above), followed by an optical processor (which controls the relative phase of the sidebands) and a PS amplifier (which combines the sidebands). Like its constituent amplifiers, this composite amplifier produces correlated sidebands [Eqs. (75) and (76) or Eqs. (77) and (78)]. The noise figure of a cascaded PS amplifier is 4 dB [Eq. (119)], which is 0.8-dB lower than that of a copier. However, it is 1-dB higher than that of a PI amplifier, because the second amplifier cannot compensate completely the noise added by the first amplifier and the processor (connecting fibers and splices). PS signal processing is obtained at only a moderate (noise) cost.

Two-mode PS links are sequences of transmission fibers (attenuators) followed by optical processors and two-mode PS amplifiers. An idealized link has two CS input sidebands. As these sidebands propagate through the link, they are periodically diminished and decorrelated by the attenuators, then augmented and re-correlated by the amplifiers [Eqs. (90) and (91)]. The noise figure of this PS link is proportional to its total loss, but is 6-dB lower than that of the corresponding PI link [Eqs. (121) and (123)]. Furthermore, the output sidebands are (almost) completely correlated (because the last element in the link is an amplifier), so one can eliminate the effects of electronic noise by detecting both sidebands. A realistic link requires only one nonzero input, because a copier placed before the first attenuator provides the required second input. Analyses show that the presence of the copier and the way in which it is configured (equalized or symmetric) have only minor effects on the noise properties of the link [Eqs. (98)–(100), and Eqs. (104) and (105)], because the (correlated) copier noise is swamped by the (uncorrelated) attenuator noise in the link. Hence, the predicted 6-dB noise advantage of a two-mode PS link is realizable.

In this paper, the number variances and correlations (moments) produced by parametric devices were described in detail. The relation between the amplitude and number moments is described in Appendix A, which also contains more examples of correlations affecting the performances of parametric devices. The main results of this paper were obtained by retaining the signal-noise contributions to the variances and correlations, and omitting the noise-noise contributions, which are usually smaller. For completeness, the latter contributions are calculated in Appendix B.

Appendices

Appendix A: Effects of correlations

If the noise-only contributions to the number variances and correlations are neglected, these second-order number moments depend only on the second-order amplitude moments (which include the amplitude correlations). Input amplitude correlations (produced by prior parametric processes) affect the output number moments and, hence, the noise properties of the current parametric process.

Let a be a mode operator, which satisfies the commutation relation [a,a ^†] = 1, and define the mode amplitude α = 〈a〉, where 〈〉 is an expectation value. Then the noise operator ν = a − 〈a〉 satisfies the commutation relation [ν,ν ^†] = 1 and has zero mean. The number operator m = a ^† a, which can be rewritten in terms of the noise operator as

m = {∣ α ∣}^{2} + α^{*} v + α v^{†} + v^{†} v .

(125)

Hence, the number mean

〈 m 〉 = {∣ α ∣}^{2} + 〈 v^{†} v 〉 .

(126)

The squared number

〈 m^{2} 〉 \approx {∣ α ∣}^{4} + 2 {∣ α ∣}^{2} 〈 v^{†} v 〉 + [{(α^{*})}^{2} 〈 v^{2} 〉 + {∣ α ∣}^{2} 〈 v^{†} v 〉 + {∣ α ∣}^{2} 〈 v v^{†} 〉 + α^{2} 〈 {(v^{†})}^{2} 〉],

(127)

where noise-only terms were neglected. The number variance 〈δm ²〉 = 〈m ²〉 − 〈m〉². By combining Eqs. (126) and (127), one finds that the number variance

〈 δ m^{2} 〉 \approx {(α^{*})}^{2} 〈 v^{2} 〉 + {∣ α ∣}^{2} 〈 v^{†} v 〉 + {∣ α ∣}^{2} 〈 v v^{†} 〉 + α^{2} 〈 {(v^{†})}^{2} 〉 .

(128)

The number variance depends on all four of the second-order self-moments, two of which are independent. (The fourth moment is the conjugate of the first, and the third differs from the second by 1.) If the mode is associated with a CS, the noise operator is a VS operator, so 〈δm ²〉 = ∣α∣² = 〈m〉, which is the standard result.

The applications considered in this paper involve two (or more) modes (sidebands), which are called the signal (1) and idler (2). Equations (125)–(128) apply to each mode separately. By combining the signal and idler versions of Eq. (125), one finds that the number product

〈 m_{1} m_{2} 〉 \approx {∣ α_{1} α_{2} ∣}^{2} + {∣ α_{1} ∣}^{2} 〈 v_{2}^{†} v_{2} 〉 + {∣ α_{2} ∣}^{2} 〈 v_{1}^{†} v_{1} 〉

+ α_{1}^{*} α_{2}^{*} 〈 v_{1} v_{2} 〉 + α_{1}^{*} α_{2} 〈 v_{1} v_{2}^{†} 〉 + α_{1} α_{2}^{*} 〈 v_{1}^{†} v_{2} 〉 + α_{1} α_{2} 〈 v_{1}^{†} v_{2}^{†} 〉 .

(129)

The number correlation 〈δm ₁ δm ₂〉 = 〈m ₁ m ₂〉 − 〈m ₁〉〈m ₂〉. By combining the signal and idler versions of Eqs. (126) with Eq. (129), one finds that the number correlation

〈 δ m_{1} δ m_{2} 〉 \approx α_{1}^{*} α_{2}^{*} 〈 v_{1} v_{2} 〉 + α_{1}^{*} α_{2} 〈 v_{1} v_{2}^{†} 〉 + α_{1} α_{2}^{*} 〈 v_{1}^{†} v_{2} 〉 + α_{1} α_{2} 〈 v_{1}^{†} v_{2}^{†} 〉 .

(130)

The number correlation depends on all four of the second-order cross-moments (correlations), two of which are independent. (The fourth moment is the conjugate of the first, and the third is the conjugate of the second, because ν ₁ and ν ₂ commute.) If the sidebands are independent CS, 〈δm ₁ δm ₂〉 = 0.

As the modes propagate through a device, the amplitudes α_i and noise operators ν_i evolve in the same way as the mode operators a_i . (Specific IO equations for a_i were stated in Sec. 3.) Equations (128) and (130) apply wherever one chooses to measure the sideband variances and correlation. Henceforth, the approximation signs will be omitted.

To illustrate how correlations develop in practice, some examples will be considered. The operation of a two-mode amplifier is governed by Eqs. (27) and (28). By combining these equations, one obtains the second-order IO equations

w_{1}^{2} = μ^{2} v_{1}^{2} + 2 μ ν v_{1} v_{2}^{†} + ν^{2} {(v_{2}^{†})}^{2},

(131)

w_{1}^{†} w_{1} = {∣ μ ∣}^{2} v_{1}^{†} v_{1} + μ^{*} ν v_{1}^{†} v_{2}^{†} + μ ν^{*} v_{1} v_{2} + {∣ ν ∣}^{2} v_{2} v_{2}^{†},

(132)

w_{1} w_{2} = μ^{2} v_{1} v_{2} + μ ν v_{1} v_{1}^{†} + μ ν v_{2}^{†} v_{2} + ν^{2} v_{1}^{†} v_{2}^{†},

(133)

w_{1} w_{2}^{†} = {∣ μ ∣}^{2} v_{1} v_{2}^{†} + μ ν^{*} v_{1}^{2} + μ^{*} ν {(v_{2}^{†})}^{2} + {∣ ν ∣}^{2} v_{1} v_{2}^{†},

(134)

where ν_i and w_i are input and output noise operators, respectively. Each of the output moments depends on input self- and cross-moments, so input correlations affect the outputs and, hence, the noise figures of two-mode amplifiers. Notice that the six moment equations decouple into two sets of three equations. The first set involves w ² ₁, w ² ₂ and w ₁ w ^† ₂, whereas the second involves w ^† ₁ w ₁, w ^† ₂ w ₂ and w ₁ w ₂. Equation (132) and its idler counterpart imply that w ^† ₁ w ₁ − w ^† ₂ w ₂ = ν ^† ₁ ν ₁ − ν ^† ₂ ν ₂. Regardless of the input conditions, noise photons are produced in pairs [19

J. M. Manley and H. E. Rowe, “Some general properties of nonlinear elements—Part I. General energy relations,” Proc. IRE 44, 904–913 (1956). [CrossRef]

, 20

M. T. Weiss, “Quantum derivation of energy relations analogous to those for nonlinear reactances,” Proc. IRE 45, 1012–1013 (1957).

If the inputs are CS, the first set of output moments are zero, whereas the other output moments

〈 w_{i}^{†} w_{i} 〉 = {∣ ν ∣}^{2}, 〈 w_{1} w_{2} 〉 = μ ν .

(135)

By combining Eqs. (128), (130) and (135), one obtains the number variances and correlation

〈 δ n_{i}^{2} 〉 = {∣ β_{i} ∣}^{2} ({∣ μ ∣}^{2} + {∣ ν ∣}^{2}),

(136)

〈 δ n_{1} δ n_{2} 〉 = β_{1}^{*} β_{2}^{*} μ ν + β_{1} β_{2} μ^{*} ν^{*},

(137)

respectively, where the β_i are output amplitudes. Equations (136) and (137) are consistent with Eqs. (31)–(33).

If two amplifiers (also labeled 1 and 2) are concatenated, the composite amplifier is governed by IO equations of the forms (27) and (28), in which the composite transfer coefficients $\bar{μ}$ = μ ₂ μ ₁ + ν ₂ ν ^* ₁ and $\bar{ν}$ = μ ₂ ν ₁ + ν ₂ μ ^* ₁. These coefficients have the squared moduli

{∣ \bar{μ} ∣}^{2} = {∣ μ_{2} μ_{1} ∣}^{2} + 2 Re (μ_{2} ν_{2}^{*} μ_{1} ν_{1}) + {∣ ν_{2} ν_{1} ∣}^{2},

(138)

{∣ \bar{ν} ∣}^{2} = {∣ μ_{2} ν_{1} ∣}^{2} + 2 Re (μ_{2} ν_{2}^{*} μ_{1} ν_{1}) + {∣ ν_{2} μ_{1} ∣}^{2},

(139)

which depend on the relative phase

θ = ϕ_{μ_{2}} - ϕ_{ν_{2}} + ϕ_{μ_{1}} + ϕ_{ν_{1}}

, and satisfy the auxiliary equation ∣ $\bar{μ}$ ∣² − ∣ $\bar{ν}$ ∣² = 1. One can control the relative phase by varying the pump phases or imposing phase shifts on the sidebands between the amplifiers. If the amplifiers combine constructively (θ = 0), ∣ $\bar{μ}$ ∣ = ∣μ ₂ μ ₁∣ + ∣ν ₂ ν ₁∣ and ∣ $\bar{ν}$ ∣ = ∣μ ₂ ν ₁∣ + ∣ν ₂ μ ₁∣, whereas if they combine destructively (θ = π), ∣ $\bar{μ}$ ∣ = ∣μ ₂ μ ₁∣ − ∣ν ₂ ν ₁∣ and ∣ $\bar{ν}$ ∣ = ∣μ ₂ ν ₁∣ − ∣ν ₂ μ ₁∣. The noise figure of the composite amplifier is determined solely by the transmission coefficient ∣ $\bar{ν}$ ∣ (or ∣ $\bar{μ}$ ∣). Suppose that the amplifiers have equal and high PI gains. If the amplifiers combine constructively, ∣ $\bar{μ}$ ∣ = ∣μ∣² + ∣ν∣² ≫ 1 and ∣ $\bar{ν}$ ∣ = 2∣μν∣ ≫ 1: The noise figures of the first amplifier and the composite amplifier both are 3 dB. This result implies that the effective noise figure of the second amplifier is 0 dB, which is not possible if the inputs to the second amplifier are uncorrelated. (The effective noise figure is the input SNR divided by the output SNR, without the proviso that the input is a CS.) Conversely, if the amplifiers combine destructively, ∣ $\bar{μ}$ ∣ = 1 and ∣ $\bar{ν}$ ∣ = 0: The noise added by the first amplifier is removed by the second amplifier, which also is not possible if the inputs to the second amplifier are uncorrelated.

According to Eqs. (135), the nonzero moments produced by the first amplifier are ∣ν ₁∣² and μ ₁ ν ₁. By combining these results with Eq. (132), one finds that the nonzero output moment

〈 w_{i}^{†} w_{i} 〉 = {∣ μ_{2} ∣}^{2} {∣ v_{1} ∣}^{2} + 2 Re (μ_{2} v_{2}^{*} μ_{1} v_{1}) + {∣ v_{2} ∣}^{2} {∣ μ_{1} ∣}^{2},

(140)

which equals ∣ $\bar{ν}$ ∣² [Eq. (139)]. Thus, to predict correctly the properties of the composite amplifier, one is required to account for the correlation between the inputs to the second amplifier. If this correlation were absent, Eqs. (128) and (132) would underestimate the output variance in the constructive regime (by a factor of 2) and overestimate it in the destructive regime (because both nonzero contributions would be positive).

The operation of a two-mode attenuator (frequency convertor) is governed by Eqs. (41) and (42). By combining these equations, one obtains the second-order IO equations

w_{1}^{2} = τ^{2} v_{1}^{2} + 2 τ ρ v_{1} v_{3} + ρ^{2} v_{3}^{2},

(141)

w_{1}^{†} w_{1} = {∣ τ ∣}^{2} v_{1}^{†} v_{1} + τ^{*} ρ v_{1}^{†} v_{3} + τ ρ^{*} v_{1} v_{3}^{†} + {∣ ρ ∣}^{2} v_{3}^{†} v_{3},

(142)

w_{1} w_{3} = - τ ρ^{*} v_{1}^{2} + ({∣ τ ∣}^{2} - {∣ ρ ∣}^{2}) v_{1} v_{3} + τ^{*} ρ v_{3}^{2},

(143)

w_{1} w_{3}^{†} = - τ ρ v_{1} v_{1}^{†} + τ^{2} v_{1} v_{3}^{†} - ρ^{2} v_{1}^{†} v_{3} + τ ρ v_{3} v_{3}^{†},

(144)

where mode 3 is the loss mode (idler). Once again, each of the output moments depends on the input self- and cross-moments, so input correlations affect the outputs and, hence, the noise figures of two-mode attenuators (frequency convertors). Notice that the six moment equations decouple into two sets of three equations. The first set involves w ² ₁, w ² ₃ and w ₁ w ₃, whereas the second involves w ^† ₁ w ₁, w ^† ₃ w ₃ and w ₁ w ^† ₃. (This decomposition differs from the previous one.) Equation (142) and its idler counterpart imply that w ^† ₁ w ₁ + w ^† ₃ w ₃ = ν ^† ₁ ν ₁ + ν ^† ₃ v ₃. Regardless of the input conditions, the total number of noise photons is conserved [19

J. M. Manley and H. E. Rowe, “Some general properties of nonlinear elements—Part I. General energy relations,” Proc. IRE 44, 904–913 (1956). [CrossRef]

, 20

M. T. Weiss, “Quantum derivation of energy relations analogous to those for nonlinear reactances,” Proc. IRE 45, 1012–1013 (1957).

If both inputs are CS, all six normally-ordered output moments are zero. The only nonzero moments are the anti-normally-ordered moments 〈w_iw ^† _i〉 = 1. By combining these results with Eqs. (128) and (130), one obtains the number variances and correlation

〈 δ n_{i}^{2} 〉 = {∣ β_{i} ∣}^{2},

(145)

〈 δ n_{1} δ n_{3} 〉 = 0,

(146)

respectively. Equations (145) and (146) reflect the fact that input CS are converted into different output CS. They are consistent with Eqs. (45) and (46).

If the input signal is an amplified CS (and the input loss mode is a VS), the nonzero input moment 〈ν ^† ₁ ν ₁〉 = ∣ν∣². By combining this result with Eqs. (141)–(144), one obtains the nonzero output moments

〈 w_{1}^{†} w_{1} 〉 = {∣ τ v ∣}^{2}, 〈 w_{3}^{†} w_{3} 〉 = {∣ ρ v ∣}^{2}, 〈 w_{1} w_{3}^{†} 〉 = - τ ρ {∣ v ∣}^{2} .

(147)

In contrast to the previous case, the outputs contain noise photons and are correlated. By combining Eqs. (128), (130) and (147), one obtains the number variances and correlation

〈 δ n_{1}^{2} 〉 = ({∣ τ μ ∣}^{2} + {∣ τ v ∣}^{2} + {∣ ρ ∣}^{2}) {∣ τ μ α_{1} ∣}^{2},

(148)

〈 δ n_{3}^{2} 〉 = ({∣ ρμ ∣}^{2} + {∣ ρv ∣}^{2} + {∣ τ ∣}^{2}) {∣ ρμ α_{1} ∣}^{2},

(149)

〈 δ n_{1} δ n_{3} 〉 = 2 {∣ τ ρ μ v α_{1} ∣}^{2},

(150)

respectively. By combining Eqs. (148)–(150), one finds that

〈 {(δ n_{1} + δ n_{3})}^{2} 〉 = ({∣ μ ∣}^{2} + {∣ ν ∣}^{2}) {∣ μ α_{1} ∣}^{2} .

(151)

The right side of Eq. (151) is the variance of the input signal: Because the attenuator preserves the number of photons, the total output variance must equal the input variance, which would not be possible if the outputs were uncorrelated. (The variance of the input loss mode is zero.)

Now consider a parametric amplifier followed by a frequency convertor that operates on the same sidebands (1 and 2). By combining Eqs. (27) and (28) with Eqs. (41) and (42), one obtains the composite IO equations

b_{1} = (τ μ) a_{1} + (ρ ν) a_{1}^{†} + (ρ μ) a_{2} + (τ ν) a_{2}^{†},

(152)

b_{2} = - (ρ^{*} μ) a_{1} + (τ^{*} ν) a_{1}^{†} + (τ^{*} μ) a_{2} - (ρ^{*} ν) a_{2}^{†} .

(153)

For simplicity, suppose that τ = 1/2^1/2 = ρ. Then Eqs. (152) and (153) can be rewritten in the simple forms

b_{1} = μ a_{+} + ν a_{+}^{†},

(154)

b_{2} = - μ a_{-} + ν a_{-}^{†},

(155)

respectively, where the superposition modes a _± = (a ₁ ± a ₂)/2^1/2. If the input signal and idler are independent CS with amplitudes α ₁ and 0, respectively, the superposition modes are independent CS with (common) amplitude α ₁/2^1/2. Equations (154) and (155) show that the outputs of the composite device are independent squeezed CS with amplitudes (∓μα ₁ + να ^* ₁)/2^1/2 [26

C. J. McKinstrie, S. Radic, M. G. Raymer, and M. V. Vasilyev, “Quantum mechanics of phase-sensitive amplification in fibers,” Opt. Commun. 257, 146–163 (2006). [CrossRef]

]: The frequency convertor decorrelates the amplifier outputs. Squeezed CS are discussed in [6–8

C. J. McKinstrie, M. Yu, M. G. Raymer, and S. Radic, “Quantum noise properties of parametric processes,” Opt. Express 13, 4986–5012 (2005). [CrossRef] [PubMed]

]. Their output numbers and variances depend on the phases of the input amplitudes and transfer coefficients.

As stated previously, the nonzero moments produced by the amplifier are 〈ν ^† _i ν_i 〉 = ∣ν∣² and 〈ν ₁ ν ₂〉 = μν. By combining these results with Eqs. (141)–(144), one finds that the nonzero moments produced by the composite device are

〈 w_{i}^{2} 〉 = \pm μ ν, 〈 w_{i}^{†} w_{i} 〉 = {∣ ν ∣}^{2},

(156)

where the + and − signs apply to modes 1 and 2, respectively. The frequency convertor destroys the cross-moment 〈ν ₁ ν ₂〉 and creates the self-moments 〈w ² _i〉, which would not be possible if the inputs were uncorrelated. It is these self-moments that make the noise properties of the device PS. By combining Eqs. (128), (130) and (156), one obtains the number variances and correlation

〈 δ n_{i}^{2} 〉 = {∣ \pm β^{*} μ + β ν^{*} ∣}^{2},

(157)

〈 δ n_{1} δ n_{2} 〉 = 0 .

(158)

respectively. Equation (157) is consistent with Eq. (16) and Eq. (158) confirms the statement that the outputs are uncorrelated.

Appendix B: Noise-noise contributions

The results of this paper were obtained by retaining the signal-noise contributions to the quantities of interest and omitting the noise-noise contributions. It was shown in [6

C. J. McKinstrie, M. Yu, M. G. Raymer, and S. Radic, “Quantum noise properties of parametric processes,” Opt. Express 13, 4986–5012 (2005). [CrossRef] [PubMed]

] and Sec. 2.2 that the noise-noise contributions to the output variances and correlations are

{〈 δ n_{i}^{2} 〉}_{nn} = 2 Σ_{k} {∣ μ_{i k} ν_{i k} ∣}^{2} + Σ_{k} Σ_{l > k} {∣ μ_{i k} ν_{i l}^{*} + μ_{i l} ν_{i k}^{*} ∣}^{2},

(159)

{〈 δ n_{i} δ n_{j} 〉}_{nn} = 2 Σ_{k} (μ_{i k} ν_{i k}^{*}) (μ_{j k}^{*} ν_{j k}) + Σ_{k} Σ_{l > k} (μ_{i k} ν_{i l}^{*} + μ_{i l} ν_{i k}^{*}) (μ_{j k}^{*} ν_{j l} + μ_{j l}^{*} ν_{j k}),

(160)

respectively. If i = j, Eq. (160) reduces to Eq. (159), as it must do. For the applications considered herein, μ_ik and ν_ik are not nonzero simultaneously. For such applications, Eqs. (159) and (160) can be simplified. The first terms on the right sides are zero. Consider the second term on the right side of Eq. (159). If k is like i (both odd or even) and l is unlike i (one odd and the other even), the contribution to the summation is ∣μ_ikν_il ∣² with k < l. Conversely, if k is unlike and l is like, the contribution is ∣μ_ilν_ik ∣² with l > k, which is equivalent to the first type of contribution with k > l. Hence, the output variance

〈 δ n_{i}^{2} 〉 = (Σ_{k} {∣ μ_{i k} ∣}^{2}) (Σ_{l} {∣ ν_{i l} ∣}^{2}),

(161)

where k is like i, l is unlike i and the subscript nn was omitted. The right side of Eq. (161) is manifestly real. Now consider the second term on the right side of Eq. (160). If i and j are like, the contributions are μ_ikν ^* _il μ ^* _jk ν_jl = μ_ikμ ^* _jkν ^* _ilν_jl with k < l or μ_ilν ^* _ikμ ^* _jlν_jk = μ_ilμ ^* _jlν ^* _ikν_jk with l > k. Hence, the output correlation

〈 δ n_{i} δ n_{j} 〉 = (Σ_{k} μ_{i k} μ_{j k}^{*}) (Σ_{l} ν_{i l}^{*} ν_{j l}),

(162)

where k is like and l is unlike. Equation (3) ensures that the right side of Eq. (162) is real. If i = j, Eq. (162) reduces to Eq. (161), as it must do. Conversely, if i and j are unlike, the contributions are μ_ikν ^* _ilμ ^* _jlν_jk = μ_ikν_jkμ ^* _jlν ^* _il with k < l or μ_ilν ^* _ikμ ^* _jkν_jl = μ_ilν_jlμ ^* _jkν ^* _ik with l > k. Hence, the output correlation

〈 δ n_{i} δ n_{j} 〉 = (Σ_{k} μ_{i k} ν_{j k}) (Σ_{l} μ_{j l}^{*} ν_{i l}^{*}),

(163)

where k is like i and l is unlike i. Equation (2) ensures that the right side of Eq. (163) is real.

It is instructive to consider some examples. Two-mode amplification is governed by Eqs. (27) and (28), in which the signal operator a ₁ is coupled to the conjugate of the idler operator a ^† ₂. The nonzero transfer coefficients μ ₁₁ = μ = μ ₂₂ and ν ₁₂ = ν = ν ₂₁, where ∣μ∣² = G is the PI gain and ∣ν∣² = G − 1. It follows from these facts, and Eqs. (161) and (163), that

〈 δ n_{i}^{2} 〉 = {∣ μ ν ∣}^{2} = 〈 δ n_{1} δ n_{2} 〉 .

(164)

Because the amplitude fluctuations associated with a CS are the same as those associated with a VS, the noise-only contributions to the sideband variances are equal. The correlation equals the (common) variance because sideband photons are created in pairs (by the destruction of pump photons).

Two-mode attenuation (frequency conversion) is governed by Eqs. (41) and (42), in which the signal operator a ₁ is coupled to the loss-mode (idler) operator a ₃. The nonzero transfer coefficients μ ₁₁ = τ = μ ^* ₃₃ and μ ₁₃ = ρ = − μ ^* ₃₁. It follows from these facts, and Eqs. (161) and (162), that

〈 δ n_{i}^{2} 〉 = 0 = 〈 δ n_{1} δ n_{3} 〉 .

(165)

These results reflect the fact that attenuation (frequency conversion) is a stable process. Unlike amplification, which is an unstable process, there is no mechanism to convert input vacuum fluctuations into output number fluctuations. (In frequency conversion, signal and pump photons are destroyed, and idler and other pump photons are created.)

Now consider four-mode processes, of which amplification followed by attenuation (Sec. 3.3), attenuation followed by amplification (Sec. 3.4) and cascaded PS amplification (Sec. 3.5) are examples. It follows from Eqs. (161)–(163) that

〈 δ n_{1}^{2} 〉 = ({∣ μ_{11} ∣}^{2} + {∣ μ_{13} ∣}^{2}) ({∣ ν_{12} ∣}^{2} + {∣ ν_{14} ∣}^{2}),

(166)

〈 δ n_{2}^{2} 〉 = ({∣ μ_{22} ∣}^{2} + {∣ μ_{24} ∣}^{2}) ({∣ ν_{21} ∣}^{2} + {∣ ν_{23} ∣}^{2}),

(167)

〈 δ n_{1} δ n_{2} 〉 = (μ_{11} ν_{21} + μ_{13} ν_{23}) (μ_{22}^{*} ν_{12}^{*} + μ_{24}^{*} ν_{14}^{*}),

(168)

〈 δ n_{1} δ n_{3} 〉 = (μ_{11} μ_{31}^{*} + μ_{13} μ_{33}^{*}) (ν_{12}^{*} ν_{32} + ν_{14}^{*} ν_{34}) .

(169)

Equations (166)–(169) are stated explicitly because they apply to all of the aforementioned examples. It is easy to verify that they are consistent with Eqs. (159) and (160).

Amplification followed by an attenuation is governed by Eqs. (48) and (49), in which ν ₁₄ = 0 = ν ₂₃. It follows from these facts and Eqs. (166)–(168) that

〈 δ n_{1}^{2} 〉 = [T_{1} G + (1 - T_{1})] T_{1} (G - 1),

(170)

〈 δ n_{2}^{2} 〉 = [T_{2} G + (1 - T_{2})] T_{2} (G - 1),

(171)

〈 δ n_{1} δ n_{2} 〉 = T_{1} T_{2} G (G - 1),

(172)

where G is the PI gain, and T ₁ and T ₂ are the signal and idler transmissions, respectively. Equations (170) and (171) are consistent with Eq. (52) of [18

C. J. McKinstrie, S. Radic, R. M. Jopson, and A. R. Chraplyvy, “Quantum noise limits on optical monitoring with parametric devices,” Opt. Commun. 259, 309–320 (2006). [CrossRef]

], which were obtained by a different method. As the transmissions T_j tend to zero, the noise-only contributions to the variances tend to zero as the first power of T_j , whereas the noise-only contributions to the correlation tends to zero as the second power. Hence, attenuation decorrelates the sideband fluctuations completely.

Attenuation followed by amplification is governed by Eqs. (59) and (60). It follows from Eqs. (166)–(168) that

〈 δ n_{i}^{2} 〉 = G (G - 1) = 〈 δ n_{1} δ n_{2} 〉 .

(173)

Equations (173) are identical to Eqs. (164) because attenuators replace CS and VS by diminished CS and identical VS, respectively, all of which have the same amplitude fluctuations, so prior attenuation does not change the noise-only contributions to the output variances and correlation.

Cascaded PS amplification is governed by Eqs. (69) and (70), from which it follows that

〈 δ n_{1}^{2} 〉 = [T H_{0} + (1 - T) G_{2}] [T (H_{0} - 1) + (1 - T) (G_{2} - 1)],

(174)

〈 δ n_{2}^{2} 〉 = [T H_{0} + (1 - T) G_{2}] [T (H_{0} - 1) + (1 - T) (G_{2} - 1)],

(175)

〈 δ n_{1} δ n_{2} 〉 = {T {[H_{0} (H_{0} - 1)]}^{\frac{1}{2}} + (1 - T) {[G_{2} (G_{2} - 1)]}^{\frac{1}{2}}}^{2} .

(176)

In Eqs. (174)–(176), the T terms represent fluctuations that were transmitted through both amplifiers (H ₀), whereas the 1 − T terms represent fluctuations that were added by the attenuators and transmitted through the second amplifier (G ₂).

A noteworthy feature of the preceding equations is the occurrence of the signal-like terms H ₀ and G ₂, the idler-like terms H ₀ − 1 and G ₂ − 1 and the (symmetric) product terms [H ₀(H ₀ − 1)]^1/2 and [G ₂(G ₂ − 1)]^1/2. One can explain why these terms occur by writing the sideband IO equations in the compact forms

b_{1} = c_{1 o} + c_{1 e}^{†}, b_{2} = c_{2 o}^{†} + c_{2 e},

(177)

where the effective-input-mode operators c _1o = μ ₁₁ a ₁ + μ ₁₃ a ₃, c _1e = ν ^* ₁₂ a ₂ + ν * ₁₄ a ₄, c _2o = ν ^* ₂₁ a ₁ + ν * ₂₃ a ₃ and c _2e = μ ₂₂ a ₂ + μ ₂₄ a ₄. These effective-mode operators satisfy the commutation relations listed in Table 1. As an example of how to read the table, the entry ∣μ ₁₁∣² + ∣μ ₁₃∣² in the first row and first column is the value of [c _1o,c ^† _1o]. All commutators of the form [x,y] are zero. The effective-mode operators also have the property that [x,y ^†] = 〈xy ^†〉, where the expectation value is associated with a (four-mode) vacuum state.

Table 1. Commutators for effective modes 1 and 2

[x,y ^†]	c _1o	c _1e	c _2o	c _2e
c _1o	∣μ ₁₁∣² + ∣μ ₁₃∣²	0	μ ₁₁ ν ₂₁ + μ ₁₃ ν ₂₃	0
c _1e	0	∣ν ₁₂∣² + ∣ν ₁₄∣²	0	ν ^* ₁₂ μ ^* ₂₂ + ν ^* ₁₄ μ ^* ₂₄
c _2o	ν ^* ₂₁ μ ^* ₁₁ + ν ^* ₂₃ μ ^* ₁₃	0	∣ν ₂₁∣² + ∣ν ₂₃∣²	0
c _2e	0	μ ₂₂ ν ₁₂ + μ ₂₄ ν ₁₄	0	∣μ ₂₂∣² + ∣μ ₂₄∣²

Equation (177) converts a four-mode process into a two-effective-mode process, the properties of which are well known. By proceeding in the standard way, one finds that the signal number 〈n ₁〉 = 〈c _1e c ^† _1e〉 and the squared number

〈 n_{1}^{2} 〉 = 〈 c_{1 o} c_{1 e} c_{1 o}^{†} c_{1 e}^{†} + c_{1 o} c_{1 e} c_{1 e} c_{1 e}^{†} + c_{1 e} c_{1 e}^{†} c_{1 o}^{†} c_{1 e}^{†} + c_{1 e} c_{1 e}^{†} c_{1 e} c_{1 e}^{†} 〉

= 〈 c_{1 o} c_{1 o}^{†} 〉 〈 c_{1 e} c_{1 e}^{†} 〉 + {〈 c_{1 e} c_{1 e}^{†} 〉}^{2} .

(178)

Terms that are obviously zero (because they involve 〈0∣c ^† or c∣0〉) were omitted from the first version of Eq. (178). One obtains the second version from the first by applying the effective-mode commutation relations described above. By combining the formulas for the first and second powers of the signal number, one finds that the signal variance

〈 δ n_{1}^{2} 〉 = 〈 c_{1 o} c_{1 o}^{†} 〉 〈 c_{1 e} c_{1 e}^{†} 〉 .

(179)

Equation (179) is consistent with Eq. (166). The IO equations (177) are related by the interchanges 1 ↔ 2 and o ↔ e. Hence, the idler number 〈n ₂〉 = 〈c _2o c ^† _2o〉 and the idler variance

〈 δ n_{2}^{2} 〉 = 〈 c_{2 o} c_{2 o}^{†} 〉 〈 c_{2 e} c_{2 e}^{†} 〉 .

(180)

Equation (180) is consistent with Eq. (167). By proceeding in a similar way, one finds that the number product

〈 n_{1} n_{2} 〉 = 〈 c_{1 e} c_{1 o} c_{2 e}^{†} c_{2 o}^{†} + c_{1 e} c_{1 o} c_{2 o} c_{2 o}^{†} + c_{1 e} c_{1 e}^{†} c_{2 e}^{†} c_{2 o}^{†} + c_{1 e} c_{1 e}^{†} c_{2 o} c_{2 o}^{†} 〉

= 〈 c_{1 o} c_{2 o}^{†} 〉 〈 c_{1 e} c_{2 e}^{†} 〉 + 〈 c_{1 e} c_{1 e}^{†} 〉 〈 c_{2 o} c_{2 o}^{†} 〉,

(181)

from which it follows that the correlation

〈 δ n_{1} δ n_{2} 〉 = 〈 c_{1 o} c_{2 o}^{†} 〉 〈 c_{1 e} c_{2 e}^{†} 〉 .

(182)

Equation (182) is consistent with Eq. (168).

In the context of a cascaded PS amplifier, the only modes of interest are the unlike modes 1 and 2, so the preceding analysis is sufficient. However, in other applications, such as four-mode amplification [5

C. J. McKinstrie, S. Radic, and M. G. Raymer, “Quantum noise properties of parametric amplifiers driven by two pump waves,” Opt. Express 12, 5037–5066 (2004). [CrossRef] [PubMed]

], the properties of modes 3 and 4 also of interest. The IO equation for mode 3 is related to the first of Eqs. (177) by the interchange 1 ↔ 3, so the variance of mode 3 and the correlation between modes 2 and 3 follow from Eqs. (179) and (182), respectively. For the like modes 1 and 3, a separate calculation is required. By following the procedure described above, and using the commutation relations listed in Table 2, one finds that the number product

〈 n_{1} n_{3} 〉 = 〈 c_{1 e} c_{1 o} c_{3 o}^{†} c_{3 e}^{†} + c_{1 e} c_{1 o} c_{3 e} c_{3 e}^{†} + c_{1 e} c_{1 e}^{†} c_{3 o}^{†} c_{3 e}^{†} + c_{1 e} c_{1 e}^{†} c_{3 e} c_{3 e}^{†} 〉

= 〈 c_{1 o} c_{3 o}^{†} 〉 〈 c_{1 e} c_{3 e}^{†} 〉 + 〈 c_{1 e} c_{1 e}^{†} 〉 〈 c_{3 e} c_{3 e}^{†} 〉 .

(183)

Table 2. Commutators for effective modes 1 and 3

[x,y ^†]	c _1o	c _1e	c _3o	c _3e
c _1o	∣μ ₁₁∣² + ∣μ ₁₃∣²	0	μ ₁₁ ^μ ^* ₃₁ + μ ₁₃ μ ^* ₃₃	0
c _1e	0	∣ν ₁₂∣2 + ∣ν ₁₄∣²	0	ν ^* ₁₂ ν ₃₂ + ν ^* ₁₄ ν ₃₄
c _3o	μ ₃₁ μ ^* ₁₁ + μ ₃₃ μ ^* ₁₃	0	∣μ ₃₁∣² + ∣μ ₃₃∣²	0
c _3e	0	ν ^* ₃₂ ν ₁₂ + ν ^* ₃₄ ν ₁₄	0	∣ν ₃₂∣² + ∣ν ₃₄∣²

Hence, the correlation

〈 δ n_{1} δ n_{3} 〉 = 〈 c_{1 o} c_{3 o}^{†} 〉 〈 c_{1 e} c_{3 e}^{†} 〉 .

(184)

Equation (184) is consistent with Eq. (169). One can deduce the properties of mode 4 from those of mode 2. Thus, the variances and correlations produced by four-mode (and multiple-mode) processes are determined by the commutators of the odd and even effective-mode operators, as are the variances and correlations produced by the prototypical two-mode processes of amplification (Sec. 3.1) and frequency conversion (Sec. 3.2).

Acknowledgments

CJM acknowledges useful discussions with S. Radic and M. G. Raymer.

References and links

1.	J. Hansryd, P. A. Andrekson, M. Westland, J. Li, and P. O. Hedekvist, “Fiber-based optical parametric amplifiers and their applications,” IEEE J. Sel. Top. Quantum Electron. 8, 506–520 (2002). [CrossRef]
2.	S. Radic and C. J. McKinstrie, “Optical amplification and signal processing in highly nonlinear optical fiber,” IEICE Trans. Electron. E88-C, 859–869 (2005). [CrossRef]
3.	P. A. Andrekson and M. Westlund, “Nonlinear optical fiber based all-optical waveform sampling,” Laser Photon. Rev. 1, 231–248 (2007). [CrossRef]
4.	S. Radic, “Parametric amplification and processing in optical fibers,” Laser Photon. Rev. 2, 489–513 (2008). [CrossRef]
5.	C. J. McKinstrie, S. Radic, and M. G. Raymer, “Quantum noise properties of parametric amplifiers driven by two pump waves,” Opt. Express 12, 5037–5066 (2004). [CrossRef] [PubMed]
6.	C. J. McKinstrie, M. Yu, M. G. Raymer, and S. Radic, “Quantum noise properties of parametric processes,” Opt. Express 13, 4986–5012 (2005). [CrossRef] [PubMed]
7.	R. Loudon and P. L. Knight, “Squeezed light,” J. Mod. Opt. 34, 709–759 (1987). [CrossRef]
8.	S. M. Barnett and P. M. Radmore, Methods in Theoretical Quantum Optics (Oxford University Press, 1997).
9.	N. Christensen, R. Leonhardt, and J. D. Harvey, “Noise characteristics of cross-phase modulation instability light,” Opt. Commun. 101, 205–212 (1993). [CrossRef]
10.	J. E. Sharping, M. Fiorentino, and P. Kumar, “Observation of twin-beam-type quantum correlation in optical fiber,” Opt. Lett. 26, 367–369 (2001). [CrossRef]
11.	R. Tang, P. Devgan, P. L. Voss, V. S. Grigoryan, and P. Kumar, “In-line frequency-nondegenerate phase-sensitive fiber-optical parametric amplifier,” IEEE Photon. Technol. Lett. 17, 1845–1847 (2005). [CrossRef]
12.	R. Tang, J. Lasri, P. S. Devgan, V. Grigoryan, and P. Kumar, “Gain characteristics of a frequency nondegenerate phase-sensitive fiber-optic parametric amplifier with phase self-stabilized input,” Opt. Express 13, 10483–10493 (2005). [CrossRef] [PubMed]
13.	C. Lundström, J. Kakande, P. A. Andrekson, Z. Tong, M. Karlsson, P. Petropoulos, F. Parmigiani, and D. J. Richardson, “Experimental comparison of gain and saturation characteristics of a parametric amplifier in phase-sensitive and phase-insensitive mode,” ECOC 2009, paper 1.1.1.
14.	J. Kakande, C. Lundström, P. A. Andrekson, Z. Tong, M. Karlsson, P. Petropoulos, F. Parmigiani, and D. J. Richardson, “Detailed characterization of a fiber-optic parametric amplifier in phase-sensitive and phase-insensitive operation,” Opt. Express 18, 4130–4137 (2010). [CrossRef] [PubMed]
15.	R. Loudon, “Theory of noise accumulation in linear optical-amplifier chains,” J. Quantum Electron. 21, 766–773 (1985). [CrossRef]
16.	R. E. Slusher and B. Yurke, “Squeezed light for coherent communications,” J. Lightwave Technol. 8, 466–477 (1990). [CrossRef]
17.	M. Vasilyev, “Distributed phase-sensitive amplification,” Opt. Express 13, 7563–7571 (2005). [CrossRef] [PubMed]
18.	C. J. McKinstrie, S. Radic, R. M. Jopson, and A. R. Chraplyvy, “Quantum noise limits on optical monitoring with parametric devices,” Opt. Commun. 259, 309–320 (2006). [CrossRef]
19.	J. M. Manley and H. E. Rowe, “Some general properties of nonlinear elements—Part I. General energy relations,” Proc. IRE 44, 904–913 (1956). [CrossRef]
20.	M. T. Weiss, “Quantum derivation of energy relations analogous to those for nonlinear reactances,” Proc. IRE 45, 1012–1013 (1957).
21.	C. J. McKinstrie, J. D. Harvey, S. Radic, and M. G. Raymer, “Translation of quantum states by four-wave mixing in fibers,” Opt. Express 13, 9131–9142 (2005). [CrossRef] [PubMed]
22.	M. G. Raymer, S. J. van Enk, C. J. McKinstrie, and H. J. McGuinness, “Interference of two photons of different color,” Opt. Commun. 283, 747–752 (2010). [CrossRef]
23.	Z. Tong, A. Bogris, C. Lundström, C. J. McKinstrie, M. Vasilyev, M. Karlsson, and P. A. Andrekson, “Modeling and measurement of the noise figure of a cascaded non-degenerate phase-sensitive parametric amplifier,” Opt. Express 18, 14820–14835 (2010). [CrossRef] [PubMed]
24.	Z. Tong, C. J. McKinstrie, C. Lundström, M. Karlsson, and P. A. Andrekson, “Noise performance of optical fiber transmission links that use non-degenerate cascaded phase-sensative amplifiers,” Opt. Express 18, 15426–15439 (2010). [CrossRef] [PubMed]
25.	L. A. Krivitsky, U. L. Andersen, R. Dong, A. Huck, C. Wittmann, and G. Leuchs, “Electronic noise-free measurements of squeezed light,” Opt. Lett. 33, 2395–2397 (2008). [CrossRef] [PubMed]
26.	C. J. McKinstrie, S. Radic, M. G. Raymer, and M. V. Vasilyev, “Quantum mechanics of phase-sensitive amplification in fibers,” Opt. Commun. 257, 146–163 (2006). [CrossRef]

OCIS Codes
(060.2320) Fiber optics and optical communications : Fiber optics amplifiers and oscillators
(190.4380) Nonlinear optics : Nonlinear optics, four-wave mixing
(270.2500) Quantum optics : Fluctuations, relaxations, and noise

ToC Category:
Quantum Optics

History
Original Manuscript: June 9, 2010
Revised Manuscript: August 26, 2010
Manuscript Accepted: August 28, 2010
Published: September 1, 2010

Citation
C. J. McKinstrie, M. Karlsson, and Z. Tong, "Field-quadrature and photon-number correlations produced by parametric processes," Opt. Express 18, 19792-19823 (2010)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-18-19-19792

Sort: Author | Year | Journal | Reset

References

J. Hansryd, P. A. Andrekson, M. Westland, J. Li, and P. O. Hedekvist, “Fiber-based optical parametric amplifiers and their applications,” IEEE J. Sel. Top. Quantum Electron. 8, 506–520 (2002). [CrossRef]
S. Radic, and C. J. McKinstrie, “Optical amplification and signal processing in highly nonlinear optical fiber,” IEICE Trans. Electron,” E 88-C, 859–869 (2005). [CrossRef]
P. A. Andrekson, and M. Westlund, “Nonlinear optical fiber based all-optical waveform sampling,” Laser Photon. Rev. 1, 231–248 (2007). [CrossRef]
S. Radic, “Parametric amplification and processing in optical fibers,” Laser Photon. Rev. 2, 489–513 (2008). [CrossRef]
C. J. McKinstrie, S. Radic, and M. G. Raymer, “Quantum noise properties of parametric amplifiers driven by two pump waves,” Opt. Express 12, 5037–5066 (2004). [CrossRef] [PubMed]
C. J. McKinstrie, M. Yu, M. G. Raymer, and S. Radic, “Quantum noise properties of parametric processes,” Opt. Express 13, 4986–5012 (2005). [CrossRef] [PubMed]
R. Loudon, and P. L. Knight, “Squeezed light,” J. Mod. Opt. 34, 709–759 (1987). [CrossRef]
S. M. Barnett, and P. M. Radmore, Methods in Theoretical Quantum Optics (Oxford University Press, 1997).
N. Christensen, R. Leonhardt, and J. D. Harvey, “Noise characteristics of cross-phase modulation instability light,” Opt. Commun. 101, 205–212 (1993). [CrossRef]
J. E. Sharping, M. Fiorentino, and P. Kumar, “Observation of twin-beam-type quantum correlation in optical fiber,” Opt. Lett. 26, 367–369 (2001). [CrossRef]
R. Tang, P. Devgan, P. L. Voss, V. S. Grigoryan, and P. Kumar, “In-line frequency-nondegenerate phase-sensitive fiber-optical parametric amplifier,” IEEE Photon. Technol. Lett. 17, 1845–1847 (2005). [CrossRef]
R. Tang, J. Lasri, P. S. Devgan, V. Grigoryan, and P. Kumar, “Gain characteristics of a frequency nondegenerate phase-sensitive fiber-optic parametric amplifier with phase self-stabilized input,” Opt. Express 13, 10483–10493 (2005). [CrossRef] [PubMed]
C. Lundström, J. Kakande, P. A. Andrekson, Z. Tong, M. Karlsson, P. Petropoulos, F. Parmigiani and D. J. Richardson, “Experimental comparison of gain and saturation characteristics of a parametric amplifier in phase-sensitive and phase-insensitive mode,” ECOC 2009, paper 1.1.1.
J. Kakande, C. Lundström, P. A. Andrekson, Z. Tong, M. Karlsson, P. Petropoulos, F. Parmigiani, and D. J. Richardson, “Detailed characterization of a fiber-optic parametric amplifier in phase-sensitive and phase-insensitive operation,” Opt. Express 18, 4130–4137 (2010). [CrossRef] [PubMed]
R. Loudon, “Theory of noise accumulation in linear optical-amplifier chains,” J. Quantum Electron. 21, 766–773 (1985). [CrossRef]
R. E. Slusher, and B. Yurke, “Squeezed light for coherent communications,” J. Lightwave Technol. 8, 466–477 (1990). [CrossRef]
M. Vasilyev, “Distributed phase-sensitive amplification,” Opt. Express 13, 7563–7571 (2005). [CrossRef] [PubMed]
C. J. McKinstrie, S. Radic, R. M. Jopson, and A. R. Chraplyvy, “Quantum noise limits on optical monitoring with parametric devices,” Opt. Commun. 259, 309–320 (2006). [CrossRef]
J. M. Manley, and H. E. Rowe, “Some general properties of nonlinear elements–Part I. General energy relations,” Proc. IRE 44, 904–913 (1956). [CrossRef]
M. T. Weiss, “Quantum derivation of energy relations analogous to those for nonlinear reactances,” Proc. IRE 45, 1012–1013 (1957).
C. J. McKinstrie, J. D. Harvey, S. Radic, and M. G. Raymer, “Translation of quantum states by four-wave mixing in fibers,” Opt. Express 13, 9131–9142 (2005). [CrossRef] [PubMed]
M. G. Raymer, S. J. van Enk, C. J. McKinstrie, and H. J. McGuinness, “Interference of two photons of different color,” Opt. Commun. 283, 747–752 (2010). [CrossRef]
Z. Tong, A. Bogris, C. Lundström, C. J. McKinstrie, M. Vasilyev, M. Karlsson, and P. A. Andrekson, “Modeling and measurement of the noise figure of a cascaded non-degenerate phase-sensitive parametric amplifier,” Opt. Express 18, 14820–14835 (2010). [CrossRef] [PubMed]
Z. Tong, C. J. McKinstrie, C. Lundström, M. Karlsson, and P. A. Andrekson, “Noise performance of optical fiber transmission links that use non-degenerate cascaded phase-sensative amplifiers,” Opt. Express 18, 15426–15439 (2010). [CrossRef] [PubMed]
L. A. Krivitsky, U. L. Andersen, R. Dong, A. Huck, C. Wittmann, and G. Leuchs, “Electronic noise-free measurements of squeezed light,” Opt. Lett. 33, 2395–2397 (2008). [CrossRef] [PubMed]
C. J. McKinstrie, S. Radic, M. G. Raymer, and M. V. Vasilyev, “Quantum mechanics of phase-sensitive amplification in fibers,” Opt. Commun. 257, 146–163 (2006). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper

Figures


Fig. 1.	Fig. 2.	Fig. 3.


Fig. 4.	Fig. 5.	Fig. 6.


Fig. 7.	Fig. 8.	Fig. 9.


Fig. 10.	Fig. 11.	Fig. 12.


Fig. 13.	Fig. 14.	Fig. 15.