The relationship between Raman signal intensity, SERS-BG signal intensity and the SERS-BG integral ƩI BG . For our experiments we utilised a commercial SERS substrate called Mato (Integrated Optics), which was covered with a 100 µmolar ATP solution. As the overall integral of a SERS spectrum is dominated by the SERS-BG IBG and only a negligible contribution is due to the narrow Raman signals, we refer in the following to the overall integral SERS spectrum as summed SERS-BG signal ƩIBG.
In Fig. 1a SERS spectra of different enhancement are shown. The intensities of the three most intense Raman signals of ATP are plotted against ƩIBG for a total of 504 spectra. The signals originate from the adenine ring breathing vibration at 731 cm− 1, the PO32−-stretching vibration at 1006 cm− 1 and the C5-N7-stretching vibration of the adenine ring at 1334 cm− 1 (for detailed spectral assignment see Figure S1 and Table S1). We observe strong linear relationships with high coefficients of determination (R2) of 0.94, 0.94 and 0.99, respectively. However, at high intensities spectra tend to deviate from the linear model, most noticeably for the signal at 731 cm− 1. A comparison of the R² values of wavenumber-wise linear models of the SERS spectra signal intensities vs ƩIBG (Fig. S2) reveals the influence of the Raman signals. At positions of Raman signals from ATP the R2 values decrease in correspondence to the analyte signals and follow the band shape, proving that the decreased linear relationship is indeed due to variations of the pronounced Raman signals. At positions away from the Raman signals, the intensity originates only from SERS-BG and correlates highly with ƩIBG.
We therefore separated the spectral information into the narrow Raman signals and an underlying smooth baseline by an asymmetric least squares baseline fitting. Since the baseline is assumed to be mainly due to the SERS-BG signal, we use it in this work as an estimation for the SERS-BG. As expected, the correlation between the SERS-BG and ƩIBG is very high (Fig. 1b). These results show that ƩIBG is a very good predictor for the intensity of SERS-BG signals in a SERS spectrum at individual relative wavenumbers.
However, much more relevant in the context of the SERS-BG imaging method introduced here is the relationship between ƩIBG and the intensity of the Raman signals, which both are LSPR modulated processes. For the Raman signals at 731 cm− 1 and (to a lesser extend) at 1006 cm− 1 the determined values differ from the linear model fit of signal intensity vs ƩIBG. This might originate from slightly different signal generation processes for SERS and SERS-BG (see introduction) or the analyte/substrate interaction. While both the SERS and SERS-BG signal intensities highly depend on the plasmonic properties of the hot spots (e.g. due to size and shape of the nanostructure), the effective number of analyte molecules and their orientation within the hot spots might have a more pronounced impact just on the SERS signal. Therefore, a high correlation of intense SERS signals with ƩIBG is not imperative. And indeed, the correlation of the Raman signals at 731 cm− 1 and 1006 cm− 1 with ƩIBG is relatively low (Fig. 1c). Nevertheless, it is still high enough to enable the differentiation of hot spot regions with high Raman signal enhancement from regions of low Raman signal enhancement by utilizing ƩIBG. The big advantage of using ƩIBG instead of a single Raman signal is its higher signal intensity by a factor of ~ 103. This allows an easy transfer to fast imaging processes by detecting ƩIBG as an integrated signal intensity along the whole SERS spectrum. This approach forms the basis for our proposed SERS-BG imaging based hot spot localisation method.
Principle of SERS-BG based SERS-Analysis. Finding the hot spots on a SERS substrate is usually difficult, because they are typically very rare and occupy just a fraction (less than 1%) of the surface10. Additionally, for high signal enhancement it is crucial to measure at the exact hot spot positions, not just in the close vicinity. For example, all 504 Raman spectra acquired in our study were taken directly at or in the vicinity of hot spots, but only 3% of them exhibited strong enhancement. Conventional SERS experiments (schematically depicted in Fig. 2a), where a laser spot is randomly focused on different positions of the SERS substrate surface until a position of high signal enhancement is found, are therefore quite challenging and time consuming procedures. Moreover, due to the rareness of hot spots, there is a high probability not to find any of the “hottest” hot spots and thus to not fully exploit the signal enhancement capabilities of the SERS substrate. This results last but not least in a decreased reproducibility and reliability. A procedure often applied to circumvent this problem, especially when dealing with heterogeneous substrates, is to use objectives with a low numerical aperture. This approach allows to focus on a large area of the substrate and thus obtain an averaged, more reproducible SERS signal and a higher probability to hit a hot spot, but it reduces the total signal intensity due to the decreased collection efficiency of low NA objectives.
By contrast, our SERS-BG imaging based method (schematically depicted in Fig. 2b) aims to acquire SERS spectra only of regions of highest signal enhancement with high spatial resolution. This is achieved by laser-scanning the red-shifted light emission IImg of the SERS substrate using fast and sensitive HyD detectors (HYbrid Detectors of photomultiplier tubes and avalanche diodes) before the actual SERS spectral acquisition. The optical window for detection of the emitted light is set to the wavelength range of the Raman signals and thus should collect the same spectral range as ƩIBG. As shown in the previous section, this should allow to localise hot spots by a high intensity of IImg. As hot spots result from structures smaller than the wavelength of the incident light, they appear as a diffraction-limited spot in these images. Using an objective of high numerical aperture allows to localise these spots with high precision. The SERS measurement is subsequently performed by positioning the laser focus at the centre of highest intensity of these usually round-shaped spots and wavelength-dispersed detection of the Raman scattered light. This should distinctly increase the chance to obtain a SERS signal of high intensity compared to the random procedure. Thanks to the high speed of the imaging method also large areas of the SERS substrate can be measured in short time and the area can be further increased by stitching images. This significantly enhances the chance to localise several hot spots or – at best – the ‘hottest’ hot spots on the substrate.
SERS-BG based hot spot localisation and differentiation. To validate that a SERS-BG measurement images real hot spots applicable for enhancing Raman signals, we spectrally investigated the two hot spots with highest intensity of IImg pixel-wise after SERS-BG imaging (Fig. 3a). They both show a gradual increase of IImg towards the centre of the hot spot, even reaching the detector limit there (Fig. 3b and c, top). The intensities of the spectral signals (shown below the enlarged SERS-BG image) increase towards the centre of the hot spots in correspondence to IImg. For hot spot #1 we observe enhancement of ATP signals at 731 cm− 1, 1006 cm1 and 1334 cm− 1 (Fig. 3b) while for hot spot #2 strong signals are due to carbonisation apparent by the so-called D- and G-band of amorphous carbon (AC) at 1313 cm− 1 and 1592 cm− 1, respectively11–14 (Fig. 3c). Such carbonisation is a commonly observed phenomenon in SERS experiments and is most often attributed to contaminations12,13 or due to photophysical reactions of the analyte13. For quantitative evaluation we integrated the intensity data of the whole spectrum, the ATP signal at 731 cm− 1 and the G-band at 1592 cm− 1 to obtain ƩIBG, ƩIATP and ƩIAC, respectively (Fig. S3). The integrals, as the spectral signals, increase towards the centre of the hot spots (Fig. 3b and c bottom).
Comparison of the integration results also reveals the correlation between IImg, ƩIBG and the spectral signal of ATP or AC, respectively (Fig. S4). In hot spot #1 the profile of ATP enhancement matches with ƩIBG (R2 = 0.97) as well as with the pixel intensities IImg obtained from the SERS-BG image (R2 = 0.93). The integrated ATP signal increases from 0.1∙104 at the edges of the hot spot to 3.4∙104 in the centre, corresponding to an additional signal enhancement factor of 34.
In hot spot #2 we observe a similar correlation of the integrated AC band at 1592 cm− 1, ƩIBG and the SERS-BG pixel intensities IImg. However, a linear relationship is only found between ƩIAC and ƩIBG. This is mainly due to saturation of the HyD detector at the three central pixels, which prevents a correlation of IImg with ƩIAC or ƩIBG (Fig. S4). The saturation of the detector is most likely caused by surface-enhanced resonance Raman-scattering of amorphous carbon species leading to an enormous signal enhancement13. This also results in a four times higher ƩIBG at the AC-enhancing hot spot #2 when compared to the ATP-enhancing hot spot #1.
To differentiate between the desired ATP signal enhancing hot spots and the unwanted AC-enhancing hot spots we applied a multimodal SERS-BG imaging approach. In addition to SERS-BG imaging in the full range from 640–720 nm we also performed SERS-BG imaging restricted to the spectral region of the intense AC D- and G-band from 680–720 nm (SERS-BG-AC). This method assesses rough variations of the spectral shape. Figure 4 exemplarily compares the resulting SERS-BG images for two hot spots (#3 and #4). Hot spot #3 shows a higher intensity in the SERS-BG channel (green) compared to the SERS-BG-AC channel (red), indicating that the main spectral intensity is found in the range of 640–680 nm. Hot spot #4 shows similar intensities in both channels, hinting towards high signal intensities in the AC signal region. Spectral acquisition reveals that the increased intensities in the SERS-BG-AC channel of hot spot #4 are due to the AC bands. In contrast the spectrum in hot spot #3 shows mostly ATP signals.
Our observations show that the diffraction-limited spots of highest SERS-BG intensity IImg on the SERS substrate found by SERS-BG imaging are indeed real hot spots and allow collection of well-enhanced SERS spectra. The exact position of the hot spot is (as expected) at the centre of the diffraction-limited spot, where SERS-BG pixel intensity as well as the intensity of spectral bands in the SERS spectrum are highest in our experiments. This proves, that the pixel intensities IImg of SERS-BG imaging assess the sum of the signal intensities analogous to ƩIBG from SERS spectra, which correlate to the analyte signal intensities (see above). The SERS-BG pixel intensity profile thus follows the signal enhancement of the SERS-BG as well as of the analyte signals and is hence well suited for hot spot localisation. By using a multimodal SERS-BG imaging approach we are furthermore able to easily discriminate between hots spots enabling the measurement of useful SERS spectra and hot spots showing unwanted carbonisation effects.
Multimodal SERS-BG based SERS substrate characterisation. In order to perform more extensive optical characterisation of the SERS substrate we combined SERS-BG imaging with reflection and transmission imaging of the incident light. Additional micro- and nanostructure imaging by SEM yielded complementary information on the surface structure of the SERS substrate.
Figure 5: Multimodal characterisation of SERS substrate. (a) Reflection (Ex 633 nm, Em 630–636 nm, 1% laser intensity), transmission (Ex 633 nm, Em 560– 750 nm, 60% laser intensity) and SERS-background images (Ex 633 nm, Em 640–720 nm, 60% laser intensity) as well as a SEM recording from the same area. SEM imaging reveals a highly chaotic structure consisting of brighter protruding hill-like regions and darker valley-like regions, which end in hole-like structures at the centre. (b) Binarised images of (a). For the optical images features of high intensity are shown (I > 0.4∙Imax). For reflection, transmission and SERS-BG image 354, 114 and 44 features are determined with an average size of 634 ± 328 nm, 938 ± 691 nm and 525 ± 255 nm, respectively. For the SEM image the hole-like structures of low intensity are shown (I < 0.3∙Imax). 1753 features are present in the image with a mean size of 162 ± 105 nm. (c) Overlay of all four measurement images (top) and binarised images (bottom). (d) Correlation between the optical parameters plotted in a multimodal room spanned by reflection, transmission and SERS-BG measurements. Each circle represents one pixel on the SERS substrate. Data are normalised to Imax for each optical modality and color-coded according to SERS-BG intensity.
We further tried to correlate the optical properties with specific nanostructural elements determined by SEM (Fig. 5a, right). However due to the chaotic structure of the SERS substrate and its complex influence on the generation of SERS signals a simple separation of the images into features of a specific intensity range was not sufficient to link structural elements to optical behaviour (Fig. 5c), which is in line with Bell et al.4.
Assessment of sampling error in conventional SERS experiments by Monte-Carlo simulations. When performing SERS experiments the measured area and measurement position have a very high influence on the uniformity and repeatability of the signal enhancement and thus on the obtained SERS spectrum itself. To estimate this influence we performed Monte-Carlo simulations on the SERS-BG image of the commercial SERS substrate. We assess the impact of the measurement area size on the uniformity of the signal and the number of hot spots within the measurement area.
For our simulations, we randomly positioned a circular mask on the SERS-BG image mimicking the conventional random SERS analysis (see also Fig. 2a; for information on handling of the image edges see Fig. S5). For the area within the mask we calculated the mean SERS-BG intensity and determined the number of hot spots present. The deviation of the mean SERS-BG intensity within the mask from that of the whole image is taken as the SERS-BG sampling error (SERS-BG SE). The number of hot spots NHS is determined by counting the hot spots present within the mask area. From 1000 simulation repetitions, we calculated the average and standard deviation (SD) of both, the SERS-BG SE and NHS. Simulations were performed for different sizes of the circular mask thus mimicking different laser spot sizes.
Figure 6 shows the influence of the mask size on the average SERS-BG SE and on NHS. The average SERS-BG signal (red dots in Fig. 6b) is close to the average intensity of the whole image independent of the spot size and only shows minor uncorrelated fluctuations. The SDs for the 1000 repeated simulations however clearly increase with decreasing mask size reaching a maximum of 46% at a minimum mask diameter of 0.29 µm revealing a high variability of the SERS-BG signals across the investigated area. Thus large spot sizes are required for reproducible measurements. The average number of hot spots within the mask (Fig. 6d) increases with the square of the mask diameter, which is due to a linear relationship between the hot spot number and the area of the mask and approaches zero at small spot sizes. This emphasises the importance of the chosen laser focus size when performing a SERS measurement in the random way described above. An objective with high numerical aperture of 1.2 for example theoretically creates a laser spot size of approximately 0.64 µm in diameter at 633 nm resulting in a SD of around 39% and a number of detected hot spots close to zero and thus is completely unsuitable for reproducible SERS measurements. Applying an objective of lower numerical aperture of 0.3 leads to a theoretical laser spot size of 2.57 µm in diameter, which almost halves the SD to around 20%. The probability to catch at least one hot spot (Figure S6a) is just around 2.3% – only in 23 out of the 1000 simulations a hot spot is assessed. Therefore, even an objective of low numerical aperture is not very suitable. Reproducible and reliable measurements are only possible for very large measurement spot sizes. This is also mentioned by the manufacturer of the substrate, who recommends a laser spot size of at least 20 µm48. From our simulations this would result in a SD of around 9% for the SERS-BG SE and a mean hot spot number of 1.3 ± 1.1. This allows to obtain a SERS signal from at least one or two hot spots in 78% or 37% of the measurements, respectively (Figure S6a). However, as already mentioned by Crawford et al.28 increasing the laser spot size above the objective-defined size requires specific technical equipment to widen the laser spot and concomitantly requires an increase in the laser power, which is not possible in all commercially available Raman spectrometers. Also, when using objectives with lower numerical aperture the collection efficiency of the scattered light is decreased, thus annihilating in part the SERS enhancement of the signals. Therefore, this approach would not be the method of choice for many applications.
An alternative approach to increase the reliability of the measurement is to increase the number of measurements in an experiment, i.e. performing multiple measurements at different positions. This allows to increase the sampled area without changing the laser spot size as the sampled area sums up with each additional measurement, assuming the selected positions do not overlap. Using the hypothetical 0.3 NA objective for example requires 60 randomly positioned measurements to cover the same sample area as with a laser spot of 20 µm in diameter. However, as mentioned above the probability to obtain an average signal of at least two hot spots (e.g. for more reproducible quantitative measurements) then is still quite low. At a laser spot size of 40 µm in diameter one measurement assesses an average hot spot number of 5.1 ± 2.2 and a probability of almost 100% to hit at least two hot spots. Covering an equivalent area with this 0.3 NA objective would require 242 repeat measurements. These numbers impressively highlight that without technical widening of the laser spot the random procedure requires very time consuming experiments with a large number of repetitions to obtain reliable averaged SERS signals. By contrast, our SERS-BG based imaging method is far less time consuming while simultaneously allowing for targeted and thus reliable hot spot measurements with objectives of any NA.
The advantages of our method over the random procedure have even more impact when an exactly defined number of hot spot is of interest. From our Monte-Carlo simulations we estimate that in a random procedure the maximum probabilities to obtain exactly 1 to 5 hot spots would be found at 19 µm, 25 µm, 34 µm, 37 µm and 41 µm laser spot diameter with probabilities of 42%, 38%, 24%, 21% and 17%, respectively (Fig. 6b). Due to the high SERS substrate heterogeneity precise SERS analytics on a defined number of hot spots is hence not achievable by the standard random procedure. Also, for some experiments it may be useful to not have averaged signals, e.g. for single molecule SERS measurements, which are ideally based on just a single hot spot. Our SERS-BG based imaging method provides a suitable solution for both tasks. Upon imaging the substrate and localising the hot spots, suitable measurement positions can be chosen that provide exactly the desired number of hot spots.