A hybrid atomistic electrodynamics-quantum mechanical approach for simulating surface-enhanced Raman scattering.

Surface-enhanced Raman scattering (SERS) is a technique that has broad implications for biological and chemical sensing applications by providing the ability to simultaneously detect and identify a single molecule. The Raman scattering of molecules adsorbed on metal nanoparticles can be enhanced by many orders of magnitude. These enhancements stem from a twofold mechanism: an electromagnetic mechanism (EM), which is due to the enhanced local field near the metal surface, and a chemical mechanism (CM), which is due to the adsorbate specific interactions between the metal surface and the molecules. The local field near the metal surface can be significantly enhanced due to the plasmon excitation, and therefore chemists generally accept that the EM provides the majority of the enhancements. While classical electrodynamics simulations can accurately simulate the local electric field around metal nanoparticles, they offer few insights into the spectral changes that occur in SERS. First-principles simulations can directly predict the Raman spectrum but are limited to small metal clusters and therefore are often used for understanding the CM. Thus, there is a need for developing new methods that bridge the electrodynamics simulations of the metal nanoparticle and the first-principles simulations of the molecule to facilitate direct simulations of SERS spectra. In this Account, we discuss our recent work on developing a hybrid atomistic electrodynamics-quantum mechanical approach to simulate SERS. This hybrid method is called the discrete interaction model/quantum mechanics (DIM/QM) method and consists of an atomistic electrodynamics model of the metal nanoparticle and a time-dependent density functional theory (TDDFT) description of the molecule. In contrast to most previous work, the DIM/QM method enables us to retain a detailed atomistic structure of the nanoparticle and provides a natural bridge between the electronic structure methods and the macroscopic electrodynamics description. Using the DIM/QM method, we have examined in detail the importance of the local environment on molecular excitation energies, enhanced molecular absorption, and SERS. Our results show that the molecular properties are strongly dependent not only on the distance of the molecule from the metal nanoparticle but also on its orientation relative to the nanoparticle and the specific local environment. Using DIM/QM to simulate SERS, we show that there is a significant dependence on the adsorption site. Furthermore, we present a detailed comparison between enhancements obtained from DIM/QM simulations and those from classical electrodynamics simulations of the local field. While we find qualitative agreement, there are significant differences due to the neglect of specific molecule-metal interactions in the classical electrodynamics simulations. Our results highlight the importance of explicitly considering the specific local environment in simulations of molecule-plasmon coupling.

Direct test of the equivalency of dynamic IR and dynamic raman spectroscopies as techniques for observing ultrafast molecular dynamics.

We report the temperature-dependent infrared (IR) and Raman spectra of Fe(CO)3(η(4)-norbornadiene). This molecule undergoes carbonyl ligand site exchange on the vibrational time scale, and the effect of this exchange is observable as coalescence of the carbonyl bands in both the IR and Raman spectra. We outline a theory that we used to account for these effects and report simulations of the experimental spectra. We used these simulations to extract the carbonyl ligand exchange rates at various temperatures from the IR and Raman data. This data was used to calculate the activation energy for carbonyl exchange, yielding activation energies of 1.2 ± 0.2 and 1.4 ± 0.1 kcal/mol from the IR and Raman data, respectively. These activation energies are statistically identical and are consistent with previously reported values. This constitutes the first direct comparison between dynamic IR and Raman spectroscopies, and we find them to give identical results.

Tuning surface-enhanced Raman scattering from graphene substrates using the electric field effect and chemical doping.

Graphene recently has been demonstrated to support surface-enhanced Raman scattering. Here, we show that the enhancement of the Raman signal of methylene blue on graphene can be tuned by using either the electric field effect or chemical doping. Both doping experiments show that hole-doped graphene yields a larger enhancement than one which is electron-doped; however, chemical doping leads to a significantly larger modulation of the enhancements. The observed enhancement correlates with the changes in the Fermi level of graphene, indicating that the enhancement is chemical in nature, as electromagnetic enhancement is ruled out by hybrid electrodynamical and quantum mechanical simulations.

Photoreaction of matrix-isolated dihydroazulene-functionalized molecules on Au{111}.

Dihydroazulenes are photochromic molecules that reversibly switch between two distinct geometric and conductivity states. Molecular design, surface attachment, and precise control over the assembly of such molecular machines are critical in order to understand molecular function and motion at the nanoscale. Here, we use surface-enhanced Raman spectroscopy on special atomically flat, plasmonically enhanced substrates to measure the photoreaction kinetics of isolated dihydroazulene-functionalized molecules assembled on Au{111}, which undergo a ring-opening reaction upon illumination with UV light and switch back to the initial isomer via thermal relaxation. Photokinetic analyses reveal the high efficiency of the dihydroazulene photoreaction on solid substrates compared to other photoswitches. An order of magnitude decrease in the photoreaction cross section of surface-bound dihydroazulenes was observed when compared with the cross sections of these molecules in solution.

A discrete interaction model/quantum mechanical method for simulating surface-enhanced Raman spectroscopy.

We have derived and implemented analytical gradients for the discrete interaction model/quantum mechanics (DIM/QM) method. DIM/QM combines an atomistic electrodynamics model with time-dependent density functional theory and thus enables modeling of the optical properties for a molecule while taking into account the local environment of a nanoparticle's surface. The DIM/QM analytical gradients allow for geometry optimizations, vibrational frequencies, and Raman spectra to be simulated for molecules interacting with metal nanoparticles. We have simulated the surface-enhanced Raman scattering (SERS) spectra for pyridine adsorbed on different sites of icosahedral nanoparticles with diameters between 1 and 8 nm. To describe the adsorption of the pyridine molecule onto the metal surface, we have implemented a coordination-dependent force field to differentiate the various local surface environments. We find that the DIM/QM method predicts geometries and frequencies that are in good agreement with full QM simulations and experiments. For the simulated SERS spectra of pyridine, we find a significant dependence on the adsorption site and the size of the metal nanoparticle. This illustrates the importance of accounting for the local environment around the molecule. The Raman enhancement factors are shown to roughly mirror the magnitude of the nanoparticle's local field about the molecule. Because the simulated nanoparticles are small, the plasmon peaks are quite broad which results in weak local electric fields and thus modest Raman enhancement factors.

Importance of Correctly Describing Charge-Transfer Excitations for Understanding the Chemical Effect in SERS.

The enhancement mechanism due to the molecule-surface chemical coupling in surface-enhanced Raman scattering (SERS) is governed to a large extent by the energy difference between the highest occupied molecular orbital (HOMO) of the metal and the lowest unoccupied molecular orbital (LUMO) of the molecule. Here, we investigate the importance of correctly describing charge-transfer excitations, using time-dependent density functional theory (TDDFT), when calculating the chemical coupling in SERS. It is well-known that TDDFT, using traditional functionals, underestimates the position of charge-transfer excitations. Here, we show that this leads to a significant overestimation of the chemical coupling mechanism in SERS. Significantly smaller enhancements are found using long-range corrected (LC) functionals as compared with a traditional generalized gradient approximation (GGA) and hybrid functionals. Enhancement factors are found to be smaller than 530 and typically less than 50. Our results show that it is essential to correctly describe charge-transfer excitations for predicting the chemical enhancement in SERS.

A discrete interaction model/quantum mechanical method to describe the interaction of metal nanoparticles and molecular absorption.

A frequency-dependent quantum mechanics/molecular mechanics method for the calculation of response properties of molecules adsorbed on metal nanoparticles is presented. This discrete interaction model/quantum mechanics (DIM/QM) method represents the nanoparticle atomistically, thus accounting for the local environment of the nanoparticle surface on the optical properties of the adsorbed molecule. Using the DIM/QM method, we investigate the coupling between the absorption of a silver nanoparticle and of a substituted naphthoquinone. This system is chosen since it shows strong coupling due to a molecular absorption peak that overlaps with the plasmon excitation in the metal nanoparticle. We show that there is a strong dependence not only on the distance of the molecule from the metal nanoparticle but also on its orientation relative to the nanoparticle. We find that when the transition dipole moment of an excitation is oriented towards the nanoparticle there is a significant increase in the molecular absorption as a result of coupling to the metal nanoparticle. In contrast, we find that the molecular absorption is decreased when the transition dipole moment is oriented parallel to the metal nanoparticle. The coupling between the molecule and the metal nanoparticle is found to be surprisingly long range and important on a length scale comparable to the size of the metal nanoparticle. A simple analytical model that describes the molecule and the metal nanoparticle as two interacting point objects is found to be in excellent agreement with the full DIM/QM calculations over the entire range studied. The results presented here are important for understanding plasmon-exciton hybridization, plasmon enhanced photochemistry, and single-molecule surface-enhanced Raman scattering.

A discrete interaction model/quantum mechanical method for describing response properties of molecules adsorbed on metal nanoparticles.

A new polarizable quantum mechanics/molecular mechanics method for the calculation of response properties of molecules adsorbed on metal nanoparticles is presented. This method, which we denote the discrete interaction model/quantum mechanics (DIM/QM) method, represents the nanoparticle atomistically which enables the modeling of the influence of the local environment of a nanoparticle surface on the optical properties of a molecule. Using DIM/QM, we investigate the excitation energies of rhodamine-6G (R6G) and crystal violet (CV) adsorbed on silver and gold nanoparticles of different quasispherical shapes and sizes. The metal nanoparticle is characterized by its static total polarizability, a reasonable approximation for frequencies far from the plasmon resonance. We observe that for both R6G and CV, the presence of the nanoparticle shifts the strongest excitation to the red approximately 40 nm and also increases the oscillator strength of that excitation. The shifts in excitation energies due to the nanoparticle surface are found to be comparable to those due to solvation. We find that these shifts decay quickly as the molecule is moved away from the surface. We also find that the wavelength shift is largest when the transition dipole moment is aligned with the edges of the nanoparticle surface where the electric field is expected to be the largest. These results show that the molecular excitations are sensitive to the local environment on the nanoparticle as well as the specific orientation of the molecule relative to the surface.

Controlling the non-resonant chemical mechanism of SERS using a molecular photoswitch.

In this work we present a detailed investigation of the Raman properties of a dithienylethene photoswitch interacting with a small gold cluster (Au(19)(+)) using time-dependent density functional theory (TD-DFT). The enhancement mechanism (CHEM) due to the molecule-surface chemical coupling in surface-enhanced Raman scattering (SERS) has been characterized for this system. We demonstrate that it is possible to control the CHEM enhancement by switching the photoswitch from its closed form to its open form. The open form of the photoswitch is found to be the strongest Raman scatterer when adsorbed on the surface whereas the opposite is found for the free molecule. This trend is explained using a simple two-state approximation. In this model the CHEM enhancement scales roughly as (omega(X)/omega(e)(4)), where omega(X) is the HOMO-LUMO gap of the free molecule and omega(e) is an average excitation between the HOMO of the photoswitch and the LUMO of the metal. We propose that the ability of this photoswitch to switch reversibly from open to closed will make it an excellent probe to control the CHEM enhancement of SERS.

Understanding the molecule-surface chemical coupling in SERS.

The enhancement mechanism due to the molecule-surface chemical coupling in surface-enhanced Raman scattering (SERS) has been characterized using time-dependent density functional theory. This has been achieved with a systematical study of the chemical enhancement of meta- and para-substituted pyridines interacting with a small silver cluster (Ag(20)). Changing the functional groups on pyridine enabled us to modulate the direct chemical interactions between the pyridine ring and the metal cluster. Surprisingly, we find that the enhancement does not increase as more charge is transferred from the pyridine ring to the cluster. Instead, we find that the magnitude of chemical enhancement is governed to a large extent by the energy difference between the highest occupied energy level (HOMO) of the metal and the lowest unoccupied energy level (LUMO) of the molecule. The enhancement scales roughly as (omega(X)/omega(e))(4), where omega(e) is an average excitation energy between the HOMO of the metal and the LUMO of the molecule and omega(X) is the HOMO-LUMO gap of the free molecule. The trend was verified by considering substituted benzenethiols, small molecules, and silver clusters of varying sizes. The results imply that molecules that show significant stabilization of the HOMO-LUMO gaps (such as those that readily accept pi-backbonding) would be likely to have strong chemical enhancement. The findings presented here provide the framework for designing new molecules which exhibit high chemical enhancements. However, it remains a challenge to accurately describe the magnitude of the Raman enhancements using electronic structure methods, especially density functional theory, because they often underestimate the energy gap.