Semiclassical Real-Time Nuclear-Electronic Orbital Dynamics for Molecular Polaritons: Unified Theory of Electronic and Vibrational Strong Couplings

Nonadiabatic Hydrogen Tunneling Dynamics for Multiple Proton Transfer Processes with Generalized Nuclear-Electronic Orbital Multistate Density Functional Theory

Proton transfer and hydrogen tunneling play key roles in many processes of chemical and biological importance. The generalized nuclear-electronic orbital multistate density functional theory (NEO-MSDFT) method was developed in order to capture hydrogen tunneling effects in systems involving the transfer and tunneling of one or more protons. The generalized NEO-MSDFT method treats the transferring protons quantum mechanically on the same level as the electrons and obtains the delocalized vibronic states associated with hydrogen tunneling by mixing localized NEO-DFT states in a nonorthogonal configuration interaction scheme. Herein, we present the derivation and implementation of analytical gradients for the generalized NEO-MSDFT vibronic state energies and the nonadiabatic coupling vectors between these vibronic states. We use this methodology to perform adiabatic and nonadiabatic dynamics simulations of the double proton transfer reactions in the formic acid dimer and the heterodimer of formamidine and formic acid. The generalized NEO-MSDFT method is shown to capture the strongly coupled synchronous or asynchronous tunneling of the two protons in these processes. Inclusion of vibronically nonadiabatic effects is found to significantly impact the double proton transfer dynamics. This work lays the foundation for a variety of nonadiabatic dynamics simulations of multiple proton transfer systems, such as proton relays and hydrogen-bonding networks.

Nuclear-Electronic Orbital Time-Dependent Configuration Interaction Method

Combining real-time electronic structure with the nuclear-electronic orbital (NEO) method has enabled the simulation of complex nonadiabatic chemical processes. However, accurate descriptions of hydrogen tunneling and double excitations require multiconfigurational treatments. Herein, we develop and implement the real-time NEO time-dependent configuration interaction (NEO-TDCI) approach. Comparison to NEO-full CI calculations of absorption spectra for a molecular system shows that the NEO-TDCI approach can accurately capture the tunneling splitting associated with the electronic ground state as well as vibronic progressions corresponding to double electron-proton excitations associated with excited electronic states. Both of these features are absent from spectra obtained with single reference real-time NEO methods. Our simulations of hydrogen tunneling dynamics illustrate the oscillation of the proton density from one side to the other via a delocalized, bilobal proton wave function. These results indicate that the NEO-TDCI approach is highly suitable for studying hydrogen tunneling and other inherently multiconfigurational systems.

Real-Time Time-Dependent Density Functional Theory for Simulating Nonequilibrium Electron Dynamics

The explicit real-time propagation approach for time-dependent density functional theory (RT-TDDFT) has increasingly become a popular first-principles computational method for modeling various time-dependent electronic properties of complex chemical systems. In this Perspective, we provide a nontechnical discussion of how this first-principles simulation approach has been used to gain novel physical insights into nonequilibrium electron dynamics phenomena in recent years. Following a concise overview of the RT-TDDFT methodology from a practical standpoint, we discuss our recent studies on the electronic stopping of DNA in water and the Floquet topological phase as examples. Our discussion focuses on how RT-TDDFT simulations played a unique role in deriving new scientific understandings. We then discuss existing challenges and some new advances at the frontier of RT-TDDFT method development for studying increasingly complex dynamic phenomena and systems.

Squeezed Protons and Infrared Plasmonic Resonance Energy Transfer

Unusual nuclear quantum effects may emerge near noble metal nanostructures such as squeezed vibrational states in molecular junctions and plasmonic resonance energy transfer in the infrared domain. Herein, nuclear quantum effects near heavy metals are studied by nuclear-electronic orbital density functional theory (NEO-DFT) with an effective core potential. For a quantum proton sandwiched between a pair of gold tips modeled by two Au6 clusters, NEO-DFT calculations suggest that the quantum proton density can be squeezed as the tip distance decreases. For an HF molecule placed near a one-dimensional Au nanowire composed of up to 34 Au atoms, real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) shows that the infrared plasmonic motion within the Au nanowire may resonantly transfer electronic energy to the HF proton vibrational stretch mode. Overall, these calculations illustrate the advantages of the NEO approach for probing nuclear quantum effects, such as squeezed proton vibrational states and infrared plasmonic resonance energy transfer.

Simulating the Excited-State Dynamics of Polaritons with Ab Initio Multiple Spawning

Over the past decade, there has been a growth of interest in polaritonic chemistry, where the formation of hybrid light-matter states (polaritons) can alter the course of photochemical reactions. These hybrid states are created by strong coupling between molecules and photons in resonant optical cavities and can even occur in the absence of light when the molecule is strongly coupled with the electromagnetic fluctuations of the vacuum field. We present a first-principles model to simulate nonadiabatic dynamics of such polaritonic states inside optical cavities by leveraging graphical processing units (GPUs). Our first implementation of this model is specialized for a single molecule coupled to a single-photon mode confined inside the optical cavity but with any number of excited states computed using complete active space configuration interaction (CASCI) and a Jaynes-Cummings-type Hamiltonian. Using this model, we have simulated the excited-state dynamics of a single salicylideneaniline (SA) molecule strongly coupled to a cavity photon with the ab initio multiple spawning (AIMS) method. We demonstrate how the branching ratios of the photodeactivation pathways for this molecule can be manipulated by coupling to the cavity. We also show how one can stop the photoreaction from happening inside of an optical cavity. Finally, we also investigate cavity-based control of the ordering of two excited states (one optically bright and the other optically dark) inside a cavity for a set of molecules, where the dark and bright states are close in energy.

Quantum to Classical Cavity Chemistry Electrodynamics

Polaritonic chemistry has ushered in new avenues for controlling molecular dynamics. However, two key questions remain: (i) Can classical light sources elicit the same effects as certain quantum light sources on molecular systems? (ii) Can semiclassical treatments of light-matter interactions capture nontrivial quantum effects observed in molecular dynamics? This work presents a quantum-classical approach addressing issues of realizing cavity chemistry effects without actual cavities. It also highlights the limitations of the standard semiclassical light-matter interaction. It is demonstrated that classical light sources can mimic quantum effects up to the second order of light-matter interaction provided that the mean-field contribution, the symmetrized two-time correlation function, and the linear response function are the same in both situations. Numerical simulations show that the quantum-classical method aligns more closely with exact quantum molecular-only dynamics for quantum light states such as Fock states, superpositions of Fock states, and vacuum squeezed states than does the conventional semiclassical approach.

First-Principles Approach for Coupled Quantum Dynamics of Electrons and Protons in Heterogeneous Systems

The coupled quantum dynamics of electrons and protons is ubiquitous in many dynamical processes involving light-matter interaction, such as solar energy conversion in chemical systems and photosynthesis. A first-principles description of such nuclear-electronic quantum dynamics requires not only the time-dependent treatment of nonequilibrium electron dynamics but also that of quantum protons. Quantum mechanical correlation between electrons and protons adds further complexity to such coupled dynamics. Here we extend real-time nuclear-electronic orbital time-dependent density functional theory (RT-NEO-TDDFT) to periodic systems and perform first-principles simulations of coupled quantum dynamics of electrons and protons in complex heterogeneous systems. The process studied is an electronically excited-state intramolecular proton transfer of o-hydroxybenzaldehyde in water and at a silicon (111) semiconductor-molecule interface. These simulations illustrate how environments such as hydrogen-bonding water molecules and an extended material surface impact the dynamical process on the atomistic level. Depending on how the molecule is chemisorbed on the surface, excited-state electron transfer from the molecule to the semiconductor surface can inhibit ultrafast proton transfer within the molecule. This Letter elucidates how heterogeneous environments influence the balance between the quantum mechanical proton transfer and excited electron dynamics. The periodic RT-NEO-TDDFT approach is applicable to a wide range of other photoinduced heterogeneous processes.

Theory and modeling of light-matter interactions in chemistry: current and future

Light-matter interaction not only plays an instrumental role in characterizing materials' properties via various spectroscopic techniques but also provides a general strategy to manipulate material properties via the design of novel nanostructures. This perspective summarizes recent theoretical advances in modeling light-matter interactions in chemistry, mainly focusing on plasmon and polariton chemistry. The former utilizes the highly localized photon, plasmonic hot electrons, and local heat to drive chemical reactions. In contrast, polariton chemistry modifies the potential energy curvatures of bare electronic systems, and hence their chemistry, via forming light-matter hybrid states, so-called polaritons. The perspective starts with the basic background of light-matter interactions, molecular quantum electrodynamics theory, and the challenges of modeling light-matter interactions in chemistry. Then, the recent advances in modeling plasmon and polariton chemistry are described, and future directions toward multiscale simulations of light-matter interaction-mediated chemistry are discussed.

Polaritonic Chemistry: Hindering and Easing Ground State Polyenic Isomerization via Breakdown of σ-π Separation

The ground state conformational isomerization in polyenes is a symmetry allowed process. Its low energy barrier is governed by electron density transfer from the formal single bond that is rotated to the nearby formal double bonds. Along the reaction pathway, the transition state is therefore destabilized. The rules of polaritonic chemistry, i.e., chemistry in a nanocavity with reflecting windows, are barely beginning to be laid out. The standing electric field of the nanocavity couples strongly with the molecular wave function and modifies the potential energy curve in unexpected ways. A quantum electrodynamics approach, applied to the torsional degree of freedom of the central bond of butadiene, shows that formation of the polariton mixes the σ-π frameworks thereby stabilizing/destabilizing the planar, reactant-like conformations. The values of the fundamental mode of the cavity field used in the absence of the cavity do not trigger this mechanism.

Cavity Born-Oppenheimer Hartree-Fock Ansatz: Light-Matter Properties of Strongly Coupled Molecular Ensembles

Experimental studies indicate that optical cavities can affect chemical reactions through either vibrational or electronic strong coupling and the quantized cavity modes. However, the current understanding of the interplay between molecules and confined light modes is incomplete. Accurate theoretical models that take into account intermolecular interactions to describe ensembles are therefore essential to understand the mechanisms governing polaritonic chemistry. We present an ab initio Hartree-Fock ansatz in the framework of the cavity Born-Oppenheimer approximation and study molecules strongly interacting with an optical cavity. This ansatz provides a nonperturbative, self-consistent description of strongly coupled molecular ensembles, taking into account the cavity-mediated dipole self-energy contributions. To demonstrate the capability of the cavity Born-Oppenheimer Hartree-Fock ansatz, we study the collective effects in ensembles of strongly coupled diatomic hydrogen fluoride molecules. Our results highlight the importance of the cavity-mediated intermolecular dipole-dipole interactions, which lead to energetic changes of individual molecules in the coupled ensemble.

Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics

When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.

Nuclear-Electronic Orbital Quantum Dynamics of Plasmon-Driven H2 Photodissociation

Leveraging localized surface plasmon resonances of metal nanoparticles to trigger chemical reactions is a promising approach for heterogeneous catalysis. First-principles modeling of such processes is challenging due to the large number of electrons and electronic excited states as well as the significance of nuclear quantum effects when hydrogen is involved. Herein, the nonadiabatic nuclear-electronic quantum dynamics of plasmon-induced H2 photodissociation near an Al13- cluster is simulated with real-time nuclear-electronic orbital time-dependent density functional theory (RT-NEO-TDDFT). This approach propagates the nonequilibrium quantum dynamics of both electrons and protons. The plasmonic oscillations are shown to inject hot electrons into the antibonding orbital of H2, thereby inducing H2 dissociation. The quantum mechanical treatment of the hydrogen nuclei leads to faster H2 photodissociation and slightly larger isotope effects. Analysis of the nonequilibrium electronic density suggests that these findings stem from enhanced excited-state electronic coupling between the plasmonic mode and the H2 antibonding orbital due to proton delocalization or zero-point energy effects. Given the low computational overhead for including nuclear quantum effects with the RT-NEO-TDDFT approach, this work paves the way for simulating nonadiabatic nuclear-electronic quantum dynamics in other plasmonic systems.

QED Theory for Controlling the Molecule-Cavity Interaction: From Solvable Analytical Models to Realistic Ones

The study of the interactions of chemical systems in a cavity and the ability to control the reactions inside the cavities become an evolving and hot field of research. Despite that, there is still a significant gap between experiment and theory. Herein, we aim to bridge this gap by starting with the analysis of solvable analytical models for reactions inside a cavity, then continuing to realistic models for many molecules inside a single mode and in a multimode cavity. In addition, we investigate different ways to control the strength of the molecule-cavity coupling term, which in turn allows controlling chemical reactions. Our analysis can benefit the development of ab initio computational methods to simulate molecular systems in polariton cavities; in addition, we show how to parameterize the model Hamiltonians in order to simulate a specific molecular system. Finally, we demonstrate the possibility of achieving isomerization, in case it is prohibited out of the cavity, by placing the reaction inside a cavity.

Chiral Polaritonics: Analytical Solutions, Intuition, and Use

Preferential selection of a given enantiomer over its chiral counterpart has become increasingly relevant in the advent of the next era of medical drug design. In parallel, cavity quantum electrodynamics has grown into a solid framework to control energy transfer and chemical reactivity, the latter requiring strong coupling. In this work, we derive an analytical solution to a system of many chiral emitters interacting with a chiral cavity similar to the widely used Tavis-Cummings and Hopfield models of quantum optics. We are able to estimate the discriminating strength of chiral polaritonics, discuss possible future development directions and exciting applications such as elucidating homochirality, and deliver much needed intuition to foster the newly flourishing field of chiral polaritonics.

Electronic Born-Oppenheimer approximation in nuclear-electronic orbital dynamics

Within the nuclear-electronic orbital (NEO) framework, the real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) approach enables the simulation of coupled electronic-nuclear dynamics. In this approach, the electrons and quantum nuclei are propagated in time on the same footing. A relatively small time step is required to propagate the much faster electronic dynamics, thereby prohibiting the simulation of long-time nuclear quantum dynamics. Herein, the electronic Born-Oppenheimer (BO) approximation within the NEO framework is presented. In this approach, the electronic density is quenched to the ground state at each time step, and the real-time nuclear quantum dynamics is propagated on an instantaneous electronic ground state defined by both the classical nuclear geometry and the nonequilibrium quantum nuclear density. Because the electronic dynamics is no longer propagated, this approximation enables the use of an order-of-magnitude larger time step, thus greatly reducing the computational cost. Moreover, invoking the electronic BO approximation also fixes the unphysical asymmetric Rabi splitting observed in previous semiclassical RT-NEO-TDDFT simulations of vibrational polaritons even for small Rabi splitting, instead yielding a stable, symmetric Rabi splitting. For the intramolecular proton transfer in malonaldehyde, both RT-NEO-Ehrenfest dynamics and its BO counterpart can describe proton delocalization during the real-time nuclear quantum dynamics. Thus, the BO RT-NEO approach provides the foundation for a wide range of chemical and biological applications.

Nonadiabatic Wave Packet Dynamics with Ab Initio Cavity-Born-Oppenheimer Potential Energy Surfaces

Strong coupling of molecules with quantized electromagnetic fields can reshape their potential energy surfaces by forming dressed states. In such a scenario, it is possible to manipulate the dynamics of the molecule and open new photochemical reaction pathways. A theoretical approach to describe such coupled molecular-photon systems is the Cavity-Born-Oppenheimer (CBO) approximation. Similarly to the standard Born-Oppenheimer (BO) approximation, the system is partitioned and the electronic part of the system is treated quantum mechanically. This separation leads to CBO surfaces that depend on both nuclear and photonic coordinates. In this work, we demonstrated, for two molecular examples, how the concept of the CBO approximation can be used to perform nonadiabatic wave packet dynamics of a coupled molecular-cavity system. The light-matter interaction is incorporated in the CBO surfaces and the associated nonadiabatic coupling elements. We show that molecular and cavity contributions can be treated on the same numerical footing. This approach gives a new perspective on the description of light-matter coupling in molecular systems.

Interference between Molecular and Photon Field-Mediated Electron Transfer Coupling Pathways in Cavities

Cavity polaritonics creates novel opportunities to direct chemical reactions. Electron transfer (ET) reactions are among the simplest reactions, and they underpin energy conversion. New strategies to manipulate and direct electron flow at the nanoscale are of particular interest in biochemistry, energy science, bioinspired materials science, and chemistry. We show that optical cavities can modulate electron transfer pathway interferences and ET rates in donor-bridge-acceptor (DBA) systems. We derive the rate for DBA electron transfer when the molecules are coupled to cavity modes, emphasizing novel cavity-induced pathway interferences with the molecular electronic coupling pathways, as these interferences allow a new kind of ET rate tuning. The interference between the cavity-induced coupling pathways and the intrinsic molecular coupling pathway is dependent on the cavity properties. Thus, manipulating the interference between the cavity-induced DA coupling and the bridge-mediated coupling offers an approach to direct and manipulate charge flow at the nanoscale.

Nuclear Gradient Expressions for Molecular Cavity Quantum ElectrodynamicsSimulations using Mixed Quantum-Classical Methods

We derive a rigorous nuclear gradient for a molecule-cavity hybrid system using the Quantum Electrodynamics Hamiltonian. We treat the electronic-photonic DOFs as the quantum subsystem, and the nuclei as the classical subsystem. Using the adiabatic basis for the electronic DOF and the Fock basis for the photonic DOF, and requiring the total energy conservation of this mixed quantum-classical system, we derived the rigorous nuclear gradient for the molecule-cavity hybrid system, which is naturally connected to the approximate gradient under the Jaynes-Cummings approximation. The nuclear gradient expression can be readily used in any mixed quantum-classical simulations and will allow one to perform the non-adiabatic on-the-fly simulation of polariton quantum dynamics. The theoretical developments in this work could significantly benefit the polariton quantum dynamics community with a rigorous nuclear gradient of the molecule-cavity hybrid system and have a broad impact on the future non-adiabatic simulations of polariton quantum dynamics.