Advances in OptoElectronicsVolume 2011 (2011), Article ID 368327, 6 pagesdoi:10.1155/2011/368327

Research Article

Modelling of Output Statistics of Single and -Mode Straight and Curved :Ti:LiNb Waveguide Amplifiers

Niculae N. Puscas

Physics Department I, Faculty of Applied Science, University “Politehnica” of Bucharest, Splaiul Independentei 313, 060042 Bucharest, Romania

Received 29 April 2011; Accepted 20 July 2011

Academic Editor: Xian Cao

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.

Abstract.

A theoretical analysis of some statistical parameters which characterize the -doped Ti:LiNbO₃ 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 -doped fibres and waveguide amplifiers [1–12].

A novel technique for evaluating the output statistics of single and -mode straight and curved -doped Ti:LiNbO₃ 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 ions incorporated into the LiNbO₃ 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 is achieved after the optical pumping to a pump band (an excited state) followed by a fast relaxation to the metastable level . 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 -doped waveguide amplifiers in the wavelength range 1440 nm < < 1640 nm.

The interaction of pump/signal photons with 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] 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 is the steady state solution of the rate equations for the upper populations of the two levels which is written in the form

In (2) the ions distribution in the waveguide has been written in the form , where 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 being the fluorescence lifetime of the ions. The absorption and emission rates of the amplified spontaneous emission are and, 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 and , 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]: 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 where are the input and output reflectivities and are the optical powers injected at or at . We assume that the signal is injected only from the section and that no noise power is entering the amplifier (i.e., for the ASE beams). By properly choosing the values of and we can account for the natural reflectivities of the LiNbO₃-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 , respectively,

The effects of gain saturation due to the amplified signal or to the ASE power are accounted for in the expression for the population : as the signal and ASE power grows the emission rate increases and decreases, leading to a reduction of and consequently of the gain factor in (1).

Considering a collection of atoms with population density distribution for the ground state and 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]: where

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 : where represent the spectral gain and the ASE photon number, respectively.

For coherent input signals with Poisson’s statistics being well known that . In the case of uniform inversion, when the population densities of the upper and lower levels (and , resp.) do not depend on the waveguide coordinate , the variance takes the canonical form:

As a measure of the gain and the noise characteristics of the waveguide amplifiers, we adopted a quality factor where 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]: where represents the spontaneous emission factor and the Plank's constant. The amplifier output statistics can also be characterised by the Fano factor and the statistical fluctuation [1]. For Poisson statistics (coherent light), the Fano factor and the statistical fluctuation correspond to and , while for Bose-Einstein statistics (incoherent light), they are given by and .

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 where the spontaneous emission factor is given by

In (16) the factor 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] where are the propagation constants which are different for different modes because of their velocity [11] and are the overlap coefficients between the pump optical field, , and the mode fields, , excited in the waveguide, which satisfy a normalization condition . 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 to the equations for and assuming in the first iteration; (2) integrate from to 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 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 . 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 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 -doped LiNbO₃ 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: (1484 nm) = 5.61 × 10−25 m²,(1484 nm) = 3.46 × 10−25 m², (1484 nm) = 1.92 × 10−25 m²,(1484 nm) = 1.105 × 10−25 m²,(1532 nm) = 17.24 × 10−25 m²,(1532 nm)= 12.15 × 10−25 m², (1532 nm) = 16.36 × 10−25 m², (1532 nm) = 11.53 × 10−25 m². The Er-profile has been considered Gaussian in depth and constant in width, with a surface concentration of about 7 × 10²⁵ m⁻³ 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⁻¹) for TE, = 4.8  dB·m⁻¹ for TM in the case of straight waveguides, and = 0.4 dB·cm⁻¹ 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 ) in the case of single-pass (signal output at with ) pass configuration, high pump regime (100 mW incident pump power) for single and for the -mode operation.

Figure 1: The spectral evolution of the quality factor pass configuration for single (S) and for the -mode (M) 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).

Figure 2: The power spectral density versus wavelength as in Figure 1.

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

Figure 3: The Fano facto versus wavelength as in Figure 1.

Figure 4: The spectral dependence of the statistical fluctuation as in Figure 1.

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.

Figure 5: The Fano factor versus the pump power in the single (S) and -mode (M) operation.

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 -doped Ti:LiNbO₃ 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 -doped LiNbO₃ 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 -doped LiNbO₃ 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.

References

E. Desurvire, Erbium-Doped Fiber Amplifiers, Wiley-Interscience, New York, NY, USA, 1994.

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.

B. E. A. Saleh and M. C. Teich, Fundamentals of Photonics, John Wiley & Sons, New York, NY, USA, 2nd edition, 2007.

M. F. Saleh, G. Giuseppe, B. E. A. Saleh, and M. C. Teich, “Modal and polarization qubits in Ti:LiNbO₃ photonic circuits for a universal quantum logic gate,” Optics Express, vol. 18, no. 19, pp. 20475–20490, 2010.

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.

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.

G. Lifante, Integrated Photonics: Fundamentals, John Wiley & Sons, West Sussex, England, 2003.

N. N. Puscas, D. Scarano, R. Girardi, and I. Montrosset, “Analysis of output statistics of single and double pass Er-doped LiNbO₃ waveguide amplifiers,” Optical and Quantum Electronics, vol. 29, no. 8, pp. 799–809, 1997.

N. N. Puscas, B. Wacogne, A. Ducariu, and B. Grappe, “Modelling the output statistics of Er-doped LiNbO₃ curved waveguide amplifiers,” Journal of Modern Optics, vol. 46, no. 6, pp. 1017–1030, 1999.

N. N. Puscas, “Modelling the spectral noise of single and double pass Er³⁺-doped Ti:LiNbO₃M-mode straight waveguide amplifiers,” Journal of Optoelectronics and Advanced Materials, vol. 4, no. 4, pp. 911–921, 2002.

N. N. Puşcaş, “Analysis of the gain and photon ststistics in ³⁺Ti:LiNbO₃M-mode straight waveguide amplifiers,” University Politehnica of Bucharest, Scientific Bulletin A, vol. 65, no. 1, pp. 73–82, 2003.

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.

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.

D. Marcuse, Light Transmission Optics, Van Nostrand Reinhold, New York, NY, USA, 1972.

R. Brinkmann, I. Baumann, M. Dinand, W. Sohler, and H. Suche, “Er-doped single and double pass Ti:LiNbO₃ waveguide amplifiers,” IEEE Journal of Quantum Electronics, vol. 30, no. 10, pp. 2356–2360, 1994.

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.

G. Jain, A. Kapoor, and E. K. Sharma, “Er-LiNbO₃ 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.