Investigating New Reactivities Enabled by Polariton Photochemistry

Cavity Controlled Upconversion in CdSe Nanoplatelet Polaritons

Exciton-polaritons provide a versatile platform for investigating quantum electrodynamics effects in chemical systems, such as polariton-altered chemical reactivity. However, using polaritons in chemical contexts will require a better understanding of their photophysical properties under ambient conditions, where chemistry is typically performed. Here, we used cavity quality factor to control strong light-matter interactions and in particular the excited state dynamics of colloidal CdSe nanoplatelets (NPLs) coupled to a Fabry-Pérot optical cavity. With increasing cavity quality factor, we observe significant population of the upper polariton (UP) state, exemplified by the rare observation of substantial UP photoluminescence (PL). Excitation of the lower polariton (LP) states results in upconverted PL emission from the UP branch due to efficient exchange of population between the LP, UP and the reservoir of dark states present in collectively coupled polaritonic systems. In addition, we measure time scales for polariton dynamics ∼100 ps, implying great potential for NPL based polariton systems to affect photochemical reaction rates. State-of-the-art quantum dynamical simulations show outstanding quantitative agreement with experiments, and thus provide important insight into polariton photophysical dynamics of collectively coupled nanocrystal-based systems. These findings represent a significant step toward the development of practical polariton photochemistry platforms.

Exciton energy transfer inside cavity-A benchmark study of polaritonic dynamics using the surface hopping method

Strong coupling between the molecular system and photon inside the cavity generates polaritons, which can alter reaction rates by orders of magnitude. In this work, we benchmark the surface hopping method to simulate non-adiabatic dynamics in a cavity. The comparison is made against a numerically exact method (the hierarchical equations of motion) for a model system investigating excitonic energy transfer for a broad range of parameters. Surface hopping captures the effects of the radiation mode well, both at resonance and off-resonance. We have further investigated parameters that can increase or decrease the rate of population transfer, and we find that surface hopping in general can capture both effects well. Finally, we show that the dipole self-energy term within our parameter regime does not significantly affect the system's dynamics.

Coherent control from quantum commitment probabilities

We introduce a general definition of a quantum committor in order to clarify reaction mechanisms and facilitate control in processes where coherent effects are important. With a quantum committor, we generalize the notion of a transition state to quantum superpositions and quantify the effect of interference on the progress of the reaction. The formalism is applicable to any linear quantum master equation supporting metastability for which absorbing boundary conditions designating the reactant and product states can be applied. We use this formalism to determine the dependence of the quantum transition state on coherences in a polaritonic system and optimize the initialization state of a conical intersection model to control reactive outcomes, achieving yields of the desired state approaching 100%. In addition to providing a practical tool, the quantum committor provides a conceptual framework for understanding reactions in cases when classical intuitions fail.

Coupling polyatomic molecules to lossy nanocavities: Lindblad vs Schrödinger description

The use of cavities to impact molecular structure and dynamics has become popular. As cavities, in particular plasmonic nanocavities, are lossy and the lifetime of their modes can be very short, their lossy nature must be incorporated into the calculations. The Lindblad master equation is commonly considered an appropriate tool to describe this lossy nature. This approach requires the dynamics of the density operator and is thus substantially more costly than approaches employing the Schrödinger equation for the quantum wave function when several or many nuclear degrees of freedom are involved. In this work, we compare numerically the Lindblad and Schrödinger descriptions discussed in the literature for a molecular example where the cavity is pumped by a laser. The laser and cavity properties are varied over a range of parameters. It is found that the Schrödinger description adequately describes the dynamics of the polaritons and emission signal as long as the laser intensity is moderate and the pump time is not much longer than the lifetime of the cavity mode. Otherwise, it is demonstrated that the Schrödinger description gradually fails. We also show that the failure of the Schrödinger description can often be remedied by renormalizing the wave function at every step of time propagation. The results are discussed and analyzed.

Impact of Cavity on Molecular Ionization Spectra

Ionization phenomena have been widely studied for decades. With the advent of cavity technology, the question arises how quantum light affects molecular ionization. As the ionization spectrum is recorded from the neutral ground state, it is usually possible to choose cavities which exert negligible effect on the neutral ground state, but have significant impact on the ion and the ionization spectrum. Particularly interesting are cases where the ion exhibits conical intersections between close-lying electronic states, which gives rise to substantial nonadiabatic effects. Assuming single-molecule strong coupling, we demonstrate that vibrational modes irrelevant in the absence of a cavity play a decisive role when the molecule is in the cavity. Here, dynamical symmetry breaking is responsible for the ion-cavity coupling and high symmetry enables control of the coupling via molecular orientation relative to the cavity field polarization. Significant impact on the spectrum by the cavity is found and shown to even substantially increase for less symmetric molecules.

Formulation of transition dipole gradients for non-adiabatic dynamics with polaritonic states

A general formulation of the strong coupling between photons confined in a cavity and molecular electronic states is developed for the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham method. The light-matter interaction is included in the Jaynes-Cummings model, which requires the derivation and implementation of the analytical derivatives of the transition dipole moments between the molecular electronic states. The developed formalism is tested in the simulations of the nonadiabatic dynamics in the polaritonic states resulting from the strong coupling between the cavity photon mode and the ground and excited states of the penta-2,4-dieniminium cation, also known as PSB3. Comparison with the field-free simulations of the excited-state decay dynamics in PSB3 reveals that the light-matter coupling can considerably alter the decay dynamics by increasing the excited state lifetime and hindering photochemically induced torsion about the C=C double bonds of PSB3. The necessity of obtaining analytical transition dipole gradients for the accurate propagation of the dynamics is underlined.

Generalized Born-Huang expansion under macroscopic quantum electrodynamics framework

Born-Huang expansion is the cornerstone for studying potential energy surfaces and non-adiabatic couplings (NACs) in molecular systems. However, the traditional approach is insufficient to describe the molecular system, which strongly interacts with quantum light. Inspired by the work by Schäfer et al., we develop the generalized Born-Huang expansion theory within a macroscopic quantum electrodynamics (QED) framework. The theory we present allows us to describe electromagnetic vacuum fluctuations in dielectric media and incorporate the effects of dressed photons (or polaritons) into NACs. With the help of the generalized Born-Huang expansion, we clearly classify electronic nuclear NACs, polaritonic nuclear NACs, and polaritonic electronic NACs. Furthermore, to demonstrate the advantage of the macroscopic QED framework, we estimate polaritonic electronic NACs without any free parameter, such as the effective mode volume, and demonstrate the distance dependence of the polaritonic electronic NACs in a silver planar system. In addition, we take a hydrogen atom in free space as an example and derive spontaneous emission rates from photonic electronic NACs (polaritonic electronic NACs are reduced to photonic electronic NACs). We believe that this work not only provides an avenue for the theoretical exploration of NACs in a nucleus-electron-polariton coupled system but also offers a more comprehensive understanding for molecules coupled with quantum light.

Controlling Molecular Photoisomerization in Photonic Cavities through Polariton Funneling

Strong coupling between photonic modes and molecular electronic excitations, creating hybrid light-matter states called polaritons, is an attractive avenue for controlling chemical reactions. Nevertheless, experimental demonstrations of polariton-modified chemical reactions remain sparse. Here, we demonstrate modified photoisomerization kinetics of merocyanine and diarylethene by coupling the reactant's optical transition with photonic microcavity modes. We leverage broadband Fourier-plane optical microscopy to noninvasively and rapidly monitor photoisomerization within microcavities, enabling systematic investigation of chemical kinetics for different cavity-exciton detunings and photoexcitation conditions. We demonstrate three distinct effects of cavity coupling: first, a renormalization of the photonic density of states, akin to a Purcell effect, leads to enhanced absorption and isomerization rates at certain wavelengths, notably red-shifting the onset of photoisomerization. This effect is present under both strong and weak light-matter couplings. Second, kinetic competition between polariton localization into reactive molecular states and cavity losses leads to a suppression of the photoisomerization yield. Finally, our key result is that in reaction mixtures with multiple reactant isomers, exhibiting partially overlapping optical transitions and distinct isomerization pathways, the cavity resonance can be tuned to funnel photoexcitations into specific reactant isomers. Thus, upon decoherence, polaritons localize into a chosen isomer, selectively triggering the latter's photoisomerization despite initially being delocalized across all isomers. This result suggests that careful tuning of the cavity resonance is a promising avenue to steer chemical reactions and enhance product selectivity in reaction mixtures.

Electron Transfer at Molecule-Metal Interfaces under Floquet Engineering: Rate Constant and Floquet Marcus Theory

Electron transfer (ET) at molecule-metal or molecule-semiconductor interfaces is a fundamental reaction that underlies all electrochemical processes and substrate-mediated surface photochemistry. In this study, we show that ET rates near a metal surface can be significantly manipulated by periodic driving (e.g., Floquet engineering). We employ the Floquet surface hopping and Floquet electronic friction algorithms developed previously to calculate the ET rates near the metal surface as a function of driving amplitudes and driving frequencies. We find that ET rates have a turnover effect when the driving frequencies increase. A Floquet Marcus theory is further formulated to analyze such a turnover effect. We then benchmark the Floquet Marcus theory against Floquet surface hopping and Floquet electronic friction methods, indicating that the Floquet Marcus theory works in the strong nonadiabatic regimes but fails in the weak nonadiabatic regimes. We hope these theoretical tools will be useful to study ET rates in the plasmonic cavity and plasmon-assisted photocatalysis.

Machine Learning Framework for Modeling Exciton Polaritons in Molecular Materials

A light-matter hybrid quasiparticle, called a polariton, is formed when molecules are strongly coupled to an optical cavity. Recent experiments have shown that polariton chemistry can manipulate chemical reactions. Polariton chemistry is a collective phenomenon, and its effects increase with the number of molecules in a cavity. However, simulating an ensemble of molecules in the excited state coupled to a cavity mode is theoretically and computationally challenging. Recent advances in machine learning (ML) techniques have shown promising capabilities in modeling ground-state chemical systems. This work presents a general protocol to predict excited-state properties, such as energies, transition dipoles, and nonadiabatic coupling vectors with the hierarchically interacting particle neural network. ML predictions are then applied to compute the potential energy surfaces and electronic spectra of a prototype azomethane molecule in the collective coupling scenario. These computational tools provide a much-needed framework to model and understand many molecules' emerging excited-state polariton chemistry.

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 Bottleneck in Colloidal Quantum Dots

Controlling the relaxation dynamics of excitons is key to improving the efficiencies of semiconductor-based applications. Confined semiconductor nanocrystals (NCs) offer additional handles to control the properties of excitons, for example, by changing their size or shape, resulting in a mismatch between excitonic gaps and phonon frequencies. This has led to the hypothesis of a significant slowing-down of exciton relaxation in strongly confined NCs, but in practice due to increasing exciton-phonon coupling and rapid multiphonon relaxation channels, the exciton relaxation depends only weakly on the size or shape. Here, we focus on elucidating the nonradiative relaxation of excitons in NCs placed in an optical cavity. We find that multiphonon emission of carrier governs the decay, resulting in a polariton-induced phonon bottleneck with relaxation time scales that are slower by orders of magnitude compared to the cavity-free case, while the photon fraction plays a secondary role.

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.

Understanding Polaritonic Chemistry from Ab Initio Quantum Electrodynamics

In this review, we present the theoretical foundations and first-principles frameworks to describe quantum matter within quantum electrodynamics (QED) in the low-energy regime, with a focus on polaritonic chemistry. By starting from fundamental physical and mathematical principles, we first review in great detail ab initio nonrelativistic QED. The resulting Pauli-Fierz quantum field theory serves as a cornerstone for the development of (in principle exact but in practice) approximate computational methods such as quantum-electrodynamical density functional theory, QED coupled cluster, or cavity Born-Oppenheimer molecular dynamics. These methods treat light and matter on equal footing and, at the same time, have the same level of accuracy and reliability as established methods of computational chemistry and electronic structure theory. After an overview of the key ideas behind those ab initio QED methods, we highlight their benefits for understanding photon-induced changes of chemical properties and reactions. Based on results obtained by ab initio QED methods, we identify open theoretical questions and how a so far missing detailed understanding of polaritonic chemistry can be established. We finally give an outlook on future directions within polaritonic chemistry and first-principles QED.

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.

On the Mechanism of Polaritonic Rate Suppression from Quantum Transition Paths

Polariton chemistry holds promise for facilitating mode-selective chemical reactions, but the underlying mechanism behind the rate modifications observed under strong vibrational coupling is not well understood. Using the recently developed quantum transition path theory, we have uncovered a mechanism of resonant suppression of a thermal reaction rate in a simple model polaritonic system consisting of a reactive mode in a bath confined to a lossless microcavity with a single photon mode. We observed the formation of a polariton during rate-limiting transitions on reactive pathways and identified the concomitant rate suppression as being due to hybridization between the reactive mode and the cavity mode, which inhibits bath-mediated tunneling. The transition probabilities that define the quantum master equation can be directly translated into a visualization of the corresponding polariton energy landscape. This landscape exhibits a double funnel structure with a large barrier between the initial and final states.

Etalon-Assisted Study of the Strong CO Ligand Vibrations of the fac-[Re(CO)3(bpy)(CH3CN)]+ Octahedral Complex

The strong CO ligand vibrations of an octahedral complex, fac-[Re (CO)₃(bpy)(CH₃CN)]⁺, in acetonitrile are observed at 2040 and 1932 cm^-1. Facial rhenium tricarbonyl systems offer very strong and isolated CO vibrations with the potential for interactions between these vibrations. This work first identifies the dominant ion-pair species using attenuated total reflection infrared (ATR-IR) absorption spectra on a dilution series and then determines the strength of these CO ligand vibrations (as isolated vibrations) with a combination of ATR-IR and etalon-based measurements that determine the absolute complex index of refraction of the solution. Finally, the etalon experiments are modeled to study the interaction between vibrations, which is a property not embedded in the solution's complex index of refraction. The ATR-IR spectra are accomplished on a dilution series as well as a larger set of spectra as these solutions evaporated. The A'(1) CO ligand band at 2040 cm^-1 is fit with a sum of three Lorentzian functions characterizing the distribution of free, solvent-separated, and contact ion pairs of this octahedral complex vs concentration. The other CO ligand band at 1932 cm^-1 is broader and complicated by the dynamics of vibrational interactions, the unresolved splitting of the A'(2) and A″ CO vibrations, and ion-pair speciation. The etalon transmission measurements vs angle were on a 0.029 M solution, and Rabi splittings of 19 and 38 cm^-1 were observed for the A'(1) CO vibration and the unresolved A'(2) + A″ CO vibrations, respectively. The great strength of the CO ligand vibrations is evident despite the use of a dilute solution. Integrated band intensities are reported in comparison to hybrid density functional calculations for isolated vibrations. Then, the observed Rabi splittings are modeled to obtain the coupling strength of the CO ligand vibration with etalon cavity modes and with each other. In summary, this work develops a method to determine the concentration of these solutions from the ATR-IR spectrum, characterizes the ion-pairing, shows that the index of refraction is not constant in the IR spectral region of interest, and develops an interaction Hamiltonian that characterizes cavity-vibration and vibration-vibration coupling.

Manipulating Attosecond Charge Migration in Molecules by Optical Cavities

The ultrafast electronic charge dynamics in molecules upon photoionization while the nuclear motions are frozen is known as charge migration. In a theoretical study of the quantum dynamics of photoionized 5-bromo-1-pentene, we show that the charge migration process can be induced and enhanced by placing the molecule in an optical cavity, and can be monitored by time-resolved photoelectron spectroscopy. The collective nature of the polaritonic charge migration process is investigated. We find that, unlike spectroscopy, molecular charge dynamics in a cavity is local and does not show many-molecule collective effects. The same conclusion applies to cavity polaritonic chemistry.

Ab Initio Molecular Cavity Quantum Electrodynamics Simulations Using Machine Learning Models

We present a mixed quantum-classical simulation of polariton dynamics for molecule-cavity hybrid systems. In particular, we treat the coupled electronic-photonic degrees of freedom (DOFs) as the quantum subsystem and the nuclear DOFs as the classical subsystem and use the trajectory surface hopping approach to simulate non-adiabatic dynamics among the polariton states due to the coupled motion of nuclei. We use the accurate nuclear gradient expression derived from the Pauli-Fierz quantum electrodynamics Hamiltonian without making further approximations. The energies, gradients, and derivative couplings of the molecular systems are obtained from the on-the-fly simulations at the level of complete active space self-consistent field (CASSCF), which are used to compute the polariton energies and nuclear gradients. The derivatives of dipoles are also necessary ingredients in the polariton nuclear gradient expression but are often not readily available in electronic structure methods. To address this challenge, we use a machine learning model with the Kernel ridge regression method to construct the dipoles and further obtain their derivatives, at the same level as the CASSCF theory. The cavity loss process is modeled with the Lindblad jump superoperator on the reduced density of the electronic-photonic quantum subsystem. We investigate the azomethane molecule and its photoinduced isomerization dynamics inside the cavity. Our results show the accuracy of the machine-learned dipoles and their usage in simulating polariton dynamics. Our polariton dynamics results also demonstrate the isomerization reaction of azomethane can be effectively tuned by coupling to an optical cavity and by changing the light-matter coupling strength and the cavity loss rate.

Non-Hermitian Hamiltonians for linear and nonlinear optical response: A model for plexcitons

In polaritons, the properties of matter are modified by mixing the molecular transitions with light modes inside a cavity. Resultant hybrid light-matter states exhibit energy level shifts, are delocalized over many molecular units, and have a different excited-state potential energy landscape, which leads to modified exciton dynamics. Previously, non-Hermitian Hamiltonians have been derived to describe the excited states of molecules coupled to surface plasmons (i.e., plexcitons), and these operators have been successfully used in the description of linear and third order optical response. In this article, we rigorously derive non-Hermitian Hamiltonians in the response function formalism of nonlinear spectroscopy by means of Feshbach operators and apply them to explore spectroscopic signatures of plexcitons. In particular, we analyze the optical response below and above the exceptional point that arises for matching transition energies for plasmon and molecular components and study their decomposition using double-sided Feynman diagrams. We find a clear distinction between interference and Rabi splitting in linear spectroscopy and a qualitative change in the symmetry of the line shape of the nonlinear signal when crossing the exceptional point. This change corresponds to one in the symmetry of the eigenvalues of the Hamiltonian. Our work presents an approach for simulating the optical response of sublevels within an electronic system and opens new applications of nonlinear spectroscopy to examine the different regimes of the spectrum of non-Hermitian Hamiltonians.

Dissociation slowdown by collective optical response under strong coupling conditions

We consider an ensemble of diatomic molecules resonantly coupled to an optical cavity under strong coupling conditions at normal incidence. Photodissociation dynamics is examined via direct numerical integration of the coupled Maxwell-Schrödinger equations with molecular rovibrational degrees of freedom explicitly taken into account. It is shown that the dissociation is significantly affected (slowed down) when the system is driven at its polaritonic frequencies. The observed effect is demonstrated to be of transient nature and has no classical analog. An intuitive explanation of the dissociation slowdown at polaritonic frequencies is proposed.

Etalon-Assisted Determination of the Complex Index of Refraction of a Solution for the Study of Strong Cavity-Vibrational Coupling of PF6- in Acetonitrile

A new method is established using an etalon cavity to assist in the determination of the wavelength-dependent complex index of refraction of a solution throughout the mid-infrared range. The results are used to study the cavity-vibration polaritons of PF₆^- in acetonitrile. Mixed states are formed by placing solution inside a pair of parallel plate mirrors with a wavelength-scale spacing, i.e., within an etalon, such that there are cavity states that are angle-tuned into resonance with the strong P-F vibrations. The dominant ν₃ vibrations of PF₆^- consist of nearly triply degenerate oscillations of the partial-positively charged phosphorous against antisymmetric concerted motions of different sets of fluorine atoms with partial negative charges. These vibrations are dominant even though the solute is 29 times less concentrated than the solvent on a molar basis. The first part of the paper describes the method of determining the complex index of refraction of the solution from a combination of etalon transmission maxima and the attenuated total reflection (ATR) absorption spectrum of the solution. The results are presented as an analytical function including a sum of 37 vibrational contributions. Absolute integrated isolated band intensities were determined to be 463 ± 4, 462 ± 7, and 266 ± 4 km/mol for the three ν₃ PF₆^- vibrations at 841.4, 847.4, and 854.0 cm^-1, respectively, which sum to 1191 ± 9 km/mol for the ν₃ band. Then, the results are used to simulate the measured etalon transmission using the transfer matrix (TM) method with and without the ν₃ target vibrations. The etalon transmission simulations reconstruct the position of cavity modes in the absence of target vibrations. They provide input data for the testing of simple quantum mechanical models for the interaction of vibrations with cavity modes and the interactions of vibrations with other vibrations within the molecule and between solute and solvent. The model shows that the nearly degenerate ν₃ vibrations interact with each other with a vibration-vibration coupling of 33 ± 5 cm^-1. This is comparable to the cavity-vibration coupling of 30.4 ± 2.9 cm^-1 of the two strongest vibrations of PF₆^-.

Polariton enhanced free charge carrier generation in donor-acceptor cavity systems by a second-hybridization mechanism

Cavity quantum electrodynamics has been studied as a potential approach to modify free charge carrier generation in donor-acceptor heterojunctions because of the delocalization and controllable energy level properties of hybridized light-matter states known as polaritons. However, in many experimental systems, cavity coupling decreases charge separation. Here, we theoretically study the quantum dynamics of a coherent and dissipative donor-acceptor cavity system, to investigate the dynamical mechanism and further discover the conditions under which polaritons may enhance free charge carrier generation. We use open quantum system methods based on single-pulse pumping to find that polaritons have the potential to connect excitonic states and charge separated states, further enhancing free charge generation on an ultrafast timescale of several hundred femtoseconds. The mechanism involves polaritons with optimal energy levels that allow the exciton to overcome the high Coulomb barrier induced by electron-hole attraction. Moreover, we propose that a second-hybridization between a polariton state and dark states with similar energy enables the formation of the hybrid charge separated states that are optically active. These two mechanisms lead to a maximum of 50% enhancement of free charge carrier generation on a short timescale. However, our simulation reveals that on the longer timescale of picoseconds, internal conversion and cavity loss dominate and suppress free charge carrier generation, reproducing the experimental results. Thus, our work shows that polaritons can affect the charge separation mechanism and promote free charge carrier generation efficiency, but predominantly on a short timescale after photoexcitation.

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.

Collective response in light-matter interactions: The interplay between strong coupling and local dynamics

Strong molecule-radiation field coupling is often reached when a large number N of molecules respond collectively to the radiation field. In electronic strong coupling, molecular nuclear dynamics following polariton excitation reflects (a) the timescale separation between the fast electronic and photonic dynamics and the slow nuclear motion on one hand, and (b) the interplay between the collective nature of the molecule-field coupling and the local nature of the molecules nuclear response on the other. The first implies that the electronic excitation takes place, in the spirit of the Born approximation, at an approximately fixed nuclear configuration. The second can be rephrased as the intriguing question, can the collective nature of the optical excitation lead to collective nuclear motion following polariton formation, resulting in so-called polaron decoupled dynamics. We address this issue by studying the dynamical properties of a simplified Holstein-Tavis-Cummings type model, in which boson modes representing molecular vibrations are replaced by two-level systems while the boson frequency and the vibronic coupling are represented by the coupling between these levels (that induces Rabi oscillations between them) and electronic state dependence of this coupling. We investigate the short-time behavior of this model following polariton excitation as well as its response to CW driving and its density of states spectrum. We find that, while some aspects of the dynamical behavior appear to adhere to the polaron decoupling picture, the observed dynamics mostly reflect the local nature of the nuclear configuration of the electronic polariton rather than this picture.

Incorporating Lindblad Decay Dynamics into Mixed Quantum-Classical Simulations

We derive the $\mathcal{L}$-MFE method to incorporate Lindblad jump operator dynamics into the mean-field Ehrenfest (MFE) approach. We map the density matrix evolution of Lindblad dynamics onto pure state coefficients using trajectory averages. We use simple assumptions to construct the $\mathcal{L}$-MFE method that satisfies this exact mapping. This establishes a method that uses independent trajectories which exactly reproduces Lindblad decay dynamics using a wavefunction description, with deterministic changes of the magnitudes of the quantum expansion coefficients, while only adding on a stochastic phase. We further demonstrate that when including nuclei in the Ehrenfest dynamics, the $\mathcal{L}$-MFE method gives semi-quantitatively accurate results, with the accuracy limited by the accuracy of the approximations present in the semiclassical MFE approach. This work provides a general framework to incorporate Lindblad dynamics into semiclassical or mixed quantum-classical simulations.

On the characteristic features of ionization in QED environments

The ionization of molecular systems is important in many chemical processes, such as electron transfer and hot electron injection. Strong coupling between molecules and quantized fields (e.g., inside optical cavities) represents a new promising way to modify molecular properties in a non-invasive way. Recently, strong light-matter coupling has shown the potential to significantly improve the rates of hot electron driven processes, for instance, in water splitting. In this paper, we demonstrate that inside an optical cavity, the residual interaction between an outgoing free electron and the vacuum field is significant. We further show that since the quantized field is also interacting with the ionized molecule, the free electron and the molecular system are correlated. We develop a theoretical framework to account for the field induced correlation and show that the interaction between the free electron and the field, free electron-field interaction, has sizable effects on the ionization potential of typical organic molecules.

Coupling, lifetimes, and "strong coupling" maps for single molecules at plasmonic interfaces

The interaction between excited states of a molecule and excited states of a metal nanostructure (e.g., plasmons) leads to hybrid states with modified optical properties. When plasmon resonance is swept through molecular transition frequency, an avoided crossing may be observed, which is often regarded as a signature of strong coupling between plasmons and molecules. Such strong coupling is expected to be realized when 2|⟨U⟩|/ℏΓ > 1, where ⟨U⟩ and Γ are the molecule-plasmon coupling and the spectral width of the optical transition, respectively. Because both ⟨U⟩ and Γ strongly increase with decreasing distance between a molecule and a plasmonic structure, it is not obvious that this condition can be satisfied for any molecule-metal surface distance. In this work, we investigate the behavior of ⟨U⟩ and Γ for several geometries. Surprisingly, we find that if the only contributions to Γ are lifetime broadenings associated with the radiative and nonradiative relaxation of a single molecular vibronic transition, including effects on molecular radiative and nonradiative lifetimes induced by the metal, the criterion 2|⟨U⟩|/ℏΓ > 1 is easily satisfied by many configurations irrespective of the metal-molecule distance. This implies that the Rabi splitting can be observed in such structures if other sources of broadening are suppressed. Additionally, when the molecule-metal surface distance is varied keeping all other molecular and metal parameters constant, this behavior is mitigated due to the spectral shift associated with the same molecule-plasmon interaction, making the observation of Rabi splitting more challenging.

Non-Hermitian Cavity Quantum Electrodynamics - Configuration Interaction Singles Approach for Polaritonic Structure with ab initio Molecular Hamiltonians

We combine ab initio molecular electronic Hamiltonians with a cavity quantum electrodynamics model for dissipative photonic modes and apply mean-field theories to the ground- and excited-states of resulting polaritonic systems. In particular, we develop a non-Hermitian configuration interaction singles theory for mean-field ground- and excited-states of the molecular system strongly interacting with a photonic mode, and apply these methods to elucidating the phenomenology of paradigmatic polaritonic systems. We leverage the Psi4Numpy framework to yield open-source and accessible reference implementations of these methods.

Cavity quantum-electrodynamical time-dependent density functional theory within Gaussian atomic basis. II. Analytic energy gradient

Following the formulation of cavity quantum-electrodynamical time-dependent density functional theory (cQED-TDDFT) models [Flick et al., ACS Photonics 6, 2757–2778 (2019) and Yang et al., J. Chem. Phys. 155, 064107 (2021)], here, we report the derivation and implementation of the analytic energy gradient for polaritonic states of a single photochrome within the cQED-TDDFT models. Such gradient evaluation is also applicable to a complex of explicitly specified photochromes or, with proper scaling, a set of parallel-oriented, identical-geometry, and non-interacting molecules in the microcavity.

An Accurate Linearized Semiclassical Approach for Calculating Cavity-Modified Charge Transfer Rate Constants

We show that combining the linearized semiclasscial approximation with Fermi's golden rule (FGR) rate theory gives rise to a general-purpose cost-effective and scalable computational framework that can accurately capture the cavity-induced rate enhancement of charge transfer reactions that occurs when the molecular system is placed inside a microcavity. Both partial linearization with respect to the nuclear and photonic degrees of freedom and full linerization with respect to nuclear, photonic, and electronic degrees of freedom (the latter within the mapping Hamiltonian approach) are shown to be highly accurate, provided that the Wigner transforms of the product (WoP) of operators at the initial time is not replaced by the product of their Wigner transforms. We also show that the partial linearization method yields the quantum-mechanically exact cavity-modified FGR rate constant for a model system in which the donor and acceptor potential energy surfaces are harmonic and identical except for a shift in the equilibrium energy and geometry, if WoP is applied.

Resolving ambiguities of the mode truncation in cavity quantum electrodynamics

This work provides the fundamental theoretical framework for few-mode cavity quantum electrodynamics by resolving the gauge ambiguities between the Coulomb gauge and the dipole gauge Hamiltonians under the photonic mode truncation. We first propose a general framework to resolve ambiguities for an arbitrary truncation in a given gauge. Then, we specifically consider the case of mode truncation, deriving gauge invariant expressions for both the Coulomb and dipole gauge Hamiltonians that naturally reduce to the commonly used single-mode Hamiltonians when considering a single-mode truncation. We finally provide the analytical and numerical results of both atomic and molecular model systems coupled to the cavity to demonstrate the validity of our theory.

Molecular orbital theory in cavity QED environments

Coupling between molecules and vacuum photon fields inside an optical cavity has proven to be an effective way to engineer molecular properties, in particular reactivity. To ease the rationalization of cavity induced effects we introduce an ab initio method leading to the first fully consistent molecular orbital theory for quantum electrodynamics environments. Our framework is non-perturbative and explains modifications of the electronic structure due to the interaction with the photon field. In this work, we show that the newly developed orbital theory can be used to predict cavity induced modifications of molecular reactivity and pinpoint classes of systems with significant cavity effects. We also investigate electronic cavity-induced modifications of reaction mechanisms in vibrational strong coupling regimes.

Trajectory Ensemble Methods Provide Single-Molecule Statistics for Quantum Dynamical Systems

The emergence of experiments capable of probing quantum dynamics at the single-molecule level requires the development of new theoretical tools capable of simulating and analyzing these dynamics beyond an ensemble-averaged description. In this article, we present an efficient method for sampling and simulating the dynamics of the individual quantum systems that make up an ensemble and apply it to study the nonequilibrium dynamics of the ubiquitous spin-boson model. We generate an ensemble of single-system trajectories, and we analyze this trajectory ensemble using tools from classical statistical mechanics. Our results demonstrate that the dynamics of quantum coherence is highly heterogeneous at the single-system level due to variations in the initial bath configuration, which significantly affects the transient exchange of coherence between the system and its bath. We observe that single systems tend to retain coherence over time scales longer than that of the ensemble. We also compute a novel thermodynamic entanglement entropy that quantifies a thermodynamic driving force favoring system-bath entanglement.

Probing Light-Induced Conical Intersections by Monitoring Multidimensional Polaritonic Surfaces

The interaction of a molecule with the quantized electromagnetic field of a nanocavity gives rise to light-induced conical intersections between polaritonic potential energy surfaces. We demonstrate for a realistic model of a polyatomic molecule that the time-resolved ultrafast radiative emission of the cavity enables following both nuclear wavepacket dynamics on, and nonadiabatic population transfer between, polaritonic surfaces without applying a probe pulse. The latter provides an unambiguous (and in principle experimentally accessible) dynamical fingerprint of light-induced conical intersections.

Polariton ring currents and circular dichroism of Mg-porphyrin in a chiral cavity

By placing Mg-porphyrin molecules in a chiral optical cavity, time reversal symmetry is broken, and polariton ring currents can be generated with linearly polarized light, resulting in a circular dichroism signal. Since the electronic state degeneracy in the molecule is lifted by the formation of chiral polaritons, this signal is one order of magnitude stronger than the bare molecule signal induced by circularly polarized light. Enantiomer-selective photochemical processes in chiral optical cavities is an intriguing future possibility.

Placing aromatic molecules in a chiral optical cavity can break time-reversal symmetry and generate polariton ring currents with a linearly polarized pump. Such currents can be probed by circular dichroism, with one order of magnitude enhancement.

A partially linearized spin-mapping approach for simulating nonlinear optical spectra

We present a partially linearized method based on spin-mapping for computing both linear and nonlinear optical spectra. As observables are obtained from ensembles of classical trajectories, the approach can be applied to the large condensed-phase systems that undergo photosynthetic light-harvesting processes. In particular, the recently derived spin partially linearized density matrix method has been shown to exhibit superior accuracy in computing population dynamics compared to other related classical-trajectory methods. Such a method should also be ideally suited to describing the quantum coherences generated by interaction with light. We demonstrate that this is, indeed, the case by calculating the nonlinear optical response functions relevant for the pump-probe and 2D photon-echo spectra for a Frenkel biexciton model and the Fenna-Matthews-Olsen light-harvesting complex. One especially desirable feature of our approach is that the full spectrum can be decomposed into its constituent components associated with the various Liouville-space pathways, offering a greater insight beyond what can be directly obtained from experiments.

Molecular Polaritonics: Chemical Dynamics Under Strong Light-Matter Coupling

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. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Ultrastrong Exciton-Photon Coupling in Broadband Solar Absorbers

The recent development of organic polaritonic solar cells, in which sunlight absorbers and photon modes of a resonator are hybridized as a result of their strong coupling, has revealed the potential this interaction offers to control and enhance the performance of these devices. In this approach, the photovoltaic cell is built in such a way that it also behaves as an optical cavity supporting spectrally well-defined resonances, which match the broad absorption bands of the dyes employed. Herein we focus on the experimental and theoretical analysis of the specific spectral and angular optical absorption characteristics of a broadband light harvester, namely a subphthalocyanine, when operating in the ultrastrong coupling regime. We discuss the implications of having a broad distribution of oscillator strengths and demonstrate that rational design of the layered structure is needed to optimize both the spectral and the angular response of the sunlight harvester dye.

Selective isomer emission via funneling of exciton polaritons

Polaritons in organic systems has shown the potential to modify chemical properties and to mediate long-range energy transfer between individual chromophores, among other capabilities. Here, we demonstrate that strong coupling and formation of organic exciton-polaritons can be used to selectively tune the isomer emission of organic molecules. By taking advantage of their delocalized and hybrid character, polaritons emerging in the strong coupling regime open a new relaxation pathway that allows for an efficient funneling of the excitation between the molecular isomers. We implement this by strong coupling to trans-DCS (E-4-dimethylamino-4′cyanostilbene)molecules, which present two isomers in different amounts when immersed in a polymer matrix. Thanks to this new relaxation pathway, the photoexcitation that is first shared by the common polaritonic mode is then selectively funneled to the excited states of one of the isomers, recognizing pure emission from the isomeric states that do not contribute to emission under normal conditions.

Quantum-electrodynamical time-dependent density functional theory within Gaussian atomic basis

Inspired by the formulation of quantum-electrodynamical time-dependent density functional theory (QED-TDDFT) by Rubio and co-workers [Flick et al., ACS Photonics 6, 2757-2778 (2019)], we propose an implementation that uses dimensionless amplitudes for describing the photonic contributions to QED-TDDFT electron-photon eigenstates. This leads to a Hermitian QED-TDDFT coupling matrix that is expected to facilitate the future development of analytic derivatives. Through a Gaussian atomic basis implementation of the QED-TDDFT method, we examined the effect of dipole self-energy, rotating-wave approximation, and the Tamm-Dancoff approximation on the QED-TDDFT eigenstates of model compounds (ethene, formaldehyde, and benzaldehyde) in an optical cavity. We highlight, in the strong coupling regime, the role of higher-energy and off-resonance excited states with large transition dipole moments in the direction of the photonic field, which are automatically accounted for in our QED-TDDFT calculations and might substantially affect the energies and compositions of polaritons associated with lower-energy electronic states.

Polariton induced conical intersection and berry phase

We investigate the Polariton induced conical intersection (PICI) created from coupling a diatomic molecule with the quantized photon mode inside an optical cavity, and the corresponding Berry Phase effects. We use the rigorous Pauli-Fierz Hamiltonian to describe the quantum light-matter interactions between a LiF molecule and the cavity, and use the exact quantum propagation to investigate the polariton quantum dynamics. The molecular rotations relative to the cavity polarization direction play a role as the tuning mode of the PICI, resulting in an effective CI even within a diatomic molecule. To clearly demonstrate the dynamical effects of the Berry phase, we construct two additional models that have the same Born-Oppenheimer surface, but the effects of the geometric phase are removed. We find that when the initial wavefunction is placed in the lower polaritonic surface, the Berry phase causes a π phase-shift in the wavefunction after the encirclement around the CI, indicated from the nuclear probability distribution. On the other hand, when the initial wavefunction is placed in the upper polaritonic surface, the geometric phase significantly influences the couplings between polaritonic states and therefore, the population dynamics between them. These BP effects are further demonstrated through the photo-fragment angular distribution. PICI created from the quantized radiation field has the promise to open up new possibilities to modulate photochemical reactivities.

Changing Vibration Coupling Strengths of Liquid Acetonitrile with an Angle-Tuned Etalon

This work is the first report on nonzero molecular vibration-vibration coupling in an infrared cavity-vibration experiment. Vibration-vibration coupling strength is determined as a cavity mode of parallel spaced mirrors (etalon mode or fringe) is angle-tuned in the region between two vibrations of liquid acetonitrile which are Fermi coupled, namely, a CN stretch dominated vibration and a nearby combination band dominated by the symmetric CH₃ bend and C-C stretch. All other infrared cavity-vibration work to date involving more than one vibration has used a value of zero for vibration-vibration coupling; however, this work starts with Fermi coupled vibrations and reveals that there are changes in the vibration-vibration coupling and cavity-vibration couplings as the cavity mode is angle-tuned between the interacting vibrations. The ability to change fundamental vibrational dynamics within a cavity is an exciting result which helps to build a foundation for understanding molecular vibrational dynamics in parallel plate etalon cavities.

Molecular Polaritons Generated from Strong Coupling between CdSe Nanoplatelets and a Dielectric Optical Cavity

We demonstrate the formation of CdSe nanoplatelet (NPL) exciton-polaritons in a distributed Bragg reflector (DBR) cavity. The molecule-cavity hybrid system is in the strong coupling regime with an 83 meV Rabi splitting, characterized from angle-resolved reflectance and photoluminescence measurements. Mixed quantum-classical dynamics simulations are used to investigate the polariton photophysics of the hybrid system by treating the electronic and photonic degrees of freedom (DOF) quantum mechanically and the nuclear phononic DOF classically. Our numerical simulations of the angle-resolved photoluminescence (PL) agree extremely well with the experimental data, providing a fundamental explanation of the asymmetric intensity distribution of the upper and lower polariton branches. Our results also provide mechanistic insights into the importance of phonon-assisted nonadiabatic transitions among polariton states, which are reflected in the various features of the PL spectra. This work proves the feasibility of coupling nanoplatelet electronic states with the photon states of a dielectric cavity to form a hybrid system and provides a new platform for investigating cavity-mediated physical and chemical processes.

Benchmarking Quasiclassical Mapping Hamiltonian Methods for Simulating Cavity-Modified Molecular Dynamics

Recent experimental realizations of strong coupling between optical cavity modes and molecular matter placed inside the cavity have opened exciting new routes for controlling chemical processes. Simulating the cavity-modified dynamics of complex chemical systems calls for the development of accurate, flexible, and cost-effective approximate numerical methods that scale favorably with system size and complexity. In this Letter, we test the ability of quasiclassical mapping Hamiltonian methods to serve this purpose. We simulated the spontaneous emission dynamics of an atom confined to a microcavity via five different variations of the linearized semiclassical (LSC) method. Our main finding is that recently proposed LSC-based methods which use a modified form of the identity operator are reasonably accurate and perform significantly better than the Ehrenfest and standard LSC methods, without significantly increasing computational costs. These methods are therefore highly promising as a general purpose tool for simulating cavity-modified dynamics of complex chemical systems.