Search

We present an efficient set of methods for propagating excited-state dynamics involving a large number of configuration interaction singles (CIS) or Tamm-Dancoff approximation (TDA) single-reference excited states. Specifically, (i) following Head-Gordon et al., we implement an exact evaluation of the overlap of singly-excited CIS/TDA electronic states at different nuclear geometries using a biorthogonal basis and (ii) we employ a unified protocol for choosing the correct phase for each adiabat at each geometry. For many-electron systems, the combination of these techniques significantly reduces the computational cost of integrating the electronic Schrodinger equation and imposes minimal overhead on top of the underlying electronic structure calculation. As a demonstration, we calculate the electronic excited-state dynamics for a hydrogen molecule scattering off a silver metal cluster, focusing on high-lying excited states, where many electrons can be excited collectively and crossings are plentiful. Interestingly, we find that the high-lying, plasmon-like collective excitation spectrum changes with nuclear dynamics, highlighting the need to simulate non-adiabatic nuclear dynamics and plasmonic excitations simultaneously. In the future, the combination of methods presented here should help theorists build a mechanistic understanding of plasmon-assisted charge transfer and excitation energy relaxation processes near a nanoparticle or metal surface.

The accurate description of large molecular systems in complex environments remains an ongoing challenge for the field of computational chemistry. This problem is even more pronounced for photo-induced processes, as multiple excited electronic states and their corresponding non-adiabatic couplings must be taken into account. Multiscale approaches such as hybrid quantum mechanics/molecular mechanics (QM/MM) offer a balanced compromise between accuracy and computational burden. Here, we introduce an open-source software package (INAQS) for non-adiabatic QM/MM simulations that bridges the sampling capabilities of the GROMACS MD package and the excited-state infrastructure of the Q-CHEM electronic structure software. The interface is simple and can be adapted easily to other MD codes. The code supports a variety of different trajectory based molecular dynamics, ranging from Born-Oppenheimer to surface hopping dynamics. To illustrate the power of this combination, we simulate electronic absorption spectra, free energy surfaces along a reaction coordinate, and the excited state dynamics of 1,3-cyclohexadiene in solution., Comment: 45 pages, 7 figures

Nuclear Berry curvature effects emerge from electronic spin degeneracy and canlead to non-trivial spin-dependent (nonadiabatic) nuclear dynamics. However, such effects are completely neglected in all current mixed quantum-classical methods such as fewest switches surface-hopping. In this work, we present a phase-space surface-hopping (PSSH) approach to simulate singlet-triplet intersystem crossing dynamics. We show that with a simple pseudo-diabatic ansatz, a PSSH algorithm can capture the relevant Berry curvature effects and make predictions in agreement with exact quantum dynamics for a simple singlet-triplet model Hamiltonian. Thus, this approach represents an important step towards simulating photochemical and spin processes concomitantly, as relevant to intersystem crossing and spin-lattice relaxation dynamics.

We analyze a model problem representing a multi-electronic molecule sitting on a metal surface. Working with a reduced configuration interaction Hamiltonian, we show that one can extract very accurate ground state wavefunctions as compared with the numerical renormalization group theory (NRG) -- even in the limit of weak metal-molecule coupling strength but strong intramolecular electron-electron repulsion. Moreover, we extract what appear to be meaningful excitation energies as well. Our findings should lay the groundwork for future {\em ab initio} studies of charge transfer processes and bond making/breaking processes on metal surfaces.

We investigate a spin-boson inspired model of electron transfer, where the diabatic coupling is given by a position-dependent phase, eiWx. We consider both equilibrium and nonequilibrium initial conditions. We show that, for this model, all equilibrium results are completely invariant to the sign of W (to infinite order). However, the nonequilibrium results do depend on the sign of W, suggesting that photo-induced electron transfer dynamics with spin–orbit coupling can exhibit electronic spin polarization (at least for some time).

A quantum simulation of vibrational strong coupling (VSC) in the collective regime via thermostatted ring-polymer molecular dynamics (TRPMD) is reported. For a collection of liquid-phase water molecules resonantly coupled to a single lossless cavity mode, the simulation shows that, as compared with a fully classical calculation, the inclusion of nuclear and photonic quantum effects does not lead to a change in the Rabi splitting but does broaden polaritonic linewidths roughly by a factor of two. Moreover, under thermal equilibrium, both quantum and classical simulations predict that the static dielectric constant of liquid water is largely unchanged inside versus outside the cavity. This result disagrees with a recent experiment demonstrating that the static dielectric constant of liquid water can be resonantly enhanced under VSC, suggesting either limitations of our approach or perhaps other experimental factors that have not yet been explored., manuscript (18 pages) + supporting information (8 pages)

We simulate vibrational strong coupling (VSC) and vibrational ultrastrong coupling (V-USC) for liquid water with classical molecular dynamics simulations. When the cavity modes are resonantly coupled to the O−H stretch mode of liquid water, the infrared spectrum shows asymmetric Rabi splitting. The lower polariton (LP) may be suppressed or enhanced relative to the upper polariton (UP) depending on the frequency of the cavity mode. Moreover, although the static properties and the translational diffusion of water are not changed under VSC or V-USC, we do find the modification of the orientational autocorrelation function of H(2)O molecules especially under V-USC, which could play a role in ground-state chemistry.

We investigate two well-known approaches for extending the fewest switches surface hopping (FSSH) algorithm to periodic time-dependent couplings. The first formalism acts as if the instantaneous adiabatic electronic states were standard adiabatic states, which just happen to evolve in time. The second formalism replaces the role of the usual adiabatic states by the time-independent adiabatic Floquet states. For a set of modified Tully model problems, the Floquet FSSH (F-FSSH) formalism gives a better estimate for both transmission and reflection probabilities than the instantaneous adiabatic FSSH (IA-FSSH) formalism, especially for slow nuclear velocities. More importantly, only F-FSSH predicts the correct final scattering momentum. Finally, in order to use Floquet theory accurately, we find that it is crucial to account for the interference between wavepackets on different Floquet states. Our results should be of interest to all those interested in laser-induced molecular dynamics.

Dynamics at molecule-metal interfaces are a subject of intense current interest and come in many different flavors of experiments: gas-phase scattering, chemisorption, electrochemistry, nanojunction transport, and heterogeneous catalysis, to name a few. These dynamics involve nuclear degrees of freedom entangled with many electronic degrees of freedom (in the metal), and as such there is always the possibility for nonadiabatic phenomena to appear: the nuclei do not necessarily need to move slower than the electrons to break the Born-Oppenheimer (BO) approximation. In this Feature Article, we review a set of dynamical methods developed recently to deal with such nonadiabatic phenomena at a metal surface, methods that serve as alternatives to Tully's independent electron surface hopping (IESH) model. In the weak molecule-metal coupling regime, a classical master equation (CME) can be derived and a simple surface hopping approach is proposed to propagate nuclear and electronic dynamics stochastically. In the strong molecule-metal interaction regime, a Fokker-Planck equation can be derived for the nuclear dynamics, with electronic DoFs incorporated into the overall friction and random force. Lastly, a broadened classical master equation (BCME) can interpolate between the weak and strong molecule-metal interactions. Here, we briefly review these methods and the relevant benchmarking data, showing in particular how the methods can be used to calculate nonequilibrium transport properties. We highlight several open questions and pose several avenues for future study.

We investigate spin-dependent electron transfer in the presence of a Duschinskii rotation. In particular, we propagate dynamics for a two-level model system for which spin-orbit coupling introduces an interstate coupling of the form $e^{iWx}$, which is both position(x)-dependent and complex-valued. We demonstrate that two-level systems coupled to Brownian oscillators with Duschinskii rotations (and thus entangled normal modes) can produce marked increases in transient spin polarization relative to two-level systems coupled to simple shifted harmonic oscillators. These conclusions should have significant relevance for modeling the effect of nuclear motion on chiral induced spin selectivity.

The interplay between static and dynamic disorder and collective optical response in molecular ensembles is an important characteristic of nanoplasmonic and nanophotonic molecular systems. Here we investigate the cooperative superradiant response of a molecular ensemble of quantum emitters under the influence of environmental disorder, including inhomogeneous broadening (as induced by static random distribution of the molecular transition frequencies) and motional narrowing (as induced by stochastic modulation of these excitation energies). The effect of inhomogeneous broadening is to destroy the coherence of the collective molecular excitation and suppress superradiant emission. However, fast stochastic modulation of the molecular excitation energy can effectively restore the coherence of the quantum emitters and lead to a recovery of superradiant emission, which is an unexpected manifestation of motional narrowing. For a light scattering process as induced by an off-resonant incident pulse, stochastic modulation leads to inelastic fluorescence emission at the average excitation energy at long times and suggests that dynamic disorder effects can actually lead to collective excitation of the molecular ensemble.

We present a preliminary surface-hopping approach for modeling intersystem crossing (ISC) dynamics between four electronic states: one singlet and one (triply degenerate) triplet. In order to incorporate all Berry force effects, the algorithm requires that, when moving along an adiabatic surface associated with the triplet manifold, \mycomment{one must also keep track of a quasi-diabatic index (akin to a "$m_s$" quantum number) for each trajectory. For a simple model problem, we find that a great deal of new physics can be captured by our algorithm, setting the stage for larger, more realistic (or perhaps even {\em ab initio}) simulations in the future.

Chemical relaxation phenomena, including photochemistry and electron transfer processes, form a vigorous area of research in which nonadiabatic dynamics plays a fundamental role. However, for electronic systems with spin degrees of freedom, there are few if any applicable and practical quasiclassical methods. Here, we show that for nonadiabatic dynamics with two electronic states and a complex-valued Hamiltonian that does not obey time-reversal symmetry (as relevant to many coupled nuclear-electronic-spin systems), the optimal semiclassical approach is to generalize Tully’s surface hopping dynamics from coordinate space to phase space. In order to generate the relevant phase-space adiabatic surfaces, one isolates a proper set of diabats, applies a phase gauge transformation, and then diagonalizes the total Hamiltonian (which is now parameterized by both R and P). The resulting algorithm is simple and valid in both the adiabatic and nonadiabatic limits, incorporating all Berry curvature effects. Most importantly, the resulting algorithm allows for the study of semiclassical nonadiabatic dynamics in the presence of spin–orbit coupling and/or external magnetic fields. One expects many simulations to follow as far as modeling cutting-edge experiments with entangled nuclear, electronic, and spin degrees of freedom, e.g., experiments displaying chiral-induced spin selectivity.

We investigate a simple and robust scheme for choosing the phases of adiabatic electronic states smoothly (as a function of geometry) so as to maximize the performance of ab initio non-adiabatic dynamics methods. Our approach is based upon consideration of the overlap matrix (U) between basis functions at successive points in time and selecting the phases so as to minimize the matrix norm of log(U). In so doing, one can extend the concept of parallel transport to cases with sharp curve crossings. We demonstrate that this algorithm performs well under extreme situations where dozens of states cross each other either through trivial crossings (where there is zero effective diabatic coupling), or through non-trivial crossings (when there is a non-zero diabatic coupling), or through a combination of both. In all cases, we compute the time-derivative coupling matrix elements (or equivalently non-adiabatic derivative coupling matrix elements) that are as smooth as possible. Our results should be of interest to all who are interested in either non-adiabatic dynamics, or more generally, parallel transport in large systems.

Chemical manifestations of strong light–matter coupling have recently been a subject of intense experimental and theoretical studies. Here we review the present status of this field. Section 1 is an introduction to molecular polaritonics and to collective response aspects of light–matter interactions. Section 2 provides an overview of the key experimental observations of these effects, while Section 3 describes our current theoretical understanding of the effect of strong light–matter coupling on chemical dynamics. A brief outline of applications to energy conversion processes is given in Section 4. Pending technical issues in the construction of theoretical approaches are briefly described in Section 5. Finally, the summary in Section 6 outlines the paths ahead in this exciting endeavor.

We explore the effects of spin-orbit coupling on nuclear wave packet motion near an out-of-equilibrium molecular junction, where nonzero Berry curvature emerges as the antisymmetric part of the electronic friction tensor. The existence of nonzero Berry curvature mandates that different nuclear wave packets (associated with different electronic spin states) experience different nuclear Berry curvatures, i.e. different pseudo-magnetic fields. Furthermore, for a generic, two-orbital two-lead model (representing the simplest molecular junction), we report significant spin polarization of the {\em electronic} current with decaying and oscillating signatures in the large voltage limit -- all as a result of {\em nuclear} motion. These results are consistent with magnetic AFM chiral-induced spin selectivity experiments. Altogether, our results highlight an essential role for Berry curvature in condensed phase dynamics, where spin separation survives dissipation to electron-hole pair creation and emerges as one manifestation of nuclear Berry curvature., Main text: 6 pages, 3 figures; Supplementary Material: 10 pages, 3 figures

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design. This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

We isolate and dissect exactly how and why Tully's fewest switches surface hopping (FSSH) algorithm fails when two wave packets come together at a crossing and "recohere". Using two different one-dimensional avoided crossing models and an initial wave function, which is a superposition of wave packets on different adiabats, we show that the key failures pertain to asymptotic nuclear momenta rather than electronic populations. Moreover, these FSSH failures stem from the fundamental assumption of independent trajectories with time-local hopping. As such, there is no possible means to correct FSSH without either (i) introducing time-nonlocal dynamics (i.e., allowing trajectories to move forward and backward in time) or (ii) requiring that trajectories interact.

A simple combination of density functional theory/time-dependent density functional theory (DFT/TDDFT) and configuration interaction is presented to fix the incorrect topology of the S0- S1 conical intersection (CI) and allow a description of bond making and bond breaking in photoinduced dynamics. The proposed TDDFT-1D method includes one lone optimized doubly excited configuration in addition to the DFT/TDDFT singly excited states within the context of a large configuration interaction Hamiltonian. Results for ethylene and stilbene are provided to demonstrate that this ansatz can yield physically meaningful potential energy surfaces near S0- S1 avoided crossings without changing the vertical excitation energies far from the relevant crossings. We also investigate the famous linear water example to show that the algorithm calculates the correct topology of the S0- S1 CI and yields the correct geometric phase.

We analyze different stochastic approaches for simulating electron transfer and potential sweep experiments in diffusive multidimensional electrochemical systems. In particular, we focus on a simple two-dimensional system, where one dimension is a traditional mass diffusion coordinate moving reactants from bulk solution to an electrode, and a second dimension represents a reorganization coordinate capturing solvent motion. We find that this multidimensional system can indeed be reduced to a simpler one-dimensional model, provided certain circumstances pertaining to separability between the two coordinates are met. Our results should begin to bridge the gap between continuum models of electrochemical dynamics/transport and ab initio models of surface electronic structure.

We show that the Berry force as computed by an approximate, mean-field electronic structure can be meaningful if properly interpreted. In particular, for a model Hamiltonian representing a molecular system with an even number of electrons interacting via a two-body (Hubbard) interaction and a spin–orbit coupling, we show that a meaningful nonzero Berry force emerges whenever there is spin unrestriction—even though the Hamiltonian is real-valued and formally the on-diagonal single-surface Berry force must be zero. Moreover, if properly applied, this mean-field Berry force yields roughly the correct asymptotic motion for scattering through an avoided crossing. That being said, within the context of a ground-state calculation, several nuances do arise as far interpreting the Berry force correctly, and as a practical matter, the Berry force diverges near the Coulson–Fischer point (which can lead to numerical instabilities). We do not address magnetic fields here.

We present a computational approach to simulate linear sweep and cyclic voltammetry experiments that does not require a discretized grid in space to quantify diffusion. By using a Green’s function solution coupled to a standard implicit ordinary differential equation solver, we are able to simulate current and redox species concentrations using only a small grid in time. As a result, where benchmarking is possible, we find that the current method is faster than (and quantitatively identical to) established techniques. The present algorithm should help open the door for studying adsorption effects in inner sphere electrochemistry.

We examine the many open questions that arise for nonadiabatic dynamics in the presence of degenerate electronic states, e.g., for singlet-to-triplet intersystem crossing where a minimal Hamiltonian must include four states (two of which are always degenerate). In such circumstances, the standard surface hopping approach is not sufficient as the algorithm does not include Berry force. Yet, we hypothesize that such a Berry force may be crucial as far as creating chiral induced spin separation, which is now a burgeoning field of study. Thus, this Perspective highlights the fact that if one can generate a robust and accurate semiclassical approach for the case of degenerate states, one will take a big step forward toward merging chemical physics with spintronics.

For a small fraction of hot CO2 molecules immersed in a liquid-phase CO2 thermal bath, classical cavity molecular dynamics simulations show that forming collective vibrational strong coupling (VSC) between the C=O asymmetric stretch of CO2 molecules and a cavity mode accelerates hot-molecule relaxation. The physical mechanism underlying this acceleration is the fact that polaritons, especially the lower polariton, can be transiently excited during the nonequilibrium process, which facilitates intermolecular vibrational energy transfer. The VSC effects on these rates (i) resonantly depend on the cavity mode detuning, (ii) cooperatively depend on molecular concentration or Rabi splitting, and (iii) collectively scale with the number of hot molecules, which is similar to Dicke's superradiance. For larger cavity volumes, due to a balance between this superradiant-like behavior and a smaller light-matter coupling, the total VSC effect on relaxation rates can scale slower than $1/N$, and the average VSC effect per molecule can remain meaningful for up to $N \sim10^4$ molecules forming VSC. Moreover, we find that the transiently excited lower polariton prefers to relax by transferring its energy to the tail of the molecular energy distribution rather than equally distributing it to all thermal molecules. Finally, we highlight the similarities of parameter dependence between the current finding with VSC catalysis observed in Fabry-Perot microcavities., Comment: 10 pages, 7 figures

Selectively exciting target molecules to high vibrational states is inefficient in the liquid phase, which restricts the use of IR pumping to catalyze ground-state chemical reactions. Here, we demonstrate that this inefficiency can sometimes be solved by confining the liquid to an optical cavity under vibrational strong coupling conditions. For a liquid solution of 13CO2 solute in a 12CO2 solvent, cavity molecular dynamics simulations show that exciting a polariton (hybrid light-matter state) of the solvent with an intense laser pulse, under suitable resonant conditions, may lead to a very strong (> 3 quanta) and ultrafast (< 1 ps) excitation of the solute, even though the solvent ends up being barely excited. By contrast, outside a cavity the same input pulse fluence can excite the solute by only half a vibrational quantum and the selectivity of excitation is low. Our finding is robust under different cavity volumes, which may lead to observable cavity enhancement on IR photochemical reactions in Fabry-Pérot cavities., 10 pages of manuscript + 23 pages of SI

We analytically calculate the electronic friction tensor for a molecule near a metal surface in the case that the electronic Hamiltonian is complex-valued, e.g. the case that there is spin-orbit coupling and/or an external magnetic field. In such a case, even at equilibrium, we show that the friction tensor is not symmetric. Instead, the tensor is the real-valued sum of one positive definite tensor (corresponding to dissipation) plus one antisymmetric tensor (corresponding to a Berry pseudomagnetic force). Moreover, we find that this Berry force can be much larger than the dissipational force, suggesting the possibility of strongly spin-polarized chemicurrents or strongly spin-dependent rate constants for systems with spin-orbit coupling., Comment: 6 pages, 2 figures (main text); 15 pages, 15 figures (supplemental material)

We present an improved approach for generating a set of optimized frontier orbitals (HOMO and LUMO) that minimizes the energy of one double configuration. We further benchmark the effect of including such a double within a rigorous configuration interaction singles or a parameterized semi-empirical time-dependent density functional theory Hamiltonian for a set of test cases. Although we cannot quite achieve quantitative accuracy, the algorithm is quite robust and routinely delivers an enormous qualitative improvement to standard single-reference electronic structure calculations.

A new scheme is proposed for modeling molecular nonadiabatic dynamics near metal surfaces. The charge-transfer character of such dynamics is exploited to construct an efficient reduced representation for the electronic structure. In this representation, the fewest switches surface hopping (FSSH) approach can be naturally modified to include electronic relaxation (ER). The resulting FSSH-ER method is valid across a wide range of coupling strengths as supported by tests applied to the Anderson-Holstein model for electron transfer. Future work will combine this scheme with ab initio electronic structure calculations.

Within the context of very simple avoided crossing, we investigate the investigate the effect of a complex diabatic coupling in determining spin-dependent rate constants and scattering states. We find that, if the molecular geometry is not linear and the Berry force is not zero, one can find significant spin polarization of the products. This study emphasizes that, when analyzing nonadiabatic reactions with spin orbit coupling (and a complex Hamiltonian), one must consider how Berry force affects nuclear motion -- at least in the context of gas phase reactions. Work is currently ongoing as far as extrapolating these conclusions to the condensed phase where interesting spin selection has been observed in recent years., 52 pages, 11 figures

Though the concept of Berry force was proposed thirty years ago, little is known about the practical consequences of this force as far as chemical dynamics are concerned. Here, we report that when molecular dynamics pass near a conical intersection, a massive Berry force can appear as a result of even a small amount of spin-orbit coupling (, Spin polarization is at the basis of quantum information and underlies some natural processes, but many aspects still need to be explored. Here, the authors, by quantum mechanical computations, show that even a weak spin-orbit coupling near a conical intersection can induce large spin selection, with consequences for spin manipulation in photochemical or electrochemical reactions.

We present an algorithm for efficient calculation of analytic nonadiabatic derivative couplings between spin-adiabatic, time-dependent density functional theory states within the Tamm-Dancoff approximation. Our derivation is based on the direct differentiation of the Kohn-Sham pseudowavefunction using the framework of Ou et al. Our implementation is limited to the case of a system with an even number of electrons in a closed shell ground state, and we validate our algorithm against finite difference at an S1/T2 crossing of benzaldehyde. Through the introduction of a magnetic field spin-coupling operator, we break time-reversal symmetry to generate complex valued nonadiabatic derivative couplings. Although the nonadiabatic derivative couplings are complex valued, we find that a phase rotation can generate an almost entirely real-valued derivative coupling vector for the case of benzaldehyde. This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article appeared in N. Bellonzi et al., J. Chem. Phys. 152, 044112 (2020) and may be found at https://doi.org/10.1063/1.5126440.

Recent experiments suggest that vibrational strong coupling (VSC) may significantly modify ground-state chemical reactions and their rates even without external pumping. The intrinsic mechanism of this "vacuum-field catalysis" remains largely unclear. Generally, modifications of thermal reactions in the ground electronic states can be caused by equilibrium or non-equilibrium effects. The former are associated with modifications of the reactant equilibrium distribution as expressed by the transition state theory of chemical reaction rates, while the latter stem from the dynamics of reaching and leaving transition state configurations. Here, we examine how VSC can affect chemical reactions rates in a cavity environment according to transition state theory. Our approach is to examine the effect of coupling to cavity mode(s) on the potential of mean force (PMF) associated with the reaction coordinate. Within the context of classical nuclei and classical photons and also assuming no charge overlap between molecules, we find that while the PMF can be affected by the cavity environment, this effect is negligible for the usual micron-length cavities used to examine VSC situations.

Recent experiments have observed that the chemical and photophysical properties of molecules can be modified inside an optical Fabry-Perot microcavity under collective vibrational strong coupling (VSC) conditions, and such modification is currently not well understood by theory. In an effort to understand the origin of such cavity induced phenomena, some recent studies have focused on the effect of the cavity environment on the nonlinear optical response of the molecular subsystem. Here, we use a recently proposed protocol for classical cavity molecular dynamics (CavMD) simulations to numerically investigate the linear and nonlinear response of liquid carbon dioxide under such VSC conditions following an optical pulse excitation. We find that applying a strong pulse of excitation to the lower hybrid light-matter state, i.e., the lower polariton (LP), can lead to an overall molecular nonlinear absorption which is enhanced by up to two orders of magnitude relative to the excitation outside the cavity. This polariton-enhanced multiphoton absorption also causes an ultrashort LP lifetime (0.2 ps) under strong illumination. Unlike usual polariton relaxation processes -- whereby polaritonic energy transfers directly to the manifold of singly excited vibrational dark states -- under the present mechanism, the LP transfers energy directly to the manifold of higher vibrationally excited dark states; these highly excited dark states subsequently relax to the manifold of singly excited states with a lifetime of tens of ps. Because the present mechanism is generic in nature, we expect these numerical predictions to be experimentally observed in different molecular systems and in cavities with different volumes.

We calculate the excitation energies of finite 1D Hubbard chains with a variety of different site energies from two perspectives: (i) the physics-based Bethe-Salpeter equation (BSE) method and (ii) the chemistry-based configuration interaction (CI) approach. Results obtained from all methods are compared against the exact values for three classes of systems: metallic, impurity-doped, and molecular (semiconducting/insulating) systems. While in a previous study we showed that the GW method holds comparative advantages versus traditional quantum chemistry approaches for calculating the ionization potentials and electron affinities across a large range of Hamiltonians, we show now that the BSE method outperforms CI approaches only for metallic and semiconducting systems. For insulating molecular systems, CI approaches generate better results.

We investigate rates of electron transfer for generalized Anderson-Holstein models in the limit of weak molecule-metal coupling, using both surface hopping and electronic friction dynamics in one and two dimensions. Overall, provided there is an external source of friction, electronic friction can sometimes perform well even in the limit of small metal-molecule coupling and capture nonadiabatic effects. However, we show that electronic friction dynamics is likely to fail if there is a competition between nonequivalent pathways. Our conclusions provide further insight into the recent observation by Ouyang et al., [J. Chem. Theory Comput., 2016, 12, 4178] regarding the applicability of Kramer's theory in the adiabatic limit to recover Marcus theory in the nonadiabatic limit.

The standard fewest-switches surface hopping (FSSH) approach fails to model nonadiabatic dynamics when the electronic Hamiltonian is complex-valued and there are multiple nuclear dimensions; FSSH does not include geometric magnetic effects and does not have access to a gauge independent direction for momentum rescaling. In this paper, for the case of a Hamiltonian with two electronic states, we propose an extension of Tully's FSSH algorithm, which includes geometric magnetic forces and, through diabatization, establishes a well-defined rescaling direction. When combined with a decoherence correction, our new algorithm shows satisfying results for a model set of two-dimensional single avoided crossings.

It is well-known that under conditions of fast electronic equilibration and weak nonadiabaticity, nonadiabatic effects induced by electron-hole pair excitations can be partly incorporated through a frictional force. However, ab initio computation of the electronic friction tensor suffers from numerical instability and usually demands a convergence check. In this study, we present an efficient and accurate interpolation method for computing the electronic friction tensor in a nearly black-box manner as appropriate for molecular dynamics. In almost all cases, our method agrees quite well with the exact friction tensor which is available for several quadratic Hamiltonians. As such, we outperform more conventional approaches that are based on the introduction of a broadening parameter. Future work will implement this interpolation approach within ab initio software packages.

We investigate detailed balance for a quantum system interacting with thermal radiation within mixed quantum-classical theory. For a two-level system coupled to classical radiation fields, three semiclassical methods are benchmarked: (1) Ehrenfest dynamics overestimate the excited-state population at equilibrium due to the failure of capturing vacuum fluctuations. (2) The coupled Maxwell-Bloch equations, which supplement Ehrenfest dynamics by damping at the full golden rule rate, underestimate the excited state population due to double-counting of the self-interaction effect. (3) $\mathrm{Ehrenfest}+\mathrm{R}$ dynamics recover detailed balance and the correct thermal equilibrium by enforcing the correct balance between the optical excitation and spontaneous emission of the quantum system. These results highlight the fact that, when properly designed, mixed quantum-classical electrodynamics can simulate thermal equilibrium in the field of nanoplasmonics.

We present a preliminary extension of the fewest switches surface hopping (FSSH) algorithm to the case of complex Hamiltonians as appropriate for modeling the dynamics of photoexcited molecules in magnetic fields. We make ansatze for the direction of momentum rescaling and we account for Berry's phase effects through "magnetic" forces as applicable in the adiabatic limit. Because Berry's phase is a nonlocal, topological characteristic of a set of entangled potential energy surfaces, we find that Tully's local FSSH algorithm can only partially capture the correct physics., Comment: 33 pages, 10 figures

We re-investigate the famous Mollow triplet and show that most of the well-known quantum characteristics of the Mollow triplet--including incoherent emission and a non-standard dependence of the sidebands on detuning--can be recovered quantitatively using semiclassical dynamics with a classical light field. In fact, by not relying on the rotating wave approximation, a semiclassical model predicts some quantum effects beyond the quantum optical Bloch equation, including higher order scattering and asymmetric sideband features. This letter highlights the fact that, with strong intensities, many putatively quantum features of light-matter interactions arise from a simple balance of mean-field electrodynamics and elementary spontaneous emission which requires minimal computational cost. Our results suggest that the application of semiclassical electrodynamics to problems with strong light-matter coupling in the fields of nanophotonics and superradiance are likely to yield a plethora of new information., 3 figures

A highly efficient scheme is proposed and benchmarked to compute 2D optical spectra. This scheme is ideally designed for electronic spectroscopy; however, the method can be applied in a straightforward way to vibrational spectroscopy as well. Our scheme performs dynamics only for the t2 duration, eliminating explicit t1 and t3 coherent dynamics and thus can achieve dramatic improvements in efficiency. To gain this efficiency, we assume the system is in the inhomogeneous regime and that there is no significant nonadiabatic transfer of population during the t1 and t3 coherence times. Preliminary results are presented for the Frenkel Hamiltonian. We obtain excellent agreement with numerically exact results (which are possible for this simplistic model Hamiltonian), capturing all relevant trends at least qualitatively (and sometimes quantitatively).

Within a generalized Anderson-Holstein model, we investigate electron transfer rates using two different surface hopping algorithms: a broadened classical master equation (BCME) and independent electron surface hopping (IESH). We find that for large enough bandwidth and density of one electron states, and in the presence of external friction, the IESH results converge to the BCME results for impurity-bath model systems, recovering both relaxation rates and equilibrium populations. Without external friction, however, the BCME and IESH results can strongly disagree, and preliminary evidence suggests that IESH does not always recover the correct equilibrium state. Finally, we also demonstrate that adding an electronic thermostat to IESH does help drive the metallic substrate to the correct equilibrium state, but this improvement can sometimes come at the cost of worse short time dynamics. Overall, our results should be of use for all computational chemists looking to model either gas phase scattering or electrochemical dynamics at a metal interface.

The dynamics of an electronic system interacting with an electromagnetic field is investigated within mixed quantum–classical theory. Beyond the classical path approximation (where we ignore all feedback from the electronic system on the photon field), we consider all electron–photon interactions explicitly according to Ehrenfest (i.e., mean-field) dynamics and a set of coupled Maxwell–Liouville equations. Because Ehrenfest dynamics cannot capture certain quantum features of the photon field correctly, we propose a new Ehrenfest+R method that can recover (by construction) spontaneous emission while also distinguishing between electromagnetic fluctuations and coherent emission.

Interfacial electrochemistry and photo(electro)catalysis are key processes that convert the energy of photons or electrons to chemical bonds in many energy conversion and storage technologies. Achieving a molecular level understanding of the funda- mental interfacial structure, energetics, dynamics, and reaction mechanisms that govern these processes represents a broad frontier for chemical physics and physical chemistry. This Special Topic contains a collection of articles that range from the development of new experimental and computational techniques to the novel application of those techniques for mechanistic studies, as the principal investigators seek a fundamental molecular understanding of both electrode/electrolyte interfaces and the relevant electrocatalytic, photocatalytic, and photoelectrochemical reactions taking place thereabout. Altogether, this col- lection of articles captures the current state of this very active, frontier research field and highlights the current and remaining key scientific challenges and opportunities.

Semiclassical electrodynamics is an appealing approach for studying light-matter interactions, especially for realistic molecular systems. However, there is no unique semiclassical scheme. On the one hand, intermolecular interactions can be described instantaneously by static two-body interactions connecting different molecules plus a classical transverse E-field; we will call this Hamiltonian #I. On the other hand, intermolecular interactions can also be described as effects that are mediated exclusively through a classical one-body E-field without any quantum effects at all (assuming we ignore electronic exchange); we will call this Hamiltonian #II. Moreover, one can also mix these two Hamiltonians into a third, hybrid Hamiltonian, which preserves quantum electron-electron correlations for lower excitations but describes higher excitations in a mean-field way. To investigate which semiclassical scheme is most reliable for practical use, here we study the real-time dynamics of a pair of identical two-level systems (TLSs) undergoing either resonance energy transfer (RET) or collectively driven dynamics. While all approaches perform reasonably well when there is no strong external excitation, we find that no single approach is perfect for all conditions. Each method has its own distinct problems: Hamiltonian #I performs best for RET but behaves in a complicated manner for driven dynamics. Hamiltonian #II is always stable, but obviously fails for RET at short distances. One key finding is that, under externally driving, a full configuration interaction description of Hamiltonian #I strongly overestimates the long-time electronic energy, highlighting the not obvious fact that, if one plans to merge quantum molecules with classical light, a full, exact treatment of electron-electron correlations can actually lead to worse results than a simple mean-field treatment., Comment: 17 pages, 7 figures

The optical response of an electronic two-level system (TLS) coupled to an incident continuous wave (cw) electromagnetic (EM) field is simulated explicitly in one dimension by the following five approaches: (i) the coupled Maxwell-Bloch equations, (ii) the optical Bloch equation (OBE), (iii) Ehrenfest dynamics, (iv) the Ehrenfest+R approach and (v) classical dielectric theory (CDT). Our findings are as follows: (i) standard Ehrenfest dynamics predict the correct optical signals only in the linear response regime where vacuum fluctuations are not important; (ii) both the coupled Maxwell-Bloch equations and CDT predict incorrect features for the optical signals in the linear response regime due to a double-counting of self-interaction; (iii) by exactly balancing the effects of self-interaction versus the effects of quantum fluctuations (and insisting on energy conservation), the Ehrenfest+R approach generates the correct optical signals in the linear regime and slightly beyond, yielding, e.g., the correct ratio between the coherent and incoherent scattering EM fields. As such, Ehrenfest+R dynamics agree with dynamics from the quantum OBE, but whereas the latter is easily applicable only for a single TLS in vacuum, the former should be applicable to large systems in environments with arbitrary dielectrics. Thus, this benchmark study suggests that the Ehrenfest+R approach may be very advantageous for simulating light-matter interactions semiclassically., 3 figures