Fiber-Optic Ring Resonator Interferometer

Ramón José Pérez Menéndez

## Abstract

In this chapter, theoretical analysis and electro-optical characterization of a fiber optic passive ring resonator interferometer (FOPRRI) were realized. First, a theoretical description and analysis of all-pass and add-drop filter configurations were performed, hence obtaining their respective transfer function, taking into account, the physical properties of a standard single-mode optical fiber as resonator constitutive material. Second, computational numerical simulations of a typical FOPRRI were performed with different analysis methods: effective index method (EIM), coupled mode theory (CMT), finite element method (FEM) and finite difference time domain (FDTD). Third, a comparative analysis of the results obtained with these methods is realized, showing the most accurate and appropriate for filter electro-optical characterization as a consequence. Last, and as an ultimate step, taking as reference the previously obtained electro-optical characterization for both configurations, the main applications are derived.

Keywords: passive optical ring resonator interferometer, single-mode optical fiber, all-pass filter, add-drop filter, EIM (effective index method), CMT (coupled mode theory), FEM (finite element method), FDTD (finite difference time domain)

## 1. Introduction

Optical interferometry constitutes an important technique used in high number of measurement processes for multiple physical magnitudes and quantitative phenomenon [1]. Particularly, fiber-optic waveguides can act as very useful and efficient transmission medium for light guidance in many interferometric sensor devices. The issue of fiber-optic interferometric devices has received an extensive treatment in the literature [2, 3]. Also, the general theory of optical resonators was extensively treated by authors, for example, Yariv and Yeh [4]. However, quantitative analytic treatment of a fiber optic passive ring resonator interferometer (FOPRRI) can be rather found in singular applications, particularly, in gyroscopes. A clear proof of it is that one of the pioneering and singular descriptions of fiberoptic ring interferometer devices was in 1976 when Vali and Shorthill designed the first fiber-optic gyroscope model [5–7]. Later, many other authors studied in detail this first design and improved some of its electro-optical parameters [8, 9]. Considered from a broader point of view, the FOPRRI device has now evolved toward what we could call the integrated optical passive ring resonator interferometer (IOPRRI). In this case, the fiber-optic ring resonator is replaced by a ring-integrated optical waveguide, made of high refractive index material, such as lithium niobate,

indium phosphide, silicon nitride, or silicon on insulator (SOI) platform. This chapter focuses on the analysis and testing of fiber-optic ring resonator. Therefore, the integrated ring resonator study, its investigations and challenges are beyond the scope of it [10].

these two frequencies gives a measure of gyro output signal, which, in turn, is proportional to rotation rate to which the fiber-optic ring resonator system is

2. Analysis of the fiber-optic all-pass filter ring resonator configuration

The all-pass ring resonator configuration is created by taking one output of a generic 2 � 2 directional coupler and feedback into one input, as it can be seen in Figure 2. Such a device exhibits a periodic cavity resonance when light traversing the ring acquires a phase shift corresponding to an integer multiple of 2π radians. The resonator can mathematically be formulated from two components: a coupling strength (from coupling coefficient κ) and a feedback path (via transmission coefficient t), in contrast to the infinite sum derivations made for Fabry-Perot and Gires-Tournois interferometers. Therefore, following the text of Rabus [20], we will derive the basic spectral properties by assuming the steady-state operation and the matching fields. Although both methods are equally valid, the field-matching method has the advantage of simplicity. The basic relations among the incident Ei1, transmitted Et1, and circulating Ei2, Eit<sup>2</sup> fields of a single resonator are derived by combining the relations for the coupler with that of the feedback path. In the frequency domain, the fields exiting the coupling region are related to the input

> Et<sup>1</sup> ¼ tEi<sup>1</sup> þ i κEi<sup>2</sup> (2) Et<sup>2</sup> ¼ tEi<sup>2</sup> � i κ Ei<sup>1</sup> (3)

> > �<sup>1</sup> <sup>p</sup> ) that takes into account, the

� � (4)

subjected.

fields via the following equations:

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

the following unitary matrix:

Figure 2.

51

where factor <sup>i</sup> is the imaginary unit (<sup>i</sup> <sup>¼</sup> ffiffiffiffiffiffi

Et<sup>1</sup> Et<sup>2</sup> � �

All-pass fiber-optic ring resonator configuration ð Þ κ<sup>1</sup> ≈ κ<sup>2</sup> ¼ κ and t<sup>1</sup> ≈ t<sup>2</sup> ¼ t .

phase shift caused by the cross-coupling of the waveguides. For simplicity, we assume the total symmetry of cross-coupling and bar-coupling coefficients for both

waveguides, and therefore, the two approximate equations are ðκ<sup>1</sup> ≈ κ<sup>2</sup> ¼ κ and t<sup>1</sup> ≈ t<sup>2</sup> ¼ tÞ. The two Eqs. (2) and (3) can be written in a compact way via

> <sup>¼</sup> t i<sup>κ</sup> �iκ t

� � Ei<sup>1</sup>

Ei<sup>2</sup>

An all-fiber passive ring resonator interferometer is formed by a directional coupler and a fiber-optic loop, as shown in Figure 1. Two main geometrical configurations are possible to reach such an optical coupling: (a) direct-coupled fiber ring resonator or (b) cross-coupled fiber ring resonator, as it can be seen in figure. Superimposition of successive beams generated by re-coupling the first coupled beam into the fiber ring results in multi-beam optical interference in a similar way to a Fabry-Pérot interferometer. Resonance will occur when the increase of the phase of the beam around the ring equals 2π n � <sup>π</sup>=<sup>2</sup> � �, <sup>π</sup>=<sup>2</sup> being the phase shift due to cross-coupling and n a positive integer (n = 1, 2, 3, …), so that the light reentering the ring after a single circulation around it is in phase with the light entering the ring for the first time. Following L. F. Stokes et al. [12, 13], the transfer function of the system for the optimum amplitude resonance condition is given by:

$$\frac{I\_T}{I\_0} = (\mathbf{1} - \boldsymbol{\chi}) \left\{ \mathbf{1} - \left[ \frac{\left(\mathbf{1} - \boldsymbol{\kappa}\right)^2}{\mathbf{1} + \boldsymbol{\kappa}^2 + 2\boldsymbol{\kappa}\sin\left(\beta L\right)} \right] \right\} \tag{1}$$

where γ is the insertion loss of the coupler, κ ¼ ð Þ 1 � γ exp ð Þ �2αL is the optical resonance coupling coefficient, α is the fiber attenuation factor, L is the total length of the fiber ring and β is the propagation constant of the re-circulating mode. The transfer function is very similar to a Fabry-Pérot interferometer operating in back-reflection mode. A finesse F of more of 1000 has been demonstrated for very low excess loss coupler and low-loss fibers [11, 12]. At resonance, significant optical power is recirculating in the ring and nonlinear effects as stimulated Brillouin scattering (SBS) can be observed for input powers of less than 100 μW [13]. As it will be seen later, fiber-optic ring resonators have been used in a number of applications, including fiber ring lasers [8], sensors [14], fiber laser gyroscopes [15, 16], optical spectrum analyzers [17, 18] and optical delay lines [19], among others. One of the first and main sensing applications of fiber-optic ring resonator is as a gyroscope and the major attraction of this kind of resonator gyro is that the cavity length can be taken much shorter than the convectional gyro in the ratio of 2/F. In order to detect the phase shift induced by rotation, a highly coherent laser light source is required. A dual frequency modulation scheme is used to control the laser frequency and tracking the resonance condition [8, 9]. The difference between

#### Figure 1.

Fiber-optic-passive-ring-resonator-interferometer (FOPRRI) all-pass configuration: (a) direct-coupled fiber ring-resonator and (b) cross-coupled fiber ring-resonator.

indium phosphide, silicon nitride, or silicon on insulator (SOI) platform. This chapter focuses on the analysis and testing of fiber-optic ring resonator. Therefore, the integrated ring resonator study, its investigations and challenges are beyond the

Interferometry - Recent Developments and Contemporary Applications

phase of the beam around the ring equals 2π n � <sup>π</sup>=<sup>2</sup>

IT I0

An all-fiber passive ring resonator interferometer is formed by a directional coupler and a fiber-optic loop, as shown in Figure 1. Two main geometrical configurations are possible to reach such an optical coupling: (a) direct-coupled fiber ring resonator or (b) cross-coupled fiber ring resonator, as it can be seen in figure. Superimposition of successive beams generated by re-coupling the first coupled beam into the fiber ring results in multi-beam optical interference in a similar way to a Fabry-Pérot interferometer. Resonance will occur when the increase of the

to cross-coupling and n a positive integer (n = 1, 2, 3, …), so that the light reentering the ring after a single circulation around it is in phase with the light entering the ring for the first time. Following L. F. Stokes et al. [12, 13], the transfer function of

where γ is the insertion loss of the coupler, κ ¼ ð Þ 1 � γ exp ð Þ �2αL is the optical resonance coupling coefficient, α is the fiber attenuation factor, L is the total length of the fiber ring and β is the propagation constant of the re-circulating mode. The transfer function is very similar to a Fabry-Pérot interferometer operating in back-reflection mode. A finesse F of more of 1000 has been demonstrated for very low excess loss coupler and low-loss fibers [11, 12]. At resonance, significant optical power is recirculating in the ring and nonlinear effects as stimulated Brillouin scattering (SBS) can be observed for input powers of less than 100 μW [13]. As it will be seen later, fiber-optic ring resonators have been used in a number of applications, including fiber ring lasers [8], sensors [14], fiber laser gyroscopes [15, 16], optical spectrum analyzers [17, 18] and optical delay lines [19], among others. One of the first and main sensing applications of fiber-optic ring resonator is as a gyroscope and the major attraction of this kind of resonator gyro is that the cavity length can be taken much shorter than the convectional gyro in the ratio of 2/F. In order to detect the phase shift induced by rotation, a highly coherent laser light source is required. A dual frequency modulation scheme is used to control the laser frequency and tracking the resonance condition [8, 9]. The difference between

Fiber-optic-passive-ring-resonator-interferometer (FOPRRI) all-pass configuration: (a) direct-coupled fiber

2

1 þ κ<sup>2</sup> þ 2 κ sin ð Þ βL

( ) " #

the system for the optimum amplitude resonance condition is given by:

<sup>¼</sup> ð Þ <sup>1</sup> � <sup>γ</sup> <sup>1</sup> � ð Þ <sup>1</sup> � <sup>κ</sup>

� �, <sup>π</sup>=<sup>2</sup> being the phase shift due

(1)

scope of it [10].

Figure 1.

50

ring-resonator and (b) cross-coupled fiber ring-resonator.

these two frequencies gives a measure of gyro output signal, which, in turn, is proportional to rotation rate to which the fiber-optic ring resonator system is subjected.

## 2. Analysis of the fiber-optic all-pass filter ring resonator configuration

The all-pass ring resonator configuration is created by taking one output of a generic 2 � 2 directional coupler and feedback into one input, as it can be seen in Figure 2. Such a device exhibits a periodic cavity resonance when light traversing the ring acquires a phase shift corresponding to an integer multiple of 2π radians. The resonator can mathematically be formulated from two components: a coupling strength (from coupling coefficient κ) and a feedback path (via transmission coefficient t), in contrast to the infinite sum derivations made for Fabry-Perot and Gires-Tournois interferometers. Therefore, following the text of Rabus [20], we will derive the basic spectral properties by assuming the steady-state operation and the matching fields. Although both methods are equally valid, the field-matching method has the advantage of simplicity. The basic relations among the incident Ei1, transmitted Et1, and circulating Ei2, Eit<sup>2</sup> fields of a single resonator are derived by combining the relations for the coupler with that of the feedback path. In the frequency domain, the fields exiting the coupling region are related to the input fields via the following equations:

$$E\_{i1} = tE\_{i1} + i\,\kappa E\_{i2} \tag{2}$$

$$E\_{t2} = tE\_{i2} - i\,\kappa E\_{i1} \tag{3}$$

where factor <sup>i</sup> is the imaginary unit (<sup>i</sup> <sup>¼</sup> ffiffiffiffiffiffi �<sup>1</sup> <sup>p</sup> ) that takes into account, the phase shift caused by the cross-coupling of the waveguides. For simplicity, we assume the total symmetry of cross-coupling and bar-coupling coefficients for both waveguides, and therefore, the two approximate equations are ðκ<sup>1</sup> ≈ κ<sup>2</sup> ¼ κ and t<sup>1</sup> ≈ t<sup>2</sup> ¼ tÞ. The two Eqs. (2) and (3) can be written in a compact way via the following unitary matrix:

$$
\begin{pmatrix} E\_{t1} \\ E\_{t2} \end{pmatrix} = \begin{pmatrix} t & i\kappa \\ -i\kappa & t \end{pmatrix} \begin{pmatrix} E\_{i1} \\ E\_{i2} \end{pmatrix} \tag{4}
$$

Figure 2. All-pass fiber-optic ring resonator configuration ð Þ κ<sup>1</sup> ≈ κ<sup>2</sup> ¼ κ and t<sup>1</sup> ≈ t<sup>2</sup> ¼ t .

where sign and i-factor in iκ elements of first matrix on the right side of equation are due to cross-coupling phase shift caused by coupler. Also, since this matrix is an unitary matrix, we find that accomplishes the following relationship between bar- (t) and cross- (κ) coupling coefficients (assuming both are independent of frequency):

$$t^2 + \kappa^2 = 1\tag{5}$$

Then, the transmitted power to through port (3) becomes 0. The critical coupling is due to destructive interference between the direct beam and the crosscoupled beam of light, both after passing their respective coupling zone. A condition termed undercoupling does occur when the ring loss exceeds the coupling

<sup>2</sup> < αL, also κ<sup>2</sup> > αL) and is the conventional mode of operation for an

place, when the round-trip loss does not exceed the coupling strength (t < a or

all-pass resonator. Figure 3(a) and (b) displays the effect of attenuation on the transmission and build-up factor for the overcoupled, critically coupled and undercoupled regimes. Experimental demonstration of the transmission

characteristics in each of these regimes for a fiber all-pass filter ring resonator can

A good idea of the behavior of the system, when resonance condition accomplishes, can be obtained when one represents the T transmission function defined as the ratio Pt1/Pi<sup>1</sup> (i.e., the relationship between the through-port power and inputport power) versus normalized detuning ϕ, Figure 3(a). Starting from input field expressions and after some derivations, the following equation can be found:

Figure 3(a) represents the T transmission function for an all-pass ring resonator with t = 0.75 and varying loss a as a parameter (from a = 0.00 up to a = 0.99). It can be seen that, in this case, critical coupling occurs when a = t = 0.75; therefore, the T

<sup>2</sup> > αL, also κ<sup>2</sup> < αL). An overcoupling condition takes

<sup>¼</sup> <sup>a</sup><sup>2</sup> � <sup>2</sup>j j<sup>t</sup> <sup>a</sup> cos <sup>ϕ</sup> <sup>þ</sup> j j<sup>t</sup> <sup>2</sup>

<sup>1</sup> � <sup>2</sup>j j<sup>t</sup> <sup>a</sup> cos <sup>ϕ</sup> <sup>þ</sup> j j<sup>t</sup> <sup>2</sup> <sup>a</sup><sup>2</sup> (11)

strength (t > a or equally: t

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

> <sup>T</sup> <sup>¼</sup> Pt<sup>1</sup> Pi<sup>1</sup>

<sup>¼</sup> Et<sup>1</sup> Ei<sup>1</sup> 

 

2

transmission function drops to zero when ϕ ¼ 0 for this a-factor value. As a

T transmission function (a) and B build-up factor (b) of ring-resonator, both as a function of normalized

detuning ϕ, (curves extracted from [21], with permission).

equally: t

Figure 3.

53

be found in [22].

The feedback path of length 2 π R (R being the ring radius) connects the output from port 4 back into input port 2, where the field can be expressed as:

$$E\_{i2} = e^{-a2\pi R} e^{i\beta 2\,\pi R} E\_{t2} \equiv a e^{i\phi} E\_{t2} \tag{6}$$

where α is the attenuation factor of the fiber, a-factor represents the ring singlepass field amplitude transmission (taking into account the fiber attenuation), β ¼ ð Þ <sup>2</sup>π=<sup>λ</sup> neff is the propagation constant of propagating mode in the fiber (where in time, neff is the effective index of the fiber) and ϕ represents the single-pass ring round-trip phase shift.

Since that adding or subtracting an integer number m of 2 π radians from singlepass phase shift does not change the value of the function, the single-pass phase shift for all resonances is defined such that its value is zero for a local resonance of interest. Furthermore, because the single-pass phase shift is directly related to the radian frequency as ϕ ¼ ωTR, where TR is the transit time of the round-trip of ring resonator and then ϕ is clearly representative of a normalized frequency detuning. Following the previous derivations, one can demonstrate that the transmission power Pt<sup>1</sup> in the output port (3) of the straight waveguide can be expressed as:

$$P\_{t1} = \left| E\_{t1} \right|^2 = \frac{a^2 + \left| t \right|^2 - 2\left| a \left| t \right| \cos \left( \phi + \phi\_c \right)}{1 + a^2 \left| t \right|^2 - 2\left| a \left| t \right| \cos \left( \phi + \phi\_c \right)} \tag{7}$$

where t ¼ j jt exp jϕ<sup>c</sup> ð Þ, j jt representing the coupling losses and ϕ<sup>c</sup> represents the phase of the coupler. Making a similar deduction, we can reach the circulating power in the ring, which is:

$$P\_{i2} = \left| E\_{i2} \right|^2 = \frac{a^2 \left( 1 - \left| t \right|^2 \right)}{1 + a^2 \left| t \right|^2 - 2\left| a \left| t \right| \cos \left( \phi + \phi\_c \right)} \tag{8}$$

On resonance, condition ϕ þ ϕ<sup>c</sup> ð Þ¼ 2π m is fulfilled, where m is an integer, and then Eqs. (7) and (8) becomes:

$$P\_{t1} = \left| E\_{t1} \right|^2 = \frac{\left( a - \left| t \right| \right)^2}{\left( 1 - a \left| t \right| \right)^2} \tag{9}$$

$$P\_{i2} = \left| E\_{i2} \right|^2 = \frac{a^2 \left( 1 - \left| t \right|^2 \right)}{\left( 1 - a \left| t \right| \right)^2} \tag{10}$$

where, as we have said above, a is the loss coefficient of the ring (for zero loss <sup>α</sup> = 0, then <sup>a</sup> = 1). The special case that occurs when <sup>a</sup> <sup>¼</sup> j j<sup>t</sup> (or <sup>κ</sup><sup>2</sup> <sup>¼</sup> <sup>α</sup>L) is named critical coupling and is of high interest in applications. This happens when internal losses around the ring are equal to the coupling losses.

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

where sign and i-factor in iκ elements of first matrix on the right side of equation are due to cross-coupling phase shift caused by coupler. Also, since this matrix is an unitary matrix, we find that accomplishes the following relationship between bar- (t) and cross- (κ) coupling coefficients (assuming both are indepen-

t

from port 4 back into input port 2, where the field can be expressed as:

�<sup>α</sup> <sup>2</sup> <sup>π</sup> <sup>R</sup> e

pass field amplitude transmission (taking into account the fiber attenuation), β ¼ ð Þ <sup>2</sup>π=<sup>λ</sup> neff is the propagation constant of propagating mode in the fiber (where in time, neff is the effective index of the fiber) and ϕ represents the single-pass ring

Ei<sup>2</sup> ¼ e

Interferometry - Recent Developments and Contemporary Applications

Pt<sup>1</sup> ¼ j j Et<sup>1</sup>

Pi<sup>2</sup> ¼ j j Ei<sup>2</sup>

2 ¼

Pt<sup>1</sup> ¼ j j Et<sup>1</sup>

Pi<sup>2</sup> ¼ j j Ei<sup>2</sup>

losses around the ring are equal to the coupling losses.

2 ¼

where, as we have said above, a is the loss coefficient of the ring (for zero loss <sup>α</sup> = 0, then <sup>a</sup> = 1). The special case that occurs when <sup>a</sup> <sup>¼</sup> j j<sup>t</sup> (or <sup>κ</sup><sup>2</sup> <sup>¼</sup> <sup>α</sup>L) is named critical coupling and is of high interest in applications. This happens when internal

The feedback path of length 2 π R (R being the ring radius) connects the output

where α is the attenuation factor of the fiber, a-factor represents the ring single-

Since that adding or subtracting an integer number m of 2 π radians from singlepass phase shift does not change the value of the function, the single-pass phase shift for all resonances is defined such that its value is zero for a local resonance of interest. Furthermore, because the single-pass phase shift is directly related to the radian frequency as ϕ ¼ ωTR, where TR is the transit time of the round-trip of ring resonator and then ϕ is clearly representative of a normalized frequency detuning. Following the previous derivations, one can demonstrate that the transmission power Pt<sup>1</sup> in the output port (3) of the straight waveguide can be expressed as:

<sup>2</sup> <sup>¼</sup> <sup>a</sup><sup>2</sup> <sup>þ</sup> j j<sup>t</sup> <sup>2</sup> � <sup>2</sup> a tj j cos <sup>ϕ</sup> <sup>þ</sup> <sup>ϕ</sup><sup>c</sup> ð Þ

<sup>a</sup><sup>2</sup> <sup>1</sup> � j j<sup>t</sup> <sup>2</sup>

where t ¼ j jt exp jϕ<sup>c</sup> ð Þ, j jt representing the coupling losses and ϕ<sup>c</sup> represents the phase of the coupler. Making a similar deduction, we can reach the circulating

On resonance, condition ϕ þ ϕ<sup>c</sup> ð Þ¼ 2π m is fulfilled, where m is an integer, and

<sup>2</sup> <sup>¼</sup> ð Þ <sup>a</sup> � j j<sup>t</sup> <sup>2</sup>

<sup>a</sup><sup>2</sup> <sup>1</sup> � j j<sup>t</sup> <sup>2</sup>

<sup>1</sup> <sup>þ</sup> <sup>a</sup><sup>2</sup> j j<sup>t</sup> <sup>2</sup> � <sup>2</sup> a tj j cos <sup>ϕ</sup> <sup>þ</sup> <sup>ϕ</sup><sup>c</sup> ð Þ (7)

<sup>1</sup> <sup>þ</sup> <sup>a</sup><sup>2</sup> j j<sup>t</sup> <sup>2</sup> � <sup>2</sup> a tj j cos <sup>ϕ</sup> <sup>þ</sup> <sup>ϕ</sup><sup>c</sup> ð Þ (8)

ð Þ <sup>1</sup> � a tj j <sup>2</sup> (9)

ð Þ <sup>1</sup> � a tj j <sup>2</sup> (10)

<sup>2</sup> <sup>þ</sup> <sup>κ</sup><sup>2</sup> <sup>¼</sup> <sup>1</sup> (5)

<sup>i</sup> <sup>β</sup><sup>2</sup> <sup>π</sup> <sup>R</sup> Et<sup>2</sup> � a ei<sup>ϕ</sup> Et<sup>2</sup> (6)

dent of frequency):

round-trip phase shift.

power in the ring, which is:

then Eqs. (7) and (8) becomes:

52

Then, the transmitted power to through port (3) becomes 0. The critical coupling is due to destructive interference between the direct beam and the crosscoupled beam of light, both after passing their respective coupling zone. A condition termed undercoupling does occur when the ring loss exceeds the coupling strength (t > a or equally: t <sup>2</sup> > αL, also κ<sup>2</sup> < αL). An overcoupling condition takes place, when the round-trip loss does not exceed the coupling strength (t < a or equally: t <sup>2</sup> < αL, also κ<sup>2</sup> > αL) and is the conventional mode of operation for an all-pass resonator. Figure 3(a) and (b) displays the effect of attenuation on the transmission and build-up factor for the overcoupled, critically coupled and undercoupled regimes. Experimental demonstration of the transmission characteristics in each of these regimes for a fiber all-pass filter ring resonator can be found in [22].

A good idea of the behavior of the system, when resonance condition accomplishes, can be obtained when one represents the T transmission function defined as the ratio Pt1/Pi<sup>1</sup> (i.e., the relationship between the through-port power and inputport power) versus normalized detuning ϕ, Figure 3(a). Starting from input field expressions and after some derivations, the following equation can be found:

$$T = \frac{P\_{t1}}{P\_{i1}} = \left| \frac{E\_{t1}}{E\_{i1}} \right|^2 = \frac{a^2 - 2|t|a\cos\phi + |t|^2}{1 - 2|t|a\cos\phi + |t|^2a^2} \tag{11}$$

Figure 3(a) represents the T transmission function for an all-pass ring resonator with t = 0.75 and varying loss a as a parameter (from a = 0.00 up to a = 0.99). It can be seen that, in this case, critical coupling occurs when a = t = 0.75; therefore, the T transmission function drops to zero when ϕ ¼ 0 for this a-factor value. As a

#### Figure 3.

T transmission function (a) and B build-up factor (b) of ring-resonator, both as a function of normalized detuning ϕ, (curves extracted from [21], with permission).

increases (up to 1) or decreases (to 0.00) from 0.75, the corresponding transmission curve arises its minimum point from 0. For a = 0.00, T transmission function reaches a constant value of 0.75<sup>2</sup> = 0.5625 for all values of ϕ. Another function that allows to obtain a good idea of the behavior of the system is the so-called B intensity build-up factor as a function of ϕ (normalized detuning), defined as:

$$\mathcal{B} = \frac{P\_{i2}}{P\_{i1}} = \left| \frac{E\_{i2}}{E\_{i1}} \right|^2 = \frac{\left(1 - |t|^2\right)a^2}{1 - 2|t|a\cos\phi + |t|^2a^2} \Rightarrow \lim\_{\phi = 2\pi m, a = 1} \tag{12}$$

$$\frac{\left(1 - |t|^2\right)a^2}{1 - 2|t|a\cos\phi + |t|^2a^2} = \frac{1 + |t|}{1 - |t|} \approx \frac{4}{\kappa^2}$$

B function can be thought as the ratio of the circulating light intensity inside the ring to incident light intensity. Figure 3(b) represents the B function (or intensity build-up factor) versus normalized detuning ϕ for the same all-pass ring than the one in Figure 3(a), when varying loss a is taken as a parameter (from a = 0.20 up to a = 0.99). The last result in Eq. (12) refers to case in which incident light is resonant with the ring (ϕ ¼ 2 mπ) and attenuation is negligible (a = 1). This happens when constructive interference is produced at the coupler port (4) between the coupled light entering the ring coming from waveguide and the recirculation light coming from input port (3). This way, the circulating intensity of light around the ring results reinforced. However, for this constructive interference to take place, a coherent source is essential at input port (1). A passive ring resonator under these conditions attains the maximum ratio of circulating power to incident power that can be achieved. Figure 4 shows two successive resonances versus normalized detuning with peak values given by right side of Eq. (12).

than in the bus, fiber ring resonators can be used for nonlinear optics applications

wavelengths corresponding to two successive resonant conditions). When

neff , approximate calculation of FSR Δλ is very simple [20], as follows:

FSR <sup>¼</sup> <sup>Δ</sup><sup>λ</sup> ¼ � <sup>2</sup><sup>π</sup>

Spectral response of through port of an all-pass fiber ring resonator (notch-type filter).

Ring resonator filters can be characterized by four parameters, which are generally used to describe their performance, namely, the free spectral range (FSR), the full width at half maximum (FWHM), the finesse F and the quality factor Q. The FSR Δλ is the distance between two successive resonance peaks (difference between the

neglecting the wavelength dependence of the fiber-optic effective refractive index

L

FSR <sup>¼</sup> <sup>Δ</sup><sup>λ</sup> <sup>¼</sup> <sup>λ</sup><sup>2</sup>

The second parameter in importance for describing the performance of the optical filter is the FWHM. It is defined as the full spectral width at half maximum (or 3 dB bandwidth 2 δλ). In the case of all-pass filter, FWHM must be redefined as full width at half depth (FWHD) and its value can be calculated as follows (in

> 1 � j jt ffiffiffiffiffi

Third parameter is the finesse F. The finesse is a figure of merit of the filter that can be defined as the ratio between FSR and FWHD. So that implementing this

j j<sup>t</sup> <sup>p</sup> ����! j j <sup>t</sup> <sup>≈</sup><sup>1</sup> <sup>≈</sup> <sup>λ</sup><sup>2</sup> <sup>κ</sup><sup>2</sup>

2π<sup>2</sup>ngR

This approximation is derived from the resonant condition next to a resonance found for the used propagation constant. A more accurate calculation is obtained taking into account the chromatic dispersion of the fiber. Then, the group refractive index of the fiber can be used instead of the effective index, so that the new

∂β ∂β � ��<sup>1</sup> <sup>≈</sup> <sup>λ</sup><sup>2</sup>

neff <sup>L</sup> (13)

ng <sup>L</sup> (14)

<sup>4</sup>π<sup>2</sup>ngR (15)

with low input intensities.

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

Figure 5.

expression for FRS to be:

55

absence of ring internal losses, that is, a = 1):

FWHD <sup>¼</sup> <sup>2</sup> δλ <sup>¼</sup> <sup>λ</sup><sup>2</sup>

The spectral characteristics of the all-pass fiber-ring resonator configuration (one waveguide, one resonator, as represented in Figure 2) as optical filter can better be viewed when spectral response of through port versus vacuum wavelength is represented. A notch type filter is obtained from observation of this resonance curve (Figure 5). For cross-coupling value of 10% (κ<sup>2</sup> <sup>¼</sup> <sup>0</sup>:10), the optical intensity in the ring can be 40 times higher than the intensity in the input port of the waveguide. This way, since the intensity in the fiber ring can be much higher

Figure 4. Build-up factor versus normalized detuning for all-pass fiber ring resonator (from [21], with permission).

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

#### Figure 5.

increases (up to 1) or decreases (to 0.00) from 0.75, the corresponding transmission curve arises its minimum point from 0. For a = 0.00, T transmission function reaches a constant value of 0.75<sup>2</sup> = 0.5625 for all values of ϕ. Another function that allows to obtain a good idea of the behavior of the system is the so-called B intensity

a2

a2 <sup>1</sup> � <sup>2</sup>j j<sup>t</sup> <sup>a</sup> cos <sup>ϕ</sup> <sup>þ</sup> j j<sup>t</sup> <sup>2</sup> <sup>a</sup><sup>2</sup> <sup>¼</sup> <sup>1</sup> <sup>þ</sup> j j<sup>t</sup>

B function can be thought as the ratio of the circulating light intensity inside the ring to incident light intensity. Figure 3(b) represents the B function (or intensity build-up factor) versus normalized detuning ϕ for the same all-pass ring than the one in Figure 3(a), when varying loss a is taken as a parameter (from a = 0.20 up to a = 0.99). The last result in Eq. (12) refers to case in which incident light is resonant with the ring (ϕ ¼ 2 mπ) and attenuation is negligible (a = 1). This happens when constructive interference is produced at the coupler port (4) between the coupled light entering the ring coming from waveguide and the recirculation light coming from input port (3). This way, the circulating intensity of light around the ring results reinforced. However, for this constructive interference to take place, a coherent source is essential at input port (1). A passive ring resonator under these conditions attains the maximum ratio of circulating power to incident power that can be achieved. Figure 4 shows two successive resonances versus normalized

The spectral characteristics of the all-pass fiber-ring resonator configuration (one waveguide, one resonator, as represented in Figure 2) as optical filter can better be viewed when spectral response of through port versus vacuum wavelength is represented. A notch type filter is obtained from observation of this resonance curve (Figure 5). For cross-coupling value of 10% (κ<sup>2</sup> <sup>¼</sup> <sup>0</sup>:10), the optical intensity in the ring can be 40 times higher than the intensity in the input port of the waveguide. This way, since the intensity in the fiber ring can be much higher

Build-up factor versus normalized detuning for all-pass fiber ring resonator (from [21], with permission).

) lim<sup>ϕ</sup>¼2πm,a¼<sup>1</sup>

<sup>1</sup> � j j<sup>t</sup> <sup>≈</sup> <sup>4</sup>

κ2

(12)

build-up factor as a function of ϕ (normalized detuning), defined as:

Interferometry - Recent Developments and Contemporary Applications

<sup>1</sup> � j j<sup>t</sup> <sup>2</sup>

<sup>1</sup> � j j<sup>t</sup> <sup>2</sup>

<sup>1</sup> � <sup>2</sup>j j<sup>t</sup> <sup>a</sup> cos <sup>ϕ</sup> <sup>þ</sup> j j<sup>t</sup> <sup>2</sup> <sup>a</sup><sup>2</sup>

<sup>B</sup> <sup>¼</sup> Pi<sup>2</sup> Pi<sup>1</sup>

Figure 4.

54

<sup>¼</sup> Ei<sup>2</sup> Ei<sup>1</sup> 

 

2 ¼

detuning with peak values given by right side of Eq. (12).

Spectral response of through port of an all-pass fiber ring resonator (notch-type filter).

than in the bus, fiber ring resonators can be used for nonlinear optics applications with low input intensities.

Ring resonator filters can be characterized by four parameters, which are generally used to describe their performance, namely, the free spectral range (FSR), the full width at half maximum (FWHM), the finesse F and the quality factor Q. The FSR Δλ is the distance between two successive resonance peaks (difference between the wavelengths corresponding to two successive resonant conditions). When neglecting the wavelength dependence of the fiber-optic effective refractive index neff , approximate calculation of FSR Δλ is very simple [20], as follows:

$$F \text{SR} = \Delta \lambda = -\frac{2\pi}{L} \left(\frac{\partial \beta}{\partial \beta}\right)^{-1} \approx \frac{\lambda^2}{n\_{\text{eff}}L} \tag{13}$$

This approximation is derived from the resonant condition next to a resonance found for the used propagation constant. A more accurate calculation is obtained taking into account the chromatic dispersion of the fiber. Then, the group refractive index of the fiber can be used instead of the effective index, so that the new expression for FRS to be:

$$F\text{SR} = \Delta\lambda = \frac{\lambda^2}{n\_\text{g}L} \tag{14}$$

The second parameter in importance for describing the performance of the optical filter is the FWHM. It is defined as the full spectral width at half maximum (or 3 dB bandwidth 2 δλ). In the case of all-pass filter, FWHM must be redefined as full width at half depth (FWHD) and its value can be calculated as follows (in absence of ring internal losses, that is, a = 1):

$$\text{FWHD} = 2\,\delta\lambda = \frac{\lambda^2}{2\pi^2 n\_\text{g}R} \xrightarrow[\sqrt{|t|}]{} \frac{|t|\approx 1}{\sqrt{|t|}} \approx \frac{\lambda^2 \kappa^2}{4\pi^2 n\_\text{g}R} \tag{15}$$

Third parameter is the finesse F. The finesse is a figure of merit of the filter that can be defined as the ratio between FSR and FWHD. So that implementing this

definition and in absence of ring internal losses, the resulting expression for the finesse of an all-pass fiber resonator is [23]:

$$\mathcal{F}\_{\;\;\;\;\;\;F\;\;\;\;\;\;\;\;\;\;\mathcal{F}\;\;\;\;\;\;\mathcal{F}\;\;\;\;\;\;\;\;\;\mathcal{F}\tag{16}$$

Et<sup>1</sup> ¼ t<sup>1</sup> þ

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

phase, respectively. It is <sup>a</sup> <sup>¼</sup> <sup>a</sup><sup>2</sup>

rms amplitude value is given by:

drop port is given by:

port) is given by:

<sup>T</sup> through <sup>¼</sup> Et<sup>1</sup>

input port), will be:

<sup>T</sup> drop <sup>¼</sup> Et<sup>2</sup>

57

Ei<sup>1</sup> � � � �

� � � �

2

<sup>¼</sup> <sup>1</sup> � <sup>t</sup>

Ei<sup>1</sup> � � � �

� � � �

2 <sup>¼</sup> <sup>t</sup> 2

2 4

t<sup>2</sup> κ<sup>1</sup> κ<sup>2</sup> a<sup>2</sup>

<sup>1</sup>=<sup>2</sup> <sup>e</sup><sup>j</sup> <sup>θ</sup> j j <sup>t</sup><sup>1</sup> <sup>2</sup>

<sup>1</sup> � <sup>t</sup><sup>1</sup> <sup>t</sup><sup>2</sup> <sup>a</sup><sup>2</sup>

þj j <sup>κ</sup><sup>1</sup> <sup>2</sup> ð Þ

3

<sup>5</sup>Ei<sup>1</sup> <sup>¼</sup> <sup>t</sup><sup>1</sup> � <sup>t</sup><sup>2</sup> a e<sup>j</sup> <sup>θ</sup>

<sup>1</sup>=<sup>2</sup> and θ ¼ 2θ1<sup>=</sup>2. The mode amplitude in the ring has

2 1 � � <sup>1</sup> � <sup>t</sup>

ð Þ 1 � t<sup>1</sup> t<sup>2</sup> a

<sup>1</sup> � <sup>t</sup><sup>1</sup> <sup>t</sup><sup>2</sup> a ej <sup>θ</sup> (19)

2 2 � �a

<sup>2</sup> (20)

4 t 2 <sup>1</sup>�t<sup>2</sup> ð Þ<sup>2</sup> sin <sup>2</sup> <sup>ϕ</sup>

<sup>1</sup> <sup>þ</sup> <sup>4</sup> <sup>t</sup><sup>2</sup>

<sup>1</sup> <sup>þ</sup> <sup>4</sup> <sup>t</sup><sup>2</sup>

(21)

2 � �

> 2 � �

(22)

2 � �

(23)

<sup>1</sup>�t<sup>2</sup> ð Þ<sup>2</sup> sin <sup>2</sup> <sup>ϕ</sup>

1

<sup>1</sup>�t<sup>2</sup> ð Þ<sup>2</sup> sin <sup>2</sup> <sup>ϕ</sup>

<sup>1</sup> � <sup>t</sup><sup>1</sup> <sup>t</sup><sup>2</sup> a ej <sup>θ</sup> (18)

<sup>1</sup>=<sup>2</sup> ej <sup>θ</sup>

In this equation, a1=<sup>2</sup> and θ1=<sup>2</sup> are used, which are the half round-trip loss and

to pass the second coupler to become the new dropped mode rms amplitude Et2. Its

Et<sup>2</sup> <sup>¼</sup> <sup>κ</sup><sup>1</sup> <sup>κ</sup><sup>2</sup> <sup>a</sup>1=2e<sup>j</sup> <sup>θ</sup>1=<sup>2</sup>

At resonance condition (θ ¼ m2π, m integer), the output intensity from the

The through port mode amplitude Et<sup>1</sup> will be zero at resonance for identical symmetrical couplers (t<sup>1</sup> = t2), when a = 1 (zero loss). This means that on resonance, all input optical power will be coupled to resonator and could be extracted by drop port. The value a = 1 (zero loss) can only be achieved by a given value of gain incorporated in the resonator ring to compensate the waveguide losses. In a pure passive ring resonator with defined waveguide constituents, the loss coefficient a is constant. By adjusting coupling parameters t<sup>1</sup> and t<sup>2</sup> to the loss coefficient a, minimum intensity can be achieved (ideally, Pt<sup>1</sup> ≈0). From Eq. (18), one can obtain:

> <sup>a</sup> <sup>¼</sup> <sup>t</sup><sup>1</sup> t2

Taking into account, Eq. (19), the transfer function of drop port (relative to

This way, the transmission spectral response of through port corresponds to an approximately Lorentzian curve centered on a resonance, where a narrow frequency band has been extracted (see Figure 7). These extracted band do appear (in inverted shape) in the spectrum response of drop port, as it can also be seen in same figure. In the case of add-drop filter, if attenuation is negligible (a≈ 1), critical coupling occurs at symmetric coupling (κ<sup>1</sup> ¼ κ<sup>2</sup> ¼ κ). For non-negligible single-pass resonator loss, critical coupling occurs when single-pass ring losses match the

2 2 � �a

<sup>2</sup> <sup>a</sup><sup>2</sup> � <sup>2</sup> <sup>t</sup><sup>1</sup> <sup>t</sup><sup>2</sup> <sup>a</sup> cos <sup>ϕ</sup> <sup>þ</sup> <sup>t</sup>

1 � 2t<sup>1</sup> t<sup>2</sup> a cos ϕ þ ð Þ t<sup>1</sup> t<sup>2</sup> a

2 1 � � <sup>1</sup> � <sup>t</sup>

1 � 2 t<sup>1</sup> t<sup>2</sup> a cos ϕ þ ð Þ t<sup>1</sup> t<sup>2</sup> a

According with Eq. (18), the transfer function of through port (relative to input

2 1

<sup>2</sup> ��������! <sup>a</sup>¼<sup>1</sup>, <sup>t</sup>1¼t2¼<sup>t</sup>

<sup>2</sup> ��������! <sup>a</sup>¼<sup>1</sup>, <sup>t</sup>1¼t2¼<sup>t</sup>

Pt<sup>2</sup>�resonance <sup>¼</sup> Et<sup>2</sup>�resonance j j<sup>2</sup> <sup>¼</sup> <sup>1</sup> � <sup>t</sup>

Finally, the fourth parameter is the quality factor Q. It represents a measure of the sharpness of the resonance. It can be defined as the ratio of the operation wavelength and the resonance width (FWHD). Its calculation is as follows:

$$\mathcal{Q} = \frac{\lambda}{\text{FWHM}} \approx \frac{\lambda}{\frac{\lambda^2 \kappa^2}{4\pi^2 n\_\text{fR}}} = \frac{\lambda}{\frac{\lambda^2 \kappa^2}{2\pi n\_\text{fL}}} = \frac{2\pi n\_\text{gL}}{\lambda \kappa^2} = \frac{n\_\text{gL}}{\lambda} \,\,\, \mathcal{F} \tag{17}$$

The quality factor Q can also be regarded as the total energy stored in the ring divided by the energy lost per optical cycle.

## 3. Analysis of the fiber-optic add-drop filter ring resonator configuration

The basic fiber-optic add-drop filter consists of one input-waveguide, one output-waveguide and one fiber-ring resonator such as the one represented in Figure 6. Both of the waveguides are side-coupled to ring resonator. The four ports of structure are defined as input port, throughput port, drop port and add port. Viewed from drop port acting as the optical output, this configuration behaves as a narrow-band amplitude filter that can add or drop a spectral band from an optical input. This configuration is mathematically equivalent to a well-known Fabry-Pérot interferometer. Therefore, the equations for transmission coefficients are easy to obtain. For simplification, Ei<sup>1</sup> is defined with unitary rms amplitude. Then, the throughput Et<sup>1</sup> mode rms amplitude in the input-waveguide is given by:

Figure 6. Add-drop fiber-optic ring resonator configuration (for symmetrical coupling: κ<sup>1</sup> ≈ κ<sup>2</sup> ¼ κ and t<sup>1</sup> ≈ t<sup>2</sup> ¼ t).

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

definition and in absence of ring internal losses, the resulting expression for the

λ2 2 π ng R λ<sup>2</sup> κ<sup>2</sup> 4π<sup>2</sup>ngR

<sup>¼</sup> <sup>2</sup><sup>π</sup>

<sup>¼</sup> <sup>2</sup><sup>π</sup> ng <sup>L</sup>

λκ<sup>2</sup> <sup>¼</sup> ng <sup>L</sup>

<sup>κ</sup><sup>2</sup> (16)

<sup>λ</sup> <sup>F</sup> (17)

FWHD <sup>≈</sup>

the sharpness of the resonance. It can be defined as the ratio of the operation wavelength and the resonance width (FWHD). Its calculation is as follows:

> λ<sup>2</sup> κ<sup>2</sup> 4π<sup>2</sup>ngR

3. Analysis of the fiber-optic add-drop filter ring resonator

Finally, the fourth parameter is the quality factor Q. It represents a measure of

<sup>¼</sup> <sup>λ</sup> λ<sup>2</sup> κ<sup>2</sup> 2 π ng L

The quality factor Q can also be regarded as the total energy stored in the ring

The basic fiber-optic add-drop filter consists of one input-waveguide, one output-waveguide and one fiber-ring resonator such as the one represented in Figure 6. Both of the waveguides are side-coupled to ring resonator. The four ports of structure are defined as input port, throughput port, drop port and add port. Viewed from drop port acting as the optical output, this configuration behaves as a narrow-band amplitude filter that can add or drop a spectral band from an optical input. This configuration is mathematically equivalent to a well-known Fabry-Pérot interferometer. Therefore, the equations for transmission coefficients are easy to obtain. For simplification, Ei<sup>1</sup> is defined with unitary rms amplitude. Then, the throughput Et<sup>1</sup> mode rms amplitude in the input-waveguide is given by:

Add-drop fiber-optic ring resonator configuration (for symmetrical coupling: κ<sup>1</sup> ≈ κ<sup>2</sup> ¼ κ and t<sup>1</sup> ≈ t<sup>2</sup> ¼ t).

<sup>F</sup> <sup>¼</sup> FSR

Interferometry - Recent Developments and Contemporary Applications

finesse of an all-pass fiber resonator is [23]:

<sup>Q</sup> <sup>¼</sup> <sup>λ</sup>

divided by the energy lost per optical cycle.

configuration

Figure 6.

56

FWHD <sup>≈</sup> <sup>λ</sup>

$$E\_{t1} = \left[t\_1 + \frac{t\_2 \kappa\_1 \kappa\_2 a\_{1/2}^2 e^{j\theta \left(\left|t\_1\right|^2 + \left|\kappa\_1\right|^2\right)}}{\mathbf{1} - t\_1 t\_2 a\_{1/2}^2 e^{j\theta}}\right] E\_{i1} = \frac{t\_1 - t\_2 a\_1 e^{j\theta}}{\mathbf{1} - t\_1 t\_2 a\_2 e^{j\theta}}\tag{18}$$

In this equation, a1=<sup>2</sup> and θ1=<sup>2</sup> are used, which are the half round-trip loss and phase, respectively. It is <sup>a</sup> <sup>¼</sup> <sup>a</sup><sup>2</sup> <sup>1</sup>=<sup>2</sup> and θ ¼ 2θ1<sup>=</sup>2. The mode amplitude in the ring has to pass the second coupler to become the new dropped mode rms amplitude Et2. Its rms amplitude value is given by:

$$E\_{t2} = \frac{\kappa\_1 \kappa\_2 a\_{1/2} e^{j\theta\_{1/2}}}{1 - t\_1 t\_2 a\_1 e^{j\theta}} \tag{19}$$

At resonance condition (θ ¼ m2π, m integer), the output intensity from the drop port is given by:

$$P\_{t2-reonance} = \left| E\_{t2-reonance} \right|^2 = \frac{\left(\mathbf{1} - t\_1^2\right)\left(\mathbf{1} - t\_2^2\right)a}{\left(\mathbf{1} - t\_1t\_2a\right)^2} \tag{20}$$

The through port mode amplitude Et<sup>1</sup> will be zero at resonance for identical symmetrical couplers (t<sup>1</sup> = t2), when a = 1 (zero loss). This means that on resonance, all input optical power will be coupled to resonator and could be extracted by drop port. The value a = 1 (zero loss) can only be achieved by a given value of gain incorporated in the resonator ring to compensate the waveguide losses. In a pure passive ring resonator with defined waveguide constituents, the loss coefficient a is constant. By adjusting coupling parameters t<sup>1</sup> and t<sup>2</sup> to the loss coefficient a, minimum intensity can be achieved (ideally, Pt<sup>1</sup> ≈0). From Eq. (18), one can obtain:

$$a = \frac{t\_1}{t\_2} \tag{21}$$

According with Eq. (18), the transfer function of through port (relative to input port) is given by:

$$\left| \mathcal{T}\_{\text{through}} = \frac{\left| E\_{t1} \right|^2}{\left| E\_{i1} \right|} = \frac{t\_2^2 a^2 - 2 \, t\_1 t\_2 a \cos \phi + t\_1^2}{1 - 2t\_1 t\_2 a \cos \phi + \left( t\_1 t\_2 a \right)^2} \frac{\frac{4t^2}{\left( 1 - t^2 \right)^2} \sin^2 \left( \frac{\phi}{2} \right)}{1 + \frac{4t^2}{\left( 1 - t^2 \right)^2} \sin^2 \left( \frac{\phi}{2} \right)} \tag{22}$$

Taking into account, Eq. (19), the transfer function of drop port (relative to input port), will be:

$$\left| \mathcal{T}\_{drop} = \left| \frac{E\_{t2}}{E\_{t1}} \right|^2 = \frac{\left(1 - t\_1^2\right)\left(1 - t\_2^2\right)a}{1 - 2t\_1t\_2a\cos\phi + \left(t\_1t\_2a\right)^2} \xrightarrow[a = \text{l\_1}\ t\_1 = t\_2 = t} \quad \frac{1}{1 + \frac{4t^2}{\left(1 - t^2\right)^2}\sin^2\left(\frac{\theta}{2}\right)}\tag{23}$$

This way, the transmission spectral response of through port corresponds to an approximately Lorentzian curve centered on a resonance, where a narrow frequency band has been extracted (see Figure 7). These extracted band do appear (in inverted shape) in the spectrum response of drop port, as it can also be seen in same figure. In the case of add-drop filter, if attenuation is negligible (a≈ 1), critical coupling occurs at symmetric coupling (κ<sup>1</sup> ¼ κ<sup>2</sup> ¼ κ). For non-negligible single-pass resonator loss, critical coupling occurs when single-pass ring losses match the

4. Numerical simulation of fiber-optic all-pass filter ring resonator

As announced at the beginning of this chapter, this section will deal specifically with the numerical simulation of a typical FOPRRI in all-pass filter configuration. First, effective index method (EIM) [25, 26] will be applied to this structure in which main constituents are made of a standard single-mode optical fiber (coupler and ring resonator). Secondly, coupled mode theory (CMT) [27–33] will be

performed to the same structure. In third time, finite element method (FEM) [34] is performed. Finally, 3D-FDTD (finite difference time domain) [35, 36] was applied to same structure. The results of all these numerical simulations are compared to

Optical layout scheme for EIM and CMT numerical simulation is represented in Figure 8. OptSim® Circuit tool (from Synopsys™ Software) was applied to obtain the transmission power spectrum of this resonant structure, which takes the effective index method (EIM) as a numerical simulation algorithm. A side-polished fiber-optic 4 4 directional coupler is used to couple the waveguide fiber to fiber ring resonator (here implemented by a series of two waveguides, 1 and 2, respectively). A standard SMF-28 (Corning™ Glass) single-mode optical fiber is assumed for both, waveguide and ring resonator as constituent materials. Optical parameters of this fiber are: basis material = fused-silica (core = boron-silicate, cladding = silica), step refractive index, core diameter = 8.2 μm, cladding diameter = 125 μm, effective refractive indexes: ncore = 1.467, ncladding = 1.460, numerical aperture (NA) = 0.143, maximum attenuation = 0.35 dB/km (at 1310 nm wavelength) and h-

Normalized Gaussian profile continuous wave (CW) optical field of 1310 nm vacuum wavelength was launched to the system from a wideband tunable laser source. To obtain the pass-through port power transmission spectrum, wavelength sweeping of laser source was made in order to fix each wavelength value for its respective measured optical power. The computation is iterated 30 times to obtain the steadystate solution. Obtained transmission power spectrum is represented in Figure 9.

All-pass FOPRRI configuration: optical layout scheme for EIM and CMT numerical simulations.

. Gap between waveguide and ring resonator is 200 nm.

each other to see which one best match each specific application.

configuration

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

parameter = 2 <sup>10</sup><sup>6</sup> <sup>m</sup><sup>1</sup>

Figure 8.

59

Figure 7. Spectral response of a fiber-optic add-drop filter (adapted from [24], with permission).

coupling, so it verified the following condition: t<sup>2</sup> a ¼ t1. Transfer functions, FWHM and extinction ratio (ER) relative to all-pass filter, through and dropt ports of adddrop filter are plotted.

It is worth to make an interpretation and comparison of the add-drop filter configuration with respect to the all-pass filter one. In the upper coupler (Figure 6), interference takes place between the coupled input field and ring circulating field. However, the lower coupler does not experience any interference (provided that excitation is from the upper waveguide only). The lower coupler may thus be viewed simply as a tap waveguide that leaks power out of the ring into the drop port. For an add-drop filter, it is desirable to operate at critical coupling for achieving the complete extinguishment of a frequency band in the through port. In this case, the sum of all losses incurred in the resonator, even the loss at the out-coupled drop port, must be taken into account. Calculations of FSR, FWHM, B build-up factor, finesse F and quality factor Q are given below according to the results found in Refs. [21, 22].

$$F \text{SR} = \Delta\lambda = \frac{\lambda^2}{n\_{\text{eff}}L} \tag{24}$$

$$F \text{WHD} \approx \frac{\lambda^2 \kappa^2}{\pi L n\_{\text{eff}}} \tag{25}$$

$$\mathcal{F} = \frac{F \text{SR}}{F \text{WHD}} \approx \frac{\pi}{\kappa^2} \tag{26}$$

$$\mathcal{Q} = \frac{\lambda}{\text{FWHM}} \approx \frac{\pi n\_{\text{g}} L}{\lambda \kappa^2} = \frac{n\_{\text{g}} L}{\lambda} \mathcal{F} \tag{27}$$

where the approximations again apply to symmetrically coupled, low-loss and weakly coupled microrings. Note that, comparing Eqs. (27) and (17), with the same optical parameters, add-drop filter has the half value of Q-factor (and also the finesse F) than that the all-pass one. This result could be expected since an all-pass filter configuration has half the loss due to external coupling than an add-drop one.

## 4. Numerical simulation of fiber-optic all-pass filter ring resonator configuration

As announced at the beginning of this chapter, this section will deal specifically with the numerical simulation of a typical FOPRRI in all-pass filter configuration. First, effective index method (EIM) [25, 26] will be applied to this structure in which main constituents are made of a standard single-mode optical fiber (coupler and ring resonator). Secondly, coupled mode theory (CMT) [27–33] will be performed to the same structure. In third time, finite element method (FEM) [34] is performed. Finally, 3D-FDTD (finite difference time domain) [35, 36] was applied to same structure. The results of all these numerical simulations are compared to each other to see which one best match each specific application.

Optical layout scheme for EIM and CMT numerical simulation is represented in Figure 8. OptSim® Circuit tool (from Synopsys™ Software) was applied to obtain the transmission power spectrum of this resonant structure, which takes the effective index method (EIM) as a numerical simulation algorithm. A side-polished fiber-optic 4 4 directional coupler is used to couple the waveguide fiber to fiber ring resonator (here implemented by a series of two waveguides, 1 and 2, respectively). A standard SMF-28 (Corning™ Glass) single-mode optical fiber is assumed for both, waveguide and ring resonator as constituent materials. Optical parameters of this fiber are: basis material = fused-silica (core = boron-silicate, cladding = silica), step refractive index, core diameter = 8.2 μm, cladding diameter = 125 μm, effective refractive indexes: ncore = 1.467, ncladding = 1.460, numerical aperture (NA) = 0.143, maximum attenuation = 0.35 dB/km (at 1310 nm wavelength) and hparameter = 2 <sup>10</sup><sup>6</sup> <sup>m</sup><sup>1</sup> . Gap between waveguide and ring resonator is 200 nm. Normalized Gaussian profile continuous wave (CW) optical field of 1310 nm vacuum wavelength was launched to the system from a wideband tunable laser source. To obtain the pass-through port power transmission spectrum, wavelength sweeping of laser source was made in order to fix each wavelength value for its respective measured optical power. The computation is iterated 30 times to obtain the steadystate solution. Obtained transmission power spectrum is represented in Figure 9.

Figure 8. All-pass FOPRRI configuration: optical layout scheme for EIM and CMT numerical simulations.

coupling, so it verified the following condition: t<sup>2</sup> a ¼ t1. Transfer functions, FWHM and extinction ratio (ER) relative to all-pass filter, through and dropt ports of add-

Spectral response of a fiber-optic add-drop filter (adapted from [24], with permission).

Interferometry - Recent Developments and Contemporary Applications

It is worth to make an interpretation and comparison of the add-drop filter configuration with respect to the all-pass filter one. In the upper coupler (Figure 6), interference takes place between the coupled input field and ring circulating field. However, the lower coupler does not experience any interference (provided that excitation is from the upper waveguide only). The lower coupler may thus be viewed simply as a tap waveguide that leaks power out of the ring into the drop port. For an add-drop filter, it is desirable to operate at critical coupling for achieving the complete extinguishment of a frequency band in the through port. In this case, the sum of all losses incurred in the resonator, even the loss at the out-coupled drop port, must be taken into account. Calculations of FSR, FWHM, B build-up factor, finesse F and quality factor Q are given below according to the results found

FSR <sup>¼</sup> Δλ <sup>¼</sup> <sup>λ</sup><sup>2</sup>

FWHD <sup>≈</sup> <sup>λ</sup><sup>2</sup> <sup>κ</sup><sup>2</sup>

<sup>F</sup> <sup>¼</sup> FSR

FWHD <sup>≈</sup> <sup>π</sup> ng <sup>L</sup>

where the approximations again apply to symmetrically coupled, low-loss and weakly coupled microrings. Note that, comparing Eqs. (27) and (17), with the same optical parameters, add-drop filter has the half value of Q-factor (and also the finesse F) than that the all-pass one. This result could be expected since an all-pass filter configuration has half the loss due to external coupling than an add-drop one.

<sup>Q</sup> <sup>¼</sup> <sup>λ</sup>

π L neff

<sup>λ</sup>κ<sup>2</sup> <sup>¼</sup> ng <sup>L</sup>

FWHD <sup>≈</sup> <sup>π</sup>

neff <sup>L</sup> (24)

<sup>κ</sup><sup>2</sup> (26)

<sup>λ</sup> <sup>F</sup> (27)

(25)

drop filter are plotted.

Figure 7.

in Refs. [21, 22].

58

Figure 9. All-pass FOPRRI configuration: transmission power spectrum.

also takes into account the effective index of the ring resonator material, so the expected results from it will be improved compared with those obtained with the EIM method alone. Further in-depth readings on BPM method can be found in Refs. [37–42]. Figure 11 shows the result of numerical simulation of all-pass ring resona-

FEM numerical simulation (from FemSIM® tool) for all-pass FOPRRI configuration (resonance at 1285 nm).

CMT method numerical simulation result (from BeamPROP® tool) for all-pass FOPRRI configuration

Finite element method (FEM) was applied to this structure by means of FemSIM® tool (Synopsys™ Software) [34]. This method also takes into account the effective index of the ring resonator material, so the expected results from it will be improved compared with those obtained with the EIM method alone. Additional treatment on the FEM method can be found in the Refs. [43–47]. The result of FemSIM® of all-pass ring resonator analyzed here is represented in Figure 12.

FullWAVE® (Synopsys™ Software) [36] is a full-vectorial numerical computation tool that allows the simulation of light wave propagation in an arbitrary geometry system. The simulation engine is based on the finite difference time domain (FDTD) method. The FDTD algorithm implements rigorous and complete timespace solution to Maxwell's equations. This method offers an efficient field wave propagation solution in a photonic-integrated structure when other methods such as beam propagation method (BPM) or finite element method (FEM) cannot cope with the structure geometry or do not provide a precise solution. Additional treatment on the FDTD method can be found in Refs. [48–52]. Figure 13 represents the result of FDTD numerical simulation of all-pass ring resonator studied here.

tor considered here.

Figure 11.

Figure 12.

61

(resonance at 1285 nm).

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

Figure 10. All-pass FOPRRI configuration: expanded transmission power spectrum around 1295 nm.

As it can be seen on expanded transmission spectrum showed in Figure 10, a resonance condition is achieved for 1285 nm wavelength. Since the effective refractive index of ring resonator is 1.467 and its optical path length is 51.420 μm, a FSR of 21.89 nm is obtained as it can also be observed in this same figure. This value agrees accurately with theoretical value as calculated from Eq. (13).

On the other hand, measured value of FWHM is approximately 2 nm. This value is also in frank agreement with value calculated analytically by means of Eq. (15), taking into account 8.184 μm as medium ring resonator radius. Then, a calculation from Eq. (16) gives the value of 10.95 for F finesse. As a consequence of last, calculation from Eq. (17) yields 642.80 for the Q quality-factor of the system.

Coupled mode theory (CMT), which uses BPM method, was applied to this structure by means of BeamPROP® tool (Synopsys™ Software) [32]. This method Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

Figure 11.

CMT method numerical simulation result (from BeamPROP® tool) for all-pass FOPRRI configuration (resonance at 1285 nm).

#### Figure 12.

FEM numerical simulation (from FemSIM® tool) for all-pass FOPRRI configuration (resonance at 1285 nm).

also takes into account the effective index of the ring resonator material, so the expected results from it will be improved compared with those obtained with the EIM method alone. Further in-depth readings on BPM method can be found in Refs. [37–42]. Figure 11 shows the result of numerical simulation of all-pass ring resonator considered here.

Finite element method (FEM) was applied to this structure by means of FemSIM® tool (Synopsys™ Software) [34]. This method also takes into account the effective index of the ring resonator material, so the expected results from it will be improved compared with those obtained with the EIM method alone. Additional treatment on the FEM method can be found in the Refs. [43–47]. The result of FemSIM® of all-pass ring resonator analyzed here is represented in Figure 12.

FullWAVE® (Synopsys™ Software) [36] is a full-vectorial numerical computation tool that allows the simulation of light wave propagation in an arbitrary geometry system. The simulation engine is based on the finite difference time domain (FDTD) method. The FDTD algorithm implements rigorous and complete timespace solution to Maxwell's equations. This method offers an efficient field wave propagation solution in a photonic-integrated structure when other methods such as beam propagation method (BPM) or finite element method (FEM) cannot cope with the structure geometry or do not provide a precise solution. Additional treatment on the FDTD method can be found in Refs. [48–52]. Figure 13 represents the result of FDTD numerical simulation of all-pass ring resonator studied here.

As it can be seen on expanded transmission spectrum showed in Figure 10, a resonance condition is achieved for 1285 nm wavelength. Since the effective refractive index of ring resonator is 1.467 and its optical path length is 51.420 μm, a FSR of 21.89 nm is obtained as it can also be observed in this same figure. This value agrees

On the other hand, measured value of FWHM is approximately 2 nm. This value is also in frank agreement with value calculated analytically by means of Eq. (15), taking into account 8.184 μm as medium ring resonator radius. Then, a calculation from Eq. (16) gives the value of 10.95 for F finesse. As a consequence of last, calculation from Eq. (17) yields 642.80 for the Q quality-factor of the system. Coupled mode theory (CMT), which uses BPM method, was applied to this structure by means of BeamPROP® tool (Synopsys™ Software) [32]. This method

accurately with theoretical value as calculated from Eq. (13).

All-pass FOPRRI configuration: expanded transmission power spectrum around 1295 nm.

Figure 9.

Figure 10.

60

All-pass FOPRRI configuration: transmission power spectrum.

Interferometry - Recent Developments and Contemporary Applications

#### Figure 13. 2D-FTDT numerical simulation results for an all-pass FOPRRI (resonance at 1285 nm).


lower waveguides and ring resonator as constituent material. Gap between waveguide and ring resonator is 200 nm. Normalized Gaussian profile continuous wave (CW) optical field of 1310 nm vacuum wavelength was launched to the system from a wideband tunable laser source. To obtain the power transmission spectra of pass-through and drop ports, wavelength sweeping of laser source was made in order to fix each wavelength value for its respective measured optical power. The computation is iterated 30 times to obtain the steady-state solution. The resulting transmission power spectra are depicted in Figure 15(a)–(c), respectively. As it can be seen on expanded drop-port transmission spectrum showed in Figure 10(c), a resonance condition is achieved for 1285 nm wavelength. Since the effective refractive index of ring resonator is 1.467 and its optical path length is 51.420 μm, a FSR of 21.89 nm is obtained as it can also be observed in this same figure. This value agrees accurately with theoretical value as calculated from Eq. (24). On the other hand, measured value of FWHM is approximately 4 nm. This value is also in frank

Add-drop FOPRRI configuration: optical layout scheme for EIM and CMT numerical simulations.

Figure 14.

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

63

agreement with value calculated analytically by means of Eq. (25).

for add-drop configuration must be half than those for the all-pass ones.

Then, a calculation from Eq. (26) gives the value of 5.475 for F finesse. As a consequence of last, calculation from Eq. (27) yields 321.40 for the Q quality-factor of the system. These results are the expected, taking into account that for the same coupling coefficients, waveguide and ring resonator geometry, the values F and Q

Comments about fundamentals and particularities of BPM, FEM and FDTD methods are given in the following. For an optical waveguide as cylindrical optical fiber, where light travels almost completely confined to its core medium, propagation mechanism can be simply described by means of the four Maxwell's equations. Then, BPM, FEM and FDTD methods assume that the solution to Maxwell's equations follows the slowly varying envelope approximation. This means that the field's amplitude variation is very small as compared to the variation in the phase and propagation distance. BPM is a fast method, which is applicable to a system whose physical structure and properties vary slowly with the geometry. The FDTD method uses the standard finite difference method to solve the Maxwell's equations. It is also

#### Table 1.

Comparison of numerical simulation methods (EIM, CMT, FEM and FDTD) for all-pass/add-drop filters.

Table 1 shows a comparison of the results of the four methods used (EIM, CMT, FEM and FDTD) for all-pass filter and add-drop filter configurations, taking into account their suitability for accuracy and calculation time. In view of transmission spectrum of the all-pass filter configuration, this can find applications such as optical band-pass filter, optical delay line, optical switch and active ring modulator.

## 5. Numerical simulation of fiber-optic add-drop filter ring resonator configuration

As realized on all-pass configuration, this section will be dedicated to numerical simulation of a typical FOPRRI in add-drop filter configuration. Then, successive numerical simulations by applying effective index method (EIM), coupled mode theory (CMT), finite element method (FEM) and 3D finite difference time domain (FDTD) will be performed to this configuration. Finally, the comparison of all these numerical simulations will be carried out to verify, which is the most appropriate for each type of application.

Optical layout scheme for EIM and CMT numerical simulation is represented in Figure 14. OptSim® Circuit tool (from Synopsys™ Software) with the effective index method (EIM) as computational algorithm was applied here to obtain the transmission power spectrum of resonant structure. The configuration analyzed can be obtained from the one considered in Section 4 by adding a second fiber 4 4 directional coupler and a second side-coupled fiber waveguide. Also, standard SMF-28 (Corning™ Glass) single-mode optical fiber is assumed for both, upper and

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

Figure 14. Add-drop FOPRRI configuration: optical layout scheme for EIM and CMT numerical simulations.

lower waveguides and ring resonator as constituent material. Gap between waveguide and ring resonator is 200 nm. Normalized Gaussian profile continuous wave (CW) optical field of 1310 nm vacuum wavelength was launched to the system from a wideband tunable laser source. To obtain the power transmission spectra of pass-through and drop ports, wavelength sweeping of laser source was made in order to fix each wavelength value for its respective measured optical power. The computation is iterated 30 times to obtain the steady-state solution. The resulting transmission power spectra are depicted in Figure 15(a)–(c), respectively. As it can be seen on expanded drop-port transmission spectrum showed in Figure 10(c), a resonance condition is achieved for 1285 nm wavelength. Since the effective refractive index of ring resonator is 1.467 and its optical path length is 51.420 μm, a FSR of 21.89 nm is obtained as it can also be observed in this same figure. This value agrees accurately with theoretical value as calculated from Eq. (24). On the other hand, measured value of FWHM is approximately 4 nm. This value is also in frank agreement with value calculated analytically by means of Eq. (25).

Then, a calculation from Eq. (26) gives the value of 5.475 for F finesse. As a consequence of last, calculation from Eq. (27) yields 321.40 for the Q quality-factor of the system. These results are the expected, taking into account that for the same coupling coefficients, waveguide and ring resonator geometry, the values F and Q for add-drop configuration must be half than those for the all-pass ones.

Comments about fundamentals and particularities of BPM, FEM and FDTD methods are given in the following. For an optical waveguide as cylindrical optical fiber, where light travels almost completely confined to its core medium, propagation mechanism can be simply described by means of the four Maxwell's equations. Then, BPM, FEM and FDTD methods assume that the solution to Maxwell's equations follows the slowly varying envelope approximation. This means that the field's amplitude variation is very small as compared to the variation in the phase and propagation distance. BPM is a fast method, which is applicable to a system whose physical structure and properties vary slowly with the geometry. The FDTD method uses the standard finite difference method to solve the Maxwell's equations. It is also

Table 1 shows a comparison of the results of the four methods used (EIM, CMT, FEM and FDTD) for all-pass filter and add-drop filter configurations, taking into account their suitability for accuracy and calculation time. In view of transmission spectrum of the all-pass filter configuration, this can find applications such as optical band-pass filter, optical delay line, optical switch and active ring modulator.

Comparison of numerical simulation methods (EIM, CMT, FEM and FDTD) for all-pass/add-drop filters.

2D-FTDT numerical simulation results for an all-pass FOPRRI (resonance at 1285 nm).

Interferometry - Recent Developments and Contemporary Applications

Comparison EIM CMT (BPM) FEM FDTD All-pass filter Good Very good Outstanding Excellent Add-drop filter Good Very good Outstanding Excellent

5. Numerical simulation of fiber-optic add-drop filter ring resonator

As realized on all-pass configuration, this section will be dedicated to numerical simulation of a typical FOPRRI in add-drop filter configuration. Then, successive numerical simulations by applying effective index method (EIM), coupled mode theory (CMT), finite element method (FEM) and 3D finite difference time domain (FDTD) will be performed to this configuration. Finally, the comparison of all these numerical simulations will be carried out to verify, which is the most appropriate

Optical layout scheme for EIM and CMT numerical simulation is represented in Figure 14. OptSim® Circuit tool (from Synopsys™ Software) with the effective index method (EIM) as computational algorithm was applied here to obtain the transmission power spectrum of resonant structure. The configuration analyzed can be obtained from the one considered in Section 4 by adding a second fiber 4 4 directional coupler and a second side-coupled fiber waveguide. Also, standard SMF-28 (Corning™ Glass) single-mode optical fiber is assumed for both, upper and

configuration

Figure 13.

Table 1.

for each type of application.

62

FDTD method allow as broad as the Maxwell's equations they embody in discretized form, the basic formulation using centered finite difference expressions, which can treat virtually any material type and geometry and the basic resources needed to implement the technique. The FDTD method treats transients in the time domain (e.g., pulses), and it is applicable over the computationally unpredictable resonance region in which a wavelength is comparable with to the object interaction size. The FDTD method offers many advantages as an electromagnetic modeling, simulation and analysis tool. Its capabilities include: 1) broadband response predictions centered about the system resonances, 2) arbitrary 2-D or 3-D model geometries, 3) interaction with an object of any conductivity from that of a perfect conductor, to that of a real metal, to that of low or zero conductivity, 4) frequency-dependent constitutive parameters for modeling most materials (lossy dielectrics, magnetic materials, anisotropic plasmas and magnetized ferrites) and 5) any type of response, including far fields derived from near fields, such as scattered fields, antenna patterns, radar-cross sections, surface response, current and power density and penetration/coupling. FTDT is an appropriate technique for time-domain analysis of passive microwave and millimeter wave structures. The strength of FDTD method lies in its capability to model any volumetric structure with what are typically rectangular cells. No requirements for symmetry or smooth surfaces are needed. Also, FDTD has been successfully used in modeling micro-structured

Coupled mode theory (CMT), which uses BPM method, was applied to this structure by means of BeamPROP® tool (Synopsys™ Software). This method also takes into account the effective index of the ring resonator material, so the expected results from it will be improved compared with those obtained with the EIM method alone. BPM numerical simulation results can be seen in Figure 16. Finite element method (FEM) was applied to this structure by means of FemSIM® tool (Synopsys™ Software). This method also takes into account the effective index of the ring resonator material, so the expected results from it will be improved compared with those obtained with the EIM method alone. FEM

Result of a FDTD method numerical simulation applied to the add-drop filter configuration considered here is represented on Figure 18 (resonance takes place at 1285 nm wavelength). As it was said above, for this structure Eq. (25) gives the value of 5.475 for F finesse and Eq. (26) yields 321.40 for the Q quality-factor of the

As a consequence of transmission spectrum of add-drop filter configuration,

this can find applications such as optical band-pass filter, optical delay line,

CMT method numerical simulation result (from BeamPROP® tool) for all-pass FOPRRI configuration

numerical simulation of this case is represented in Figure 17.

micro-strip structures.

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

system.

Figure 16.

65

(resonance at 1285 nm).

#### Figure 15.

Transmission spectra of add-drop configuration: (a) pass-through port, (b) drop port and (c) drop-port (expanded).

known as the full wave method. It discretizes the equations and solves them using a linear algebraic algorithm.

To accelerate the finite difference method simulation, standard approaches such as Krylov subspace-based model order reduction or multi-grid analysis are used. It can be used for all types of structures. As compared with BPM, FDTD is at least an order of magnitude slower. The basic elements that are of application in the case of

## Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

FDTD method allow as broad as the Maxwell's equations they embody in discretized form, the basic formulation using centered finite difference expressions, which can treat virtually any material type and geometry and the basic resources needed to implement the technique. The FDTD method treats transients in the time domain (e.g., pulses), and it is applicable over the computationally unpredictable resonance region in which a wavelength is comparable with to the object interaction size. The FDTD method offers many advantages as an electromagnetic modeling, simulation and analysis tool. Its capabilities include: 1) broadband response predictions centered about the system resonances, 2) arbitrary 2-D or 3-D model geometries, 3) interaction with an object of any conductivity from that of a perfect conductor, to that of a real metal, to that of low or zero conductivity, 4) frequency-dependent constitutive parameters for modeling most materials (lossy dielectrics, magnetic materials, anisotropic plasmas and magnetized ferrites) and 5) any type of response, including far fields derived from near fields, such as scattered fields, antenna patterns, radar-cross sections, surface response, current and power density and penetration/coupling. FTDT is an appropriate technique for time-domain analysis of passive microwave and millimeter wave structures. The strength of FDTD method lies in its capability to model any volumetric structure with what are typically rectangular cells. No requirements for symmetry or smooth surfaces are needed. Also, FDTD has been successfully used in modeling micro-structured micro-strip structures.

Coupled mode theory (CMT), which uses BPM method, was applied to this structure by means of BeamPROP® tool (Synopsys™ Software). This method also takes into account the effective index of the ring resonator material, so the expected results from it will be improved compared with those obtained with the EIM method alone. BPM numerical simulation results can be seen in Figure 16.

Finite element method (FEM) was applied to this structure by means of FemSIM® tool (Synopsys™ Software). This method also takes into account the effective index of the ring resonator material, so the expected results from it will be improved compared with those obtained with the EIM method alone. FEM numerical simulation of this case is represented in Figure 17.

Result of a FDTD method numerical simulation applied to the add-drop filter configuration considered here is represented on Figure 18 (resonance takes place at 1285 nm wavelength). As it was said above, for this structure Eq. (25) gives the value of 5.475 for F finesse and Eq. (26) yields 321.40 for the Q quality-factor of the system.

As a consequence of transmission spectrum of add-drop filter configuration, this can find applications such as optical band-pass filter, optical delay line,

#### Figure 16.

CMT method numerical simulation result (from BeamPROP® tool) for all-pass FOPRRI configuration (resonance at 1285 nm).

known as the full wave method. It discretizes the equations and solves them using a

Transmission spectra of add-drop configuration: (a) pass-through port, (b) drop port and (c) drop-port

Interferometry - Recent Developments and Contemporary Applications

To accelerate the finite difference method simulation, standard approaches such as Krylov subspace-based model order reduction or multi-grid analysis are used. It can be used for all types of structures. As compared with BPM, FDTD is at least an order of magnitude slower. The basic elements that are of application in the case of

linear algebraic algorithm.

Figure 15.

(expanded).

64

Interferometry - Recent Developments and Contemporary Applications

taking into account three criteria: spent computational time (seconds), memory

A conventional model of a fiber optic passive ring resonator interferometer (FOPRRI) was theoretically analyzed in two main configurations: all-pass filter and add-drop filter. Then, four numerical simulation methods (EIM, CMT, FEM and 3D-FDTD) are successively applied on each of these two configurations. In both cases, results showed that EIM method is the most suitable when no very accuracy of results is needed (error < 7%) and low time of computation is desirable (<70 s). However, when high accuracy is considered the most important issue, 3D-FDTD is the best method (error < 1%) at the expense of a much larger computing time (≈1200 sec) than others. CMT (BPM) method provides substantial performance results with computational time similar than that EIM one, but an error < 1%. Finally, FEM is suitable for simulations that can be approximated with 2D techniques since it provides substantial performance results (error < 1%) with compu-

The author wishes to thank the Integrated Optics and Optoelectronics Laboratory of Faculty of Sciences of University of Oviedo-Spain for technical support and

resources applied (MByte) and error found (%).

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

tational time only an order of magnitude above EIM or CMT.

\*Address all correspondence to: ramonjose.perez@lugo.uned.es

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

assistance provided in the realization of this work.

7. Conclusions

Acknowledgements

Author details

UNED, Lugo, Spain

67

Ramón José Pérez Menéndez

provided the original work is properly cited.

Figure 17. FEM numerical simulation (from FemSIM® tool) for all-pass FOPRRI configuration (resonance at 1285 nm).

Figure 18. 2D-FTDT numerical simulation results for an add-drop FOPRRI (resonance at 1285 nm).

optical switch, inertial rotation sensor, photonic biosensor, tailored response filter and active ring modulator.

## 6. Performance comparison of numerical simulation methods used

Table 2 shows performance comparison of all four methods used (EIM, CMT, FEM and FDTD) for both the configurations: all-pass filter and add-drop filter and


Table 2.

Comparison of computational time, memory resources and error methods (EIM, CMT, FEM and FDTD).

taking into account three criteria: spent computational time (seconds), memory resources applied (MByte) and error found (%).

## 7. Conclusions

A conventional model of a fiber optic passive ring resonator interferometer (FOPRRI) was theoretically analyzed in two main configurations: all-pass filter and add-drop filter. Then, four numerical simulation methods (EIM, CMT, FEM and 3D-FDTD) are successively applied on each of these two configurations. In both cases, results showed that EIM method is the most suitable when no very accuracy of results is needed (error < 7%) and low time of computation is desirable (<70 s). However, when high accuracy is considered the most important issue, 3D-FDTD is the best method (error < 1%) at the expense of a much larger computing time (≈1200 sec) than others. CMT (BPM) method provides substantial performance results with computational time similar than that EIM one, but an error < 1%. Finally, FEM is suitable for simulations that can be approximated with 2D techniques since it provides substantial performance results (error < 1%) with computational time only an order of magnitude above EIM or CMT.

## Acknowledgements

The author wishes to thank the Integrated Optics and Optoelectronics Laboratory of Faculty of Sciences of University of Oviedo-Spain for technical support and assistance provided in the realization of this work.

## Author details

optical switch, inertial rotation sensor, photonic biosensor, tailored response

2D-FTDT numerical simulation results for an add-drop FOPRRI (resonance at 1285 nm).

FEM numerical simulation (from FemSIM® tool) for all-pass FOPRRI configuration (resonance at 1285 nm).

Interferometry - Recent Developments and Contemporary Applications

6. Performance comparison of numerical simulation methods used

Table 2 shows performance comparison of all four methods used (EIM, CMT, FEM and FDTD) for both the configurations: all-pass filter and add-drop filter and

Criterion/Method EIM CMT (BPM) FEM 2D-FDTD/3D-FDTD

Comparison of computational time, memory resources and error methods (EIM, CMT, FEM and FDTD).

Computational time (sec) 2D < 0.05 <0.10 1.7 4.3 Computational time (sec) 3D 68 100 547 1170 Memory (MByte) 16 17 17 17 Error (%) <7 <1 <1 <1

filter and active ring modulator.

Figure 17.

Figure 18.

Table 2.

66

Ramón José Pérez Menéndez UNED, Lugo, Spain

\*Address all correspondence to: ramonjose.perez@lugo.uned.es

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## References

[1] Hariharan P. Optical Interferometry. 2nd ed. London: Elsevier, Academic Press; 2003

[2] Sharma U, Wei X. Fiber optic interferometric devices. Chapter 2. In: Kang JU, editor. Fiber Optic Sensing and Imaging. New York: Springer; 2013. pp. 29-53

[3] Yin S, Ruffin PB, Yu FTS, editors. Fiber Optic Sensors. 2nd ed. Boca Raton FL (USA): CRC Press; 2008. pp. 333-366

[4] Yariv A, Yeh P. Optical Resonators from Photonics. Chapter 4. 6th ed. Oxford University Press, New York; 2007. pp. 156-210

[5] Vali V, Shorthill RW. Fiber ring interferometer. Applied Optics. 1976; 15(5):1099-2000

[6] Vali V, Shorthill RW. Fiber laser gyroscopes. Proceedings of SPIE. 1976; 76:110

[7] Vali V, Shorthill RW. Ring interferometer 950 m long. Applied Optics. 1977;16:290

[8] Ezequiel S, Balsamo SR. Passive ring resonator laser gyroscope. Applied Physics Letters. 1977;30:478

[9] Meyer RE, Ezequiel S. Passive fiberoptic ring resonator for rotation sensing. Optics Letters. 1984;8:644-646

[10] Rao YJ, Jackson DA. Principles of fiber-optic interferometry. In: Grattan LS, Meggitt BT, editors. Optical Fiber Sensor Technology. Boston, MA: Springer; 2000. pp. 167-191

[11] Yve CY, Peng JD, Liao YB, Zhou BK. Fibre ring resonator with finesse of 1260. Electronics Letters. 1988;24:622-623

[12] Stokes LF, Chodorow M, Shaw HJ. All-single-mode fiber resonator. Optics Letters. 1982;7:288-290

characteristics of fiber ring resonators. IEEE Journal of Quantum Electronics.

Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

> [34] RSoft Products-Synopsys Optical Solutions. URL: http://optics.synopsys.

> [36] RSoft Products-Synopsys Optical Solutions. URL: http://optics.synopsys.

[38] Thylén L. The beam propagation method: An analysis of its applicability. Optics Quantum Electron. 1983;15:

Nakano H. Beam propagation method using Padé approximant operators. Transactions IEICE Japan. 1994;J-77-C-I:

[40] Hadley GR. Wide-angle beam propagation using Padé approximant operators. Optics Letters. 1992;17:

[41] Liu PL, Li BJ. Study on form birefringence in waveguide devices using the semi-vectorial beam propagation method. IEEE

Photonics Technology Letters. 1991;3:

[42] Koshiba M, Tsuji Y. A wide-angle finite element beam propagation method. IEEE Photonics Technology

Letters. 1996;8:1208-1210

McGraw-Hill; 1973

[43] Zienckiewitz OC. The Finite Element Method. 3rd ed. New York:

[37] Feit MD, Freck JA Jr. Light propagation in graded-index optical fibers. Applied Optics. 1978;17:

[39] Yamauchi J, Shibayama J,

[35] FDTD Solutions-Lumerical's Nanophotonic FDTD Simulation Software. URL: http://lumerical.com/

com/RSoft/FemSIM

com/RSoft/FullWAVE

fdtd

3990-3998

433-439

490-494

1426-1428

913-915

[24] Bogaerts W et al. Silicon micro-ring

[25] Knox RM, Thoulios PP. Integrated circuits for the millimeter through optical frequency range. In: Proceedings of Symposium on Submillimeter Waves. Vol. 20; Brooklyn-New York; 1970.

[26] Marcatili EA. Dielectric rectangular waveguide and directional coupler for integrated optics. The Bell System Technical Journal. 1969;48:2071-2102

[27] Mode Solutions-Waveguide Mode Solver and Propagation Simulator. URL: http://www.lumerical.com/mode

[28] WGMODES-Photonics Research Laboratory. URL: http://www.

photonics.umd.edu/software/wgmodes/

[29] COMSOL Multiphysics. URL: http://

[30] Phoenix Software. URL: http:// www.phoenixbv.com/FieldDesigner

[31] Yariv A. Coupled-mode theory for guided-wave optics. IEEE Journal of Quantum Electronics. September 1973;

[32] RSoft Products-Synopsys Optical Solutions. URL: http://optics.synopsys.

[33] Phoenix Software. URL: http:// www.phoenixbv.com/OptoDesigner

comsol.com/mode

QE-9(9):919-9333

69

com/RSoft/BeamPROP

[23] Vien V. Optical Micro-Ring Resonators: Theory, Techniques and Applications. Boca Raton, FL: CRC

resonators. Laser and Photonics Reviews. 2012;6(1):47-53

2004;40:736-730

Press; 2017

pp. 497-515

[13] Stokes LF, Chodorow M, Shaw HJ. All-fiber stimulated Brillouin ring laser with sub-milliwatt pump threshold. Optics Letters. 1982;7:509-511

[14] Jackson DA, Jones JDC. Fiber-optic sensors. Optica Acta. 1986;33:1469-1503

[15] Zarinetchi F, Smith SP, Ezekiel S. Stimulated Brillouin fiber-optic laser gyroscope. Optics Letters. 1991;16: 229-231

[16] Hill KO, Kawasaki BS, Johnson DC. CW Brillouin laser. Applied Physics Letters. 1976;28:608-609

[17] Kalli K, Jackson DA. Ring resonator optical spectrum analyzer with 20 kHz resolution. Optics Letters. 1992;17: 1090-1092

[18] Tai S, Kyuma K, Nakayama T. Novel measuring method for spectral linewidth of laser diodes using fiberoptic ring resonators. Electronics Letters. 1985;21:91-93

[19] Newton SA, Howland RS, Jackson KP, Shaw HJ. High speed pulse-train generation using single-mode fiber recirculating delay lines. Electronics Letters. 1983;19:757-758

[20] Rabus D. Integrated Ring Resonators. New York: Springer; 2007

[21] Heebner JE et al. Optical Micro-Resonators: Theory, Fabrication and Applications. London: Springer; 2008

[22] Heebner JE, Wong V, Schweingsberg A, Boyd RW, Jackson DJ. Optical transmission Fiber-Optic Ring Resonator Interferometer DOI: http://dx.doi.org/10.5772/intechopen.80569

characteristics of fiber ring resonators. IEEE Journal of Quantum Electronics. 2004;40:736-730

References

Press; 2003

pp. 29-53

2007. pp. 156-210

15(5):1099-2000

Optics. 1977;16:290

76:110

68

[1] Hariharan P. Optical Interferometry. 2nd ed. London: Elsevier, Academic

Interferometry - Recent Developments and Contemporary Applications

[12] Stokes LF, Chodorow M, Shaw HJ. All-single-mode fiber resonator. Optics

[13] Stokes LF, Chodorow M, Shaw HJ. All-fiber stimulated Brillouin ring laser with sub-milliwatt pump threshold. Optics Letters. 1982;7:509-511

[14] Jackson DA, Jones JDC. Fiber-optic sensors. Optica Acta. 1986;33:1469-1503

[15] Zarinetchi F, Smith SP, Ezekiel S. Stimulated Brillouin fiber-optic laser gyroscope. Optics Letters. 1991;16:

[16] Hill KO, Kawasaki BS, Johnson DC. CW Brillouin laser. Applied Physics

[17] Kalli K, Jackson DA. Ring resonator optical spectrum analyzer with 20 kHz resolution. Optics Letters. 1992;17:

[18] Tai S, Kyuma K, Nakayama T. Novel

[19] Newton SA, Howland RS, Jackson KP, Shaw HJ. High speed pulse-train generation using single-mode fiber recirculating delay lines. Electronics

Resonators. New York: Springer; 2007

[21] Heebner JE et al. Optical Micro-Resonators: Theory, Fabrication and Applications. London: Springer; 2008

measuring method for spectral linewidth of laser diodes using fiberoptic ring resonators. Electronics

Letters. 1985;21:91-93

Letters. 1983;19:757-758

[20] Rabus D. Integrated Ring

[22] Heebner JE, Wong V, Schweingsberg A, Boyd RW, Jackson DJ. Optical transmission

Letters. 1976;28:608-609

Letters. 1982;7:288-290

229-231

1090-1092

[2] Sharma U, Wei X. Fiber optic interferometric devices. Chapter 2. In: Kang JU, editor. Fiber Optic Sensing and Imaging. New York: Springer; 2013.

[3] Yin S, Ruffin PB, Yu FTS, editors. Fiber Optic Sensors. 2nd ed. Boca Raton FL (USA): CRC Press; 2008. pp. 333-366

[4] Yariv A, Yeh P. Optical Resonators from Photonics. Chapter 4. 6th ed. Oxford University Press, New York;

[5] Vali V, Shorthill RW. Fiber ring interferometer. Applied Optics. 1976;

[6] Vali V, Shorthill RW. Fiber laser gyroscopes. Proceedings of SPIE. 1976;

[8] Ezequiel S, Balsamo SR. Passive ring resonator laser gyroscope. Applied Physics Letters. 1977;30:478

[9] Meyer RE, Ezequiel S. Passive fiberoptic ring resonator for rotation sensing.

[10] Rao YJ, Jackson DA. Principles of fiber-optic interferometry. In: Grattan LS, Meggitt BT, editors. Optical Fiber Sensor Technology. Boston, MA: Springer; 2000. pp. 167-191

[11] Yve CY, Peng JD, Liao YB, Zhou BK. Fibre ring resonator with finesse of 1260. Electronics Letters. 1988;24:622-623

Optics Letters. 1984;8:644-646

[7] Vali V, Shorthill RW. Ring interferometer 950 m long. Applied [23] Vien V. Optical Micro-Ring Resonators: Theory, Techniques and Applications. Boca Raton, FL: CRC Press; 2017

[24] Bogaerts W et al. Silicon micro-ring resonators. Laser and Photonics Reviews. 2012;6(1):47-53

[25] Knox RM, Thoulios PP. Integrated circuits for the millimeter through optical frequency range. In: Proceedings of Symposium on Submillimeter Waves. Vol. 20; Brooklyn-New York; 1970. pp. 497-515

[26] Marcatili EA. Dielectric rectangular waveguide and directional coupler for integrated optics. The Bell System Technical Journal. 1969;48:2071-2102

[27] Mode Solutions-Waveguide Mode Solver and Propagation Simulator. URL: http://www.lumerical.com/mode

[28] WGMODES-Photonics Research Laboratory. URL: http://www. photonics.umd.edu/software/wgmodes/

[29] COMSOL Multiphysics. URL: http:// comsol.com/mode

[30] Phoenix Software. URL: http:// www.phoenixbv.com/FieldDesigner

[31] Yariv A. Coupled-mode theory for guided-wave optics. IEEE Journal of Quantum Electronics. September 1973; QE-9(9):919-9333

[32] RSoft Products-Synopsys Optical Solutions. URL: http://optics.synopsys. com/RSoft/BeamPROP

[33] Phoenix Software. URL: http:// www.phoenixbv.com/OptoDesigner [34] RSoft Products-Synopsys Optical Solutions. URL: http://optics.synopsys. com/RSoft/FemSIM

[35] FDTD Solutions-Lumerical's Nanophotonic FDTD Simulation Software. URL: http://lumerical.com/ fdtd

[36] RSoft Products-Synopsys Optical Solutions. URL: http://optics.synopsys. com/RSoft/FullWAVE

[37] Feit MD, Freck JA Jr. Light propagation in graded-index optical fibers. Applied Optics. 1978;17: 3990-3998

[38] Thylén L. The beam propagation method: An analysis of its applicability. Optics Quantum Electron. 1983;15: 433-439

[39] Yamauchi J, Shibayama J, Nakano H. Beam propagation method using Padé approximant operators. Transactions IEICE Japan. 1994;J-77-C-I: 490-494

[40] Hadley GR. Wide-angle beam propagation using Padé approximant operators. Optics Letters. 1992;17: 1426-1428

[41] Liu PL, Li BJ. Study on form birefringence in waveguide devices using the semi-vectorial beam propagation method. IEEE Photonics Technology Letters. 1991;3: 913-915

[42] Koshiba M, Tsuji Y. A wide-angle finite element beam propagation method. IEEE Photonics Technology Letters. 1996;8:1208-1210

[43] Zienckiewitz OC. The Finite Element Method. 3rd ed. New York: McGraw-Hill; 1973

[44] Koshiba M, Saitoh H, Eguchi M, Hirayama K. Simple scalar finite element approach to optical waveguides. IEEE Proceedings Journal. 1992;139:166-171

[45] Koshiba M. Optical Waveguide Theory by the Finite Element Method. Dordecht, Holland: KTK Scientific Publishers and Kluwer Academic Publishers; 1992

[46] Meunier G. The Finite Element Method for Electromagnetic Modeling. River Street-Hoboken NJ: John Wiley & Sons; 2008

[47] Kwon YW, Bang H. The Finite Element Method Using Matlab. CRC Press; 1977

[48] Stern M. Semivectorial polarized finite difference method optical waveguides with arbitrary index profiles. IEEE Proceedings Journal. Piscataway NJ (USA); 1988;135:56-63

[49] Stern M. Semivectorial polarized Hfield solutions for dielectric waveguides with arbitrary index profiles. IEEE Proceedings Journal. Piscataway NJ (USA); 1988;135:333-338

[50] Kawano K, Kitoh T, Kohtuku M, Takeshita T, Hasumi Y. 3-D semivectorial analysis to calculate facet reflectivities of semiconductor optical waveguides based on the bi-directional method of line BPM (MoL-BPM). IEEE Photonics Technology Letters. 1998;10: 108-110

[51] Gerdes B, Lunitz B, Benish D, Pregla R. Analysis of slab waveguide discontinuities including radiation and absorption effects. Electronics Letters. 1992;28:1013-1015

[52] Taflove A. Computational Electrodynamics: Finite Difference Time Domain Method. Artech House: Norwood, MA; 1995

**71**

**Chapter 5**

**Abstract**

Detection

*and Pavel V. Minkovich*

interferometry, photonic crystal fiber

**1. Introduction**

Modal Interferometers Based

on a Tapered Special Photonic

*Vladimir (or Uladzimir) P. Minkovich, Joel Villatoro*

visualized as a shift of the peaks in the output spectrum pattern.

**Keywords:** modal interferometer, optical fiber devices, optical fiber sensors,

Photonic crystal fibers (PCFs), also known as microstructured or holey optical fibers [1–4], consist of a waveguiding core surrounded by a system of air channels in glass cladding that run along the length of the fiber and arranged in a hexagonal structure around the core. The structure of PCFs enables new possibilities for optical sensing in comparison with standard optical fibers. The most common approach consists of making a sample to interact with the evanescent fields of the PCF guided modes [5–7]. To do so one has to fill the holes with the sample, a gas or liquid, for example, and then the analysis or detection is carried out. In some situation such a process may be inconvenient or impractical. Mach-Zehnder modal interferometers (MZMIs) based on no adiabatic tapered silica PCFs, first fabricated in our facilities [8–10], are attractive for sensing application because of their intrinsic advantages, such as simplicity of fabrication and practical using, high sensitivity, small size, and immunity to electromagnetic interference and their aptitude for remote measurements. Initially, our MZMIs were used for high-resolution refractive index sensing of liquids with indices ranging from 1.41 to 1.45 [8]. In this work, it is presented

Crystal Fiber for Highly Sensitive

The use of a tapered special photonic crystal fiber (PCF) with collapsed air holes in the waist (the thinnest part of a taper) for highly sensitive detection of strain, high temperature, and fast detection of hydrogen with concentrations between 1.2 and 5.6 vol.% and biosensing is demonstrated. In the tapered PCF, a fundamental core mode couples to a few modes of the solid taper waist. Owing to the beating between the waist modes, the transmission spectra of the tapered PCF exhibit several interference peaks, which are sensitive to refractive index changes of a medium that surrounds the taper and also to changes of a taper length. The changes can be

## **Chapter 5**

[44] Koshiba M, Saitoh H, Eguchi M, Hirayama K. Simple scalar finite element approach to optical

Interferometry - Recent Developments and Contemporary Applications

waveguides. IEEE Proceedings Journal.

[45] Koshiba M. Optical Waveguide Theory by the Finite Element Method. Dordecht, Holland: KTK Scientific Publishers and Kluwer Academic

[46] Meunier G. The Finite Element Method for Electromagnetic Modeling. River Street-Hoboken NJ: John Wiley &

[47] Kwon YW, Bang H. The Finite Element Method Using Matlab. CRC

[48] Stern M. Semivectorial polarized finite difference method optical waveguides with arbitrary index profiles. IEEE Proceedings Journal. Piscataway NJ (USA); 1988;135:56-63

[49] Stern M. Semivectorial polarized Hfield solutions for dielectric waveguides with arbitrary index profiles. IEEE Proceedings Journal. Piscataway NJ

[50] Kawano K, Kitoh T, Kohtuku M,

semivectorial analysis to calculate facet reflectivities of semiconductor optical waveguides based on the bi-directional method of line BPM (MoL-BPM). IEEE Photonics Technology Letters. 1998;10:

[51] Gerdes B, Lunitz B, Benish D, Pregla

discontinuities including radiation and absorption effects. Electronics Letters.

(USA); 1988;135:333-338

Takeshita T, Hasumi Y. 3-D

R. Analysis of slab waveguide

[52] Taflove A. Computational Electrodynamics: Finite Difference Time Domain Method. Artech House:

1992;28:1013-1015

Norwood, MA; 1995

1992;139:166-171

Publishers; 1992

Sons; 2008

Press; 1977

108-110

70
