1. Introduction
Gas sensors have wide applications in toxic monitoring, explosive detection and pollution control. One method for detecting dangerous gas is based on surface plasmon resonance (SPR) [
1,
2,
3]. When the surface plasmonic polaritons (SPPs) along the metal/dielectric interface are excited with wave vector that matches SPR, sensitive gas sensing can be achieved. On the other hand, due to properties such as low-loss and availability for scaling up and down, bound states in the continuum (BICs) or “dark optical states” sustained within the all-dielectric micro/nanostructures can provide another platform for sensing applications, in which surface polaritons that propagate at the interface of dielectric gratings can be utilized for trapping and manipulating photons [
4,
5,
6].
THz wave, covering the electromagnetic spectrum from 1 to 10 THz, has wide applications in prospects of studying various chemical gases since the rotational or transitional energy of most gaseous molecules lies in such range [
7]. Flammable gases such as N
2O, HCN, H
2 have strong absorption peaks in terahertz bands which are also termed as “fingerprints.” For example, according to rotation-vibration spectra, the absorption peak of hazardous gas molecular HCN lies at 1.24 THz [
8].
Gas sensors based on graphene have been extensively studied since SPPs can be supported in grapheme in the range of THz wave. Compared with noble metals, graphene shows distinguished properties such as less intrinsic loss and the flexible tunability that realized by the chemical and electrical doping methods. Due to the relative large refractive indices of graphene in the THz regime, the figure of merit (FOM) in THz gas sensor is typically large (e.g., above 1150 /RIU in
Table 1); however, the sensitivity in angle (the angle shift as in response to the change of refractive index) is typically small (e.g., 34.11°/RIU in
Table 1). In addition, since the imaginary part of permittivity of graphene is large in the THz range, prism with very high refractive index or an air gap in between the graphene and prism is usually required [
9], which makes the whole sensing system bulky and may hinder the potential application. In this paper, by introducing all polymeric gratings instead of prism on the 1D-PhC, higher FOM and high sensitivity can be numerically achieved at the same time.
Bloch surface wave (BSW) could be classified as one kind of BICs, it has applications in refractive index sensing area [
10,
11,
12]. Quite similar to SPPs [
13], BSWs propagating along the interface between the dielectric layer and the one-dimensional photonic crystals (1DPhCs), strong enhancement in electric field can be produced at the interface due to the bandgap effect. Since there are a variety of dielectric materials that can be selected to construct 1DPhCs, the dispersion of BSW can be arbitrarily tuned, as compared with the case in SPP that provided by limited novel metals [
14]. In addition, because of the low loss property of dielectric materials, relatively narrower resonant peaks can be obtained in BSW-based structures.
In this paper, by combing with the merits of THz wave and BSW, a new kind of THz gas sensor is designed with polymeric materials that are inexpensive to fabricate. According to the simulation results based on finite element method (FEM), the proposed gas sensor is sensitive to the change of refractive index real part of HCN, as indicated by the high FOM measured by frequency ~147.3 /RIU and FOM measured by the incident angle ~301.3°/RIU. The dispersion of BSW has been discussed to explain the sensing mechanism. The above mentioned sensing mechanism is based on the change of the real part of the refractive index, in which the analyte does not show absorption to the incident THz wave. In the following, we present discussions on another way of gas sensing which is sensitive to the change of the imaginary part of the refractive index, in such cases the analyte shows absorption to THz wave. Such as N2O, for gases with an absorption spectrum in the 10–20 THz band, gas sensing applications can be realized. At the end of this paper, comparisons in terms of sensitivity, fabrication materials, operation frequency and gas contents have been conducted and the potential applications of the proposed design have been discussed.
2. Structure and Analysis
The schematic of BSW-based THz gas sensor is shown in
Figure 1. The proposed dielectric multilayer sustains transverse electric (TE)-polarized (transverse magnetic (TM)-polarized light can also be supported) BSWs which has N = 9 alternating dielectric layers of polyvinylidene fluoride (PVDF) (denoted with light green color, with a thickness of
tPVDF = 7 μm) and polycarbonate (PC) (denoted with gray color, with a thickness of
tPC = 5 μm). The buffering layer (with the thickness of
tbuf = 24 μm), which is composed of a high refractive index PC material, is introduced at the bottom of the PhC structure and a coupling grating is introduced under the buffering layer (with grating constant Λ = 10 μm, height of grating
hg = 10 μm and aspect ratio AR = 0.28). Both the buffering layer and the grating are made of PC material. The sensor is surrounded by gas, with refractive index
ng. During operation, when the sensor is placed in the atmosphere of the detected gas, a terahertz wave is incident from the bottom in a direction at an angle θ to the normal. Because the change in the refractive index of the gas produces different degrees of perturbation, the resonance frequency shifts and the parameters of the excited BSW are also changed. The electromagnetic waves are exponentially elapsed in the vertical direction in the photonic crystal layer.
In the terahertz region, the PVDF dielectric function exhibits a resonance,
where ε
opt = 2.0, ε
dc = 50.0, ω
TO = 0.3 THz and γ = 0.1 THz [
15]. In comparison with PVDF, the dielectric response of PC polymer is frequency independent, having a purely real dielectric constant ε
PC = 2.56. The high and low refractive index layers are PC and PVDF materials, respectively and their dielectric constants are 2.56 and 1.9570-0.0004
i at 10 THz.
When material losses are negligible and the number of reflector periods is infinite, it indicates that there exists a frequency range for complete TE polarization reflection for a given angle of incidence, according to the theory of planar periodic reflector SPR. The formula of the wave vector is given as follows:
where
k0 =
2π
/λ and
kBSW are the magnitudes of the free space wave vector and grating BSW wave vector respectively,
n0 is the refractive index of the sensing gas, θ is the incident angle, Λ is the grating period and m is an integer [
16,
17], the negative sign “−” means the propagation direction of the BSW is opposite that of incident
k0 vector that horizontally projected. As THz wave is incident on the grating, it produces the reflected and BSW waves. As according to References [
16,
17,
18], the direction of the BSW wave is opposite to the direction of the x component of the incident wave, it thus originates the negative sign as shown in Equation (2). In this paper,
m = 1 and the incident angle is 0° < θ < 90°. The wavelength range corresponding to a certain grating constant can be obtained as follow:
when Λ = 10 μm, the effective index of refraction 1.4 <
neff < 1.6, from Equation (3). Therefore, the wave numbers of BSW is in the range of 1.4
k0 to 1.6
k0. We can get the operation frequency
f for producing BSW, which lies at the range: 11.54 THz <
f < 21.43 THz, such analytical result is close to the simulation result as indicated by the frequency range of the purple line that shown in
Figure 2a.
Figure 2a,b shows the band diagrams of the periodic PVDF/PC dielectric structure for TE and TM polarized modes. The numerical simulations were performed by means of the finite element method [
19]. The white regions represent the 1D PhC forbidden band. By introducing the buffer layer that destroys the periodicity of the 1D PhC, it is possible to excite the BSW which only appears within the band gap (see the pink line). Compared with TE mode, TM mode has a much narrower band gap due to the Brewster effect; in addition, the frequency range for observing BSW in the same structure is quite limited (see the blue line in
Figure 2b).
Figure 2c shows that for the TE case, when the refractive index of the background
n0 changes, the
kx of the excited BSW is varied. According to Equations (2) and (3), the
kx is directly connected with the incident angle and the operation frequency (please refer it to the dispersion diagram shown in
Figure 2a), by measuring the incident angle or operation frequency, the aim of sensing
n0 may be achieved.
When
kx and frequency satisfy certain condition (please refer it to the pink line shown in
Figure 2a), it is possible to produce BSW. The electric field distribution of the BSW mode in the
y direction is shown in
Figure 3. It can be seen that the electric field is mainly concentrated in the buffer layer and it gradually decreases along +
y directions. The electric field distribution of TE mode should satisfy the following formula:
where,
kx,inc and
kz,inc are the
x and
z components of
kinc Λ
x and Λ
z are the grating period in the corresponding direction.
Sm,n(
y) is the Fourier coefficients, which can be computed by solving Maxwell’s equations in Fourier space [
20].
By tuning the concentration and (or) pressure of the measured gas, the dielectric constant can be varied. Such phenomenon may be described by the equation as given as [
8]:
where we select the HCN as the measured gas, note that
is the background dielectric constant.
is the angular frequency corresponding to the absorption peak.
MHz/hPa is the damping rate.
is the absorption strength factor, which is a function of pressure
P and concentration
C. Given that the frequency of the incident THz wave is 14THz, the refractive index of HCN is:
n0 = 1.000 + 2.14 × 10
−7i. By increasing pressure or concentration, the refractive index of the gas can be tuned; here we assume the real part of the
n0 changes from 1.000 to 1.003, apparent shifts in the spectra can be found in
Figure 4a,b, where
Sθ decreases with the increase of
n0, that is,
Sθ and FOM can reach 123°/RIU and 147.3 /RIU, respectively, which can come up with following formula:
The change is because the angle sensitivity is proportional to the reciprocal of the incident cosine value. FOM has the same but very slow change of
Sθ, which is shown in
Figure 4c. It can be seen that as
n0 increases by 0.01, FOM only decreases 4.1/RIU. Similarly, the phenomena can also be observed in
Figure 4b,d, when wave is incident at a fixed 60°. By observing the shift of frequency spectrum, the change of
Sf is not obvious in
Figure 4b and FOM increases with the increase of
n0, that is,
Sf and FOM can reach 4.7 THz/RIU and 301.3 /RIU, respectively. This sensor can simultaneously reach a high and stable sensitivity and FOM in both angular and frequency domain. It needs to be pointed out that, since the contents of the top layer and substrate are PC, it is easy to fabricate micro-porous structure at the interface of buffer layer [
21]. By doing so, more polar molecules may be adsorbed to the relatively large surface, which may enhance the sensitivity [
22]; similarly, as more un-polar molecules are filled in the porous space of PC surfaces, the effective index of refraction of the top layer can also be adjusted, which in turn help enhance the sensitivity.
At present, most researches on THz gas sensor focus on sensing the change of real part of the analyte, see References [
1,
2,
3,
7]. The structure based on BSW as proposed here may also be applied for detecting the minor change of the imaginary part of the refractive index of gas. In the area of atmospheric pollution and dangerous gas detection, gases such as N
2O, HCN, H
2 and CO
2 have obvious absorption peaks in the THz band. From Equation (8), it can be seen that, when the operation frequency of the sensor is consistent with the absorption frequency of gas molecules, in which ω = ω
0i, different concentration
C and (or) pressure
P will make the imaginary part of the refractive index of analyte be varied; by measuring the change of transmission in the spectra, the sensing to the imaginary part may also be achieved. Here we select the gas N
2O with the absorption peak appeared in 17.664 THz (which is equivalent to the 588.726 cm
−1) to detect the change of imaginary part of the N
2O gas [
23]. Since there is no reference that can show the relationship between
and concentration C (or pressure P), in the following we consider four cases where the imaginary part is 0, 0.00002, 0.0002 and 0.002 respectively. The transmission spectrum as simulated with finite element method is shown in
Figure 5. It can be found that the transmittance is sensitive to the minor change of the imaginary part of the refractive index (e.g., the imaginary part changes from 0 to 0.002, the transmission changes from nearly 0 to 0.15) and the maximum sensitivity can reach 664 (transmittance unit/ one unit of the imaginary part of the refractive index). Further optimization in terms of increasing the number of photonic crystal layers, replacing PVDF with another polymer film that shows higher refractive index may also increase the sensitivity. It turns out that our design is not only sensitive for sensing the real part of the refractive index of the gas but also be applied as monitoring the minor change of imaginary part of the refractive index of the analyte. Consider that hazardous gases such as N
2O, HCN, H
2 and CO
2 show “fingerprints” in the THz band, the proposed design may find certain applications in sensing the concentration or pressure of such gases.
We discussed the real-part-refractive-index sensing of HCN at 14 THz and the imaginary-part-refractive-index sensing of N
2O near the absorption spectrum. This proves the applicability of the proposed structure. In recent theoretical terahertz sensing articles, only a few have explored sensing applications for practical gases. Performance comparison between the proposed BSW based design and other reported THz gas sensors is summarized in
Table 1. References [
1,
2,
3] all discuss the sensing sensitivity of different refractive indices and only Reference [
7] discusses the concentration sensing of H
2. It is difficult to find the relationship between the corresponding dielectric constant and concentration (pressure) of the corresponding gas, so it is hard to do practical inquiries into specific gases. For the comparison of sensing performance, angular sensitivity and the value of
Sf/
f can be applied. The
Sθ has a good performance compared to References [
1,
2,
3]. However compared to Reference [
7] in terms of frequency sensitivity, the value of
Sf/
f can only reach 0.34 /RIU. It is noted that since the refractive index contrast between PC and PVDF is quite limited ~1, which is much smaller than the case in silicon and air (3.4:1) as reported in Reference [
7], the normalized sensitivity of frequency
Sf/
f in the proposed design is small. However, the sensitivity in angle as shown in the BSW-based THz sensor is relatively high; in addition, considering the availability for sensing the minor change of both real and imaginary parts of refractive index of gases, our design may find potential applications in sensing hazardous gases in the terahertz band.