## Abstract

Noise is one of the basic characteristics of optical amplifiers. Whereas there are various noise sources, the intrinsic one is quantum noise that originates from Heisenberg’s uncertainty principle. This chapter describes quantum noise in optical amplifiers, including population-inversion–based amplifiers such as an Erbium-doped fiber amplifier and a semiconductor optical amplifier, and optical parametric amplifiers. A full quantum mechanical treatment is developed based on Heisenberg equation of motion for quantum mechanical operators. The results provide the quantum mechanical basis for a classical picture of amplifier noise widely used in the optical communication field.

### Keywords

- quantum noise
- noise figure
- quantum mechanics
- population-inversion–based amplifier
- optical parametric amplifier

## 1. Introduction

Noise is one of the important properties in optical amplifiers [1]. The intrinsic noise characteristic is determined by quantum mechanics, especially Heisenberg’s uncertainty principle. This chapter describes quantum noise in optical amplifiers in terms of quantum mechanics. After brief introduction of a classical treatment usually used in the optical communication field, properties of an optically amplified light, such as the mean amplitude, the mean photon number, and their variances, are derived based on first principles of quantum mechanics. Two kinds of optical amplifiers are treated: amplifiers based on two-level interaction in a population-inverted medium, i.e., an Erbium-doped fiber amplifier and a semiconductor optical amplifier, and those based on parametric interaction in an optical nonlinear medium. The results presented here provide the quantum mechanical basis to a phenomenological classical treatment conventionally employed for describing amplifier noise.

## 2. Classical treatment

A classical treatment of amplifier noise is widely employed in the optical communication field [1, 2], whereas it originates from quantum mechanics. Before presenting a quantum mechanical treatment, we briefly introduce the classical treatment. We first consider the light intensity or the photon number outputted from an amplifier. A photon-number rate equation for light propagating through a population-inverted medium can be expressed as

where * n* is the number of photons,

N

_{2}and

N

_{1}are the numbers of atoms at the higher and lower energy states in the medium, respectively, and

*is a constant representing the photon emission/absorption efficiency. The first, second, and third terms represent stimulated emission, absorption, and spontaneous emission, respectively. The efficiency*g

*is common in this phenomenon in a simple two-level model [3]. Assuming that*g

N

_{1}and

N

_{2}are uniform along the medium length, the photon number at the output is calculated from Eq. (1) as

where * L* is the medium length. The first term represents amplified signal photons with a signal gain of exp[

*(*g

N

_{2}–

N

_{1})

*] ≡*L

*. The second term represents amplified spontaneous emission (ASE) photons, which can be rewritten as*G

n

_{sp}(

*– 1) with*G

n

_{sp}≡

N

_{2}/(

N

_{2}–

N

_{1}). The parameter

n

_{sp}is called “population inversion parameter” or “noise factor.” Note that the above equations are for the photon number in one mode in terms of the frequency and the polarization.

Eq. (2) shows that the output photon number is composed of amplified signal photons and ASE photons. Accordingly, the output amplitude is supposed to be a summation of amplitudes of amplified signal and ASE lights as

where _{out}, _{in}, and _{ASE} are the amplitudes of the output light, the input signal light, and ASE light, respectively. The second term provides the ASE power as <|_{ASE}|^{2}> = _{sp}(* G* – 1)

*Δ*hf

*, where < > denotes the mean value,*f

*is one photon energy, and Δ*hf

*is the ASE bandwidth. Regarding the ASE phase, it is supposed to be completely random because spontaneous emission occurs randomly. Thus, the average of the ASE amplitude is supposed to be zero: <*f

E

_{ASE}> = 0. Here, we decompose

E

_{ASE}into the real and imaginary parts, which are supposed to be isotropic because the phase is random: <{Re[

E

_{ASE}]}

^{2}> = <{Im[

E

_{ASE}]}

^{2}> = <|

E

_{ASE}|

^{2}>/2 =

n

_{sp}(

*–1)*G

*Δ*hf

*/2. Subsequently, the variance of each amplitude component is <{Re[*f

E

_{ASE}]}

^{2}> − <Re[

E

_{ASE}]>

^{2}= <{Im[

E

_{ASE}]}

^{2}> − <Im[

E

_{ASE}]>

^{2}=

n

_{sp}(

*– 1)*G

*Δ*hf

*/2.*f

Intensity noise is evaluated using Eq. (3). The output intensity is given by _{out} = |_{out}|^{2}, and its fluctuation is evaluated by the variance of _{out} as

where the postulate of the ASE light phase being random is used in averaging. The first term represents 2 × <signal output intensity> × <ASE intensity>, which is called the “signal-spontaneous beat noise.” The second and third terms represent the intensity variance of the ASE light, which is called the “spontaneous-spontaneous beat noise.”

As an indicator for the amplifier noise performance, the “noise figure (NF)” is usually used. It is defined as the ratio of the signal-to-noise ratios (SNRs) at the input and output of an amplifier in terms of the optical intensity: NF ≡ (SNR)_{in}/(SNR)_{out} where SNR ≡ (mean intensity)^{2}/(variance of the intensity) in the signal mode. The square of the mean intensity at the output is calculated from Eq. (3) as <|_{out}|^{2}> = * G*|

E

_{in}|

^{2}+ <|

E

_{ASE}|

^{2}>, and the output variance is expressed as Eq. (4); thus the output SNR is expressed as

On the other hand, the input SNR is evaluated for pure monochromatic light in the definition of the noise figure. In quantum mechanics, such a light is called “coherent state,” whose photon-number variance is equal to the mean photon number: <^{2}> − <* n*>

^{2}= <

*>. Thus, <|*n

E

_{in}|

^{4}> − <|

E

_{in}|

^{2}>

^{2}= |

E

_{in}|

^{2}

*where*hf

*is attached for the dimension to be matched. Subsequently, the input SNR is given by <|*hf

E

_{in}|

^{2}>

^{2}/(<|

E

_{in}|

^{4}> − <|

E

_{in}|

^{2}>

^{2}) = |

E

_{in}|

^{2}/

*. Therefore, the NF is expressed as*hf

where * G*|

E

_{in}|

^{2}> > |

E

_{ASE}|

^{2}is assumed, and <|

E

_{ASE}|

^{2}> =

n

_{sp}(

*– 1)*G

*is substituted because the signal mode is considered here.*hf

The above-mentioned classical treatment is widely used for noise in optical amplifiers. However, it is based on phenomenological assumptions. (i) Eq. (3) is phenomenologically provided. Though the solution of the photon-number rate equation indicates that the output photon number is composed of amplified signal photons and ASE photons (Eq. (2)), this result does not logically conclude that the output amplitude is a linear summation of the amplified signal and the ASE amplitudes as Eq. (3). (ii) The phase of ASE light is assumed to be random, which is a phenomenological postulate, not logically derived from first principles. Although the above classical treatment is correct and useful in fact, we need quantum mechanics for theoretically confirming its validity, which is presented in the following sections.

## 3. Quantum mechanics

In this section, we briefly review quantum mechanics, especially the Heisenberg picture [4]. The basic concept of quantum mechanics is that a physical state is probabilistic and the theory only provides mean values of physical quantities, which is given by a quantum mechanical inner product of a physical quantity operator

where _{0}> is an initial state.

The most important operator in discussing quantum mechanical properties of light is the “annihilation operator,” * α*> and an eigenvalue

*:*α

*> is called “coherent state” and corresponds to pure monochromatic light in classical optics, and the eigenvalue*α

*corresponds to its complex amplitude in the photon-number unit. The Hermitian conjugate of the annihilation operator*α

*> is a photon number state having*n

*photons. Another important operator is the photon-number operator defined as*n

We discuss quantum noise of optical amplifiers in the following sections, using the above-mentioned framework of quantum mechanics. Note that the above operator

## 4. Population-inversion–based amplifiers

Erbium-doped fiber amplifiers (EDFAs) are widely used in optical communications. Optical semiconductor amplifiers are also being developed for compact and integrated amplifying devices. They amplify signal light through interaction between light and a two-level atomic system with population inversion. This section discusses quantum noise in population-inversion-based amplifiers [5].

### 4.1. Heisenberg equation

The Hamiltonian for a light-atom interacting system can be expressed as [4]

The first and second terms are the Hamiltonians of light and atoms without interaction, respectively, where * ω* is the angular frequency of light;

*and |1>*

_{j}

*denoting the upper and lower energy states, respectively, satisfying <1|1>*

_{j}

_{j,k}= <2|2>

_{j,k}=

δ

_{j,k}and <1|2>

_{j,k}= <2|1>

_{j,k}= 0;

*indicates a specific atom. The third term is the interaction Hamiltonian between light and atoms, which represents energy exchange such that a photon is created while an atom transits from the upper to the lower states and vice versa, with*j

*being the coupling coefficient.*κ

_{j}

Applying the above Hamiltonian to the Heisenberg equations for

Employing the variable translations

where Δ* ω* ≡

_{j}

*–*ω

ω

_{a}

^{(j)}, i.e., the frequency detuning between light and a two-level system.

We solve Eq. (10) by an iterative approximation. First, the first-order solutions are derived by substituting the initial values {

The solutions of these equations are

Next, these first-order solutions are substituted into the right-hand side of Eq. (10a), and the second-order solution is calculated as

We regard Eq. (13) as the time evolution of the field operator during a short time * τ*, and rewrite it as

where

Eq. (14) is the basic expression for discussing quantum properties of light that travels through an amplifier. For the discussion, we also need an initial state of the system at _{0}. It can be expressed as

where |* Ψ*(

t

_{0})> denotes the initial state of light, |Ψ

_{a}(

t

_{0})> denotes that of atoms, and

c

_{1}and

c

_{2}are the probability amplitudes of an atom being in the lower and upper states, respectively, satisfying |

c

_{1}|

^{2}+ |

c

_{2}|

^{2}= 1. We use Eqs. (14) and (16) in the following calculations.

### 4.2. Mean amplitude

We first discuss the mean amplitude. The mean amplitude, denoted as

The average of the transition operator

Assuming that the energy levels of each atom are densely distributed, the summation in this equation can be replaced by an integral in the frequency domain, and the real part of Eq. (18) is further calculated as

where Ω is the frequency detuning; * ρ* is the density of atoms in the frequency domain and is assumed to be constant around a resonant frequency Ω = 0 as

ρ

_{0}; {

*, |*κ

c

_{2}|

^{2}, and |

c

_{1}|

^{2}} are assumed to be identical for any atom; and

N

_{1}≡

ρ

_{0}|

c

_{1}|

^{2}and

N

_{2}≡

ρ

_{0}|

c

_{2}|

^{2}are the numbers of atoms at the lower and upper energy states, respectively. On the other hand, the imaginary part of Eq. (18) can be rewritten as

In this expression, the contents of the integral is an odd function around the resonant frequency Ω = 0. Thus, the imaginary part equals 0. Regarding the average of _{1}^{*}_{2}) = 0 on average. Subsequently,

_{j}

Eq. (21) describes the time evolution of the mean amplitude of light traveling through an amplifying medium, i.e., the time evolution in a frame moving along with the light, during a short time. This expression can be translated to the spatial evolution along the medium length as

where * g* ≡ 2

κ

^{2}

*/*π

*with*v

*being the light velocity. Applying a Taylor expansion*v

*(*x

z

_{0}+ Δ

*) =*z

*(*x

z

_{0}) + [

*/*dx

*](*dz

z

_{0}) × Δ

*to this equation, we obtain the following differential equation:*z

Eq. (23) includes (_{2} – _{1}), which depends on the atoms’ state at a local position. Here, we assume that the atoms’ state is uniform along the medium length independent of * z*, and this condition is satisfied in uniformly pumped amplifiers with no gain saturation. With this assumption, Eq. (23) can be analytically solved, and the mean amplitude at the amplifier output is expressed as

where * G* ≡ exp[

*(*g

N

_{2}–

N

_{1})

*] with*L

*being the amplifier length. This result, derived from the Heisenberg equation, is equivalent to the classical expression, i.e., Eq. (3) with <*L

E

_{ASE}> = 0, and confirms that the mean amplitude of ASE light is zero.

### 4.3. Mean photon number

We next discuss the mean photon number. The short-time evolution of the photon-number operator is expressed from Eq. (14) as

from which the short-time evolution of the mean photon number is obtained as

In deriving Eq. (26), higher-order interaction terms are neglected, because the short-time evolution is considered here. This short-time evolution is translated to the short-length evolution along the medium length as

from which the following spatial differential equation is obtained:

This equation is equivalent to the photon-number rate equation given in Eq. (1). Therefore, similar to Eq. (1), the output of the mean photon number is calculated as

The first and second terms represent amplified signal photons and ASE photons, respectively. It is noted that ASE photons appear at the output even though there is no such light in the mean amplitude as shown in Eq. (24).

### 4.4. Amplitude fluctuation

We next discuss amplitude fluctuations or noise. The light amplitude has two quadratures, i.e., the real and imaginary components. The operators representing each component are

From Eq. (14), the short-time evolution of the mean square of the real component is expressed as

where * O*(

τ

^{2}) terms are neglected. On the other hand, the square of the average of the real component is expressed as

From Eqs. (30) and (31), the short-time evolution of the variance of the real component is obtained as

This equation is translated to the short-length evolution as

from which the following differential equation is obtained:

From Eq. (34), the variance of the real component at the output is calculated as

The variance of the imaginary component is similarly calculated as _{x2}^{2}(* L*) =

Gσ

_{x2}

^{2}(0) + (1/4)(2

n

_{sp}– 1)(

*– 1). The first term in Eq. (35) represents amplified fluctuations from the incident light, with a gain*G

*whose square root equals the amplitude gain [Eq. (24)], and the second term represents additional fluctuations that are superimposed onto the amplified fluctuation through the amplification process. This input and output relationship of amplitude fluctuations can be schematically illustrated in the complex amplitude space (constellation) as shown in Figure 1.*G

For a coherent incident state, i.e., whose amplitude variances is _{x1}^{2} (0) = _{x2}^{2} (0) = 1/4 [4], Eq. (35) is rewritten as

The first term 1/4 corresponds to the inherent quantum noise of a coherent state, and the second term represents amplitude fluctuations at the amplifier output in a classical picture. Recalling that the mean amplitude at the amplifier output is that amplified from the incident light with no addition mean field, as indicated in Eq. (24), Eq. (36) suggests that the amplifier output can be regarded as a summation of a clean signal light (i.e., coherent state), displaced from the initial mean amplitude position, and fluctuating light, whose mean value and variance are 0 and _{sp}(* G* – 1)/2, respectively, in one quadrature. Figure 2 illustrates this output condition in the complex amplitude space. Noting that the variance of the fluctuating light equals half of the spontaneous photon number indicated in the second term in Eq. (28), we can say that the above picture illustrated in Figure 2 is equivalent to the classical picture of amplitude noise described in Section 2, where the ASE power is given by <|

E

_{ASE}|

^{2}> =

n

_{sp}(

*– 1)*G

*Δ*hf

*and the variance of the real component of ASE light is given by <{Re[*f

E

_{ASE}]}

^{2}> − <Re[

E

_{ASE}]> =

n

_{sp}(

*– 1)*G

*Δ*hf

*/2. Therefore, the classical picture introduced in Section 2 is confirmed by the quantum mechanical treatment presented here, except for the inherent quantum noise of 1/4. This noise 1/4 is sometimes called “vacuum fluctuation” or “zero-point fluctuation,” that appears owing to quantum mechanics.*f

### 4.5. Photon-number fluctuation

We next discuss photon-number fluctuations. These fluctuations are evaluated employing the variance of the photon number as

where higher-order interaction terms are neglected as before. This expression can be translated to the short-length evolution as

from which the following differential equation is obtained:

with * g*(

N

_{2}–

N

_{1}), and the second and third terms represent the number of photons generated at a local position, which propagate and reach the medium end while being amplified by the first term. Then, the solution of Eq. (39) can be expressed as

Here,

From this equation and Eq. (29), the photon-number variance at the amplifier output is expressed as

Recalling that the mean photon numbers of the amplified signal and the spontaneous emission are ^{2}, corresponding to the spontaneous-spontaneous beat noise represented by the second and third terms in Eq. (4). The third and fourth terms denote the mean photon numbers of the amplified signal and spontaneous emission, respectively, corresponding to the inherent quantum noises of the amplified signal light and the spontaneous emission, respectively. In the fifth term,

The first and second terms in Eq. (42) correspond to the classical intensity noise represented by Eq. (4), as described above, supporting the classical treatment. In addition, the inherent quantum noises are included in Eq. (42), owing to the full quantum mechanical treatment, and the amplified excess noise is simultaneously included as well. Sometimes in the classical treatment, the inherent quantum noise is phenomenologically added as the shot noise arising at the electrical stage after direct detection [2]. In fact, however, it exists in the optical stage as derived above. Therefore, the inherent quantum noise is sometimes called “optical shot noise.”

### 4.6. Noise figure

The noise figure, defined as the ratio of the signal-to-noise ratios (SNRs) at the input and output of an amplifier in terms of the light intensity or the photon number, is usually used as an indicator for the noise performance of an amplifier. Based on the above results, we describe the noise figure of population-inversion-based amplifiers in this subsection. The output SNR is obtained from Eqs. (29) and (42) as

where only the signal power and the signal-spontaneous beat noise are taken into account, assuming that the amplified signal is sufficiently larger than the spontaneous emission. On the other hand, the input SNR is evaluated for a coherent state, according to the definition of the noise figure, which is

This expression equals the classical result given by Eq. (6).

The noise figure is proportional to the population inversion parameter _{sp} = _{2}/(_{2} – _{1}), as shown above. The minimum value of _{sp} is 1, which is achieved when _{1} = 0, i.e., the fully inverted condition where all atoms are in the upper states. Under this condition, NF = 3 dB for * G* >> 1. This is the quantum-limited noise figure of population-inversion-based optical amplifiers. Near-quantum-limited noise figure has been demonstrated experimentally in Erbium-doped fiber amplifiers [6, 7].

The fact that the noise performance is determined by the population inversion parameter can be intuitively understood as follows. The source of amplifier noise is spontaneous emission. A small amount of spontaneous emission suggests a good noise performance. However, spontaneous emission is roughly proportional to the signal gain (Eq. (29)), which is desired to be high as an amplifier. Thus, the amount of spontaneous emission normalized to the signal gain, (ASE power)/(signal gain), can be an indicator for the noise performance. The spontaneous emission rate is proportional to the number of atoms in the upper energy level _{2}, i.e., (ASE power)∝_{2}, and the signal gain is determined by the balance between stimulated emission and absorption and thus is proportional to the difference between the numbers of atoms in the upper and lower states, roughly speaking, i.e., (signal gain)∝(_{2} – _{1}). Subsequently, (ASE power)/(signal gain)∝_{2}/(_{2} – _{1}) = _{sp}, which suggests that the amplifier noise performance is determined by the population inversion parameter _{sp}.

## 5. Optical parametric amplifiers

Whereas population-inversion–based amplifiers are widely used, there is another type of optical amplifiers, that is an optical parametric amplifier (OPA) based on optical nonlinearity [8]. When signal light is incident onto a nonlinear medium along with intense pump light, a signal and idler photons are created from one pump photon in case of second-order nonlinearity, satisfying the energy conservation of ℏ_{s} + ℏ_{i} = ℏ_{p} (_{s}, _{i}, and _{p} are the angular frequencies of the signal, idler, and pump lights, respectively), or they are created from two pump photons in case of third-order nonlinearity, satisfying ℏ_{s} + ℏ_{i} = ℏ_{p1} + ℏ_{p2}. Through this photon exchange phenomenon, the signal light is amplified. This signal amplification scheme also offers optical signal processing functions such as wavelength conversion and generation of phase-conjugated light [9]. This section describes quantum noise in OPAs [10].

### 5.1. Heisenberg equation

The Hamiltonian for parametric interaction between signal and idler via pump light(s) can be expressed as [11]

The first and second terms are the Hamiltonians of signal and idler lights without interaction, respectively, where * χ*. Since the pump light is so intense that its quantum properties do not matter here, the pump light is treated classically, whose amplitude

E

_{p}is included in the coupling coefficient as

*∝*χ

E

_{p}or

E

_{p1}

E

_{p2}for the second- or third-order nonlinear interaction, respectively.

From the Heisenberg equation with the above Hamiltonian, temporal differential equations for the field operators are obtained as

These temporal differential equations can be translated to spatial ones as

where * β* =

*(*n

*/*ω

*) is the propagation constant (*c

*: the refractive index,*n

*: the light velocity in the vacuum). The above equations can be simplified by the variable translation*c

Here, we consider the propagation phase of the right-hand term in the above equations. The coefficient * χ* includes the pump light amplitude as

*∝*χ

E

_{p}or

E

_{p1}

E

_{p2}, and the pump amplitude can be expressed as

E

_{p}=

E

_{p}(0)exp(−

iβ

_{p}

*) under no pump-depletion condition (*z

β

_{p}is the propagation constant of the pump light). Subsequently,

*∝*χ

E

_{p}(0)exp(−

iβ

_{p}

*) or*z

E

_{p1}(0)

E

_{p2}(0) exp[−

*(*i

β

_{p1}+

β

_{p2})

*]. From these considerations, Eq. (48) can be rewritten as*z

where * κ* is the coupling coefficient in the spatial domain, excluding exp(−

iβ

_{p}

*) or exp[−*z

*(*i

β

_{p1}+

β

_{p2})

*], and Δ*z

*≡*β

β

_{s}+

β

_{i}–

β

_{p}or

β

_{s}+

β

_{i}–

β

_{p1}–

β

_{p2}. This parameter Δ

*is called “phase mismatch,” and determines the signal gain of an OPA as shown later. As for*β

*, its absolute value is |*κ

*| =*κ

*|*dγ

E

_{p1}| |

E

_{p2}| when an optical fiber is used as a nonlinear medium, where

*is the nonlinear coefficient, and*γ

*is the degeneracy factor that takes 2 and 1 for*d

f

_{p1}≠

f

_{p2}and

f

_{p1}=

f

_{p2}, respectively. Regarding the phase of

*, it is determined by the incident phase(s) of the pump light(s).*κ

From Eq. (49), the signal field operator at the output is calculated as

where * g* ≡ {|

*|*κ

^{2}– (Δ

*/2)*β

^{2}}

^{1/2},

*is the medium length, and*L

*≡ arg(*φ

*). Eq. (50a) includes the field operators of signal and idler lights, i.e., the signal and the idler are treated separately. For particular frequency conditions such as 2*κ

ω

_{s}=

ω

_{p}for the second-order nonlinearity or 2

ω

_{s}=

ω

_{p1}+

ω

_{p2}for the third-order nonlinearity, the idler frequency equals to the signal frequency,

ω

_{i}=

ω

_{s}, and the signal and the idler are degenerate. Under such conditions, Eq. (50a) is rewritten as

As shown later, degenerate and nondegenerate OPAs have definitely different characteristics. For simplifying mathematical expressions, hereafter, we rewrite Eq. (50) as

with

The mean values of the physical quantities after amplification can be evaluated using Eq. (51). In the evaluation, we need the initial state in addition. Here, we assume that only signal light is incident to an OPA, and express the initial state as * Ψ*>

_{s}and |0>

_{i}denote the signal and idler states, respectively. When the initial state is a coherent state as |

*>*Ψ

_{s}= |

*>, we have*α

*are the mean photon number and the phase of the incident signal light, respectively. On the other hand,*θ

### 5.2. Mean amplitude, photon number, and signal gain

The mean amplitude and photon number at the output are evaluated by

Eq. (53) indicates that the signal field is simply amplified while preserving the phase state, with no additional field on average. On the other hand, Eq. (54) shows that the output photons consist of two components. The first term is proportional to the incident photon number, which corresponds to the amplified signal photons with a gain of

It is noted in this expression that parameter * g* = {|

*|*κ

^{2}– (Δ

*/2)*β

^{2}}

^{1/2}is equivalent to the gain coefficient. When Δ

*= 0,*β

*is maximum and the signal gain is maximum. Therefore, it is important to satisfy the condition Δ*g

*= 0, that is called the “phase matching condition,” in implementing an OPA [9]. The second term in Eq. (54) is independent on the signal input, and represents spontaneously emitted photons. The above results, i.e., the spontaneous light does not appear in the mean amplitude while it does in the photon number, suggest that the amplitude of the spontaneous emission is completely random. However, we do not know how random it is at this stage. The photon number of the spontaneous emission is expressed from the second term in Eqs. (54) and (52c) as*β

where Eq. (55) is applied. This expression is equivalent to the spontaneous photon number in population-inversion-based amplifiers indicated in Eq. (29) with _{sp} = 1. This correspondence suggests that nondegenerate OPAs can offer the ideal noise performance achievable in EDFAs, which is shown later.

Regarding degenerate OPA, on the other hand, its mean output amplitude is calculated as

where _{0} is the phase of the incident signal light. The mean output amplitude does not have a simple form as in nondegenerate OPA (Eq. (53)). Under the condition where Δ* β* = 0 and the gain coefficient

*is so large as cosh(*g

*) ≈ sinh(*gL

*) ≈*gL

*/2, Eq. (57) is approximated as*e

^{gL}

where Δ ≡ _{0} + (* φ* –

*)/2 is introduced. This expression indicates that the phase state of the incident signal light is not transferred to the output.*ϕ

The mean photon number in degenerate OPA is calculated from Eq. (51) as

Unfortunately, this equation cannot be further developed, because we cannot readily calculate * Ψ*>. However, for a coherent incident state, these quantities can be evaluated using

In this expression, the first term represents amplified signal photons, and the second term represents spontaneous emission whose mean amplitude is zero as indicated in Eq. (57).

From the first term in Eq. (60), the signal gain is expressed as

which is dependent on the relative phase Δ = _{0} + (* φ*−

*)/2. Hence, degenerate OPA is called “phase-sensitive amplifier (PSA).” The maximum gain is obtained when Δ = 0 as*ϕ

### 5.3. Amplitude fluctuation

Next, we evaluate the amplitude noise in OPAs. For the evaluation, the light amplitude is decomposed into two quadratures and the variance of each quadrature is calculated, as in Section 4.3. In case of OPAs, the output field operator is phase-shifted by (* φ* +

*)/2, as indicated by Eq. (51). Accordingly, we introduce a phase-shifted field operator defined as*ϕ

_{0}|

_{0}> − <Ψ

_{0}|

_{0}>

^{2}.

The calculation result for nondegenerate OPA is expressed as

where Eqs. (55) and (56) are applied. The first term represents noise amplified from the incident light, and the second term represents additional noise superimposed via OPA. Note that Eq. (63) is equivalent to the amplitude variance of population-inversion-based amplifiers shown by Eq. (35) with _{sp} = 1, suggesting that the ideal noise performance achievable in EDFAs can be obtained in OPA. Similar to Eq. (35), Eq. (63) is rewritten for a coherent incident state as

The first term is the inherent quantum noise of a coherent state, and the second term represents noise superimposed via OPA in a classical picture. The sum of the second terms of the two quadratures equals the photon number of spontaneous emission light indicated in Eq. (56). This consideration supports the classical noise treatment, described in Section 2, where spontaneous emission with random phase is superimposed onto signal light at the amplifier output.

For degenerate OPA, on the other hand, the amplitude variances are calculated as

where _{0} ≡ * G*(Δ = 0) is the phase-synchronized gain introduced in Eq. (62). The results show that the output variances are unequal in the two quadratures, such that

σ

_{x1}

^{2}is enhanced while

σ

_{x2}

^{2}is depressed. Since

*+*φ

*)/2 is enhanced and that along the orthogonal axis is depressed, in the complex amplitude space. It is noted in Eq. (65a) that it only has fluctuations amplified from incident light with no additional fluctuation, unlike population-inversion-based amplifiers (Eq. (35)) and nondegenerate OPA (Eq. (63)). This result suggests that even the noise-enhanced quadrature*ϕ

*) ≈ sinh(*gL

*) ≈*gL

*/2 and then*e

^{gL}

*≈*A

*, is rewritten as*B

σ

_{x1}

^{2}≈ 0. This consideration indicates that the amplitude noise along the phase axis orthogonal to the signal output phase approaches zero as the signal gain increases. The above-mentioned characteristics of amplitude fluctuation in degenerate OPA can be illustrated in the complex amplitude space as shown in Figure 3.

### 5.4. Photon-number fluctuation and noise figure

Next, we discusses photon-number fluctuations in OPAs, which are evaluated through the photon-number variance as _{n}^{2} = <Ψ_{0}|_{0}> − <Ψ_{0}|_{0}>^{2}. Using Eq. (51), the average of the square of the photon-number operator for nondegenerate OPA is calculated as

Subsequently, the photon-number variance is obtained as

where Eqs. (55) and (56) are applied. Recalling that the mean photon number of the amplified signal is * G* – 1), as indicated by Eqs. (54)–(56), each term in Eq. (68) can be interpreted as follows. The first term is equivalent to 2 × (signal light intensity) × (spontaneous light intensity), which corresponds to the signal-spontaneous beat noise. The second term is equivalent to (spontaneous light intensity)

^{2}, which corresponds to the spontaneous-spontaneous beat noise. The third and fourth terms are equal to the mean photon numbers of the amplified signal and of spontaneous emission, respectively, and correspond to the optical shot noises of the amplified signal light and the spontaneous emission, respectively.

Next, we consider degenerate OPA. As indicated by Eq. (59), properties of the photon number in degenerate OPA are hard to evaluate for an arbitrary initial state. Thus, we assume a coherent incident state here. From Eq. (51), the average of the square of the photon-number operator for the initial state |Ψ_{0}> = |* α*>

_{s}⊗ |0>

_{i}is calculated as

Subsequently, the photon-number variance at the output is

where Eq. (61) is applied. This expression cannot be decomposed and interpreted as that of nondegenerate OPA indicated by Eq. (68), which could be because the amplitude distribution is not simply isotropic in two quadratures, unlike nondegenerate OPA.

The noise figure can be evaluated from the results obtained above. For nondegenerate OPA, it is obtained as

where only the signal-spontaneous beat noise is considered for the output SNR, according to the definition on the noise figure. This noise figure equals that of ideal population-inversion–based amplifiers indicated by Eq. (44) with _{sp} = 1. For the degenerate case, on the other hand, it is expressed as

For Δ = 0, NF = 1 (0 dB), suggesting no SNR degradation in phase-synchronized degenerate OPA. In fact, a noise figure of less than 3 dB in a phase-sensitive amplifier has been experimentally demonstrated [12, 13].

## 6. Conclusion

This chapter describes quantum noise of optical amplifiers. Full quantum mechanical treatment based on the Heisenberg equation for physical quantity operators was presented, by which quantum properties of optical amplifiers were derived from first principles. The obtained results are consistent with a conventional classical treatment, except for the inherent quantum noise or the zero-point fluctuation, providing the theoretical base to the conventional phenomenological treatment.

## References

- 1.
Yamamoto Y, Inoue K. Noise in amplifiers. Journal of Lightwave Technology. 2003; 21 :2895-2915. DOI: 10.1109/JLT.2003.816887 - 2.
Olsson N. Lightwave system with optical amplifiers. Journal of Lightwave Technology. 1989; 7 :1071-1082 - 3.
Yariv A. Quantum Electronics. 3rd ed. USA: Wiley; 1989. 676 p - 4.
Loudon R. The Quantum Theory of Light. 3rd ed. New York: Oxford; 2000. 438 p - 5.
Inoue K. Quantum mechanical treatment of optical amplifiers based on population inversion. IEEE Journal of Quantum Electronics. 2014; 50 :563-567. DOI: 10.1109/JQE.2014.2325906 - 6.
Smart RG, Zyskind JL, Sulhoff JW, DiGiovanni DJ. An investigation of the noise figure and conversion efficiency of 0.98 μm pumped erbium-doped fiber amplifiers under saturated conditions. IEEE Photonics Technology Letters. 1992; 4 :1261-1264 - 7.
Laming RI, Zechael MN, Payne DN. Erbium-doped fiber amplifier with 54 dB gain and 3.1 dB noise figure. IEEE Photonics Technology Letters. 1992; 4 :1345-1347 - 8.
Marhic M. Fiber Optical Parametric Amplifiers, Oscillators and Related Devices. 1st ed. New York: Cambridge; 2008. 366 p - 9.
Hansryd J, Andrekson PA, Westlund M, Li J, Hedekvist P. Fiber-based optical parametric amplifiers and their applications. IEEE Journal of Selected Topics in Quantum Electronics. 2002; 8 :506-520 - 10.
Inoue K. Quantum noise in parametric amplification under phase-mismatched conditions. Optics Communication. 2016; 366 :71-76. DOI: 10.1016/j.optcom.2015.12.034 - 11.
Walls D, Milburn G. Quantum Optics. 2nd ed. Berlin: Springer; 2008. 425 p - 12.
Imajuku W, Takada A, Yamabayashi Y. Inline coherent optical amplifier with noise figure lower than 3 dB quantum limit. Electronics Letters. 2000; 36 :63-64 - 13.
Tong Z, Bogris A, Lundström C, McKinstrie CJ, Vasilyev M, Karosson M, Andrekson A. Modeling and measurement of the noise figure of a cascaded non-degenerate phase-sensitive parametric amplifier. Optics Express. 2010; 18 :14820-14835