Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
Early Diagnosis of Bronchopulmonary Dysplasia with E-Nose: A Pilot Study in Preterm Infants
Previous Article in Journal
Advanced 3D Face Reconstruction from Single 2D Images Using Enhanced Adversarial Neural Networks and Graph Neural Networks
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Frequency Asymmetric Absorption–Transmission Metastructures–Photonic Crystals and Their Application as a Refractive Index Sensor

College of Electronic and Optical Engineering & College of Flexible Electronics (Future Technology), Nanjing University of Posts and Telecommunications, Nanjing 210023, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(19), 6281; https://doi.org/10.3390/s24196281
Submission received: 12 August 2024 / Revised: 17 September 2024 / Accepted: 18 September 2024 / Published: 28 September 2024
(This article belongs to the Section Communications)
Figure 1
<p>Schematic diagram of the MPC: (<b>a</b>) ordinary periodic structure; (<b>b</b>) introducing the analyte layer.</p> ">
Figure 2
<p>Forward absorption and backward transmission curves in ordinary periodic MPC.</p> ">
Figure 3
<p>(<b>a</b>) Forward absorption and backward transmission curves of the MPC introduced into the analyte layer. The black solid line represents the forward absorption curve, the red solid line indicates the backward transmission curve, and the black dashed arrow marks the resonance frequency. (<b>b</b>) Energy distribution of forward and backward incident electric fields at resonance frequency points.</p> ">
Figure 4
<p>Forward absorption curves corresponding to different RI: (<b>a</b>) frequency range from 12 to 18 GHz, (<b>b</b>) frequency range from 17.00 to 17.25 GHz, and (<b>c</b>) relationship between wavelength and forward absorption curves.</p> ">
Figure 5
<p>When <span class="html-italic">N</span> = 3, (<b>a</b>) shows the relationship between the resonance wavelength and RI, and (<b>b</b>) shows the linear fitting of the resonance frequency <span class="html-italic">f</span> and resonance wavelength <span class="html-italic">λ</span> to RI.</p> ">
Figure 6
<p>When <span class="html-italic">N</span> = 4, (<b>a</b>) shows the change in the forward absorption curve with <span class="html-italic">n<sub>c</sub></span>, and (<b>b</b>) shows the linear fitting of the resonance frequency with RI. The blue dots represent the sampling points at various RI, while the black dashed lines indicate the fitting curves.</p> ">
Figure 7
<p>When <span class="html-italic">N</span> = 5, (<b>a</b>) shows the change in the forward absorption curve with <span class="html-italic">n<sub>c</sub></span>, and (<b>b</b>) shows the linear fitting of the resonance frequency with RI. The orange dots represent the sampling points at various RI, while the black dashed lines indicate the fitting curves.</p> ">
Figure 8
<p>Resonance frequency and period number <span class="html-italic">N</span> = 3, <span class="html-italic">N</span> = 4, and <span class="html-italic">N</span> = 5 are related in various cases.</p> ">
Figure 9
<p>The linear fitting of the resonance frequency with RI at (<b>a</b>) <span class="html-italic">φ</span> = 30° and (<b>b</b>) <span class="html-italic">φ</span> = 50°. The green and blue stars represent the sampling points for the different RI, while the black dashed lines indicate the fitted curves.</p> ">
Figure 10
<p>Resonance frequency and tilt angles <span class="html-italic">φ</span> = 30°, <span class="html-italic">φ</span> = 40°, and <span class="html-italic">φ</span> = 50° are related in various cases.</p> ">
Versions Notes

Abstract

:
In this paper, a kind of metastructure–photonic crystal (MPC) with multi-frequency asymmetric absorption–transmission properties is proposed. It is composed of various dielectric layers arranged in a periodically tilting pattern. When electromagnetic waves (EMWs) enter from the opposite direction, MPC shows an obvious asymmetry. EMWs are absorbed at 13.71 GHz, 14.37 GHz, and 17.10 GHz in forward incidence, with maximum absorptions of 0.919, 0.917, and 0.956, respectively. In the case of backward incidence, transmission above 0.877 is achieved. Additionally, the MPC is utilized for refractive index (RI) sensing, allowing for wide RI range detection. The refractive index unit is denoted as RIU. The RI detection range is 1.4~3.0, with the corresponding absorption peak variation range being 17.054~17.194 GHz, and a sensitivity of 86 MHz/RIU. By adjusting the number of MPC cycles and tilt angle, the sensing performance and operating frequency band can be tailored to meet various operational requirements. This MPC-based RI sensor is simple to fabricate and has the potential to be used in the development of high-performance and compact sensing devices.

1. Introduction

Asymmetric devices refer to electromagnetic waves (EMWs) that exhibit different propagation behaviors when incident from opposite directions, such as absorption, reflection, and transmission [1,2,3]. This bidirectional transmission disparity enables asymmetric devices to find applications in sensing [4], all-optical computing [5], optical signal processing [6], and multiplexing [7]. Researchers have explored various methods to create asymmetry, with building artificial composite structures like asymmetric gratings [8], photonic crystals [9], and metamaterials [10,11] being widely recognized as effective in breaking spatial inversion symmetry [12]. In recent years, the potential use of asymmetric propagation in integrated photonic systems has turned it into a thriving research area. Xu et al. [11] utilized a pair of asymmetric subwavelength gratings and passive hyperbolic metamaterials to achieve visible frequency asymmetry. By coupling the grating in the material’s propagation mode, they achieved asymmetric light transmission with a contrast exceeding 14 dB at central wavelengths of 532 nm and 633 nm. Liu et al. [13] introduced a dual-frequency asymmetric device based on a double-sided composite metal grating, created by connecting a silver grating on the outer surface with an inward inverted T-shaped grating with a Ag film. The device achieved high transmittances of 0.67 and 0.82 at operating wavelengths of 627 nm and 1238 nm, respectively. However, current research primarily focuses on broadband or single-frequency cases, where the opposite direction exhibits similar propagation characteristics. In some other applications, achieving multi-frequency different propagation characteristics of the asymmetric function is crucial.
As a new material design concept, metastructures consist of artificially designed subwavelength units [14] that can acquire properties not found in traditional materials through special multi-scale configuration designs [15], such as a negative refractive index and negative dielectric constant [16]. Due to their unique capability to control the propagation direction, amplitude, and polarization state of EMWs, metastructures can potentially replace traditional bulky optical devices, facilitating the integration and miniaturization of optical devices [17]. Photonic crystals (PCs) can be categorized into one-dimensional, two-dimensional, and three-dimensional structures based on the spatial distribution of dielectric materials [18]. The dielectric constant of PCs exhibits periodicity in space, resulting in a periodic distribution of light refraction. This unique structure gives PCs distinct physical properties like photonic band gaps and photon locality [19], enhancing the interaction between light and matter through effects such as slow optical response [20] and resonator resonance [21]. PCs possess unique capabilities such as reflection, conduction, beam splitting, and coupling to EMWs [22], making them essential materials for micro-nano optoelectronic integration. However, most reported PCs currently suffer from poor structural adjustability and low reflectivity due to their densely packed structure [23]. Metastructures–photonic crystals (MPCs) are a new type of artificial electromagnetic material that combines metastructures with the concept of traditional PCs. By precisely adjusting the geometric and material properties of MPC, its flexibility and energy localization can be enhanced to achieve a specific response to EMWs.
Because of their small size, high sensitivity, and strong resistance to electromagnetic interference, optical sensors have garnered significant attention from researchers, making them a focal point in refractive index (RI) measurements in recent years [24,25,26]. Currently, various physical mechanisms have been utilized to improve sensing performance, such as surface plasmon resonance [27], Fano resonance [28], and PC defect mode resonance [29,30]. In 2016, Chen et al. [27] first proposed a liquid-core optical fiber surface plasmon resonance sensor with an adjustable nanoporous silica coating. By adjusting the RI of the nanoporous silica coating, the sensor can be used for different RI detection ranges. In 2019, Elnaz et al. [26] proposed a tunable polarimetric sensitive RI sensor with a sensitivity of up to 10,000 nm/RIU and a sensitivity that can be adjusted based on the surface conductivity of graphene-based chiral structures. In 2023, Li et al. [31] proposed a sensor based on hyperbolic metamaterials, which offers an adjustable index interval of analytes for detecting biochemical molecules in the terahertz band. In addition, the specialized design of superstructures to enhance absorption or transmission has been extensively studied for their application in refractive index detection. Li et al. [32] utilized a graphene-based silicon carbide grating and Tamm plasmon structure to achieve a dual-band absorption peak that is highly sensitive to the refractive index of the air layer. Hu et al. [33] excited the Tamm plasmon polariton resonance based on a graphene-based hybrid Tamm plasmon structure to develop a multi-channel absorption and refractive index sensor. These studies play a crucial role in enhancing absorption or transmission, providing significant advancements in the field of optical sensing. Due to the rapid development of the biochemical industry, the demand for sensing and detection continues to rise. Therefore, it is necessary to promote the use of sensing materials and diversify device design.
In this paper, a kind of plasma-based MPC is proposed to achieve the asymmetric propagation of EMWs and wide-range RI detection, enriching the design possibilities for asymmetric devices and sensors. The MPC consists of materials with varying RIs arranged periodically at an oblique angle. Due to the asymmetrical spatial distribution of the MPC, EMW propagation characteristics differ significantly depending on the incident direction, exhibiting a noticeable asymmetric absorption–transmission (AAT) phenomenon. When EMWs are incident in the forward direction, a high absorptivity of 0.94 is achieved at the operating frequency of 24.92 GHz, while in the backward case, a high transmittance of 0.98 is achieved. Even after adding an analyte layer, the AAT phenomenon persists and exhibits multi-frequency behavior. At frequencies of 13.71 GHz, 14.37 GHz, and 17.10 GHz, the forward absorptivity values are 0.919, 0.917, and 0.956, respectively, with backward transmittance above 0.877. To showcase the superiority of the proposed MPC, RI sensing is conducted, demonstrating a wide RI range of 1.4 to 3.0. The absorption peak corresponding to the RI values shifts from 17.194 GHz to 17.054 GHz, with a sensitivity of 86 MHz/RIU. Furthermore, this paper discusses in detail the effects of cycle number and tilt angle on MPC performance, considering potential errors in MPC fabrication and practical usage. The proposed MPC offers advantages such as multifunctionality, compact size, large measurement range, and suitability for asymmetric device design and biosensing applications.

2. The Theoretical Model

The schematic diagram of the proposed MPC is shown in Figure 1. We use M to represent the plasma layer, F to represent the porous silicon layer, and N = 3 to represent the number of periods. In Figure 1a, in the ordinary structure, M and F are arranged periodically, and the arrangement rule is (MF)2N. In Figure 1b, the analyte layer (C) is introduced into the original structure, breaking the periodicity of the structure, and its arrangement rule is (MF)NC(MF)N. Two coordinate systems (x, y, z) and (x’, y’, z’) are defined in Figure 1. The coordinate system (x’, y’, z’) is obtained by rotating the coordinate system (x, y, z) clockwise by an angle of φ centered on the y-axis. Two arrows in different directions at the top and bottom indicate the forward and backward incidence of EMWs. In this paper, the initial condition set by incident EMWs is transverse magnetic (TM) wave [34], and TM mode propagation means that there is no magnetic field component in the propagation direction. The background is air, and the permittivity of layers M, F, and C is expressed by εm, εf, and εc, respectively, where εf = 1.3 and εc = 3.7 [35]. The thickness is indicated by dm, df, and dc. In Figure 1a, dm = 0.4 mm, df = 1.0 mm, and the length of the hypotenuse l = 20 mm. In Figure 1b, dm = 0.6 mm, df = 1.0 mm, dc = 1.0 mm, l = 20 mm. φ is the inclination angle of the MPC with respect to the z-axis, θ is the incidence angle of the TM wave, and the specific values are φ = 40°and θ = 50°. Related parameter settings are shown in Table 1.
Due to the introduction of the magnetic field, the plasma is seen as an anisotropic medium, and its dielectric function is described as follows [36]:
ε m = ε 1 2 ε 2 2 ε 1
ε 1 = 1 ω p 2 ( ω + i ν c ) ω [ ( ω + i ν c ) 2 ω c 2 ]
ε 2 = ω p 2 ω c ω [ ( ω + i ν c ) 2 ω c 2 ]
where ωp = (e2ne/0)1/2 is the plasma frequency, ω is the frequency of the incident TM wave, νc is the collision frequency of electrons, ωc = (eB/m) is the plasma electron oscillation frequency, m is the mass of the electron, e is the charge of the electron, and ne is the plasma density. Their respective values are as follows: B = 0.1 T, ne = 1018 m−3, and νc = 0.1ωp [36].
In the ordinary periodic structure, depicted in Figure 1a, the effective relative permittivity tensor of the MPC in the coordinate system xyz’ is [37]
ε = ε x x 0 0 0 ε y y 0 0 0 ε z z
According to the effective medium model [38], MPC can be described as a uniaxial anisotropic medium. εxx’ represents the component of the permittivity tensor parallel to the x’-axis, while εyy’ and εzz’ are the components of the permittivity tensor parallel to the y’-axis. Refer to the following formula for the specific form [38]:
1 ε x x = 1 d m + d f ( d m ε m + d f ε f )
ε y y = ε z z = d m ε m + d f ε f d m + d f
The relation between the coordinate system xyz and the coordinate system xyz′ is as follows [37]
x y z = cos φ 0 sin φ 0 1 0 sin φ 0 cos φ x y z
The effective relative permittivity tensor of the MPC in the coordinate system xyz is
ε = ε x x 0 ε x z 0 ε y y 0 ε z x 0 ε z z
where εxx, εyy, and εzz can be represented as
ε x x = ε x x cos 2 φ + ε y y sin 2 φ ε z z = ε x x sin 2 φ + ε y y cos 2 φ ε x z = ε z x = ε y y ε x x sin φ cos φ
In the coordinate system xyz′, the electric displaces vector D′, and electric filed vector E’ is
D = ε 0 ε E
In the coordinate system xyz, it is
D = ε 0 ε E
The magnetic field has only y-directional components for both the forward and backward incidence of the TM waves, so [37]
H y = H y + y + H y y = H y 0 + e i k z + z + k x x ω t y + H y 0 e i k z z + k x x ω t y
where kz+ and kz indicate the +z and −z components of the TM wave vector.
So, we can obtain [37]
D x 0 D z = ε 0 ε x x 0 ε x z 0 ε y y 0 ε z x 0 ε z z E x 0 E z
According to Maxwell’s equations, the following equation can be derived
k z + H y + + k z H y = ω ε 0 ε x x E x + ε x z E z k x H y + + H y = ω ε 0 ε z x E x + ε z z E z
E x = ε z z k z + + ε x z k x ω ε 0 ε x x ε z z ε x z 2 H y + + ε z z k z + ε x z k x ω ε 0 ε x x ε z z ε x z 2 H y
From × E = B t , for the forward incident wave, we have [37]
E x z E z x = i ω μ 0 μ H y +
The following result can be obtained by combining the above equation:
ε z z k z 2 + ε x z k x + k x ε z x k z + ε x x k x 2 = k 0 2 ε x x ε z z ε x z 2
The two results of solving Equation (17) are as follows
k z 1 = ( ε x z k x + ε x z 2 ε z z ε x x k x 2 k 0 2 ε z z ) / ε z z
k z 2 = ( ε x z k x ε x z 2 ε z z ε x x k x 2 k 0 2 ε z z ) / ε z z
where k0 = ω/c and kx = k0 sinθ [37].
Equations (18) and (19) reveal the root cause of AAT characteristics. They represent the forward incidence and backward incidence of EMWs, respectively. For ordinary dielectric materials, kz1 and kz2 are inversely related, and the structure can be viewed as reversible. However, in the MPC, kx is a non-zero constant because the spatial symmetry is broken. Therefore, the values of kz1 and kz2 are no longer negative, resulting in a significant difference. This leads to different propagation modes of EMWs in the forward and reverse incidents [37].
In addition, the quantitative evaluation of RI sensors involves various indicators, such as sensitivity (S), figure of merit (FOM), and quality factor (Q) [39].
One of the metrics used to evaluate the performance of a proposed sensor is FOM. This measurement is calculated by the ratio of the half-peak width to the peak resonance wavelength [39].
F O M = S F W H M
Among it, FWHM stands for full width at half maximum, which represents the bandwidth of the RI sensor response and indicates the range of signals that the sensor can effectively detect. A narrower FWHM generally indicates higher sensitivity and a higher FOM.
Q is used to measure the energy stored in the MPC relative to the energy lost through dissipative mechanisms such as absorption and scattering. A high Q value indicates that the MPC performs well in efficiently detecting RI changes in the surrounding environment. Mathematically, Q can be defined as follows [39]:
Q = f 0 F W H M
where f0 is the resonance frequency.

3. Analysis and Discussion

The simulation results presented in this paper were obtained using the HFSS method and calculated through the finite element method. In the simulation, the x- and y-directions are subject to periodic boundary conditions, while the z-direction is treated as an open space.
Figure 2 illustrates the EMW propagation characteristic curves corresponding to the ordinary periodic MPC in Figure 1a at frequencies ranging from 22 GHz to 28 GHz. It is evident from the figure that there is a significant difference in the propagation of the MPC when EMWs are incident in the forward and backward directions. Upon forward incidence, EMWs exhibit absorption peaks at 23.04 GHz, 24.92 GHz, and 27.84 GHz, respectively. At the three distinct resonance frequencies, the forward absorption rate can reach 0.34, 0.79, and 0.94, respectively. When the incident light is in the backward case, as Equations (18) and (19), the propagation mode will change. During this scenario, EMWs are predominantly transmitted, with a transmittance exceeding 0.9 within 22~28 GHz. Notably, at f = 24.92 GHz, the backward transmittance is 0.98, showcasing a very high forward absorption and backward transmittance, demonstrating excellent AAT characteristics.
After introducing the analyte layer into the MPC, the structure diagram is depicted in Figure 1b. Figure 3 illustrates the forward absorption and backward transmission curves in this scenario, revealing the continued presence of the AAT phenomenon, which is significantly enhanced. Under forward incidence, EMWs still trigger three absorption peaks at 13.71 GHz, 14.37 GHz, and 17.10 GHz, respectively. However, unlike Figure 2, the absorption rates corresponding to the three resonant frequencies are all above 0.9, specifically 0.919, 0.917, and 0.956, respectively. The backward transmittance rates are 0.877, 0.896, and 0.953, respectively, maintaining high levels. The root cause of the asymmetric propagation illustrated in Figure 3a is explained by Equations (18) and (19). Due to the non-uniform distribution of the MPC dielectric function, the components of the dielectric constant tensor vary in different directions. This variation alters the values of kz1 and kz2, causing EMWs to experience different propagation modes when incident from opposite directions. Additionally, the MPC possesses a specific structural eigenfrequency. When the frequency of the incident EMW aligns with this intrinsic frequency, the absorption and transmission of the EMW by the MPC change significantly. For instance, at the resonance frequency, the electric field energy reaches its maximum, making the EMW more readily absorbed and resulting in a peak in the absorption curve. Figure 3b clearly illustrates the energy distribution of EMWs within the MPC, revealing that the area of the concentrated electric field energy corresponds to the trough in the absorption curve. These outcomes demonstrate that the proposed MPC exhibits multi-frequency AAT characteristics following the introduction of the defect layer, offering a fresh approach for designing asymmetric devices across various frequencies.
Currently, RI sensors have increasingly extensive applications in biomedicine, environmental monitoring, and other fields. By detecting changes in the RI of the medium being measured, the resonance condition changes, leading to variations in wavelength or frequency, enabling the detection of RI. The forward absorption depicted in Figure 3 exhibits sharp spikes. Upon calculation, the Q values of the absorption peaks at the resonance frequencies of 13.71 GHz, 14.37 GHz, and 17.10 GHz are 93.3, 131.8, and 229.8, respectively. Notably, the absorption peak at 17.10 GHz boasts the highest Q value, indicating potential for sensor applications.
To evaluate the performance of the proposed MPC-based RI sensor, the impact of the different RI of the analyte layer on the forward absorption curves is calculated while keeping other parameters constant (as shown in Table 1). The relevant calculation results are presented in Figure 4. The RI of the analyte layer is denoted by nc, with a measurement range of 1.2 RIU to 3.8 RIU. To observe the influence of different RI on forward absorption, values of nc = 1.4 RIU, 1.8 RIU, 2.2 RIU, 2.6 RIU, and 3.0 RIU are selected for analysis. Figure 4a illustrates the variation in the forward absorption curves with changing RI within the range of 12 to 18 GHz. For clearer observation, Figure 4b displays the changes in the value corresponding to the forward absorption peak in our study of 17~17.25 GHz. It is evident that as nc increases, the peak of forward absorption will redshift, and the absorption rate will gradually increase. Specifically, when nc = 1.4, the forward absorption peaks at 17.194 GHz with an absorption rate of 0.12. As nc increases to 1.8, 2.2, and 2.6, the absorption peaks shift to 17.150 GHz, 17.114 GHz, and 17.082 GHz, respectively, with corresponding absorption rates of 0.43, 0.74, and 0.92. At nc = 3.0, the peak absorption rate reaches 0.99 at 17.054 GHz. Figure 4c illustrates the relationship between wavelength and the absorption curves. As the value of nc increases, the forward absorption peak shifts toward longer wavelengths, demonstrating a redshift phenomenon. When nc equals 3.0, the peak absorption rate reaches 0.99 at a wavelength of 1.76 mm.
Figure 5a illustrates the resonance frequency ranging from 17.220 GHz to 17.008 GHz as the RI of the analyte layer C varies from 1.2 RIU to 3.8 RIU. The results indicate that the resonance frequency of MPC varies with the RI of the analyte layer C and maintains excellent linearity within the range of 1.4~3.0 RIU. Figure 5b depicts the relationship between RI and resonance frequency, showing a negative correlation as the resonance frequency decreases with the increase in nc. The RI detection range is 1.4~3.0 RIU, with the frequency changing from 17.194 GHz to 17.054 GHz. The linear fitting equation is f = −0.086nc + 17.31, and the S is 86 MHz/RIU. Simultaneously, its linearity R2 is 0.9924, representing the degree of fit. When nc = 3.0 RIU, the Q value is 315.8, and the FOM value is 1.66 RIU−1. Of course, when considering the relationship between wavelength and RI, the corresponding values are given by the equation λ = 0.0088nc + 1.73, 8.8 µm/RIU and 0.9923.
In order to evaluate the effect of the number of periods N on the performance of the sensor, we calculate the forward absorption curves for different refractive indices at N = 4 and N = 5. When N = 4, it is shown in Figure 6. Figure 6a illustrates the change in the absorption peak as nc increases. Overall, it can be observed that as nc increases, the absorption peak shifts to lower frequencies, and the absorption rate increases. When nc = 1.4 RIU, the forward absorption peaks at 18.035 GHz with an absorption rate of 0.09. With an increase in nc to 3.0 RIU, there is a peak absorption rate of 0.96 at 17.844 GHz. In Figure 6b, the resonance frequency and RI are linearly fitted. The detection range of RI remains at 1.4~3.0 RIU, and the resonance frequency shows a negative correlation with RI. The linear fitting equation is f = −0.119nc + 18.20, with an S of 119 MHz/RIU and a linearity R2 of 0.9974. When nc = 3.0 RIU, the Q value is 302.4, and the FOM value is 2.02 RIU−1. When N = 5, the absorption curves are depicted in Figure 7. In Figure 7a, when nc = 1.4 RIU, the peak of forward absorption occurs at 18.639 GHz, with an absorption rate of 0.12. With an increase in nc to 3.0 RIU, there is a peak absorption rate of 0.98 at 18.514 GHz. In Figure 7b, RI detection ranges from 1.4 RIU to 3.0 RIU, and the frequency ranges from 18.639 GHz to 18.514 GHz. The linear fitting equation is f = −0.078 nc + 18.74, with an S of 78 MHz/RIU and a linearity R2 of 0.9944. When nc = 3.0, the Q value is 363.0, and the FOM value is 1.53 RIU−1.
Under different period numbers N, the correlation between resonance frequency and RI is shown in Figure 8, where N = 3, 4, and 5 are represented by red, blue, and orange, respectively. It is evident that as N increases, the resonance frequency also increases. Based on the analysis, altering the period number N will impact the resonance frequency and RI values. Therefore, selecting the appropriate N is essential for enhancing sensor performance.
In addition, we also studied the influence of structural tilt angle φ on sensor performance. The relationship between resonance frequency and RI under different tilt angles is shown in Figure 9. In Figure 9a, when φ = 30°, the RI detection range is 1.4~3.0 RIU, and the frequency range is 16.006~16.392 GHz. The linear fitting equation is f = −0.242nc + 16.73, the S value is 242 MHz/RIU, and the linearity R2 is 0.9999. When nc = 3.0, the Q value is 76.2 and the FOM value is 1.15 RIU−1, indicating poor sensor performance. In Figure 9b, with φ = 50°, the frequency range is 18.09~18.43 GHz. The linear fitting equation is f = −0.212 nc + 18.72, S is 212 MHz/RIU, and the linearity R2 is 0.9983. When nc = 3.0, the Q value is 335 and the FOM value is 3.93 RIU−1, showing an improvement in performance compared to φ = 30°.
The relationship between resonance frequency and RI at different tilt angles is illustrated in Figure 10, where φ = 30°, 40°, and 50° are denoted by green, red, and light blue, respectively. As φ increases from 30° to 40°, there is a significant rise in the resonance frequency. At φ = 40° and φ = 50°, the frequencies are more closely aligned. The sensor’s performance can be manipulated by adjusting the structural tilt angle.
Table 2 presents a comparison with other studies. It is evident that while the MPC-based RI sensor in this paper may not excel in performance, its standout feature lies in its versatility as an asymmetric device, a unique characteristic not observed in prior research findings.

4. Conclusions

In summary, this paper presents a novel kind of MPC design that enables the asymmetric propagation of EMWs and wide-range RI detection. The MPC demonstrates unique electromagnetic responses at positive and negative incidences, with a forward absorptivity of 0.94 and a backward transmittance exceeding 0.98 at a resonance frequency of 24.92 GHz. This distinct feature is further evident in multi-frequency AAT after the introduction of the analyte layer, offering detection capabilities. In RI sensing applications, the MPC achieves a broad detection range from 1.4 to 3.0 RIU with a significant absorption peak frequency shift, validating its potential for RI detection. The MPC not only showcases advantages such as versatility, compact size, and a wide measurement range but also indicates promising applications in asymmetric optical components, biosensing, and other advanced technology fields, paving the way for the development of high-performance optical devices.

Author Contributions

L.L.: investigation, data curation, formal analysis, writing—original draft, and visualization. X.L.: software and validation. H.Z.: conceptualization, methodology, supervision, and writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Samples of the compounds are available from the authors.

Acknowledgments

This work was supported by the National College Students Innovation and Entrepreneurship Training Program (Grant No. 202310293015Z) and College Student Innovation Training Program of Nanjing University of Posts and Telecommunications.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Li, Z.; Liu, W.; Cheng, H.; Chen, S.; Tian, J. Tunable dual-band asymmetric transmission for circularly polarized waves with graphene planar chiral metasurfaces. Opt. Lett. 2016, 41, 3142–3145. [Google Scholar] [CrossRef] [PubMed]
  2. Li, J.S.; Bai, F.Q. Dual-band terahertz polarization converter with high-efficiency asymmetric transmission. Opt. Mater. Express 2020, 10, 1853–1861. [Google Scholar] [CrossRef]
  3. Liu, D.J.; Xiao, Z.Y.; Wang, Z.H. Multi-Band Asymmetric Transmission and 90° Polarization Rotator Based on Bi-Layered Metasurface with F-Shaped Structure. Plasmonics 2017, 12, 445–452. [Google Scholar] [CrossRef]
  4. Tao, H.; Chieffo, L.R.; Brenckle, M.A.; Siebert, S.M.; Liu, M.; Strikwerda, A.C.; Fan, K.; Kaplan, D.L.; Zhang, X.; Averitt, R.D.; et al. Metamaterials on Paper as a Sensing Platform. Adv. Mater. 2011, 23, 3197–3201. [Google Scholar] [CrossRef]
  5. Wong, R.Y.M.; Tso, C.Y.; Chao, C.Y.H.; Huang, B.; Wan, M.P. Ultra-broadband asymmetric transmission metallic gratings for subtropical passive daytime radiative cooling. Sol. Energy Mater. Sol. Cells 2018, 186, 330–339. [Google Scholar] [CrossRef]
  6. Stephen, L.; Yogesh, N.; Subramanian, V. Broadband asymmetric transmission of linearly polarized electromagnetic waves based on chiral metamaterial. J. Appl. Phys. 2018, 123, 033103. [Google Scholar] [CrossRef]
  7. Shang, G.; Guan, C.; Zhang, K.; Wu, Q.; Liu, J.; Ding, X.; Li, H.; Burokur, S.N.; Ding, X. Design of a frequency-multiplexed metasurface with asymmetric transmission. Opt. Lett. 2022, 47, 4504–4507. [Google Scholar] [CrossRef]
  8. Li, S.; Huang, L.R.; Ling, Y.H.; Liu, W.B.; Ba, C.F.; Li, H.H. High-performance asymmetric optical transmission based on coupled complementary subwavelength gratings. Sci. Rep. 2019, 9, 17117. [Google Scholar] [CrossRef]
  9. Mandatori, A.; Bertolotti, M.; Sibilia, C. Asymmetric transmission of some two-dimensional photonic crystals. J. Opt. Soc. Am. B 2007, 24, 685–690. [Google Scholar] [CrossRef]
  10. Chen, K.; Feng, Y.; Cui, L.; Zhao, J.; Jiang, T.; Zhu, B. Dynamic control of asymmetric electromagnetic wave transmission by active chiral metamaterial. Sci. Rep. 2017, 7, 42802. [Google Scholar] [CrossRef]
  11. Xu, T.; Lezec, H.J. Visible-frequency asymmetric transmission devices incorporating a hyperbolic metamaterial. Nat. Commun. 2014, 5, 4141. [Google Scholar] [CrossRef]
  12. Sheikh Ansari, A.; Iyer, A.K.; Gholipour, B. Asymmetric transmission in nanophotonics. Nanophotonics 2023, 12, 2639–2667. [Google Scholar] [CrossRef]
  13. Liu, H.; Zhang, Y.; Chen, C.; Bai, B.; Zhang, X.; Zhang, H.; Shao, Q.; Wu, B.; Zhao, B. Dynamically adjustable and high-contrast asymmetric optical transmission based on bilateral compound metallic gratings. Opt. Laser Technol. 2021, 140, 107033. [Google Scholar] [CrossRef]
  14. Ebbesen, T.W.; Lezec, H.J.; Ghaemi, H.F.; Thio, T.; Wolff, P.A. Extraordinary optical transmission through sub-wavelength hole arrays. Nature 1998, 391, 667–669. [Google Scholar] [CrossRef]
  15. Fan, R.H.; Xiong, B.; Peng, R.W.; Wang, M. Constructing Metastructures with Broadband Electromagnetic Functionality. Adv. Mater. 2020, 32, 1904646. [Google Scholar] [CrossRef] [PubMed]
  16. Chen, H.R.; Cheng, Y.Z.; Zhao, J.C.; Mao, X.S. Multi-band terahertz chiral metasurface with giant optical activities and negative refractive index based on T-shaped resonators. Mod. Phys. Lett. B 2018, 32, 1850366. [Google Scholar] [CrossRef]
  17. An, Q.; Li, D.; Liao, W.; Liu, T.; Joralmon, D.; Li, X.; Zhao, J. A Novel Ultra-Wideband Electromagnetic-Wave-Absorbing Metastructure Inspired by Bionic Gyroid Structures. Adv. Mater. 2023, 35, 2300659. [Google Scholar] [CrossRef]
  18. Tang, G.J.; He, X.T.; Shi, F.L.; Liu, J.W.; Chen, X.D.; Dong, J.W. Topological photonic crystals: Physics, designs, and applications. Laser Photonics Rev. 2022, 16, 2100300. [Google Scholar] [CrossRef]
  19. Butt, M.; Khonina, S.N.; Kazanskiy, N. Recent advances in photonic crystal optical devices: A review. Opt. Laser Technol. 2021, 142, 107265. [Google Scholar] [CrossRef]
  20. Baba, T. Slow light in photonic crystals. Nat. Photonics 2008, 2, 465–473. [Google Scholar] [CrossRef]
  21. Wang, J.; Liu, C.; Wang, F.; Su, W.; Yang, L.; Lv, J.; Fu, G.; Li, X.; Liu, Q.; Sun, T. Surface plasmon resonance sensor based on coupling effects of dual photonic crystal fibers for low refractive indexes detection. Results Phys. 2020, 18, 103240. [Google Scholar] [CrossRef]
  22. Gao, Y.F.; Sun, J.P.; Xu, N.; Jiang, Z.; Hou, Q.C.; Song, H.; Jin, M.C.; Zhang, C. Manipulation of topological beam splitter based on honeycomb photonic crystals. Opt. Commun. 2021, 483, 126646. [Google Scholar] [CrossRef]
  23. Liu, X.; Liu, J.; Wei, B.; Yang, D.; Luo, L.; Ma, D.; Huang, S. Bio-Inspired Highly Brilliant Structural Colors and Derived Photonic Superstructures for Information Encryption and Fluorescence Enhancement. Adv. Sci. 2023, 10, 2302240. [Google Scholar] [CrossRef] [PubMed]
  24. Qi, L.; Zhao, C.L.; Yuan, J.; Ye, M.; Wang, J.; Zhang, Z.; Jin, S. Highly reflective long period fiber grating sensor and its application in refractive index sensing. Sens. Actuators B 2014, 193, 185–189. [Google Scholar] [CrossRef]
  25. Rakhshani, M.R.; Mansouri-Birjandi, M.A. High sensitivity plasmonic refractive index sensing and its application for human blood group identification. Sens. Actuators B 2017, 249, 168–176. [Google Scholar] [CrossRef]
  26. Shokati, E.; Asgari, S.; Granpayeh, N. Dual-Band Polarization-Sensitive Graphene Chiral Metasurface and its Application as a Refractive Index Sensor. IEEE Sens. J. 2019, 19, 9991–9996. [Google Scholar] [CrossRef]
  27. Chen, Y.; Li, X.; Zhou, H.; Hong, X.; Geng, Y. Refractive index detection range adjustable liquid-core fiber optic sensor based on surface plasmon resonance and a nano-porous silica coating. J. Phys. D Appl. Phys. 2016, 49, 355102. [Google Scholar] [CrossRef]
  28. Zhang, Y.; Liu, W.; Li, Z.; Li, Z.; Cheng, H.; Chen, S.; Tian, J. High-quality-factor multiple Fano resonances for refractive index sensing. Opt. Lett. 2018, 43, 1842–1845. [Google Scholar] [CrossRef]
  29. Wan, B.F.; Xu, Y.; Zhou, Z.W.; Zhang, D.; Zhang, H.F. Theoretical Investigation of a Sensor Based on One-Dimensional Photonic Crystals to Measure Four Physical Quantities. IEEE Sens. J. 2021, 21, 2846–2853. [Google Scholar] [CrossRef]
  30. Kumar, N.; Suthar, B.; Nayak, C.; Bhargava, A. Analysis of a gas sensor based on one-dimensional photonic crystal structure with a designed defect cavity. Phys. Scr. 2023, 98, 065506. [Google Scholar] [CrossRef]
  31. Li, Z.; Chen, S.; Zhang, Y.; Tong, C.; Guo, X.; Shen, J.; Li, C. Terahertz band range adjustable hyperbolic metamaterial refractive index sensor. Results Phys. 2023, 49, 106477. [Google Scholar] [CrossRef]
  32. Li, Z.; Hu, J.; Wang, Z.; Chen, Y.; Li, M.; Zhou, J.; Wu, J.; Wang, J. Perfect absorption and strong phonon polariton with a graphene-based silicon carbide grating and Tamm plasmonic structures. J. Opt. Soc. Am. B 2022, 39, 2164–2171. [Google Scholar] [CrossRef]
  33. Hu, J.; Huang, Y.; Chen, Y.; Hu, Z.-D.; Wu, J.; Wang, J. High-sensitivity multi-channel refractive-index sensor based on a graphene-based hybrid Tamm plasmonic structure. Opt. Mater. Express 2021, 11, 3833–3843. [Google Scholar] [CrossRef]
  34. Jiang, W.; Liu, J.; Xia, T.; Liu, N.; Liu, Q.H. A Necessary and Sufficient Condition for Having Independent TE and TM Modes in an Anisotropic Waveguide. IEEE Trans. Microwave Theory Tech. 2017, 65, 3660–3670. [Google Scholar] [CrossRef]
  35. Zaky, Z.A.; Ahmed, A.M.; Shalaby, A.S.; Aly, A.H. Refractive index gas sensor based on the Tamm state in a one-dimensional photonic crystal: Theoretical optimisation. Sci. Rep. 2020, 10, 9736. [Google Scholar] [CrossRef] [PubMed]
  36. Wang, Y.; Liu, S.; Zhong, S. Tunable multichannel terahertz filtering properties of dielectric defect layer in one-dimensional magnetized plasma photonic crystal. Opt. Commun. 2020, 473, 125985. [Google Scholar] [CrossRef]
  37. Fang, Y.T.; Zhang, Y.C. Perfect Nonreciprocal Absorption Based on Metamaterial Slab. Plasmonics 2018, 13, 661–667. [Google Scholar] [CrossRef]
  38. Wangberg, R.; Elser, J.; Narimanov, E.E.; Podolskiy, V.A. Nonmagnetic nanocomposites for optical and infrared negative-refractive-index media. J. Opt. Soc. Am. B 2006, 23, 498–505. [Google Scholar] [CrossRef]
  39. Sampath, D.; Narasimhan, V. One-Dimensional Defect Layer Photonic Crystal Sensor for Purity Assessment of Organic Solvents. ACS Omega 2024, 9, 9625–9632. [Google Scholar] [CrossRef]
  40. Sui, P.; Zhang, A.; Pan, F.; Chang, P.; Pan, H.; Liu, F.; Wang, J.; Cao, C. High sensitivity refractive index sensor with wide detection range and high linearity based on LSPR in hollow-core anti-resonance fiber. Opt. Laser Technol. 2022, 155, 108427. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the MPC: (a) ordinary periodic structure; (b) introducing the analyte layer.
Figure 1. Schematic diagram of the MPC: (a) ordinary periodic structure; (b) introducing the analyte layer.
Sensors 24 06281 g001
Figure 2. Forward absorption and backward transmission curves in ordinary periodic MPC.
Figure 2. Forward absorption and backward transmission curves in ordinary periodic MPC.
Sensors 24 06281 g002
Figure 3. (a) Forward absorption and backward transmission curves of the MPC introduced into the analyte layer. The black solid line represents the forward absorption curve, the red solid line indicates the backward transmission curve, and the black dashed arrow marks the resonance frequency. (b) Energy distribution of forward and backward incident electric fields at resonance frequency points.
Figure 3. (a) Forward absorption and backward transmission curves of the MPC introduced into the analyte layer. The black solid line represents the forward absorption curve, the red solid line indicates the backward transmission curve, and the black dashed arrow marks the resonance frequency. (b) Energy distribution of forward and backward incident electric fields at resonance frequency points.
Sensors 24 06281 g003
Figure 4. Forward absorption curves corresponding to different RI: (a) frequency range from 12 to 18 GHz, (b) frequency range from 17.00 to 17.25 GHz, and (c) relationship between wavelength and forward absorption curves.
Figure 4. Forward absorption curves corresponding to different RI: (a) frequency range from 12 to 18 GHz, (b) frequency range from 17.00 to 17.25 GHz, and (c) relationship between wavelength and forward absorption curves.
Sensors 24 06281 g004
Figure 5. When N = 3, (a) shows the relationship between the resonance wavelength and RI, and (b) shows the linear fitting of the resonance frequency f and resonance wavelength λ to RI.
Figure 5. When N = 3, (a) shows the relationship between the resonance wavelength and RI, and (b) shows the linear fitting of the resonance frequency f and resonance wavelength λ to RI.
Sensors 24 06281 g005
Figure 6. When N = 4, (a) shows the change in the forward absorption curve with nc, and (b) shows the linear fitting of the resonance frequency with RI. The blue dots represent the sampling points at various RI, while the black dashed lines indicate the fitting curves.
Figure 6. When N = 4, (a) shows the change in the forward absorption curve with nc, and (b) shows the linear fitting of the resonance frequency with RI. The blue dots represent the sampling points at various RI, while the black dashed lines indicate the fitting curves.
Sensors 24 06281 g006
Figure 7. When N = 5, (a) shows the change in the forward absorption curve with nc, and (b) shows the linear fitting of the resonance frequency with RI. The orange dots represent the sampling points at various RI, while the black dashed lines indicate the fitting curves.
Figure 7. When N = 5, (a) shows the change in the forward absorption curve with nc, and (b) shows the linear fitting of the resonance frequency with RI. The orange dots represent the sampling points at various RI, while the black dashed lines indicate the fitting curves.
Sensors 24 06281 g007
Figure 8. Resonance frequency and period number N = 3, N = 4, and N = 5 are related in various cases.
Figure 8. Resonance frequency and period number N = 3, N = 4, and N = 5 are related in various cases.
Sensors 24 06281 g008
Figure 9. The linear fitting of the resonance frequency with RI at (a) φ = 30° and (b) φ = 50°. The green and blue stars represent the sampling points for the different RI, while the black dashed lines indicate the fitted curves.
Figure 9. The linear fitting of the resonance frequency with RI at (a) φ = 30° and (b) φ = 50°. The green and blue stars represent the sampling points for the different RI, while the black dashed lines indicate the fitted curves.
Sensors 24 06281 g009
Figure 10. Resonance frequency and tilt angles φ = 30°, φ = 40°, and φ = 50° are related in various cases.
Figure 10. Resonance frequency and tilt angles φ = 30°, φ = 40°, and φ = 50° are related in various cases.
Sensors 24 06281 g010
Table 1. Main parameter settings of the model.
Table 1. Main parameter settings of the model.
ParameterSymbolQuantity
Ordinary Periodic MPCMPC Containing the Analyte Layer
Thickness of plasma layerdm0.4 mm0.6 mm
Thickness of porous silicon layerdf1.0 mm1.0 mm
Length of hypotenusel20.0 mm20.0 mm
Angle of inclinationφ40°40°
The incidence angle of EMWsθ50°50°
Thickness of analyte layer dc/1.0 mm
Table 2. Comparison with other studies.
Table 2. Comparison with other studies.
Ref.RI Detection RangeSensitivityFOMMultifunctional Device
[25]1.0~1.12320 nm/RIUNoNo
[29]1.05~2.0534.3 GHz/RIU634No
[31]1.000~1.45125.07 THz/RIU30.81No
[40]1.28~1.4427,890 nm/RIU139.78No
This work1.4~3.086 MHz/RIU1.66Yes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lei, L.; Li, X.; Zhang, H. Multi-Frequency Asymmetric Absorption–Transmission Metastructures–Photonic Crystals and Their Application as a Refractive Index Sensor. Sensors 2024, 24, 6281. https://doi.org/10.3390/s24196281

AMA Style

Lei L, Li X, Zhang H. Multi-Frequency Asymmetric Absorption–Transmission Metastructures–Photonic Crystals and Their Application as a Refractive Index Sensor. Sensors. 2024; 24(19):6281. https://doi.org/10.3390/s24196281

Chicago/Turabian Style

Lei, Lei, Xiang Li, and Haifeng Zhang. 2024. "Multi-Frequency Asymmetric Absorption–Transmission Metastructures–Photonic Crystals and Their Application as a Refractive Index Sensor" Sensors 24, no. 19: 6281. https://doi.org/10.3390/s24196281

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop