The coupling between propagating surface plasmon polaritons (SPPs) in Au nanostructures and the exciton transitions of cyanine dye J-aggregates has been examined using leakage radiation microscopy. Real space images of the nanostructures give the propagation lengths of the leaky SPP modes, and Fourier space images yield their dispersion curves. The dispersion curves show an avoided crossing when the structures are coated with J-aggregates, with a Rabi splitting of approximately 30 meV. The lifetimes of the coupled states were calculated by combining the measured propagation lengths with the group velocities obtained from the dispersion curves. The SPP lifetimes decrease from ∼50 fs for the bare Au nanostructures, to ∼10 fs in the avoided crossing region for the coupled J-aggregate/Au nanostructure system. Analytical Holstein-Tavis-Cummings model calculations and finite element simulations of the coupled system show that the decrease in lifetime is primarily due to energy dissipation into dark states associated with the J-aggregates.

Exciton transport in molecular aggregates is a fundamental process governing the performance of organic optoelectronics and light-harvesting systems. While most theoretical studies have emphasized long-time transport behavior, recent advances in ultrafast spectroscopy have brought into focus the short-time regime, in which exciton motion remains ballistic on femtosecond-to-picosecond timescales. In this work, we develop an analytical framework for short-time exciton dynamics in a one-dimensional lattice subject to both on-site energetic (diagonal) disorder and intermolecular coupling (off-diagonal) fluctuations. Utilizing the reciprocal-space analysis, we derive closed-form expressions for the first and second spatial moments considering both localized excitation and moving Gaussian initial conditions. Our analytical and numerical results show that, while the long-time dynamics are influenced by diagonal disorder, the short-time ballistic expansion is governed primarily by off-diagonal disorder. Crucially, we reveal a synergistic interplay between the average intermolecular coupling and the off-diagonal coupling disorder strength, demonstrating that they contribute equivalently to short-time exciton transport. Moreover, we integrate this generic disorder model with a realistic molecular system within the framework of macroscopic quantum electrodynamics, thereby providing a theoretical foundation for characterizing and optimizing ultrafast energy flow of disordered molecular aggregates in complex dielectric media.

We present numerically exact inchworm quantum Monte Carlo results for the real-time dynamics of the spin polarization in the sub-Ohmic spin-boson model at finite temperature. We focus in particular on the localization and coherence behavior of the model, extending our previous study at low temperature [Goulko et al., Phys. Rev. Lett. 134, 056502 (2025)]. As the temperature increases, the system becomes less localized and less coherent. The loss of coherence, which is controlled by two independent mechanisms-a smooth damping-driven crossover and a sharp frequency-driven transition-exhibits a nontrivial temperature dependence. While both types of coherence loss occur at lower coupling in the high temperature regime, the frequency exhibits a sharper drop at high temperatures and this drop is observed for all values of the sub-Ohmic exponent, in contrast to the zero-temperature case. We discuss the full temperature-dependent dynamical phase diagram of the system and the interplay between coherence and localization across a wide range of physical parameters.

Simulating the dynamics of molecular excitons in complex nanophotonic environments requires integrating rigorous electromagnetic simulations with accurate treatments of open quantum system dynamics. In this work, we develop MQED-QD (Macroscopic Quantum Electrodynamics for Quantum Dynamics), a robust computational package for simulating exciton dynamics in arbitrary dielectric and plasmonic environments. Based on the MQED framework, the package offers a unified workflow for constructing the dyadic Green's functions from classical electromagnetic solvers, parametrizing quantum master equations, and propagating the time evolution to determine the molecular subsystem's dynamical properties. To demonstrate the package's capabilities, we simulate exciton transport within a one-dimensional molecular chain near a silver nanostructure, including benchmarking against planar surfaces and exploring the influence of silver nanorods. Our results reveal that surface plasmon polaritons on nanorods drastically enhance long-range dipole-dipole interactions, accelerating exciton delocalization and yielding higher participation ratios compared to planar geometries. MQED-QD provides a powerful, open-source package that facilitates the rational design of nanoscale architectures by elucidating accurate molecular exciton dynamics in conjunction with nanophotonics and plasmonics.

Strong light-matter interactions between molecules and quantized electromagnetic fields inside an optical cavity open up novel possibilities, although inevitably influenced by disorder, an inherent attribute of realistic molecular systems. Here, we explore the steady-state optical response of molecular emitters within a lossy cavity, with a focus on the combined effects of static and dynamic disorders, as frequently observed in solution-phase experiments. By analyzing the transmission spectra, molecular energy change, and dark state population, we uncover the contrasting effects of static and dynamic disorders on the dark state population and its interplay with polariton states. We find that the Rabi splitting exhibits an inversion with increasing disorder strength, where the maximum splitting is determined by the interplay of static and dynamic disorders. Furthermore, we identify a dark state-induced polariton linewidth narrowing, revealing a mechanism distinct from motional narrowing induced by frequency fluctuations. These mechanistic insights highlight the critical role of dark states, establishing a foundation for future developments in the fields of polariton chemistry and strong coupling spectroscopy.

We investigate the real-time dynamics of the sub-Ohmic spin-boson model across a broad range of coupling strengths, using the numerically exact inchworm quantum Monte Carlo algorithm. From short- and intermediate-time dynamics starting from an initially decoupled state, we extract signatures of the zero-temperature quantum phase transition between localized and delocalized states. We show that the dynamical phase diagram thus obtained differs from the equilibrium phase diagram in both the values of critical couplings and the associated values of the critical exponent. We also identify and quantitatively analyze two competing mechanisms for the crossover between coherent oscillations and incoherent decay. Deep in the sub-Ohmic regime, the crossover is driven by the damping of the oscillation amplitude, while closer to the Ohmic regime the oscillation frequency itself drops sharply to zero at large coupling.

The notion of strong light-matter coupling is typically associated with the observation of Rabi splitting, corresponding to the formation of the hybrid light-matter states known as polaritons. However, this relationship is derived from the assumption that disorder can be ignored or acts as a perturbative effect. Contrary to the conventional treatment of disorder effects, we investigated the impact of strong disorder on the absorption spectrum by developing a nonperturbative effective model combined with classical electrodynamics simulations. Intriguingly, we find that strong disorder leads to an enhanced spectral splitting that closely resembles Rabi splitting, yet originates from a fundamentally different mechanism as induced by the dark modes. Specifically, we examine a disordered molecular ensemble in proximity to a plasmonic nanodisk and demonstrate disorder-induced spectral splitting in the absorption spectrum. This conclusion raises a controversial issue, suggesting that both polaritons (which dominate in the strong coupling regime) and dark modes (which dominate in the strong disorder regime) can lead to spectral splitting, and one cannot distinguish them solely based on the steady-state absorption spectrum.

Exciton transport in molecular aggregates with magic-angle orientation is expected to be strongly suppressed due to their negligible dipole-dipole interactions. However, recent reports show that light-matter interactions can significantly enhance exciton transport, attributed to the effective long-range coupling mediated by the photonic fields. To elucidate their interplay, we employ the macroscopic quantum electrodynamics framework to simulate exciton transport within a chromophore array arranged in a magic-angle configuration in proximity to a silver surface. Our results show a significant enhancement of the exciton mean square velocity that is robust across variations in chromophore-surface separation, intermolecular distance, and molecular transition frequency. Furthermore, based on the image-dipole method, we derive analytical expressions that agree well with numerical simulations, revealing the enhancement's origin in the near-field coupling term as induced by the radiative scattering at the metallic surface. More importantly, we observe nontrivial differences in the power-law scaling of the mean square velocity near metallic surfaces compared to free space. Our findings highlight the potential to control exciton transport by designing coupled exciton-photon systems and engineering the dielectric environments.

Fabry-Perot microcavities can strongly enhance interactions between light and molecules, leading to the formation of hybrid light-matter states known as polaritons. Polaritons possess much smaller effective masses and much larger group velocities when the molecules are resonant with cavity modes that have finite (non-zero) in-plane wavevectors, giving rise to the possibilities of long-range and ultrafast ballistic transport. In this paper, we present results of numerical simulations of the ultrafast ballistic transport phenomenon in real space and time during and after initialization with a short, spatially localized pulse. We find that the transport of the molecular excitons as induced by the external light field is synchronized with the evolution of the enhanced and localized electromagnetic field inside the cavity. Moreover, the synchronized transport rate is in good agreement with the group velocities predicted from a calculated dispersion relation across a wide range of frequencies. These simulations provide an intuitive tool for understanding the collective motion of light and excitons and helps to better understand how experimental observations of polaritons should be interpreted.

Strong light-matter coupling within an optical cavity leverages the collective interactions of molecules and confined electromagnetic fields, giving rise to the possibilities of modifying chemical reactivity and molecular properties. While collective optical responses, such as enhanced Rabi splitting, are often observed, the overall effect of the cavity on molecular systems remains ambiguous for a large number of molecules. In this paper, we investigate the non-adiabatic electron transfer process in electron donor-acceptor pairs influenced by collective excitation and local molecular dynamics. Using the timescale difference between reorganization and thermal fluctuations, we derive analytical formulas for the electron transfer rate constant and the polariton relaxation rate. These formulas apply to any number of molecules (N) and account for the collective effect as induced by cavity photon coupling. Our findings reveal a non-monotonic dependence of the rate constant on N, which can be understood by the interplay between electron transfer and polariton relaxation. As a result, the cavity-induced quantum yield increases linearly with N for small N (as predicted by a simple Dicke model) but shows a turnover and suppression for large N. We also interrelate the thermal bath frequency and the number of molecules, suggesting the optimal number for maximizing enhancement. The analysis provides an analytical insight for understanding the collective excitation of light and electron transfer, helping to predict the optimal condition for effective cavity-controlled chemical reactivity.

We investigate the effect of a cavity on nonlinear two-photon transitions of a molecular system and we analyze how such an effect depends on the cavity quality factor, the field enhancement, and the possibility of dephasing. We find that the molecular response to strong light fields in a cavity with a variable quality factor can be understood as arising from a balance between (i) the ability of the cavity to enhance the field of an external probe and promote multiphoton transitions more easily and (ii) the fact that the strict selection rules on multiphoton transitions in a cavity support only one resonant frequency within the excitation range. Although our simulations use a classical level description of the radiation field (i.e., we solve Maxwell-Bloch or Maxwell-Liouville equations within the Ehrenfest approximation for the field-molecule interaction), based on experience with this level of approximation in the past studies of plasmonic and polaritonic systems, we believe that our results are valid over a wide range of external probing.

The interplay between collective optical response and molecular static and dynamic disorder is studied using simple effective Hamiltonians for an ensemble of two-level emitters inside and outside a single-mode cavity. We model environmental disorder by randomly modulating the molecular transition frequencies and the coupling between the emitters and the electromagnetic field. We also consider effects of intermolecular interactions and orientational disorder. We investigate how these effects lead to new features in the steady state absorption (outside the cavity), transmission spectra (inside the cavity) and the yield of local molecular processes such as a unimolecular reaction. Outside the cavity, the collective behavior is manifested in the linewidth of the steady state absorption, the emission spectrum, and the local chemical yield. Inside the cavity, however, the collective behavior primarily determines the Rabi splitting. Effects of intermolecular interactions under orientational disorder are also studied. For the most part, for all types of disorder, if we increase disorder, we find a reduction in the collective nature of the molecular response (smaller effective N), and therefore, the Rabi splitting contraction occurs with orientational disorder. Moreover, we find that static disorder is more destructive to collective behavior than dynamic disorder.

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

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

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

We model a collection of N two-level systems (TLSs) coupled to a multimode cavity via Meyer-Miller-Stock-Thoss (MMST) dynamics, sampling both electronic and photonic zero-point energies (ZPEs) and propagating independent trajectories in Wigner phase space. By investigating the ground-state stability of a single TLS, we use MMST dynamics to separately study both electronic ZPE effects (which would naively lead to the breakdown of the electronic ground state) and photonic ZPE effects (which would naively lead to spontaneous absorption). By contrast, including both effects, i.e., sampling both electronic and photonic ZPEs, leads to the dynamical stability of the electronic ground state. Therefore, MMST dynamics provides a practical way to identify the contributions of self-interaction and vacuum fluctuations. More importantly, we find that MMST dynamics can predict accurate quantum dynamics for both electronic populations and electromagnetic field intensity in the high saturation limit. For a single TLS in a cavity, MMST dynamics correctly predicts the initial exponential decay of spontaneous emission, Poincare recurrences, and the positional dependence of a spontaneous emission rate. For an array of N equally spaced TLSs with only one TLS excited initially, MMST dynamics correctly predicts the modification of spontaneous emission rate as a function of the spacing between TLSs. Finally, MMST dynamics also correctly models Dicke's superradiance and subradiance, i.e., the dynamics when all TLSs are excited initially, including the correct quantum statistics for the delay time (as found by counting trajectories, for which a full quantum simulation is hard to achieve). Therefore, this work raises the possibility of simulating large-scale collective light-matter interactions with methods beyond mean-field theory.

We investigate different approaches to derive the proper Floquet-based quantum-classical Liouville equation (F-QCLE) for laser-driven electron-nuclear dynamics. The first approach projects the operator form of the standard QCLE onto the diabatic Floquet basis and then transforms to the adiabatic representation. The second approach directly projects the QCLE onto the Floquet adiabatic basis. Both approaches yield a form that is similar to the usual QCLE with two modifications: (1) The electronic degrees of freedom are expanded to infinite dimension and (2) the nuclear motion follows Floquet quasi-energy surfaces. However, the second approach includes an additional cross derivative force due to the dual dependence on time and nuclear motion of the Floquet adiabatic states. Our analysis and numerical tests indicate that this cross derivative force is a fictitious artifact, suggesting that one cannot safely exchange the order of Floquet state projection with adiabatic transformation. Our results are in accord with similar findings by Izmaylov et al., [J. Chem. Phys. 140, 084104 (2014)] who found that transforming to the adiabatic representation must always be the last operation applied, although now we have extended this result to a time-dependent Hamiltonian. This paper and the proper derivation of the F-QCLE should lay the basis for further improvements of Floquet surface hopping.

Quantum nondemolition (QND) measurements of photons is a much pursued endeavor in the field of quantum optics and quantum information processing. Here, we propose a hybrid optoelectromechanical platform that integrates a cavity system with a hybrid electromechanical probe for QND photon counting. Building upon a mechanical-mode-mediated nonperturbative electro-optical dispersive coupling, our protocol performs the QND photon counting measurement by means of the current-voltage characteristics of the probe. In particular, we show that the peak voltage shift of the differential conductance is linearly dependent on the photon occupation number, thus providing a sensitive measure of the photon number, especially in the strong optomechanical coupling regime. Given that our proposed hybrid system is compatible with state-of-the-art experimental techniques, we discuss its implementations and anticipate applications in quantum optics and polariton physics.