1. Introduction
Integrated photonic sensors hold significant importance in various fields due to their ability to accurately and efficiently detect and measure physical, chemical, and biological parameters[
1,
2]. These sensors utilize light to sense changes in their surroundings, offering advantages such as high sensitivity, fast response times, and miniaturization[
3]. Their integration into compact and scalable platforms enables the development of portable, low-cost, and highly sensitive sensing devices for applications ranging from environmental monitoring and healthcare diagnostics to telecommunications and industrial process control[
4,
5,
6]. By leveraging the unique properties of light, integrated photonic sensors contribute to advancements in fields such as biomedical research, environmental sustainability, and the optimization of industrial processes, ultimately fostering innovation and improving quality of life[
7,
8,
9].
The Mach-Zehnder interferometer (MZI) structure holds immense significance in sensing applications due to its ability to efficiently manipulate and detect light signals[
10,
11]. This structure utilizes the principle of interference to measure changes in optical properties, making it highly sensitive to external stimuli such as temperature, pressure, and refractive index variations[
12,
13,
14]. Its compact size, low power consumption, and high sensitivity make it ideal for a wide range of sensing applications, including environmental monitoring, biomedical diagnostics, and industrial process control[
16,
17,
18]. Additionally, its versatility allows for integration with various materials and technologies, paving the way for innovative sensing solutions with improved performance and reliability.
The Sagnac loop serves multiple purposes in optical systems. In the context of MZI, like the loop-terminated Mach–Zehnder interferometer (LT-MZI), the Sagnac loop acts as a crucial component to enhance performance and functionality. One of the key functions of the Sagnac loop in an LT-MZI is to reduce the length of the interferometer by a factor of two. This reduction in length is advantageous for compactness and integration into miniaturized optical devices, allowing for space-efficient designs without compromising performance. Moreover, the Sagnac loop enhances the interferometer's phase sensitivity. By circulating light through the loop, it increases the effective path length travelled by the light, thereby amplifying the interference effects[
19]. In interferometric sensors, including those used for refractive index sensing, common mode rejection is crucial for eliminating unwanted noise and environmental perturbations. The Sagnac loop configuration inherently provides excellent common mode rejection, as any disturbance affecting both arms of the interferometer equally will not affect the interference pattern[
20].
Carbon dioxide (CO
2) monitoring is of paramount importance due to its significant implications for both environmental and human health. As a major greenhouse gas, CO
2 plays a crucial role in regulating the earth's climate, influencing temperature patterns, weather phenomena, and ultimately the balance of ecosystems[
21]. Moreover, in indoor environments, monitoring CO
2 is essential for ensuring air quality and occupant comfort. Elevated CO
2 levels can indicate poor ventilation, which can lead to a buildup of other pollutants and potential health risks, including headaches, fatigue, and respiratory problems[
22,
23]. By monitoring CO
2 concentrations, it is possible to maintain optimal indoor air quality, promote productivity, and safeguard the well-being of occupants in various settings, including homes, schools, offices, and healthcare facilities. In this work, a novel design of a CO
2 gas sensor is proposed which is based on an LT-MZI embedded with polyhexamethylene biguanide (PHMB) functional polymer in the sensing arm to determine the minute change in the gas concentration.
2. Device Design and Analysis
In this paper, a novel approach for developing a highly compact CO
2 gas sensor is proposed by leveraging an SOI-based LT-MZI and PHMB polymer which is deposited onto the sensing arm as shown in
Figure 1. The inset shows the E-field distribution at 1550 nm in the cross-sectional view of the waveguide covered with 300 nm thick PHMB polymer. The height and width of the waveguide are considered at 220 nm and 250 nm, respectively. The LT-MZI configuration incorporates a Sagnac loop at the end of the conventional MZI, effectively halving the interferometer's length. This advancement holds promise for enhanced sensor performance, scalability, and versatility in various applications requiring precise refractive index measurements. The length of the directional coupler waveguide and arms of the MZI segment are designated by L
DC and L, respectively. The waveguide maintains a consistent width, denoted as w, across its entirety. Introducing the gap, represented as g, between the two waveguides of the directional coupler segment facilitates efficient coupling. Moreover, the fixed bending radius, denoted as R, set at 2 µm, ensures the design's compactness, a crucial aspect in crafting effective sensing devices.
Table 1 presents the range of the geometric parameters used in the device design.
In an accurate endeavor, rigorous simulations are conducted to dissect transmittance patterns and field distributions. Employing the 2D-finite element method (2D-FEM) via COMSOL Multi-physics software, the intricacies of light propagation are determined. Leveraging the rich materials library within COMSOL, refractive indices for pivotal materials like silicon, SiO2 and air are obtained, ensuring alignment with real-world parameters. To maintain precision, the device design undergoes precise segmentation into subdomains, further subdivided into triangular mesh elements. Each element features a grid size precisely set at λ/100, striking an optimal balance between simulation accuracy and computational efficiency. In the pursuit of realistic electromagnetic wave phenomena analysis, an open-bounded domain is established as a fundamental imperative. To emulate this environment faithfully, scattering boundary conditions (SBC) are judiciously implemented along the outer perimeters of the FEM simulation window.
Through these systematic approaches and the utilization of computing resources— an 8-core processor (3.8 GHz) paired with 128 GB of RAM—a suite of simulations characterized by their high fidelity and reliability are generated. These simulations offer invaluable insights into transmittance characteristics and field distributions, serving as a cornerstone for the optimization and advancement of photonic device designs.
3. Design of a 3 dB Directional Coupler
An integrated photonic 3dB directional coupler is a crucial component in optical communication systems, facilitating efficient power splitting and combining between two waveguides on a chip[
24]. It achieves a balanced power division (50:50) through precise design and fabrication techniques, where half of the input power is transferred to each output port. This device operates based on the principle of evanescent field coupling, enabling seamless transmission of optical signals with minimal loss and crosstalk. In general, the coupling length (L
DC) can be calculated using the formula:
Where βeff is the effective propagation constant of the coupled modes, the effective propagation constant βeff depends on the waveguide parameters, such as the waveguide’s refractive indices, the waveguide’s width and height, and the wavelength of light.
The transmission spectrum is plotted for a directional coupler based on an SOI platform with L
DC varying between 0.5 µm and 12 µm. The transmission is calculated as, P
out/P
in, where P
out and P
in are the output and input power, respectively, While designing with COMSOL Multiphysics tool, it considers several types of losses of the light propagation in optical waveguides. These include absorption losses, where the waveguide material absorbs light energy. Software also accounts for radiation losses from waveguide bends and mode mismatches, as well as coupling losses at interfaces and junctions. Additionally, mode conversion losses occur when higher-order modes are present or when unintended mode coupling happens. By incorporating these factors, accurate modeling of waveguide performance and efficiency is obtained. Incremental changes in the width (w) of the waveguide core from 250 nm to 280 nm are made to examine their influence on the coupling length, as depicted in
Figure 2 (a-d). Notably, the beat length of the directional coupler demonstrates an increase as w rises from 250 nm to 280 nm. In optical directional couplers, the beat length signifies the length over which the phase difference between the coupled modes changes by 2π radians (one full cycle). This metric is influenced by factors including the wavelength of light, refractive indices of the materials, and the physical dimensions of the coupler. Consequently, for w=250 nm, achieving a 3 dB power distribution remains feasible even at an L
DC as small as 3 µm.
Figure 2 (e) and
Figure 2 (f) show the normalized H-field distribution in the directional coupler for L
DC=3 µm and L
DC= 4 µm, respectively.
4. Optimization of LT-MZI Structure
In this section, the LT-MZI structure undergoes optimization to ascertain the spectral characteristics of the device. At this stage, the LT-MZI structure is examined in its bare form, devoid of any PHMB polymer covering on the sensing arm and lacks an upper cladding in this particular configuration. Initially, the impact of variations in LDC on both interference dips and the free spectral range (FSR) is thoroughly investigated which is a fundamental parameter that defines the spacing between consecutive resonance frequencies within the device. It represents the spectral range over which the resonator can operate without encountering spectral overlap between adjacent resonances. Essentially, the FSR determines the resolution or spacing capability of the resonator's spectral response. It is typically influenced by various factors such as the physical dimensions of the resonator, including its length, width, and height, as well as the refractive index of the material.
Interference dips, a key aspect in MZI structures, signify instances where the output intensity of the interferometer diminishes due to interference between the two optical paths. LDC is varied between 3 µm and 6 µm whereas the other geometric parameters such as w, g and L are fixed at 250 nm, 100 nm and 5 µm, respectively. The transmission of the device is plotted by utilizing the expression: ; where Pout and Pin are the output and input power, respectively.
From
Figure 3 (a-g), it's evident that as L
DC increases from 3 µm to 6 µm, the FSR of the device decreases from 90 nm to 57 nm, accompanied by a noticeable blueshift in the dips. The extinction ratio (ER) of the device ranges from 19 dB to 27 dB, indicating a robust destructive interference resulting from the interaction of two light waves with opposite phases. For further investigation, L
DC is fixed at 6 µm which provides a good transmission characteristic in the desired wavelength range.
The width (w) of the waveguide in integrated photonic sensors plays an important role in defining the sensitivity of the device. As w approaches cut-off, even slight variations in factors like material properties, dimensions, or the presence of external elements can cause noticeable alterations in its transmission characteristics. This heightened sensitivity is advantageous for sensing applications because it enables the detection of subtle changes in the environment with greater precision. Consequently, designing a waveguide near cut-off can enhance the sensitivity of sensors based on waveguide principles, making them more effective for detecting changes in parameters such as temperature, pressure, or chemical composition.
The impact of w on the FSR of the device, spanning from 250 nm to 270 nm is determined, as illustrated in
Figure 4 (a-c). While keeping other geometric parameters such as L
DC, g, and L constant at 6 µm, 100 nm, and 5 µm respectively, a notable trend is observed. Specifically, the FSR exhibited a discernible rise from 57 nm to 92.5 nm with an increment in w from 250 nm to 270 nm. This observation underscores the pivotal role that waveguide dimensions, which in turn contribute to the total path length of the device, play in shaping its spectral characteristics.
In the design of an LT-MZI, the gap between the two waveguides of a 3 dB directional coupler significantly impacts the power division across the wavelength spectrum. The spacing between the waveguides determines the coupling strength between them, which in turn affects the transfer of power from one waveguide to the other. In the context of a wavelength spectrum, variations in the g can lead to changes in the interference pattern observed within the coupler. When g is small, i.e. 100 nm, the FSR of the LT-MZI is 57 nm, which increases beyond 80 nm as g approaches 175 nm (see
Figure 4 (d-g)). Therefore, precise control over the waveguide separation is crucial for optimizing the power division and maintaining spectral integrity in the LT-MZI structure.
5. Results and Discussion
CO
2, a colorless, odorless, and non-flammable greenhouse gas, is primarily generated through the combustion of carbonaceous materials or animal metabolism. Its escalating concentration and growth rate pose severe threats to the environment, including global warming and unsustainable environmental challenges[
25]. Consequently, there is a pressing need for highly sensitive CO
2 gas sensors to monitor greenhouse gas concentrations effectively. With increasing CO
2 levels in the surrounding medium, the refractive index of functionalization layers undergoes discernible changes as shown in
Table 2.
After fabricating the primary sensing structure using electron-beam lithography (EBL) and reactive ion etching (RIE), a second lithography process is required. In this process, a photoresist layer is spin-coated over the entire sample. Subsequently, a small sensing window is left open on the sensing arm of the device. PHMB is then spin-coated within this window, and the thickness of the PHMB layer is precisely controlled by adjusting the spin-coating speed. For effective CO
2 gas sensing applications, a PHMB thickness of approximately 240-250 nm is optimal[
26].
Polyhexamethylene biguanide (PHMB), featuring basic amide-bearing functional groups and belonging to the guanidine polymer family, emerges as a prime candidate for functionalization polymer layers due to its sensitivity and selectivity to low CO
2 concentrations at room temperature and atmospheric pressure[
27]. Notably, unlike previous functional materials, PHMB does not necessitate the presence of water vapor for CO
2 absorption and release. Gas sensors employing PHMB exhibit a linear correlation between the resonance wavelength blueshift and CO
2 concentration[
22,
29]. Remarkably, there are no reported instances of gases other than CO
2 interfering with the sensitivity of the PHMB polymer layer, rendering it highly suitable for CO
2 gas sensing applications. In a study analyzing PHMB's selectivity as a functional layer for CO
2 gas against hydrogen gas, it was confirmed that the sensor displayed no response to hydrogen or nitrogen gas, ensuring accurate detection of CO
2 without false positives[
26].
The concentration of CO2 gas deemed dangerous for humans typically begins to pose health risks when it exceeds certain thresholds. While CO2 is naturally present in the atmosphere at concentrations around 400 ppm, indoor environments can experience elevated levels due to poor ventilation or specific activities such as combustion. Generally, concentrations above 1000 ppm are considered to cause discomfort and may lead to symptoms such as headaches, dizziness, and drowsiness. At higher concentrations, exceeding 5000 ppm, CO2 can induce more severe symptoms, including difficulty breathing, increased heart rate, and impaired cognitive function. Extremely high concentrations, reaching tens of thousands of ppm, are life-threatening and can result in unconsciousness or even death. Therefore, maintaining CO2 levels below 1000 ppm in indoor environments is advisable to ensure the safety and well-being of occupants.
Initially, the MZI structure without an end-loop is investigated for its CO
2 gas sensing capabilities. The schematic diagram of this structure is depicted in
Figure 5(a). The intensity interrogation method in a MZI structure involves monitoring changes in the intensity of light at the output_1 port to detect variations in external conditions, such as gas concentration. In this method, the interference pattern of light, which results from the phase difference between the two arms of the interferometer, is analyzed. As the external condition (e.g., gas concentration) changes, it affects the refractive index of the PHMB polymer resulting in the shift in the effective refractive index in the sensing arm of the MZI, leading to a shift in the interference pattern. This shift causes a change in the intensity of light at specific wavelengths.
The transmission spectrum of the device is plotted over a wavelength range of 1450 nm to 1580 nm, with varying gas concentrations, as shown in
Figure 5(b). The geometric parameters of the device are fixed, with w = 250 nm, g = 100 nm, L
DC = 6 µm, and L = 5 µm. At the operational wavelength of 1535.5 nm, destructive interference occurs, resulting in a dip in the transmission spectrum that becomes more pronounced as the gas concentration increases.
Figure 5(c) illustrates the power measured at Output_1 at the wavelength of 1535.5 nm under different gas concentrations. It is observed that the power decreases from -7.15 dB to -15.38 dB as the gas concentration rises from 0 ppm to 366 ppm. However, beyond 366 ppm, a slight increase in power is noted, which may be attributed to weak destructive interference. Additionally, the H-field distribution within the device is presented for gas concentrations of 0 ppm and 366 ppm in
Figure 5(d) and
Figure 5(e), respectively. These figures provide a visual representation of the field changes due to varying gas concentrations.
The transmission spectrum of the LT-MZI structure under varying CO
2 gas concentrations is methodically depicted for two distinct arm lengths: L=5 µm and L=10 µm, illustrated in
Figure 6 (a) and
Figure 6 (b) respectively. Analyzing these figures reveals a noticeable trend: as the CO
2 gas concentration escalates from 0 to 524 ppm, an observable red shift in the interference dip occurs. However, the magnitude of this wavelength shift is intrinsically linked to the length of the MZI arms.
The sensitivity of a MZI structure increases with longer arm lengths (L=15 µm) compared to shorter ones (L=5 µm) as shown in
Figure 6 (c). Firstly, longer arms provide a greater interaction length for the propagating optical signal within the device, allowing more opportunities for the signal to experience phase changes induced by external stimuli, such as refractive index variations induced by the change in the properties of PHMB deposited on the sensing arm. This increased interaction length enhances the sensitivity of the MZI to these external perturbations. Additionally, longer arm lengths typically result in a larger optical path length imbalance between the two arms of the interferometer. This imbalance amplifies the phase difference between the two arms, leading to more significant interference pattern changes in response to minute changes in the external conditions. Furthermore, longer arms inherently possess a higher spatial resolution, enabling the detection of smaller changes in the environment or sample properties. The sensitivity of the proposed device estimated at 7.63 pm/ppm, 34.46 pm/ppm and 74.78 pm/ppm for the arm length (L) of 5 µm, 10 µm and 15 µm, respectively. Furthermore,
Table 3 showcases not only the footprint of the device but also the sensitivity it offers.
The normalized H-field distribution within the LT-MZI structure is vividly illustrated in
Figure 6. It delineates the wavelength where both destructive interference and constructive interference occur, as depicted in
Figure 7 (a) and
Figure 7 (b) respectively. Key geometric parameters of the device are as follows: L
DC=6 µm, L=5 µm, w=250 nm, g=100 nm. Additionally, in this particular scenario, the CO
2 gas concentration is assumed to be negligible.
Table 4 showcases a compilation of recent publications focusing on different types of photonic sensors utilizing PHMB polymer tailored for diverse applications, including subwavelength grating slot waveguides, metasurfaces, and optical fibre interferometers, specifically engineered for CO
2 gas sensing. Notably, the analysis reveals a notable trend: the sensor device proposed based on the LT-MZI structure emerges as a standout performer in terms of its unparalleled responsiveness to CO
2 gas within the ambient environment, surpassing the efficacy of previous works in this domain.
Figure 1.
Schematic representation of LT-MZI structure embedded with PHMB polymer on the sensing arm for CO2 gas sensing. Inset presents the H-field distribution at 1550 nm in the cross-sectional view of the waveguide covered with PHMB polymer.
Figure 1.
Schematic representation of LT-MZI structure embedded with PHMB polymer on the sensing arm for CO2 gas sensing. Inset presents the H-field distribution at 1550 nm in the cross-sectional view of the waveguide covered with PHMB polymer.
Figure 2.
Transmission of a directional coupler segment versus LDC for waveguide width (w) of (a) 250 nm, (b) 260 nm, (c) 270 nm, (d) 280 nm. The normalized H-field distribution in the directional coupler for, (e) LDC=3 µm, and (f) LDC= 4 µm. Note: The gap between the waveguides is fixed at 100 nm.
Figure 2.
Transmission of a directional coupler segment versus LDC for waveguide width (w) of (a) 250 nm, (b) 260 nm, (c) 270 nm, (d) 280 nm. The normalized H-field distribution in the directional coupler for, (e) LDC=3 µm, and (f) LDC= 4 µm. Note: The gap between the waveguides is fixed at 100 nm.
Figure 3.
Transmission spectrum of the LT-MZI structure versus LDC, (a) 3 µm, (b) 3.5 µm, (c) 4 µm, (d) 4.5 µm, (e) 5 µm, (f) 5.5 µm, (g) 6 µm. Note: w=250 nm, g=100 nm, L=5 µm.
Figure 3.
Transmission spectrum of the LT-MZI structure versus LDC, (a) 3 µm, (b) 3.5 µm, (c) 4 µm, (d) 4.5 µm, (e) 5 µm, (f) 5.5 µm, (g) 6 µm. Note: w=250 nm, g=100 nm, L=5 µm.
Figure 4.
Transmission spectrum of LT-MZI versus w, (a) 250 nm, (b) 260 nm, (c) w=270 nm. Note: LDC=6 µm, g=100 nm, L=5 µm. Transmission spectrum of LT-MZI versus g, (d) 100 nm, (e) 125 nm, (f) 150 nm, (g) 175 nm. Note. LDC=6 µm, L= 5 µm, w=250 nm.
Figure 4.
Transmission spectrum of LT-MZI versus w, (a) 250 nm, (b) 260 nm, (c) w=270 nm. Note: LDC=6 µm, g=100 nm, L=5 µm. Transmission spectrum of LT-MZI versus g, (d) 100 nm, (e) 125 nm, (f) 150 nm, (g) 175 nm. Note. LDC=6 µm, L= 5 µm, w=250 nm.
Figure 5.
(a) Schematic representation of MZI structure without end loop, (b) Transmission spectrum versus CO2 gas concentration, (c) Power decay at Output_1 for the operational wavelength of 1535.5 nm versus CO2 gas concentration. Normalized H-field distribution in the device at the wavelength of 1535.5 nm for the gas concentration of, (d) 0 ppm, and (e) 366 ppm.
Figure 5.
(a) Schematic representation of MZI structure without end loop, (b) Transmission spectrum versus CO2 gas concentration, (c) Power decay at Output_1 for the operational wavelength of 1535.5 nm versus CO2 gas concentration. Normalized H-field distribution in the device at the wavelength of 1535.5 nm for the gas concentration of, (d) 0 ppm, and (e) 366 ppm.
Figure 6.
Transmission spectrum of LT-MZI structure in the presence of CO2 gas with different concentrations. The length (L) of the sensing arm is, (a) 5 µm, (b) 10 µm. (c) Sensitivity of the device for L=5 µm, 10, 15 µm. Note: The remaining geometric parameters are, LDC=6 µm, w= 250 nm, and g=100 nm.
Figure 6.
Transmission spectrum of LT-MZI structure in the presence of CO2 gas with different concentrations. The length (L) of the sensing arm is, (a) 5 µm, (b) 10 µm. (c) Sensitivity of the device for L=5 µm, 10, 15 µm. Note: The remaining geometric parameters are, LDC=6 µm, w= 250 nm, and g=100 nm.
Figure 7.
Normalized H-field distribution in the LT-MZI structure at, (a) λ=1559.8 nm, (b) 1566 nm.
Figure 7.
Normalized H-field distribution in the LT-MZI structure at, (a) λ=1559.8 nm, (b) 1566 nm.
Table 1.
Geometric parameters of LT-MZI structure.
Table 1.
Geometric parameters of LT-MZI structure.
Variable |
Expression |
Value |
w |
Width of the waveguide |
250,260,270 nm |
g |
The gap between waveguides of directional coupler |
100,125,150,175 nm |
R |
Radius |
2 µm |
LDC
|
Length of directional coupler |
0.5 to 12 µm |
L |
Length of MZI structure |
1,3,5,7,10 µm |
Table 2.
Refractive index of PHMB versus CO
2 gas concentration[
27,
28].
Table 2.
Refractive index of PHMB versus CO
2 gas concentration[
27,
28].
PHMB refractive index |
Gas concentration (ppm) |
1.55 |
0 |
1.54 |
215 |
1.53 |
262 |
1.52 |
328 |
1.51 |
366 |
1.49 |
434 |
1.48 |
524 |
Table 3.
Device footprint versus sensitivity.
Table 3.
Device footprint versus sensitivity.
L (µm) |
Device footprint (µm×µm) |
Sensitivity (pm/ppm) |
5 |
45×16 |
7.63 |
10 |
55×16 |
34.46 |
15 |
65×16 |
74.78 |
Table 4.
Sensitivity comparison of the proposed sensor design with previously published works.
Table 4.
Sensitivity comparison of the proposed sensor design with previously published works.
Device design |
Sensitivity |
Ref. |
Silicon subwavelength grating slot waveguide |
12.9 pm/ppm |
[30] |
Silicon nano-cylinders metasurface |
17.3 pm/ppm |
[29] |
Microbubble resonator |
0.49 pm/ppm |
[28] |
Silicon microring refractometric sensor |
6×10-9 RIU/ppm |
[26] |
Optical fibre Fabry-Perot Interferometer |
12.2 pm/ppm |
[26] |
Plasmonic perfect absorber |
17.32 pm/ppm |
[31] |
Loop-terminated MZI structure |
7.63 pm/ppm for L=5 µm; 34.46 pm/ppm for L=10 µm; 74.78 pm/ppm for L=15 µm |
This work |