Not dark yet for strong light-matter coupling to accelerate singlet fission dynamics.

Polaritons are unique hybrid light-matter states that offer an alternative way to manipulate chemical processes. In this work, we show that singlet fission dynamics can be accelerated under strong light-matter coupling. For superexchange-mediated singlet fission, state mixing speeds up the dynamics in cavities when the lower polariton is close in energy to the multiexcitonic state. This effect is more pronounced in non-conventional singlet fission materials in which the energy gap between the bright singlet exciton and the multiexcitonic state is large ( > 0.1 eV). In this case, the dynamics is dominated by the polaritonic modes and not by the bare-molecule-like dark states, and, additionally, the resonant enhancement due to strong coupling is robust even for energetically broad molecular states. The present results provide a new strategy to expand the range of suitable materials for efficient singlet fission by making use of strong light-matter coupling.

Introduction

The implications that strong light-matter coupling can have in chemistry have recently raised a lot of interest.1, 2, 3, 4, 5, 6 This is because it offers an unconventional way to manipulate chemical processes by modifying the energy landscape as well as the dynamics., When an ensemble of molecules interacts with a confined light mode, new eigenstates of the system emerge in the strong coupling regime., This happens when the light-matter interaction exceeds the intrinsic decay rates of both the molecular excitations and the cavity photons. The new eigenstates of the system consist of two hybrid light-matter states known as polaritons and a manifold of dark states, superpositions of the molecular excitations that do not couple to the photon mode. The main difference between both sets of states is that polaritons are delocalized thanks to their cavity photon contribution, while dark states usually behave similarly to uncoupled single-molecule excitons.

In the past decade, there has been a lot of effort devoted to understanding the effects of electronic strong coupling in molecular systems, e.g., the modification of potential energy surfaces,,,, conical intersections,11, 12, 13, 14, 15 and electron and energy-transfer phenomena.16, 17, 18, 19, 20, 21, 22 The finite lifetime of the cavity photons23, 24, 25, 26 and the presence of a dense dark-state manifold27, 28, 29, 30, 31, 32, 33, 34, 35 are key to understanding polaritonic chemistry phenomena. As already noted, dark states may wash out polaritonic effects in setups with collective, i.e., many-molecule, strong coupling, as there is a macroscopic number of them compared with only two polaritons per cavity mode.36, 37, 38 For instance, there is certain controversy on strong light-matter coupling effects on thermally activated delayed fluorescence (TADF), in which a triplet state repopulates the singlet state responsible for the delayed emission., Some authors have suggested that there are limitations on polaritonic effects on reverse intersystem crossing., This is because the initial state is localized on a single molecule while a polariton is delocalized over the entire molecular ensemble that couples to the cavity photon. As further discussed below, the rate for population relaxation from a localized to a delocalized state is penalized by a factor , where N is the number of entities participating in the delocalized state.,,, In polariton-assisted TADF, direct population transfer from the triplet state to the lower polariton () then competes with the much faster process of relaxation between localized single-molecule states; that is, from the triplet state to the singlet dark-state manifold., However, the reverse process is not penalized when there is a manifold of N localized final states, as the factor in the rate to any single state is compensated by the number of states N.

Singlet fission is a downconversion photophysical reaction in which a spin-singlet exciton splits into two independent spin-triplet states (Equation 1).44, 45, 46 The recent interest in this unique phenomenon has been driven by its potential capacity to overcome the Shockley-Queisser limit for the efficiency of single junction solar cells.48, 49, 50 It is well accepted that the singlet fission process involves the generation of a multiexcitonic intermediate, the zero-spin triplet-pair state (),, which eventually splits into two uncoupled triplet states. Commonly, the formation of the state (first step in Equation 1) is the rate-limiting process in singlet fission and is thus the main subject of study in the field.

In contrast to TADF under strong coupling, a polaritonic mode could potentially be the initial state instead of the final one in a singlet fission process. Previous theoretical works investigating singlet fission under strong coupling have mainly focused on the single-molecule case,53, 54, 55 while studies of the collective situation are still scarce. There are no clear guidelines yet on if, how, and when singlet fission can benefit from the presence of polaritons. This is the main motivation for the current work exploring cavity-modified singlet fission. We focus on the situation of many-molecule strong coupling, which is experimentally much easier to achieve than few-molecule strong coupling (only possible in strongly subwavelength plasmonic nanocavities), and is simultaneously more promising for leading to practical light-harvesting devices in the strong coupling regime.57, 58, 59, 60 On the experimental side, there has been some initial work investigating singlet fission or related processes such as triplet-triplet annihilation in optical cavities or with plasmonic nanostructures.,,61, 62, 63 Most of them have focused on the long-time (nanoseconds to microseconds) dynamics., In a recent study, transient optical spectroscopy was employed to monitor the early dynamics of a TIPS-pentacene film placed in a Fabry-Perot resonator. The rate constants extracted from fitting the experimental data did not show significant differences between cavity and non-cavity situations. These results suggested that the dark-state manifold dominated the dynamics and thus collective strong coupling was unable to modify the singlet fission process.

In this work, we focus on singlet-exciton fission dynamics under collective strong coupling and explore how the dynamics is affected by the presence of dark states and state broadening due to both the natural linewidth of the vibronic peaks and energetic disorder. The central question we aim to answer is whether the state mixing induced by strong light-matter coupling can enhance the singlet fission rate in prototypical organic materials, and which materials are most suitable for such an application. Our results show that singlet fission indeed becomes faster when the is spectrally tuned to be on resonance with the multiexcitonic state. Importantly, this mechanism is also operative for non-conventional singlet fission materials that present a large singlet-multiexcitonic state gap. We also find that the enhancement mechanism is not significantly affected by dark-state-induced dephasing in these compounds. Moreover, when energetic state broadening is considered, the enhancement in the rate is more robust against disorder for strongly exothermic materials than for conventional ones. This combination of properties opens up a whole range of opportunities for materials that have not been explored for singlet fission to date. Throughout the manuscript, we use italic characters to denote diabatic electronic states, while eigenstates of the system are indicated in boldface and labeled according to their main diabatic contribution.

Results and discussion

Singlet fission dynamics in an optical cavity

To address the feasibility of singlet fission under strong light-matter coupling, we first treat a single-molecule model and then extend it to the many-molecule case in the following section. Our model describes several molecular electronic excitations, a bath of intramolecular vibrations, and a cavity photon. We fully consider the coupling between the electronic excitations and the photonic mode, and treat the electronic-vibrational interaction perturbatively by relying on a master equation approach based on Bloch-Redfield theory. There are two main reasons behind this choice. First, it allows us to naturally incorporate two essential ingredients that have a leading role in light-matter strong coupling: the effect of many molecules as well as the finite cavity lifetime (through an additional Lindblad term). Second, in contrast to the commonly employed Fermi-Golden rule approaches that rely on a perturbative treatment of the electronic couplings,64, 65, 66, 67, 68 Bloch-Redfield theory treats them exactly and can properly describe singlet fission dynamics when the coupling to the vibrational bath is not too strong.69, 70, 71 Note that the physics of the dark polaritons described in the Holstein-Tavis-Cummings model commonly used to describe molecular polaritons,, which treats a single vibrational mode explicitly, are also represented in our Bloch-Redfield simulations in a slightly different form. For instance, X-type polaritons, which are effectively dark states in the Holstein-Tavis-Cummings model, can be represented within the single-excitation subspace by considering energetic broadening of the states.

The Bloch-Redfield master equation for the time evolution of the density matrix describing the electronic and photonic degrees of freedom can be written as:where , and d indices run over the eigenstates of the system Hamiltonian, are the eigenfrequency differences, and is the Bloch-Redfield tensor accounting for the system-bath interaction; i.e., the coupling between electronic and vibrational modes. Relaxation rates are expressed in terms of the spectral density representing a bilinear electron-phonon interaction that we approximate by an Ohmic function with a Lorentz-Drude cutoff,, for which we have chosen characteristic parameters of representative singlet fission materials. Although the specifics of the spectral density of the environment can be relevant for quantitative predictions for a given molecular species, here we focus on providing insight and general guidelines on the circumstances under which collective strong coupling modifies the singlet fission dynamics, and thus use a general molecule-independent spectral density. Note that we have checked that our results are not strongly affected by the spectral density parameters; that is, the cutoff and reorganization energy.

The system Hamiltonian for the molecule interacting with the cavity mode is given bywith the electronic Hamiltonianwhere and are the energies and interstate couplings, respectively, for the (diabatic) electronic states involved in singlet fission. We use a four-state model to represent a system with the ground state (), the optically active singlet exciton (), the triplet-pair state (), and a charge-transfer () state, as suggested by Reichman and coworkers. Specifically, we consider first the case of a slightly exothermic formation of the state ( meV, eV) and a higher-lying state ( meV). These energetics are representative of efficient singlet fission compounds in which the state mediates the population transfer from the to the multiexcitonic state via a superexchange mechanism. In many singlet fission materials, the couplings between the and the and states are (at least) one order of magnitude larger than the direct interaction, since the former contain one-electron terms while the latter can be approximated as the difference between two bielectronic integrals. Here we take meV, whereas we disregard the direct coupling ( meV). Note that previous work on singlet fission under collective strong coupling focused on a direct mechanism where no intermediates are involved.

The cavity term in (Equation 3) takes the form , where is the frequency of the cavity mode, which we take in resonance with the optical exciton (), and and are the bosonic creation and annihilation operators, respectively. Finally, the interaction between the cavity photon and the electronic states can be expressed by the Jaynes-Cummings Hamiltonian for the case of a single singlet fission site ():with being the coupling strength between the cavity photon and the i-th electronic excited state, which is half the Rabi frequency (). Unless otherwise indicated, we consider strong coupling to the bright state, with meV. We precisely choose this value because it is roughly equal to the energy gap between the and states, thus placing the close to the eigenstate. On the other hand, we consider both and to be optically non-active transitions, hence, their coupling to the cavity vanishes. In this work we are interested in the linear-response regime, such that simulations can be restricted up to the single-excitation subspace.

We account for the finite lifetime of the cavity photon due to radiative and nonradiative decay by including a Lindblad term in our simulations. In this work, we have considered a 50-ps cavity lifetime. We have checked that the minimum cavity lifetime necessary for our results to hold is around 10 ps. Note that, for cavity photon frequencies around 2 eV, as those we consider here, this would correspond to a quality factor of . While such Q factors are beyond those used in the context of molecular strong coupling, we note that they can be achieved experimentally at optical frequencies., In particular, in Najer et al., a value of was reported.

We would like to stress that, when planning to manipulate molecular photophysics with cavities, it is important to be aware of this additional deactivation channel, which is absent in the non-cavity situation, and might compete with the intrinsic molecular processes. In the situation we explore in this work, the non-cavity singlet fission rate is faster than the cavity decay. Bearing in mind that formation in pentacene, one of the most efficient singlet fission materials, occurs within 80 fs, one should aim for cavity lifetimes beyond the femtosecond range in order to modify singlet fission dynamics via polariton formation.

The singlet fission dynamics within and outside the optical cavity is represented in Figure 1. For the non-cavity case (Figure 1A), the initially populated state relaxes to the state with a mean time of 5 ps, as obtained from the Bloch-Redfield rate from the adiabatic to the eigenstate,expressed in terms of the eigenstate expansion coefficients and , and the power spectrum , which characterizes the environment’s ability to absorb or release the energy required for the transition between the two eigenstates to happen at a given finite temperature (in our case, room temperature). According to this expression, the singlet fission rate will be non-zero under two conditions: (1) both states must have a common diabatic contribution, and (2) the environment of the molecular vibrations is able to meet the energetic requirements for the transition to occur. As shown in Figure 1A, the very small and contributions to the adiabatic state are sufficient to drive the population relaxation. Note that the time evolution of the photophysical reaction proceeds with no significant population of the state, which acts only as a mediator agent.

Figure 1

Energy spectra and population dynamics for N = 1

(A and B) Energy spectrum (left), initial eigenstate diabatic composition (middle), and population dynamics (right) for (A) bare case, and strong coupling with (B) ( meV). Note that the eigenstates are colored according to their main diabatic contribution. Color code: cavity photon (yellow), (dark green), (red), (dark blue), (lime green), (indigo), and (orange).

Next, we focus on how strong light-matter coupling can affect the singlet fission dynamics. In this case, we explore the time evolution of the system starting from the . Since the population transfer rate is analogous to that for the non-cavity case in Equation 6 by replacing with , then a straightforward strategy to enhance the singlet fission rate might be to increase the mixing between states. For that, we couple the cavity photon with the state such that the is placed close to resonance with the state. In this scenario, the acquires some character and the state acquires some polaritonic, i.e., and cavity photon, character, thus enhancing the population transfer rate. This is indeed what happens, as shown in Figure 1B, where the singlet fission dynamics is accelerated at earlier times compared with the non-cavity situation, with a considerably shorter mean time (0.4 ps). Note that state mixing was invoked in Polak et al. to qualitatively interpret the reverse process, triplet-triplet annihilation, in an optical cavity. Although the rate populating the state is modified under strong coupling, the -mediated singlet fission mechanism is maintained, as shown by the nearly null population of the state during the entire photophysical reaction, and the fact that the multiexcitonic state is not populated without the presence of the state. Notice also that, because of the finite cavity lifetime (50 ps), the ground state () becomes populated due to the non-zero cavity photon character of the and states. It is important to highlight that, since the is practically degenerate with the state, it is pure dephasing that promotes the population transfer; i.e., it is the bath spectral density evaluated at zero frequency that determines the Bloch-Redfield rate and, therefore, energy-dependent details in this spectral density are not too critical here. Finally, note that although we have not taken into account the cavity dispersion in our simulations, it could be detrimental to the polariton-accelerated singlet fission rate we observe because an additional relaxation channel, namely, intra-band scattering in the branch, would be present.

Collective strong light-matter coupling

In the following we extend our model to the case of N equivalent (non-interacting) singlet fission sites, that is, multiple electronic systems, interacting with a cavity photon. For that, we expand the electronic terms of the system Hamiltonian following the so-called Tavis-Cummings model:,where corresponds to an electronic excitation ( or ) at the k-th singlet fission site, and the light-matter coupling strength is the same for all k sites and is chosen so that the collective Rabi splitting stays constant (). In the following we also restrict the system to one excitation at most. Note that we do not consider permanent dipole moments in our model since we are interested in the collective regime where these are not enhanced by a factor, as opposed to the transition dipole moment. In the following, we show results for . Note that the relevant parameter here is the collective Rabi splitting and that we have checked that our results are converged for this value.

The most important difference between the single singlet fission site previously discussed and the collective case is the presence of the dark singlet-exciton manifold; i.e., linear combinations of the molecular states that do not couple to the cavity photon. These states, together with the two polaritons and the and manifolds, constitute the eigenstates of the collective system. As shown in Figure 2, the fast relaxation is maintained, and because of detailed balance, the saturation population increases since for larger N the many states outnumber the and act as a population sink. Therefore, this polariton-enhanced mechanism does not suffer from the issues arising in polariton-assisted TADF, previously discussed, in which a single-molecule excitation has to transfer to a collective polariton, incurring a penalty in the transition rate.

Figure 2

Energy spectra and population dynamics for N = 20

Energy spectrum (left), initial eigenstate diabatic composition (middle), and population dynamics (right) for ( meV). Note that the eigenstates are colored according to their main diabatic contribution. Color code: cavity photon (yellow), (dark green), (red), (dark blue), (lime green), (indigo), (orange), and dark states (royal blue).

It is straightforward to understand why the rapid initial relaxation from the to the manifold still holds in the collective case. The collective strong coupling version of Equation 6 is given bywhere labels the specific state of the manifold. These states are basically linear combinations of the non-cavity eigenstate on each site, , therefore, the rate constant can be approximated as . Since the amplitude of the state on each k-site is equal to the amplitude scaled by , i.e., , and since , then . Therefore, in the macroscopic limit, the rate from the to the set of N states is independent of N and is dictated by the amount of character the acquires in the case (). According to this analysis, collective strong coupling will affect the early singlet fission dynamics when the is close to resonance with the manifold, such that . Notice that, in the collective case, the ground state is barely populated, in contrast to the single-site results (Figure 1B). This is because only one eigenstate of the manifold has a non-vanishing cavity contribution.

Promoting strongly exothermic singlet fission

An important factor at play here that we have not addressed yet is the role of the {} dark states. In the example we have discussed, since this manifold is relatively close in energy to the , it is significantly populated during the dynamics, as shown in Figure 2. From the dark states, singlet fission then proceeds with essentially the same rate as in the non-cavity situation, then diminishing the rate enhancement due to strong coupling. However, the further these {} dark states are spectrally separated from the and the manifold, the less populated they will be after excitation to the and, as a consequence, they will be less detrimental to the cavity-based singlet fission rate enhancement. Therefore, we hypothesize that the ideal candidates for singlet fission under collective strong coupling are those compounds with a large enough gap, such that, when the is excited, the {} dark states are energetically too high to be populated. Standard singlet fission materials exhibit an state that ideally lies slightly above the multiexcitonic singlet state, i.e., weakly exothermic singlet fission, which severely limits the pool of potentially efficient singlet fission compounds. Also, prototypical efficient singlet fission materials, such as pentacene and rubrene, have poor photochemical stabilities,78, 79, 80 preventing their use in real-world devices. Therefore, according to our hypothesis, strong coupling could enhance the singlet fission rate in materials that have usually been ignored since excessive exoergicity is known to be detrimental in the non-cavity situation.,,

To test our hypothesis, in Figure 3 we plot the population dynamics dependence on the diabatic gap while keeping the gap constant. For the cavity case, as the gap varies, the Rabi splitting is chosen such that the lies close to resonance with the eigenstates ( meV above), just like in Figures 1B and 2. For the non-cavity case, the singlet fission dynamics is greatly slowed down as the separation between the two states increases as the amplitude contribution in Equation 6 decreases due to reduced state mixing. In contrast, within the cavity, the state is populated much faster and independently of the gap, as long as the is brought close to resonance with the eigenstates. These results therefore indicate that strong light-matter coupling can also accelerate the singlet fission dynamics of materials with large gaps.

Figure 3

Adiabatic population for the non-cavity ( initial state) and cavity ( initial state and ) situations

The single-molecule coupling constant meV as the gap increases since a larger Rabi splitting is needed for the to be close to resonance with the eigenstates.

Therefore, within the strong coupling regime, it becomes possible to efficiently generate triplet-pair states not only for slightly exothermic energetics; e.g., in pentacene. Relaxing the near-degeneracy criterion () should allow fast singlet fission processes to be obtained in a wider range of molecular systems with respect to those identified to date. In order to exemplify the potential impact in the search for singlet fission chromophores, in the following, we consider a set of organic molecules with suitable properties: (1) exothermic singlet fission (); (2) non-vanishing transition dipole moment between ground and lowest excited singlet, so strong coupling to the cavity mode could be achieved; and (3) sufficiently large ( eV) -to- gap in order to ensure molecular stability for practical applications. Our molecular test set contains 262 organic molecules obtained from the Cambridge Structural Database (CSD) selected following the protocol designed by Padula and collaborators. In order for singlet fission to compete with other relaxation pathways, such as internal conversion, intersystem crossing, or radiationless decay to the ground state, it should take place on a picosecond (or even sub-picosecond) timescale. Therefore, here we consider that good candidates for efficient singlet fission are those able to reach a 50% population of the state within 5 ps after photoexcitation. Of course, this is a somewhat arbitrary limit, but it reasonably serves our purpose to compare singlet fission efficiency with and without strong coupling. For the non-cavity scenario, it occurs for eV (Figure 3), which is in the order of the gap measured in crystalline pentacene,, and slightly higher than the value employed in the previous sections (0.08 eV). About 26% of the considered chromophores exhibit eV (green bars in Figure 4) and are expected to be the most promising singlet fission compounds. In contrast, if we consider eV as the upper limit for the Rabi splitting,, the number of potential molecules able to undergo singlet fission efficiently when placed in an optical cavity increases to 66% of the total set of studied molecules (orange bars in Figure 4). These results demonstrate how strong coupling can substantially increase the pool of suitable molecular candidates, which represents one of the main current challenges for the practical implementation of singlet fission in optoelectronic devices. Also note that the values reported here can be seen as rather conservative, since in the cavity model the energy difference remains constant while varying the gap, debilitating the efficiency of the CT-mediated mechanism (for large gaps, the CT states lie rather high with respect to and ). Moreover, we have not considered here the case of direct singlet fission, i.e., direct coupling, which would also benefit from the energy-level alignment between and states.

Figure 4

Histogram (in eV) of a list of 262 organic molecules

Experimental structures from Montalti et al. Singlet and triplet energies obtained from Padula et al.

Effect of state broadening

So far we have assumed equivalent singlet fission sites with well-defined discrete , , and energies. Since our conclusions rely on the being close to resonance with the manifold, it is not clear whether they will still hold when taking into account the energetic broadening of the states. In the following, we account for the linewidth of the Franck-Condon vibronic transition that couples to the cavity mode and also inhomogeneous broadening due to energetic disorder because of different local environments. To explore this effect, we select the molecular energies by sampling a Gaussian distribution centered at with . We consider a maximum full-width half-maximum (FWHM) value of 0.1 eV, which is characteristic of the vibronic absorption peaks of polycyclic aromatic hydrocarbons. We explore two limiting cases, one in which the energies of the three molecular states of each singlet fission site are jointly sampled, that is, the same random number is used to sample the Gaussian distribution for all three states of a given site, and another one for which the three states are sampled independently by using different random numbers. Note that, in the former case, the energy gaps between states remain constant, while in the latter they may vary.

In Figure 5, we plot the strong coupling results (red) and compare them with the non-cavity ones (blue) with (light) and without (dark) state broadening. Figures 5A and 5B show the results for a system with an gap of 0.2 eV, representative of the class of materials we propose for cavity-mediated singlet fission. Like in the previous sections, the coupling strength is chosen such that the is close to resonance with the eigenstates ( meV). In general, the singlet fission dynamics in the cavity is slowed down once state broadening is considered (light versus dark), but it is still much faster than the non-cavity situation (red versus blue). Therefore, our proposal is still operative in a realistic situation where the states have a finite width, and strong coupling can enhance singlet fission dynamics of materials that present a large gap even in the presence of state broadening.

Figure 5

(diabatic) population dynamics with (light) and without (dark) state broadening

(A–D) Non-cavity results are plotted in blue and strong coupling ones in red. The shaded area includes the mean SD of 50 realizations. Results are shown for the case. Coupling strength in (A) and (B) meV, and in (C) and (D) meV. The initial state for the cavity simulations is the . In the presence of state broadening, it is taken to be , where is the cavity photon coefficient of the r-th eigenstate , and only those eigenstates below the cavity photon are considered. For the non-cavity simulations, the initial state is taken as the eigenstate with the largest amplitude.

We explore next the situation in which the gap is smaller (0.08 eV) and, accordingly, the Rabi splitting is also reduced (0.15 eV) to match the spectral location with that of the eigenstate. Our results are shown in Figures 5C and 5D for the two sampling cases. Note first that, when comparing the dark lines (the same in both panels), in which state broadening is not taken into account, the cavity and non-cavity dynamics (red versus blue) are more alike than for the case with a larger gap (Figures 5A and 5B). As discussed before, this is due to the rapid relaxation to the dark-state manifold when the gap is small. Moreover, when state broadening is incorporated in the simulations, cavity dynamics (light red line) turns out to be very similar to the non-cavity situation (light blue line) and the shaded areas even overlap. These results provide insight as to why, in a recent experimental work, no significant differences were found when comparing the cavity and non-cavity singlet fission dynamics of TIPS-pentacene. In the experimental setup, both the Rabi splitting and the FWHM of the vibronic peak that is coupled to the cavity mode were eV. Notably, the gap in TIPS-pentacene is also expected to be close to this value. Therefore, the experimental conditions resemble those illustrated in Figures 5C and 5D, which predict very moderate changes in singlet fission dynamics. This further confirms that singlet fission rates in slightly exothermic systems are not expected to be largely boosted by strong coupling. We also note that, although the cavity lifetime of the experimental work was much shorter that the one we consider here ( 20 fs versus 50 ps), the fact that singlet fission dynamics inside and outside the cavity took place at practically the same rate is an indication that, in the cavity experiments, the polaritons populated the dark-state manifold before decaying due to their photonic contribution. Therefore, the cavity lifetime is not a relevant parameter since it did not play an active role.

Finally, regarding the 50-ps cavity lifetime we consider in this work, although a long cavity lifetime is associated with a narrow linewidth, under strong coupling with organic molecules, when the polariton overlaps with (the tail of) the broad absorption band of the molecules, its effective linewidth will be increased, while decay of the excitation will not. We thus expect that even very long-lived cavity modes will lead to relatively broadband polaritons that can capture a sufficient fraction of the incoming light to be potentially useful for singlet fission photovoltaic applications.

To conclude, we have demonstrated that singlet fission dynamics can be accelerated under collective strong coupling via state mixing when the is almost resonant with the state. We have also shown that this effect is much more beneficial for compounds with large gaps, thus reducing both the population transfer from the to the dark-state manifold and the detrimental impact that energetic broadening has on the resonant mechanism behind the rate acceleration. Given that the main characteristic of conventional singlet fission materials is a small gap, and that systems such as polycyclic aromatic hydrocarbons are unstable and pose significant challenges for practical applications, our results provide a new perspective and suggest a new paradigm for cavity-mediated singlet fission with materials that have not previously been considered for singlet fission. Our results also highlight the fact that cavities with long-lived photons beyond a picosecond lifetime are required for singlet fission to benefit from strong light-matter coupling. We hope that our results can serve as a guide and inspire future experiments to realize efficient singlet fission with unconventional compounds in optical cavities.

Experimental procedures

Resource availability

Lead contact

Further information and requests for resources should be directed to and will be fulfilled by the lead contact, Clàudia Climent (ccliment@sas.upenn.edu).

Materials availability

This study did not generate new unique materials.

Methods

In our simulations, each molecular excited state is coupled to identical and independent baths. We consider Ohmic spectral densities with a Lorentizian cutoff , where is the cutoff frequency and λ the reorganization energy. Following Berkelbach et al., we take the values meV and meV, which are characteristic of organic aromatic molecules.

The power spectrum is given bywith the Bose occupation factor , and we take K.

To account for the finitie cavity lifetime, we include a Lindblad term in our simulations, , where , where is the bosonic destruction operator of the cavity photon, and κ is the lifetime.

Note that we do not rely on the secular approximation of the Bloch-Redfield equation. This is because this approximation fails when there are eigenstates close in energy. The master equation has been solved with the Qutip package.,