Abdulkarem H. M. Almawgania,
Hussein A. Elsayedb,
Ahmed Mehaney*b,
T. A. Tahacd,
Ziyad Awadh Alrowailic,
Ghassan Ahmed Alie,
Walied Sabrab,
Sayed Asaduzzamanfg and
Ashour M. Ahmedbh
aElectrical Engineering Department, College of Engineering, Najran University, Najran, Kingdom of Saudi Arabia
bPhysics Department, Faculty of Science, Beni-Suef University, Beni-Suef, 62512, Egypt. E-mail: ahmed011236@science.bsu.edu.eg; ahalmawgani@nu.edu.sa
cPhysics Department, College of Science, Jouf University, P.O. Box 2014, Sakaka, Saudi Arabia
dPhysics and Engineering Mathematics Department, Faculty of Electronic Engineering, Menoufia University, Menouf 32952, Egypt
eInformation Systems Department, College of Computer Sciences and Information Systems, Najran University, Najran, Saudi Arabia
fDepartment of CSE, Rangamati Science and Technology University, Rangamati, Bangladesh
gDepartment of Computer Science and Engineering, Daffodil International University, Dhaka, Bangladesh
hPhysics Department, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11623, Saudi Arabia
First published on 27th February 2023
In this research, we have a theoretical simple and highly sensitive sodium chloride (NaCl) sensor based on the excitation of Tamm plasmon resonance through a one-dimensional photonic crystal structure. The configuration of the proposed design was, [prism/gold (Au)/water cavity/silicon (Si)/calcium fluoride (CaF2)10/glass substrate]. The estimations are mainly investigated based on both the optical properties of the constituent materials and the transfer matrix method as well. The suggested sensor is designed for monitoring the salinity of water by detecting the concentration of NaCl solution through near-infrared (IR) wavelengths. The reflectance numerical analysis showed the Tamm plasmon resonance. As the water cavity is filled with NaCl of concentrations ranging from 0 g l−1 to 60 g l−1, Tamm resonance is shifted towards longer wavelengths. Furthermore, the suggested sensor provides a relatively high performance compared to its photonic crystal counterparts and photonic crystal fiber designs. Meanwhile, the sensitivity and detection limit of the suggested sensor could reach the values of 24700 nm per RIU (0.576 nm (g l)−1) and 0.217 g l−1, respectively. Therefore, the suggested design could be of interest as a promising platform for sensing and monitoring NaCl concentrations and water salinity as well.
Various techniques are used to determine the salinity of the water where it can be measured by some methods such as physical and chemical methods.5 Researchers used refractivity phenomena, density, and conductivity for salinity measurements in the case of using physical methods. While in chemical methods, directly determine the concentration of both NaCl and magnesium (Mg) dissolved in seawater.5 It is worth mentioning that the physical methods are more convenient and faster, but they are affected by electrical interference due to their association with electrical measurements.5,6 From this point, the optical techniques are used instead to determine and monitor the water salinity.
From the optical techniques, the most popular one is the photonic crystal periodic structure.7–12 Photonic crystals (PCs) have the possibility to control the electromagnetic waves propagation with different frequencies. One of the most important features which exhibited by the PC designs is the photonic band gap (PBG). Within this feature, specific frequencies of the electromagnetic waves cannot propagate through the PC design. By breaking the periodicity of PC designs, the confinement of light can be achieved where some resonant modes have appeared in the resulted optical spectra.
In the literature, the salinity level in salt solutions are determined using the photonic crystal optical fiber.13–17 Besides, many reports devoted attention towards the usage of the defective PC designs for detecting the salinity of seawater. For instance, in ref. 18, the authors reported about the salinity sensing based on a defective PC design. In their work, the transmittance spectrum is presented based on using a one dimensional (1D)-PC composed of fused silica and seawater. Also, the authors in ref. 19 have reported 1D-defective PC design for simultaneous detection of both seawater salinity and temperature in the vicinity of the 2 × 2 transfer matrix method (TMM).
Furthermore, the physical Tamm plasmon (TP) resonance phenomena was demonstrated in PC designs for many sensing purposes. Basically, Tamm resonance results from the interaction of the incident electromagnetic waves at the interface between a metal and Bragg mirror.20,21 There have been many previous reports about the demonstration of the TP resonance feature in PC designs such as in ref. 20–23. TP resonance showed high sensitivity and quality factor in the optical sensing techniques of PC structures.24–28 For the above-mentioned reasons, we try in this research to make a TP resonance/1D PC sensor to detect the concentrations of salt solutions in water. It is known that the salt level in the saltwater can be specified and indicated by the level of NaCl which, in turn, indicates the most predominant substance in salt water. Accordingly, in our study, we theoretically present the principle of the salinity sensing of NaCl solution based on the TP resonance feature. The current optical salinity sensor is constructed on a defective PC design. Various concentrations of NaCl samples will be injected into the respective cavity in the proposed PC design. Through simulation, we can obtain a resonance (Tamm resonance) peak, which has the tunability feature by adjusting the concentration of the NaCl samples. Finally, the performance of current salinity sensor is presented by investigating some important factors such as the sensitivity and the quality factor.
Fig. 1 Three dimensions representation of the suggested NaCl sensor that configured as, [prism/Au/water cavity/(Si/CaF2)N/glass substrate]. |
The NaCl solution samples will be injected into a cavity inside the proposed PC design where it is set between a metallic layer of gold (Au) and Si/CaF2 multilayers. The metallic layer located on the design top surface is necessary for appearing the Tamm resonance in the resulted PBG as we mentioned previously in the introduction. In the front, a prism receives the incident electromagnetic waves. Now, the whole proposed salinity sensor will have a configuration, [prism/Au/water cavity/(Si/CaF2)N/glass substrate], where N describes the structure's periodicity.
By using the well-known TMM, the optical parameters of the current salinity sensor are estimated as follows.29–31 According to this method, the electric and magnetic components in the Z-direction for the TE mode take the following form:
Ej(z) = Gjexp(−ikjz) + Hjexp(ikjz) = Ey+ + Ey− | (1) |
(2) |
Then, the electric and magnetic components can be described by the matrix formalization as follows.
(3) |
(4) |
This can be done when the incident waves interact with a distinct layer having a thickness of dj = Z2 − Z1. According to eqn (3) and (4) we find that,
(5) |
In these equations, fj gives the relation between the electric and magnetic fields components along the two sides of the layer j, θj describes the incident angle, and Kj represents the wave vector. Then, the total characteristic matrix for all the structure layers that interact with the incident waves can be described as follows.
(6) |
After that, the estimated reflectance of the incident electromagnetic waves is given as follows,
(7) |
In our calculations, the Drude model formulates the permittivity of Au layer as follows.32,33
(8) |
(9) |
For the fabrication of the proposed sensor, an Au layer was deposited on the prism by using the RF sputtering method.37 Next, a polyimide (PI) polymer layer was prepared by using the chemical vapor deposition (CVD) technique on the prism/Au.38 Then, Si/CaF2 multilayer films were grown on prism/Au/PI. The periodic Si/CaF2 multilayer structures can be deposited at room temperature by using the molecular beam epitaxy (MBE).39,40 The Si/CaF2 multilayer structures were studied in many previous experimental works.41,42 The whole multilayer structure is now configured as prism/Au/PI/(Si/CaF2)10/air. The cavity in a stack of multilayers structure can be formed by using chemical etching for the PI polymer layer. The final structure is prism/Au/cavity/(Si/CaF2)10/air.
Then, Fig. 3 displays the reflectance spectrum of the PC that configured as, [prism/(Si/CaF2)10/substrate] at normal incidence. Here, the thickness of Si and CaF2 layers are proposed as 100 nm and 900 nm. Furthermore, the substrate is designed from glass having a refractive index of 1.52. Fig. 3 demonstrates the formation of the PBG with a width of 425 nm where their edges from the left side to the right one are located at 1456 nm and 1881 nm, respectively. However, a PBG with this width has a little interest specifically TP resonance modes are highly sensitive to variation of refractive indices due to NaCl concentration variations.
To simulate the detection of NaCl, we have introduced a 20 nm thin metallic layer and a cavity layer of water with a thickness of 5 μm on the PC top surface (see Fig. 1). Hence, the last visualization of our design was, [prism/Au/water cavity/(Si/CaF2)10/substrate]. In this regard, the Au layer refractive index could be investigated from eqn (8) in which ωP = 1.323 × 1016 and γ = 1.26 × 1014.45 Fig. 3B shows the appearance of three sharp dips with intensities 4.3%, 1.07%, and 0.04% in the reflectance spectrum at wavelength positions of 1483.33 nm, 1651.77 nm, and 1859.34 nm, respectively. These sharp reflectance dips could be due to TP resonance. Such resonance is introduced through the reflectance spectrum of multilayer structures or PCs because the free electrons' oscillations are quantized at the interface between the thin metallic layer and cavity.46 Physically, the sharp TP resonant modes (the three sharp dips) with a high resonance in the water cavity is due to the coupling between two fields. One field is a narrow mode and this is due to the plasmon resonance in the metallic layer and the other one through the evanescent fields present in the Bragg mirror.47 To sum up, it means that the narrow TP modes that appear in the liquid cavity come from the hybridization between the narrow plasmon mode of the metal layer and evanescent modes of the dielectric layers, producing slow-guided modes (three narrow dips) of low group velocity.48,49 Therefore, these resonance modes could be sensitive against the refractive index variations with concentration. Therefore, the detection of NaCl solution is correlated with the shift of these dips with the variations in the concentration of the used analyte. However, the presence of only one resonance dip through the reflectance spectrum could offer more simplicity besides overcoming any conflict during the detection procedure unlike the case of two or three dips.
Thus, we aimed to investigate only one TP resonance through the reflectance spectrum. In this context, the role of the angle of incidence could be crucial. Thus, we consider TE mode of polarization at the oblique incidence in which the incidence angle is equivalent to 50°. Fig. 4A shows the intensity of the reflection for prism/(Si/CaF2)10/substrate structure, which is about 100% within the spectrum range at an angle of 50°. This is because the angle is now close to the critical angle and a total internal reflection (TIR) occurred. Physically, the optical path length (OPL) of a light ray moving through a medium is given by the product of the refractive index of the medium and the geometric length of the ray. The OPL is a very important in determining the phase shift of the light and governs diffraction and interference of light as it propagates. Based on Bragg's law, the bandgap of the PC structure could be due to the destructive interference of multiple reflections of light at each interface between the layers of the crystal. As the angle incident angle increase from 0 to 50°, the OPL changes, and hence the number of constructive interference increase and the number of peaks that can be localized in the cavity layer increases as well. This means that the number of localized resonance states increases. At an angle of 50°, only one dip starts to appear in the spectrum, and this is logic because the angle now is close to the critical angle. Then, all dips disappeared, and the intensity of the reflection is about 100% due to the total reflection.
On the other hand, for [prism/Au/water cavity/(Si/CaF2)10/substrate] designed structure, only one dip started to appear in the spectrum, as shown in Fig. 4B. When there is a light beam incident on the prism with an incident angle greater than the critical angle, a very small portion of light can penetrate the Au-prism interface and generate an evanescent wave [https://scitation.aip.org/content/contributor/AU0025156].50 This evanescent wave energized the surface plasmon wave (SPW). Therefore, the reflectivity of the incident beam will drop to zero, thereby leading to a single TP resonance.
The choice of such condition is coming after the optimization of the angle of incidence. Therefore, at 50°, a single TP resonance appears through the wavelengths of interest as shown in Fig. 4. From this figure, the dip showed a reflectivity of 1.82% where it is located at λTP = 1681.55 nm. Moreover, the change in optical path length produces a change in both the number and position of the reflectance dips with the increase of incidence angle. As a result, this dip will be the key axis during the detection procedure of our study.
The thickness of the Au layer has a strong effect on the TP resonance mode and sensor performance. The optimum value of the Au layer is shown in Fig. 5. Fig. 5A displays that the resonant TP dip is shifted downwards the lower wavelengths with increasing the Au layer's thickness from 5 nm to 40 nm. The intensity of dip reflectivity of the TP mode has a low value at the Au layer with a thickness of 20 nm, as shown in Fig. 5B. This condition corresponds to an excellent coupling of the TP resonance by reducing the electron energy loss and the narrow width of the dip. Therefore, the optimum value of the Au layer's thickness is 20 nm, which provides a relatively high sensor performance.
Fig. 5 The effect of change the thickness of the Au layer on the(A) reflection spectrum (B) intensity of reflection at of TP resonance. |
Meanwhile, there a large difference between the defective PC designs such as the results presented in ref. 18 and our PC design. Notably, Chaves et al. discuss the role of 1D PCs in the detection and monitoring of water salinity based on the inclusion of a defect layer in the middle of the designed PCs.18 By introducing the defect layer, a resonant mode is introduced through the PBG due to the broken periodicity of the designed PC structure. Then, its position can be shifted with the variation of the saline solution concentration as a result of the changes in the OPL of the incident radiation with the refractive index of the defect layer. In contrast, our study is essentially based on the localization of TP resonance due to the deposition of a thin metallic layer on the top surface of a 1D PCs. In fact, we believe that the suggested design could be of a potential interest compared to the conventional 1D defective PCs due to its relatively high performance. Then, the inclusion of a cavity layer between the thin metallic layer and the 1D PCs does not provide any chance for introducing a resonant and its role is limited for the inclusion of the analyte. In particular, this layer is classified as a cap layer that could provide some effects on the path length of incident radiation.
Thus, the position of TP resonance is expected to change significantly by immersing different concentrations of the NaCl through the cavity as shown in Fig. 6. In this regard, the proposed structure has been designed for detecting the NaCl concentration in water at room temperature (25 °C). Notably, within a small range of temperatures from 20 to 60 °C, the thermo-optic (TO) coefficients for Si, CaF2, and Au are very small. Moreover, some of the experimental verifications of PC structures on Si are insensitive to the little temperature changes especially that is close to room temperature.51 This means that the thermal effect on the sensor's behavior can be negligible in the temperature regime of interest for the working device.
Thus, at 10 g l−1 of NaCl, the resonance dip (TP resonance) changed from 1681.55 to 1685.75 nm. A small increase in the refractive index of the cavity with increasing the NaCl concentration produces the shift in the wavelength position of TP resonance. Such an effect could lead to the change in the optical path length of the incident radiation. For more increase in the concentration of NaCl to 20 g l−1, 30 g l−1, 40 g l−1, 50 g l−1 and 60 g l−1, the resonance dip is shifted upwards to wavelength positions of 1690.55 nm, 1696.05 nm, 1702.1 nm, 1708.84 nm, and 1716.13 nm, respectively, as shown in Fig. 6 and the color map in Fig. 7. Moreover, the reflectance values of these dips is slightly increasing with the increase of NaCl concentration. Meanwhile, the reflectivity increases from 1.82% to 2.42% as the concentration increases from 0 g l−1 to 60 g l−1, respectively. Here, our design is characterized by a relative high sensitivity to the concentration variations despite the very little increase in the refractive index of the cavity with the increments of NaCl concentrations as shown in Fig. 8. Such a result could make our design very promising in the detection and sensing of NaCl concentration and salinity as well compared to other physical methods and techniques.
Fig. 7 A color map for the reflectivity of TP resonance at different concentrations of NaCl solution. |
After that, the TP resonance position regarding NaCl concentration was plotted in Fig. 9. The figure indicates that the increase in the NaCl concentration leads to a significant shift in the position of TP resonance towards the higher wavelengths. Accordingly, the TP resonance position increases at increasing the NaCl concentration. This response can be described by fitting the numerical results as depicted in Fig. 9. The quadratic fitting leads to a relatively high coincidence with the investigated numerical results in Fig. 9. Such fitting could be mathematically described according to the following equation:
TP position (nm) = 0.003125C2 + 0.38918C + 1681.55 | (10) |
Fig. 9 The resulted relation between the spectral position of TP resonance and the NaCl concentration according to a quadratic fitting. |
This equation reveals that the concentration of NaCl solution can be simply obtained based on the TP resonance position. In other words, at a distinct wavelength of TP resonance, eqn (10) can be solved to obtain the exact value of NaCl concentration.
Then, the performance of the designed structure is mandatory to determine the validity of such structure as an effective sensor for the monitoring and detecting of NaCl concentration and salinity as well. Thus, there are many different parameters that can highlight the sensor performance such as sensitivity (S), figure of merit (FM), detection limit (DL), signal to noise ratio (SNR) and quality factor (QF). These parameters are widely described and analyzed analytically in many literatures.42–46 The mathematical formulas of these parameters can be calculated from the following relations:52–56
(11) |
(12) |
(13) |
(14) |
Then, the performance of the designed structure is mandatory to determine the validity of such structure as an effective sensor for the monitoring and detecting of NaCl concentration and salinity as well. In this regard, Fig. 10 describes the sensitivity of our design against the concentration of NaCl solution. The value of sensitivity is equivalent to the TP resonance shift due to the change in NaCl concentration. Fig. 10 shows a linear increase of sensitivity with the NaCl concentration. Here, the value of S can reach 0.576 nm (g l−1) as the concentration of the NaCl solution is equivalent to 60 g l−1. Such value is seemed to be very promising through the detection procedure especially in the monitoring of low concentrations. In addition to that, it could give a positive indication towards other parameters like the detection limit and quality factor as well. Therefore, we have plotted in Fig. 11 the change in the quality factor against the concentration variations. As you see from eqn (11)–(14), the values of sensitivity and QF in Fig. 10 and 11, respectively, increases due to the changes in the spectral position of TP resonance. Therefore, the displacement of the TP mode (ΔλTP) increases with NaCl concentration, which, in turn, makes the sensitivity and QF values increase with NaCl concentration as well. Here, QF describes the ratio between the wavelength of TP resonance relative to the full width at half maximum (FWHM). In this context, the presence of sharp dips with small values of FWHM leads to high values of QF. Fig. 11 indicates that the values of QF have a nonlinear increase with the increments of NaCl concentration. Such response is mainly connected with the values of the spectral position of the TP resonance and FWHM as well at a specified concentration of the NaCl solution. In particular, the changes in the concentration values leads to a significant effect on the refractive index of the cavity layer. Meanwhile, a different response of the OPL of the incident waves in the vicinity of the FWHM and TP spectral position is expected. In this regard, the QF has the value of 686.452 at NaCl concentration of 60 g l−1. Then, Fig. 12 shows the detection limit of the proposed sensor versus the variations of NaCl concentration. Such parameter has an important role on the sensors' performance. DL is strongly depending on the position of TP resonance, S, and QF. The Fig. 12 clarifies that the values of DL are approximately decreasing linearly with the concentration. DL provides a minimum value of 0.217 g l−1 at concentration = 60 g l−1. In addition, DL records relatively low values as a result of the high values of S and QF as well. In particular, low values of DL support the high performance of sensors.
Finally, to the best of our knowledge, the proposed design could be of significant interest as a promising platform for the detection and monitoring of NaCl concentration and salinity as well. Specifically, our designed sensor provides a relatively high performance compared to its counterparts in photonic and optical techniques as listed in Table 1. As illustrated in Table 1, the proposed sensor in this study has achieved a good reasonable performance compared to other designs in the literature according to the obtained performance parameter values. Based on Table 1, the introduced design has a high sensitivity value of 24700 (nm/RIU), which is considered higher than the values of previous designs. Also, the proposed sensor demonstrated a DL value of 5.91 × 10−6, which is better than that of previous designs. In addition, the proposed sensor showed a QF value of 686.5, which is an acceptable value compared to other designs.
The designed structure | Sensitivity (nm/RIU) | QF | DL (RIU) | Reference |
---|---|---|---|---|
Photonic crystal biosensor | 1118 | Not mentioned | 10−3 | 57 |
Photonic crystal cavity | 3300 | 103 | 10−5 | 58 |
1D-porous silicon photonic crystals | 6770 | 662 | Not mentioned | 59 |
Photonic crystal cavity and fiber loop ring-down technique | 450 | Not mentioned | ∼1.6 × 10−4 | 60 |
Photonic crystal containing graphene | 1179 | Not mentioned | ∼2.2 × 10−5 | 61 |
Photonic crystal covered with a perforated gold film | 17 | 3 × 104 | Not mentioned | 62 |
Cavity photonic crystal | 777 | 2576 | Not mentioned | 63 |
Photonic crystal cavities | 656 | 2719 | Not mentioned | 64 |
[Prism/Au/water cavity/(Si/CaF2)10/substrate] | 24700 | 686.5 | 5.91 × 10−6 | Our design |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ra00308f |
This journal is © The Royal Society of Chemistry 2023 |