Generalization of the Tavis-Cummings model for multi-level anharmonic systems: Insights on the second excitation manifold

Polariton spectra under the collective coupling regime. II. 2D non-linear spectra

In our previous work [Mondal et al., J. Chem. Phys. 162, 014114 (2025)], we developed several efficient computational approaches to simulate exciton-polariton dynamics described by the Holstein-Tavis-Cummings (HTC) Hamiltonian under the collective coupling regime. Here, we incorporated these strategies into the previously developed Lindblad-partially linearized density matrix (L-PLDM) approach for simulating 2D electronic spectroscopy (2DES) of exciton-polariton under the collective coupling regime. In particular, we apply the efficient quantum dynamics propagation scheme developed in Paper I to both the forward and the backward propagations in the PLDM and develop an efficient importance sampling scheme and graphics processing unit vectorization scheme that allow us to reduce the computational costs from O(K2)O(T3) to O(K)O(T0) for the 2DES simulation, where K is the number of states and T is the number of time steps of propagation. We further simulated the 2DES for an HTC Hamiltonian under the collective coupling regime and analyzed the signal from both rephasing and non-rephasing contributions of the ground state bleaching, excited state emission, and stimulated emission pathways.

Polariton spectra under the collective coupling regime. I. Efficient simulation of linear spectra and quantum dynamics

We outline two general theoretical techniques to simulate polariton quantum dynamics and optical spectra under the collective coupling regimes described by a Holstein-Tavis-Cummings (HTC) model Hamiltonian. The first one takes advantage of sparsity of the HTC Hamiltonian, which allows one to reduce the cost of acting polariton Hamiltonian onto a state vector to the linear order of the number of states, instead of the quadratic order. The second one is applying the well-known Chebyshev series expansion approach for quantum dynamics propagation and to simulate the polariton dynamics in the HTC system; this approach allows us to use a much larger time step for propagation and only requires a few recursive operations of the polariton Hamiltonian acting on state vectors. These two theoretical approaches are general and can be applied to any trajectory-based non-adiabatic quantum dynamics methods. We apply these two techniques with our previously developed Lindblad-partially linearized density matrix approach to simulate the linear absorption spectra of the HTC model system, with both inhomogeneous site energy disorders and dipolar orientational disorders. Our numerical results agree well with the previous analytic and numerical work.

Non-equilibrium rate theory for polariton relaxation dynamics

We derive an analytic expression of the non-equilibrium Fermi's golden rule (NE-FGR) expression for a Holstein-Tavis-Cumming Hamiltonian, a universal model for many molecules collectively coupled to the optical cavity. These NE-FGR expressions capture the full-time-dependent behavior of the rate constant for transitions from polariton states to dark states. The rate is shown to be reduced to the well-known frequency domain-based equilibrium Fermi's golden rule (E-FGR) expression in the equilibrium and collective limit and is shown to retain the same scaling with the number of sites in non-equilibrium and non-collective cases. We use these NE-FGR to perform population dynamics with a time-non-local and time-local quantum master equation and obtain accurate population dynamics from the initially occupied upper or lower polariton states. Furthermore, NE-FGR significantly improves the accuracy of the population dynamics when starting from the lower polariton compared to the E-FGR theory, highlighting the importance of the non-Markovian behavior and the short-time transient behavior in the transition rate constant.

Probing the anharmonicity of vibrational polaritons with double-quantum two-dimensional infrared spectroscopy

Strong coupling between the molecular vibrations and electromagnetic fields of light confined to an infrared cavity leads to the formation of vibro-polaritons - quasi-particles thought to provide the means to control the rates of chemical reactions inside a dark cavity. Despite the mechanisms indicating how vibrational coupling to the vacuum fields can affect the reaction rates are still not well understood, it has been recently demonstrated that the formation of the polariton states alters the ultrafast relaxation dynamics of the strongly coupled system. The relaxation dynamics in molecules, which is known to be important for the chemical reactivity, is directed by anharmonic couplings involving multiple intra- and inter-molecular vibrational degrees of freedom. However, the impact of the molecular anharmonicity on the polariton states remains elusive. Some theoretical models, employed to interpret the experimental observations, assume that vibrational polaritons are harmonic. Others assume a certain anharmonicity of vibro-polaritons; however, to date, it has not been experimentally determined. Herein, we performed double-quantum two-dimensional third-order nonlinear infrared spectroscopy of the carbonyl stretching (C=O) vibrational modes in a thin film of polymethyl methacrylate polymer (PMMA) strongly coupled to the surface lattice resonances of the periodic arrays of half-wavelength infrared disk antennas. We found that, indeed, the mechanical anharmonicity of polaritons is very small. Quantitatively, our results place an upper bound on a polariton mechanical anharmonicity of 2 cm-1, compared with that of the C=O mode in a PMMA film of 15 cm-1. Thus, our results support previous assumptions regarding the harmonic character of vibro-polaritons.

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.

Modeling the near-field effect on molecular excited states using the discrete interaction model/quantum mechanical method

Strong light-matter interactions significantly modify the optical properties of molecules in the vicinity of plasmonic metal nanoparticles. Since the dimension of the plasmonic cavity approaches that of the molecules, it is critical to explicitly describe the nanoparticle junctions. In this work, we use the discrete interaction model/quantum mechanical (DIM/QM) method to model the coupling between the plasmonic near-field and molecular excited states. DIM/QM is a combined electrodynamics/quantum mechanical model that uses an atomistic description of the nanoparticle. We extend the DIM/QM method to include the local field effects in the sum-over-state formalism of time-dependent density functional theory. As a test of the method, we study the interactions between small organic chromophores and metal nanoparticles. In particular, we examine how the inclusion of multiple electronic transitions and intermolecular interactions modify the coupling between molecules and nanoparticles. Using the sum-over-state formalism of DIM/QM, we show that two-state models break down when the plasmon excitation is detuned from the molecular excitations. To gain further insight, we compare the simple coupled-dipole model (CDM) with the DIM/QM model. We find that CDM works well for simple systems but fails when going beyond the single molecule or single nanoparticle cases. We also find that the coupling depends strongly on the site of the nanoparticle in which the chromophore couples to. Our work suggests the importance of explicitly describing the cavity to capture the atomistic level local field environment in which the molecule strongly couples to.

Singly-excited resonant open quantum system Tavis-Cummings model with quantum circuit mapping

Tavis-Cummings (TC) cavity quantum electrodynamical effects, describing the interaction of N atoms with an optical resonator, are at the core of atomic, optical and solid state physics. The full numerical simulation of TC dynamics scales exponentially with the number of atoms. By restricting the open quantum system to a single excitation, typical of experimental realizations in quantum optics, we analytically solve the TC model with an arbitrary number of atoms with linear complexity. This solution allows us to devise the Quantum Mapping Algorithm of Resonator Interaction with N Atoms (Q-MARINA), an intuitive TC mapping to a quantum circuit with linear space and time scaling, whose N+1 qubits represent atoms and a lossy cavity, while the dynamics is encoded through 2N entangling gates. Finally, we benchmark the robustness of the algorithm on a quantum simulator and superconducting quantum processors against the quantum master equation solution on a classical computer.

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.

Dissociation dynamics of a diatomic molecule in an optical cavity

We study the dissociation dynamics of a diatomic molecule, modeled as a Morse oscillator, coupled to an optical cavity. A marked suppression of the dissociation probability, both classical and quantum, is observed for cavity frequencies significantly below the fundamental transition frequency of the molecule. We show that the suppression in the probability is due to the nonlinearity of the dipole function. The effect can be rationalized entirely in terms of the structures in the classical phase space of the model system.

Cavity-enabled enhancement of ultrafast intramolecular vibrational redistribution over pseudorotation

Vibrational strong coupling (VSC) between molecular vibrations and microcavity photons yields a few polaritons (light-matter modes) and many dark modes (with negligible photonic character). Although VSC is reported to alter thermally activated chemical reactions, its mechanisms remain opaque. To elucidate this problem, we followed ultrafast dynamics of a simple unimolecular vibrational energy exchange in iron pentacarbonyl [Fe(CO) 5 ] under VSC, which showed two competing channels: pseudorotation and intramolecular vibrational-energy redistribution (IVR). We found that under polariton excitation, energy exchange was overall accelerated, with IVR becoming faster and pseudorotation being slowed down. However, dark-mode excitation revealed unchanged dynamics compared with those outside of the cavity, with pseudorotation dominating. Thus, despite controversies around thermally activated VSC modified chemistry, our work shows that VSC can indeed alter chemistry through a nonequilibrium preparation of polaritons.