Comparison of numerical simulation methods (EIM, CMT, FEM and FDTD) for all-pass/add-drop filters.
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.
- 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)
Optical interferometry constitutes an important technique used in high number of measurement processes for multiple physical magnitudes and quantitative phenomenon . 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 . 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 fiber-optic ring interferometer devices was in 1976 when Vali and Shorthill designed the first fiber-optic gyroscope model [5, 6, 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 .
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 , being the phase shift due to cross-coupling and
where γ is the insertion loss of the coupler, is the optical resonance coupling coefficient,
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
where sign and
The feedback path of length 2
Since that adding or subtracting an integer number
where , representing the coupling losses and represents the phase of the coupler. Making a similar deduction, we can reach the circulating power in the ring, which is:
where, as we have said above,
Then, the transmitted power to through port (3) becomes 0. The critical coupling is due to destructive interference between the direct beam and the cross-coupled beam of light, both after passing their respective coupling zone. A condition termed
A good idea of the behavior of the system, when resonance condition accomplishes, can be obtained when one represents the T
Figure 3(a) represents the T transmission function for an all-pass ring resonator with
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
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% (), 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 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,
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:
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 ). 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,
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 :
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:
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,
In this equation, and are used, which are the half round-trip loss and phase, respectively. It is and . The mode amplitude in the ring has to pass the second coupler to become the new dropped mode rms amplitude
At resonance condition (,
The through port mode amplitude
According with Eq. (18), the transfer function of through port (relative to input port) is given by:
Taking into account, Eq. (19), the transfer function of drop port (relative to input port), will be:
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 (),
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].
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, 28, 29, 30, 31, 32, 33] will be performed to the same structure. In third time, finite element method (FEM)  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 h-parameter = 2 10−6 m−1. 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 steady-state solution. Obtained transmission power spectrum is represented in Figure 9.
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) . 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. Further in-depth readings on BPM method can be found in Refs. [37, 38, 39, 40, 41, 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) . 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, 44, 45, 46, 47]. The result of FemSIM® of all-pass ring resonator analyzed here is represented in Figure 12.
FullWAVE® (Synopsys™ Software)  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 time-space 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, 49, 50, 51, 52]. Figure 13 represents the result of FDTD numerical simulation of all-pass ring resonator studied here.
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.
|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 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 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 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 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, 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 taking into account three criteria: spent computational time (seconds), memory resources applied (MByte) and error found (%).
|Computational time (sec) 2D||< 0.05||<0.10||1.7||4.3|
|Computational time (sec) 3D||68||100||547||1170|
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.
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.