An Investigation of High-Performance Pressure Sensor by Employing Polymer Defect-Based 1D Annular Photonic Crystal

1. Introduction

Photonic crystal (PhC) is a brand-new class of optical material, which comprises a periodic arrangement of dielectric constant and can be realized in 1D, 2D, and 3D configurations [1,2,3]. In the last decade, 1D PhC attracted much attention within the scientific community because of its unique capability of controlling the flow of electromagnetic waves inside the crystal. The trapping of light inside the PhC structure is governed by the photonic bandgap (PBG) effect, where a particular spectral band is prohibited to propagate through the structure [4]. Remarkably, by inserting a defect layer in such a structure, its periodicity breaks, leading to the formation of strongly localized allowed states inside the PBG [5]. The PBG acts as a reflection band and, therefore, the defect mode appears to be a sharp passband in the PBG. The characteristics of the defect mode can be controlled by changing the defect layer’s refractive index and thickness. This tunable property of the defect mode makes the defect-based 1D PhC the most researched topic due to its broad-ranging application domains like sensing, filtering, waveguiding, agriculture, food safety, and bio-medical applications [6,7,8,9].

With the ongoing remarkable advancements in nanotechnology, it is feasible to fabricate various geometrical configurations of PhC [10,11,12]. Recently, a new configuration of 1D PhC known as 1D cylindrical multilayer structure has fascinated a lot of researchers worldwide due to its remarkable performance, particularly in sensing and filtering applications [13,14,15]. The cylindrical PhC is referred to as an annular photonic crystal (APC) or circular photonic crystal (CPC). The APC structure can be constructed by arranging alternative layers of different materials in a cylindrical ring fashion [16]. The flow of EM waves in APC was first investigated by M.A. Kaliteevkki and his team in 1999. After that, C.A. Hu et al. proposed a modified transfer matrix method (TMM) using cartesian coordinates for the analysis of transmittance and reflection characteristics of the APC [17]. Additionally, the authors thoroughly explained the fabrication feasibility of APC structure and studied the liquid sensing application. Furthermore, J. Scheuer et al. [18] introduced a ring defect into the APC structure to make annular Bragg lasers. Moreover, a wide variety of materials, such as metal, dielectric, superconductor, metamaterial, etc. can be used to develop APC structures [19,20,21]. Recently, APC structures have been extensively studied for sensing applications. For example, M.M. Abadla et al. theoretically demonstrated a temperature sensor using 1D APC, where the authors studied the shift in the transmittance characteristics of the structure with a temperature rise. They achieved a significant temperature sensitivity of 0.033 nm/^oC [22]. S.K. Srivastava et al. designed a 1D ternary annular PhC, where a unit cell consists of a semiconductor, superconductor, and air. The authors investigated the reflectance spectrum for both the TE and TM modes by changing the starting radius and temperature [23]. A defect-based 1D APC with a configuration of (Si/SiO²)^N/2 (TiO₂) (Si/SiO²)^N/2 is proposed in [24]. With the help of TMM, the authors studied the transmission spectrum to analyze the shift in the defect mode frequency for varying temperatures across the defect layer. In the same vein, A. Mehany et al. designed a porous silicon-based multilayer 1D APC, where different layers contain a different percentage of porosity. The authors studied the shifting nature of the Tamm resonance mode by varying the analyte refractive index from 1 to 1.5 and accomplished a remarkable sensitivity of 6770 nm/RIU [25]. S.K. Srivastava et al. designed a 1D APC by using double negative metamaterials and studied its reflectance spectrum at different azimuthal mode numbers [26]. Similarly, S. Gandhi et al. reported a porous Si-based annular 1D PhC to detect the creatinine concentrations in blood, where they observed a maximum sensitivity of 640.29 nm/RIU and a quality factor in the order of 10⁵. S. Jena et al. demonstrated a hydrostatic pressure sensor using different configurations of regular plannar 1D PhC [27]. In addition, F. Segovia et al. [28] and A. Herrera et al. [29] studied the shifting nature of the defect mode and band gap in relation to a rising hydrostatic pressure by using plannar 1D PhC. In another related study, Liu et al. carried out an experimental investigation on the dual-core photonic crystal fiber for pressure sensing from 0 MPa to 45 MPa [30]. The regulation and detection of such high pressure are quite difficult, which imposes a great challenge in the experimental realization of the device. Although the pressure sensing application has been explored by using plannar 1D PhC and photonic crystal fiber, there is still plenty of improvement potential for sensing performance by investigating the effect of pressure using the APC.

Various materials can be used for creating different layers and defect layers of the PhC. These materials can be semiconductors, dielecctrics, metals, superconductors, metamaterials, etc. Recently, polymer materials have emerged as a potential candidate for designing 1D PhC owing to their exceptional tunable optical properties. By stretching and compressing the polymer materials [31,32], it is convenient to modify their properties, which is beneficial in numerous applications, such as sensors, filters, reflectors, etc. [33,34].

This study is significant research as it provides a couple of novelties in the current research scenario. A thorough analysis of a novel polymer-based 1D APC is provided in the study to be used in the pressure sensing application field. The primary benefit of using a polymer material is that it shows an elasto-optic effect so that the optical properties can be easily modified by applying external pressure. The mainstay of this research is to study the tunable characteristics of the defect mode with pressure variations. Moreover, a complete optimization of different geometrical parameters has been performed to achieve a remarkable sensing performance. Additionally, a 3D plot of the reflectance spectrum is studied to enhance the novelty of this study. It is believed that this work can be a torchbearer in the pressure-sensing application field.

2. Proposed Structure and Mathematical Modelling

The cross-sectional view of the proposed binary defect-based APC is illustrated in Figure 1. The design comprises an inner-most circular region (a starting section), which has a refractive index (RI) of $n_0$ and radius of ${\rho }_0$. Also, alternating layers of Si and SiO₂ are deposited, whose RI is $n_1$ and $n_2$, respectively, and the radius is ${\rho }_1$ and ${\rho }_2$, respectively. A defect layer of RI $n_d$ and radius ${\rho }_d$ is inserted in the proposed APC to realize a defect mode in the transmittance spectrum. The period of APC before and after the defect layer is $\raisebox{1ex}{$N$}\!\left/ \!\raisebox{-1ex}{$2$}\right.$, where N denotes the total period of the structure. Here, the value of N is taken as 10. The entire configuration of the designed structure is ${\left(Si/Si{O}_2\right)}^{\raisebox{1ex}{$N$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}/Defect/{\left(Si/Si{O}_2\right)}^{\raisebox{1ex}{$N$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}$. In this study, the proposed structure is analyzed by considering the defect layer as polymer materials like PS and PMMA. The outer region (the final medium) is the substrate layer, which has RI $n_f$ and radius ${\rho }_f$.

It is assumed that the cylindrical waves emerged from the core and diverged further toward the higher layers. In the cylindrical coordinates, the Helmholtz wave equations are used to mathematically explain the E-type mode (Ez, Hϕ, H$\gamma$), which can be written as follows [35,36]:

$$\frac{1}{\gamma }\frac{\partial }{\partial \gamma }\left(\rho \frac{\partial E_z}{\partial \rho }\right)+\frac{1}{{\gamma }^2}\frac{{\partial }^2{E}_z}{\partial {\varnothing }^2}=-{\epsilon }_0{\mu }_0ϵ\mu w^2{E}_z$$

(1)

By stating $E_z$ as $E_z=F\left(r\right)\theta (\varnothing )$, and readjusting, we can write [37],

$$E_z\left(\gamma ,\varnothing ,z\right)=\left[A{W}_v\left(q\gamma \right)+B{Z}_v\left(q\gamma \right)\right]e^{iv\varnothing }$$

(2)

where $W_v$ and $Z_v$denote the 1st and 2nd order of the Bessel’s function, respectively. The symbol q represents the propagation constant, which is expressed as $q^2={\epsilon }_0{\mu }_0ϵ\mu w^2={k_0}^2{n}^2$. The tangential component is stated as follows:

$$H_{\varnothing }=\frac{i}{\omega \mu }\frac{\partial E_z}{\partial \gamma }=\frac{i}{\omega \mu }\left[A{W^{\text{'}}}_v\left(q\gamma \right)+B{Z^{\prime }}_v\left(q\gamma \right)\right]e^{iv\varnothing }$$

(3)

To simplify the calculation, we consider only the fundamental mode (i.e., v=0), so the simplified field equations are written as follows:

$$E_z=\left[A{W}_0\left(q\gamma \right)+B{Z}_0\left(q\gamma \right)\right]$$

(4)

$$H_{\varnothing }=-\frac{iq}{\omega \mu }\left[A{W}_1\left(q\gamma \right)+B{Z}_1\left(q\gamma \right)\right]=-ip=\left[A{W}_1\left(q\gamma \right)+B{Z}_1\left(q\gamma \right)\right]$$

(5)

The fields can be represented in a matrix form as follows [22]:

$$\left(\begin{array}{c}{E}_z\\ H_{\varnothing }\end{array}\right)=\left(\begin{array}{cc}{W}_0\left(qp\right)& Z_0\left(qp\right)\\ -ip{W}_1\left(qp\right)& -ip{Z}_1\left(qp\right)\end{array}\right)\left(\begin{array}{c}A\\ B\end{array}\right)$$

(6)

where $p=\raisebox{1ex}{$q$}\!\left/ \!\raisebox{-1ex}{$\omega \mu $}\right..$The amplitude factors A and B can be expressed as,

$$\left(\begin{array}{c}A\\ B\end{array}\right)={\left(\begin{array}{cc}{W}_0\left(qp\right)& Z_0\left(qp\right)\\ -ip{W}_1\left(qp\right)& -ip{Z}_1\left(qp\right)\end{array}\right)}^{-1}\left(\begin{array}{c}{E}_z\\ H_{\varnothing }\end{array}\right)$$

(7)

By an equivalent field at $\gamma ={\gamma }_0$ is given as:

$${\left(\begin{array}{c}{E}_0\\ H_0\end{array}\right)}_{\gamma ={\gamma }_0}={\left(\begin{array}{c}{E}_1\\ H_1\end{array}\right)}_{\gamma ={\gamma }_0}=\left(\begin{array}{cc}{W}_0\left(q_1{\gamma }_0\right)& Z_0\left(q_1{\gamma }_0\right)\\ -i{p}_1{W}_1\left(q_1{\gamma }_0\right)& -i{p_{1}Z}_1\left(q_1{\gamma }_0\right)\end{array}\right)\left(\begin{array}{c}{A}_1\\ B_1\end{array}\right)$$

(8)

As the amplitude remains the same within a particular layer, we can write,

$$\left(\begin{array}{c}{A}_1\\ B_1\end{array}\right)={\left(\begin{array}{cc}{W}_0\left(q_1{\gamma }_1\right)& Z_0\left(q_1{\gamma }_1\right)\\ -i{p}_1{W}_1\left(q_1{\gamma }_1\right)& -i{p_{1}Z}_1\left(q_1{\gamma }_1\right)\end{array}\right)}^{-1}{\left(\begin{array}{c}{E}_2\\ H_2\end{array}\right)}_{\gamma ={\gamma }_1}$$

(9)

By putting eq. 9 in eq. 8, we get,

$${\left(\begin{array}{c}{E}_0\\ H_0\end{array}\right)}_{\gamma ={\gamma }_0}=\left(\begin{array}{cc}{W}_0\left(q_1{\gamma }_0\right)& Z_0\left(q_1{\gamma }_0\right)\\ -i{p}_1{W}_1\left(q_1{\gamma }_0\right)& -i{p_{1}Z}_1\left(q_1{\gamma }_0\right)\end{array}\right){\left(\begin{array}{cc}{W}_0\left(q_1{\gamma }_1\right)& Z_0\left(q_1{\gamma }_1\right)\\ -i{p}_1{W}_1\left(q_1{\gamma }_1\right)& -i{p_{1}Z}_1\left(q_1{\gamma }_1\right)\end{array}\right)}^{-1}{\left(\begin{array}{c}{E}_2\\ H_2\end{array}\right)}_{\gamma ={\gamma }_1}=Q_1{\left(\begin{array}{c}{E}_2\\ H_2\end{array}\right)}_{\gamma ={\gamma }_1}$$

(10)

where $Q_1$ is the transfer matrix. Upon simplifying $Q_1$, we can get the following equations:

$$Q_{11}=-\frac{\pi q_1}{2}\left[W_0\left(q_1{\gamma }_0\right)Z_1\left(q_1{\gamma }_1\right)-W_1\left(q_1{\gamma }_1\right)Z_0\left(q_1{\gamma }_0\right)\right]$$

(11a)

$$Q_{12}=\frac{i\pi q_1}{2p_1}\left[W_0\left(q_1{\gamma }_0\right)Z_0\left(q_1{\gamma }_1\right)-W_0\left(q_1{\gamma }_1\right)Z_0\left(q_1{\gamma }_0\right)\right]$$

(11b)

$$Q_{21}=\frac{i\pi p_1{q}_1}{2}\left[W_1\left(q_1{\gamma }_0\right)Z_1\left(q_1{\gamma }_1\right)-W_1\left(q_1{\gamma }_1\right)Z_1\left(q_1{\gamma }_0\right)\right]$$

(11c)

$$Q_{22}=\frac{\pi q_1}{2}\left[W_1\left(q_1{\gamma }_0\right)Z_0\left(q_1{\gamma }_1\right)-W_0\left(q_1{\gamma }_1\right)Z_1\left(q_1{\gamma }_0\right)\right]$$

(11d)

By performing a similar type of analysis for other layers, the electric field at the 1st layer and the Nth layer can be related as,

$${\left(\begin{array}{c}{E}_0\\ H_0\end{array}\right)}_{\gamma ={\gamma }_0}=Q_1{Q}_2{Q}_3....Q_N{\left(\begin{array}{c}{E}_f\\ H_f\end{array}\right)}_{\gamma ={\gamma }_N}=M{\left(\begin{array}{c}{E}_f\\ H_f\end{array}\right)}_{\gamma ={\gamma }_N}$$

(12)

where the symbol M denotes the overall transfer matrix, which can be obtained by multiplying the transfer matrix of the individual layers.

2.1 Pressure-Dependent Refractive Index of Constituent Materials

The photoelasticity effect is the method, which is utilized to investigate the sensor’s pressure, where the material’s refractive index (RI) differs in accordance with different stress/pressure values. Specifically, the stress is associated with field components through the photoelastic constant as the fourth power of the refractive index. This variation in the RI of the material for different pressure values can be studied by using a photoelastic constant. By using the photo-elasticity effect, the dependence of the refractive index on the stress effect can be expressed as follows [38]:

$$\nabla \left(\frac{1}{n^2}\right)=\sum s_{ab}{\sigma }_b$$

(13)

where $s_{ab}$ represents the strain-optic coefficient, n denotes the RI, and ${\sigma }_b$ represents the stress component. For an isotropic medium, the pressure-dependent RI can be written as [39]:

$$n=n_0-\left(c_1+c_2\right)\sigma$$

(14)

where $n_0$ is the material’s RI, which corresponds to zero applied pressure, $c_1$ and $c_2$ can be specified as follows:

$$c_1=\frac{{n_0}^3(s_{11}-2\vartheta s_{12})}{2E}$$

(15)

$$c_2=\frac{{n_0}^3\left[\right(s_{11}-\vartheta \left(s_{11}-s_{12}\right)]}{2E}$$

(16)

where $s_{11}$ and $s_{12}$ signify the strain-optic coefficients, $\vartheta$ represents the Poisson’s ratio and E denotes the Young’s modulus. According to the proposed structure, the values of different constants for the materials Si, SiO₂, PS, and PMMA are obtained from [40,41].

3. Results and Discussion

The suggested 1D APC is designed by using the successive layers of Si and SiO₂ with a period N=10, where a defect layer is introduced in the middle of the APC. The structure is analyzed by separately considering the polymer materials (i.e., PS and PMMA) as the defect layer. Therefore, the proposed two different configurations are ${\left(Si/Si{O}_2\right)}^{\raisebox{1ex}{$N$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}/PS/{\left(Si/Si{O}_2\right)}^{\raisebox{1ex}{$N$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}$ and ${\left(Si/Si{O}_2\right)}^{\raisebox{1ex}{$N$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}/PMMA/{\left(Si/Si{O}_2\right)}^{\raisebox{1ex}{$N$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}$. The thickness of the Si and SiO₂ layer is taken as ${\rho }_1=108nm$ and ${\rho }_2=254nm$, respectively, whereas the thickness of the starting medium is ${\rho }_0=905nm$. It is observed that the thickness of the starting medium (${\rho }_0$) is considerably larger compared to the thickness of the constituent layers (${\rho }_{1}and{\rho }_2$) of the APC. A greater value of ${\rho }_0$ confirmed that the designed APC analogously acts to the plannar 1D PhC. Initially, the thickness of the defect layer is considered as ${\rho }_d=234nm$, and the period (N) is taken as 10. As considered, light incidents normally onto the structure, and the azimuthal number (m) is assumed to be zero in order to avoid the computational complexity. The modified transfer matrix method (TMM) is employed to study the reflectance spectrum of the proposed APC. The polymer materials (PMMA and PS) are comparatively less rigid than the traditional semiconductor materials like Si and SiO₂. Therefore, they possess highly tunable RI characteristics concerning pressure change.

To begin with, the reflectance spectrum of the structure is studied by separately considering the defect layer as PS and PMMA, which is shown in Figure 2. It can be observed that a band gap is obtained from 1145 nm to 2071 nm in the reflectance. Remarkably, a sharp dip (i.e., a defect mode) within the band gap is perceived at 1475 nm and 1450 mm for the PS and PMMA as the defect layer or layers, respectively. The mainstay of this study is to investigate the defect mode’s characteristics concerning various applied pressures across the defect layer.

Next, the pressure-sensing application of the proposed APC has been analyzed by selecting PS as the defect medium. Figure 3 demonstrates the shift in the defect mode inside the band gap by increasing the pressure across the defect layer from 0 GPa to 3 GPa. It is observed that with changes in the pressure, the band gap range remains fixed, whereas the defect mode position undergoes a noticeable shift. In particular, the position (wavelength) of the defect mode is blue-shifted from 1475 nm to 1427 nm to 1382 nm to 1341 nm as the pressure increases from 0 GPa to 1 GPa to 2 GPa to 3 GPa, respectively. Also, it is observed that the intensity of the defect mode decreases with an increase in the applied pressure.

After that, a shift in the position of the defect mode occurred with different applied pressures by considering PMMA as the defect layer. As depicted in Figure 4, the wavelength of the defect mode is moved from 1450 nm to 1414 nm to 1379 nm to 1348 nm, as the pressure increases from 0 GPa to 1 GPa to 2 GPa to 3 GPa, respectively. The primary reason for such a shift (as shown in Figure 3 and Figure 4) in the defect mode position is that the RI of the PS and PMMA layer changes with a change in pressure, which modifies the overall transfer matrix, and results in a definite shift in the wavelength of the defect mode. This phenomenon can be well explained by the following mathematical equation:

$$\sigma =m\lambda =n_{eff}\vartheta$$

(17)

where $\sigma$ and $\vartheta$ denote the optical and geometrical path difference, respectively. $\lambda$ represents the wavelength, $n_{eff}$ is the effective refractive index, and m is an integer. When the pressure across the defect layer increases, the value of $n_{eff}$ increases, thereby making the $\lambda$ shift toward a higher value in order to keep$\sigma$ fixed.

The variation in the position of the defect mode is plotted in relation to changes in the pressure parameter, as shown in Figure 5. It is observed that for the PS and PMMA together as a defect layer, the defect mode wavelength decreases almost linearly with an increase in pressure. For a low pressure (0-2 MPa), the PS-based APC, the defect mode position or positions took place at relatively a higher wavelength compared to the PMMA-based APC. By contrast, for high pressure (i.e., 3 MPa), the PS-based APC showed a higher defect mode wavelength than the PMMA-based APC. The fitting equation for the variation of the defect peak position for constant changes in pressure can be mathematically represented as follows:

$$\begin{gathered} \mathrm{For}\mathrm{PS}:{\lambda }_{defect}=1.94{\sigma }^2-50.96\sigma +1476 \\ \mathrm{For}\mathrm{PMMA}:{\lambda }_{defect}=1.06{\sigma }^2-37.43\sigma +1450 \end{gathered}$$

where ${\lambda }_{defect}$ is in $nm$ and $\sigma$ is in $GPa$.

The variation in the intensity of the defect mode peak is investigated with changes in pressure, as shown in Figure 6. A lower intensity value of the defect mode indicated accurate detection of the defect mode’s position. As shown in the figure, for both the PS and PMMA-based APC, the intensity decreases with a rise in the applied pressure. In particular, the intensity decreases from 0.1699 to 0.104 to 0.0402 to 0.0038 for the PS-based APC and 0.1305 to 0.0825 to 0.0374 to 0.0081 for PMMA-based APC, as the pressure increases from 0 GPa to 1 GPa to 2 GPa to 3 GPa, respectively.

After that, the sensor’s sensitivity, which is the most significant parameter to assess the sensor’s performance, is calculated. The sensitivity is measured by the geometrical parameters at ${\rho }_1=108nm$, ${\rho }_2=254nm$, ${\rho }_d=234nm$ and N=10. As depicted in Figure 7, the sensitivity of the sensor decreases with a rise in pressure. For the PS-based APC configuration, a maximum sensitivity of 49.31 nm/GPa has been achieved at $\sigma$=0 GPa, whereas for the PMMA-based APC structure, a maximum sensitivity of 36.73 nm/GPa is obtained at $\sigma$=0 GPa. Accordingly, such a significant sensitivity of the proposed structure can make it a robust candidate in the pressure sensing application field. In addition to the pressure sensitivity, other important sensing parameters, such as quality factor, FWHM, and figure of merit for both the PS defect-based APC and the PMMA defect-based APC, were also computed, as illustrated in Table 1 and Table 2, respectively.

To provide a deeper understanding of the shifting nature of the defect mode, a 3D plot of the reflectance has been examined for the pressure and wavelength. Figure 8 and Figure 9 represent the reflectance spectrum concerning pressure variation for PS defect and PMMA defect-based APC, respectively. As illustrated, the blue color line indicates the lower and upper edge of the band gap. It can be observed that the defect mode appeared within the band gap, and the wavelength of this defect mode continued to decrease with an increase in the applied pressure. This nature of variation is exactly similar to the variation explained in Fg.5.

To optimize the thickness of the defect layer (${\rho }_d$), the sensitivity and FoM were computed by changing the value of ${\rho }_d$ from 180 nm to 260 nm, and at a fixed ${\rho }_1=108nm$, ${\rho }_2=254nm$, and N=10. It is observed that for 260 nm<${\rho }_d$<180 nm, both the defect mode intensity and full-width half maximum (FWHM) are high, leading to the realization of low sensing parameters. Thus, ${\rho }_d$ is considered within the range of 180 nm and 260 nm only. For each value of ${\rho }_d$, the sensitivity and FoM were measured for the applied pressure within the range of 0 GPa and 3 GPa. As shown in Figure 10, it is observed that with an increase in thickness, sensitivity increases for both PS and PMMA defect-based APC structures. A maximum sensitivity of 51.29 nm/GPa and 38.92 nm/GPa is obtained for the PS and PMMA defect-based APC configuration, respectively. As shown in Figure 11, it is obvious that FoM increases with an increase in ${\rho }_d$. It is also observed that the PS defect-based APC achieved a higher FoM versus the PMMA defect-based APC structure. A maximum FoM of 326.2 1/GPa and 250.6 1/GPa is accomplished for the PS and PMMA defect-based APC configuration, respectively. Based on the above discussion of results, it is confirmed that an optimum sensitivity and FoM are obtained at ${\rho }_d=260nm$. Thus, the optimized value of ${\rho }_d$ is obtained as 260 nm.

Subsequently, the effect of the number of periods of APC on sensitivity and FoM is studied, as illustrated in Figure 12 and Figure 13, respectively. This simulation is carried out by setting ${\rho }_1=108nm$ and ${\rho }_2=254nm$, ${\rho }_d=260nm$, and changing the period (N) from 8 to 12. As shown in Figure 12, there is very little variation in sensitivity with an increase in the value of N. For the PS defect-based APC, the sensitivity is significantly higher compared to the PMMA defect-based APC. On the other hand, the value of N has a substantial impact on the FoM. As depicted in Figure 13, the FoM increases with an increase in the value of N. Like sensitivity, the PS defect-based APC obtained superior FoM compared to the PMMA defect-based APC. Consequently, sensitivity is the most significant sensing parameter, and it remains almost unaffected by changes in the value of N. Thus, the optimized value of N is obtained as 10.

4. Conclusions

In this study, a novel polymer defect-based 1D annular photonic crystal sensor is designed and examined. The proposed sensor is effectively capable of sensing pressure from 0 MPa to 3 MPa. The designed structure comprises alternate layers of Si and SiO₂. This structure is studied by separately considering two polymer defect layers (i.e., PS and PMMA). The modified transfer matrix method is employed to examine the reflectance characteristics of the proposed structure. The sensing principle is based on the tuning of the defect mode wavelength inside the band gap of the reflection spectra with different applied pressures. A 3D plot of the reflectance spectrum of the proposed APC is studied with pressure variation. Under the case of normal incidence, and at optimized structure parameters of ${\rho }_1=108nm$ and ${\rho }_2=254nm$, ${\rho }_d=260nm$, and N=10, the proposed PS defect-based APC provided a maximum sensitivity and FoM of 51.29 nm/GPa and 301.7 GPa^-1, respectively. On the other hand, the PMMA defect-based APC provided a maximum sensitivity and FoM of 38.92 nm/GPa and 231.7 GPa^-1, respectively. Therefore, it is concluded that for the designed structure, PMMA as a polymer material outperformed compared to PS polymer in the pressure sensing application. The fabrication’s feasibility and stability of the reflectance characteristics of the proposed structure made the designed structure in this study a suitable candidate for the pressure sensing application.