1. Introduction
In industry, sensors/transducers have applications that include healthcare and medicine, air quality, food safety, and fuel storage. A problem in the aviation industry is the contamination of fuel with foreign organisms such as yeast [
1,
2,
3]. The accumulation of bio-organisms, which feed on the carbon in fuels, reduces fuel stability, corrodes storage tanks, and degrades the functions of valves, pumps, and other mechanisms in a fuel system; all of these problems are termed fuel biofouling [
4,
5]. Portable transducers that are inexpensive, collect data for long time intervals, and generate data quickly after a binding event, can be used to monitor the status of jet fuel stored in tankers.
When detecting chemicals or organisms, it can be important to quantify the thickness and index of refraction of an adhered layer of analyte on the transducer surface, as well as the refractive index of the background/bulk media [
6,
7]. Quantifying these three variables requires multiparametric transducer input; in the case of a GMR sensor, these are three or more resonant peaks that shift due to these sensor variables [
6]. A GMR sensor produces resonant peaks when interrogated with white light. As the sensor variables change, the wavelengths of the resonant peaks shift, and these shifts can be correlated to the magnitude that the sensor variables have changed [
6,
8].
Past work has determined sensor variable values using a sensor’s sensitivity, and the magnitudes of resonant wavelength shifts are monitored as one sensor variable is fixed and the other sensor variables are solved for [
9]. This method usually requires that one biolayer variable is held constant as the other varies [
10,
11]. The sensitivities can be expressed as:
These expressions denote the biolayer index of refraction sensitivity {
S(
nbio)}, biolayer thickness sensitivity {
S(
dbio)}, and bulk index of refraction sensitivity {
S(
nbulk)}, respectively. In the above expressions, Δλ, Δ
nbio, Δ
dbio, and Δ
nbulk are the change in resonance wavelength, change in biolayer index of refraction, change in biolayer thickness, and change in bulk (background) index of refraction, respectively. Using sensitivity to determine biolayer sensor values falls short because the biolayer thickness (biolayer index of refraction) sensitivity is dependent on the value of the biolayer index of refraction (biolayer thickness). For instance, a greater biolayer thickness yields a greater absolute value of biolayer index of refraction sensitivity. This is because the evanescent tail of a resonant mode sees a larger change for a thicker layer. The analogue of this concept applies to the magnitude of biolayer thickness sensitivity: a greater biolayer index of refraction yields a greater absolute value of biolayer thickness sensitivity. Restricting one biolayer variable to a constant value, during biolayer growth, to calculate the other biolayer variable is a method to recon with the interdependence of the biolayer variables’ sensitivities; we avoid this restriction using our method of analysis. Our method expands on past work that used a GMR sensor with two resonant modes and simulation tools to create a lookup table [
6]. In the work performed by Magnusson et al., the lookup table is used to solve for the biolayer index of refraction and the bulk index of refraction, while biolayer thickness is held constant [
6].
There have been methods proposed to interpret transducer output and deconvolve the biolayer/adlayer thickness and biolayer/adlayer index of refraction [
12,
13,
14]. These works include using a surface plasmon resonance sensor to conduct two experiments to obtain two data sets, then interpreting the data to determine the dielectric constant and thickness of a dielectric layer [
15]. Another related work discusses the theory of using a combination GMR-SPR (guided mode resonance-surface plasmon resonance) sensor with three resonant modes; among the three resonant modes, the biolayer sensitivities differ by orders of magnitude and the bulk sensitivities also differ, to a lesser extent [
16]. Using the values of sensitivity, matrix methods are proposed to calculate the biolayer thickness, biolayer index of refraction, and bulk index of refraction [
16]. In addition to a sensor being multiparametric, such as those mentioned above, it is also of great importance that a sensor does not require the analyte to be tagged or altered for detection (label-free) [
17].
A bio-selective layer can be applied to the sensor surface to capture the desired analyte and reject other substances in the environment [
4,
6,
18,
19]. Thus, a GMR sensor can be tailored to detect an analyte of choice without the use of labels, as in label-free sensing. Multiparametric and label-free sensing is in high demand in industry, and it can be performed using a GMR sensor [
6,
20].
Other label-free sensors include surface plasmon resonance sensors, integrated interferometers, MEMS-based sensors, nano-sensors (rods and particles), Bragg grating sensors, photonic crystal-based sensors, ring-resonator sensors, ellipsometry, and grating coupled sensors [
6,
21,
22]. Sensor schemas that utilize labels include immunomagnetic separation, polymerase chain reaction, and standard immunoassay; these sensor types use luminescent, radioactive, absorptive, and fluorescent labels [
6]. Sensors that utilize labels require the extra step of altering the analyte for detection. In contrast, the GMR sensor surface is altered to selectively capture and detect the analyte of interest.
Surface plasmon resonance (SPR) sensors, which are most like GMR sensors, utilize a resonance effect at the interface between a dielectric and metal [
23,
24]. When a TM polarized electromagnetic wave reaches a dielectric and metal interface at a specific angle of incidence, the electromagnetic wave becomes evanescent at the interface while interacting with the free electrons in the metal [
25]. This phenomenon produces an absorption minimum in the spectrum that has high angular and spectral sensitivity [
26]. The linewidth of SPR sensors is large, and only a TM mode produces an SPR response. In contrast, GMR sensors have a smaller linewidth, preferred over a large one, and the resonance effect can be produced using both the TE and TM modes; this allows the monitoring of the changes in a greater number of sensor variables with mixed polarization states.
The GMR sensor produces resonant modes by diffracting incident broadband light into leaky waveguide modes, allowing standing waves to form in the sensor at specific wavelengths (frequencies) as eigenmodes [
27]. The GMR properties that admit these selected wavelengths into the sensor (coupled in), also allow these wavelengths out of the sensor (coupled out). Because standing waves have allowed multiple photons to constructively interfere, the efficiency of the light coupled out of the sensor is high [
27]. The light coupled out at high efficiency is referred to as being resonant. These resonant wavelength spectra are narrow and sensitive to changes on the surface of the GMR structure, for example, chemical reactions or the presence of a biomaterial [
7,
19,
27,
28].
Magnusson and Wang proposed the use of the GMR effect for sensor applications due to the GMR filter’s tunable properties, based on the resonance structure parameters and refractive indices [
27]. Tibuleac et al. and Wawro et al. introduced new GMR biosensor devices, in addition to applications of the sensors integrated with optical fibers [
7,
29]. Utilization of modal and polarization differentiation for multiparametric biosensors is a pivotal attribute of this technology [
6].
GMR sensors are highly sensitive to their resonance parameters, which is innate in the fundamental resonance effect [
6]. The resonant wavelength values of the GMR device are perturbed as the structure parameters change due to the attachment of a biomolecular layer on the device. A bio-selective layer on the GMR device can preferentially bind with a target analyte; this avoids additional data processing and foreign tags [
6]. The GMR sensor has attributes including enriched data sets, label free sensing, and economic fabrication. These are qualities that will lead to the continued application of this sensor technology in several fields [
30].
The sensor used in this work is required to have three or more resonant modes and it must be easily fabricated. To achieve the goal of three or more resonant modes, a sensor with a relatively thick homogeneous layer is developed. The thick homogeneous layer supports the resonant modes well. Here, we fabricate a 1-D grating with a two-part period. In addition, an aspect ratio is chosen so that the pillars in the grating are easily formed with a low probability of collapsing.
The sensor model shown in
Figure 1a is simple and easily fabricated with a grating aspect ratio of 0.81. In
Figure 1b,c the TM and TE spectra are shown for a simple grating (d
h = 0 nm) and for d
h values from 300 to 500 nm; these simulations use the rigorous coupled-wave analysis (RCWA) module in the RSoft DiffractMOD software (Synopsys, Inc., Mountain View, USA) [
31,
32]. The grating without a homogenous layer has no resonant peaks in the spectrum of interest. At a d
h value of 500 nm, two resonances occur in each polarization state for a total of four resonant peaks. The thicker homogeneous layer results in more resonant peaks that are used to detect changes in the multiple sensor variables. To this end, a homogeneous layer of 500 nm is chosen for this work to allow for multiparametric data collection from the sensor (
Figure 1b,c green line).
In this work we use simulation tools, with our multiparametric sensor, to generate an extensive library of biolayer thicknesses, biolayer indices of refraction, bulk indices of refraction, and the associated resonant wavelength shifts produced by our sensor. We use our library (also termed lookup table) and an inversion algorithm to take an input of the measured shifts of three resonant wavelengths (measured in one spectrometer reading), and then output the shifts in the value of three sensor variables (biolayer thickness, biolayer index of refraction, and bulk index of refraction). This method of determining biolayer and bulk sensor value shifts is novel, specifically because we simultaneously determine shifts in the values of two biolayer variables. To add to the usefulness of our method, we also determine the shifts in value of one bulk variable. The creation of a lookup table through a quick automated process and the utilization of an inversion algorithm expands what can be executed with multiparametric sensors/transducers.
2. Experiments and Methods
2.1. Fabricated Sensor Specifications
The guided-mode resonance (GMR) sensor consisted of a quartz substrate overlayed with silicon nitride (Si
3N
4) via plasma enhanced chemical vapor deposition (PECVD). The Si
3N
4 was deposited at a rate of 290 Angstrom/min. After Si
3N
4 was deposited on the substrate, it was patterned using laser interference lithography (LIL) and reactive ion etching (RIE) to achieve the fill factor and grating depth needed (
Figure 2).
The sensor described above was used to quantify Concanavalin A (Con A) and Yarrowia yeasts that bind to the sensor surface. The sensor surface was functionalized with Con A because it attracts the glycoproteins found on the cell wall of the yeasts [
4]. The reflection response of the sensor was measured with an Ocean Optics USB4000-VIS-NIR spectrum analyzer via an Ocean Optics bifurcated optical fiber. The layout of the sensor setup is illustrated in
Figure 3.
2.2. Yeast Cell Preparation
The yeast cells used were in a suspension media including 20% glycerol. The yeast cells were of the fungal genus Yarrowia in the Ascomycota phylum family. There are 4 to 6 chromosomes in the yeast cells, and the genome encodes 6448 genes [
4]. The Yarrowia species of yeast characteristically form biofilms [
33]. The yeast cells were reduced from a concentration of
to
cells/mL by diluting them in phosphate buffer saline (PBS).
2.3. Sensor Surface Functionalization
The methanol used in this work was from Sigma-Aldrich (St. Louis, MO, USA); Concanavalin A (Con A) and glutaraldehyde were from Santa Cruz Biotechnology, Inc. (Dallas, TX, USA); and (3-Aminopropyl)triethoxysilane (APTES) was from Acros Organics (Carlsbad, CA, USA).
The sensor surface was functionalized in accordance with the work carried out by Abdallah et al. [
4]. The molecules used on the sensor surface were meant to promote specific binding to the desired analyte (yeast cells) and simultaneously reduce non-specific binding. The silane in the APTES solution was used to capture the protein Con A. Con A has a high affinity to bind to the polysaccharides that are readily found on the cell wall of yeasts—the analyte of choice [
34].
First, the sensor was placed in 3% (3-Aminopropyl)triethoxysilane (APTES) and 97% methanol for 30 min to salinize the surface. The sensor was then gently agitated for 60 s in a 60:40 mixture of methanol and deionized (DI) water; it was then dried in a vacuum furnace for 30 min at a temperature of 95 °C and a pressure of 15 mm Hg. A mixture consisting of 0.7% glutaraldehyde (GA) and 99.3% DI water was prepared and the sensor was submerged in the mixture for 60 min. Subsequently, the sensor was agitated in DI water for 60 s. The sensor was then incubated in 1 mg/mL of Con A in PBS for 120 min; for part of this time, the reflection spectrum was measured, and the resonant wavelength shifts at the end of the period were recorded as the total shifts. Then, the sensor was placed in a DI water bath and gently agitated for 60 s. A dilution of yeast cells in PBS with a concentration of 8.25 × 105 cells/mL was used to soak the sensor surface for 120 min; the spectrum was measured for part of the incubation period. The wavelength shifts at the end of the recorded portion of this interval were recorded as the total shifts.
2.4. Inversion: Translating Resonant Shifts to Sensor Variables
A lookup table of simulated resonance shifts, and corresponding sensor variable value changes was used to invert resonant shift values to changes in the biolayer and bulk environment of the sensor. The RCWA module in RSoft DiffractMOD software (Synopsys, Inc., Mountain View, USA) was used for all simulations in this work [
32]. A valuable feature of RSoft is the ability to set a value range of GMR structure parameters (like depth/thickness and index of refraction) and step sizes through the value range. For example, the grating depth range can be set from 50 nm to 250 nm and the step size can be chosen as 10 nm. This would result in producing 21 simulations that iterate through grating depth values from 50 nm to 250 nm (i.e., 50, 60, 70, … 250 nm) with all other GMR parameters held constant. Each simulation consists of reflectance over a wavelength spectrum, with the wavelength range and resolution selected by the user. The RSoft software can also perform iterative RCWA calculations over multiple sensor parameters based on ranges of values and step sizes through the ranges set by the user.
To produce the lookup table in this work, the sensor’s structural parameters (period, fill factor, grating depth, and homogeneous layer depth) and indices of refraction of materials were held constant, and a biolayer was modeled on the sensor surface. The term, sensor variables, is used to refer to the values of the modeled biolayer and a value of the bulk environment. Iterative simulations were set to cycle through values of biolayer thickness, biolayer index of refraction, and the bulk index of refraction; each simulation generates the reflectance spectrum for a set of sensor variable values.
In our work, the iteration settings were as follows: the biolayer thickness range was 2 nm to 20 nm and the step size was 3 nm (i.e., 2, 5, 8, … 20 nm), the biolayer index of refraction range was 1.334 RIU to 1.430 RIU and the step size was 0.006 RIU (i.e., 1.334, 1.340, 1.346, … 1.430 RIU), and the bulk index of refraction range was 1.334 RIU to 1.430 and the step size was 0.006 RIU (i.e., 1.334, 1.340, 1.346, … 1.430 RIU). This resulted in 7 values of biolayer thickness, 17 values of biolayer index of refraction, and 17 values of bulk index of refraction. A simulation was performed for every permutation of the biolayer and bulk sensor variable values set by the user; this resulted in a total of 2023 reflectance spectra produced via RSoft simulation. We developed a MATLAB program to quickly evaluate the peak wavelength values for each of the 4 resonant modes by finding the maximum reflectance value in a wavelength range corresponding to a given mode. For a single simulation, a set of 4 resonant wavelengths, biolayer thickness, biolayer index of refraction, and bulk index of refraction were grouped together and inserted into a digital file referred to as the lookup table; this was all performed within our MATLAB program. The lookup table consisted of all permutations of the sensor variables set by the user and the corresponding resonant wavelength values for the 4 resonant modes. Our use of the iterative ability of RSoft and the MATLAB code we developed were key to taking experimental resonance wavelength shifts as input and producing GMR biolayer thickness, biolayer index of refraction, and bulk index of refraction as output. For the results reported herein, out of the set of 4 available modes (see
Figure 1), we used the 3 modes with the highest sensitivity.
The experimentally measured wavelength shifts were taken as input and compared to the wavelength shifts from the lookup table using the formula . Here, ΔλTi is the change in resonant wavelength (Δλ) from the lookup table (T) for the resonant mode (i). ΔλEi is the change in resonant wavelength (Δλ) from the experimental input (E) for the resonant mode (i). The set of resonant wavelength shifts from the lookup table for which S is minimized correlates to a set of sensor variables that were taken as the output.
Let us detail how Smin was calculated. We determined the difference between the lookup table resonance wavelength shift for resonant mode TM0 and the experimental resonance wavelength shift for resonant mode TM0; this value was then squared. This computation was also performed for the TM1 and TE1 modes. These 3 values (one for each mode used) were summed and then square rooted to get S. This calculation was performed for all 2023 lookup table entries.
Of the 2023 lookup table entries, the 10 with the lowest values of S were kept (S1 to S10) and all others were discarded. The 2nd lowest value of S (S2) was used to gauge the significance of the 8 S values that were greater (S3 to S10). Any S3 to S10 value that had a percent difference from S2 that was greater than ~30% was discarded. Statistical analysis of our data led us to using ~30% as the best cutoff point. As stated above, each S value corresponded to an entry in the lookup table. The biolayer thicknesses, biolayer indices of refraction, and the bulk indices of refraction corresponding to the remaining values of S were averaged and taken as the output biolayer thickness, biolayer index of refraction, and bulk index of refraction. Finally, to verify, these 3 physical values, inserted into the RCWA simulation software, generated the set of Δλ(simulation) that was approximately equal to the set of Δλ(experimental); this fact supports our use of this method.
4. Conclusions
Shifts in the wavelength of three resonant modes produced by a GMR sensor are used to simultaneously quantify the change in biolayer thickness, biolayer index of refraction, and bulk index of refraction. Extracting shifts in value for these three sensor variables using one spectral measurement, is an improvement in the field of biosensing. In past work, GMR sensors were used to monitor the bulk index of refraction and only one biolayer variable at a time.
The sensor used in this work is functionalized, with the protein Con A, for detection of Yarrowia yeast cells. During analyte accumulation, a set of three measured resonant wavelength shifts are used as inputs in a lookup table, and the outputs are three sensor variable value shifts: bulk index of refraction, biolayer thickness, and biolayer index of refraction. A label-free process that uses a single spectral measurement to determine changes in biolayer and bulk sensor variables is novel in practice. The multiparametric sensor, demonstrated via experiment and simulation in this work, can be used to quickly obtain detailed information about the sensor environment and the accumulation of a desired analyte over time. The key to our study is the use of a lookup table and an inversion algorithm to simultaneously monitor three sensor variables (two biolayer and one bulk) as biomolecules are captured on a sensor surface; this is a feat that has not been accomplished in past research—based on our review of the literature.
The simulation of the 2023 reflection spectra, each with a different variable value set, is an automated process using rigorous numerical models. This process can be applied to other sensors of interest and the number of variable value sets can be chosen freely, greatly exceeding the modest set count applied here. Implementing our inversion method with sensors that have a high sensitivity is expected to lead to a more accurate inversion; this would be important future work.