Magnetic dipolar interaction between correlated triplets created by singlet fission in tetracene crystals.

Singlet fission can potentially break the Shockley-Queisser efficiency limit in single-junction solar cells by splitting one photoexcited singlet exciton (S1) into two triplets (2T1) in organic semiconductors. A dark multiexciton state has been proposed as the intermediate connecting S1 to 2T1. However, the exact nature of this multiexciton state, especially how the doubly excited triplets interact, remains elusive. Here we report a quantitative study on the magnetic dipolar interaction between singlet-fission-induced correlated triplets in tetracene crystals by monitoring quantum beats relevant to the multiexciton sublevels at room temperature. The resonances of multiexciton sublevels approached by tuning an external magnetic field are observed to be avoided, which agrees well with the theoretical predictions considering a magnetic dipolar interaction of ∼ 0.008 GHz. Our work quantifies the magnetic dipolar interaction in certain organic materials and marks an important step towards understanding the underlying physics of the multiexciton state in singlet fission.

The spin-allowed singlet fission (SF) process is highly efficient in some organic semiconductors123456789. By creating two triplet excitons from a single photoexcited singlet exciton, the SF process of carrier multiplication can be implanted in multiple device architectures to improve the efficiency of solar conversion1011121314. However, the intrinsic mechanism responsible for the fast fission process is still under intense debate, despite remarkable progresses have recently been made23456781516171819202122. The essential role played by the intermediate multiexciton (ME) state (1(TT)) has been realized12345151720 with the SF process being described as1,

The doubly excited ME state, also referred to as the correlated pair of triplet excitons, which are entangled with spin coherence, is created when a singlet-excited molecule shares its energy with a neighbouring molecule at its ground state. The magnetic interaction between the correlated triplets is insightful for revealing the nature of the ME state by providing valuable information about their spatial separation and the effect of spin coherence172023. However, it is challenging to quantify the weak magnetic dipolar interaction between the correlated triplet excitons at the ME state with a short lifetime1820.

Magnetic dipolar interaction is pivotal for many technically significant processes in organic materials including exciton fission/fusion1, organic magneto-resistance24 and organic photovoltaics25. In a SF sensitizer, a conceivable approach to quantify the magnetic dipolar interaction is to investigate the interaction-induced energy shift of the ME state. In crystalline tetracene, the energy differences between the ME sublevels can be monitored through quantum beating signals in the singlet population172026. As depicted in Fig. 1a, exciton fusion from two ME sublevels induces quantum beats in population of the S1 state, manifesting themselves as an oscillation in the time-resolved fluorescence (TRFL) spectrum. After first introduced in 1980s26, this quantum beat phenomenon in tetracene has recently been re-examined and comprehensively explained by Bardeen's group172027. Being a direct evidence of the ME state, the quantum beating signal has been regarded as a fingerprint of SF in tetracene151720. However, the signature of magnetic dipolar interaction has never been directly captured in tetracene2026, and only an upper bound of ∼0.2 GHz was given due to the limitations in previous experiment20.

Figure 1

Interaction-induced anti-crossing among ME sublevels.

(a) The scenario of the quantum beats related to the ME sublevels in tetracene crystals at zero field. Exciton fusion from sublevels of ME state induces the delayed fluorescence where quantum beats manifest themselves as oscillations in the TRFL spectrum. (b,c) The field dependences of two near-resonant sublevels for the cases of non-interaction and with interaction between correlated triplets, respectively. With magnetic interaction, the emergence of mixed states results in an avoided crossing with separated upper and lower branches (c).

Here we exploit a new scenario of the interaction-induced anti-crossing to quantify the interaction strength in crystalline tetracene. We manipulate the ME sublevels by applying an external magnetic field to approach the level-crossing resonance. The magnetic dipolar interaction between the correlated triplets causes an avoided level crossing, resulting in two branches of levels with an energy gap. The interaction strength is proportional to the gap size that can be precisely evaluated by measuring the frequencies of quantum beats relevant to the resultant upper and lower levels with TRFL spectroscopy. With this approach, we have succeeded in quantitatively determining the magnetic dipolar interaction between the correlated triplet excitons at the ME state to be ∼0.008 GHz. The unexpected weak interaction strength indicates that it is necessary to consider the exciton size for a comprehensive understanding of the SF process.

Results

Theoretical consideration

Figure 1 schematically shows the scenario of the quantum beats relevant to the ME sublevels involving in the SF process. We describe the sublevels of ME state with a spin-dependent Hamiltonian of two correlated triplets. The Hamiltonian includes two isolated triplets (for example, α and β) and their mutual interaction (Hint)20, that is,

The isolated triplet dipole (for example, α or β, as represented by ‘Tri') can be described by the Hamiltonian as2028

Here the first item represents the Zeeman shift due to the applied magnetic field, where S, B, μ and g are the spin operator, the external magnetic field, the Bohr magneton and the Lande g-factor, respectively. The latter two terms give zero-field Hamiltonian with parameters of D* and E* characterized by ESR experiments29. The Hamiltonian for the magnetic dipolar interaction between the correlated triplet dipoles can be written as24

where S and S are the spin operators of the two triplets, R is the displacement vector between the two triplets with the magnitude of R0 and X is the parameter of interaction strength. As described in Supplementary Note 1, we calculate the energies of the ME sublevels by solving the Hamiltonian of equation (2) with a basis of nine eigenstates . In spite of weak magnitude, the interaction part can cause substantial differences from the energy levels predicted by the interaction-free model (Supplementary Figs 1 and 2). Those sublevels with non-zero mapping to the S1 state involve in the exciton fission/fusion processes and directly contribute to the quantum beats (Fig. 1a)20.

Theoretical calculation predicts a sublevel crossing with a magnetic field of ∼420 G applied along the x axis (Supplementary Fig. 1). The mutual interaction opens a gap of Δ≈2X at the level-crossing resonance (Fig. 1c, Supplementary Fig. 1b,d). Since the gap size is small, it is difficult to directly measure the beat frequencies in time domain due to damping of the oscillations20. To overcome this obstacle, we slightly tilt the external magnetic field with a small angle (θ) relative to the x axis in the xy plane (Fig. 2a) to introduce an extra perturbation (δEy) in the Hamiltonian. The gap size then increases to be Δ≈2X+2δEy (Supplementary Note 2), so the evaluation becomes more accurate because the visibility of quantum beat is significantly enhanced in the TRFL trace.

Figure 2

Magnetic-field-dependent quantum beats.

(a) Schematic diagram of the external magnetic field applied at the xy plane with a small angle (θ) tilted relative to the x axis. (b) The decay curves of fluorescence dynamics with three different external magnetic field value tilted at θ=2° near the resonance of level crossing. Inset shows the oscillation part obtained by subtracting the multi-exponential decay components from the raw data recorded at 420 G. (c) The relative amplitudes of quantum beats are plotted as functions of the beat frequency and the field magnitude. The dashed lines indicate the expected beat frequencies calculated with interaction-free model.

Level-crossing resonance

Figure 2b plots the TRFL traces recorded under different magnetic fields at θ∼2°. An anomalous oscillation with the frequency of ∼0.05 GHz emerges due to the newly opened gap when the field is at 420 G (inset, Fig. 2b). The field-dependent amplitudes and frequencies of the measured quantum beats (Fig. 2c) are compared with theoretical calculations (Supplementary Fig. 3) to further confirm the above assignments. The interaction-free model can well explain the observed high-frequency beats (>0.1 GHz) (Supplementary Fig. 3), but fails to account for the anomalous part (∼420 G) without considering the opened energy gap (Δ). A control experiment has been performed with the magnetic field applied at z direction (Supplementary Note 3). Under this configuration, the anomalous beat no longer exists since no level crossing is found (Supplementary Figs 6 and 7), verifying that the observed low-frequency quantum beat is indeed a result of gap opened by the magnetic dipolar interaction and perturbation.

Next, we analyse the θ-dependent oscillatory components to quantify the interaction strength (Fig. 3a). We assume the displacement vector (R/R0) to be along the nearest neighbour orientation316. The resonant energies of all nine ME sublevels are calculated by exactly diagonalizing Hamiltonian (2) and indexed by their zero-field energies. When the magnetic dipolar interaction is included (Supplementary Fig. 1), the avoided crossing occurs between the three and four sublevels (Fig. 3c) as well as the sixth, seventh and ninth sublevels (Fig. 3d) at ∼420 G. The mapping products of for selected sublevel pairs, which directly reflect the beating amplitudes, are compared in Fig. 3e,f, suggesting the dominant contribution from the third and fourth sublevels (Supplementary Note 2; Supplementary Fig. 5). In Fig. 3b, the measured peak oscillation frequency of this anomalous beating signal as a function of the tilt angle θ is compared with the theoretical curves of energy separation between three and four sublevels calculated with different interaction strengths. The perturbation induced by field tilting corresponds to the zero-interaction curve (X=0) in Fig. 3b with δEy∝θ (Supplementary Fig. 4). The disparity of experimental data from the linear dependence on the tilt angle near θ≈0 is an evidence of existing magnetic dipolar interaction, which can be best reproduced by calculation with the interaction strength of X=0.008 GHz (Fig. 3b).

Figure 3

Quantum beats near level-crossing resonances.

(a) The oscillation components are plotted versus the delay time and tilt angle (θ). The external field is set at 420 G. (b) Experimental results of θ-dependent beat frequencies are compared with theoretical calculated curves considering different strengths (X) of magnetic dipolar interaction. (c–f) Theoretical considerations of the interaction-induced anti-crossing of various ME sublevels (for details, see Supplementary Note 1). At the resonance field of ∼420 G, the anti-crossing occurs in two cases, resulting in the mixed states of three and four sublevels (c) and the mixed states of sixth, seventh, and ninth sublevels (d). The field-dependent amplitudes of low-frequency quantum beats calculated approximately as the mapping products of selected sublevel pairs for these two cases are plotted in e and f, respectively.

Strong magnetic field limit

Furthermore, we have also examined the quantum beats at the strong magnetic field limit. A field of ∼3,000 G is applied in the xz plane with an angle of Φ to the x axis as shown in Fig. 4a. In this regime, multiple beat frequencies collapse into one frequency, because only two sublevels dominate the mappings to the S1 state (Supplementary Fig. 2; Supplementary Note 4)27. The energy alignments of the two sublevels calculated using models with and without the interaction are shown in Fig. 4c,b, respectively. The two sublevels approach to degeneracy at Φ≈69° (Fig. 4b), which is avoided due to the magnetic dipolar interaction (Fig. 4c). The oscillatory amplitude is plotted as a function of delay time and the angle deviated from the degeneracy value (ΔΦ) in Fig. 4d. The measured ΔΦ-dependent beat frequencies shown in Fig. 4e agree well with the calculated curve having X=0.008 GHz. The value agrees well with that obtained at 420 G, which is significantly larger than the bounds of errors and spurious effects including the hyperfine coupling as discussed in Supplementary Methods (Supplementary Table 3). The strength measured here is at the same order of magnitude as that adopted to interpret the data of optically detected magnetic resonance of germinate triplet pair at the T1+T1 state in bis(triisopropylsilylethynyl)-tetracene films23. The weak magnitude measured here is also in consistence with the upper bound estimated in literature20 and theoretical analysis of other magnetic effects30. It is over two orders of magnitude smaller than the energy difference between the sublevels at zero field, which might explain why the signature of this magnetic dipolar interaction was absent in previous studies2026.

Figure 4

Quantum beats at the strong magnetic field limit.

(a) Schematic diagram of a strong field applied at the xz plane having an angle Φ with respect to the x axis. In this strong field limit, only two dominant sublevels are involved in the fission and fusion processes. (b,c) The calculated energy alignments of two sublevels for the cases of non-interaction and with interaction between correlated triplets, respectively. The level crossing occurs at Φ≈69°. (d) The oscillation components are plotted versus the delay time and the angle (ΔΦ) deviated from the resonant value. The magnitude of the external field is ∼3,000 G. (e) The experimentally measured beating frequencies are plotted as a function of ΔΦ under the strong field limit. The curves are theoretically calculated results with different strengths of magnetic interaction.

Discussion

The orientation of displacement vector (R/R0) with respect to the z axis of magnetic tensor may affect the evaluation of the interaction strength. We have considered all four possible configurations for two neighbouring molecules in ab plane (Supplementary Fig. 8). The interaction strength X that best reproduces the experimental data (Figs 3f and 4d) is far below the theoretical values for two nearest neighbour molecules (Supplementary Tables 1 and 2). The separation distance between two triplets at the ME state equivalent to the measured value of X can be estimated to be ∼1.7 nm, which is about two to four times of the intermolecular distance (Supplementary Note 5). The unexpected large separation between the correlated triplets can be plausibly caused by the delocalization effect of photoexcited singlet excitons31. Previous PL study has suggested the size of singlet excitons to be in the order of tens of molecules at low temperature31, which might also affect the exciton dynamics in certain extent at room temperature7. The pair of triplets resulted from the delocalized singlet excitons can distribute over a distance of multiple molecules, which can naturally explain the measured weak magnetic interaction.

Moreover, though the diffusion of triplet excitons at the ME state remains poorly understood, triplet diffusion at the T1+T1 state has been shown to be quite efficient in tetracene crystals32. In principle, the triplet diffusion may increase the separation distance between the triplet excitons over time and cause a time-dependent X. Such effect has not been clearly seen in this study, which could be ascribed as the consequence of quantum beat spectroscopy used in this study. As depicted in Fig. 1a, only the correlated triplet pairs that recombine through the exciton fusion process contribute significantly to the beating signal. The triplet pairs that diffuse apart may lose the spin coherence and/or the coherence needed for exciton fusion. Those triplet pairs generate delayed fluorescence but do not make significant contribution to the beating signal. Consequently, the time dependence of interaction strength is not manifested in the quantum beat spectroscopy here. Additional work will be needed to fully uncover the effect of exciton diffusion on the interaction between the triplet excitons at the ME state.

Overall, the spins of correlated triplets are weakly coupled with the environment, making our method reliable in evaluating the magnetic dipolar interaction at room temperature. This work represents the first successful quantification of the magnetic dipolar interaction between the correlated triplet excitons at the ME state in crystalline tetracene. The weak interaction strength observed here provides a direct evidence that the initial step (S0+S1↔1(TT)) of SF process involves molecules that are not nearest to each other. The entangled triplets at the ME state delocalized over many molecules arising from the delocalization of initial singlet excitons could be a plausible origin of the weak interaction strength. Such phenomena indicate that an alternative scenario different from the available models considering two proximately coupled molecules may be necessary to fully understand the SF process in crystalline tetracene. Including the effect of exciton delocalization to elucidate SF may be enlightening to address the debates between a charge-transfer mechanism and coherent superposition model responsible for the generation of intermediate ME states2161920333435. The quantum beating behaviour in crystalline tetracene due to spin coherence of the correlated triplet pairs is sensitive to the magnitude and direction of the magnetic field. These features deserve more in-depth investigations and can stimulate further adventures in applying SF to many new areas of research such as spintronics and quantum information science363738.

Methods

Sample preparation and magnetic field alignment

Tetracene single crystals with the thickness of ∼1 μm and size up to 5 × 5 mm2 were prepared by the method of physical vapour deposition. The crystallographic axes were determined by X-ray diffraction and polarization microscopy (Supplementary Fig. 9), which were further employed to determine the magnetic axes with a transform matrix as established by electron spin resonance measurements (Supplementary Methods). A rotatable magnetic coil was used for magnetic-field-dependent experiments. The samples were mounted on a multi-axis platform to realize the desired alignment with respect to the magnetic field.

Optical characterizations

The second harmonic field (400 nm) of the pulses emitted from a Ti:Sapphire oscillator (800 nm, Vitara, Coherent) was chosen as the excitation source. The repetition rate was reduced down to 4 MHz by a pulse picker and the excitation flux was kept at a low level (∼4 nJ cm−2) to avoid the effect of exciton–exciton annihilation. The emission light was collected by a fibre and routed to a spectrograph. The TRFL spectra at 530 nm were recorded with the technique of time-correlated single-photon counting by an avalanche photodiode having a temporal resolution of better than 50 ps. The multi-exponential decay components were subtracted from TRFL traces before Fourier transform to extract the frequencies and amplitudes of quantum beating signals (Supplementary Fig. 10).

Theoretical simulation

The resonant energies of ME sublevels are calculated by exact diagonalization of the full spin-dependent Hamiltonian. Quantum beat in the detected TRFL signal is modelled by the dynamics of the density matrix elements governed by the quantum Liouville equations (Supplementary Note 1).

Additional information

How to cite this article: Wang, R. et al. Magnetic dipolar interaction between correlated triplets created by singlet fission in tetracene crystals. Nat. Commun. 6:8602 doi: 10.1038/ncomms9602 (2015).