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

Development of discrete interaction models for ultra-fine nanoparticle plasmonics

Plasmonics serves as a most outstanding feature of nanoparticle technology and is nowadays used in numerous applications within imaging, sensing and energy harvesting, like plasmonically enhanced solar cells, nanoparticle bioimaging, plasmon-controlled fluorescence for molecular tracking in living cells, plasmon-controlled electronic molecular devices and surface enhanced Raman spectroscopy for single molecular detection. Although plasmonics has been utilized since ancient times, the understanding of its basic interactions has not been fully achieved even under the emergence of modern nanoscience. In particular, it has been difficult to address the "ultra-fine" 1-10 nm regime, important for applications especially in bioimaging and biomedical areas, where neither classical nor quantum based theoretical methods apply. Recently, new approaches have been put forward to bridge this size gap based on semi-empirical discrete interaction models where each atom makes a difference. A primary aim of this perspective article is to review some of the most salient features of these models, and in particular focus on a recent extension - the extended discrete interaction model (Ex-DIM), where the geometric and environmental features are extended - and highlight a set of benchmark studies using this model concerning size, shape, material, temperature dependence and other characteristics of ultra-fine plasmonic nanoparticles. We also analyze new possibilities offered by the model for designing ultra-fine plasmonic particles for applications in the areas of bioimaging, biosensing, photothermal therapy, infrared light harvesting and photodetection. We foresee that future modelling activities will be closely connected to collaborative experimental work including synthesis, device fabrication and measurements with feedback and validation in a systematic fashion. With this strategy we can expect that modelling of ultra-fine plasmonics particles can be integrated in the development of novel plasmonic systems with unprecedented performance and applicability.

Time-dependent Kohn-Sham electron dynamics coupled with nonequilibrium plasmonic response via atomistic electromagnetic model

Computational modeling of plasmon-mediated molecular photophysical and photochemical behaviors can help us better understand and tune the bound molecular properties and reactivity and make better decisions to design and control nanostructures. However, computational investigations of coupled plasmon-molecule systems are challenging due to the lack of accurate and efficient protocols to simulate these systems. Here, we present a hybrid scheme by combining the real-time time-dependent density functional theory (RT-TDDFT) approach with the time-domain frequency dependent fluctuating charge (TD-ωFQ) model. At first, we transform ωFQ in the frequency-domain, an atomistic electromagnetic model for the plasmonic response of plasmonic metal nanoparticles (PMNPs), into the time-domain and derive its equation-of-motion formulation. The TD-ωFQ introduces the nonequilibrium plasmonic response of PMNPs and atomistic interactions to the electronic excitation of the quantum mechanical (QM) region. Then, we combine TD-ωFQ with RT-TDDFT. The derived RT-TDDFT/TD-ωFQ scheme allows us to effectively simulate the plasmon-mediated "real-time" electronic dynamics and even the coupled electron-nuclear dynamics by combining them with the nuclear dynamics approaches. As a first application of the RT-TDDFT/TD-ωFQ method, we study the nonradiative decay rate and plasmon-enhanced absorption spectra of two small molecules in the proximity of sodium MNPs. Thanks to the atomistic nature of the ωFQ model, the edge effect of MNP on absorption enhancement has also been investigated and unveiled.

A recursive cell multipole method for atomistic electrodynamics models

For large plasmonic nanoparticles, retardation effects become important once their length becomes comparable to the wavelength of light. However, most models do not incorporate retardation effects due to the high computational cost of solving for the optical properties of large atomistic electrodynamics systems. In this work, we derive and implement a recursive fast multipole method (FMM) in Cartesian coordinates that includes retardation effects. In this method, higher-order electrodynamic interaction tensors used for the FMM are calculated recursively, thus greatly reducing the implementation complexity of the model. This method allows for solving of the optical properties of large atomistic nanoparticles with controlled accuracy; in practice, taking the expansion to the fifth order provides a good balance of accuracy and computational time. Finally, we study the effects retardation has on the near- and far-field properties of large plasmonic nanoparticles with over a million atoms using this method. We specifically focus on nanorods and their dimers, which are known to generate highly confined fields in their junctions. In the future, this method can be applied to simulations in which accurate near-field properties are required, such as surface-enhanced Raman scattering.

Molecular Simulation Study of Surface-Enhanced Raman Scattering of Liquid Water

We developed in our previous study [J. Chem. Phys., 2021, 155, 174118, J. Phys. Chem. A, 2022, 126, 4762] a classical electronic and molecular dynamics simulation method to describe the optical response of metal material in solution based on an atomistic model by incorporating the classical equation of motion for free electrons under an applied electric field. To show further usefulness of the method, in the present study, we apply it to surface-enhanced Raman scattering of liquid water to examine the signal enhancement of the solution system caused by plasmon resonance effects of a silver nanoparticle.

Modeling the near-field effect on molecular excited states using the discrete interaction model/quantum mechanical method

Strong light-matter interactions significantly modify the optical properties of molecules in the vicinity of plasmonic metal nanoparticles. Since the dimension of the plasmonic cavity approaches that of the molecules, it is critical to explicitly describe the nanoparticle junctions. In this work, we use the discrete interaction model/quantum mechanical (DIM/QM) method to model the coupling between the plasmonic near-field and molecular excited states. DIM/QM is a combined electrodynamics/quantum mechanical model that uses an atomistic description of the nanoparticle. We extend the DIM/QM method to include the local field effects in the sum-over-state formalism of time-dependent density functional theory. As a test of the method, we study the interactions between small organic chromophores and metal nanoparticles. In particular, we examine how the inclusion of multiple electronic transitions and intermolecular interactions modify the coupling between molecules and nanoparticles. Using the sum-over-state formalism of DIM/QM, we show that two-state models break down when the plasmon excitation is detuned from the molecular excitations. To gain further insight, we compare the simple coupled-dipole model (CDM) with the DIM/QM model. We find that CDM works well for simple systems but fails when going beyond the single molecule or single nanoparticle cases. We also find that the coupling depends strongly on the site of the nanoparticle in which the chromophore couples to. Our work suggests the importance of explicitly describing the cavity to capture the atomistic level local field environment in which the molecule strongly couples to.

Accelerated constant-voltage quantum mechanical/molecular mechanical method for molecular systems at electrochemical interfaces

The structure and electronic properties of a molecule at an electrochemical interface are changed by interactions with the electrode surface and the electrolyte solution, which can be significantly modulated by an applied voltage. We present an efficient self-consistent quantum mechanics/molecular mechanics (QM/MM) approach to study a physisorbed molecule at a metal electrode-electrolyte interface under the constant-voltage condition. The approach employs a classical polarizable double electrode model, which enables us to study the QM/MM system in the constant-voltage ensemble. A mean-field embedding approximation is further introduced in order to overcome the difficulties associated with statistical sampling of the electrolyte configurations. The results of applying the method to a test system indicate that the adsorbed molecule is no less or slightly more polarized at the interface than in the bulk electrolyte solution. The geometry of the horizontally adsorbed molecule is modulated by their electrostatic interactions with the polarizable electrode surfaces and also the interactions with cations attracted toward the interface when the adsorbate is reduced. We also demonstrate that the approach can be used to quantitatively evaluate the reorganization energy of a one electron reduction reaction of a molecule in an electrochemical cell.

Analytic high-order energy derivatives for metal nanoparticle-mediated infrared and Raman scattering spectra within the framework of quantum mechanics/molecular mechanics model with induced charges and dipoles

This work is devoted to deriving and implementing analytic second- and third-order energy derivatives with respect to the nuclear coordinates and external electric field within the framework of the hybrid quantum mechanics/molecular mechanics method with induced charges and dipoles (QM/DIM). Using these analytic energy derivatives, one can efficiently compute the harmonic vibrational frequencies, infrared (IR) and Raman scattering (RS) spectra of the molecule in the proximity of noble metal clusters/nanoparticles. The validity and accuracy of these analytic implementations are demonstrated by the comparison of results obtained by the finite-difference method and the analytic approaches and by the full QM and QM/DIM calculations. The complexes formed by pyridine and two sizes of gold clusters (Au18 and Au32) at varying intersystem distances of 3, 4, and 5 Å are used as the test systems, and Raman spectra of 4,4'-bipyridine in the proximity of Au2057 and Ag2057 metal nanoparticles (MNP) are calculated by the QM/DIM method and compared with experimental results as well. We find that the QM/DIM model can well reproduce the IR spectra obtained from full QM calculations for all the configurations, while although it properly enhances some of the vibrational modes, it artificially overestimates RS spectral intensities of several modes for the systems with very short intersystem distance. We show that this could be improved, however, by incorporating the hyperpolarizability of the gold metal cluster in the evaluation of RS intensities. Additionally, we address the potential impact of charge migration between the adsorbate and MNPs.

Do We Really Need Quantum Mechanics to Describe Plasmonic Properties of Metal Nanostructures?

Optical properties of metal nanostructures are the basis of several scientific and technological applications. When the nanostructure characteristic size is of the order of few nm or less, it is generally accepted that only a description that explicitly describes electrons by quantum mechanics can reproduce faithfully its optical response. For example, the plasmon resonance shift upon shrinking the nanostructure size (red-shift for simple metals, blue-shift for d-metals such as gold and silver) is universally accepted to originate from the quantum nature of the system. Here we show instead that an atomistic approach based on classical physics, ωFQFμ (frequency dependent fluctuating charges and fluctuating dipoles), is able to reproduce all the typical "quantum" size effects, such as the sign and the magnitude of the plasmon shift, the progressive loss of the plasmon resonance for gold, the atomistically detailed features in the induced electron density, and the non local effects in the nanoparticle response. To support our findings, we compare the ωFQFμ results for Ag and Au with literature time-dependent DFT simulations, showing the capability of fully classical physics to reproduce these TDDFT results. Only electron tunneling between nanostructures emerges as a genuine quantum mechanical effect, that we had to include in the model by an ad hoc term.

Evaluation of molecular photophysical and photochemical properties using linear response time-dependent density functional theory with classical embedding: Successes and challenges

Time-dependent density functional theory (TDDFT) based approaches have been developed in recent years to model the excited-state properties and transition processes of the molecules in the gas-phase and in a condensed medium, such as in a solution and protein microenvironment or near semiconductor and metal surfaces. In the latter case, usually, classical embedding models have been adopted to account for the molecular environmental effects, leading to the multi-scale approaches of TDDFT/polarizable continuum model (PCM) and TDDFT/molecular mechanics (MM), where a molecular system of interest is designated as the quantum mechanical region and treated with TDDFT, while the environment is usually described using either a PCM or (non-polarizable or polarizable) MM force fields. In this Perspective, we briefly review these TDDFT-related multi-scale models with a specific emphasis on the implementation of analytical energy derivatives, such as the energy gradient and Hessian, the nonadiabatic coupling, the spin-orbit coupling, and the transition dipole moment as well as their nuclear derivatives for various radiative and radiativeless transition processes among electronic states. Three variations of the TDDFT method, the Tamm-Dancoff approximation to TDDFT, spin-flip DFT, and spin-adiabatic TDDFT, are discussed. Moreover, using a model system (pyridine-Ag20 complex), we emphasize that caution is needed to properly account for system-environment interactions within the TDDFT/MM models. Specifically, one should appropriately damp the electrostatic embedding potential from MM atoms and carefully tune the van der Waals interaction potential between the system and the environment. We also highlight the lack of proper treatment of charge transfer between the quantum mechanics and MM regions as well as the need for accelerated TDDFT modelings and interpretability, which calls for new method developments.

State-averaged CASSCF with polarizable continuum model for studying photoreactions in solvents: Energies, analytical nuclear gradients, and non-adiabatic couplings

This paper presents state-averaged complete active space self-consistent field in polarizable continuum model (PCM) for studies of photoreactions in solvents. The wavefunctions of the solute and the PCM surface charges of the solvent are optimized simultaneously such that the state-averaged free energy is variationally minimized. The method supports both fixed weights and dynamic weights where the weights are automatically adjusted based on the energy gaps. The corresponding analytical nuclear gradients and non-adiabatic couplings are also derived. Furthermore, we show how the new method can be entirely formulated in terms of seven basic operations, which allows the implementation to benefit from existing high-performance libraries on graphical processing units. Results demonstrating the accuracy and performance of the implementation are presented and discussed. We also apply the new method to the study of minimal conical intersection search and photoreaction energy pathways in solvents. Effects from the polarity of the solvents and different formulas of dynamic weights are compared and discussed.

Local-Field Effects in Linear Response Properties within a Polarizable Frozen Density Embedding Method

In this work, we present a polarizable frozen density embedding (FDE) method for calculating polarizabilities of coupled subsystems. The method (FDE-pol) combines a FDE method with an explicit polarization model such that the expensive freeze/thaw cycles can be bypassed, and approximate nonadditive kinetic potentials are avoided by enforcing external orthogonality between the subsystems. To describe the polarization of the frozen environment, we introduce a Hirshfeld partition-based density-dependent method for calculating the atomic polarizabilities of atoms in molecules, which alleviates the need to fit the atomic parameters to a specific system of interest or to a larger general set of molecules. We show that the Hirshfeld partition-based method predicts molecular polarizabilities close to the basis set limit, and thus, a single basis set-dependent scaling parameter can be introduced to improve the agreement against the reference polarizability data. To test the model, we characterized the uncoupled and coupled response of small interacting molecular complexes. Here, the coupled response properties include the perturbation of the frozen system due to the external perturbation which is ignored in the uncoupled response. We show that FDE-pol can accurately reproduce both the exact uncoupled polarizability and the coupled polarizabilities of the supermolecular systems. Using damped response theory, we also demonstrate that the coupled frequency-dependent polarizability can be described by including local field effects. The results emphasize the necessity of including local-field effects for describing the response properties of coupled subsystems, as well as the importance of accurate atomic polarizability models.

Role of metal-nanostructure features on tip-enhanced photoluminescence of single molecules

Tip-enhanced photoluminescence (TEPL) experiments have recently reached the ability to investigate single molecules exploiting resolution at the submolecular level. Localized surface plasmon resonances of metallic nanostructures have the capability of enhancing an impinging electromagnetic radiation in the proximity of their surface, with evident consequences both on absorption and emission of molecules placed in the same region. We propose a theoretical analysis of these phenomena in order to interpret TEPL experiments on single molecules, including a quantum mechanical description of the target molecule equilibrated with the presence of two nanostructures representative of the nanocavity usually employed in STMs. The approach has been applied to the zinc phthalocyanine molecule, previously considered in recent TEPL experiments [Yang et al., Nat. Photonics 14, 693-699 (2020)]. This work has the aim of providing a comprehensive theoretical understanding of the experimental results, particularly focusing on the investigation of the tip features that majorly influence the excitation and fluorescence processes of the molecule, such as the geometry, the dielectric function, and the tip-molecule distance.

Going Beyond the Limits of Classical Atomistic Modeling of Plasmonic Nanostructures

Theoretical modeling of plasmonic phenomena is of fundamental importance for rationalizing experimental measurements. Despite the great success of classical continuum modeling, recent technological advances allowing for the fabrication of structures defined at the atomic level require to be modeled through atomistic approaches. From a computational point of view, the latter approaches are generally associated with high computational costs, which have substantially hampered their extensive use. In this work, we report on a computationally fast formulation of a classical, fully atomistic approach, able to accurately describe both metal nanoparticles and graphene-like nanostructures composed of roughly 1 million atoms and characterized by structural defects.

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

Since surface-enhanced Raman scattering (SERS) is of considerable interest for sensing applications in aqueous solution, the role that solvent plays in the spectroscopy must be understood. However, these efforts are hindered due to a lack of simulation approaches for modeling solvent effects in SERS. In this work, we present an atomistic electrodynamics-quantum mechanical method to simulate SERS in aqueous solution based on the discrete interaction model/quantum mechanical method. This method combines an atomistic electrodynamics model of the nanoparticle with a time-dependent density functional theory description of the molecule and a polarizable embedding method for the solvent. The explicit treatment of solvent molecules and nanoparticles results in a large number of polarizable dipoles that need to be considered. To reduce the computational cost, a simple cut-off based approach has been implemented to limit the number of dipoles that need to be treated without sacrificing accuracy. As a test of this method, we have studied how solvent affects the SERS of pyridine in the junction between two nanoparticles in aqueous solution. We find that the solvent leads to an enhanced SERS due to an increased local field at the position of the pyridine. We further demonstrate the importance of both image field and local field effects in determining the enhancements and the spectral signatures. Our results show the importance of describing the local environment due to the solvent molecules when modeling SERS.

A polarizable three-layer frozen density embedding/molecular mechanics approach

We present a novel multilayer polarizable embedding approach in which the system is divided into three portions, two of which are treated using density functional theory and their interaction is based on frozen density embedding (FDE) theory, and both also mutually interact with a polarizable classical layer described using an atomistic model based on fluctuating charges (FQ). The efficacy of the model is demonstrated by extending the formalism to linear response properties and applying it to the simulation of the excitation energies of organic molecules in aqueous solution, where the solute and the first solvation shell are treated using FDE, while the rest of the solvent is modeled using FQ charges.

Hybrid theoretical models for molecular nanoplasmonics

The multidisciplinary nature of the research in molecular nanoplasmonics, i.e., the use of plasmonic nanostructures to enhance, control, or suppress properties of molecules interacting with light, led to contributions from different theory communities over the years, with the aim of understanding, interpreting, and predicting the physical and chemical phenomena occurring at molecular- and nano-scale in the presence of light. Multiscale hybrid techniques, using a different level of description for the molecule and the plasmonic nanosystems, permit a reliable representation of the atomistic details and of collective features, such as plasmons, in such complex systems. Here, we focus on a selected set of topics of current interest in molecular plasmonics (control of electronic excitations in light-harvesting systems, polaritonic chemistry, hot-carrier generation, and plasmon-enhanced catalysis). We discuss how their description may benefit from a hybrid modeling approach and what are the main challenges for the application of such models. In doing so, we also provide an introduction to such models and to the selected topics, as well as general discussions on their theoretical descriptions.

Plasmonic Resonances of Metal Nanoparticles: Atomistic vs. Continuum Approaches

The fully atomistic model, ωFQ, based on textbook concepts (Drude theory, electrostatics, quantum tunneling) and recently developed by some of the present authors in Nanoscale, 11, 6004-6015 is applied to the calculation of the optical properties of complex Na, Ag, and Au nanostructures. In ωFQ, each atom of the nanostructures is endowed with an electric charge that can vary according to the external electric field. The electric conductivity between nearest atoms is modeled by adopting the Drude model, which is reformulated in terms of electric charges. Quantum tunneling effects are considered by letting the dielectric response of the system arise from atom-atom conductivity. ωFQ is challenged to reproduce the optical response of metal nanoparticles of different sizes and shapes, and its performance is compared with continuum Boundary Element Method (BEM) calculations.

The chemical effect goes resonant - a full quantum mechanical approach on TERS

Lately, experimental evidence of unexpectedly extremely high spatial resolution of tip-enhanced Raman scattering (TERS) has been demonstrated. Theoretically, two different contributions are discussed: an electromagnetic effect, leading to a spatially confined near field due to plasmonic excitations; and the so-called chemical effect originating from the locally modified electronic structure of the molecule due to the close proximity of the plasmonic system. Most of the theoretical efforts have concentrated on the electromagnetic contribution or the chemical effect in case of non-resonant excitation. In this work, we present a fully quantum mechanical description including non-resonant and resonant chemical contributions as well as charge-transfer phenomena of these molecular-plasmonic hybrid systems at the density functional and the time-dependent density functional level of theory. We consider a surface-immobilized tin(ii) phthalocyanine molecule as the molecular system, which is minutely scanned by a plasmonic tip, modeled by a single silver atom. These different relative positions of the Ag atom to the molecule lead to pronounced alterations of the Raman spectra. These Raman spectra vary substantially, both in peak positions and several orders of magnitude in the intensity patterns under non-resonant and resonant conditions, and also, depending on, which electronic states are addressed. Our computational approach reveals that unique - non-resonant and resonant - chemical interactions among the tip and the molecule significantly alter the TERS spectra and are mainly responsible for the high, possibly sub-Angstrom spatial resolution.

Present and Future of Surface-Enhanced Raman Scattering

The discovery of the enhancement of Raman scattering by molecules adsorbed on nanostructured metal surfaces is a landmark in the history of spectroscopic and analytical techniques. Significant experimental and theoretical effort has been directed toward understanding the surface-enhanced Raman scattering (SERS) effect and demonstrating its potential in various types of ultrasensitive sensing applications in a wide variety of fields. In the 45 years since its discovery, SERS has blossomed into a rich area of research and technology, but additional efforts are still needed before it can be routinely used analytically and in commercial products. In this Review, prominent authors from around the world joined together to summarize the state of the art in understanding and using SERS and to predict what can be expected in the near future in terms of research, applications, and technological development. This Review is dedicated to SERS pioneer and our coauthor, the late Prof. Richard Van Duyne, whom we lost during the preparation of this article.

Resolving Molecular Structures with High-Resolution Tip-Enhanced Raman Scattering Images

Vibrational modes of a single molecule can be visualized by tip-enhanced Raman spectroscopy with atomic resolution. However, the exact vibrations associated with these Raman scattering images are still in debate due to the lack of theoretical interpretation. In this work, we systematically study the Raman scattering images of a single Co(II)-tetraphenylporphyrin molecule. The stable structure whose Raman scattering images consistently match experimental results is discovered. Furthermore, we elucidate the effects of near-field localizations and field gradient on the resolution in Raman scattering images. The approach of locally integrated Raman polarizability density employed in this work provides an intuitive explanation of the origin of the experimental Raman scattering images.

A classical picture of subnanometer junctions: an atomistic Drude approach to nanoplasmonics

The description of optical properties of subnanometer junctions is particularly challenging. Purely classical approaches fail, because the quantum nature of electrons needs to be considered. Here we report on a novel classical fully atomistic approach, ωFQ, based on the Drude model for conduction in metals, classical electrostatics and quantum tunneling. We show that ωFQ is able to reproduce the plasmonic behavior of complex metal subnanometer junctions with quantitative fidelity to full ab initio calculations. Besides the practical potentialities of our approach for large scale nanoplasmonic simulations, we show that a classical approach, in which the atomistic discretization of matter is properly accounted for, can accurately describe the nanoplasmonics phenomena dominated by quantum effects.

Theoretical Study of the Charge Transfer Exciton Binding Energy in Semiconductor Materials for Polymer:Fullerene-Based Bulk Heterojunction Solar Cells

Recent efforts and progress in polymer solar cell research have boosted the photovoltaic efficiency of the technology. This efficiency depends not only on the device architecture but also on the material properties. Thus, insight into the design of novel semiconductor materials is vital for the advancement of the field. This paper looks from a theoretical viewpoint into two of the factors for the design of semiconductor materials with applications to bulk heterojunction solar cells: the charge transfer exciton binding energy and the nanoscale arrangement of donor and acceptor molecules in blend systems. Being aware that the exciton dissociation of local excitons in charge transfer states initiates the charge generation process, the excited state properties of four oligomers (one donor-type: PEO-PPV; and three donor-acceptor-types: PTFB, PTB7, and PTB7-Th) and two fullerene derivatives ([60]PCBM and [70]PCBM), previously reported in the literature as having high electrical conductance, are studied. With such a study, the donor molecules, either of donor-type or donor-acceptor type, are screened as candidates for [60]PCBM- and/or [70]PCBM-based bulk heterojunctions. The charge transfer energy and charge transfer exciton binding energy of suitable donor:acceptor bulk heterojunctions, some of them not yet fabricated, are studied. Further, the charge transfer exciton binding energies of [60]PCBM- and [70]PCBM-based blends are compared. A combination of molecular dynamics simulations with calculations based on Kohn-Sham density functional theory (KS-DFT) and its time-dependent extension (KS-TDDFT) is used. An important feature of this work is that it incorporates the effect of the environment of the quantum chemical system in KS-DFT or KS-TDDFT calculations through a polarizable discrete reaction field (DRF). Our predictions in terms of the influence of the nanoscale arrangement of donor and acceptor molecules on the performance of organic solar cells indicate that bulk heterojunction morphologies for donor-acceptor-type oligomers lead to their lowest excited states having charge transfer character. Further, we find that in terms of favorable charge transfer exciton binding energy, the PTB7-Th:[70]PCBM blends outperform the other blends.

A Discrete Interaction Model/Quantum Mechanical Method for Simulating Plasmon-Enhanced Two-Photon Absorption

In this work, we extend the discrete interaction model/quantum mechanical (DIM/QM) method to simulate plasmon-enhanced two-photon absorption (PETPA). The metal nanoparticle is treated atomistically by means of electrodynamics, while the molecule is described using damped cubic response theory within a time-dependent density functional theory framework. Using DIM/QM, we study the PETPA of para-nitroaniline ( p-NA) with a focus on the local and image field effects, the molecular orientation effects, and the molecule-nanoparticle distance effects. Our findings show that the enhancement is more complex than the simple | E|⁴ enhancement mechanism, where | E| is the local field at the position of the molecule. Because of specific interactions with the nanoparticle, we find that a TPA dark state of p-NA can be significantly enhanced through a coupling with the plasmon excitation. The results presented in this work illustrate that the coupling between molecular excitations and plasmons can give rise to unusual and complex behavior in nonlinear spectroscopy that cannot simply be understood by considering the optical properties of the individual molecules and nanoparticles separately. The method presented here provides detailed insights into the enhancement of nonlinear optical properties of molecules coupled to plasmonic nanoparticles.

Real-Time Description of the Electronic Dynamics for a Molecule Close to a Plasmonic Nanoparticle

The optical properties of molecules close to plasmonic nanostructures greatly differ from their isolated molecule counterparts. To theoretically investigate such systems from a quantum-chemistry perspective, one has to take into account that the plasmonic nanostructure (e.g., a metal nanoparticle-NP) is often too large to be treated atomistically. Therefore, a multiscale description, where the molecule is treated by an ab initio approach and the metal NP by a lower level description, is needed. Here we present an extension of one such multiscale model [Corni, S.; Tomasi, J. J. Chem. Phys.2001, 114, 3739], originally inspired by the polarizable continuum model, to a real-time description of the electronic dynamics of the molecule and of the NP. In particular, we adopt a time-dependent configuration interaction (TD CI) approach for the molecule, the metal NP is described as a continuous dielectric of complex shape characterized by a Drude-Lorentz dielectric function, and the molecule-NP electromagnetic coupling is treated by an equation-of-motion (EOM) extension of the quasi-static boundary element method (BEM). The model includes the effects of both the mutual molecule-NP time-dependent polarization and the modification of the probing electromagnetic field due to the plasmonic resonances of the NP. Finally, such an approach is applied to the investigation of the light absorption of a model chromophore, LiCN, in the presence of a metal-NP of complex shape.

Conformational Contrast of Surface-Mediated Molecular Switches Yields Ångstrom-Scale Spatial Resolution in Ultrahigh Vacuum Tip-Enhanced Raman Spectroscopy

Tip-enhanced Raman spectroscopy (TERS) combines the ability of scanning probe microscopy (SPM) to resolve atomic-scale surface features with the single-molecule chemical sensitivity of surface-enhanced Raman spectroscopy (SERS). Here, we report additional insights into the nature of the conformational dynamics of a free-base porphyrin at room temperature adsorbed on a metal surface. We have interrogated the conformational switch between two metastable surface-mediated isomers of meso-tetrakis(3,5-ditertiarybutylphenyl)-porphyrin (H₂TBPP) on a Cu(111) surface. At room temperature, the barrier between the porphyrin ring buckled up/down conformations of the H₂TBPP-Cu(111) system is easily overcome, and a 2.6 Å lateral resolution by simultaneous TERS and STM analysis is achieved under ultrahigh vacuum (UHV) conditions. This work demonstrates the first UHV-TERS on Cu(111) and shows TERS can unambiguously distinguish the conformational differences between neighboring molecules with Ångstrom-scale spatial resolution, thereby establishing it as a leading method for the study of metal-adsorbate interactions.

Quantum mechanics capacitance molecular mechanics modeling of core-electron binding energies of methanol and methyl nitrite on Ag(111) surface

We study a newly devised quantum mechanics capacitance molecular mechanics (QMCMM) method for the calculation of core-electron binding energies in the case of molecules adsorbed on metal surfaces. This yet untested methodology is applied to systems with monolayer of methanol/methyl nitrite on an Ag(111) surface at 100 K temperature. It was found out that the studied C, N, and O 1s core-hole energies converge very slowly as a function of the radius of the metallic cluster, which was ascribed to build up of positive charge on the edge of the Ag slab. Further analysis revealed that an extrapolation process can be used to obtain binding energies that deviated less than 0.5 eV against experiments, except in the case of methanol O 1s where the difference was as large as 1.8 eV. Additional QM-cluster calculations suggest that the latter error can be connected to the lack of charge transfer over the QM-CMM boundary. Thus, the results indicate that the QMCMM and QM-cluster methods can complement each other in a holistic picture of molecule-adsorbate core-ionization studies, where all types of intermolecular interactions are considered.

Theory of Linear and Nonlinear Surface-Enhanced Vibrational Spectroscopies

The vibrational spectroscopy of molecules adsorbed on metal nanoparticles can be enhanced by many orders of magnitude so that the detection and identification of single molecules are possible. The enhancement of most linear and nonlinear vibrational spectroscopies has been demonstrated. In this review, we discuss theoretical approaches to understanding linear and nonlinear surface-enhanced vibrational spectroscopies. A unified description of enhancement mechanisms classified as either electromagnetic or chemical in nature is presented. Emphasis is placed on understanding the spectral changes necessary for interpretation of linear and nonlinear surface-enhanced vibrational spectroscopies.

Plasmon-enhanced two-photon-induced isomerization for highly-localized light-based actuation of inorganic/organic interfaces

Two-photon initiated photo-isomerization of an azobenzene moiety adsorbed on silver nanoparticles (Ag NPs) is demonstrated. The azobenzene is linked to a materials-binding peptide that brings it into intimate contact with the Ag NP surface, producing a dramatic enhancement of its two-photon absorbance. An integrated modeling approach, combining advanced conformational sampling with Quantum Mechanics/Capacitance Molecular Mechanics and response theory, shows that charge transfer and image charges in the Ag NP generate local fields that enhance two-photon absorption of the cis isomer, but not the trans isomer, of adsorbed molecules. Moreover, dramatic local field enhancement is expected near the localized surface plasmon resonance (LSPR) wavelength, and the LSPR band of the Ag NPs overlaps the azobenzene absorbance that triggers cis to trans switching. As a result, the Ag NPs enable two-photon initiated cis to trans isomerization, but not trans to cis isomerization. Confocal anti-Stokes fluorescence imaging shows that this effect is not due to local heating, while the quadratic dependence of switching rate on laser intensity is consistent with a two-photon process. Highly localized two-photon initiated switching could allow local manipulation near the focal point of a laser within a 3D nanoparticle assembly, which cannot be achieved using linear optical processes.

Toward an Enhancement of the Photoactivity of Multiphotochromic Dimers Using Plasmon Resonance: A Theoretical Study

Building dimers of organic photochromic compounds paves the way to multifunctional switches, but such architectures often undergo partial photoreactivity only. Combining photochromism of molecules and plasmon resonance of gold nanoparticles (NPs) is known to affect the photochromism of monomers, yet the impact on multimers remains unknown. Here we propose a theoretical study of dimers of dithienylethenes by the mean of a hybrid calculation scheme (discrete-interaction model/quantum mechanics). We aim to assess how the optical properties of multiphotochromes are tuned by the influence of the plasmon resonances. We show that, for a typical chemisorption orientation on the NP, the absorption bands responsible for the photochromism are significantly enhanced for both the doubly open and mixed closed-open isomers of the dyad, hinting that plasmon resonance could be used to boost the generally poor photoactivity of dithienylethene dyads.

Plasmonic circular dichroism of 310- and α-helix using a discrete interaction model/quantum mechanics method

Plasmonic circular dichroism (CD) of chiral molecules in the near field of plasmonic nanoparticles (NPs) may be used to enhance molecular CD signatures or to induce a CD signal at the plasmon resonance. A recent few-states theory explored these effects for model systems and showed an orientation dependence of the sign of the induced CD signal for spherical NPs. Here, we use the discrete interaction model/quantum mechanical (DIM/QM) method to simulate the CD and plasmonic CD of the 310- and α-helix conformations of a short alanine peptide. We find that the interactions between the molecule and the plasmon lead to significant changes in the CD spectra. In the plasmon region, we find that the sign of the CD depends strongly on the orientation of the molecule as well as specific interactions with the NP through image dipole effects. A small enhancement of the CD is found in the molecular region of the spectrum, however, the molecular signatures may be significantly altered through interactions with the NP. We also show that the image dipole effect can result in induced plasmonic CD even for achiral molecules. Overall, we find that the specific interactions with the NP can lead to large changes to the CD spectrum that complicates the interpretation of the results.

Dynamical coupling of plasmons and molecular excitations by hybrid quantum/classical calculations: time-domain approach

The presence of plasmonic material influences the optical properties of nearby molecules in untrivial ways due to the dynamical plasmon-molecule coupling. We combine quantum and classical calculation schemes to study this phenomenon in a hybrid system that consists of a Na(2) molecule located in the gap between two Au/Ag nanoparticles. The molecule is treated quantum-mechanically with time-dependent density-functional theory, and the nanoparticles with quasistatic classical electrodynamics. The nanoparticle dimer has a plasmon resonance in the visible part of the electromagnetic spectrum, and the Na(2) molecule has an electron-hole excitation in the same energy range. Due to the dynamical interaction of the two subsystems the plasmon and the molecular excitations couple, creating a hybridized molecular-plasmon excited state. This state has unique properties that yield e.g. enhanced photoabsorption compared to the freestanding Na(2) molecule. The computational approach used enables decoupling of the mutual plasmon-molecule interaction, and our analysis verifies that it is not legitimate to neglect the back coupling effect when describing the dynamical interaction between plasmonic material and nearby molecules. Time-resolved analysis shows nearly instantaneous formation of the coupled state, and provides an intuitive picture of the underlying physics.

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.

A periodic charge-dipole electrostatic model. II. A kinetic-exchange-correlation correction

We extend the periodic charge-dipole electrostatic model, see I. V. Bodrenko, M. Sierka, E. Fabiano, and F. Della Sala, J. Chem. Phys. 137, 134702 (2012), to include a kinetic-exchange-correlation (KXC) correction. The KXC correction is approximated by means of an extended-Hückel-type formula, it is exact in the infinite jellium model and it is also computationally efficient as it requires only the computation of overlap integrals. Tests on the linear response of silver slabs to an external electrostatic perturbation show that the KXC correction yields a very accurate description of induced dipole and of the whole induced charge density profile. We also show that the KXC parameters are quite transferable and related to the atomic polarizability.