Polaron-Polariton Assisted Thermally Activated Superradiance

Introduction.— Molecular aggregates have attracted significant attention due to their unique optical properties and potential applications in light-emitting devices [1, 2, 3, 4, 5]. One of their most striking features is superradiance—the enhancement of collective spontaneous emission rate beyond that of an individual monomer [6]. Pioneering studies of molecular superradiance trace back to Kasha’s seminal works on excitons in aggregates, which revealed that Coulombic dipole-dipole interactions among monomers lead to distinct photophysical behaviors depending on the packing configuration of the monomers [7, 8, 9, 10]. Subsequent studies have extended Kasha’s framework to include exciton-phonon interactions [11, 12, 13, 3] and disorder [14, 15], both of which typically diminish superradiance. In addition, superradiance in molecular aggregates is generally suppressed as the temperature increases [16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26]. To the best of our knowledge, thermally activated superradiance has rarely been observed and is typically anticipated only in molecular aggregates where the bright state lies energetically above dark states [27, 28], as reported in H aggregates [29] and two-dimensional aggregates [30]. In contrast, this phenomenon is not expected in one-dimensional J aggregates, where such high-lying bright-state configurations are absent.

Recently, coupling molecules to confined electromagnetic fields, such as cavity photons or polaritons (photons dressed by dielectric environments), has emerged as a novel approach for modifying molecular processes, including energy transfer [31, 32, 33, 34, 35] and chemical reactions [36, 37, 38, 39, 40, 41]. In particular, polariton-assisted superradiance in a J aggregate has been theoretically predicted to yield anomalous giant rate enhancement, exceeding that in free space by roughly an order of magnitude [42]. However, the previous study has not accounted for the vibrational and thermal effects, which are known to suppress superradiance in molecular aggregates. To overcome these limitations, we develop a general theory of single-photon superradiance that describes molecular aggregates coupled to both phonons and polaritons at finite temperatures. This theory combines macroscopic quantum electrodynamics (MQED) [43, 44, 45, 46]—a first-principles quantization scheme for electromagnetic fields in inhomogeneous, dispersive and absorbing media—with the polaron quantum master equation approach [47, 48]. We show that, in contrast to free space, superradiance in molecular aggregates coupled to polaritons can be enhanced by increasing either the exciton-phonon coupling or the temperature. Our results reveal an unconventional thermally activated superradiance mechanism in one-dimensional J aggregates and open new directions for controlling molecular emission via dielectric engineering.

Model Hamiltonian.— We consider molecular aggregates, composed of a set of identical electronic two-level systems (TLSs), interacting simultaneously with phonons (originate from intramolecular and local environmental vibrational modes) and polaritons (photons dressed by the dielectric environment). Based on the multipolar coupling framework of MQED, the total Hamiltonian under the electric-dipole approximation can be expressed as

The first term describes the TLSs constituting the aggregate, where represents the transition frequency of each monomer, and () denotes the raising (lowering) operator for monomer , with and corresponding to the electronic ground and excited states of monomer , respectively. The second term describes the phonon (vibrational) modes, where and denote the frequency and the annihilation operator of the -th phonon mode that couples to monomer . The third term describes the interactions between the monomers and the phonon modes, where denotes the Huang-Rhys factor that characterizes the coupling between monomer and its phonon mode . We consider that each monomer couples individually to a separate phonon bath and that both the phonon bath properties and exciton-phonon coupling strengths are identical across all monomers; therefore, the notation can be further simplified as and , and the exciton-phonon coupling is determined by the phonon (vibrational) spectral density . The fourth term describes the polaritons (photons dressed by the dielectric environment), where denotes the annihilation operator of the bosonic vector fields (polariton fields). The fifth term describes the interactions between the monomers and polaritons, where denotes the transition dipole operator of monomer , and represents its corresponding transition dipole moment. The field operator is defined as , where the tensor is explicitly given by . The symbols , , and denote the permittivity of free space, the relative permittivity, and the speed of light in free space, respectively; the tensor denotes the dyadic Green’s function that satisfies macroscopic Maxwell’s equations (Maxwell’s equations in matter) , where is the three-dimensional identity matrix. The information regarding the exciton-polariton coupling is fully encoded within the generalized polariton spectral density , with elements .

Generalized theory of superradiance.— Our goal is to explore the combined effect of exciton-phonon and exciton-polariton couplings on the collective decay rates of the aggregate (excited) states. Following previous studies [49, 42], the aggregate states are defined as , where the local excited state of monomer is expressed as , with . The coefficients correspond to elements of the unitary matrix that diagonalizes the free-space dipole-dipole interaction matrix , whose matrix elements are explicitly given by for and for , with denoting the free-space dyadic Green’s function [50]. The state for which the coefficients exhibit the same sign for all is further referred to as the superradiant or bright state, owing to its maximal collective decay rate.

We focus on molecular aggregates that are weakly coupled to polaritons, and the interactions among the monomers are dominated by free-space dipole-dipole interactions. To derive the superradiance rate of a molecular aggregate, we first transform the total Hamiltonian in Eq. (1) into the polaron frame, simplifying the treatment of exciton-phonon interactions, and derive the equation of motion from the quantum Liouville equation. Note that the superradiance rate of a molecular aggregate cannot be obtained by a standard second-ordered polaron quantum master equation due to the strong intermolecular dipole-dipole interactions. Instead, we transform the quantum Liouville equation into a modified rotating-frame. Then, we treat phonons and polaritons as baths, apply the Born-Markov and secular approximations (weak coupling between molecules and polaritons ), derive a Lindblad-type quantum master equation, and finally identify the collective transition rate from the aggregate excited state to the ground state . A detailed derivation and discussion can be found in Supplemental Material (SM) [51], and the final result of is given by:

The modified polariton spectral densities are defined as and , which correspond to the polariton emission and absorption processes, respectively. The thermal occupation number of polaritons is given by , where and denote the Boltzmann constant and the temperature of the polariton bath, respectively. The Franck-Condon weighted density of states (FCWD) is the Fourier transform of the phonon time correlation function , where the phonon phase function is given by . Here, denotes the temperature of all the independent phonon baths. The transition energy of aggregate state is expressed as , where .

The collective transition rate recovers the results of several previous studies under various limiting conditions, including (i) the single-monomer case [52], (ii) the absence of exciton-phonon coupling [42], and (iii) the free-space scenario without exciton-phonon coupling [49]. This consistency supports the validity of our approach, and detailed demonstrations are provided in SM [51]. We note that Eq. (2) provides a robust description in the weak light-matter coupling regime, which is characterized by Markovian dynamics. Under strong light-matter coupling, the dynamics of the aggregate are instead expected to exhibit Rabi oscillations [53, 50].

System setup.— We consider a minimal model composed of two monomers with an energy gap of eV and a transition dipole moment magnitude of Debye, positioned above a silver surface in a head-to-tail (J aggregate) configuration, as depicted in Fig. 1(a). The monomers are separated by 1 nm, and the molecule-surface distance is set to 10 nm to prevent the intrinsic loss mechanisms in the metal from dominating the decay process [54] and to ensure the system remains within the weak exciton-polariton coupling regime [42]. A more detailed analysis of the exciton-polariton coupling regime can be found in the previous work [55]. In this setup, the monomers are coupled to a global polariton bath, while each monomer is individually coupled to an independent phonon (vibrational) bath, as illustrated in Fig. 1(b). The dielectric function of silver is modeled using an analytical expression [56] fitted to experimental data [57]. For the phonon spectral density, we employ a super-Ohmic form, i.e., , where is a dimensionless exciton-phonon coupling constant (not related to the exciton-polariton coupling strength), and denotes the cutoff frequency, fixed at cm^-1 throughout this study. Note that our theoretical framework applies to arbitrary dielectric environments and is not restricted to the super-Ohmic form of the phonon spectral density. This setup is chosen to reveal an anomalous superradiance behavior. We focus on the transition rate of the superradiant state, i.e., , where the expansion coefficients satisfy for all . This state corresponds to the lowest-energy excited state of a J aggregate. Moreover, we set the temperature .

Polaron-polariton effect on superradiance.— First, to explore the effect of polarons and polaritons on molecular superradiance, we vary the exciton-phonon coupling strength and place molecules above a silver surface to introduce exciton-polariton coupling. Fig. 2(a) presents the rate enhancement, defined as the transition rate of the dimer normalized by that of the monomer , as a function of at for systems with exciton-polariton coupling (silver surface) and without it (free space). In free space, the rate enhancement decreases monotonically as increases; in contrast, when above a silver surface, the enhancement initially increases before eventually decreasing. This initial rise in rate enhancement highlights a key distinction between the polariton-free and polariton-assisted superradiance, and appears to contradict the conventional understanding that exciton-phonon coupling generally suppresses superradiance. To investigate the origin of this anomalous behavior, we plot the transition rates of the dimer and monomer separately in Fig. 2(b). It is shown that the trend observed in the rate enhancement closely follows that of the dimer’s superradiance rate .

To better understand the influence of phonons on the superradiance rate, we decompose the FCWD in Eq. (2) into the zero-phonon line (ZPL) and the phonon sideband (PSB) [58, 59, 60, 61] as , with

where , reflecting the absence of correlations between different phonon baths that each molecule couples to. Note that the ZPL and PSB correspond to the zero-phonon and phonon-assisted transitions, respectively. For the treatment of , the translation symmetry of each monomer above the surface and the close proximity of the monomers allow us to define [42]. Therefore, at zero temperature, where , the collective transition rate in Eq. (2) simplifies to . We further define the zero-phonon and phonon-assisted transition rates as , such that the total transition rate can be expressed as

The absence of a prefactor in the phonon-assisted transition rate arises from the fact that each molecule is coupled to its own independent local phonon bath. For a dimer, . The overall behavior of in Fig. 2(b) thus reflects the combining effect of and , and we plot them separately in Fig. 2(c).

Figure 2(c) shows that, in free space, increasing naturally enhances the phonon-assisted transition rate but suppresses the zero-phonon transition rate , and the relative changes in these two processes are nearly equal in magnitude. Since contributes twice in , the overall superradiance rate is suppressed as increases. On the other hand, above the surface, and both exhibit non-monotonic behavior. Furthermore, by decomposing , Fig. 2(d) shows that the key distinction between the monomer and dimer responses lies in the monotonic decrease of with increasing when above the surface. This can be understood by the analytical solution , which shows that its magnitude decays exponentially as increases.

To elucidate the origin of the non-monotonic trend of , we examine the components of , including and , as shown in Figs. 2(e) and 2(f). As increases, the peak of shifts to higher energy due to the reduction in the free-space dipole-dipole interaction between the monomers, i.e., . This shift initially improves the resonance between and but eventually results in off-resonance. Although the magnitude of decreases monotonically with , the resonance-induced enhancement dominates at small , leading to the initial rise of . The rate in Fig. 2(c) exhibits a similar non-monotonic dependence with a slight difference that the enhancement persists over a broader range of . This distinction is explained by the profile of in Fig. 2(g): first, is broadband and thus overlaps a wider resonance window with than the sharply peaked ; second, the amplitude of initially grows as rises from to , while that of decreases monotonically.

The above results show that predicting the pattern of the polaron-polariton assisted superradiance is nontrivial. Nevertheless, we offer a transparent and physically intuitive way for analysis by evaluating: (i) the resonance condition between the FCWD and the polariton spectral density and (ii) the relative magnitude of the FCWD. In addition, since the dipole-dipole interaction , a larger transition dipole moment of the monomer typically leads to a more significant activation effect. Note that our theoretical framework can naturally accommodate disorder [15, 62, 63], e.g., by varying monomer transition frequencies to model static disorder. The observed phenomenon is expected to persist as long as the energetic disorder remains small relative to the intermolecular dipole-dipole interaction. Conversely, comparable disorder and intermolecular coupling may lead to richer phenomena worthy of further investigation.

Temperature control of superradiance.— While we have demonstrated the possibility of polaron-polariton assisted superradiance, directly controlling the exciton-phonon coupling strength in experiments remains challenging, hindering the verification of the predicted behavior. To overcome this difficulty, we propose an experimentally feasible alternative: tuning the temperature , which can also effectively modulates the exciton-phonon coupling strength. Specifically, the phonon phase function can be rewritten as , where serves as an effective and dependent exciton-phonon coupling strength function. Since increases monotonically with , raising effectively enhances the exciton-phonon coupling strength.

In Fig. 3, we examine the effect of on superradiance with . The trend observed in the rate enhancement in Fig. 3(a) is similar to that of the -dependent case. The initial rise in the rate enhancement for a dimer above a silver surface can be understood by , where we have assumed as . Both and initially increase with due to the thermally induced blue shifts of and [upper and lower panels of Fig. 3(b)]. This spectral shift first brings FCWD into resonance with [Fig. 2(e)], enhancing the superradiance rate, but eventually drives it out of resonance, resulting in a turnover behavior.

The non-monotonic temperature dependence presented in Fig. 3 reveals an anomalous thermally activated superradiance in one-dimensional J aggregates. Remarkably, this thermally induced enhancement becomes more pronounced and persists to higher temperatures, even up to room temperature, with increasing aggregate size N (see SM [51]). In addition, the emission rate can exhibit superlinear or sublinear scaling with N, or even decrease with increasing N, depending on the alignment between the monomer transition energy and the surface plasmon polariton resonance [42]. Moreover, the mechanism of this phenomenon differs fundamentally from that observed in earlier studies [29, 30]. In our case, the behavior stems from the interplay between phonons and polaritons, whereas the previously reported effect originates solely from the monomer packing configuration. Notably, our theoretical framework is sufficiently general to capture both mechanisms.

Conclusion.— By combining macroscopic quantum electrodynamics and a modified polaron quantum master equation approach, we develop a general theory of superradiance in molecular excitons coupled to both phonons and polaritons at finite temperature. According to our theoretical approach, we demonstrate that superradiance can be enhanced by increasing either the exciton-phonon coupling strength or the temperature. This finding predicts an anomalous thermally activated superradiance in one-dimensional J aggregates, which is uncommon in free space and unexplained by previous theoretical frameworks. Open questions remain, including quantifying the fraction of superradiant decay through radiative versus nonradiative channels in lossy media [64, 65] and maximizing the radiative fraction. Furthermore, it is also intriguing to understand how polaritons influence superradiance under external driving fields and multiphoton excitation [66, 67] through quantum master equations, as well as in large aggregates where spatial inhomogeneity becomes significant [68, 69]. We believe that this work highlights the potential of dielectric environment engineering for tailoring collective light emission in molecular aggregates.