# Modeling of Mid-IR Amplifier Based on an Erbium-Doped Chalcogenide Microsphere

Authors: P. Bia; A. Di Tommaso; M. De Sario
Journal: International Journal of Optics
 (2012)

---

## Abstract

An optical amplifier based on a tapered fiber and an Er3+-doped chalcogenide microsphere is designed and optimized. A dedicated 3D numerical model, which exploits the coupled mode theory and the rate equations, is used. The main transitions among the erbium energy levels, the amplified spontaneous emission, and the most important secondary transitions pertaining to the ion-ion interactions have been considered. Both the pump and signal beams are efficiently injected and obtained by a suitable design of the taper angle and the fiber-microsphere gap. Moreover, a good overlapping between the optical signals and the rare-earth-doped region is also obtained. In order to evaluate the amplifier performance in reduced computational time, the doped area is partitioned in sectors. The obtained simulation results highlight that a high-efficiency midinfrared amplification can be obtained by using a quite small microsphere.


---

## Body

## 1. Introduction

In recent years, the micro- and nanospherical resonators have attracted great interest for their high quality Q-factor, material versatility, manufacturing easiness, and dopant hosting flexibility for active devices. In rare-earth-doped microspheres, Whispering Gallery Modes (WGMs) can strongly enhance the light-matter interaction, because of their very high Q-factor and small mode volume. A number of application fields could exploit the intriguing WGM properties, such as those involving polarization transmission, coupled-resonator-induced transparency, biosensor analysis [1, 2], nonlinear optics, cavity quantum electrodynamics (QED), and quantum information processing [3].Rare-earth-doped microspheres based on silica, phosphate, tellurite, and ZBLAN glass host materials [3–7] show ultralow lasing thresholds and very narrow emission linewidths. Chalcogenide glass has recently attracted significant interest as a material for the manufacturing of active microsphere resonators. In fact, they exhibit lower modal volumes, higher refractive indices, and high absorption and emission cross-sections [8]. Moreover, these glasses allow realizing efficient mid-IR amplifiers, thanks to their low phonon energy and to their high capability to host rare earth ions.In this work, the mathematical model of a microsphere amplifier is described. The signal wavelength isλs=2.76μm and the pump wavelength is λp=980nm. The model appears more complete than those reported in the literature for active [9–11] and passive [12, 13] devices and it improves the previous work [14]. It involves the coupled mode theory and the rate equations and it allows the simulation of a tapered fiber coupled to the rare-earth-doped chalcogenide glass microsphere. The following features are considered: (i) the radiative and nonradiative rates, at both pump and signal wavelengths, (ii) the stimulated emission at the signal wavelength, (iii) the amplified spontaneous emission noise (ASE), (iv) the lifetime of the considered energy levels, and (v) the ion-ion energy transfers. The model is implemented with a homemade 3D numerical code and the numerical results are shown.The paper is structured as follows: Section2 includes the mathematical model, Section 3 the numerical results, and finally Section 4 the conclusions.


## 2. Mathematical Model

A sketch of the considered system is shown in Figure1. It consists of a microsphere and a tapered fiber both in Ga5Ge20Sb10S65 chalcogenide glass. Ppin and Psin stand for the input pump and signal powers, respectively, Ppout and Psout for the output pump and signal powers, respectively, δ is the angle of the taper, g is the gap between the cavity and the fiber, R0 is the sphere radius, and a is the waist fiber radius and S is the thickness of doped area.Figure 1
Schematic of a microsphere coupled to a tapered fiber.By assuming that the input signals are narrowband, each one gives rise to a singlel, m, n WGM within the microsphere. Moreover, on the assumption that the microsphere is only doped in the outer layer, the most amplified WGM is the fundamental one (n=1).Figure2 shows the transitions among the energy levels considered in the developed numerical code. In particular, the population inversion is possible between 4I11/2 and 4I13/2 energy levels, by pumping at the wavelength λp=980nm. As a consequence, the amplification can occur at the signal wavelength λs=2760nm. Moreover, the energy transitions due to cross-relaxation and cooperative upconversion effects have also to be considered in the simulations, because of the high dopant concentration and the comparable lifetime of 4I13/2, 4I11/2, and 4I9/2 energy levels. The numerical values of absorption, emission cross-sections, upconversion, and cross-relaxation coefficients as well as the energy level lifetimes are reported in Table 1.Table 1
Spectroscopic parameters ofEr3+-doped Ga5Ge20Sb10S65 chalcogenide glass.
Parameters
Values
σ
13 at 980 nm
1.32
×
10
-
24
m
2
σ
32 at 2760 nm
1.60
×
10
-
24
m
2
C
33
20
×
10
-
24
m
2
s
-
1
C
22
30
×
10
-
24
m
2
s
-
1
C
41
5
×
10
-
24
m
2
s
-
1
τ
21
1.83
×
10
-
3
s
τ
32
1.37
×
10
-
3
s
τ
43
1.08
×
10
-
3
s
τ
54
0.13
×
10
-
3
sFigure 2
Diagram of energy levels.The active behavior is modeled by calculating the coupling coefficients, the quality factor, the mode volume, the transition rates. These parameters are obtained with an accurate evaluation of electromagnetic field mode properties in both fiber and microsphere [15]. The following conditions are considered:(i)
the energy gaps between the Stark levels of each manifold are small when compared to the energy separation between different manifolds;(ii)
the rates between the Stark levels are much faster than those between two different manifolds;(iii)
the population in any manifold is in local thermal equilibrium with the glass lattice.The following differential equations describe the time-variation of the field amplitudeAs (signal) and Ap (pump) inside the microsphere, considering both the dopant ions and the fiber coupling:
(1)dAs(t)dt=-12(1τexts+1τ0s)As(t)+j1τextsTcsAsin(t)-jΔωsAs(t)+gsAs(t),dAp(t)dt=-12(1τextp+1τ0p)Ap(t)+j1τextpTcpApin(t)-jΔωpAp(t)+gpAp(t),
where Asin and Apin are the field amplitudes at the input end of the optical fiber; being a={s,p}, Δωa=ωWGMa-ωina is the frequency detuning of the fiber input signal from the WGM frequency ωWGM, Tca=2πR0neffa/c is the circulating time inside the microsphere (round trip time), c is the speed of light in vacuum, and neffa is the WGM effective refractive index. More precisely, the effective index is almost constant in the resonance bandwidth, because of its narrowness. As a consequence, the modal dispersion can be neglected, and the group index results are almost equal to the effective index. The intrinsic lifetime is τ0=1/(κ02)=Q0/ωWGM, where Q0 is the intrinsic quality factor. κ0 is the intrinsic cavity decay rate, which depends on the total losses due to material absorption, surface scattering losses, radiative losses, and whispering gallery losses. Both the radiative and surface scattering losses can be neglected in the microspheres having the diameter around a few tens of microns and a perfect surface. The coupling lifetime τext=1/κext2=mπ/(ωK2), where κext is the cavity decay rate or coupling coefficient, denotes the coupling phenomenon between microsphere and optical fiber. The field overlapping, K, is calculated according to the overlap integral [15]:
(2)K=∫V‍k02(ns2-n02)2βfEf⋅Es*dV,
where k0=2π/λ; ns and n0 are, respectively, the wave vector in the vacuum, the microsphere, and background refractive indices; βf is the fiber propagation constant; Ef is the electric field of the fiber fundamental mode; Es is the WGM field both normalized on the r^·θ^ plane. Moreover, in (1) the gain/attenuation due to the dopant ions, evaluated over a round-trip, is defined as follows:
(3)ga=c2neffaΓa∑i=25‍∑j=1i-1‍(Niσij-Njσji)a={s,p},
where Ni is ion population of the i-th energy level, σji is the absorption cross-section, and σij is the emission cross-section, Γa is the overlap factor of each WGM with the rare-earth-doped region. The Ni values are obtained by using the rate equation model, arising from the energy level diagrams in Figure 2 and reported in [14].The evolution of both the pump and signal field amplitudes at the fiber output end is given by:(4)Aaout(t)=1-TcaτextaAain(t)+jTcaτextaAa(t)a={s,p},
where Aa is evaluated at the nearest point between the microsphere and the optical fiber (r=R0, θ=0, ϕ=π/2). Finally, the overall gain (transmittance) of the amplifier, is calculated by using [16]:
(5)G=|AsoutAsin|2=|1-Tcsτexts+jTcsτextsAsAsin|2.


## 3. Numerical Results

The developed numerical code has a low computational cost (reduced computational time and memory) compared with FDTD and FEM-based algorithms. Moreover, it is flexible and it can be easily used to evaluate the amplifier performance in several configurations. In fact, the geometrical parameters (e.g.,R0, δ, a, g), the operational ones (e.g., modulation, frequency, power of the input signals) as well as the physical ones (the rare earth concentration, the thickness of the doped region, and the refractive indices) can be varied. Moreover, the numerical code can be easily extended to the analysis of more complex rare earth and lasing systems.A number of simulations have been performed to demonstrate the feasibility of2760nm signal amplification. In particular, a parametric investigation is carried out to evaluate the amplifier performances. In the simulations, the following parameters are used: input signal power Psin=10nW, input pump power Ppin=100mW, thickness of doped region S=3μm, for which the corresponding overlap factors are Γs=0.89 for the signal and Γp=0.99 for the pump.Figure3 depicts the signal gain as a function of fiber-microsphere gap g, for three different microsphere radii R0=20 μm, 30 μm and 40 μm, with the taper angle δ=0.03rad, and dopant concentration equal to 0.5% in weight percentage (Ntot=5.77×1025ions/m3). It can be observed that the gain increases with the gap until a value beyond which the transmittance drastically decreases. This is caused by the undercoupling condition for the power signal, being low the power coupled to the microsphere. Undercoupling condition occurs for smaller gap by increasing the microsphere radius, since the structure better confines the evanescent wave.Figure 3
Signal amplification as a function of fiber-microsphere gapg for three different microsphere radii.Figure4 illustrates the signal gain as a function of the taper angle δ, for three different fiber radii a=500nm, 650nm, and 800nm. The choice of the range of the taper angle takes into account the structure durability. The geometrical and physical parameters are R0=30μm, g=540nm, and Ntot=5.77×1025ions/m3. If the fiber radius a=500nm, the gain increases by decreasing the taper angle and also the decreasing of the coupling factor occurs. Moreover, the amplification occurs for high coupling factor values, corresponding to high taper angle values. In particular, the minimum δ value required for the signal amplification depends on the fiber radius. In fact, by using a small fiber radius, the increasing of the coupling factor occurs due to the enhancement of the output evanescent wave. In Figure 4, the undercoupling condition is not occurring for the whole δ range, if a=500nm; whereas the undercoupling condition occurs for δ<0.035rad, if a=800nm.Figure 4
Signal amplification as a function of the taper angle for three different values of the waist fiber radius.Figure5 depicts the signal gain as a function of a fiber radius, for three microsphere radii R0=20μm,30μm,40μm, with taper angle δ=0.03rad, and dopant concentration Ntot=5.77×1025ions/m3. These fiber radii have been considered because lower values determine a signal leakage due to the low confinement strength at the signal wavelength. Whereas, the greater radii induce a multimodal propagation at the pump wavelength. The maximum gain value occurs for the microsphere radius R0=40μm and fiber radius a=500nm. Moreover, by increasing the fiber radius, the coupling decreases leading to the undercoupling condition and reducing the amplifier performances.Figure 5
Signal amplification as a function of the taper angle for three different values of the waist fiber radius.Figure6 shows the signal gain as a function of the erbium concentration NEr, for three different microsphere-fiber gaps g=500mm, 550mm, 600mm, with δ=0.03rad, R0=30μm, a=700nm. By increasing the ion concentration, the enhancement of the signal amplification occurs until a threshold value. In fact, if high erbium concentrations are used, the population inversion occurs as long as the available power is sufficient. The threshold value can be suitably increased by enhancing the coupling (e.g., by reducing the gap) or decreased by reducing the coupling (e.g., by increasing the gap). Low gap values induce the cavity power reduction, reaching the value not allowing the doped ion population inversion, thus limiting the amplifier performances.Figure 6
Signal amplification as a function of the erbium concentration for three different values of fiber gap.


## 4. Conclusions

In this paper, anEr3+-doped chalcogenide microsphere amplifier evanescently coupled with a tapered optical fiber has been designed by means of a homemade 3D numerical code. It includes both the rate equation and coupled mode theory models. Spherical coordinates are used to find the solution of the scalar Helmholtz equation needed for the electromagnetic analysis of the microsphere. The homemade numerical code allows to evaluate the amplifier performance by changing several parameters: fiber-microsphere gap, thickness of erbium doped region, fiber taper angle, erbium concentration, and operative parameters such as pump and signal power.The microsphere resonator shows low threshold power (≈80mW) and high gain (≈8dB) at the signal wavelength λs=2760nm. Moreover the total size of the device is few tens of microns. As future developments, slow light phenomenon could be obtained, exploiting the coherent population oscillations in order to realize buffers for optical telecommunications and optical logic gates. In addition, this device could be used to generate low threshold pump powers. Finally, the model can be easily extended to take into account both the mode degeneracy and the multimodal pump and signal, in order to match the simulation dynamics with the realistic device one. As a drawback, a more complete calculation will require more computational time.



---

*Source: 808679-2012-07-15.xml*