Abstract

We have employed capillary force deposition on ion beam patterned substrates to fabricate chains of 60 nm gold nanospheres ranging in length from 1 to 9 nanoparticles. Measurements of the surface-averaged SERS enhancement factor strength for these chains were then compared to the numerical predictions. The SERS enhancement conformed to theoretical predictions in the case of only a few chains, with the vast majority of chains tested not matching such behavior. Although all of the nanoparticle chains appear identical under electron microscope observation, the extreme sensitivity of the SERS enhancement to nanoscale morphology renders current nanofabrication methods insufficient for consistent production of coupled nanoparticle chains. Notwithstanding this fact, the aggregate data also confirmed that nanoparticle dimers offer a large improvement over the monomer enhancement while conclusively showing that, within the limitations imposed by current state-of-the-art nanofabrication techniques, chains comprising more than two nanoparticles provide only a marginal signal boost over the already considerable dimer enhancement.

1. Introduction

Scientific interest in surface-enhanced Raman spectroscopy (SERS) has increased dramatically since the detection of single molecules was first reported in 1997 [1]. While enhancement factors measured from bulk SERS-active materials typically hover around 106 [2, 3], researchers have succeeded in designing “hot spots” on substrates capable of producing values over 108 (on the order of the level necessary for detecting single molecules) [47]. The prospect of being able to boost sensitivity of SERS measurements at the single molecule level coupled with the fingerprint-like uniqueness of the Raman spectrum makes SERS a sought-after candidate for incorporation into sensor technologies; particularly, for those who are interested in detecting chemical compounds at the trace level.

Although researchers have discovered many nanoscale-textured metallic surfaces capable of producing the strong electromagnetic fields that give rise to SERS enhancements, metal nanoparticle aggregates have long been recognized as ideal materials for the achievement of truly large enhancement factors. Among these, the simplicity of the spherical nanoparticle dimer has made it a structure of choice in many attempts to build arrays of highly sensitive and reproducible SERS-active structures [8]. In a recent study, however, the question has been posed as to whether the two-particle linear aggregate is, indeed, the optimal configuration for generating SERS enhancements, prompting researchers to consider other morphologies [9]. One such configuration is closely spaced metal nanoparticles arranged end to end to form a chain (a natural extension of the dimer configuration). This structure was theoretically investigated by Wang et al. [9], who carried out finite integral technique (FIT) simulations on gold and silver nanoparticle chains comprised of even numbers of nanoparticles ranging in length from 2 nanoparticles to infinity. The results of the simulations indicate that the field enhancement at the gap between the two centermost particles (Wang refers to this as the “c-point”) can reach a maximum at a number larger than two, suggesting that the dimer might not always be the ideal configuration for achieving the maximum field enhancement. In contrast to these predictions, recent theoretical and experimental results by Wustholz et al. [10] suggest that, in most realistic experiments, the signal from the hottest hot spot in a cluster is often extremely dominant and, thus, renders unrealistic the possibility of enhancement due to near-field coupling along the chain as a whole. In this paper, the main aim is to test these conclusions on precisely designed nanoparticle chains. Through the combination of focused ion beam (FIB) milling and capillary forced deposition [11], closely spaced linear nanoclusters were fabricated. The SERS enhancement arising from these structures were measured and subsequently compared with theoretically predicted values.

2. Generalized Mie Theory Simulations

Because it would not be possible for us to measure the SERS enhancement factor (EF) directly at the c-point, we have carried out generalized Mie theory (GMT) [1215] simulations to map the electromagnetic strength over the entire surface of gold nanospheres chains ranging in length from 1 to 9 nanoparticles with interparticle separations of 1–1.5 nm (these separations still yield reliable results in the frame of classical electromagnetic theory).

The surface Raman enhancement is proportional to the product of the squares of the field enhancement at the incident and Raman frequencies,||𝑓𝜔𝐺=Raman||2||𝑓𝜔pump||2.(1) Provided that |𝜔Raman𝜔| is smaller than the spectral response of the metal nanostructure and the excitation source is polarized along the axis of the chain, the SERS enhancement can be approximated to scale with the fourth power of the electric field enhancement [16]:||𝐸𝐺(𝑟)loc||(𝑟)4||𝐸0||(𝑟)4.(2) Since the molecular species to be detected is assumed to form a monolayer on the surface of the nanospheres, the Raman EF must be averaged over the entire chain surface1𝐺(𝑟)=𝑁𝑖4𝜋𝑅2𝑖𝑁𝑖𝑖𝐺(𝑟)𝑑𝜎,(3) where 𝑁 is the number of nanoparticles in the chain, 𝑅𝑖 is the radius of the 𝑖th nanoparticle, and the fields were calculated at a distance of 0.5 nm from the surface which, based on the estimated thickness of the benzenethiol monolayer (the reporter molecule used in the experimental part of this study), is the most likely location of the reporter molecule. Here, it is important to note that, although the length of the benzenethiol molecule and bond length on gold are known, the interparticle distance for the structures fabricated in this experiment cannot be known exactly. Similarly, a monolayer of thiophenol molecules surrounding the spheres and the silica/silicon substrate invariably alters the index of refraction of the area immediately surrounding the metal surface. To account for this range in gap size and the effective index of refraction of the surrounding medium, 𝑛e, we calculated the surface averaged enhancement for a gap of 1 nm and 𝑛e=1.0 (Figure 1(a)) as well as for a gap of 1.5 nm and 𝑛e=1.5 (Figure 1(b)). As expected, the surface Raman enhancements of the nanochains (𝑁2) reach values on the order of 107-108 over that of the monomer at resonance. Furthermore, the plasmon resonances shift to longer wavelengths as the particle number increases. Figures 1(c) and 1(d) show the SERS enhancement at the laser line (632.8 nm) as a function of particle number, corresponding to Figures 1(a) and 1(b), respectively. Significantly, a peak is apparent in the SERS enhancement at a chain length of 3 nanoparticles for the 1 nm gap (𝑛e=1.0) and 2 nanoparticles for the 1.5 nm gap (𝑛e=1.5). It is important to clarify that this simulation approach was designed to produce results that could be compared to experimental measurements of the SERS EF where reporter molecules actually coat all metallic surfaces. This is in contrast to the FIT calculations performed by Wang et al. which by focusing in the field enhancement at the “c-point” cannot be directly compared to a far-field experimental measurement.

3. Finite Difference Time Domain (FDTD) Simulations

To identify the link between nanoparticle chain length and the location of the surface Raman enhancement peak, we also calculated the extinction spectrum for the linear chains. As our GMT code is designed for near-field calculations, we have rather turned to apply a fast and commercially available finite difference time domain (FDTD) method (Lumerical Solutions, Inc) for far-field results. While FDTD cannot compete with the accuracy of GMT for the near field calculations, it is a perfect alternative in the far field. In order to ensure accuracy, simulations were carried out with a mesh override region surrounding the entire structure with cell size limited to 0.8 nm and convergence of the results were confirmed by comparing the results of multiple runs.

Results from these simulations are shown in Figure 2. First, the redshift of the extinction peak with respect to the increased particle number is apparent, indicating as expected, the tight correlation between the optical response and the morphology of the scatterer. Most importantly, however, the location of the highest peak of the extinction with respect to wavelength coincides with the highest EF in the near-field calculations in Figure 1. Here, we note that when we narrow our observations to a fixed wavelength (633 nm in this case) small secondary peaks become relevant, making the direct comparison between the far and near fields less direct. Figures 1(c) and 1(d) plot the EF magnitude at the fixed wavelength againts the number of particles in each chain. One can observe that the EF peaks at a particular number of particles, and that number is not necessarily two (Figure 1(c)).

4. Sample Nanofabrication

To experimentally test the simulation results, a templating method was employed to create nanoparticle chains of specific lengths and orientations. FIB milling was used to create trenches in SiO2 which were approximately 80 nm in width, 70 nm in depth, and ranging from 80 to 560 nm in length. Gold nanospheres (Ted Pella) with a diameter of 60 nm were placed into these trenches using a capillary force deposition method [11]. To accomplish this, the nanopatterned substrate was cleaned in piranha solution and subsequently subjected to a brief plasma treatment to ensure the hydrophilicity of the SiO2 layer. Next, the SiO2 was immersed in an aqueous gold colloid which was allowed to evaporate at 70° overnight. As the meniscus slowly recedes across the substrate, nanoparticles concentrate at the three-phase contact line (Figure 3). Because the trenches are designed to host only a specific number and configuration of nanoparticles, this process creates arrays of nanoparticle chains of predetermined lengths. Scanning-electron micrographs of the substrates were taken to identify satisfactory chains. To remove carbonaceous contamination that often accompanies SEM imaging, substrates were subsequently cleaned with piranha solution and treated with O2 plasma for 10 minutes to remove any unwanted material that the electron beam might have deposited. Finally, nanochains were immersed in a 10 mM solution of benzenethiol overnight, removed, rinsed for 60 seconds with ethanol to remove excess benzenethiol, and then blown dry with N2.

5. Raman Measurements

Raman spectra of individual chains were obtained using a Leica microscope equipped with a confocal Raman spectroscopic system (Renishaw InVia) and a 30 mW 632.8 nm laser excitation source. All spectra were measured using a 50× objective with a numerical aperture of 0.75. The Raman signal was collected by a TE air-cooled 576×400 CCD array preceded by two notch filters (OD>12) to block the laser line. Preliminary tests concluded that, without sufficient filtering, persistent exposure to the laser beam causes the benzenethiol monolayer to degrade. To avoid photodegradation of the monolayer, the stability of the Raman signal was monitored over a range of beam powers and exposure times. These tests determined that the combination of a 100 μW laser beam power (spot size diameter ~2 μm) delivered on the sample with a 30-second acquisition time yielded sufficiently clear and stable spectra. To account for the possibility that the laser spot was not perfectly aligned over the chain and/or not perfectly uniform, a flat scanning stage was used to raster over the structure. From these measurements, the strongest, most well-defined spectrum from each of the structures was selected for data analysis. A typical spectrum for a nanoparticle heptamer is shown in Figure 4. The most pronounced peaks appear at 1000 cm−1 (ring out-of-plane deformation and C–H out-of-plane bending), 1080 cm−1 (C–C symmetric stretching and C–S stretching), and 1577 cm−1 (C–C symmetric stretching) [17]. Time-stability tests were performed in order to select the most stable Raman peak for comparison in these tests. All the peaks were clear and steady, but the 1080 cm−1 peak was found to show the best signal to noise ratio and was, thus, used for all data analyses performed in this study.

6. Summary of Results and Discussion

Experimental results are plotted over GMT simulations for a 1.5 nm interparticle gap in Figure 5. Results from two data sets are shown to delineate the range of outcomes for this particular experiment. Specifically, the first set represents the largest measured enhancements (i.e., the best enhancements obtained in this study) for each chain length (red trace), and a second set with the results falling within one standard deviation of the overall data average (blue trace). It can be observed that the second set deviates from theory especially at longer chain lengths, reaching a plateau of enhancement rather than decreasing for larger particle numbers. The large enhancement ratio between the monomer and the dimer is unsurprising, as it is well known that the tiny gaps formed by closely spaced nanoparticles are the source of large electric fields which give rise to large SERS enhancement factors. This situation becomes more complex, however, when the number of gaps in the chain is increased through the addition of nanoparticles. It has been shown that the SERS EF is extremely sensitive to interparticle separation, falling off at the rate of approximately one order of magnitude per nanometer of separation [18]. For longer nanochains, a mixture of coupled and decoupled neighboring particles might coexist. Moreover, the extreme sensitivity of the SERS enhancement to a gamut of factors, including nanoparticle size [19], shape [20], crystal face [21], surface roughness [22], and particle-particle spacing [2327] is well documented and is known to account for significant variations in EFs arising from seemingly-identical nanoparticle clusters. Taking these two items into consideration, one can see that in practice it is unlikely that the enhancement that would arise from a perfectly coupled nanochain of a specific length is being observed. Moreover, it is quite likely that a single dominant hot spot is accounting for the majority of the enhancement observed for all chain lengths >1. Significantly, these considerations qualitatively explain the plateau of the EF in Figure 5 as well as emphasize the importance of making hot spots through gap formation over additional enhancements achieved through resonances in complex geometries, thereby supporting the Wustholz et al. [10] conclusion.

At the other end of the spectrum of experimental outcomes, the largest enhancement results do show a significant agreement with theory. While it would be impossible to accurately determine whether the peak at a chain length of 3 nanoparticles arises from perfect interparticle coupling along the entire length of the chain, the measurement seems to indicate that this was the case and it is certainly relevant to include a discussion of this ideal situation. Referring back to Figure 1, despite a small redshift, the maximum SERS enhancement for the monomers depicted in Figures 1(a) and 1(b) is located around 525 nm and 570 nm, respectively, corresponding to the far-field LSPR peak. For aggregates consisting of more than two spheres, the hybridized plasmon modes appear as a result of the interaction of primary plasmons in nearby spheres [28, 29]. This coupling dramatically enhances the local electric field, creating an SERS enhancement, that is, orders of magnitude larger than that of a single sphere. In addition to lowering the bonding mode, hybridized plasmon modes of higher angular momentum are also excited by the wave polarized parallel to the long axis of the chain. As the chain length is increased, both EF peaks experience a redshift. These chains demonstrate a sufficiently strong dipole moment and are, thus, able to couple to far field radiation efficiently, thereby creating the peaks in the EF observed in both the simulation and experimental data.

7. Conclusions

In summary, GMT and FDTD simulations were performed to model the surface-averaged electromagnetic field enhancements arising from plasmonic coupling in closely spaced noble metal nanoparticle chains. The simulations suggest that there exists an “ideal” chain length capable of producing modest improvements on the typical dimer enhancement factor. To experimentally test the model, chains of 60 nm gold nanospheres were fabricated using a capillary force deposition technique. Subsequent Raman measurements made from these structures showed modest agreement with the GMT simulations in some cases and a complete lack of agreement in others. In general, the vast majority of the tested nanoparticle chains showed no more enhancement beyond the typical for a two-particle dimmer. From these results, it was concluded that it is unlikely that the individual nanoparticles comprising the chains made using the capillary force deposition method technique are consistently able to couple completely along the entire structure, thereby departing from the ideal situation modeled in the calculations performed in this study. However, comparison of experimental and theoretical results confirms that the formation of dimers provides a large improvement in signal strength over the monomer enhancement. While perfect plasmonic coupling of the ideal number of nanoparticles can provide improvement over the dimer enhancement, it is marginal under the best of circumstances and often introduces more difficulties in the fabrication process than is justifiable by the small boost in signal.

Acknowledgments

This work was supported in part by the Hartley Corporation Fellowship for Graduate Women in Science and the SPIE Graduate Student fellowship. The authors thank Dr. Richard Superfine for many useful conversations.