Based on the above parameters,
Figure 2 plots the corresponding transmittance spectrum of the structure when the blood plasma concentration changes from 0 g/L to 50 g/L in steps of 10 g/L. It can be seen from
Figure 2 that there is a complete forbidden band in the wavelength range of 1500–1600 nm. The increasing blood plasma concentration causes the whole forbidden band to slightly shift to regions of longer waves while the forbidden band width basically remained unchanged. Consequently, the central position of the forbidden band is still close to the central wavelength
λ0 = 1550 nm. For the whole PT-symmetry structure, the plasma blood sample layer is equivalent to the defect cavity of the structure, so that the corresponding defect mode (transmittance peak) will come into being in the transmittance band gap of the structure. Surprisingly, we have noticed in this figure that, different from the defect mode of the general photonic crystal structure, the transmittance of the defect mode in the forbidden band is greater than 1. Meanwhile, when the blood plasma concentration changes from 0 g/L to 50 g/L, the position of the transmittance peak in the forbidden band changes from 1529 nm to 1550 nm, and the amplification effect of the structure on the working light wave is gradually increases with the red shift of the transmittance peak. This is predominantly dependent on the strong coupling resonance effect between the PT-symmetry unit and the microcavity, the structure converts pumped energy from the outside into electromagnetic energy of incident light, and thus, the transmittance of defect mode is amplified. In combination with Equation (1), when the working light wavelength is equal to the central wavelength, the refractive index of AB layer accords strictly with the PT-symmetry condition, and the coupling resonance effect of the structure is the strongest. Therefore, the amplification effect of the structure on the working light wave is most obvious at the center wavelength of the forbidden band.
For improving the amplification effect of the structure on the working light wave and enhancing its detection performance in the whole blood plasma concentration range, it is essential to accurately optimize and calibrate the geometric parameters of the structure and material properties of the PT-symmetry dielectric layer. The following analysis focuses on the optimization process of parameters such as unit cell number N of the PT-symmetry structure and the thickness of the plasma blood sample layer dD, as well as the influence of macroscopic Lorentz oscillation intensity α in the PT-symmetry dielectric layer on the detection performance.
3.1. Optimization of Unit Cell Number N of PT-Symmetry Structure and Sample Layer Thickness dD
According to the above analysis, the amplification effect of the structure on the incident light is caused by the PT-symmetry unit. Hence, we need to investigate the influence of the unit cell number
N on the amplification effect of the defect mode. Initially, maintaining the other structural parameters as fixed, in
Figure 3, the transmittance spectrum at blood plasma concentration of 50 g/L for the structure (ACB)
ND(BCA)
N is shown with the four unit cell values given:
N = 5, 6, 7 and 8. Since the unit cell number
N has a significant difference in the influence on the transmittance, here, logarithmic transformation, with dB as the unit, is taken for the transmittance during calculation. The transformation relation can be written as:
T(dB) = 10lg
T.
As shown in
Figure 3, when the blood plasma concentration remains unchanged, the inhibition of the structure on the incident light in the forbidden band will gradually fortify with the increase in
N, while the transmittance of the defect mode nonmonotonous change. When
N = 5, 6, 7 and 8, the transmittance is 3.532 dB, 19.86 dB, −2.962 dB and −14.69 dB, respectively. It can be seen that the transmittance of the defect mode is largest when
N = 6. This indicates that the structure has the strongest amplification effect on the transmittance of defect mode in the center of the forbidden gap. In order to further study the influence of unit cell number
N of PT-symmetry structure on the transmittance of defect mode in the whole blood plasma concentration range, we calculated the transmittance of defect mode of the structure (ACB)
ND(BCA)
N in different unit cells with the blood plasma concentration changing from 0 g/L to 50 g/L (taking 10 g/L as a step). The calculation results are shown in
Figure 4.
As shown in
Figure 4, with the increase in blood plasma concentration, the position of defect mode is undergoing an obvious red shift. When
N = 5 and 6, the transmittance of defect mode increases with the increase in plasma concentration. On the contrary, when
N = 7 and 8, the transmittance of defect mode decreased with the increase in blood plasma concentration. Changing the unit cell number
N does not affect the position of defect mode in the case of the same blood plasma concentration, but significantly affects the amplification effect of structure on the transmittance of the defect mode. As the unit cell number
N increases from 5 to 8, the transmittance of defect mode firstly increases and then decreases, and reaches the maximum at the same blood plasma concentration when
N = 6. This is because the group velocity of the working light wave in the gain layer (B layer) is lower than that in the loss layer (B layer) when
N = 6. The photon travels in the gain layer for a relatively longer time, which leaves the energy of the pumped light source more converted into the electromagnetic energy of the working light through the gain layer; accordingly, the structure exhibits an amplification effect. In order to reasonably take advantage of this amplification effect to improve the detection performance of the structure, we select
N = 6 as the final optimization result.
In order to ensure the amplification effect of the structure on the transmittance is a suitable match to the whole blood plasma concentration range, so as to avoid affecting the detection performance of the structure due to the small transmittance in the low concentration range, the thickness of the sample layer has been exactly adjusted. When the blood plasma concentration changes from 0 g/L to 50 g/L, the frequency shift range of the transmittance peak can be symmetric with the central wavelength of the forbidden band. Considering that the transmittance peak position is observed in a linear relation with blood plasma concentration. Initially,
Cp = 25 g/L will be set, the thickness increment ∆
dD of the sample layer increases successively from 0 nm with an interval of 10 nm, and the transmittance spectra of structures under different ∆
dD are calculated.
Figure 5 reveals the variation law of the defect mode position and transmittance along with ∆
dD according to the calculation.
As shown in
Figure 5, when ∆
dD increased from 0 nm to 63 nm, the position of defect mode moves towards the long-wave direction and the transmittance increases gradually, and the transmittance peak is located at the wavelength of 1550 nm when ∆
dD = 63 nm. Meanwhile, we can readily observe the transmittance of defect mode reaches the maximum when its position is at the strongest coupling point, which corresponds to the position shown by the dotted line in
Figure 5. For ∆
dD greater than 63 nm, we detect that the transmittance of defect mode decreases accordingly. Therefore, we take the sample increment ∆
dD = 63 nm as the final optimization result and determine the sample layer thickness
dD = 1999 nm.
3.2. Detection Mechanism of Blood Plasma Concentration and Analysis of the Structural Performance
Based on the above optimized structural parameters, the transmittance spectrum of the structure at different blood plasma concentration has been calculated. We chose the blood plasma concentration point from 0 g/L to 50 g/L and separated by 10 g/L. The calculation results are shown in
Figure 6a. In
Figure 6a, we observe the obvious red shift of the transmittance peak motivated by defect mode with the increase in blood plasma concentration. It moves from 1539 nm corresponding, to 0 g/L, and to 1561 nm, corresponding to 50 g/L. Meanwhile, a change in the magnitude of the transmittance peak varies depending on blood plasma concentration. When blood plasma concentration is located in the center of the interval, the PT-symmetry structure results in the significant amplification of the transmittance of defect mode. When the value of the blood plasma concentration moved from the center of the interval to the high or low concentration interval, the transmittance decreases gradually, and the red shift interval of the whole transmittance peak is distributed symmetrically around the strong coupling point (
λ0 = 1550 nm). In order to facilitate the detection of sample concentration by the position of defect mode, the fitting curve between sample concentration and the position of defect mode is presented in
Figure 6. Referring to the calculation results in
Figure 6b, the fitting function can be obtained as:
where
Cp and
λ are the values of applied blood plasma concentration and transmittance peak position. From
Figure 6b, it can be seen that there is a good linear relationship between sample concentration and defect mode position, with a high correlation coefficient up to 0.9999.
The sensitivity of the structure can be expressed as the ratio of the wavelength shift amount ∆
λ to the sample concentration variation amount ∆
Cp or the ratio of the wavelength shift amount ∆
λ to the sample refractive index variation amount ∆
n, which is defined using the expression [
24]:
Since there is a linear relationship between the sample concentration and shift of the defect mode, the slope of Equation (7)’s fitting line can be regarded as the detection sensitivity of the structure in the whole blood plasma concentration range according to the definition of sensitivity. The calculation result shows that the detection sensitivity of the structure can reach 0.4409 nm/(g/L) (or 227.05 nm/RIU). Therefore, based on the good linear relationship between sample concentration and defect mode position and the high detection sensitivity of the structure, the blood plasma concentration value of the corresponding sample layer can be calculated by detecting the transmittance peak position with Equation (7), so as to realize the detection of the physicochemical properties of plasma blood samples.
In the actual detection, the performance and the efficiency of the sensor is often evaluated by comparing and analyzing multiple indexes, such as the quality factor (Q), the figure of merit (FOM) and detection limit (DL) of the structure. Therefore, according to the calculation results, multiple indexes should be considered to optimize the structural parameters. The quality factor Q is inversely proportional to loss of energy value and directly proportional to stored energy inside the cavity, the greater the Q value is, the more obvious the localization of defect mode is; therefore, a sensor has high Q values that demonstrate that the proposed sensor has very high accuracy because the defect mode has a highly narrow bandwidth [
13]. The Q is expressed by the ratio of transmittance peak position
λ to the full width at half maximum (FWHM) and can be calculated according to the following equation [
25]:
The ratio between the sensitivity (
S) and the FWHM is referred to as the figure of merit (FOM), which is obtained by:
FOM is proportional to the product of
S and Q. The larger the FOM value is, the better the detection performance of the sensor is. The DL is inversely proportional to
S and Q, according to [
26]:
The detection limit indicates the smallest detectable refractive index change, it can be seen from Equations (9) and (10) that DL is inversely proportional to FOM. Therefore, the larger the FOM value is, the smaller the detection limit (DL) of the structure is. Meanwhile, combining Equations (9) and (10), (11) can be expressed as:
When the transmission peak produces a small displacement, the FWHM of the transmission peak can be considered to remain unchanged. Therefore, it can be seen from Equation (12) that DL∙S is a fixed value, detection limit (DL) and sensitivity (S) show a trade-off relationship. The smaller the DL of the sensor, the smaller the refractive index change of the detectable blood sample, and the higher the sensitivity S. The simulation results show that adjusting the value of macroscopic Lorentz oscillation intensity α will affect the amplification effect of the structure, change the half-width of the defect mode transmittance peak and then have a great influence on the structure detection performance. The following section concentrates on the influence of α value on the quality factor, figure of merit and resolution of the structure.
Other structural parameters remain constant, and the blood plasma concentration is set as
Cp = 25 g/L to calculate the transmittance of the structure with α value increasing from 2.3 × 10
−4 to 2.9 × 10
−4 in steps of 0.1 × 10
−4. The 3D plot of the transmittance spectrum is shown in
Figure 7, where the position of transmittance peak remains unchanged, while the transmittance shows a tendency to increase first and then decrease with the increscent of
α values. When
α values range from 2.5 × 10
−4 to 2.7 × 10
−4, the amplification effect of structure is most obvious at the same blood plasma concentration. Since adjusting the value of α exclusively changes the gain and loss level of the AB layer and does not affect the position of the transmittance peak in the forbidden band, the value of α does not affect the concentration detection using the fitting result of Equation (7).
It can be seen from Equations (9) and (10) that the full width at half maximum FWHM of the transmittance peak affects the quality factor and the figure of merit of the structure, while the FWHM is related to the amplification effect of the structure on the working light, and its value also affects the resolution of the structure. In order to better analyze the influence of
α value on the overall detection performance when the blood plasma concentration is detected within the whole range, we calculated the transmittance spectrum at different blood plasma concentrations. When the blood plasma concentrations increase from 0 g/L to 50 g/L, based on
Figure 7, we take the average FWHM of the defect mode transmittance peak as a reference factor for analysis. In
Figure 8a, we present the variation curve of the average FWHM of the defect mode with
α. Considering that the α value is limited by the actual doping process, the
α is accurate to 10
−6 in the calculation.
As can be seen from
Figure 8a, the decrease in the average FWHM value of the transmittance peak occurred when
α increases from 2.3 × 10
−4 to 2.55 × 10
−4, and the minimum value of the average FWHM reaches 0.0238 nm when
α = 2.55 × 10
−4. As
α continues to increase from 2.55 × 10
−4 to 2.9 × 10
−4, the average FWHM value increases gradually. In order to further analyze the performance of the structure, we calculate the quality factor and the figure of merit FOM of the structure at different blood plasma concentrations with the same
α value, and take the average quality factor and the figure of merit as the detection performance of the structure under the same
α value.
Figure 8b represents the variation curve of the average quality factor and average figure of merit of the structure with
α. It can be seen that when
α = 2.55 × 10
−4, the average quality factor and the average figure of merit of the structure are the maximum, reaching 78,564 and 11,515 RIU
−1, respectively. Therefore, the optimal value of
α is 2.55 × 10
−4.
For expressing the detection performance parameters of PT-symmetry microcavity structure at
α = 2.55 × 10
−4 clearly, the data in
Figure 8 are extracted.
Table 1 shows the full width at half maximum FWHM, quality factor Q and figure of merit FOM at different blood plasma concentrations when
α = 2.55 × 10
−4; at the same time, Equation (11) is used to calculate the detection limit of the structure at different blood plasma concentrations, the calculation results are listed in
Table 1.We have noticed that when the blood plasma concentration is located in the center of the whole interval, the structure shows a compelling detection performance. By contrast, the detection performance decreases in the high and low concentration range. For the whole blood plasma concentration range, the minimum Q, the minimum FOM and the maximum DL of the structure also reaches 22,300, 3244 RIU
−1 and 12.6 × 10
−6 RIU, respectively.
Resolution is another indicator of the structural performance, which can be described as the minimal detectable blood plasma concentration; the smaller the value, the better the detection performance of the structure. In this paper, the position of the transmittance peak is used for concentration detection, so the resolution can be derived according to the change in blood plasma concentration and the displacement of the transmittance peak. The shift in the defect mode position occurs when the blood plasma concentration changes. For practical detection, the displacement interval of the transmittance peak should be greater than the full width at half maximum of the transmittance peak. As can be seen in the displacement of the transmittance peak of
Figure 9a, when
Cp = 0 g/L and 0.05 g/L, the interval between the two peaks is 0.0212 nm, which is exactly equal to the full width at half maximum, so the resolution of the structure at 0 g/L is 0.05 g/L. This method can be exploited to measure the structure resolution at different blood plasma concentrations. Since the mean full width at half maximum reaches the minimum value at
α = 2.55 × 10
−4, and the blood plasma concentration has a good linear relationship with the position of transmittance peak, therefore, the structure had the best resolution in the whole concentration range. According to the relationship between blood plasma concentration and transmittance peak, we plot the resolution of the structure in different blood plasma concentration in
Figure 9b. The values are measured based on the thickness of sample layer
dD = 1999 nm,
α = 2.55 × 10
−4 and
N = 6. It can be seen that the resolution of the structure varies from 0.006 g/L to 0.057 g/L in the whole range of blood plasma concentrations. When the blood plasma concentration is between 15 g/L and 25 g/L, the resolution value is relatively small, and the structure detection performance is more remarkable. In the high and low concentration ranges, the resolution is basically maintained at 0.05 g/L, and the average resolution is up to 0.038 g/L for the whole blood plasma concentration range.
Regarding the sensor structure proposed in the literature [
20], the maximum detection sensitivity is only 45.06 nm/RIU in the whole blood plasma concentration range, and the full width at half maximum FWHM of the transmittance peak is wider than the frequency shift of the transmittance peak position in the whole blood plasma concentration range, which results in the quality factor, figure of merit and resolution of the structure extremely small. As for the sensor for serum tissue detection proposed in literature [
27], although the sensitivity of the structure can reach 6857.89 nm/RIU, the transmittance of the defect mode is lower than 0.4, and the FWHM of the transmittance peak is relatively large, which not only affects the practical application of the structure, but also is not conducive to the signal processing of the detection equipment.
Table 2 shows the optimized structural parameters and the performance of the structure in blood plasma concentration detection. By adjusting and optimizing the plasma blood sample layer
dD =1999 nm, the cycle of PT-symmetry unit
N = 6, and the macroscopic Lorentz oscillation intensity
α = 2.55 × 10
−4, under these conditions, the structure is used to detect the plasma blood samples with blood plasma concentration from 0 g/L to 50 g/L. Its sensitivity, average quality factor, average figure of merit, average detection limit and average resolution are 0.4409 nm/(g/L) (or 227.05 nm/RIU), 78,564, 11,515 RIU
−1, 5.1 × 10
−6 RIU and 0.038 g/L, respectively. Compared with other research on biomedical sensors, the detection performance of the structure is comprehensively analyzed in this paper, and each performance index is improved synchronously. The reported results indicate that the structure has potential application value in blood detection.