Multiscale Molecular Dynamics Simulations of Polaritonic Chemistry

A Vision for the Future of Multiscale Modeling

The rise of modern computer science enabled physical chemistry to make enormous progresses in understanding and harnessing natural and artificial phenomena. Nevertheless, despite the advances achieved over past decades, computational resources are still insufficient to thoroughly simulate extended systems from first principles. Indeed, countless biological, catalytic and photophysical processes require ab initio treatments to be properly described, but the breadth of length and time scales involved makes it practically unfeasible. A way to address these issues is to couple theories and algorithms working at different scales by dividing the system into domains treated at different levels of approximation, ranging from quantum mechanics to classical molecular dynamics, even including continuum electrodynamics. This approach is known as multiscale modeling and its use over the past 60 years has led to remarkable results. Considering the rapid advances in theory, algorithm design, and computing power, we believe multiscale modeling will massively grow into a dominant research methodology in the forthcoming years. Hereby we describe the main approaches developed within its realm, highlighting their achievements and current drawbacks, eventually proposing a plausible direction for future developments considering also the emergence of new computational techniques such as machine learning and quantum computing. We then discuss how advanced multiscale modeling methods could be exploited to address critical scientific challenges, focusing on the simulation of complex light-harvesting processes, such as natural photosynthesis. While doing so, we suggest a cutting-edge computational paradigm consisting in performing simultaneous multiscale calculations on a system allowing the various domains, treated with appropriate accuracy, to move and extend while they properly interact with each other. Although this vision is very ambitious, we believe the quick development of computer science will lead to both massive improvements and widespread use of these techniques, resulting in enormous progresses in physical chemistry and, eventually, in our society.

Unraveling a Cavity-Induced Molecular Polarization Mechanism from Collective Vibrational Strong Coupling

We demonstrate that collective vibrational strong coupling of molecules in thermal equilibrium can give rise to significant local electronic polarizations in the thermodynamic limit. We do so by first showing that the full nonrelativistic Pauli-Fierz problem of an ensemble of strongly coupled molecules in the dilute-gas limit reduces in the cavity Born-Oppenheimer approximation to a cavity-Hartree equation for the electronic structure. Consequently, each individual molecule experiences a self-consistent coupling to the dipoles of all other molecules, which amount to non-negligible values in the thermodynamic limit (large ensembles). Thus, collective vibrational strong coupling can alter individual molecules strongly for localized "hotspots" within the ensemble. Moreover, the discovered cavity-induced polarization pattern possesses a zero net polarization, which resembles a continuous form of a spin glass (or better polarization glass). Our findings suggest that the thorough understanding of polaritonic chemistry, requires a self-consistent treatment of dressed electronic structure, which can give rise to numerous, so far overlooked, physical mechanisms.

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.

Non-Hermitian molecular dynamics simulations of exciton-polaritons in lossy cavities

The observation that materials can change their properties when placed inside or near an optical resonator has sparked a fervid interest in understanding the effects of strong light-matter coupling on molecular dynamics, and several approaches have been proposed to extend the methods of computational chemistry into this regime. Whereas the majority of these approaches have focused on modeling a single molecule coupled to a single cavity mode, changes to chemistry have so far only been observed experimentally when very many molecules are coupled collectively to multiple modes with short lifetimes. While atomistic simulations of many molecules coupled to multiple cavity modes have been performed with semi-classical molecular dynamics, an explicit description of cavity losses has so far been restricted to simulations in which only a very few molecular degrees of freedom were considered. Here, we have implemented an effective non-Hermitian Hamiltonian to explicitly treat cavity losses in large-scale semi-classical molecular dynamics simulations of organic polaritons and used it to perform both mean-field and surface hopping simulations of polariton relaxation, propagation, and energy transfer.

Collective Strong Coupling Modifies Aggregation and Solvation

Intermolecular (Coulombic) interactions are pivotal for aggregation, solvation, and crystallization. We demonstrate that the collective strong coupling of several molecules to a single optical mode results in notable changes in the molecular excitations around a single perturbed molecule, thus representing an impurity in an otherwise ordered system. A competition between short-range coulombic and long-range photonic correlations inverts the local transition density in a polaritonic state, suggesting notable changes in the polarizability of the solvation shell. Our results provide an alternative perspective on recent work in polaritonic chemistry and pave the way for the rigorous treatment of cooperative effects in aggregation, solvation, and crystallization.

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.

Catalytic Reaction Mechanism of Bacterial GH92 α-1,2-Mannosidase: A QM/MM Metadynamics Study

The catalytic mechanism of aC a + 2 ${C{a}^{+2}}$ -dependent family 92 α ${{\rm \alpha }}$ -mannosidase, which is abundantly present in human gut flora and malfunctions leading to the lysosomal storage disease α-mannosidosis, has been investigated using quantum mechanics/molecular mechanics and metadynamics methods. Computational efforts show that the enzyme follows a conformational itinerary of and theC a + 2 ${C{a}^{+2}}$ ion serves a dual purpose, as it not only distorts the sugar ring but also plays a crucial role in orchestrating the arrangement of catalytic residues. This orchestration, in turn, contributes to the facilitation of O S 2 ${{{\rm \ }}^{{\rm O}}{{\rm S}}_{2}}$ conformers for the ensuing reaction. This mechanistic insight is well-aligned with the experimental predictions of the catalytic pathway, and the computed energies are of the same order of magnitude as the experimental estimations. Hence, our results extend the mechanistic understanding of glycosidases.

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.

Tuning the Coherent Propagation of Organic Exciton-Polaritons through the Cavity Q-factor

Transport of excitons in organic materials can be enhanced through polariton formation when the interaction strength between these excitons and the confined light modes of an optical resonator exceeds their decay rates. While the polariton lifetime is determined by the Q(uality)-factor of the optical resonator, the polariton group velocity is not. Instead, the latter is solely determined by the polariton dispersion. Yet, experiments suggest that the Q-factor also controls the polariton propagation velocity. To understand this observation, the authors perform molecular dynamics simulations of Rhodamine chromophores strongly coupled to Fabry-Pérot cavities with various Q-factors. The results suggest that propagation in the aforementioned experiments is initially dominated by ballistic motion of upper polariton states at their group velocities, which leads to a rapid expansion of the wavepacket. Cavity decay in combination with non-adiabatic population transfer into dark states, rapidly depletes these bright states, causing the wavepacket to contract. However, because population transfer is reversible, propagation continues, but as a diffusion process, at lower velocity. By controlling the lifetime of bright states, the Q-factor determines the duration of the ballistic phase and the diffusion coefficient in the diffusive regime. Thus, polariton propagation in organic microcavities can be effectively tuned through the Q-factor.

Multi-scale molecular dynamics simulations of enhanced energy transfer in organic molecules under strong coupling

Exciton transport can be enhanced in the strong coupling regime where excitons hybridize with confined light modes to form polaritons. Because polaritons have group velocity, their propagation should be ballistic and long-ranged. However, experiments indicate that organic polaritons propagate in a diffusive manner and more slowly than their group velocity. Here, we resolve this controversy by means of molecular dynamics simulations of Rhodamine molecules in a Fabry-Pérot cavity. Our results suggest that polariton propagation is limited by the cavity lifetime and appears diffusive due to reversible population transfers between polaritonic states that propagate ballistically at their group velocity, and dark states that are stationary. Furthermore, because long-lived dark states transiently trap the excitation, propagation is observed on timescales beyond the intrinsic polariton lifetime. These insights not only help to better understand and interpret experimental observations, but also pave the way towards rational design of molecule-cavity systems for coherent exciton transport.

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.

Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics

The exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron-ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron-nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

The Rise and Current Status of Polaritonic Photochemistry and Photophysics

The interaction between molecular electronic transitions and electromagnetic fields can be enlarged to the point where distinct hybrid light-matter states, polaritons, emerge. The photonic contribution to these states results in increased complexity as well as an opening to modify the photophysics and photochemistry beyond what normally can be seen in organic molecules. It is today evident that polaritons offer opportunities for molecular photochemistry and photophysics, which has caused an ever-rising interest in the field. Focusing on the experimental landmarks, this review takes its reader from the advent of the field of polaritonic chemistry, over the split into polariton chemistry and photochemistry, to present day status within polaritonic photochemistry and photophysics. To introduce the field, the review starts with a general description of light-matter interactions, how to enhance these, and what characterizes the coupling strength. Then the photochemistry and photophysics of strongly coupled systems using Fabry-Perot and plasmonic cavities are described. This is followed by a description of room-temperature Bose-Einstein condensation/polariton lasing in polaritonic systems. The review ends with a discussion on the benefits, limitations, and future developments of strong exciton-photon coupling using organic molecules.

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.

Strong light-matter coupling in pentacene thin films on plasmonic arrays

Utilizing strong light-matter coupling is an elegant and powerful way to modify the energy landscapes of excited states of organic semiconductors. Consequently, the chemical and photophysical properties of these organic semiconductors can be influenced without the need of chemical modification but simply by implementing them in optical microcavities. This has so far mostly been shown in Fabry-Pérot cavities and with organic single crystals or diluted molecules in a host matrix. Here, we demonstrate strong, simultaneous coupling of the two Davydov transitions in polycrystalline pentacene thin films to surface lattice resonances supported by open cavities made of silver nanoparticle arrays. Such thin films are more easily fabricated and, together with the open architecture, more suitable for device applications.

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.

Vibrational coupling with O-H stretching increases catalytic efficiency of sucrase in Fabry-Pérot microcavity

Vibrational strong coupling (VSC) has been reported as a polariton-based method for modulating the rate of biochemical reactions. Herein, we studied how VSC modulates the sucrose hydrolysis. By monitoring the refractive index-induced shift of Fabry-Pérot microcavity, in which the catalytic efficiency of sucrose hydrolysis can be increased at least two times, as VSC was tuned to resonate with the stretching vibration of O-H bonds. This research provides new evidence for applying VSC in life sciences, which holds great promise to improving enzymatic industries.

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.

QM/MM Modeling of Vibrational Polariton Induced Energy Transfer and Chemical Dynamics

Vibrational strong coupling (VSC) provides a novel means to modify chemical reactions and energy transfer pathways. To efficiently model chemical dynamics under VSC in the collective regime, herein a hybrid quantum mechanical/molecular mechanical (QM/MM) cavity molecular dynamics (CavMD) scheme is developed and applied to an experimentally studied chemical system. This approach can achieve linear scaling with respect to the number of molecules for a dilute solution under VSC by assuming that each QM solute molecule is surrounded by an independent MM solvent bath. Application of this approach to a dilute solution of Fe(CO)₅ in n-dodecane under VSC demonstrates polariton dephasing to the dark modes and polariton-enhanced molecular nonlinear absorption. These simulations predict that strongly exciting the lower polariton may provide an energy transfer pathway that selectively excites the equatorial CO vibrations rather than the axial CO vibrations. Moreover, these simulations also directly probe the cavity effect on the dynamics of the Fe(CO)₅ Berry pseudorotation reaction for comparison to recent two-dimensional infrared spectroscopy experiments. This theoretical approach is applicable to a wide range of other polaritonic systems and provides a tool for exploring the use of VSC for selective infrared photochemistry.

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.

Polaritonic Chemistry from First Principles via Embedding Radiation Reaction

The coherent interaction of a large collection of molecules with a common photonic mode results in strong light-matter coupling, a feature that has proven highly beneficial for chemistry and has introduced the research topics polaritonic and QED chemistry. Here, we demonstrate an embedding approach to capture the collective nature while retaining the full ab initio representation of single molecules─an approach ideal for polaritonic chemistry. The accuracy of the embedding radiation-reaction ansatz is demonstrated for time-dependent density-functional theory. Then, by virtue of a simple proton-tunneling model, we illustrate that the influence of collective strong coupling on chemical reactions features a nontrivial dependence on the number of emitters and can alternate between strong catalyzing and an inhibiting effect. Bridging classical electrodynamics, quantum optical descriptions, and the ab initio description of realistic molecules, this work can serve as a guiding light for future developments and investigations in the quickly growing fields of QED chemistry and QED material design.

Controlling Exciton Propagation in Organic Crystals through Strong Coupling to Plasmonic Nanoparticle Arrays

Exciton transport in most organic materials is based on an incoherent hopping process between neighboring molecules. This process is very slow, setting a limit to the performance of organic optoelectronic devices. In this Article, we overcome the incoherent exciton transport by strongly coupling localized singlet excitations in a tetracene crystal to confined light modes in an array of plasmonic nanoparticles. We image the transport of the resulting exciton-polaritons in Fourier space at various distances from the excitation to directly probe their propagation length as a function of the exciton to photon fraction. Exciton-polaritons with an exciton fraction of 50% show a propagation length of 4.4 μm, which is an increase by 2 orders of magnitude compared to the singlet exciton diffusion length. This remarkable increase has been qualitatively confirmed with both finite-difference time-domain simulations and atomistic multiscale molecular dynamics simulations. Furthermore, we observe that the propagation length is modified when the dipole moment of the exciton transition is either parallel or perpendicular to the cavity field, which opens a new avenue for controlling the anisotropy of the exciton flow in organic crystals. The enhanced exciton-polariton transport reported here may contribute to the development of organic devices with lower recombination losses and improved performance.

Energy-efficient pathway for selectively exciting solute molecules to high vibrational states via solvent vibration-polariton pumping

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 ¹³CO₂ solute in a ¹²CO₂ 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.

Hybrid light-matter states formed in the strong light-matter coupling regime can alter the molecular ground-state reactivity. Here, Li et al. computationally demonstrate that pumping a collection of solvent molecules forming hybrid vibrational light-matter states in an optical cavity can excite solute molecules to very high excited states.

Identifying Vibrations that Control Non-adiabatic Relaxation of Polaritons in Strongly Coupled Molecule-Cavity Systems

The strong light-matter coupling regime, in which excitations of materials hybridize with excitations of confined light modes into polaritons, holds great promise in various areas of science and technology. A key aspect for all applications of polaritonic chemistry is the relaxation into the lower polaritonic states. Polariton relaxation is speculated to involve two separate processes: vibrationally assisted scattering (VAS) and radiative pumping (RP), but the driving forces underlying these two mechanisms are not fully understood. To provide mechanistic insights, we performed multiscale molecular dynamics simulations of tetracene molecules strongly coupled to the confined light modes of an optical cavity. The results suggest that both mechanisms are driven by the same molecular vibrations that induce relaxation through nonadiabatic coupling between dark states and polaritonic states. Identifying these vibrational modes provides a rationale for enhanced relaxation into the lower polariton when the cavity detuning is resonant with specific vibrational transitions.

A perspective on ab initio modeling of polaritonic chemistry: The role of non-equilibrium effects and quantum collectivity

This Perspective provides a brief introduction into the theoretical complexity of polaritonic chemistry, which emerges from the hybrid nature of strongly coupled light-matter states. To tackle this complexity, the importance of ab initio methods is highlighted. Based on those, novel ideas and research avenues are developed with respect to quantum collectivity, as well as for resonance phenomena immanent in reaction rates under vibrational strong coupling. Indeed, fundamental theoretical questions arise about the mesoscopic scale of quantum-collectively coupled molecules when considering the depolarization shift in the interpretation of experimental data. Furthermore, to rationalize recent findings based on quantum electrodynamical density-functional theory (QEDFT), a simple, but computationally efficient, Langevin framework is proposed based on well-established methods from molecular dynamics. It suggests the emergence of cavity-induced non-equilibrium nuclear dynamics, where thermal (stochastic) resonance phenomena could emerge in the absence of external periodic driving. Overall, we believe that the latest ab initio results indeed suggest a paradigmatic shift for ground-state chemical reactions under vibrational strong coupling from the collective quantum interpretation toward a more local, (semi)-classically and non-equilibrium dominated perspective. Finally, various extensions toward a refined description of cavity-modified chemistry are introduced in the context of QEDFT, and future directions of the field are sketched.

Cavity-Modified Unimolecular Dissociation Reactions via Intramolecular Vibrational Energy Redistribution

While the emerging field of vibrational polariton chemistry has the potential to overcome traditional limitations of synthetic chemistry, the underlying mechanism is not yet well understood. Here, we explore how the dynamics of unimolecular dissociation reactions that are rate-limited by intramolecular vibrational energy redistribution (IVR) can be modified inside an infrared optical cavity. We study a classical model of a bent triatomic molecule, where the two outer atoms are bound by anharmonic Morse potentials to the center atom coupled to a harmonic bending mode. We show that an optical cavity resonantly coupled to particular anharmonic vibrational modes can interfere with IVR and alter unimolecular dissociation reaction rates when the cavity mode acts as a reservoir for vibrational energy. These results lay the foundation for further theoretical work toward understanding the intriguing experimental results of vibrational polaritonic chemistry within the context of IVR.

Theoretical Challenges in Polaritonic Chemistry

Polaritonic chemistry exploits strong light-matter coupling between molecules and confined electromagnetic field modes to enable new chemical reactivities. In systems displaying this functionality, the choice of the cavity determines both the confinement of the electromagnetic field and the number of molecules that are involved in the process. While in wavelength-scale optical cavities the light-matter interaction is ruled by collective effects, plasmonic subwavelength nanocavities allow even single molecules to reach strong coupling. Due to these very distinct situations, a multiscale theoretical toolbox is then required to explore the rich phenomenology of polaritonic chemistry. Within this framework, each component of the system (molecules and electromagnetic modes) needs to be treated in sufficient detail to obtain reliable results. Starting from the very general aspects of light-molecule interactions in typical experimental setups, we underline the basic concepts that should be taken into account when operating in this new area of research. Building on these considerations, we then provide a map of the theoretical tools already available to tackle chemical applications of molecular polaritons at different scales. Throughout the discussion, we draw attention to both the successes and the challenges still ahead in the theoretical description of polaritonic chemistry.

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

Molecular polaritons have become an emerging platform for remotely controlling molecular properties through strong light-matter interactions. Herein, a semiclassical approach is developed for describing molecular polaritons by self-consistently propagating the real-time dynamics of classical cavity modes and a quantum molecular subsystem described by the nuclear-electronic orbital (NEO) method, where electrons and specified nuclei are treated quantum mechanically on the same level. This semiclassical real-time NEO approach provides a unified description of electronic and vibrational strong couplings and describes the impact of the cavity on coupled nuclear-electronic dynamics while including nuclear quantum effects. For a single o-hydroxybenzaldehyde molecule under electronic strong coupling, this approach shows that the cavity suppression of excited state intramolecular proton transfer is influenced not only by the polaritonic potential energy surface but also by the time scale of the chemical reaction. This work provides the foundation for exploring collective strong coupling in nuclear-electronic quantum dynamical systems within optical cavities.

Polariton relaxation under vibrational strong coupling: Comparing cavity molecular dynamics simulations against Fermi's golden rule rate

Under vibrational strong coupling (VSC), the formation of molecular polaritons may significantly modify the photo-induced or thermal properties of molecules. In an effort to understand these intriguing modifications, both experimental and theoretical studies have focused on the ultrafast dynamics of vibrational polaritons. Here, following our recent work [Li et al., J. Chem. Phys. 154, 094124 (2021)], we systematically study the mechanism of polariton relaxation for liquid CO₂ under a weak external pumping. Classical cavity molecular dynamics (CavMD) simulations confirm that polariton relaxation results from the combined effects of (i) cavity loss through the photonic component and (ii) dephasing of the bright-mode component to vibrational dark modes as mediated by intermolecular interactions. The latter polaritonic dephasing rate is proportional to the product of the weight of the bright mode in the polariton wave function and the spectral overlap between the polariton and dark modes. Both these factors are sensitive to parameters such as the Rabi splitting and cavity mode detuning. Compared to a Fermi's golden rule calculation based on a tight-binding harmonic model, CavMD yields a similar parameter dependence for the upper polariton relaxation lifetime but sometimes a modest disagreement for the lower polariton. We suggest that this disagreement results from polariton-enhanced molecular nonlinear absorption due to molecular anharmonicity, which is not included in our analytical model. We also summarize recent progress on probing nonreactive VSC dynamics with CavMD.

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.

Experimental characterization techniques for plasmon-assisted chemistry

Plasmon-assisted chemistry is the result of a complex interplay between electromagnetic near fields, heat and charge transfer on the nanoscale. The disentanglement of their roles is non-trivial. Therefore, a thorough knowledge of the chemical, structural and spectral properties of the plasmonic/molecular system being used is required. Specific techniques are needed to fully characterize optical near fields, temperature and hot carriers with spatial, energetic and/or temporal resolution. The timescales for all relevant physical and chemical processes can range from a few femtoseconds to milliseconds, which necessitates the use of time-resolved techniques for monitoring the underlying dynamics. In this Review, we focus on experimental techniques to tackle these challenges. We further outline the difficulties when going from the ensemble level to single-particle measurements. Finally, a thorough understanding of plasmon-assisted chemistry also requires a substantial joint experimental and theoretical effort.

A Theoretical Investigation about Photoswitching of Azobenzene Adsorbed on Ag Nanoparticles

The optical properties of hybrid systems composed of silver nanoparticles (NPs) and azobenzene molecules were systematically investigated by combining the real-time time-dependent density functional theory (RT-TDDFT) approach with the classical electrodynamics finite difference time domain (FDTD) technique for the solution of Maxwell’s equations. In order to reflect the chemical interaction between azobenzene and metal more exactly, except for adsorbed molecules, a Ag cluster separated from NP was also dealt, using RT-TDDFT. We studied the different factors affecting the surface-enhanced absorption spectra. It was found that the electric field amplified by plasmon resonance of Ag NPs can have an overall enhancement to the molecular light absorption throughout the whole energy range. The resonance between the electron and the plasmon excitation results in a larger percentage of enhancement in the absorption spectrum the closer the resonance peak is. The enhancement ratio of the resonance peak is the largest. The plasmon–exciton coupling and the optical properties of different isolate isomers influence the line shape of the absorption spectra. The dipole interaction and electronic transfer between azobenzene molecules and Ag NPs also change the shape of spectroscopy from the absorption enhancement ratio and the location of the peak. Physical and chemical factors lead to photoswitching in these hybrid systems together.

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.

Metamaterial Analogues of Strongly Coupled Molecular Ensembles

The formation of polariton modes due to the strong coupling of light and matter has led to exciting developments in physics, chemistry, and materials science. The potential to modify the properties of molecular materials by strongly coupling molecules to a confined light field is so far-reaching and so attractive that a new field known as "polaritonic chemistry" is now emerging. However, the molecular scale of the materials involved makes probing strong coupling at the individual resonator level extremely challenging. Here, we offer a complementary approach based upon metamaterials, an approach that enables us to use cm-scale structures, thereby opening a new way to explore strong coupling phenomena. As proof-of-principle, we show that metamolecules placed inside a radio frequency cavity may exhibit strong coupling and show that near-field radio frequency techniques allow us, for the first time, to probe the response of individual metamolecules under strong coupling conditions.

Nonlinear features of Fano resonance: a QM/EM study

The nonlinear Fano effects on the absorption of hybrid systems composed of a silver nanosphere and an indoline dye molecule have been systematically investigated by the hybrid approach, which combines the quantum mechanics method (QM) with the computational electromagnetic method (EM). The absorption spectra of the dye molecule in the proximity of an Ag nanoparticle have been calculated by changing the incident field intensity, the phenomenological dephasing of molecular excitation, and the enhancement ratio of the near field. The contribution of molecular nonlinear response properties and the quantum interferences of the incident and scattered fields and of resonant plasmon-molecular excitations to the spectra has been identified. It is in no doubt that Fano resonance due to the plasmon-molecular interaction can appear in both the weak and strong field regimes; however, the Fano effect is more pronounced in the strong field regime where quantum interference leads to a nonlinear Fano effect controlled by a complex field-dependent Fano factor. When the incident field is strong enough, the resonance antisymmetry structure is spectrally resolved, and it changes with the change of the field intensity. As the field intensity varies from weak to strong, the Fano lineshape's asymmetry increases with increasing intensity in the beginning, and then decreases with a further increase of the field intensity attributed to the increase of the detuning energy induced by the integrated energy shift upon field dressing during the excitation. Decreasing the enhancement ratio of the near field or the dephasing of molecular excitation can also control the spectral lineshape transformation from an asymmetric profile to a symmetric Lorentzian lineshape. These findings are consistent with previous experimental and theoretical observations arisen by quantum interferences and are expected to stimulate further work toward exploring the plasmon-molecular interplay and the applications of Fano resonance in optical switching and sensing.

Effect of molecular Stokes shift on polariton dynamics

When the enhanced electromagnetic field of a confined light mode interacts with photoactive molecules, the system can be driven into the regime of strong coupling, where new hybrid light-matter states, polaritons, are formed. Polaritons, manifested by the Rabi split in the dispersion, have shown potential for controlling the chemistry of the coupled molecules. Here, we show by angle-resolved steady-state experiments accompanied by multi-scale molecular dynamics simulations that the molecular Stokes shift plays a significant role in the relaxation of polaritons formed by organic molecules embedded in a polymer matrix within metallic Fabry-Pérot cavities. Our results suggest that in the case of Rhodamine 6G, a dye with a significant Stokes shift, excitation of the upper polariton leads to a rapid localization of the energy into the fluorescing state of one of the molecules, from where the energy scatters into the lower polariton (radiative pumping), which then emits. In contrast, for excitonic J-aggregates with a negligible Stokes shift, the fluorescing state does not provide an efficient relaxation gateway. Instead, the relaxation is mediated by exchanging energy quanta matching the energy gap between the dark states and lower polariton into vibrational modes (vibrationally assisted scattering). To understand better how the fluorescing state of a molecule that is not strongly coupled to the cavity can transfer its excitation energy to the lower polariton in the radiative pumping mechanism, we performed multi-scale molecular dynamics simulations. The results of these simulations suggest that non-adiabatic couplings between uncoupled molecules and the polaritons are the driving force for this energy transfer process.

Multi-scale dynamics simulations of molecular polaritons: The effect of multiple cavity modes on polariton relaxation

Coupling molecules to the confined light modes of an optical cavity is showing great promise for manipulating chemical reactions. However, to fully exploit this principle and use cavities as a new tool for controlling chemistry, a complete understanding of the effects of strong light-matter coupling on molecular dynamics and reactivity is required. While quantum chemistry can provide atomistic insight into the reactivity of uncoupled molecules, the possibilities to also explore strongly coupled systems are still rather limited due to the challenges associated with an accurate description of the cavity in such calculations. Despite recent progress in introducing strong coupling effects into quantum chemistry calculations, applications are mostly restricted to single or simplified molecules in ideal lossless cavities that support a single light mode only. However, even if commonly used planar mirror micro-cavities are characterized by a fundamental mode with a frequency determined by the distance between the mirrors, the cavity energy also depends on the wave vector of the incident light rays. To account for this dependency, called cavity dispersion, in atomistic simulations of molecules in optical cavities, we have extended our multi-scale molecular dynamics model for strongly coupled molecular ensembles to include multiple confined light modes. To validate the new model, we have performed simulations of up to 512 Rhodamine molecules in red-detuned Fabry-Pérot cavities. The results of our simulations suggest that after resonant excitation into the upper polariton at a fixed wave vector, or incidence angle, the coupled cavity-molecule system rapidly decays into dark states that lack dispersion. Slower relaxation from the dark state manifold into both the upper and lower bright polaritons causes observable photo-luminescence from the molecule-cavity system along the two polariton dispersion branches that ultimately evolves toward the bottom of the lower polariton branch, in line with experimental observations. We anticipate that the more realistic cavity description in our approach will help to better understand and predict how cavities can modify molecular properties.

Dipolar coupling of nanoparticle-molecule assemblies: An efficient approach for studying strong coupling

Strong light-matter interactions facilitate not only emerging applications in quantum and non-linear optics but also modifications of properties of materials. In particular, the latter possibility has spurred the development of advanced theoretical techniques that can accurately capture both quantum optical and quantum chemical degrees of freedom. These methods are, however, computationally very demanding, which limits their application range. Here, we demonstrate that the optical spectra of nanoparticle-molecule assemblies, including strong coupling effects, can be predicted with good accuracy using a subsystem approach, in which the response functions of different units are coupled only at the dipolar level. We demonstrate this approach by comparison with previous time-dependent density functional theory calculations for fully coupled systems of Al nanoparticles and benzene molecules. While the present study only considers few-particle systems, the approach can be readily extended to much larger systems and to include explicit optical-cavity modes.

Cavity molecular dynamics simulations of vibrational polariton-enhanced molecular nonlinear absorption

Recent experiments have observed that the chemical and photophysical properties of molecules can be modified inside an optical Fabry-Pérot 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 simulations to numerically investigate the linear and the 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 that 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.

Nonadiabatic phenomena in molecular vibrational polaritons

Nonadiabatic phenomena are investigated in the rovibrational motion of molecules confined in an infrared cavity. Conical intersections (CIs) between vibrational polaritons, similar to CIs between electronic polaritonic surfaces, are found. The spectral, topological, and dynamic properties of the vibrational polaritons show clear fingerprints of nonadiabatic couplings between molecular vibration, rotation, and the cavity photonic mode. Furthermore, it is found that for the investigated system, composed of two rovibrating HCl molecules and the cavity mode, breaking the molecular permutational symmetry, by changing   ³⁵Cl to   ³⁷Cl in one of the HCl molecules, the polaritonic surfaces, nonadiabatic couplings, and related spectral, topological, and dynamic properties can deviate substantially. This implies that the natural occurrence of different molecular isotopologues needs to be considered when modeling realistic polaritonic systems.

Polaritonic Chemistry: Collective Strong Coupling Implies Strong Local Modification of Chemical Properties

A fundamental question in the field of polaritonic chemistry is whether collective coupling implies local modifications of chemical properties scaling with the ensemble size. Here we demonstrate from first-principles that an impurity present in a collectively coupled chemical ensemble features such locally scaling modifications. In particular, we find the formation of a novel dark state for a nitrogen dimer chain of variable size, whose local chemical properties are altered considerably at the impurity due to its embedding in the collectively coupled environment. Our simulations unify theoretical predictions from quantum optical models (e.g., collective dark states and bright polaritonic branches) with the single molecule quantum chemical perspective, which relies on the (quantized) redistribution of charges leading to a local hybridization of light and matter. Moreover, our findings suggest that recently developed ab initio methods for strong light-matter coupling are suitable to access these local polaritonic effects and provide a detailed understanding of photon-modified chemistry.

Simulating photodissociation reactions in bad cavities with the Lindblad equation

Optical cavities, e.g., as used in organic polariton experiments, often employ low finesse mirrors or plasmonic structures. The photon lifetime in these setups is comparable to the timescale of the nuclear dynamics governing the photochemistry. This highlights the need for including the effect of dissipation in the molecular simulations. In this study, we perform wave packet dynamics with the Lindblad master equation to study the effect of a finite photon lifetime on the dissociation of the MgH⁺ molecule model system. Photon lifetimes of several different orders of magnitude are considered to encompass an ample range of effects inherent to lossy cavities.

Hybrid theoretical models for molecular nanoplasmonics

The multidisciplinary nature of the research in molecular nanoplasmonics, i.e., the use of plasmonic nanostructures to enhance, control, or suppress properties of molecules interacting with light, led to contributions from different theory communities over the years, with the aim of understanding, interpreting, and predicting the physical and chemical phenomena occurring at molecular- and nano-scale in the presence of light. Multiscale hybrid techniques, using a different level of description for the molecule and the plasmonic nanosystems, permit a reliable representation of the atomistic details and of collective features, such as plasmons, in such complex systems. Here, we focus on a selected set of topics of current interest in molecular plasmonics (control of electronic excitations in light-harvesting systems, polaritonic chemistry, hot-carrier generation, and plasmon-enhanced catalysis). We discuss how their description may benefit from a hybrid modeling approach and what are the main challenges for the application of such models. In doing so, we also provide an introduction to such models and to the selected topics, as well as general discussions on their theoretical descriptions.

Born-Oppenheimer approximation in optical cavities: from success to breakdown

The coupling of a molecule and a cavity induces nonadiabaticity in the molecule which makes the description of its dynamics complicated. For polyatomic molecules, reduced-dimensional models and the use of the Born-Oppenheimer approximation (BOA) may remedy the situation. It is demonstrated that contrary to expectation, BOA may even fail in a one-dimensional model and is generally expected to fail in two- or more-dimensional models due to the appearance of conical intersections induced by the cavity.

Virial Relations for Electrons Coupled to Quantum Field Modes

In this work, we present a set of virial relations for many electron systems coupled to both classical and quantum fields, described by the Pauli-Fierz Hamiltonian in dipole approximation and using length gauge. Currently, there is growing interest in solutions of this Hamiltonian because of its relevance for describing molecular systems strongly coupled to photonic modes in cavities and in the possible modification of chemical properties of such systems compared to the ones in free space. The relevance of such virial relations is demonstrated by showing a connection to mass renormalization and by providing an exact way to obtain total energies from potentials in the framework of quantum electrodynamical density functional theory.

Chemistry in Quantum Cavities: Exact Results, the Impact of Thermal Velocities, and Modified Dissociation

In recent years tremendous progress in the field of light-matter interactions has unveiled that strong coupling to the modes of an optical cavity can alter chemistry even at room temperature. Despite these impressive advances, many fundamental questions of chemistry in cavities remain unanswered. This is also due to a lack of exact results that can be used to validate and benchmark approximate approaches. In this work we provide such reference calculations from exact diagonalization of the Pauli-Fierz Hamiltonian in the long-wavelength limit with an effective cavity mode. This allows us to investigate the reliability of the ubiquitous Jaynes-Cummings model not only for electronic but also for the case of ro-vibrational transitions. We demonstrate how the commonly ignored thermal velocity of charged molecular systems can influence chemical properties while leaving the spectra invariant. Furthermore, we show the emergence of new bound polaritonic states beyond the dissociation energy limit.