Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical Review A

Stochastic-Chebyshev-expansion approach for the simulation of linear polariton spectroscopy in the collective-coupling regime

Braden M. Weight, Pengfei Huo

Published Jun 23, 2025• 1 citation
Publisher

Abstract

Molecular Polariton is becoming one of the leading directions to control a multitude of chemical and physical processes, such as charge transfer, selective bond breaking, and excited state dynamics. Accurately and efficiently simulating polariton properties under the collective coupling regimes (between $N$ molecules and the cavity mode) remains a central theoretical challenge. In this work, we use a stochastic resolution of the identity approach coupled with a Chebyshev expansion to compute various polariton photophysical properties, with a substantially reduced computational effort than would be needed for a direct diagonalization of the same Hamiltonian, which is often the bottleneck for such large dimensionality. Such quantities of interest are the total density of states (the eigenspectrum of the Hamiltonian) and the transmission spectrum (a probe of the photonic degrees of freedom), the latter of which is a direct observable in the experiment. We simulate the linear spectroscopy of molecule-cavity hybrid systems, specifically exploring the effects of the distribution and magnitude of molecular disorder for one, few, and many coupled molecules. We compare our numerical results to recent work, which formulated analytic expressions in the large-$N$ limit for the spectroscopic signals. We find that our results match those of the analytic results when $N= 100$, at which point we find that the collective effects for linear spectroscopy are converged.

DOI
10.1103/x5q3-tm1q

Related papers

Coherent multidimensional spectroscopy in polariton systems

Daniela Gallego-Valencia, Lars Mewes, Johannes Feist +1 more

The fast dynamics of molecular polaritonics is scrutinized theoretically through the implementation of two-dimensional spectroscopy protocols. We derive conceptually simple and computationally efficient formulas to calculate two-dimensional spectra for molecules, each of them modeled as a system of two electronic states including vibrational relaxation, immersed in an optical cavity, thus coupled to quantized radiation. Cavity photon losses and molecular relaxation are incorporated into the Hamiltonian dynamics to form an open quantum system that is solved through a master equation. In the collective case, the relaxation dynamics into dark states is revealed to be the crucial factor to explain the asymmetries in both the diagonal and cross peaks of two-dimensional spectra for long waiting times between excitation and detection, a feature shown by recent experiments. Our theoretical method provides a deeper insight in those processes that yield relevant signals in multidimensional molecular spectroscopy.

arXiv:2403.04734•Physical Review A•Jun 11, 2024•15 citations

Multidimensional semiclassical single- and double-quantum spectroscopy of anharmonic molecular polaritons

Michael Reitz, Harsh Bhakta, Wei Xiong +1 more

We present a general and efficient approach to compute phase-resolved multidimensional spectra of anharmonic molecular polaritons, based on a semiclassical evolution of the molecular Hamiltonian and cavity field in the large-$\mathcal{N}$ limit of many molecules coupled to a confined photonic mode. By systematically expanding the response in both amplitudes and phases of the input fields, our method enables a transparent and computationally simple construction of phase-cycled two-dimensional single- and double-quantum polariton spectra from the underlying nonlinear signal components. Here, phase cycling acts as an analogue of phase matching with oblique pulses, allowing for the isolation of the contributing nonlinear pathways in Liouville space. We specialize to vibrational polaritons and benchmark the method through direct comparison with experimentally measured single-quantum spectra, providing an explanation for the longstanding puzzle of the polariton bleach effect observed at short waiting times. Further, we show how the imprint of various types of anharmonicities on the double-excitation manifold can be directly probed and analyzed through double-quantum coherence spectroscopy. Taken together, our results establish a practical and powerful framework for the modeling and interpretation of nonlinear spectroscopic experiments on strongly coupled light-matter platforms and for guiding the design of cavity-enhanced molecular platforms.

quant-pharXiv:2604.21158•Apr 24, 2026

Cavity Quantum Electrodynamics of Molecular Systems : Polarization Effects, Multiple Excitations, and Vibrational Selectivity

Garcia Borges, Lucas

The field of Polaritonic Chemistry studies molecules strongly coupled to the light field of optical cavities. The resulting modification of the potential energy landscape has been demonstrated to enable control over the chemical and physical properties of materials and molecular ensembles. This thesis extends the Tavis–Cummings model to describe molecules under strong light–matter coupling, incorporating static dipole moments and dipole self-energy contributions derived from the Pauli–Fierz Hamiltonian. These additional terms are shown to be essential for reproducing cavity-modified molecular dynamics and spectra. Simulations on the nuclear wavefunction propagation of MgH+ molecules demonstrate the impact of these terms, while an effective two-level system model accurately captures the behavior of large ensembles. Beyond single excitations, the thesis also explores the higher excitation manifolds of the Tavis-Cummings model, which reveals that allowing for more than a single excitation makes the reaction of the involved polaritons entropically more favorable. Open-system simulations highlight how the coherence of the system's initial state governs decay pathways and enables long-lived energy storage. Finally, a multilevel vibrational model coupled to a two-mode cavity is introduced, where a Raman scheme enables selective mixing of vibrational states. The simulation results show that continuous pumping can drive molecules into specific vibrational states while suppressing electronic excitation, both for single molecules and small ensembles. Together, these results provide new theoretical insights into the mechanisms by which optical cavities reshape molecular energy landscapes. The results suggest approaches for controlling photonic, electronic, and vibrational dynamics in strongly coupled systems, contributing to the theoretical foundation of polariton chemistry.

•Jan 1, 2025•0 citations

Molecular Electron Transfer in Optical Cavities: From Excitonic to Vibronic Polaritons

Takumi Hidaka, Tomohiro Fukushima, Nguyen Thanh Phuc

Strong coupling between molecular excitations and quantized electromagnetic fields in optical cavities provides a powerful means to control the physical and chemical properties of molecular systems. Here, we study electron transfer (ET) dynamics in cavity-coupled molecules using the numerically exact hierarchical equations of motion (HEOM) method, which captures nonperturbative and non-Markovian effects beyond standard perturbative theories. We identify distinct resonance and collective effects associated with polariton formation and show that the ET rate saturates in the strong-coupling regime, a feature not captured by perturbative approaches. We further extend the cavity-modified ET model by incorporating the nuclear-coordinate dependence of molecular electric dipole moments, which gives rise to a three-body interaction involving molecular electronic and vibrational degrees of freedom and cavity photons. This vibronic polariton formation leads to non-monotonic, oscillatory dependencies of the ET rate on the light-matter coupling strength and cavity frequency, which we attribute to quantum interference among multiple transfer pathways. These findings establish cavity-modified electron transfer as a multichannel quantum process governed by the interplay of electronic, vibrational, and photonic degrees of freedom.

physics.chem-phCross-listedarXiv:2602.23748•Mar 2, 2026

For molecular polaritons, disorder and phonon timescales control the activation of dark states in the thermodynamic limit

Tianchu Li, Pranay Venkatesh, Qiang Shi +1 more

Collective light-matter systems host an extensive manifold of dark states whose role in the emergence of thermodynamic behavior remains poorly understood, especially in the presence of disorder and structured environments. Here, we develop a hybrid matrix product state-hierarchical equations of motion (MPS-HEOM) approach that enables numerically exact simulations of polariton dynamics from a few emitters to the thermodynamic limit under both static and dynamic disorder. This allows us, for the first time, to provide a quantitative and operational answer to the long-standing question of what is the minimum system size required to reach the thermodynamic limit in collective polaritonic systems. By introducing a convergence scale, $N_{T}$, i.e., the number of molecules required for the photonic dynamics to reach the thermodynamic limit, we show that dynamic disorder generally poses a greater computational challenge than static disorder. We attribute this behavior to the suppression of collective light-matter dynamics by disorder, which dynamically activates non-collective degrees of freedom. We further find that $N_{T}$ exhibits a turnover behavior as the bath becomes more Markovian, as the bath timescales regulate bright-to-dark energy transfer and the involvement of dark and gray states. Hence, phonon timescales control both the breakdown of collective behavior and the growth of $N_{T}$. Our results establish the suppression of collective behavior as the key mechanism governing thermodynamic convergence in disordered light-matter systems.

physics.chem-phCross-listedarXiv:2603.06868•Mar 10, 2026

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

Jorge A. Campos-Gonzalez-Angulo, Joel Y. Yuen-Zhou

Confined electromagnetic modes strongly couple to collective excitations in ensembles of quantum emitters, producing light–matter hybrid states known as polaritons. Under such conditions, the discrete multilevel spectrum of molecular systems offers an appealing playground for exploring multiphoton processes. This work contrasts predictions from the Tavis–Cummings model in which the material is a collection of two-level systems, with the implications of considering additional energy levels with harmonic and anharmonic structures. We discuss the exact eigenspectrum, up to the second excitation manifold, of an arbitrary number N of oscillators collectively coupled to a single cavity mode in the rotating-wave approximation. Elaborating on our group-theoretic approach [New J. Phys. 23, 063081 (2021)], we simplify the brute-force diagonalization of N2 × N2 Hamiltonians to the eigendecomposition of, at most, 4 × 4 matrices for arbitrary N. We thoroughly discuss the eigenstates and the consequences of weak and strong anharmonicities. Furthermore, we find resonant conditions between bipolaritons and anharmonic transitions where two-photon absorption can be enhanced. Finally, we conclude that energy shifts in the polaritonic states induced by anharmonicities become negligible for large N. Thus, calculations with a single or few quantum emitters qualitatively fail to represent the nonlinear optical response of the collective strong coupling regime. Our work highlights the rich physics of multilevel anharmonic systems coupled to cavities absent in standard models of quantum optics. We also provide concise tabulated expressions for eigenfrequencies and transition amplitudes, which should serve as a reference for future spectroscopic studies of molecular polaritons.

arXiv:2202.01433•The Journal of Chemical Physics•Apr 5, 2022•41 citations
Molecular Polariton is becoming one of the leading directions to control a multitude of chemical and physical processes, such as charge transfer, selective bond breaking, and excited state dynamics. Accurately and efficiently simulating polariton properties under the collective coupling regimes (between NNN molecules and the cavity mode) remains a central theoretical challenge. In this work, we use a stochastic resolution of the identity approach coupled with a Chebyshev expansion to compute various polariton photophysical properties, with a substantially reduced computational effort than would be needed for a direct diagonalization of the same Hamiltonian, which is often the bottleneck for such large dimensionality. Such quantities of interest are the total density of states (the eigenspectrum of the Hamiltonian) and the transmission spectrum (a probe of the photonic degrees of freedom), the latter of which is a direct observable in the experiment. We simulate the linear spectroscopy of molecule-cavity hybrid systems, specifically exploring the effects of the distribution and magnitude of molecular disorder for one, few, and many coupled molecules. We compare our numerical results to recent work, which formulated analytic expressions in the large-NNN limit for the spectroscopic signals. We find that our results match those of the analytic results when N=100N= 100N=100, at which point we find that the collective effects for linear spectroscopy are converged.
The fast dynamics of molecular polaritonics is scrutinized theoretically through the implementation of two-dimensional spectroscopy protocols. We derive conceptually simple and computationally efficient formulas to calculate two-dimensional spectra for molecules, each of them modeled as a system of two electronic states including vibrational relaxation, immersed in an optical cavity, thus coupled to quantized radiation. Cavity photon losses and molecular relaxation are incorporated into the Hamiltonian dynamics to form an open quantum system that is solved through a master equation. In the collective case, the relaxation dynamics into dark states is revealed to be the crucial factor to explain the asymmetries in both the diagonal and cross peaks of two-dimensional spectra for long waiting times between excitation and detection, a feature shown by recent experiments. Our theoretical method provides a deeper insight in those processes that yield relevant signals in multidimensional molecular spectroscopy.
We present a general and efficient approach to compute phase-resolved multidimensional spectra of anharmonic molecular polaritons, based on a semiclassical evolution of the molecular Hamiltonian and cavity field in the large-N\mathcal{N}N limit of many molecules coupled to a confined photonic mode. By systematically expanding the response in both amplitudes and phases of the input fields, our method enables a transparent and computationally simple construction of phase-cycled two-dimensional single- and double-quantum polariton spectra from the underlying nonlinear signal components. Here, phase cycling acts as an analogue of phase matching with oblique pulses, allowing for the isolation of the contributing nonlinear pathways in Liouville space. We specialize to vibrational polaritons and benchmark the method through direct comparison with experimentally measured single-quantum spectra, providing an explanation for the longstanding puzzle of the polariton bleach effect observed at short waiting times. Further, we show how the imprint of various types of anharmonicities on the double-excitation manifold can be directly probed and analyzed through double-quantum coherence spectroscopy. Taken together, our results establish a practical and powerful framework for the modeling and interpretation of nonlinear spectroscopic experiments on strongly coupled light-matter platforms and for guiding the design of cavity-enhanced molecular platforms.
The field of Polaritonic Chemistry studies molecules strongly coupled to the light field of optical cavities. The resulting modification of the potential energy landscape has been demonstrated to enable control over the chemical and physical properties of materials and molecular ensembles. This thesis extends the Tavis–Cummings model to describe molecules under strong light–matter coupling, incorporating static dipole moments and dipole self-energy contributions derived from the Pauli–Fierz Hamiltonian. These additional terms are shown to be essential for reproducing cavity-modified molecular dynamics and spectra. Simulations on the nuclear wavefunction propagation of MgH+ molecules demonstrate the impact of these terms, while an effective two-level system model accurately captures the behavior of large ensembles. Beyond single excitations, the thesis also explores the higher excitation manifolds of the Tavis-Cummings model, which reveals that allowing for more than a single excitation makes the reaction of the involved polaritons entropically more favorable. Open-system simulations highlight how the coherence of the system's initial state governs decay pathways and enables long-lived energy storage. Finally, a multilevel vibrational model coupled to a two-mode cavity is introduced, where a Raman scheme enables selective mixing of vibrational states. The simulation results show that continuous pumping can drive molecules into specific vibrational states while suppressing electronic excitation, both for single molecules and small ensembles. Together, these results provide new theoretical insights into the mechanisms by which optical cavities reshape molecular energy landscapes. The results suggest approaches for controlling photonic, electronic, and vibrational dynamics in strongly coupled systems, contributing to the theoretical foundation of polariton chemistry.
Strong coupling between molecular excitations and quantized electromagnetic fields in optical cavities provides a powerful means to control the physical and chemical properties of molecular systems. Here, we study electron transfer (ET) dynamics in cavity-coupled molecules using the numerically exact hierarchical equations of motion (HEOM) method, which captures nonperturbative and non-Markovian effects beyond standard perturbative theories. We identify distinct resonance and collective effects associated with polariton formation and show that the ET rate saturates in the strong-coupling regime, a feature not captured by perturbative approaches. We further extend the cavity-modified ET model by incorporating the nuclear-coordinate dependence of molecular electric dipole moments, which gives rise to a three-body interaction involving molecular electronic and vibrational degrees of freedom and cavity photons. This vibronic polariton formation leads to non-monotonic, oscillatory dependencies of the ET rate on the light-matter coupling strength and cavity frequency, which we attribute to quantum interference among multiple transfer pathways. These findings establish cavity-modified electron transfer as a multichannel quantum process governed by the interplay of electronic, vibrational, and photonic degrees of freedom.
Collective light-matter systems host an extensive manifold of dark states whose role in the emergence of thermodynamic behavior remains poorly understood, especially in the presence of disorder and structured environments. Here, we develop a hybrid matrix product state-hierarchical equations of motion (MPS-HEOM) approach that enables numerically exact simulations of polariton dynamics from a few emitters to the thermodynamic limit under both static and dynamic disorder. This allows us, for the first time, to provide a quantitative and operational answer to the long-standing question of what is the minimum system size required to reach the thermodynamic limit in collective polaritonic systems. By introducing a convergence scale, NTN_{T}NT​, i.e., the number of molecules required for the photonic dynamics to reach the thermodynamic limit, we show that dynamic disorder generally poses a greater computational challenge than static disorder. We attribute this behavior to the suppression of collective light-matter dynamics by disorder, which dynamically activates non-collective degrees of freedom. We further find that NTN_{T}NT​ exhibits a turnover behavior as the bath becomes more Markovian, as the bath timescales regulate bright-to-dark energy transfer and the involvement of dark and gray states. Hence, phonon timescales control both the breakdown of collective behavior and the growth of NTN_{T}NT​. Our results establish the suppression of collective behavior as the key mechanism governing thermodynamic convergence in disordered light-matter systems.
Confined electromagnetic modes strongly couple to collective excitations in ensembles of quantum emitters, producing light–matter hybrid states known as polaritons. Under such conditions, the discrete multilevel spectrum of molecular systems offers an appealing playground for exploring multiphoton processes. This work contrasts predictions from the Tavis–Cummings model in which the material is a collection of two-level systems, with the implications of considering additional energy levels with harmonic and anharmonic structures. We discuss the exact eigenspectrum, up to the second excitation manifold, of an arbitrary number N of oscillators collectively coupled to a single cavity mode in the rotating-wave approximation. Elaborating on our group-theoretic approach [New J. Phys. 23, 063081 (2021)], we simplify the brute-force diagonalization of N2 × N2 Hamiltonians to the eigendecomposition of, at most, 4 × 4 matrices for arbitrary N. We thoroughly discuss the eigenstates and the consequences of weak and strong anharmonicities. Furthermore, we find resonant conditions between bipolaritons and anharmonic transitions where two-photon absorption can be enhanced. Finally, we conclude that energy shifts in the polaritonic states induced by anharmonicities become negligible for large N. Thus, calculations with a single or few quantum emitters qualitatively fail to represent the nonlinear optical response of the collective strong coupling regime. Our work highlights the rich physics of multilevel anharmonic systems coupled to cavities absent in standard models of quantum optics. We also provide concise tabulated expressions for eigenfrequencies and transition amplitudes, which should serve as a reference for future spectroscopic studies of molecular polaritons.