Next Article in Journal
Implications of the Onset of Sweating on the Sweat Lactate Threshold
Previous Article in Journal
Segmentation of Unsound Wheat Kernels Based on Improved Mask RCNN
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sensing Mechanisms of Rough Plasmonic Surfaces for Protein Binding of Surface Plasmon Resonance Detection

by
Treesukon Treebupachatsakul
1,
Siratchakrit Shinnakerdchoke
1 and
Suejit Pechprasarn
2,*
1
Department of Biomedical Engineering, School of Engineering, King Mongkut’s Institute of Technology Ladkrabang, Bangkok 10520, Thailand
2
College of Biomedical Engineering, Rangsit University, Pathum Thani 12000, Thailand
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(7), 3377; https://doi.org/10.3390/s23073377
Submission received: 13 January 2023 / Revised: 5 March 2023 / Accepted: 16 March 2023 / Published: 23 March 2023
(This article belongs to the Section Biosensors)

Abstract

:
Surface plasmon resonance (SPR) has been utilized in various optical applications, including biosensors. The SPR-based sensor is a gold standard for protein kinetic measurement due to its ultrasensitivity on the plasmonic metal surface. However, a slight change in the surface morphology, such as roughness or pattern, can significantly impact its performance. This study proposes a theoretical framework to explain sensing mechanisms and quantify sensing performance parameters of angular surface plasmon resonance detection for binding kinetic sensing at different levels of surface roughness. The theoretical investigation utilized two models, a protein layer coating on a rough plasmonic surface with and without sidewall coatings. The two models enable us to separate and quantify the enhancement factors due to the localized surface plasmon polaritons at sharp edges of the rough surfaces and the increased surface area for protein binding due to roughness. The Gaussian random surface technique was employed to create rough metal surfaces. Reflectance spectra and quantitative performance parameters were simulated and quantified using rigorous coupled-wave analysis and Monte Carlo simulation. These parameters include sensitivity, plasmonic dip position, intensity contrast, full width at half maximum, plasmonic angle, and figure of merit. Roughness can significantly impact the intensity measurement of binding kinetics, positively or negatively, depending on the roughness levels. Due to the increased scattering loss, a tradeoff between sensitivity and increased roughness leads to a widened plasmonic reflectance dip. Some roughness profiles can give a negative and enhanced sensitivity without broadening the SPR spectra. We also discuss how the improved sensitivity of rough surfaces is predominantly due to the localized surface wave, not the increased density of the binding domain.

1. Introduction

Surface plasmon resonance (SPR) is a guided electromagnetic phenomenon that occurs when electrons oscillate between a noble metal and dielectric when illuminated by an external electric field, for example, an incident light with a matching momentum [1,2]. The coupling between the incident photon and electron oscillation leads to optical energy absorption [3], resulting in a dark band in plasmonic reflectance spectra. The coupling of SPR is a strong candidate for real-time monitoring of binding interactions [4], leading to the development of sensors based on SPR along with numerous other optical sensing technologies. Examples of applications using SPR include host–pathogen detection [5,6], biochemical interactions [7,8], ultrasonic detection [9], refractive index sensing [10,11], voltage sensing [12], drug discovery [13,14], microscopic imaging [15], and therapeutic monitoring [16].
A widely adopted SPR-based sensor consists of a glass substrate with a refractive index of n0, a noble metal film with a refractive index of nm, and a sensing region with a refractive index of ns. This structure, also known as the Kretschmann configuration [17], is depicted in Figure 1a. There are several measurement schemes for SPR, including angular interrogation detection [18], wavelength scanning detection [19], intensity detection [20], and phase detection [21].
The angular interrogation detection [18] measures reflectance while varying the incident angle θ0 at a fixed wavelength of λ, illuminating the thin plasmonic metal due to the explained light-matter interaction between photons and electrons. It manifests as a dark band in the reflectance spectra at a specific incident angle, known as the plasmonic angle θsp, as depicted in Figure 1b. Furthermore, alterations on the metal film surface, such as molecular binding at an active site on the sensor surface, resulting in a subsequent alteration in the resonant coupling condition and a shift in the plasmonic angle, as depicted in Figure 1b.
Figure 1b shows the plasmonic angle shift from the water-backing environment with a refractive index of 1.33 to a 5 nm thick bovine serum albumin (BSA) with a BSA refractive index nBSA of 1.35 for an 80 mg/mL concentration [22] adhering to the metal surface, respectively. Figure 1b was calculated using the computational method explained later in Section 2.
Multiple studies [23,24,25,26,27,28] have reported that the surface profile, including the thickness and roughness of the metal film, significantly affects the performance of the SPR-based sensor for bulk refractive index sensing. For example, Kurihara et al. [27] explained that the SPR dips could be classified into large, small, and equivalent relationships between absorption and excitation frequency, depending on the frequency relationship between the optical property of the light source and the metal layer thickness. Furthermore, Hoffman et al. [26] demonstrated a significant change in dielectric functions and SPR wave vectors as the roughness amplitudes increased. Our recent publication [28] found that the sensitivity, full width at half maximum, intensity contrast, and plasmonic dip position for measuring the refractive index using an SPR-based sensor with rough gold films leads to a degradation of the sensor’s overall quality.
Several papers [29,30] reported that the sensitivity of binding kinetic measurement increases with the roughness level due to a larger surface area and higher surface protein in contact with the plasmonic metal. It is crucial to note that branching molecules, such as dextran, have been employed in binding surface preparation to increase the binding density. Unlike in rough plasmonic surfaces, for example, Liu et al. [30] employed a Ni seed layer on a silver-based SPR sensor to diminish the surface roughness, resulting in a smaller FWHM. In addition, Byun et al. [29] analyzed the localized surface plasmon resonance (LSPR) on a nanowire-based structure with different degrees of roughness. They also concluded that there was a change in the extinction peak amplitude caused by the LSPR effect.
Various surface treatment techniques have been employed as a postprocessing step for metal film fabrication to reduce surface roughness. Conventionally, the metal layer of the SPR-based sensor was manufactured using a sputter coating, which has a root-mean-square (RMS) roughness of 1.4 nm to 1.5 nm [31]. The metal film can be further smoothed through surface treatments, such as thermal annealing [32], helium ion beam [33], laser ablation [34], chemically grown single-crystalline gold [35], mica substrate utilizing [36], and chemical polishing [37,38]. Table 1 summarizes the roughness values for the mentioned treatment methods.
This paper theoretically investigated multiple roughness levels on a uniform gold sensor to analyze roughness’s effects on the SPR protein binding performance. This study has chosen the gold material due to its chemical stability and biocompatibility for biochemical and bioscience-based applications. A theoretical framework is proposed to analyze comparative sensing performance parameters, including protein binding sensitivity (S), the plasmonic dip position (n0sinθsp), the full width at half maximum (FWHM), the intensity contrast (ΔI), the average plasmonic dip intensity (Isp), and the figure of merit (FoM). The theoretical investigation was based on rigorous coupled-wave analysis [39,40] and Monte Carlo simulation [41].
It is significant to note that the research differs from previous literature regarding protein binding to SPR sensors with different roughness levels. We systematically separate and quantify the two enhancement schemes: the LSPR and the increased binding density. The proposed method provides insight into the underlying performance enhancement of SPR sensors for binding kinetic measurement. In addition, we believe that the research provides a complete explanation of the impact of surface roughness on the molecular binding application of the SPR-based sensor performance and its limitations.

2. Materials and Methods

2.1. Rough Surface Model and Simulated Structures

The rough gold surface was constructed using numerical modeling with Gaussian random surfaces [42], which were adopted from Byun et al. [29], and characterized by varying maximum roughness height (h) and correlation length (cl) from 2 nm to 20 nm and 5 nm to 50 nm, respectively. The roughness factors were selected to cover SPR excitation and its cut-off positions. Different mode characteristics are present in the chosen roughness scope, including effects from LSPR and enhanced protein binding density, and these will be discussed later in the results section.
A function s(x) consisting of digital numbers 0 and 1 in a random sequence was first developed, where the term x indicated the spatial distance along the x-axis of the sensor to generate the rough surface. The function was then multiplied by the maximum roughness height, h. The generated profile was then converted to the frequency domain using Fourier transformation, expressed as S(ω) = F[s(x)], where S(ω) is the Fourier-transformed output, ω is the frequency domain axis, and F is the Fourier-transform operator. The S(ω) function was then multiplied with the Gaussian distribution function, G(ω) = exp(−x2/(cl2/2)). The G(ω) function essentially acts as a low-pass filter for the signal, where the term cl indicates the roughness correlation length. Next, the inverse Fourier transformation was applied to convert the signal back to the spatial domain signal and rearrange the function into a 200 columns rectangular contour. It is worth noting that higher cl values indicate a smoother surface. Therefore, a theoretical uniform surface is a Gaussian random surface with h equal to zero and cl approaching infinity.
The rough surface was then rearranged into a rectangular coordinate for rigorous coupled-wave analysis simulation by dividing the surface into multiple layers with a height of 0.1 nm each. The rough layer was then assembled with the uniform gold surface of the sensor. The rough surface structure profile consisted of three regions: (1) the yellow area beneath the surface indicating the plasmonic gold material with a refractive index of 0.18 + 3.43i [43], (2) the light green contour indicating the 5 nm-thick BSA binding layer [44] with a refractive index of 1.35, equivalent to 80 mg/mL concentration [22], and (3) the blue region indicating the water-based environment sensing area with a refractive index of 1.33. BSA was chosen as a substrate layer in this investigation due to its numerous applications, including cell culture and studies of biochemical reactions. Furthermore, BSA, including in biomedical research, has been widely adopted as a protein concentration standard in laboratory experiments.
The theoretical study was conducted using three sets of samples:
(1) no BSA protein layer case, as shown in Figure 2a;
(2) the second sample was a rough surface with only a BSA layer on top with no sidewall coating, as depicted in Figure 2b. This model indicated an equivalent amount of protein to the ideally smooth metal, allowing us to evaluate the enhancement due to the LSPR without the influence of the increased amount of protein due to the roughness;
(3) Conversely, the BSA layer in the third sample type was coated on top of the structure and its rough-surface sidewall, as illustrated in Figure 2c. The investigation of this structure aimed to replicate the protein solution binding to the entire surface, leading to an increased protein sample binding concentration. However, the amount of BSA in the third sample can be significantly denser than in the second sample when the surface is highly rough, as explained later in Section 3.3. Therefore, the increased protein and the LSPR can enhance the SPR response in the third sample, which the second model can suppress to quantify and separate the two enhancement mechanisms.
The sidewall coating behavior of proteins in nature was more complicated than in the simulation. Even though sidewall coating is typically used to prevent edges from coming into direct contact with water and to increase binding density, this assumption may not always be accurate. Even with sidewall coating, edges can remain in direct contact with water, which could affect the precision of the provided data. In such a way, the enhanced sensitivity and FoM will be between models (2) and (3). The results section will discuss that despite the absence of sidewall coating, roughness profiles can improve the SPR measurement due to the LSPR. However, many studies on polymer technology [45,46,47,48] assist in controlling protein fabrication. Saleem et al. [45] and Nuutinen et al. [46] investigate the application of polymeric material on resonant wave gratings (RWG). They utilized TiO2 as a cover layer on the gratings structure, and then the protein layer was uniformly bound to the material. Mayet et al. [47] published a similar study in which an oxide-nitride-oxide passivation layer evenly coated the sidewalls of the gratings. Zhang et al. [48] also conducted a theoretical study involving protein sidewall coating on a metal-layer-assisted double-grating (MADG) biosensor. These indicate that the experimental and theoretical information agree well and that the protein coating behavior can be controlled.
Figure 3 illustrates the overall simulation diagram, which consists of a BK7 glass substrate with a refractive index of 1.52, a uniform gold surface with a thickness of 50-h nm, a mean rough gold layer with a thickness and depth of h, a water-based environment sensing area with a refractive index of 1.33, a binding analyte (an 80 mg/mL concentrated BSA protein) with a refractive index of 1.35, and a thickness of 5 nm lining the top and sidewall of the rough surface layer, and a linearly p-polarized HeNe laser with a wavelength of 633 nm. The laser illuminates the gold film at an incident angle, scanning from the critical angle to 90°. According to the sampling theory, the unit cell length is 1 μm and consists of multiple rows with 1 nm each and 200 columns to ensure that the unit cell can provide different levels of the employed roughness with the minimum sampling feature size of 5 nm. The proposed 200 columns to represent the roughness in this study will be discussed and validated later in the results section.
As depicted in Figure 3, roughness profiles simulated in this study were simplified to a one-dimensional roughness profile, where the y-axis was treated as an infinite extension of the x-z plane. The Krestschmann configuration in this study was a prism-based Kretschmann configuration with TM polarization illumination. The TM incident electric field can be resolved in the x-z plane, with no electric field along the y-axis. That is the reason for the uniform y-axis in this study. It is essential to point out that this assumption may not hold if the SPR illumination is achieved using a high NA objective lens; there will be p-polarization components along the y-axis [49]. It will be shown in the results section that the proposed simplified simulation can give similar SPR reflectance spectra compared to the experimental results reported in the literature for the prism-based configuration and fabricated two-dimensional rough-gold sensors at different roughness levels.

2.2. Simulation of the Surface Plasmon Resonance Detection

This research employed the rigorous coupled-wave analysis (RCWA) function [40,50] and the Monte Carlo simulation to calculate the reflectance curve under a MATLAB 2022b environment, utilizing parallel computing and graphic processing toolboxes. In addition, 151 diffraction orders were applied to ensure that all results from the studied surfaces achieved convergence, which will be further discussed in Section 3.1.
Since a single pair of h and cl can generate multiple surfaces, the Monte Carlo simulation was applied to estimate an average response of 100 surface profiles for the three sample types. These responses from the rough surfaces generated by h and cl were then analyzed using the RCWA. Thus, each Monte Carlo simulation consists of performance factors’ values from 100 times RCWA simulation with a hundred different generated rough surfaces using a single value of h and cl, as illustrated in the flow chart in Figure 4. This procedure allows for a high accuracy and stability of the calculated performance parameters across all values of the roughness factors.

2.3. Quantitative Performance Parameters

The studied performance parameters of surface plasmon resonance detection for protein binding applications included sensitivity (S), full width at half maximum (FWHM), intensity contrast (ΔI), average plasmonic dip intensity (Isp), and figure of merit (FoM). This section provided definitions and descriptions of the presented performance factors; these parameters will be computed to explain the SPR sensing mechanisms for the protein binding in the results section.
(1)
The sensitivity (S) was calculated by dividing the change in the plasmonic wave vector (k), also known as the plasmonic shifting distance, by the product of the protein layer thickness and the difference in the binding region’s refractive indices, as shown in Equation (1) and Figure 5b.
S = Δ k d Δ n s = 2 π n 0 Δ s i n θ s p λ d B S A Δ n s
The terms n0, θsp, dBSA, λ, and Δns are the refractive index of the glass substrate (1.52), the plasmonic angle, the 5-nm BSA protein layer, the wavelength of the incident light (633 nm), and the difference in sample refractive indices of 0.02, respectively.
(2)
The n0sinθsp indicates the angular position at the minimum plasmonic intensity dip (Isp), which can be used to analyze the sensor’s detection range.
(3)
The FWHM in this manuscript was calculated as the average width of the two plasmonic dips, measured at 50% of the reflectance spectra’s intensity, as expressed in Equation (2). Figure 5a illustrates the FWHM of the SPR curves, depicted as a black arrow considering the unsymmetrical nature of SPR dips.
F W H M = Δ k n o r , w a t e r + Δ k n o r , B S A 2
The terms  Δ k n o r ,   w a t e r  and  Δ k n o r , B S A  are the FWHM of the normalized reflectance curves, measured at half of the maximum magnitude when the refractive index of the sensing region was 1.33 and 1.35, respectively. In addition, several techniques can be applied to measure the FWHM of the unsymmetrical SPR dip, such as Lorentzian fit [51], Green’s function [52], and Horner’s curve fitting [53].
(4)
The ΔI is one of the unique parameters which helps to determine the quality of the surface plasmon resonance-based sensor. It can be calculated as the average difference between the reflectance at the plasmonic angle and the intensity at the critical angle of the reflectance spectra, as expressed in Equation (3) and illustrated in Figure 5a.
Δ I = Δ I w a t e r + Δ I B S A 2
The terms ΔIwater and ΔIBSA are the difference in the intensity at the plasmonic angles of the reflectance spectra detected when the 5 nm thick BSA binding layer is absent and present, respectively.
(5)
The Isp is another factor for considering the sensor’s performance. It can be defined as the average of the reflectance values at the plasmonic angles as expressed in Equation (4) and shown in Figure 5b.
I s p = I s p , w a t e r + I s p , B S A 2
The terms Isp,water, and Isp,BSA are the intensity at the plasmonic angles of the reflectance spectra computed without and with the binding substrate.
(6)
FoM encapsulates an overall sensor quality, usually defined as a ratio between the S and the FWHM [54], as expressed in Equation (5).
F o M 1 = S F W H M
It is interesting to point out that the FoM depends on the detection mechanism for intensity-based measurement; the FoM is calculated based on the S and FWHM and can also be related to intensity parameters, such as ΔI and Isp, and considering the shot noise model [55], as expressed in Equation (6).
F O M 2 = S × Δ I F W H M × I s p 4

2.4. Enhancement Ratio

Later in the result section, we will discuss that the FWHM, ΔI, and Isp are comparable for both the nonsidewall and sidewall models, indicating that the additional protein from the roughness does not affect these parameters. The S and FoM are the only performance parameters that can be affected by the sidewall protein. Here, the enhancement ratio denoted by ER compares the roughness models’ relative change in S and FoM to the ideal uniform gold sensor response. Furthermore, the LSPR wave, protein increment, and roughness effect will be quantified and classified using ER as expressed in Equations (7) to (9). Equations (7) to (9) show the calculations for the roughness enhancement effect, LSPR wave effect, and protein increase, respectively.
E R S , r o u g h n e s s = S s i d e w a l l S u n i f o r m
E R S , L S P R = S n o n s i d e w a l l S u n i f o r m
E R S , p r o t e i n = S s i d e w a l l S n o n s i d e w a l l S u n i f o r m = E R S , r o u g h n e s s E R S , L S P R
The terms ERS,roughness, ERS,LSPR, and ERS,protein indicate the sensitivity enhancement ratios caused by roughness (the LSPR and the increased binding density), the LSPR wave effect, and increased protein binding density, respectively. Furthermore, Suniform, Ssidewall, and Snon-sidewall are the sensitivities obtained by the ideal uniform gold model, BSA with the sidewall model, and BSA without the sidewall model, respectively. The enhancement for the FoM can be calculated using the relative change in the FoM, similar to the sensitivity.

3. Results

3.1. Convergence Test for Extreme Cases

The convergence test on the optimal diffraction orders was employed in four extreme cases of the study, consisting of the rough surfaces generated using (a) h of 2 nm and cl of 5 nm, (b) h of 20 nm and cl of 5 nm, (c) h of 2 nm and cl of 50 nm, and (d) h of 20 nm and cl of 50 nm, as shown in Figure 6a. The test was performed using the rigorous coupled-wave analysis to calculate the p-polarized HeNe laser’s optical reflectance with multiple diffraction orders ranging from 1 to 171. The numerical fluctuation between the 149 and 151 diffraction orders was 7 × 104, 0.002, 0.020, and 0.003 for the four tested surfaces. Employing a higher number of diffraction orders will surely increase the accuracy of the results; however, it also requires longer computing time and more resources. Therefore, the diffraction orders of 151 were employed in this study’s simulation cases with a numerical stability of 0.01 for the reflectance calculations.
The convergence test also inspected the adequate number of columns representing the rough surfaces, as shown in Figure 6b. The number of columns was investigated from 10 to 220 for the average plasmonic intensity of the no BSA model, BSA without the sidewall model, and BSA with the sidewall model. The tested roughness model was constructed based on the same four extreme cases. The plasmonic intensity’s differences between the 200 and 220 columns for (a) h of 2 nm and cl of 5 nm, (b) h of 20 nm and cl of 5 nm, (c) h of 2 nm and cl of 50 nm, and (d) h of 20 nm and cl of 50 nm was 7 × 104, 0.002, 0.018, and 0.005 for the proposed three models. The following paragraphs will discuss how this roughness level can enhance protein binding. The analysis demonstrates that 200 columns are enough to simulate the roughness frameworks.

3.2. Comparison between the Proposed RCWA Simulation and Monte Carlo-Based Method and Reported Experimental Results in the Literature

One crucial question that we will discuss in this section is how accurate or how inaccurate the proposed RCWA simulation and Monte Carlo-based method can be compared to the experimental results. Yang et al. [56] reported experimental SPR reflectance measurements for different roughness levels for 632.8 nm incident wavelength and 57.5 nm gold films. These gold films were annealed using thermal annealing at four different temperatures at 100 °C to 400 °C for 15 min, allowing them to vary the roughness of the films. Their roughness parameters were measured using an atomic force microscope (AFM).
Figure 7 shows SPR reflectance spectra calculated using the proposed procedure using the roughness levels and the roughness heights reported by Yang et al. [56]. Table 2 describes the roughness levels, including the root-mean-square (RMS) roughness and the roughness heights, experimental plasmonic angles, and the plasmonic angles calculated from Figure 7. The plasmonic angles of the experimental results and the simulated SPR dips agree well, with slight discrepancies within 0.09 degrees. The slight error shown in Table 2 could be due to the gold used in the simulation may differ from the actual gold properties in the experiment.

3.3. Quantitative Performance Parameters of SPR Sensors at Different Roughness Levels

Ideally smooth gold SPR-based sensor
All quantitative performance parameters were initially calculated based on Equations (1) to (5) after the RCWA on the ideally smooth gold layer for the SPR-based sensor was simulated. As a result, the computed S, plasmonic dip position (n0sinθsp), FWHM, ΔI, Isp, FoM1, and FoM2 were 122.59 rad/µm2, 1.440, 0.41 rad∙RIU/µm, 0.90, 0.007, 299.00, and 980.66, respectively. These are the theoretical limits of the quantitative performance parameters of an ideally uniform SPR gold sensor calculated from Figure 5a.
Rough gold surfaces
Figure 8a,b illustrate the S of the SPR on a nonsidewall BSA model and a sidewall BSA model to different roughness using RCWA and Monte Carlo simulation computed using Equation (1). Furthermore, the increment of the BSA binding domain, in percent, in the sidewall model compared to the BSA-on-top model at different roughness was shown in Figure 8c. Two gratings’ models with equivalent h and cl were generated using the same principle to calculate the protein increment. The additional sidewall protein can be calculated by subtracting the nonsidewall model from the sidewall model. Later, the percentage of protein increment can be determined by dividing the sidewall protein by the nonsidewall model’s protein amount. According to the contour representation, the protein difference remained negligible at the low-to-moderately rough surface. However, the increased amount of protein due to the roughness significantly increased as the roughness height, h, escalated to over 10 nm for a short correlation length, cl.
The S can be classified into four major regions:
(1)
The “no sidewall protein enhancement” corresponds to the region where the letters ‘a’ and ‘b’ are in Figure 8. This area indicated the roughness levels in which the sidewall BSA was negligible. The protein lying on top of the rough surfaces overlapped the sidewall protein. Moreover, two types of SPR results were investigated in this region: the negative-sensitivity LSPR and the degraded-sensitivity reflectance spectra, as illustrated as the operating point ‘a’ and ‘b’, respectively. The negative plasmonic dip shifting means the plasmonic angle shifts towards a lower plasmonic resonance angle when the sample refractive index increases. It is established that this is due to the adverse diffraction orders of gratings or scattering surfaces [57]. The reflectance spectra at the operating point ‘a’ (h and cl of 3 nm and 10 nm, respectively), as shown in Figure 9a, had an S of −379.88 rad/μm2; approximately threefold of the magnitude of the S acquired from the ideal uniform sensor. Within the same region but at rougher surfaces, with h of 10 nm and cl of 30 nm (operating point ‘b’), the S returned to the positive value but degraded to 85.31 rad/μm2, indicating a 30.41% decrease in the sensitivity. The averaged reflectance spectra calculated from 100 structure profiles at these roughness levels are shown in Figure 9b;
(2)
The “positive-sensitivity LSPR” only described the S affected by the LSPR. Here, the additional protein was presented but had an insignificant effect on the S. The reflectance curves spectra at the operating point ‘c’ obtained at h of 12 nm and cl of 20 nm indicated the enhanced S due to the LSPR wave effect only. Figure 9c illustrates the reflectance spectra at this location. In addition, the two employed models’ curves in this roughness level did not look significantly different due to the small increment of the BSA amount;
(3)
The “LSPR and protein” was investigated as the roughness rose. In this region, S of the SPR spectra can improve due to the localized LSPR effect and the presence of an additional protein. The models indicate a significant difference in sensitivity values at the operating point ‘d’, located at h of 9 nm and cl of 8 nm. For the non-sidewall BSA model, the S of 1191.13 rad/μm2 was investigated. In contrast, the BSA sidewall model achieved a significantly higher S of 1608.02 rad/μm2, having a 35.00% increase from the nonsidewall BSA model, as shown in Figure 9d. The binding sensitivity enhancement due to the LSPR was 9.72 times at this roughness level.
In contrast, the BSA sidewall model enhanced the S by 13.12 times more than the ideally smooth gold sensor, indicating that the S enhancement due to the increased protein was 3.68 times. In addition, this crucially high S was due to the deterioration of a plasmonic dip structure, resulting in a lower n0sinθsp for the nonprotein coated model, which will be explained in the next part. The S enhancement factor due to additional protein agreed with the estimated protein concentration at the roughness level, as shown in Figure 8c;
(4)
At the extreme roughness levels (near the bottom right corner of the contour in Figure 8), the plasmonic dip structure deteriorated, as indicated by the “No SPR” region, due to the scattering loss [58] of the rough surface. The propagation length of surface plasmon polaritons is strongly distorted due to roughness [59], resulting in an undetectable region, as indicated in Figure 9e for the operating point ‘e’ (obtained at h of 16 nm and cl of 15 nm). The plasmonic dips have virtually no intensity contrast. The SPR reflectance dips weakly maintained their structure and distorted at the higher roughness. The increase in protein quantity due to roughness was not the only reason for the improved S but also the enhancement of the LSPR from the roughness peaks. There is a substantial tradeoff between the roughness level, the sensitivity, and the scattering loss of the SPR dip.
Figure 10a illustrates the average n0sinθsp of the SPR dip positions calculated from the BSA without the sidewall model, the BSA with the sidewall model, and the without BSA coating model. The parameter can be used to identify the sensor’s detection range. An increased surface roughness could result in a higher plasmonic angle and a shorter detection range. Moreover, a higher n0sinθsp results in a bigger numerical aperture (NA) requirement of the illumination lens for angular scanning applications. Figure 10b compares the reflectance curves from the model without the BSA coating and their plasmonic dip position at h of 0 nm (ideal smooth), 5 nm, 10 nm, and 15 nm, and the equivalent cl of 20 nm. As the theoretical uniform surface became rougher to 10 nm, the n0sinθsp shifted gradually. At this cl, the shifting distances from h of 0 nm to 5 nm and h of 5 nm to 10 nm were very similar at approximately 0.05. However, when the h is increased to 15 nm, the n0sinθsp moves slightly backward due to the deterioration of the plasmonic structure due to surface roughness.
Figure 10c shows the FWHM of the reflectance spectra at the simulated rough surfaces. The parameter was computed using Equation (3). The FWHM results for the sidewall and no sidewall cases were similar; therefore, only the FWHM obtained from the non-sidewall BSA model was presented to shorten the manuscript. The growth of roughness strongly broadens the size of the reflectance spectra. At slightly rough surfaces with an h lower than 5 nm, the FWHM was affected by h only. For instance, operating point ‘a’ had an FWHM of 0.58 rad∙RIU/µm, 41.46% higher than the ideal uniform sensor’s FWHM. However, when h was greater than 5 nm, cl significantly affected the alternation in FWHM. At rough surfaces, with a cl less than 35 nm, the reflectance spectra may have a deteriorated structure, reducing ΔI. The crucially low ΔI can result in a slightly lower FWHM than the region with a larger cl (smoother surface). These changes in ΔI will be discussed further in the following paragraphs. In addition, operating points ‘b’, ‘c’, and ‘d’ were in this region.
On the other hand, at the operating point ‘e’, there was a crucial loss in the SPR structure, resulting in an unidentified FWHM. However, the intensity depth Isp of the SPR dips considerably deteriorated due to the scattering loss compared to the smooth sensor surface. Apart from the mentioned regions, it can be increased up to 0.83 rad∙RIU/µm at an h and a cl of 18 nm and 30 nm, respectively, which was approximately two times wider than the FWHM of the uniform gold SPR detection, as shown in Figure 10d. The rougher surface leads to stronger scattering and a broader SPR dip.
Figure 11a,b indicate the average ΔI and Isp at various roughness levels. The results obtained from the BSA without the sidewall model and the BSA with the sidewall model were slightly distinct, with a maximum difference of approximately 0.1 for both parameters, as shown in Figure 11a,b. Therefore, only the average values were reported to shorten the manuscript. In addition, this indicated that the presence of a sidewall protein does not significantly affect the change in ΔI and Isp. However, the surface roughness still substantially impacted the sensor’s output. At an operating point ‘d’, for instance, the average ΔI was 0.089, which was 10.08 times less intensity contrast than ΔI achieved by the theoretical uniform gold sensor. Furthermore, the average Isp crucially rose to 0.74 (105.71 times compared to the ideally smooth gold sensor’s Isp).
The reported quantitative performance parameters were combined to compute the FoM, calculated based on Equation (5), at different roughness levels for BSA with no sidewall model and BSA with a sidewall model, as shown in Figure 12a,b, respectively. Figure 12c,d show FoM2 calculated using the intensity-based FoM expressed in Equation (6) [55]. The protein enhancement due to the roughness neither distorts nor enhances the LSPR in the range of studied rough surfaces.
The FoM1 in Equation (5) depends mostly on how far the SPR dip moves S, not the optical intensity. Moreover, the results can be interpreted as the degradation and enhancement trends due to LSPR relying on surface roughness. In contrast, the additional protein can enhance the performance parameters, including FoM, at the expense of a more demanding NA requirement. The sensor achieved its maximum FoM1 and FoM2 of 1777.81 and 2400.02 at an h of 9 nm and a cl of 8 nm (operating point ‘d’) for the nonsidewall BSA and sidewall BSA models, respectively. The reflectance spectra at this point had a significantly increased S and a slightly dropped FWHM due to both the LSPR effect and the protein increment, resulting in the crucially high FoM1. In addition, a negative FoMs region resulted from the negative S values. The FoMs in some of these reported operating positions could surpass the sensing capability of the ideal uniform gold sensor. For instance, the operating point ‘a’ had an FoM1 of −654.97, about two times higher than the theoretical smooth gold sensor’s FoM1. Other operating points’ FoMs and other performance parameters are indicated in Table 3.
The SPR detection’s FoM2, calculated involving S, FWHM, ΔI, and Isp, had similar results, as illustrated in Figure 12c,d for the BSA without the sidewall and BSA with the sidewall model, respectively. The highest FoM2 for the SPR detection was at the same operating point as the FoM1, with an FoM2 of 616.84 for the nonsidewall BSA model and 1020.33 for the sidewall BSA model, respectively. Therefore, the surface roughness can slightly improve the overall sensing performance in the intensity detection of the SPR-based sensor for the binding kinetic measurement.
Table 3 compares the seven investigated performance parameters obtained from the five operating points. It is seen that the surface roughness can result in an enhancement (operating points ‘c’ and ‘d’), a degradation (operating point ‘b’), or a complete deterioration (operating point ‘e’) of the SPR detection due to both the LSPR effect and the additional protein. The S can be significantly increased with highly rough surfaces in exchange for the substantial deterioration in the application relying on intensity detection. There is a tradeoff between the S enhancement due to the increased protein concentration, in other words, binding density, and the FWHM of the SPR reflectance dip. Despite this limitation, the negative-sensitivity detection (operating point ‘a’) could also be an efficient alternative for utilizing the SPR intensity-measurement applications due to its improved sensitivity and slightly degraded FWHM.
The ER and the increased protein density and surface roughness effects on the S and FoMs are summarized in Table 4. The ER for the roughness effect, the LSPR effect, and the protein increment were calculated based on Equations (6) to (8). The only factor affecting the sensor’s performance for operating points ‘a’ and ‘b’ was the LSPR since the protein density was negligible, as explained in the earlier section. Note that some ERs for the operating points ‘b’, ‘c’, and ‘d’ are less than 1.00, indicating that the roughness level experienced a degradation effect. On the other hand, the additional protein and the LSPR affected the operating points ‘c’ and ‘d’.
The roughness level of the h of 9 nm and the cl of 8 nm, corresponding to the highest FoM1 and FoM2 responses discussed earlier, was further analyzed to validate their sensing response. One crucial performance parameter for protein kinetics is to check whether the sensor has a linear responsivity for different protein concentrations, such as the thickness deposited on the sensor surface. Figure 13 compares the linear responsitivity of the three models for the h of 9 nm and the cl of 8 nm roughness sensor computed in this study for the different deposited thicknesses of BSA from 0 nm (bare gold) to 10 nm. The three models had linear responsivity covering the protein thickness, with R-square values for the linear equation fitting of 0.9985, 0.9990, and 0.9990 for the ideal smooth gold surface, the rough surface with no sidewall protein coating, and the rough surface with a sidewall protein coating, respectively. The slopes (Δsinθsp/dBSA) of the three linear lines were 8.13 × 10−4, 1.69 × 10−3, and 2.01 × 10−3, for the three models, which can be converted to the S in Equation (1) and it was found that the S values were 122.59 rad/μm2, 1202.43 rad/μm2, and 1595.30 rad/μm2 for the uniform gold surface, the BSA without the sidewall model, and the BSA with the sidewall model for the 5-nm thick BSA layer case. The S values agree with the S discussed in Table 3, indicating that the rough surface can provide a linear responsivity covering 0 nm to 10 nm. Having mentioned that the uniform sidewall protein coating assumption may not be accurate, Figure 13 provides a feasible approach to predicting the expected S if the protein does not fully cover the sidewall. The partial sidewall coating’s S will be between the no sidewall coating case (red line) and the uniform sidewall coating (yellow line).

4. Discussion

4.1. Effect of LSPR and Sidewall Protein from Roughness

Table 4 shows that the wave effect was responsible for most of the impact on sensor quality. For S and FoM1, all roughness levels that preserve the SPR spectrum have a percent impact of wave effect greater than 70%. The FoM2, which considers the ΔI and Isp, experienced a greater impact from the sidewall protein. Furthermore, the percentage of S and FoM1 enhancement of the operating points ‘c’ and ‘d’ in Figure 8c were close to the percentage of protein increment due to roughness. The calculations using the explained framework using the RCWA and Monte Carlo simulation agree well with the additional protein calculated directly from the grating profile. The agreement provides an independent measurement of the protein density enhancement due to the roughness and tells us that the enhancement due to the LSPR and the increased protein are independent. This investigation implies that increasing the protein amount increases the S and the FoMs at the same rate.
In this paper, we have explained and provided insight into the detection mechanisms of rough surfaces. Rough surfaces can enhance the binding domain density, S, and FoM1 with an expense of FWHM and optical intensity contrast of the SPR dip. It is essential to point out that the advantage of rough surfaces is that they provide only a single SPR dip or SPR mode from scatterings of surface plasmon polaritons (SPPs). In other words, an average effect of scattered SPPs and LSPR hotspots at the sharp edges [60]. However, rough plasmonic surfaces also have drawbacks, including fabrication repeatability, nonlocalized heat sources due to the localized hotspots at the sharp edges, and nonuniform distribution of the SPPs, making it difficult to be employed in optical imaging applications. Another alternative approach to achieve a similar enhancement of the binding domain density, S, and FoMs is plasmonic grating. Plasmonic gratings allow the SPPs to propagate in a well-defined path and support the additional sidewall protein. The plasmonic gratings can support and excite multiple modes due to the diffraction orders leading to a limited detection range due to crosstalks and overlapping modes [61].

4.2. Intensity and Phase Detection Schemes

Aside from angular interrogation, the SPR measurement can be done using the wavelength scanning method through intensity and phase detection, as in the angular interrogation investigated in this paper. In addition, these schemes are not affected by scanning wavelengths and angles. Here, the five discussed operation points and the ideal uniform surface sensor were analyzed for their intensity detection and phase detection performances, as illustrated in Figure 14. Figure 14a shows that the roughness profiles can generate a more significant optical intensity contrast than the theoretically smooth plasmonic layer case. For instance, the operating point ‘b’, which has significantly low sensitivity; on the other hand, can provide a crucial change in the plasmonic intensity of 0.27. The phase detection results have an alternative outcome, as shown in Figure 14b. Only the operating point ‘a’ (negative sensitivity SPR) has an outstanding performance in phase detection, while the other operating positions have degraded their phase responses. According to the analysis, the negative sensitivity measurement of SPR with a rough surface can be a strong candidate for binding kinetics measurement.

4.3. Field Plots on Exemplary Optical Structures

Light possesses both an electric field and a magnetic field, and this property of electromagnetic waves was used to explain the reason behind the performance improvement of the rough SPR-based sensor. For TM-polarized light, the polarization can excite the SPR effects, or the electric field will affect the metallic surface along the x and z axes, while there are no magnetic fields in the two axes. The total intensity of the electric field along the x-axis (|Ex|2) and z-axis (|Ez|2) provides a good indicator for the summation of the LSPR strengths across the surface.
Figure 15 illustrates the electric field in the x-direction and z-direction for the TM-polarization case. Figure 15a–c illustrate the total electric field intensity of the SI unit, |Ex|2 + |Ez|2, in three different SPR-based sensor surfaces, including a theoretical smooth gold film, a gold grating structure, and a rough gold surface, respectively. For the uniform SPR-based sensor, the electric field intensity was evenly distributed along the gold surface with a maximum magnitude of 11.28  ×  103, as shown in Figure 15a. On the other hand, the total intensity of the electric field for a gratings sensor can reach up to 16.08 × 105, which is a hundred times greater than the field intensity of the ideal smooth sensor case. Note that the color scales in Figure 15b,c have been adjusted to clearly show the electric field intensity across the entire surface (their maximum values could exceed the provided color bar limit). The gratings’ height, period, and fill factor of 24 nm, 500 nm, and 0.5 were selected according to the analyzed rough surface model.
Moreover, the electric field with the significantly high intensity only appears on the sharp edges of the gratings, while the intensity in the groove area ranges from 6.00  ×  103 to 30.00  ×  103. The rough surface in the operating point ‘d’ (h and cl of 12 nm) was chosen to investigate the magnetic field behavior due to the highest achieved FoM. The maximum electric field intensity obtained can be varied from approximately 14.92 × 105 to 17.51  ×  105, depending on the randomly simulated surface. Even though the gratings and the rough surfaces produced similar electric field intensities, the roughness model provided significantly higher sharp edges, resulting in more LSPR hotspots and increased sensitivity.

5. Conclusions

The SPR-based sensor was responsive to changes in the plasmonic sensor surface. Therefore, a slight change in the surface structure, including the roughness, can strongly affect the sensor’s performance. This manuscript analyzes the effect of roughness on binding kinetic measurement for SPR detection by employing the RCWA and the Monte Carlo simulation on the constructed rough surface based on the Gaussian random surfaces method. The simulation structure consisted of two models: the nonsidewall BSA model, indicating the same amount of protein as the roughness changed, and the sidewall BSA model, where the protein quantity increased with the roughness. The convergence test on the number of required diffraction orders and simulated columns to represent the roughness and the SPR phenomenon was also applied to obtain the optimum environment for the theoretical study. The calculated reflectance spectra were then observed. Then, the performance factors were observed and calculated, including S, n0sinθsp, FWHM, ΔI, Isp, and FoMs.
The surface roughness can improve, degrade, or deteriorate the SPR sensing performance due to the generated LSPR and the additional protein amount. The sensor can have an exceptionally high S with highly rough surfaces due to the distorted plasmonic dip structure, which reduced the n0sinθsp and field enhancement effect from the surface roughness. However, utilizing the rough surface shows a tradeoff between the improved S and the crucial deterioration in the intensity detection capability. One suggested approach for the SPR application is to utilize the negative-sensitivity SPR detection at slight roughness levels. The reflectance spectra obtained from this region can increase the S to threefold of the ideal uniform gold sensor’s S while maintaining the FWHM, ΔI, and Isp. The operating points for the highest FoM1 and FoM2 enhanced the FoMs by 8.03 and 1.04 times compared to the sensor with a theoretically smooth surface. In addition, the roughness level could be utilized in other SPR measurement schemes, including intensity and phase detection. Further analysis indicates that the sensing improvement and degradation were primarily affected by the LSPR effect of over 70% in all operating cases. In conclusion, the knowledge of surface roughness can be employed to enhance the protein binding sensing performance of the SPR-based sensor.

Author Contributions

Conceptualization, S.P. and S.S.; methodology, S.P.; software, S.S. and S.P.; validation: S.P., S.S. and T.T.; formal analysis, S.P.; investigation, S.P., S.S. and T.T.; resources, S.P.; data curation, S.S.; writing—original draft preparation, S.S.; writing—review and editing, S.P. and T.T.; visualization, S.S.; supervision, S.P. and T.T.; project administration, S.P. and T.T.; funding acquisition, S.P. and T.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Research Institute of Rangsit University (RSU) and the School of Engineering of King Mongkut’s Institute of Technology Ladkrabang (KMITL)—scholarship for Siratchakrit Shinnakerdchoke, funding number KREF016601.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to acknowledge the resources, fruitful discussion, and suggestions from the College of Biomedical Engineering, Rangsit University, Thailand, and the School of Engineering, KMITL, Thailand.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brongersma, M.L.; Kik, P.G. Surface Plasmon Nanophotonics; Springer: Dordrecht, The Netherlands, 2007; Volume 131. [Google Scholar]
  2. Somekh, M.G.; Pechprasarn, S. Surface plasmon, surface wave, and enhanced evanescent wave microscopy. In Handbook of Photonics for Biomedical Engineering; Springer: Dordrecht, The Netherlands, 2017; pp. 503–543. [Google Scholar]
  3. Homola, J.; Yee, S.S.; Gauglitz, G. Surface plasmon resonance sensors. Sens. Actuators B Chem. 1999, 54, 3–15. [Google Scholar] [CrossRef]
  4. Suda, Y.; Arano, A.; Fukui, Y.; Koshida, S.; Wakao, M.; Nishimura, T.; Kusumoto, S.; Sobel, M. Immobilization and clustering of structurally defined oligosaccharides for sugar chips: An improved method for surface plasmon resonance analysis of protein-carbohydrate interactions. Bioconjug. Chem. 2006, 17, 1125–1135. [Google Scholar] [CrossRef]
  5. Mariani, S.; Minunni, M. Surface plasmon resonance applications in clinical analysis. Anal. Bioanal. Chem. 2014, 406, 2303–2323. [Google Scholar] [CrossRef] [PubMed]
  6. Chen, J.; Park, B. Label-free screening of foodborne Salmonella using surface plasmon resonance imaging. Anal. Bioanal. Chem. 2018, 410, 5455–5464. [Google Scholar] [CrossRef] [PubMed]
  7. Douzi, B. Protein–protein interactions: Surface plasmon resonance. In Bacterial Protein Secretion Systems: Methods and Protocols; Humana Press: New York, NY, USA, 2017; pp. 257–275. [Google Scholar]
  8. Drescher, D.G.; Selvakumar, D.; Drescher, M.J. Analysis of protein interactions by surface plasmon resonance. Adv. Protein Chem. Struct. Biol. 2018, 110, 1–30. [Google Scholar]
  9. Sangworasil, M.; Pechprasarn, S.; Learkthanakhachon, S.; Ittipornnuson, K.; Suvarnaphaet, P.; Albutt, N. Investigation on feasibility of using surface plasmons resonance (SPR) sensor for ultrasonic detection: A novel optical detection of ultrasonic waves. In Proceedings of the 2016 9th Biomedical Engineering International Conference (BMEiCON), Laung Prabang, Laos, 7–9 December 2016; pp. 1–3. [Google Scholar]
  10. Shen, M.; Learkthanakhachon, S.; Pechprasarn, S.; Zhang, Y.; Somekh, M.G. Adjustable microscopic measurement of nanogap waveguide and plasmonic structures. Appl. Opt. 2018, 57, 3453–3462. [Google Scholar] [CrossRef] [PubMed]
  11. Suvarnaphaet, P.; Pechprasarn, S. Enhancement of long-range surface plasmon excitation, dynamic range and figure of merit using a dielectric resonant cavity. Sensors 2018, 18, 2757. [Google Scholar] [CrossRef] [Green Version]
  12. Suvarnaphaet, P.; Pechprasarn, S. Quantitative cross-platform performance comparison between different detection mechanisms in surface plasmon sensors for voltage sensing. Sensors 2018, 18, 3136. [Google Scholar] [CrossRef] [Green Version]
  13. De Mol, N.J.; Dekker, F.J.; Broutin, I.; Fischer, M.J.; Liskamp, R.M. Surface plasmon resonance thermodynamic and kinetic analysis as a strategic tool in drug design. Distinct ways for phosphopeptides to plug into Src-and Grb2 SH2 domains. J. Med. Chem. 2005, 48, 753–763. [Google Scholar] [CrossRef] [Green Version]
  14. Navratilova, I.; Hopkins, A.L. Emerging role of surface plasmon resonance in fragment-based drug discovery. Future Med. Chem. 2011, 3, 1809–1820. [Google Scholar] [CrossRef]
  15. Chow, T.W.; Pechprasarn, S.; Meng, J.; Somekh, M.G. Single shot embedded surface plasmon microscopy with vortex illumination. Opt. Express 2016, 24, 10797–10805. [Google Scholar] [CrossRef] [PubMed]
  16. Englebienne, P.; Van Hoonacker, A.; Verhas, M. Surface plasmon resonance: Principles, methods and applications in biomedical sciences. Spectroscopy 2003, 17, 255–273. [Google Scholar] [CrossRef]
  17. Kretschmann, E.; Raether, H. Radiative decay of non radiative surface plasmons excited by light. Z. Nat. A 1968, 23, 2135–2136. [Google Scholar] [CrossRef]
  18. Guo, J.; Keathley, P.D.; Hastings, J. Dual-mode surface-plasmon-resonance sensors using angular interrogation. Opt. Lett. 2008, 33, 512–514. [Google Scholar] [CrossRef] [PubMed]
  19. Zeng, Y.; Zhou, J.; Wang, X.; Cai, Z.; Shao, Y. Wavelength-scanning surface plasmon resonance microscopy: A novel tool for real time sensing of cell-substrate interactions. Biosens. Bioelectron. 2019, 145, 111717. [Google Scholar] [CrossRef] [PubMed]
  20. Ran, B.; Lipson, S. Comparison between sensitivities of phase and intensity detection in surface plasmon resonance. Opt. Express 2006, 14, 5641–5650. [Google Scholar] [CrossRef]
  21. Kabashin, A.V.; Patskovsky, S.; Grigorenko, A.N. Phase and amplitude sensitivities in surface plasmon resonance bio and chemical sensing. Opt. Express 2009, 17, 21191–21204. [Google Scholar] [CrossRef] [PubMed]
  22. Barer, R.; Tkaczyk, S. Refractive index of concentrated protein solutions. Nature 1954, 173, 821–822. [Google Scholar] [CrossRef]
  23. Crowell, J.; Ritchie, R. Surface-plasmon effect in the reflectance of a metal. JOSA 1970, 60, 794–799. [Google Scholar] [CrossRef]
  24. Braundmeier, A., Jr.; Arakawa, E. Effect of surface roughness on surface plasmon resonance absorption. J. Phys. Chem. Solids 1974, 35, 517–520. [Google Scholar] [CrossRef]
  25. Rahman, T.S.; Maradudin, A.A. Surface-plasmon dispersion relation in the presence of surface roughness. Phys. Rev. B 1980, 21, 2137. [Google Scholar] [CrossRef]
  26. Hoffmann, A.; Lenkefi, Z.; Szentirmay, Z. Effect of roughness on surface plasmon scattering in gold films. J. Phys. Condens. Matter 1998, 10, 5503. [Google Scholar] [CrossRef]
  27. Kurihara, K.; Suzuki, K. Theoretical understanding of an absorption-based surface plasmon resonance sensor based on Kretchmann’s theory. Anal. Chem. 2002, 74, 696–701. [Google Scholar] [CrossRef]
  28. Treebupachatsakul, T.; Shinnakerdchoke, S.; Pechprasarn, S. Analysis of effects of surface roughness on sensing performance of surface plasmon resonance detection for refractive index sensing application. Sensors 2021, 21, 6164. [Google Scholar] [CrossRef]
  29. Byun, K.M.; Yoon, S.J.; Kim, D. Effect of surface roughness on the extinction-based localized surface plasmon resonance biosensors. Appl. Opt. 2008, 47, 5886–5892. [Google Scholar] [CrossRef] [PubMed]
  30. Liu, H.; Wang, B.; Leong, E.S.; Yang, P.; Zong, Y.; Si, G.; Teng, J.; Maier, S.A. Enhanced surface plasmon resonance on a smooth silver film with a seed growth layer. ACS Nano 2010, 4, 3139–3146. [Google Scholar] [CrossRef]
  31. Schug, C.; Schempp, S.; Lamparter, P.; Steeb, S. Surface roughness of sputter-deposited gold films: A combined x-ray technique and AFM study. Surf. Interface Anal. 1999, 27, 670–677. [Google Scholar] [CrossRef]
  32. Zhang, J.; Irannejad, M.; Yavuz, M.; Cui, B. Gold nanohole array with sub-1 nm roughness by annealing for sensitivity enhancement of extraordinary optical transmission biosensor. Nanoscale Res. Lett. 2015, 10, 238. [Google Scholar] [CrossRef] [Green Version]
  33. Zhang, C.; Li, J.; Belianinov, A.; Ma, Z.; Renshaw, C.K.; Gelfand, R.M. Nanoaperture fabrication in ultra-smooth single-grain gold films with helium ion beam lithography. Nanotechnology 2020, 31, 465302. [Google Scholar] [CrossRef] [PubMed]
  34. Ng, D.K.; Bhola, B.S.; Bakker, R.M.; Ho, S.T. Ultrasmooth gold films via pulsed laser deposition. Adv. Funct. Mater. 2011, 21, 2587–2592. [Google Scholar] [CrossRef]
  35. Wieduwilt, T.; Kirsch, K.; Dellith, J.; Willsch, R.; Bartelt, H. Optical fiber micro-taper with circular symmetric gold coating for sensor applications based on surface plasmon resonance. Plasmonics 2013, 8, 545–554. [Google Scholar] [CrossRef]
  36. Diebel, J.; Löwe, H.; Samori, P.; Rabe, J. Fabrication of large-scale ultra-smooth metal surfaces by a replica technique. Appl. Phys. A 2001, 73, 273–279. [Google Scholar] [CrossRef]
  37. Miller, M.S.; Ferrato, M.-A.; Niec, A.; Biesinger, M.C.; Carmichael, T.B. Ultrasmooth gold surfaces prepared by chemical mechanical polishing for applications in nanoscience. Langmuir 2014, 30, 14171–14178. [Google Scholar] [CrossRef]
  38. Leandro, L.; Malureanu, R.; Rozlosnik, N.; Lavrinenko, A. Ultrathin, ultrasmooth gold layer on dielectrics without the use of additional metallic adhesion layers. ACS Appl. Mater. Interfaces 2015, 7, 5797–5802. [Google Scholar] [CrossRef] [Green Version]
  39. Gaylord, T.K.; Moharam, M. Analysis and applications of optical diffraction by gratings. Proc. IEEE 1985, 73, 894–937. [Google Scholar] [CrossRef]
  40. Moharam, M.; Grann, E.B.; Pommet, D.A.; Gaylord, T. Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings. J. Opt. Soc. Am. A 1995, 12, 1068–1076. [Google Scholar] [CrossRef]
  41. Thadson, K.; Visitsattapongse, S.; Pechprasarn, S. Deep learning-based single-shot phase retrieval algorithm for surface plasmon resonance microscope based refractive index sensing application. Sci. Rep. 2021, 11, 16289. [Google Scholar] [CrossRef] [PubMed]
  42. O’Donnell, K.; Mendez, E. Experimental study of scattering from characterized random surfaces. JOSA A 1987, 4, 1194–1205. [Google Scholar] [CrossRef] [Green Version]
  43. Johnson, P.B.; Christy, R.-W. Optical constants of the noble metals. Phys. Rev. B 1972, 6, 4370. [Google Scholar] [CrossRef]
  44. Davidson, A.M.; Brust, M.; Cooper, D.L.; Volk, M. Sensitive analysis of protein adsorption to colloidal gold by differential centrifugal sedimentation. Anal. Chem. 2017, 89, 6807–6814. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Saleem, M.R.; Ali, R. Polymer Resonant Waveguide Gratings. In Emerging Waveguide Technology; IntechOpen: London, UK, 2018. [Google Scholar]
  46. Nuutinen, T.; Karvinen, P.; Rahomäki, J.; Vahimaa, P. Resonant waveguide grating (RWG): Overcoming the problem of angular sensitivity by conical, broad-band illumination for fluorescence measurements. Anal. Methods 2013, 5, 281–284. [Google Scholar] [CrossRef]
  47. Mayet, A.S.; Cansizoglu, H.; Gao, Y.; Ghandiparsi, S.; Kaya, A.; Bartolo-Perez, C.; AlHalaili, B.; Yamada, T.; Devine, E.P.; Elrefaie, A.F. Surface passivation of silicon photonic devices with high surface-to-volume-ratio nanostructures. JOSA B 2018, 35, 1059–1065. [Google Scholar] [CrossRef]
  48. Zhang, C.; Zhou, Y.; Mi, L.; Ma, J.; Wu, X.; Fei, Y. High performance of a metal layer-assisted guided-mode resonance biosensor modulated by double-grating. Biosensors 2021, 11, 221. [Google Scholar] [CrossRef] [PubMed]
  49. Berguiga, L.; Zhang, S.; Argoul, F.; Elezgaray, J. High-resolution surface-plasmon imaging in air and in water: V (z) curve and operating conditions. Opt. Lett. 2007, 32, 509–511. [Google Scholar] [CrossRef]
  50. Moharam, M.; Gaylord, T. Rigorous coupled-wave analysis of planar-grating diffraction. JOSA 1981, 71, 811–818. [Google Scholar] [CrossRef]
  51. Chu, Y.; Schonbrun, E.; Yang, T.; Crozier, K.B. Experimental observation of narrow surface plasmon resonances in gold nanoparticle arrays. Appl. Phys. Lett. 2008, 93, 181108. [Google Scholar] [CrossRef]
  52. Zhang, J.; Pitter, M.C.; Liu, S.; See, C.; Somekh, M.G. Surface-plasmon microscopy with a two-piece solid immersion lens: Bright and dark fields. Appl. Opt. 2006, 45, 7977–7986. [Google Scholar] [CrossRef]
  53. Thadson, K.; Sasivimolkul, S.; Suvarnaphaet, P.; Visitsattapongse, S.; Pechprasarn, S. Measurement precision enhancement of surface plasmon resonance based angular scanning detection using deep learning. Sci. Rep. 2022, 12, 2052. [Google Scholar] [CrossRef]
  54. Meng, Q.-Q.; Zhao, X.; Lin, C.-Y.; Chen, S.-J.; Ding, Y.-C.; Chen, Z.-Y. Figure of merit enhancement of a surface plasmon resonance sensor using a low-refractive-index porous silica film. Sensors 2017, 17, 1846. [Google Scholar] [CrossRef] [Green Version]
  55. Treebupachatsakul, T.; Boosamalee, A.; Chaithatwanitch, K.; Pechprasarn, S. Generalized figure of merit for plasmonic dip measurement-based surface plasmon resonance sensors. Biomed. Opt. Express 2022, 13, 1784–1800. [Google Scholar] [CrossRef]
  56. Yang, Z.; Liu, C.; Gao, Y.; Wang, J.; Yang, W. Influence of surface roughness on surface plasmon resonance phenomenon of gold film. Chin. Opt. Lett. 2016, 14, 042401. [Google Scholar] [CrossRef] [Green Version]
  57. Zhao, M.; Wang, J.; Zhang, Y.; Ge, M.; Zhang, P.; Shen, J.; Li, C. Self-referenced refractive index sensor based on double-dips method with bimetal-dielectric and double-groove grating. Opt. Express 2022, 30, 8376–8390. [Google Scholar] [CrossRef]
  58. Agarwal, S.; Prajapati, Y.; Singh, V. Influence of metal roughness on SPR sensor performance. Opt. Commun. 2017, 383, 113–118. [Google Scholar] [CrossRef]
  59. Pechprasarn, S.; Chow, T.W.; Somekh, M.G. Application of confocal surface wave microscope to self-calibrated attenuation coefficient measurement by Goos-Hänchen phase shift modulation. Sci. Rep. 2018, 8, 8547. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Wy, Y.; Jung, H.; Hong, J.W.; Han, S.W. Exploiting plasmonic hot spots in Au-based nanostructures for sensing and photocatalysis. Acc. Chem. Res. 2022, 55, 831–843. [Google Scholar] [CrossRef]
  61. Sasivimolkul, S.; Pechprasarn, S.; Somekh, M.G. Analysis of open grating-based Fabry–Pérot resonance structures with potential applications for ultrasensitive refractive index sensing. IEEE Sens. J. 2021, 21, 10628–10636. [Google Scholar] [CrossRef]
Figure 1. (a) The Kretschmann configuration consisting of a water-based sensing region, BSA binding layer, noble metallic film, a glass substrate, and incident light, and (b) the SPR sensor’s outputs consisting of reflectance curves, measuring from the water-based environment sensing region shown in the blue curve and a 5 nm-thick BSA coating layer in water shown in the red curve.
Figure 1. (a) The Kretschmann configuration consisting of a water-based sensing region, BSA binding layer, noble metallic film, a glass substrate, and incident light, and (b) the SPR sensor’s outputs consisting of reflectance curves, measuring from the water-based environment sensing region shown in the blue curve and a 5 nm-thick BSA coating layer in water shown in the red curve.
Sensors 23 03377 g001
Figure 2. The Gaussian random rough surface for h 20 nm of and cl of 5 nm for (a) no BSA protein coating case (bare gold case), (b) BSA protein without sidewall coating, and (c) BSA protein with sidewall coating. The figures on the right show the zoomed-in version of rough surfaces for the three models.
Figure 2. The Gaussian random rough surface for h 20 nm of and cl of 5 nm for (a) no BSA protein coating case (bare gold case), (b) BSA protein without sidewall coating, and (c) BSA protein with sidewall coating. The figures on the right show the zoomed-in version of rough surfaces for the three models.
Sensors 23 03377 g002
Figure 3. Simulated rough SPR sensor based on the Kretschmann configuration.
Figure 3. Simulated rough SPR sensor based on the Kretschmann configuration.
Sensors 23 03377 g003
Figure 4. The flowchart of the overall simulation processes, including the rigorous coupled-wave analysis simulation, Monte Carlo simulation, and the performance parameters calculation.
Figure 4. The flowchart of the overall simulation processes, including the rigorous coupled-wave analysis simulation, Monte Carlo simulation, and the performance parameters calculation.
Sensors 23 03377 g004
Figure 5. The calculation procedures of the discussed quantitative performance parameters for the ideal uniform gold surface sensor, including (a) ΔI, FWHM, (b) S, and Isp.
Figure 5. The calculation procedures of the discussed quantitative performance parameters for the ideal uniform gold surface sensor, including (a) ΔI, FWHM, (b) S, and Isp.
Sensors 23 03377 g005
Figure 6. (a) The convergence test for the four extreme cases of the rough surfaces constructed with h of 2 nm and cl of 5 nm (solid blue line), h of 2 nm and cl of 50 nm (dashed blue line), h of 20 nm and cl of 5 nm (solid red line), and h of 20 nm and cl of 50 nm (dashed red line), tested with the number of diffractions orders up to 171, and (b) The convergence test on the simulated columns for the rough surfaces at the same four extreme cases, illustrated as average plasmonic intensities of the no BSA model, the nonsidewall BSA model, and the sidewall BSA model.
Figure 6. (a) The convergence test for the four extreme cases of the rough surfaces constructed with h of 2 nm and cl of 5 nm (solid blue line), h of 2 nm and cl of 50 nm (dashed blue line), h of 20 nm and cl of 5 nm (solid red line), and h of 20 nm and cl of 50 nm (dashed red line), tested with the number of diffractions orders up to 171, and (b) The convergence test on the simulated columns for the rough surfaces at the same four extreme cases, illustrated as average plasmonic intensities of the no BSA model, the nonsidewall BSA model, and the sidewall BSA model.
Sensors 23 03377 g006
Figure 7. SPR reflectance spectra simulated using the proposed procedure at different roughness parameters reported by Yang et al. [56].
Figure 7. SPR reflectance spectra simulated using the proposed procedure at different roughness parameters reported by Yang et al. [56].
Sensors 23 03377 g007
Figure 8. S of the (a) nonsidewall BSA model, (b) sidewall BSA model at different roughness levels in rad/µm2, and (c) the increase of the protein density due to sidewalls in percent as the surface roughness increased.
Figure 8. S of the (a) nonsidewall BSA model, (b) sidewall BSA model at different roughness levels in rad/µm2, and (c) the increase of the protein density due to sidewalls in percent as the surface roughness increased.
Sensors 23 03377 g008
Figure 9. (a) the reflectance spectra with a negative sensitivity at the operating point ‘a’, (b) degraded-sensitivity reflectance curves at the operating point ‘b’, (c) enhanced-sensitivity reflectance spectra at the operating point ‘c’, (d) highly enhanced-sensitivity reflectance spectra at the operating point ‘d’, and (e) the reflectance curves at the highly rough surface where the sensitivity cannot be computed at the operating point ‘e’. Note that the dotted curves indicate the plasmonic angles of the ideal uniform gold cases and the black arrows indicate the plasmonic angle shift direction.
Figure 9. (a) the reflectance spectra with a negative sensitivity at the operating point ‘a’, (b) degraded-sensitivity reflectance curves at the operating point ‘b’, (c) enhanced-sensitivity reflectance spectra at the operating point ‘c’, (d) highly enhanced-sensitivity reflectance spectra at the operating point ‘d’, and (e) the reflectance curves at the highly rough surface where the sensitivity cannot be computed at the operating point ‘e’. Note that the dotted curves indicate the plasmonic angles of the ideal uniform gold cases and the black arrows indicate the plasmonic angle shift direction.
Sensors 23 03377 g009
Figure 10. (a) the average plasmonic dip location, n0sinθsp, for the three models at different roughness levels, and (b) the reflectance spectra from the nonprotein coated model at h of 0 nm (ideal smooth), 5 nm, 10 nm, and 15 nm, and the equivalent cl of 20 nm. (c) The FWHM of the SPR detection at different degrees of roughness and (d) reflectance spectra for the rough gold surface of h and cl of 18 nm and 40 nm, respectively, compared to the uniform gold surface. Note that the black arrows in (b) indicate the plasmonic angle shift direction when the h increased.
Figure 10. (a) the average plasmonic dip location, n0sinθsp, for the three models at different roughness levels, and (b) the reflectance spectra from the nonprotein coated model at h of 0 nm (ideal smooth), 5 nm, 10 nm, and 15 nm, and the equivalent cl of 20 nm. (c) The FWHM of the SPR detection at different degrees of roughness and (d) reflectance spectra for the rough gold surface of h and cl of 18 nm and 40 nm, respectively, compared to the uniform gold surface. Note that the black arrows in (b) indicate the plasmonic angle shift direction when the h increased.
Sensors 23 03377 g010
Figure 11. (a) the average ΔI and (b) the average Isp of both the nonsidewall BSA and the sidewall BSA models.
Figure 11. (a) the average ΔI and (b) the average Isp of both the nonsidewall BSA and the sidewall BSA models.
Sensors 23 03377 g011
Figure 12. (a) the FoM1 of the BSA without the sidewall, (b) BSA with the sidewall model; (c) FoM2 for SPR detection [55] obtained from the BSA without the sidewall, and (d) BSA with the sidewall model.
Figure 12. (a) the FoM1 of the BSA without the sidewall, (b) BSA with the sidewall model; (c) FoM2 for SPR detection [55] obtained from the BSA without the sidewall, and (d) BSA with the sidewall model.
Sensors 23 03377 g012
Figure 13. The Δsinθsp for different deposited thicknesses of BSA from 0 nm (bare gold) to 10 nm for ideal smooth gold surface, rough gold surface with h of 9 nm, and cl of 8 nm for the no sidewall and sidewall coating cases.
Figure 13. The Δsinθsp for different deposited thicknesses of BSA from 0 nm (bare gold) to 10 nm for ideal smooth gold surface, rough gold surface with h of 9 nm, and cl of 8 nm for the no sidewall and sidewall coating cases.
Sensors 23 03377 g013
Figure 14. (a) Intensity detection and (b) phase detection of the operating point ‘a’, ‘b’, ‘c’, ‘d’, and the theoretically uniform surface SPR-based sensor.
Figure 14. (a) Intensity detection and (b) phase detection of the operating point ‘a’, ‘b’, ‘c’, ‘d’, and the theoretically uniform surface SPR-based sensor.
Sensors 23 03377 g014
Figure 15. The plots of the total electric field intensity in the SI unit in the x-direction and z-direction (|Ex|2 + |Ez|2) on the gold-based SPR sensor with (a) an ideal smooth surface; (b) gratings structure with a gratings’ height of 24 nm, a gratings’ period of 500 nm, and a fill factor of 0.5; and (c) a rough sensing surface with h and cl of 9 nm and 8 nm, respectively (operating point ‘d’).
Figure 15. The plots of the total electric field intensity in the SI unit in the x-direction and z-direction (|Ex|2 + |Ez|2) on the gold-based SPR sensor with (a) an ideal smooth surface; (b) gratings structure with a gratings’ height of 24 nm, a gratings’ period of 500 nm, and a fill factor of 0.5; and (c) a rough sensing surface with h and cl of 9 nm and 8 nm, respectively (operating point ‘d’).
Sensors 23 03377 g015
Table 1. The root mean square (RMS) roughness of metal film after different surface treatments.
Table 1. The root mean square (RMS) roughness of metal film after different surface treatments.
Surface Treatment TechniquesRMS Roughness
Sputter Coating (No Treatment) [31]1.4–2.5 nm
Thermal Annealing [32]<1 nm
Chemically Grown Single-Crystalline Gold [35]<1 nm
Chemical Polishing [37,38]0.38 ± 0.05 nm
Helium Ion Beam [33]0.267 nm
Mica Substrate Utilizing [36]0.2 nm
Laser Ablation [34]0.17 nm
Table 2. Shows the roughness parameters and plasmonic angles from the experimental results reported by Yang et al. [56] and the plasmonic angles from Figure 7.
Table 2. Shows the roughness parameters and plasmonic angles from the experimental results reported by Yang et al. [56] and the plasmonic angles from Figure 7.
Thermal Annealing Temoerature in °CRMS of Roughness (nm)Roughness Height (nm)Experimental θsp (Degrees)Simulated θsp Using the Proposed Method
(Degrees)
100 °C1.2606.144.0444.04
200 °C0.9063.743.9443.93
300 °C0.7002.943.7543.70
400 °C0.4151.543.6543.56
Table 3. The quantitative performance parameters comparison between the surface plasmon resonance detection employed with a uniform and rough gold layer.
Table 3. The quantitative performance parameters comparison between the surface plasmon resonance detection employed with a uniform and rough gold layer.
Quantitative
Performance Parameters
Ideal Uniform Surfaceh = 3 nm,
cl = 10 nm
(Operating Point ‘a’)
h = 10 nm,
cl = 30 nm
(Operating Point ‘b’)
h = 12 nm,
cl = 20 nm
(Operating Point ‘c’)
h = 9 nm,
cl = 8 nm
(Operating Point ‘d’)
h = 16 nm,
cl = 15 nm
(Operating Point ‘e’)
No Sidewall BSASidewall BSANo Sidewall BSASidewall BSANo Sidewall BSASidewall BSANo Sidewall BSASidewall BSANo Sidewall BSASidewall BSA
S
(rad/µm2)
122.59−379.88−379.8885.3185.31996.111280.001191.131608.02--
n0sinθsp (unitless)1.4401.4811.4811.5061.5061.4951.4961.4861.4901.466 *1.473 *
FWHM
(rad∙RIU
/µm)
0.410.580.580.790.790.660.650.670.670.770.75
ΔI
(unitless)
0.900.880.880.470.470.130.150.100.140.0010.003
Isp
(unitless)
0.0070.0180.0180.410.410.720.700.690.600.750.69
FoM1
(unitless)
299.00−654.97−654.97107.98107.981509.261969.231777.802400.02--
FoM2
(unitless)
980.66−1677.42−1677.4292.5292.52590.75833.81616.841020.33--
* The n0sinθsp of the operating point ‘e’ can still be determined. However, the reflectance spectrum of the model without a protein coating lost its plasmonic dip structure, resulting in an unidentified S and FoM.
Table 4. Shows the enhancement ratio of S and FoMs compared to an ideal uniform gold case and their effect in percentage due to the additional protein and wave effect from the surface roughness.
Table 4. Shows the enhancement ratio of S and FoMs compared to an ideal uniform gold case and their effect in percentage due to the additional protein and wave effect from the surface roughness.
h = 3 nm,
cl = 10 nm
(Operating Point ‘a’)
h = 10 nm,
cl = 30 nm
(Operating Point ‘b’)
h = 12 nm,
cl = 20 nm
(Operating Point ‘c’)
h = 9 nm,
cl = 8 nm
(Operating Point ‘d’)
ERS,roughness−3.100.7010.4413.12
ERS,LSPR−3.10 (100%)0.70 (100%)8.13 (77.87%)9.72 (74.07%)
ERS,protein0 (0%)0 (0%)2.31 (22.13%)3.40 (25.93%)
ERFoM1,roughness−2.190.366.598.03
ERFoM1,LSPR−2.19 (100%)0.36 (100%)5.05 (76.64%)5.95 (74.07%)
ERFoM1,protein0 (0%)0 (0%)1.54 (23.36%)2.08 (25.93%)
ERFoM2,roughness−1.710.090.851.04
ERFoM2,LSPR−1.71 (100%)0.09 (100%)0.60 (70.85%)0.63 (60.45%)
ERFoM2,protein0 (0%)0 (0%)0.25 (29.15%)0.41 (39.55%)
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

Treebupachatsakul, T.; Shinnakerdchoke, S.; Pechprasarn, S. Sensing Mechanisms of Rough Plasmonic Surfaces for Protein Binding of Surface Plasmon Resonance Detection. Sensors 2023, 23, 3377. https://doi.org/10.3390/s23073377

AMA Style

Treebupachatsakul T, Shinnakerdchoke S, Pechprasarn S. Sensing Mechanisms of Rough Plasmonic Surfaces for Protein Binding of Surface Plasmon Resonance Detection. Sensors. 2023; 23(7):3377. https://doi.org/10.3390/s23073377

Chicago/Turabian Style

Treebupachatsakul, Treesukon, Siratchakrit Shinnakerdchoke, and Suejit Pechprasarn. 2023. "Sensing Mechanisms of Rough Plasmonic Surfaces for Protein Binding of Surface Plasmon Resonance Detection" Sensors 23, no. 7: 3377. https://doi.org/10.3390/s23073377

APA Style

Treebupachatsakul, T., Shinnakerdchoke, S., & Pechprasarn, S. (2023). Sensing Mechanisms of Rough Plasmonic Surfaces for Protein Binding of Surface Plasmon Resonance Detection. Sensors, 23(7), 3377. https://doi.org/10.3390/s23073377

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