Advances in OptoElectronics

Advances in OptoElectronics / 2011 / Article

Research Article | Open Access

Volume 2011 |Article ID 368327 |

Niculae N. Puscas, "Modelling of Output Statistics of Single and ℳ -Mode Straight and Curved E r 𝟑 + :Ti:LiNb O 𝟑 Waveguide Amplifiers", Advances in OptoElectronics, vol. 2011, Article ID 368327, 6 pages, 2011.

Modelling of Output Statistics of Single and ℳ -Mode Straight and Curved E r 𝟑 + :Ti:LiNb O 𝟑 Waveguide Amplifiers

Academic Editor: Xian Cao
Received29 Apr 2011
Accepted20 Jul 2011
Published28 Sep 2011


A theoretical analysis of some statistical parameters which characterize the Er3+-doped Ti:LiNbO3 single and ℳ-mode straight and curved waveguides is presented in this paper. In the derivation and the evaluation of the spectral optical quality factor, the power spectral density, the Fano factor, the statistical fluctuation, and the spontaneous emission factor we used the small gain approximation, and the photon statistics master equation of the linear amplifier (considering that the photon number distribution is determined by the normalized mode intensity profiles which are not uniform in the transversal section of the waveguide), transposed to the case of straight and curved amplifiers. The simulation results show the evolution of the above-mentioned parameters under various pump regimes and waveguide lengths.

1. Introduction

Lithium niobate photonic circuits permit the generation, transmission, and processing of photons to be accommodated on a single chip. Compact photonic circuits with multiple components integrated on a single chip are crucial for efficiently implementing quantum information processing schemes.

The study of the noise spectral distribution and the output statistics in a system where coherent and chaotic (thermal) fields are superimposed plays an important role in obtaining integrated amplifiers with low noise and high optical gain. This is the reason why, over the last decade, a great attention has been devoted to the analysis of the statistical properties of Er3+-doped fibres and waveguide amplifiers [1–12].

A novel technique for evaluating the output statistics of single and ℳ-mode straight and curved Er3+-doped Ti:LiNbO3 amplifiers, like the spectral optical quality factor, the power spectral density, the output mean photon number, the statistical fluctuation, the Fano factor, and the spontaneous emission factor, is proposed in this paper. The reason of this analysis is related to the fact that often in the waveguides is excited not only the fundamental mode but also other high-order (ℳ) modes [11, 13], which influence the output gain, noise figure, signal-to-noise ratio, and the statistical properties of the waveguide.

The paper is organized as follows. Section 2 is devoted to the basic equations used to evaluate the population of the upper atomic energy level and the evolution of the pump, signal, power spectral density, Amplified Spontaneous Emission (ASE), the Fano factor, the statistical fluctuation, and the spontaneous emission factor which characterize the output statistics of the single and ℳ-mode waveguides. Also, in this section the model used for the calculation of optical field distribution in the bent waveguides is presented [6–11]. Section 3 deals with the discussion of the computed results, while the conclusion of this work is presented in Section 4.

2. Theoretical Considerations

The quantum theory of the statistical properties in the optical waveguide amplifiers is based on the more general quantum theories of coherent light, coherent light/matter interaction, noise, and laser oscillations [1, 8–11]. We used the two-level model to describe the above-mentioned processes [1, 8–11]. The transitions between the two lowest energetic level manifolds of the Er3+ ions incorporated into the LiNbO3 lattice determine absorption and optical amplification by stimulated emission in the range 1440 nm < 𝜆 < 1640 nm. The optical amplification of a signal wave takes place if the population inversion with respect to the ground state 4I15/2 is achieved after the optical pumping to a pump band (an excited state) followed by a fast relaxation to the metastable level 4I13/2. Using a pump radiation having 𝜆 = 1484 nm it is possible to neglect the excited state absorption (ESA) and to consider a quasi-two-level model for the simulation of the Er3+-doped Ti∶LiNbO3 waveguide amplifiers in the wavelength range 1440 nm < 𝜆 < 1640 nm.

The interaction of pump/signal photons with Er3+ ions can be represented using the rate Equations [8–11]. In the case of the steady state regime by discretizing the considered frequency spectrum in a number of intervals Δ𝜈 small enough so that the frequency-dependent quantities can be assumed constant over each of them, one obtains the equations for the pump, signal, and ASE evolution after integration over the waveguide transversal cross-section in the form [8–11] ğ‘‘ğ‘ƒğ‘‘ğ‘§Â±ğ‘˜ğ‘šî€·ğœŽ(𝑧)=±𝑒,𝑘𝑚+ğœŽğ‘Ž,𝑘𝑚𝐴𝑁2(𝑥,𝑦,𝑧)ğ‘–ğ‘šâˆ’î‚µğœŽ(𝑥,𝑦)ğ‘‘ğ‘¥ğ‘‘ğ‘¦ğ‘Ž,𝑘𝑚𝐴𝑁2(𝑥,𝑦,𝑧)𝑖𝑚(𝑥,𝑦)𝑑𝑥𝑑𝑦+𝛼𝑘𝑚×𝑃±𝑘𝑚(𝑧)+â„Žğœˆğ‘˜Î”ğœˆğœŽğ‘’,𝑘𝑚𝐴𝑁2(𝑥,𝑦,𝑧)𝑖𝑚(,𝑥,𝑦)𝑑𝑥𝑑𝑦(1) where 𝑃±𝑘𝑚(𝑧) represent the power in the spectral interval Δ𝜈 around frequency 𝜈𝑘, for the polarization 𝑚 in the forward (+) and backward (−) propagation directions, with 𝑥 and 𝑦 being the coordinates in the plane perpendicular to the waveguide axis 𝑧 and 𝑁2 is the steady state solution of the rate equations for the upper populations of the two levels which is written in the form𝑁2(𝑥,𝑦,𝑧)=𝑁𝑇0𝑑×∑(𝑥,𝑦)𝑘,𝑚𝜏/â„Žğœˆğ‘˜î€¸ğœŽğ‘Ž,𝑘𝑚𝑃+𝑘𝑚(𝑧)+𝑃−𝑘𝑚𝑖(𝑧)𝑚(𝑥,𝑦)∑1+𝑘,𝑚𝜏/â„Žğœˆğ‘˜ğœŽî€¸î€·ğ‘Ž,𝑘𝑚+ğœŽğ‘’,𝑘𝑚𝑃+𝑘𝑚(𝑧)+𝑃−𝑘𝑚𝑖(𝑧)𝑚.(𝑥,𝑦)(2)

In (2) the Er3+ ions distribution in the waveguide has been written in the form 𝑁𝑇(𝑥,𝑦)=𝑁𝑇0𝑑(𝑥,𝑦), where 𝑁𝑇0=𝑁𝑇(0,0) and 𝑑(𝑥,𝑦) is the normalized dopant distribution. Similar notations have been used for the cross-sections, ğœŽğ‘Ž,𝑚(𝜈) and ğœŽğ‘’,𝑚(𝜈) for the optical waveguide losses 𝛼𝑘𝑚. The term â„Žğœˆğ‘˜Î”ğœˆ in (1) represents an equivalent spontaneous emission input noise power in the frequency slot Δ𝜈 corresponding to the frequency 𝜈𝑘 [1].

Spontaneous transitions from the excited states to the ground state are taken into account by the rate 𝐴21=1/𝜏,𝜏 being the fluorescence lifetime of the Er3+ ions. The absorption and emission rates of the amplified spontaneous emission are 𝑊ASE12and𝑊ASE21, respectively.

The energy fine structure is taken into account by wavelength-dependent absorption and emission cross-sections. The cross-sections profiles constitute the basis for the computation of the gain, noise figure, and other parameters which characterize the amplifiers. The pump (signal) absorption, and emission rates are 𝑅12(𝑊12) and 𝑅21(𝑊21), respectively. To analyse the behaviour of the amplifier in terms of signal gain, pump absorption and ASE output power, it is convenient to distinguish in the optical power at frequency 𝜈𝑘 the coherent contribution due to injected signal or pump beams (stimulated emission or absorption) from the incoherent one due to the (amplified) spontaneous emission. We can do this by observing that the forcing term in (1) contributes only to the incoherent ASE power and by using different “𝑘” indexes for the two kinds of contribution. The stimulated transition rates due to the pump, signal, and ASE fields can be written in terms of the corresponding field intensities 𝐼𝑚(𝑥,𝑦,𝑧,𝜈), the absorption and emission cross-sections, ğœŽğ‘Ž,𝑚(𝜈) and ğœŽğ‘’,𝑚(𝜈), and the photon energies â„Žğœˆ [1, 8]:𝑅12,𝑊12,𝑊ASE12=𝑚=TE,TMî€œğœŽğ‘Ž,𝑚𝐼(𝜈)𝑚(𝑥,𝑦,𝑧,𝜈)ğ‘…â„Žğœˆğ‘‘ğœˆ,21,𝑊21,𝑊ASE21=𝑚=TE,TMî€œğœŽğ‘’,𝑚(𝐼𝜈)𝑚(𝑥,𝑦,𝑧,𝜈)â„Žğœˆğ‘‘ğœˆ.(3) In (1)–(3) the index 𝑚 allows to distinguish the two polarizations sides of the guided field. Moreover, the optical intensities can be expressed as the product of a power spectral density 𝑃𝑚(𝑥,𝑦) and the normalized field intensities 𝑖𝑚(𝑥,𝑦), which can be assumed constant in the frequency range of interest.

The boundary conditions for (1) can be written in the form 𝑃+𝑘𝑚(𝑧=0)=1−𝑅0,𝑘𝑚𝑃𝑖𝑛0,𝑘𝑚+𝑅0,𝑘𝑚𝑃−𝑘𝑚𝑃(𝑧=0),−𝑘𝑚(𝑧=𝐿)=1−𝑅𝐿,𝑘𝑚𝑃𝑖𝑛𝐿,𝑘𝑚+𝑅𝐿,𝑘𝑚𝑃+𝑘𝑚(𝑧=𝐿),(4) where 𝑅0,𝑘𝑚,𝑅𝐿,𝑘𝑚 are the input and output reflectivities and 𝑃𝑖𝑛0,𝑘𝑚,𝑃𝑖𝑛𝐿,𝑘𝑚 are the optical powers injected at 𝑧=0 or at 𝑧=𝐿. We assume that the signal is injected only from the section 𝑧=0 and that no noise power is entering the amplifier (i.e., 𝑃𝑖𝑛0,𝑘𝑚=𝑃𝑖𝑛𝐿,𝑘𝑚=0 for the ASE beams). By properly choosing the values of 𝑅0 and 𝑅𝐿 we can account for the natural reflectivities of the LiNbO3-air interface, as well as analyse the case of a wavelength-selective coating; moreover, if we set 𝑅𝐿 to a high value (approx. 100%),our model can describe the behavior of a double-pass amplifier for the pump and the signal. For the numerical integration of the system of equations (1)-(2) in our model we took into account the overlap between the normalized field intensity and the (normalized) dopant distribution and the (normalized) upper laser level population by introducing the corresponding overlap integrals Γ𝑇,𝑚 and Γ2,𝑚, respectively, Γ𝑇,𝑚(𝑧)=𝐴𝑑(𝑥,𝑦)𝑖𝑚(𝑥,𝑦)𝑑𝑥𝑑𝑦,(5)Γ2,𝑚1(𝑧)=𝑁𝑇0𝐴𝑁2(𝑥,𝑦,𝑧)𝑖𝑚(𝑥,𝑦)𝑑𝑥𝑑𝑦.(6)

The effects of gain saturation due to the amplified signal or to the ASE power are accounted for in the expression for the population 𝑁2: as the signal and ASE power grows the emission rate increases and 𝑁2 decreases, leading to a reduction of Γ2,𝑚 and consequently of the gain factor ğ‘Žğ‘˜ğ‘šÎ“2,𝑚 in (1).

Considering a collection of atoms with population density distribution 𝑁1(𝑥,𝑦,𝑧) for the ground state and 𝑁2(𝑥,𝑦,𝑧) for the excited state, the rate of change of the probability for having 𝑛 photons at frequency 𝜈 in a single longitudinal mode (𝑃𝑛) is given by the photon statistics master equation of the linear amplifier [1].

Assuming that the normalized mode intensity profiles 𝑖𝑚(𝑥,𝑦) are not uniform in the transversal section of the waveguide the photon number distribution becomes 𝑛×𝑖𝑚(𝑥,𝑦). Multiplying the photon statistics master equation by 𝑖𝑚(𝑥,𝑦) and integrating over the transversal section of the waveguide we obtain the following [9]: 𝑑𝑃𝑛(𝑧,𝜈)𝑑𝑧=𝛾𝑒(𝑧,𝜈)𝑛𝑃𝑛−1(𝑧,𝜈)+ğ›¾ğ‘Ž(𝑧,𝜈)𝑛𝑃𝑛+1−𝛾(𝑧,𝜈)𝑒(𝑧,𝜈)(𝑛+1)+ğ›¾ğ‘Žî€»ğ‘ƒ(𝑧,𝜈)𝑛𝑛(𝑧,𝜈),(7) where 𝛾𝑒=𝛾𝑒(𝑧,𝜈)=ğœŽğ‘’(𝜈)𝐴𝑁2(𝑥,𝑦,𝑧)𝑖𝑚(𝛾𝑥,𝑦)𝑑𝑥𝑑𝑦,ğ‘Ž=ğ›¾ğ‘Ž(𝑧,𝜈)=ğœŽğ‘Žî€œ(𝜈)𝐴𝑁1(𝑥,𝑦,𝑧)𝑖𝑚(𝑥,𝑦)𝑑𝑥𝑑𝑦.(8)

In (7) ğœŽğ‘Ž and ğœŽğ‘’ represent the cross-sections for the absorption and emission processes. Furthermore, multiplying (7) by the number of photons, 𝑛, and summing over 𝑛 we obtain the output expression for the photon mean value ⟨𝑛(𝑧)⟩: ⟨𝑛(𝑧)⟩=𝐺(𝑧)⟨𝑛(0)⟩+𝑁(𝑧),(9) where 𝐺(𝑧,𝜈)=exp𝑧0î€ºğ›¾ğ‘’î€·ğ‘§î…žî€¸,ğœˆâˆ’ğ›¾ğ‘Žî€·ğ‘§î…žî€¸î€»,𝜈−𝛼(𝜈)ğ‘‘ğ‘§î…žî‚¼,(10)𝑁(𝑧,𝜈)=𝐺(𝑧,𝜈)𝑧0ğ›¾ğ‘’î€·ğ‘§î…žî€¸,𝜈𝐺(ğ‘§î…žğ‘‘,𝜈)z(11) represent the spectral gain and the ASE photon number, respectively.

For coherent input signals with Poisson’s statistics 𝑃𝑛(0)=(⟨𝑛(0)⟩𝑛/𝑛!)𝑒−⟨𝑛(0)⟩ being well known that ğœŽ2(0)=⟨𝑛(0)⟩. In the case of uniform inversion, when the population densities of the upper and lower levels (𝑁1and 𝑁2, resp.) do not depend on the waveguide coordinate 𝑧, the variance takes the canonical form: ğœŽ2+(𝑧)=𝐺(𝑧)⟨𝑛(0)⟩+𝑁(𝑧)2𝐺(𝑧)𝑁(𝑧)⟨𝑛(0)⟩+𝑁2.(𝑧)(12)

As a measure of the gain and the noise characteristics of the waveguide amplifiers, we adopted a quality factor𝑄(𝑧,𝜈)=𝐺(𝑧,𝜈),𝐹(𝑧,𝜈)(13) where 𝐹(𝑧,𝜈)=1+2𝑁(𝑧,𝜈)𝐺(𝑧,𝜈)(14) represents the noise figure which quantifies the noise properties of an optical amplifier contributing to the deterioration of the signal-to-noise ratio with a purely shot noise.

As mentioned before spontaneous emission is also present in any amplifier. Small amount of this spontaneous emission gets amplified and comes out along with signal as amplified ASE noise.

The ASE noise is generally modelled as white noise with the power spectral density [12]: 𝑆(𝑧,𝜈)=𝑛sp[](𝑧,𝜈)𝐺(𝑧,𝜈)−1â‹…â„Žğœˆ,(15) where 𝑛sp(𝑧,𝜈) represents the spontaneous emission factor and ℎ the Plank's constant. The amplifier output statistics can also be characterised by the Fano factor 𝑓(𝑧)=ğœŽ2(𝑧)/⟨𝑛(𝑧)⟩ and the statistical fluctuation 𝑒(𝑧)=ğœŽ(𝑧)/⟨𝑛(𝑧)⟩ [1]. For Poisson statistics (coherent light), the Fano factor and the statistical fluctuation correspond to 𝑓=1 and 𝑒=(⟨𝑛⟩)−1/2, while for Bose-Einstein statistics (incoherent light), they are given by 𝑓=⟨𝑛⟩+1 and 𝑒=(1+1/⟨𝑛⟩)−1/2.

Assuming an input signal characterised by Poisson statistics, and in the limit of high input signals and high gains, the Fano factor and the statistical fluctuation are given by 𝑓(𝑧,𝜈)≈1+2𝑛sp[],1(𝑧,𝜈)𝐺(𝑧,𝜈)−1𝑒(𝑧,𝜈)≈√𝐺(𝑧,𝜈)⟨𝑛(0)⟩1+2𝑛sp[](𝑧,𝜈)𝐺(𝑧,𝜈)−11/2,(16) where the spontaneous emission factor is given by𝑛sp(𝑧,𝜈)=𝑁(𝑧,𝜈)=𝐺(𝑧,𝜈)−1𝐺(𝑧,𝜈)𝐺(𝑧,𝜈)−1𝑧0ğ›¾ğ‘’î€·ğ‘§î…žî€¸,𝜈𝐺(ğ‘§î…ž,𝜈)ğ‘‘ğ‘§î…ž.(17)

In (16) the factor ∫[1+2𝐺(𝐿,𝜈)𝐿0(𝛾𝑒(ğ‘§î…ž,𝜈)/𝐺(ğ‘§î…ž,𝜈))ğ‘‘ğ‘§î…ž]1/2 characterizes the deviation of the output statistics from Poisson.

Often in the waveguides is excited not only the fundamental mode but also other high-order (ℳ) modes, which influence the output gain, noise figure, and the statistical properties of the waveguide. The normalized field transversal intensity distribution can be written as [13] 𝑖(𝑥,𝑦)=𝑀𝑗=1𝜂𝑗𝑖𝑗𝑖(𝑥,𝑦)exp𝜔𝑡−𝛽𝑗𝑧,(18) where 𝛽𝑗 are the propagation constants which are different for different modes because of their velocity [11] and 𝜂𝑗=∬𝐸𝑝(𝑥,𝑦)⋅𝐸𝑗(𝑥,𝑦)𝑑𝑥𝑑𝑦∬𝐸2𝑝(𝑥,𝑦)𝑑𝑥𝑑𝑦∬𝐸2𝑗(𝑥,𝑦)𝑑𝑥𝑑𝑦(19) are the overlap coefficients between the pump optical field, 𝐸𝑝(𝑥,𝑦), and the mode fields, 𝐸𝑗(𝑥,𝑦), excited in the waveguide, which satisfy a normalization condition ∑ℳ𝑗=1𝜂𝑗=1. In our model we considered that the waveguide is pumped by a radiation having a Gaussian distribution of the field using an optical fiber. The integrals in (19) are extended over the transversal section of the waveguide, with 𝑥 representing the width and 𝑦 the depth, respectively.

In order to determine the optical field distribution in the bent waveguides we calculated first the refractive index profiles using the Fick’s diffusion law [12]. After that, the optical mode fields were calculated numerically for both TE and TM polarisation using the effective index method presented in papers [7, 14].

3. Discussion of the Simulation Results

The system of coupled first-order differential equations (1)-(2) for the optical power components and the upper population level can only be solved by numerical methods. Moreover, the presence of the boundary conditions (4) at both the extremities of the device requires an iterative procedure of integration: (1) integrate from 𝑧=0 to 𝑧=𝐿 the equations for 𝑃+𝑘𝑚(𝑧) and assuming 𝑃−𝑘𝑚(𝑧)=0 in the first iteration; (2) integrate from 𝑧=𝐿 to 𝑧=0 the equations for 𝑃−𝑘𝑚(𝑧), using the values of 𝑃+𝑘𝑚(𝑧) found previously; (3) restart from point (1), using for 𝑃−𝑘𝑚(𝑧) the value found at point (2). The iterations are stopped when the change in 𝑃±𝑘𝑚(𝑧) is smaller than a prescribed value.

A particular care should be taken for the evaluation of Γ2,𝑚(𝑧) because this quantity determines the incremental gain in (13) and depends on the actual value of the approximation for 𝑃±𝑘𝑚(𝑧); this implies that, at each step of the forward and backward integration along 𝑧, we have to evaluate the integral Γ2,𝑚(𝑧). Because of our choice of a Runge-Kutta formula (4th order, 4 stages) as the basic integration method, at each step we should evaluate four times (one for each stage) the overlap integral, within a large expense of computer time. Therefore, we decided to perform the transversal integration only at the first stage at each step, and to use the computed value also for the other three stages; this approximation appears not to have a great influence on the solution and represents a good compromise between accuracy and computation time. In our simulations, the spontaneous emission spectrum is divided into 100 slots which corresponds to a wavelength resolution Δ𝜆 = 2 nm in the region 1450–1650 nm.

The intensity profiles 𝑖𝑚(𝑥,𝑦) used in the evaluation of Γ2,𝑚(𝑧) and Γ𝑇,𝑚 are introduced in the model as a set of measured values, while the dopant distribution 𝑑(𝑥,𝑦) is approximated by Gaussian, erfc, or constant functions in depth (𝑦) and in width (𝑥). The simulation of the optical amplification in Er3+-doped LiNbO3 waveguide has been performed using parameters obtained from the literature [5–10]. Using (1), (2), (10), (11) and (13)–(17) we have calculated numerically the spectral dependence of the quality factor, the power spectral density, the Fano factor, and the statistical fluctuation for a single-pass configuration of the optical amplifier. We assumed a 1484 nm pump and a signal at 𝜆 = 1531 nm, having 1 μW input power.

We used the following values for the absorption (ğ‘Ž) and emission (𝑒) cross-sections of the pump (𝑝) and signal (𝑠) for TE and TM polarizations: ğœŽğ‘ŽTE(1484 nm) = 5.61 × 10−25 m2,ğœŽğ‘ŽTM(1484 nm) = 3.46 × 10−25 m2, ğœŽğ‘’TE(1484 nm) = 1.92 × 10−25 m2,ğœŽğ‘’TM(1484 nm) = 1.105 × 10−25 m2,ğœŽğ‘ŽTE(1532 nm) = 17.24 × 10−25 m2,ğœŽğ‘ŽTM(1532 nm)= 12.15 × 10−25 m2, ğœŽğ‘’TE(1532 nm) = 16.36 × 10−25 m2, ğœŽğ‘’TM(1532 nm) = 11.53 × 10−25 m2. The Er-profile has been considered Gaussian in depth and constant in width, with a surface concentration of about 7 × 1025 m−3 and a diffusion depth of 20 μm and of 5.12 μm in width (defined at 1/𝑒). We assumed the following values for the scattering loss and spontaneous emission lifetime: 𝛼 = 3.7 dB·m−1) for TE, 𝛼 = 4.8  dB·m−1 for TM in the case of straight waveguides, and 𝛼 = 0.4 dB·cm−1 for TE and TM in the case of the curved ones and 𝜏 = 2.6 ms, respectively. The length of the waveguide in our simulations is 𝐿 = 5.4 cm and the pump or signal assumed to be TE polarized if not explicitly stated.

Figure 1 presents the spectral evolution of the quality factor of the signal (13), the signal gain being defined as 𝐺(𝑧)=ln⌊𝑃signal(𝑧)/𝑃signal(0)⌋) in the case of single-pass (signal output at 𝑧=𝐿 with 𝑅(𝐿)=0) pass configuration, high pump regime (100 mW incident pump power) for single and for the ℳ-mode operation.

As can be seen from Figure 1 the peak values of the quality factor are greater in the single mode operation because the noise figures in the case of ℳ-mode operation are greater than those corresponding to the single one. The ℳ-mode operation in comparison with the single one determines the diminution of the gain and the enhancement of the noise figure because the overlap integral between the population of the excited level and the normalized intensity field profile is smaller in the case of ℳ-mode operation than in the case of the single one (i.e., 1.21 times in the case of Gaussian profile of the dopant for a waveguide having 5.4 cm length and for an input pump power of 100 mW).

Considering a radius of curvature of the waveguide is about 5 cm and its length is 2.5 cm (values commonly used in the laboratories when manufacturing integrated optics Mach-Zehnder interferometers) we obtained a quality factor of about 0.13 for 𝜆 = 1.531 μm.

The power spectral density (Equation (15) is presented in Figure 2). For a high input pump power in the ℳ-mode operation the spontaneous emission factor (which determines the power spectral density) is less than 2 over most of the 60 nm spectral range considered in the figure (i.e., for Gaussian profile of the dopant in the spectral ranges 1.45 μm ÷ 1.48 μm and 1.62 μm ÷ 1.65 μm, which correspond to the amplification of an equivalent input noise of one photon per unit frequency).

In Figures 3 and 4 the spectral dependences of the Fano factor (3) and the statistical fluctuation (4) in single pass configuration are presented.

The dependence of the Fano factor in the case of TE polarization for a signal at 𝜆 = 1531 nm versus the input pump power for a straight waveguide length 𝐿 = 5.4 cm in the single and ℳ-mode operation is presented in Figure 5.

As can be seen form Figure 5 the Fano factor increases rather low with the increasing of the pump power. In the same conditions the statistical fluctuation decreases when the pump power is augmented. This fact confirms that the output statistics is approximately Poissonian for pump powers about 100 mW.

The Fano factor and the statistical fluctuation increase several orders of magnitude with the waveguide length over a range ~10 cm in both single and ℳ-mode operation. In the case of straight waveguides, the output statistics is approximately Poissonian for waveguide lengths up to 6 cm.

Computing the same parameters in the case of a curved waveguide having a length of about 2.5 cm for the above mentioned conditions we concluded that the photon statistics can be assumed to be Poissonian.

Therefore, it seems that the amplifier length is responsible for the reduction of the Poissonian aspect of the amplified light. This is not the case for shorter waveguide lengths but higher pumping levels, for which the statistical properties of the light are roughly maintained. A consequence concerns the design of complex structures, when the coherence of the amplified light must be conserved. If the miniaturization of the integrated devices has to be considered, it is preferable to choose short lengths and high-level pumpings, rather than long arms and low pumping regimes.

The above obtained simulated results show that the photon statistics which characterize the Er3+-doped Ti:LiNbO3 waveguide amplifiers (the Fano factor, the statistical fluctuation and the spontaneous emission factor) are correlated with the amplifying phenomena of light (gain, noise figure, signal-to-noise ratio, quality factor).

4. Conclusions

In the small gain approximation an original analysis of the output noise statistical properties of a single and ℳ-mode Er3+-doped LiNbO3 straight and curved waveguide amplifiers has been presented. The simulations concern the quality factor, power spectral factor, Fano factor, and the statistical fluctuation. This analysis demonstrates that the Poissonian photon statistics are maintained for pump powers lower than 100 mW and waveguide lengths smaller than 5 cm for straight waveguides and 2.5 cm in the case of curved ones.

Our simulation results concerning the optical gain, noise figure, and quality factor are in agreement with other experimental and theoretical results [15–17].

The above-mentioned parameters are important when the coherence of the amplified light is an issue for the doped waveguides under investigation. In this paper we have shown that the coherence of the output amplified signal is more sensitive to the waveguide length than to the pumping power. The theoretical results of this simulation characterize the Er3+-doped LiNbO3 waveguide amplifiers from the point of view of noise statistical properties and can be used in the better understanding of the amplification process. Also, they can be used for the design of directional couplers, symmetrical and asymmetrical Mach-Zehnder interferometers, and other complex rare earth-doped integrated circuits.


  1. E. Desurvire, Erbium-Doped Fiber Amplifiers, Wiley-Interscience, New York, NY, USA, 1994.
  2. P. Diament and M. C. Teich, “Evolution of the statistical properties of photons passed through a traveling-wave laser amplifier,” IEEE Journal of Quantum Electronics, vol. 28, no. 5, pp. 1325–1334, 1992. View at: Publisher Site | Google Scholar
  3. B. E. A. Saleh and M. C. Teich, Fundamentals of Photonics, John Wiley & Sons, New York, NY, USA, 2nd edition, 2007.
  4. M. F. Saleh, G. Giuseppe, B. E. A. Saleh, and M. C. Teich, “Modal and polarization qubits in Ti:LiNbO3 photonic circuits for a universal quantum logic gate,” Optics Express, vol. 18, no. 19, pp. 20475–20490, 2010. View at: Publisher Site | Google Scholar
  5. S. Reza, H. Herrmann, R. Ricken, V. Quiring, and W. Sohler, “Spectral characteristics of an integrated tunable frequency shifted feedback laser in erbium doped lithium niobate,” in Proceedings of the 8th International Conference on Optoelectronics, Fiber Optics and Photonics (Photonics '06), p. 250, Hyderabad, India, 2006, paper FrA 5. View at: Google Scholar
  6. J. Perina and J. Perina Jr., “Photon statistics of a contradirectional nonlinear coupler,” Quantum and Semiclassical Optics, vol. 7, no. 5, pp. 849–862, 1995. View at: Publisher Site | Google Scholar
  7. G. Lifante, Integrated Photonics: Fundamentals, John Wiley & Sons, West Sussex, England, 2003.
  8. N. N. Puscas, D. Scarano, R. Girardi, and I. Montrosset, “Analysis of output statistics of single and double pass Er-doped LiNbO3 waveguide amplifiers,” Optical and Quantum Electronics, vol. 29, no. 8, pp. 799–809, 1997. View at: Google Scholar
  9. N. N. Puscas, B. Wacogne, A. Ducariu, and B. Grappe, “Modelling the output statistics of Er-doped LiNbO3 curved waveguide amplifiers,” Journal of Modern Optics, vol. 46, no. 6, pp. 1017–1030, 1999. View at: Google Scholar
  10. N. N. Puscas, “Modelling the spectral noise of single and double pass Er3+-doped Ti:LiNbO3M-mode straight waveguide amplifiers,” Journal of Optoelectronics and Advanced Materials, vol. 4, no. 4, pp. 911–921, 2002. View at: Google Scholar
  11. N. N. Puşcaş, “Analysis of the gain and photon ststistics in 3+Ti:LiNbO3M-mode straight waveguide amplifiers,” University Politehnica of Bucharest, Scientific Bulletin A, vol. 65, no. 1, pp. 73–82, 2003. View at: Google Scholar
  12. Y. N. Singh, H. M. Gupta, and V. K. Jain, “Optical amplifiers in broadcast optical networks: a survey,” IETE Technical Review, vol. 16, no. 5, 6, pp. 449–459, 1999. View at: Google Scholar
  13. I. Turek, I. Martinèek, and R. Stránsky, “Interference of modes in optical fibers,” Optical Engineering, vol. 39, no. 5, pp. 1304–1309, 2000. View at: Publisher Site | Google Scholar
  14. D. Marcuse, Light Transmission Optics, Van Nostrand Reinhold, New York, NY, USA, 1972.
  15. R. Brinkmann, I. Baumann, M. Dinand, W. Sohler, and H. Suche, “Er-doped single and double pass Ti:LiNbO3 waveguide amplifiers,” IEEE Journal of Quantum Electronics, vol. 30, no. 10, pp. 2356–2360, 1994. View at: Publisher Site | Google Scholar
  16. J. Guo-Liang, S. Gong-Wang, M. Huan, H. Li-Li, and L. Qu, “Gain and noise figure of a double-pass waveguide amplifier based on Er/Yb-doped phosphate glass,” Chinese Physics Letters, vol. 22, no. 11, pp. 2862–2869, 2005. View at: Publisher Site | Google Scholar
  17. G. Jain, A. Kapoor, and E. K. Sharma, “Er-LiNbO3 waveguide: field approximation for simplified gain calculations in DWDM application,” Journal of the Optical Society of America B, vol. 26, no. 4, pp. 633–639, 2009. View at: Publisher Site | Google Scholar

Copyright © 2011 Niculae N. Puscas. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.