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

Identifying Differences between Semiclassical and Full-Quantum Descriptions of Plexcitons

Strong light-matter coupling between molecules and plasmonic nanoparticles gives rise to new hybrid eigenstates of the coupled system, commonly referred to as polaritons or, more precisely, plexcitons. Over the past decade, it has been amply shown that molecular electron dynamics and photophysics can be drastically affected by such interactions, thus paving the way for light-induced control of molecular excited state properties and reactivity. Here, by combining the ab initio molecular description and classical or quantum modeling of arbitrarily shaped plasmonic nanostructures within the stochastic Schrödinger equation, we present two approaches, one semiclassical and one full-quantum, to follow in real time the electronic dynamics of plexcitons while realistically taking plasmonic dissipative losses into account. The full-quantum theory is compared with the semiclassical analogue under different interaction regimes, showing (numerically and theoretically) that even in the weak-field and weak-coupling limit a small-yet-observable difference arises.

GPAW: An open Python package for electronic structure calculations

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe-Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn-Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

Recent Advances in Real-Time Time-Dependent Density Functional Theory Simulations of Plasmonic Nanostructures and Plasmonic Photocatalysis

Plasmonic catalysis provides a possible means for driving chemical reactions under relatively mild conditions. Rational design of these systems is impeded by the difficulty in understanding the electron dynamics and their interplay with reactions. Real-time, time-dependent density functional theory (RT-TDDFT) can provide dynamic information on excited states in plasmonic systems, including those relevant to plasmonic catalysis, at time scales and length scales that are otherwise out of reach of many experimental techniques. Here, we discuss previous RT-TDDFT studies of plasmonic systems, focusing on recent work that gains insight into plasmonic catalysis. These studies provide insight into plasmon dynamics, including size effects and the role of specific electronic states. Further, these studies provide significant insight into mechanisms underlying plasmonic catalysis, showing the importance of charge transfer between metal and adsorbate states, as well as local field enhancement, in different systems.

A Full Quantum Mechanical Approach Assessing the Chemical and Electromagnetic Effect in TERS

Tip-enhanced Raman spectroscopy (TERS) is a valuable method for surface analysis with nanometer to angstrom-scale resolution; however, the accurate simulation of particular TERS signals remains a computational challenge. We approach this challenge by combining the two main contributors to plasmon-enhanced Raman spectroscopy and to the high resolution in TERS, in particular, the electromagnetic and the chemical effect, into one quantum mechanical simulation. The electromagnetic effect describes the sample's interaction with the strong, highly localized, and inhomogeneous electric fields associated with the plasmonic tip and is typically the thematic focus for most mechanistic studies. On the other hand, the chemical effect covers the different responses to the extremely close-range and highly position-sensitive chemical interaction between the apex tip atom(s) and the sample, and, as we could show in previous works, plays an often underestimated role. Starting from a (time-dependent) density functional theory description of the chemical model system, comprised of a tin(II) phthalocyanine sample molecule and a single silver atom as the tip, we introduce the electromagnetic effect through a series of static point charges that recreate the electric field in the vicinity of the plasmonic Ag nanoparticle. By scanning the tip over the molecule along a 3D grid, we can investigate the system's Raman response on each position for nonresonant and resonant illumination. Simulating both effects on their own already hints at the achievable signal enhancement and resolution, but the combination of both creates even stronger evidence that TERS is capable of resolving submolecular features.

Dimers of Plasmonic Nanocubes to Reach Single-Molecule Strong Coupling with High Emission Yields

Reaching reproducible strong coupling between a quantum emitter and a plasmonic resonator at room temperature, while maintaining high emission yields, would make quantum information processing with light possible outside of cryogenic conditions. We theoretically propose to exploit the high local curvatures at the tips of plasmonic nanocubes to reach Purcell factors of >10⁶ at visible frequencies, rendering single-molecule strong coupling more easily accessible than with the faceted spherical nanoparticles used in recent experimental demonstrations. In the case of gold nanocube dimers, we highlight a trade-off between coupling strength and emission yield that depends on the nanocube size. Electrodynamic simulations on silver nanostructures are performed using a realistic dielectric constant, as confirmed by scattering spectroscopy performed on single nanocubes. Dimers of silver nanocubes feature Purcell factors similar to those of gold while allowing emission yields of >60%, thus providing design rules for efficient strongly coupled hybrid nanostructures at room temperature.

DNA sequencing by Förster resonant energy transfer

We propose a new DNA sequencing concept based on nonradiative Förster resonant energy transfer (FRET) from a donor quantum dot (QD) to an acceptor molecule. The FRET mechanism combined with the nanopore-based DNA translocation is suggested as a novel concept for sequencing DNA molecules. A recently-developed hybrid quantum/classical method is employed, which uses time-dependent density functional theory and quasistatic finite difference time domain calculations. Due to the significant absorbance of DNA bases for photon energies higher than 4 eV, biocompatibility, and stability, we use Zinc-Oxide (ZnO) QD as a donor in the FRET mechanism. The most sensitivity for the proposed method to DNA is achieved for the Hoechst fluorescent-dye acceptor and 1 nm ZnO-QD. Results show that the insertion of each type of DNA nucleobases between the donor and acceptor changes the frequency of the emitted light from the acceptor molecule between 0.25 to 1.6 eV. The noise analysis shows that the method can determine any unknown DNA nucleobases if the signal-to-noise ratio is larger than 5 dB. The proposed concept and excellent results shed light on a new promising class of DNA sequencers.

Electronic Dynamics of a Molecular System Coupled to a Plasmonic Nanoparticle Combining the Polarizable Continuum Model and Many-Body Perturbation Theory

The efficiency of plasmonic metallic nanoparticles in harvesting and concentrating light energy in their proximity triggers a wealth of important and intriguing phenomena. For example, spectroscopies are able to reach single-molecule and intramolecule sensitivities, and important chemical reactions can be effectively photocatalyzed. For the real-time description of the coupled dynamics of a molecule's electronic system and of a plasmonic nanoparticle, a methodology has been recently proposed (J. Phys. Chem. C. 120, 2016, 28774-28781) which combines the classical description of the nanoparticle as a polarizable continuum medium with a quantum-mechanical description of the molecule treated at the time-dependent configuration interaction (TDCI) level. In this work, we extend this methodology by describing the molecule using many-body perturbation theory: the molecule's excitation energies, transition dipoles, and potentials computed at the GW/Bethe-Salpeter equation (BSE) level. This allows us to overcome current limitations of TDCI in terms of achievable accuracy without compromising on the accessible molecular sizes. We illustrate the developed scheme by characterizing the coupled nanoparticle/molecule dynamics of two prototype molecules, LiCN and p-nitroaniline.

Interband plasmon-enhanced optical absorption of DNA nucleobases through the graphene nanopore

We propose a novel, to the best of our knowledge, plasmonic-based methodology for the purpose of fast DNA sequencing. The interband surface plasmon resonance and field-enhancement properties of graphene nanopore in the presence of the DNA nucleobases are investigated using a hybrid quantum/classical method (HQCM), which employs time-dependent density functional theory and a quasistatic finite difference time domain approach. In the strong plasmonic-molecular coupling regime where the plasmon and DNA absorption frequencies are degenerated, the optical response of DNA molecule in the vicinity of the nanopore is enhanced. In contrast, when the plasmon and nucleobases resonances are detuned the distinct peaks and broadening of the molecular resonances represent the inherent properties of the nucleobase. Due to the different optical properties of DNA nucleobases in the ultraviolet (UV) region of light, the signal corresponding to the replacement of nucleobases in a DNA block can be determined by considering the differential absorbance. Results show the promising capability of the present mechanism for practical DNA sequencing.

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.

Time-Resolved Excited-State Analysis of Molecular Electron Dynamics by TDDFT and Bethe-Salpeter Equation Formalisms

In this work, a theoretical and computational set of tools to study and analyze time-resolved electron dynamics in molecules, under the influence of one or more external pulses, is presented. By coupling electronic-structure methods with the resolution of the time-dependent Schrödinger equation, we developed and implemented the time-resolved induced density of the electronic wavepacket, the time-resolved formulation of the differential projection density of states (ΔPDOS), and of transition contribution map (TCM) to look at the single-electron orbital occupation and localization change in time. Moreover, to further quantify the possible charge transfer, we also defined the energy-integrated ΔPDOS and the fragment-projected TCM. We have used time-dependent density-functional theory (TDDFT), as implemented in ADF software, and the Bethe-Salpeter equation, as provided by MolGW package, for the description of the electronic excited states. This suite of postprocessing tools also provides the time evolution of the electronic states of the system of interest. To illustrate the usefulness of these postprocessing tools, excited-state populations have been computed for HBDI (the chromophore of GFP) and DNQDI molecules interacting with a sequence of two pulses. Time-resolved descriptors have been applied to study the time-resolved electron dynamics of HBDI, DNQDI, LiCN (being a model system for dipole switching upon highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) electronic excitation), and Ag22. The computational analysis tools presented in this article can be employed to help the interpretation of fast and ultrafast spectroscopies on molecular, supramolecular, and composite systems.

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.

A potential sensing mechanism for DNA nucleobases by optical properties of GO and MoS2 Nanopores

We propose a new DNA sensing mechanism based on optical properties of graphene oxide (GO) and molybdenum disulphide (MoS2) nanopores. In this method, GO and MoS2 is utilized as quantum dot (QD) nanopore and DNA molecule translocate through the nanopore. A recently-developed hybrid quantum/classical method (HQCM) is employed which uses time-dependent density functional theory and quasi-static finite difference time domain approach. Due to good biocompatibility, stability and excitation wavelength dependent emission behavior of GO and MoS2 we use them as nanopore materials. The absorption and emission peaks wavelengths of GO and MoS2 nanopores are investigated in the presence of DNA nucleobases. The maximum sensitivity of the proposed method to DNA is achieved for the 2-nm GO nanopore. Results show that insertion of DNA nucleobases in the nanopore shifts the wavelength of the emitted light from GO or MoS2 nanopore up to 130 nm. The maximum value of the relative shift between two different nucleobases is achieved by the shift between cytosine (C) and thymine (T) nucleobases, ~111 nm for 2-nm GO nanopore. Results show that the proposed mechanism has a superior capability to be used in future DNA sequencers.

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.