Open Access Article
Timo
Küster
and
Geoffrey D.
Bothun
*
Department of Chemical Engineering, University of Rhode Island, 2 East Alumni Ave, Kingston, RI 02881, USA. E-mail: gbothun@uri.edu; Tel: +1-401-874-9518
First published on 12th June 2021
The accurate and fast measurement of nitrate in seawater is important for monitoring and controlling water quality to prevent ecologic and economic disasters. In this work we show that the in situ detection of nitrate in aqueous solution is feasible at nanomolar concentrations through surface enhanced Raman spectroscopy (SERS) using native nanostructured gold substrates without surface functionalization. Spectra were analyzed as collected or after standard normal variate (SNV) normalization, which was shown through Principal Component Analysis (PCA) to reduce spectral variations between sample sets and improve Langmuir adsorption model fits. An additional normalization approach based on the substrate silicon template showed that silicon provided an internal standard that accounted for the spectral variance without the need for SNV normalization. Nitrate adsorption was well-described by the Langmuir adsorption model, consistent with an adsorbed monolayer, and a limit of detection of 64 nM nitrate was obtained in ultrapure water, representing environmentally relevant concentrations. Free energy calculations based on the Langmuir adsorption constants, approximating equilibrium adsorption constants, and calculated self-energy arising from image charge, accounting for electrostatic interactions with a polarizable nanostructured substrate, suggest that nitrate adsorption was partially driven by an entropy gain presumably due to dehydration of the gold substrate and/or nitrate ion. This work is being extended to determine if similar statistical and normalization methods can be applied to nitrate detection in complex natural waters where non-target ions and molecules are expected to interfere.
If nitrate and other nutrients reach excessive levels they can contribute to eutrophication, which is a process where rapid plant growth diminishes the supply of dissolved oxygen necessary to support higher trophic marine life. In addition to creating “dead zones” for respiring species, eutrophication can promote algae growth and the formation of large algae blooms that release toxins such as domoic acid into the surrounding waters.6,7 The negative economic impacts of eutrophication tied to fishing, aquaculture, and tourism can be significant, as can be the adverse health effects of individuals who are exposed to harmful algae blooms (HABs). Affordable, deployable, and accurate tools to improving the spatiotemporal detection of inorganic nitrogen in situ will increase ecological monitoring and inform computational approaches that are being developed to predict emerging conditions that may lead to the formation of HABs. Early warning of these devastating ecological events will provide additional time to alert coastal stakeholders and enact countermeasures.
Many methods exist for the detection of nitrate in fresh and wastewater.8 However, the continuous in-field detection of nitrate in seawater is less explored. The most commonly used method for nitrate detection in seawater is UV-vis spectroscopy with, for example, commercially available instruments providing limits of detection as low as 0.5 μM.9
Surface Enhanced Raman Spectroscopy (SERS) provides an ultrasensitive platform for sensor design and has been described as a molecular fingerprinting technique capable of single molecular detection when ordered metallic nanostructured substrates are employed.10 For environmental applications, SERS measurements can be conducted with lasers in the near infrared regime where there is little interference with water11 and portable handheld Raman instruments enable in situ field measurements. One of the challenges to SERS detection is the dependency of the signal intensity on the distance from the surface (z), proportional to z−12, requiring analyte molecules to be within approximately z ≤ 4 nm for the SERS effect to be observable.11 Because of this requirement, SERS measurements are commonly taken from analytes dried from solution on a substrate. In an in situ solution phase measurement, the diffusive transport of a target analyte to the fixed surface and the adsorption affinity of the analyte are limiting factors, potentially reducing the signal strength and responsiveness of the sensor.11,12
Nitrate detection through normal Raman spectroscopy and SERS has been demonstrated by others.13–21 With advances in sensing equipment and nano-fabrication, efforts have been made to enhance detection capabilities for nutrient pollutants. Early work in the SERS field explored colloidal gold nanoparticle solutions, etched wafers, and gold sputtered nanoparticles where, for example, ions are enriched due to charge attraction and hotspots form via aggregation.11,12,22 A fixed two-dimensional substrate provides a more practical platform for continuous environmental sensing than dispersed nanoparticles.23In situ nitrate measurements were shown to be feasible at micromolar concentrations with functionalized thiol based self-assembled monolayers (SAMs) formed on commercially available substrates (LOD = 8.06 μM on Klarite™; 10 μM on Silmeco SERStrate Au)24,25 or with gold nanoparticles.14,26 Other approaches include the usage of reporter molecules on the substrate or particle surface, such as immobilized Griess reagents to form azo dyes in the presence of nitrite, which have a strong and specific Raman signal.16,27 For the detection of nitrate through this approach, where LODs between 1.52 μM (ref. 16) and 30.7 μM (ref. 27) nitrite have been reported, a reduction step from nitrate to nitrite is required and the diazotization is not easily reversible, limiting the feasibility for continuous in field measurement applications.28
An additional challenge in SERS detection is the need to overcome the diffusion barrier between solution and substrate. The application of electric fields and charged substrates has been tested to force charged analyte molecules closer to the SERS substrate.28 Surface functionalization with cationic ligands has been shown to increase anion adsorption on SERS sensors;24,28,29 however, this approach may not be feasible for field detection as the ligands are not ion selective and will attract other negatively charged organic and inorganic molecules that will interfere with the measurements. Chemometrics is a tool commonly used in the SERS community. Inclusion of advanced normalization methods and investigation of SER spectra for hidden trends through principal component analysis and machine learning approaches shows to be promising to extract formerly hidden trends from the obtained data.30,31
The selection of currently available techniques for in situ nitrate detection clearly shows the need for a modified measurement approach. Hence, we focused on non-functionalized substrates using principal component analysis, standard normal variate spectral normalization and internal silicon standards to reduce background noise and remove signal bias due to spectral features inherent to the substrates. The objective of this work is two-fold: (1) to demonstrate low level in situ SERS detection of nitrate in aqueous solution through statistical and internal normalization using non-functionalized nanostructured gold substrates and (2) to employ adsorption models and image charge theory to gain additional fundamental insight on nitrate detection.
An initial water contact angle of 129° was measured on dry, as-received substrates indicative of a hydrophobic surface. After 60 min the water droplet wetted the substrate surface yielding a contact angle of 66°. When these wetted substrates were dried, and the water contact angle was again measured, there was minimal change in the initial contact angle. This observation shows that there was consistent wetting behavior of the substrates after initial water exposure. Therefore, all in situ SERS measurements were conducted after initial water exposure.
Raman spectra were measured for as-received (dry) SERS substrates and SERS substrates in deionized water using a 3D printed ‘beaker’ with an insert that immersed the substrate (Fig. 1D and E). Three common peaks are labeled with two peaks at approximate wavenumbers of 930 cm−1 and 1170 cm−1. The peaks centered near 930 cm−1 and 1170 cm−1 are likely due to hydrocarbons, possibly short aliphatic chains, that adsorbed from the atmosphere onto the gold surface. Additional unlabeled peaks over a wavenumber range of approximately 1280 cm−1 to 1400 cm−1 are consistent with CO2 and were also observed in solution Raman spectroscopy (Fig. S1, ESI†).36 The third peak at 520 cm−1 corresponds to crystalline silicon from the SERS substrate. The position and intensity of this peak when immersed in water was consistent across the substrates, providing a potential internal standard for nitrate detection. Patze et al.39 have used the silicon peak at 521 cm−1 arising from a similar substrate as an internal standard to improve the detection of the antibiotic sulfamethoxazole.
![]() | (1) |
n(z) = no e−βw0(z)qi2 | (2) |
The calculated distribution of monovalent ions near a planar gold surface are shown in Fig. 2. At a distance of 0.2 nm from the surface, which is slightly larger than the radii of an ion based on 0.2λB and approximates contact between the ion and the surface, there is nearly a 2.5-fold increase in the ion concentration compared to the bulk. The ion concentration decreases exponentially from the surface and is close to the bulk concentration at a distance of 10 nm. This analysis suggests that the region from approximately 0.2 to 3 nm can be considered an enrichment zone for SERS detection. Ion enrichment factors near the SERS surface are independent of bulk ion concentration. As a result, in situ SERS detection of NO3− is expected to be a direct function of the bulk NO3− concentration.
The origin of the peak at 1332 cm−1 is less clear. Gajaraj et al. assign dried nitrite (NO2−) to a wavenumber of 1326 cm−1.25 Elsewhere, asymmetric stretching (ν3), which is a weak Raman band compared to symmetric stretching, of NO3− has been reported at 1345 cm−1 for dissociated NaNO3 in solution.41 However, we observed peaks in the same region for dissociated NaNO3 and water in the solution Raman spectra, which suggests this may be related to dissolved gas. In the SERS spectra (Fig. 3A and B) the peak at 1332 cm−1 is neighbored by a peak with a maximum near 1380 cm−1 that is present across the range of [NaNO3] examined. Based on this analysis, the 1332 cm−1 peak is assigned to ν3 NO3− and the 1380 cm−1 peak is attributed to water or dissolved gases also observed in the solution Raman spectra.
While the NO3− peaks identified at 1079 cm−1 and 1332 cm−1 are clearly observed the differences in peak intensity with increasing [NaNO3], particularly at low concentrations, are subtle. Principal component analysis (PCA) was conducted for time-averaged spectra at each [NaNO3] concentration over the wavenumber range from 1000 cm−1 to 1350 cm−1 to verify that these differences are statistically significant. Principal component 1 (PC1) accounted for 81.5% and 84.9% of the variance for baselined and normalized spectra, respectively. With SNV normalization, PC1 accounted for more spectral variance indicating that this method can improve NO3− detection by accounting for spectral noise between measurements. Fig. 3C and D show that each [NaNO3] concentration has a distinct PC1 value, with PC1 decreasing (baselined) or increasing (normalized) with [NaNO3].
The nature of our measurement setup allows us to observe sorption processes in real time. We show this by plotting the time-averaged intensity and normalized intensity of SERS peaks assigned to NO3− at 1079 cm−1 (ν1; Fig. 4A) and 1332 cm−1 (ν3; Fig. 4B) as a function of nitrate concentration. Data were analyzed with the Langmuir adsorption model assuming that the intensity or normalized intensity, I, were directly proportional to the concentration of NO3− bound to the surface.
![]() | (3) |
In eqn (3)Isat is the saturated intensity or normalized intensity, C is the bulk [NO3−], and K is the Langmuir constant defined as
![]() | (4) |
Based on the calculated ion distribution, this places NO3− at the gold surface with near 2.5-fold increase in concentration relative to the bulk solution. The agreement with the Langmuir model suggests the presence of a monolayer, which is supported by the image charge theory discussed above. The K values associated with the model fit, ∼0.011 nM−1 based on intensity and ∼0.007 nM−1 based on normalized intensity, were the same for each peak (1079 cm−1 and 1332 cm−1), confirming that both peaks were specific to NO3− in the SER spectra. The goodness of the fit (R2) reported in Fig. 4 shows a higher agreement of the fitting function with the SNV normalized (R1079, SNV2 = 0.94 and R1332, SNV2 = 0.93) than for the baselined data (R1079, BL2 = 0.74 and R1332, BL2 = 0.81), confirming our previous finding that SNV normalization was successfully applied to reduce spectral variability. Fitting parameters for the Langmuir model are summarized in Table S1.†
The correlation found in the data allows for the calculation of a limit of detection, by estimating the intensity at the LOD, Ib, from the mean
and standard deviation σb of the blank43 according to eqn (5).
![]() | (5) |
Following procedures described in the literature,44,45 the calculated intensity can be equated to a concentration C, through rearrangement of the Langmuir isotherm as shown in eqn (6).
![]() | (6) |
A low limit of detection of 64 nM was determined for nitrate in MilliQ water despite the measurements being conducted in situ with as received SERS substrates where there was no specific affinity for the target analyte.
Assuming that the Langmuir adsorption constants, K (Fig. 4), are equal to the equilibrium adsorption constants, the resulting adsorption free energy based on ΔG = −RT
ln(K) was approximately −4 kJ mol−1. The adsorption entropy was determined from ΔG = ΔH − TΔS assuming that the adsorption enthalpy, ΔH, can be estimated as the self-energy arising from w0(z) at z = 0.2 nm (−2.2 kJ mol−1). The entropy term TΔS was 1.8 kJ mol−1 K−1 denoting a gain in entropy upon nitrate adsorption presumably due to the dehydration of the NO3− (with a hydration number = 4; Table 1) and/or gold substrate surface.
The Langmuir adsorption model was also used to fit the ratio of intensities at 1079 cm−1 and 1332 cm−1 to the intensity at 520 cm−1 to determine if the silicon peak is a suitable internal standard (Fig. 5). Intensity or normalized intensity ratios as a function of [NO3−1] were superimposable and the Langmuir constants, K, for model fits at 1079 cm−1 and 1332 cm−1 were 0.005 nM−1 for both peaks. These values are closer to the Langmuir constants K obtained from SNV normalized data (0.007 nM−1 for both cases) than to those of the baselined data (0.011 nM−1 for both cases), showing that both normalization methods yield reasonable results.
SERS measurements were conducted with a laser intensity of 100 mW under orbital raster scanning with an integration time of 5 s. A fresh substrate was used for every concentration series. Before the experiment was started the laser was focused by continuously measuring the Raman signal with an integration time of 0.2 s while varying the distance between substrate and laser lens until the highest signal was obtained (approximately 12 mm between laser and substrate).
Measurements were conducted with increasing concentrations of NaNO3 in ultrapure water, ranging from 0 to 2213 nM. A freshly 3D-printed beaker was used for each concentration series, washed with ultrapure water and dry blown with compressed nitrogen before usage. During the measurements, the instrument was covered with a housing that blocked ambient light. The background signal of the as-received, dry substrates was measured in triplicate before the beaker was filled with measurement solution. During all measurements background light subtraction was activated. Concentration series measurements were conducted by filling the beaker with MilliQ water (liquid height = 3.5 mm above substrate) and refocusing the laser to account for the focal change introduced by the medium. Measurements were started within 30 s of solvent exposure. Each concentration was measured every 2 min for up to one hour. Concentrations were increased by emptying the beaker and rinsing it with three times the beaker volume of MilliQ water, steadily dispensed from a pipette with low pressure. After the cleaning procedure, the beaker was refilled with solution of the next higher concentration. The measurement procedure was repeated until the highest concentrated solution was reached.
![]() | (7) |
) and I(
) are the SNV modified and unmodified peak intensities, respectively, at a given wavenumber
,
is the average intensity of the spectrum, and σ is the standard deviation of the spectrum. The method produced negative intensity values of the baseline on the order of −0.4 a.u. for all spectra. To allow for the comparison of normalized intensities by a common starting point, and to allow for easier fitting the baseline for all spectra was manually shifted to zero, so that no negative values remained in the spectra.
Principal Component Analysis (PCA) was applied to all collected spectra after baselining with and without the application of SNV. To conduct PCA, spectra were grouped by concentration and analyzed with the OriginPro application “PCA for Spectroscopy” (Version 2019b, OriginLab Corporation, Northampton MA, USA).
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/d1na00156f |
| This journal is © The Royal Society of Chemistry 2021 |