The cavity electromagnetic field within the polarizable continuum model of solvation

Absorption Intensities of Organic Molecules from Electronic Structure Calculations versus Experiments: the Effect of Solvation, Method, Basis Set, and Transition Moment Gauge

Recently, we derived experimental oscillator strengths (OSs) from well-defined UV-visible absorption spectral peaks of 100 molecules in solution. Here, we focus on a subset of transitions with the highest reliability to further benchmark the OSs from several wave function methods and density functionals. We consider multiple basis sets, transition moment gauges (length, velocity, and mixed), and solvent corrections. Most transitions in the comparison set come from conjugated molecules and have π → π* character. We use an automated algorithm to assign computed transitions to experimental bands. OSs computed using the Tamm-Dancoff approximation (TDA), CIS, or EOM-CCSD exhibited a strong gauge dependence, which is diminished in linear response theories (TD-DFT, TD-HF, and to a smaller degree LR-CCSD). OSs calculated from TD-DFT with PCM solvent models are systematically larger than apparent OSs derived from experimental spectra. For example, fcomp from hybrid functionals and PCM have mean absolute errors that are ∼10% of n·fexp, where n is a solvent refractive index factor that arises from the energy flux of the radiation field in a dielectric (solvent). Theoretical cavity field corrections considering spherical cavities do not improve the agreement between computed and experimental data. Corrections that account for the molecular shape and the direction of transition dipole moments, or that explicitly account for the effect of solvent molecules on the local field, should be more appropriate.

Time-dependent surface-enhanced Raman scattering: A theoretical approach

A new procedure for computing the time-dependent Raman scattering of molecules in the proximity of plasmonic nanoparticles (NPs) is proposed, drawing inspiration from the pioneering Lee and Heller's theory. This strategy is based on a preliminary simulation of the molecular vibronic wavefunction in the presence of a plasmonic nanostructure and an incident light pulse. Subsequently, the Raman signal is evaluated through an inverse Fourier Transform of the coefficients' dynamics. Employing a multiscale approach, the system is treated by coupling the quantum mechanical description of the molecule with the polarizable continuum model for the NP. This method offers a unique advantage by providing insights into the time evolution of the plasmon-enhanced Raman signal, tracking the dynamics of the incident electric field. It not only provides for the total Raman signal at the process's conclusion but also gives transient information. Importantly, the flexibility of this approach allows for the simulation of various incident electric field profiles, enabling a closer alignment with experimental setups. This adaptability ensures that the method is relevant and applicable to diverse real-world scenarios.

Hybrid quantum-classical polarizability model for single molecule biosensing

Optical whispering gallery mode biosensors are able to detect single molecules through effects of their polarizability. We address the factors that affect the polarizability of amino acids, which are the building blocks of life, via electronic structure theory. Amino acids are detected in aqueous environments, where their polarizability is different compared to the gasphase due to solvent effects. Solvent effects include structural changes, protonation and the local field enhancement through the solvent (water). We analyse the impact of these effects and find that all contribute to an increased effective polarizability in the solvent. We also address the excess polarizability relative to the displaced water cavity and develop a hybrid quantum-classical model that is in good agreement with self-consistent calculations. We apply our model to calculate the excess polarizability of 20 proteinogenic amino acids and determine the minimum resolution required to distinguish the different molecules and their ionised conformers based on their polarizability.

Perspectives on weak interactions in complex materials at different length scales

Nanocomposite materials consist of nanometer-sized quantum objects such as atoms, molecules, voids or nanoparticles embedded in a host material. These quantum objects can be exploited as a super-structure, which can be designed to create material properties targeted for specific applications. For electromagnetism, such targeted properties include field enhancements around the bandgap of a semiconductor used for solar cells, directional decay in topological insulators, high kinetic inductance in superconducting circuits, and many more. Despite very different application areas, all of these properties are united by the common aim of exploiting collective interaction effects between quantum objects. The literature on the topic spreads over very many different disciplines and scientific communities. In this review, we present a cross-disciplinary overview of different approaches for the creation, analysis and theoretical description of nanocomposites with applications related to electromagnetic properties.

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.

Simulating Plasmon Resonances of Gold Nanoparticles with Bipyramidal Shapes by Boundary Element Methods

Computational modeling and accurate simulations of localized surface plasmon resonance (LSPR) absorption properties are reported for gold nanobipyramids (GNBs), a class of metal nanoparticle that features highly tunable, geometry-dependent optical properties. GNB bicone models with spherical tips performed best in reproducing experimental LSPR spectra while the comparison with other geometrical models provided a fundamental understanding of base shapes and tip effects on the optical properties of GNBs. Our results demonstrated the importance of averaging all geometrical parameters determined from transmission electron microscopy images to build representative models of GNBs. By assessing the performances of LSPR absorption spectra simulations based on a quasi-static approximation, we provided an applicability range of this approach as a function of the nanoparticle size, paving the way to the theoretical study of the coupling between molecular electron densities and metal nanoparticles in GNB-based nanohybrid systems, with potential applications in the design of nanomaterials for bioimaging, optics and photocatalysis.

Mediation of resonance energy transfer by two polarisable particles

The molecular quantum electrodynamics theory is employed to calculate the matrix element and Fermi golden rule rate for resonant transfer of electronic excitation energy between a donor and an acceptor in the vicinity of two neutral electric dipole polarizable particles, which play the role of bridging species. The emitter and absorber couple linearly to the electric displacement field via their electric dipole moments, while each mediator interacts quadratically with this field through its dynamic polarizability. This form of interaction Hamiltonian enables fourth-order perturbation theory to be used to compute the probability amplitude together with summation over 24 time-ordered diagrams representing a single virtual photon exchange between each pair of coupled particles. Expressions for the migration rate mediated by two inert molecules are obtained for an arbitrary arrangement of the four species that are in fixed mutual orientation or are freely tumbling. These formulae are valid for all interparticle separation distances outside the orbital overlap region. From the general result, rate equations applicable to an equidistant collinear configuration of the four bodies are evaluated. Near- and far-zone limiting forms of the transfer rate for the relay pathway are also calculated and exhibit inverse sixth and inverse square dependences on relative separation distances between pairs of particles, confirming the short-range (radiationless) and long-range (radiative) energy transfer mechanisms associated with two-body theory. The distance behavior of interference terms between two-, three-, and four-body terms is also examined, and the relative importance of each contribution to the total transfer rate is discussed.

Quantum optimal control theory for solvated systems

In this work, we generalize the quantum optimal control theory (QOCT) of molecules subject to ultrashort laser pulses to the case of solvated systems, explicitly including the solvent dielectric properties in the system's quantum Hamiltonian. A reliable description of the solvent polarization is accounted for within the polarizable continuum model (PCM). The electron dynamics for the molecules in solution is coupled with the dynamics of the surrounding polarizable environment, which affects the features of the optimized laser pulse. To illustrate such effects, numerical applications of the developed method to the study of optimal population of selected excited states of two molecular solvated systems are presented and discussed.

Nonequilibrium Solvent Polarization Effects in Real-Time Electronic Dynamics of Solute Molecules Subject to Time-Dependent Electric Fields: A New Feature of the Polarizable Continuum Model

We develop an extension of the time-dependent equation-of-motion formulation of the polarizable continuum model (EOM-TDPCM) to introduce nonequilibrium cavity field effects in quantum mechanical calculations of solvated molecules subject to time-dependent electric fields. This method has been implemented in Octopus, a state-of-the-art code for real-space, real-time time-dependent density functional theory (RT-TDDFT) calculations. To show the potential of our methodology, we perform EOM-TDPCM/RT-TDDFT calculations of trans-azobenzene in water and in other model solvents with shorter relaxation times. Our results for the optical absorption spectrum of trans-azobenzene show (i) that cavity field effects have a clear impact in the overall spectral shape and (ii) that an accurate description of the solute shape (as the one provided within PCM) is key to correctly account for cavity field effects.

Perspective: Quantum Hamiltonians for optical interactions

The multipolar Hamiltonian of quantum electrodynamics is extensively employed in chemical and optical physics to treat rigorously the interaction of electromagnetic fields with matter. It is also widely used to evaluate intermolecular interactions. The multipolar version of the Hamiltonian is commonly obtained by carrying out a unitary transformation of the Coulomb gauge Hamiltonian that goes by the name of Power-Zienau-Woolley (PZW). Not only does the formulation provide excellent agreement with experiment, and versatility in its predictive ability, but also superior physical insight. Recently, the foundations and validity of the PZW Hamiltonian have been questioned, raising a concern over issues of gauge transformation and invariance, and whether observable quantities obtained from unitarily equivalent Hamiltonians are identical. Here, an in-depth analysis of theoretical foundations clarifies the issues and enables misconceptions to be identified. Claims of non-physicality are refuted: the PZW transformation and ensuing Hamiltonian are shown to rest on solid physical principles and secure theoretical ground.

A quantum-mechanical perspective on linear response theory within polarizable embedding

We present a derivation of linear response theory within polarizable embedding starting from a rigorous quantum-mechanical treatment of a composite system. To this aim, two different subsystem decompositions (symmetric and nonsymmetric) of the linear response function are introduced and the pole structures as well as residues of the individual terms are discussed. In addition to providing a thorough justification for the descriptions used in polarizable embedding models, this theoretical analysis clarifies which form of the response function to use and highlights complications in separating out subsystem contributions to molecular properties. The basic features of the presented expressions and various approximate forms are illustrated by their application to a composite model system.

Open-ended formulation of self-consistent field response theory with the polarizable continuum model for solvation

The study of high-order absorption properties of molecules is a field of growing importance. Quantum-chemical studies can help design chromophores with desirable characteristics. Given that most experiments are performed in solution, it is important to devise a cost-effective strategy to include solvation effects in quantum-chemical studies of these properties. We here present an open-ended formulation of self-consistent field (SCF) response theory for a molecular solute coupled to a polarizable continuum model (PCM) description of the solvent. Our formulation relies on the open-ended, density matrix-based quasienergy formulation of SCF response theory of Thorvaldsen, et al., [J. Chem. Phys., 2008, 129, 214108] and the variational formulation of the PCM, as presented by Lipparini et al., [J. Chem. Phys., 2010, 133, 014106]. Within the PCM approach to solvation, the mutual solute-solvent polarization is represented by means of an apparent surface charge (ASC) spread over the molecular cavity defining the solute-solvent boundary. In the variational formulation, the ASC is an independent, variational degree of freedom. This allows us to formulate response theory for molecular solutes in the fixed-cavity approximation up to arbitrary order and with arbitrary perturbation operators. For electric dipole perturbations, pole and residue analyses of the response functions naturally lead to the identification of excitation energies and transition moments. We document the implementation of this approach in the Dalton program package using a recently developed open-ended response code and the PCMSolver libraries and present results for one-, two-, three-, four- and five-photon absorption processes of three small molecules in solution.

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.

Local electric fields and molecular properties in heterogeneous environments through polarizable embedding

In spectroscopies, the local field experienced by a molecule embedded in an environment will be different from the externally applied electromagnetic field, and this difference may significantly alter the response and transition properties of the molecule. The polarizable embedding (PE) model has previously been developed to model the local field contribution stemming from the direct molecule-environment coupling of the electromagnetic response properties of molecules in solution as well as in heterogeneous environments, such as proteins. Here we present an extension of this approach to address the additional effective external field effect, i.e., the manifestations of the environment polarization induced by the external field, which allows for the calculation of properties defined in terms of the external field. Within a response framework, we report calculations of the one- and two-photon absorption (1PA and 2PA, respectively) properties of PRODAN-methanol clusters as well as the fluorescent protein DsRed. Our results demonstrate the necessity of accounting for both the dynamical reaction field and effective external field contributions to the local field in order to reproduce full quantum chemical reference calculations. For the lowest π→π* transition in DsRed, inclusion of effective external field effects gives rise to a 1.9- and 3.5-fold reduction in the 1PA and 2PA cross-sections, respectively. The effective external field is, however, strongly influenced by the heterogeneity of the protein matrix, and the resulting effect can lead to either screening or enhancement depending on the nature of the transition under consideration.

Gauge-invariant expectation values of the energy of a molecule in an electromagnetic field

In this paper, we show that the full Hamiltonian for a molecule in an electromagnetic field can be separated into a molecular Hamiltonian and a field Hamiltonian, both with gauge-invariant expectation values. The expectation value of the molecular Hamiltonian gives physically meaningful results for the energy of a molecule in a time-dependent applied field. In contrast, the usual partitioning of the full Hamiltonian into molecular and field terms introduces an arbitrary gauge-dependent potential into the molecular Hamiltonian and leaves a gauge-dependent form of the Hamiltonian for the field. With the usual partitioning of the Hamiltonian, this same problem of gauge dependence arises even in the absence of an applied field, as we show explicitly by considering a gauge transformation from zero applied field and zero external potentials to zero applied field, but non-zero external vector and scalar potentials. We resolve this problem and also remove the gauge dependence from the Hamiltonian for a molecule in a non-zero applied field and from the field Hamiltonian, by repartitioning the full Hamiltonian. It is possible to remove the gauge dependence because the interaction of the molecular charges with the gauge potential cancels identically with a gauge-dependent term in the usual form of the field Hamiltonian. We treat the electromagnetic field classically and treat the molecule quantum mechanically, but nonrelativistically. Our derivation starts from the Lagrangian for a set of charged particles and an electromagnetic field, with the particle coordinates, the vector potential, the scalar potential, and their time derivatives treated as the variables in the Lagrangian. We construct the full Hamiltonian using a Lagrange multiplier method originally suggested by Dirac, partition this Hamiltonian into a molecular term Hm and a field term Hf, and show that both Hm and Hf have gauge-independent expectation values. Any gauge may be chosen for the calculations; but following our partitioning, the expectation values of the molecular Hamiltonian are identical to those obtained directly in the Coulomb gauge. As a corollary of this result, the power absorbed by a molecule from a time-dependent, applied electromagnetic field is equal to the time derivative of the non-adiabatic term in the molecular energy, in any gauge.

Open-ended response theory with polarizable embedding: multiphoton absorption in biomolecular systems

We present the theory and implementation of an open-ended framework for electric response properties that includes effects from the molecular environment modeled by the polarizable embedding model.

Excited states in large molecular systems through polarizable embedding

Using the polarizable embedding model enables rational design of light-sensitive functional biological materials.

Stereo-electronic, vibrational, and environmental contributions to polarizabilities of large molecular systems: a feasible anharmonic protocol

Reliable computations of linear and non-linear optical properties of molecular systems in condensed phases require a proper account of stereo-electronic, vibrational, and environmental effects. In the framework of density functional theory, these effects can be accurately introduced using second-order vibrational perturbation theory in conjunction with polarizable continuum models. We illustrate the combination of an anharmonic description of the ground-state potential energy surface with solvation effects treated with the polarizable continuum model (PCM) in the calculation of the electronic, zero-point, and pure vibrational polarizabilities of selected systems. The description of the solvation environment is enriched by taking into account the dynamical aspects of the solute-solvent interactions through the inclusion of both electronic and vibrational non-equilbrium effects, as well as the direct effect of the solvent on the electric field that generates the molecular response (local field effect). This treatment yields accurate results which can be directly compared with experimental findings without the need of empirical corrections.