Singlet fission for quantum information and quantum computing: the parallel JDE model

Singlet fission is a photoconversion process that generates a doubly excited, maximally spin entangled pair state. This state has applications to quantum information and computing that are only beginning to be realized. In this article, we construct and analyze a spin-exciton hamiltonian to describe the dynamics of the two-triplet state. We find the selection rules that connect the doubly excited, spin-singlet state to the manifold of quintet states and comment on the mechanism and conditions for the transition into formally independent triplets. For adjacent dimers that are oriented and immobilized in an inert host, singlet fission can be strongly state-selective. We make predictions for electron paramagnetic resonance experiments and analyze experimental data from recent literature. Our results give conditions for which magnetic resonance pulses can drive transitions between optically polarized magnetic sublevels of the two-exciton states, making it possible to realize quantum gates at room temperature in these systems.

Materials capable of storing and manipulating quantum data must maintain quantum coherences and entanglement over timescales that are orders of magnitude longer than the system's quantum beat period 1 . But quantum states are fragile, and most materials do not sustain quantum coherences when temperatures are in excess of a few Kelvin. This "tyranny of low temperature" is a major hurdle to realizing accessible quantum computing and information technologies.
From a quantum information perspective, the photoproducts of singlet fission can solve two outstanding problems associated with the tyranny of low temperature, and they do it in complementary ways. Each corresponds to one of the criteria for quantum computing set forth by Ref. 1 . To paraphrase: a system capable of quantum information processing must be prepared in a pure quantum state, not a mixed one 2 . Once initialized, one must be able to execute a deterministic sequence of unitary operations, or quantum gates, on that state so that it can be coaxed into collapsing on the final state, which is the solution to a computational problem, with high probability and in polynomial time 3,4 .
Because magnetic resonance experiments operate in the "strong-field" regime where the Rabi frequency, , is comparable to the transition frequency, ω (Fig. 1a), they can drive arbitrary unitary operations between quantum spin states 5 . But the gap between the ground and excited states in these experiments is small relative to k B T , and as a result, there is a great deal of thermally generated uncertainty in the initial state of the system. In the language of the density matrix, the initial state is mixed, not pure. Optical experiments, by contrast, have a large gap between ground and excited states relative to k B T (Fig. 1b), so that thermal fluctuations do not generate appreciable uncertainty in the ensemble of initial states. The Rabi frequency in optical experiments, however, is perturbatively small relative to the transition frequency. It is therefore very difficult or impossible to induce a population inversion in an optical experiment, which severely limits the ability of a purely optical experiment to perform quantum gate operations.
In more recent years, researchers have become interested in systems where it is possible to generate spin polarization through optical pumping (Fig. 1c) 6,7 . The optical field removes the uncertainty in the initial state of the system, and magnetic resonance experiments on the optically prepared photoproduct, well in the strong-field regime, perform the quantum gate operations. The nitrogen-vacancy (NV) center in diamond is an example of a system that operates on these principles. In solid-state systems, like the NV center, spin centers are implanted into the material post-synthesis. Because the defects are often randomly dispersed, these materials have problems with scalability that might be overcome in molecular systems that are synthesized from the bottom up.
In this article, we explore the phenomenon of singlet fission as a novel platform on which one might build quantum data structures and gates near room temperature. Singlet fission is a photophysical interconversion process that takes place between specifically designed organic chromophores where, following excitation with a photon ( γ ), an optically bright singlet state, S 0 S 1 , rapidly relaxes into a doubly excited, spin-singlet state, 1 TT Eqs. (1-2) 8-10 , In this notation, the superscript indicates the state's multiplicity in terms of the total spin S, and the TT designates the spatial nature of the spin wavefunction. The | 1 TT� spin wavefunction is a maximally entangled, coherent superposition over three of the nine spin sublevels belonging to the two triplets on the chromophores [11][12][13] . In the last few years, electron paramagnetic resonance (EPR) experiments have combined pulsed optical laser excitation with magnetic resonance to reveal that the 1 TT state evolves into various 2S+1 TT M exciton states 14-22 , where the subscript, M, refers to the magnetic sublevel −S ≤ M ≤ +S . Although singlet fission is a spin-conserving process (Eq. 2), the triplet pair states (Eq. 3) are not eigenfunctions of the electronic spin hamiltonian; they are non-stationary and evolve in time.
In crystalline systems, these excitons may hop to neighboring sites, becoming increasingly more distant, and eventually unpair into T+T 23 : There has not been a consistent microscopic theory that can explain the set of relaxation phenomena embodied in Eqs. (1)(2)(3)(4).
While the lion's share of attention over the past decade has focused on maximizing the conversion of the 1 TT state to T+T for solar energy applications 24,25 , we argue that for quantum information applications, one should seek to instead limit the decay into independent excitons, by designing molecules that make the conversion from 1 TT to 2S+1 TT M as state-selective as possible.
In this work we consider singlet fission between chromophores in two classes of systems that are widespread in the literature: covalently linked organic dimer molecules 26 and organic crystals comprised of chromophore pairs 14 . As is often the case, the selection rules governing the quantum relaxation phenomena depend sensitively on molecular symmetries. In particular, a pair of identical chromophores, where one molecule is related to the other by a translation (Fig. 2b), will also exhibit a permutation symmetry for the exciton triplet pair. As we show, that symmetry isolates the 3 TT triplet states so that spin relaxation only proceeds between 1 TT and 5 TT states. The most state-selective relaxation occurs between the 1 TT singlet state and the 5 TT 0 quintet state in an ordered and immobilized system of molecular dimers that have their magnetic principal axes mutually parallel to one another.

Results
The spin-exciton hamiltonian. We follow with a derivation of the spin-exciton hamiltonian for the triplet pairs in Eq. (3), exploiting approximations for the light atom molecules characteristic of singlet fission chromophores. The singlet fission process depicted in Eq. (2) is often much faster (picosecond or sub-picosecond) than the timescales on which the states in Eq. (4) evolve (nanoseconds to microseconds). The spin-orbit interaction is small and is often ignored, but it is straightforward to include 27 . Optically induced spin polarization generates states for quantum information applications near room temperature. (a) In magnetic resonance experiments, the Rabi frequency, , can be comparable to the transition (Larmor) frequency, ω . In this strong field limit, it is possible to completely manipulate a quantum state (qubit). However, the energy gap between states is small relative to k B T . This is a source of uncertainty when the state is initialized. (b) In optical transitions, the energy gap is large relative to k B T . But the Rabi frequency is much less than the transition frequency, which means that gate operations done with weak optical fields will be incomplete and noisy. (c) By coupling optical excitations to an internal conversion process, such as singlet fission (wavy arrow), one may capitalize on the advantages of both methods, provided that the relaxation is state-selective.
Scientific Reports | (2020) 10:18480 | https://doi.org/10.1038/s41598-020-75459-x www.nature.com/scientificreports/ After making the Born-Oppenheimer approximation, suppressing orbital degrees of freedom and ignoring hyperfine interactions, we begin with a general hamiltonian that is bilinear in all spin-spin interactions, where s i is the electron spin operator for orbital i, and the indices i and j enumerate the HOMO and LUMO levels in a frontier molecular orbital description of the chromophore pair (Fig. 2a). Similarly, the g-tensors in singlet fission chromophores are often isotropic 8 , though anisotropy in the g-tensors can also be included perturbatively. O ij is a rank-two tensor that accounts for the spin-spin interactions between the four electrons in the four orbitals. Within this framework, O ij depends implicitly on integrals over spatial wavefunctions. Akin to electronically nonadiabatic effects in molecules, fluctuations in nuclear coordinates and exciton hopping can make the O ij time-dependent 23,28 . These time-dependent O ij parameters can drive transitions between different spin states of the spin-exciton hamiltonian.
The interaction tensor O ij for each electron pair can be decomposed into three separate terms: a scalar "isotropic" part, an antisymmetric tensor of rank one, and a traceless, anisotropic tensor of rank two. The isotropic term yields the usual Dirac-Heisenberg exchange coupling 29 . Unlike the isotropic exchange interaction, both the antisymmetric and the anisotropic terms are formally relativistic in nature, and in light-atom molecules they are much smaller than the isotropic term 30 . For example, the isotropic exchange interaction in a single pentacene molecule splits the singlet and triplet levels by about 1 eV 8 while the anisotropic interaction splits the triplet levels by about 1 GHz or ∼ 10 −6 eV 31 . The antisymmetric term is usually negligible in aromatic hydrocarbons and so we ignore it.
To simplify the hamiltonian, we decompose it into intra-chromophore and inter-chromophore interactions (see Supplementary Information S1). Because we are not interested in modeling, for example, the transitions between the singlet and triplet state on a single chromophore, we do not include the intra-chromophore isotropic exchange interaction in the hamiltonian. We do, however, keep the intra-chromophore anisotropic coupling.
Bridged Dimer www.nature.com/scientificreports/ We also include the isotropic inter-chromophore exchange interaction but discard the much smaller interchromophore anisotropic coupling. Finally, we set all of the inter-chromophore isotropic exchange interactions equal to the same number, J, (see Supplementary Information S1). The spin-exciton hamiltonian then takes a compact form,

Spin-Spin Interactions
where S A = s 1 + s 2 and S B = s 3 + s 4 are the spin operators associated with chromophores A and B. The first term in Eq. (6) is the isotropic exchange interaction between two triplet excitons.
By ′ are the hamiltonians associated with the intra-chromophore anisotropic interactions, taken to be of the spin-dipole form, written in the canonical "zero-field splitting" (ZFS) form from EPR literature 29 . The primed coordinates denote the principal axes of the magnetic dipole tensor for chromophore A and chromophore B. The highest symmetry case is the one we analyze here, where all principal directions for chromophore A are parallel to those of chromophore B.
In small organic chromophores like pentacene and tetracene, the single quantum triplet excited state is strongly localized 32 . The triplet states, each bound to a single chromophore, form molecular (Frenkel) excitons 33 . Like the spatial wavefunctions, the single quantum spin wavefunctions should also be tightly bound and have a definite triplet multiplicity. Equation (6) describes the interactions between two triplet Frenkel spin-excitons, modeled as two spin-1 objects, each one interacting with itself through the spin-dipole interaction and coupled to one another through exchange.
Because the hamiltonian depends only on the parameters J, D and E and the principal axes of A and B are parallel, we refer to the spin-exciton hamiltonian in Eq. (6) as the parallel JDE model. As we will show, restricting the axes of A and B to be parallel imposes symmetries that make relaxation between magnetic sublevels maximally state-selective. All of the parameters in the JDE model can either be computed in electronic structure or measured in mixed optical/magnetic resonance experiments.
While hamiltonians similar to Eq. (6) have appeared more recently in literature [14][15][16][17]19,28,34 , researchers either use approximations that make Eq. (6) spin-conserving or they analyze it at weak J. Methods of arriving at a spin-conserving JDE model differ between authors, but they amount to making an effective spin approximation, such that, for example, S 2 Az + S 2 Bz ≈ S 2 z . One then discards the isotropic exchange term, which is spin-conserving anyway, and writes the hamiltonian with renormalized D and E parameters, where the A and B site-spin operators are replaced by total spin operators. The D and E zero-field splitting parameters in crystals are renormalized to the spatially averaged D * and E * 35,36 . In the "strong exchange limit", |J|/|D| ≫ 1 , D becomes D/3, but only for the quintets 15 . Reference's 37 "singlet character" approximation is not spin-conserving, but is only appropriate when the Zeeman term dominates over J. In this theory the initial and final states are adiabatically connected to one another 28 . The theory presented in this manuscript is a nonadiabatic analog that is applicable when J is large. Spin conserving approximations, like the strong or weak exchange approximation, are suitable for work that does not consider transitions between states of different multiplicity. They are manifestly incapable of describing the kind of intersystem crossing, from 1 TT to 5 TT M for example, that recent EPR experiments have observed. These newer measurements necessitate the development of the theory presented here.
Spin dynamics and selection rules at zero applied field. The exchange term JS A · S B is the largest energy scale in the hamiltonian. It is rotationally invariant and diagonal in the total spin representation. We therefore take the diagonal elements of the JDE hamiltonian in the basis of total S 2 and S z as our reference hamiltonian, H 0 , (see Supplementary Information S3). The couplings between states of different multiplicity are off-diagonal perturbations, V, that cause transitions between these levels. They depend solely on D and E. The terms that comprise V are bilinear products of S A and S B Cartesian spin operators, and in the parallel case the largest of these terms ∼ DS Az S Bz . It has the form of an effective "Zeeman" interaction, where the z-component of the magnetic field produced by exciton spin-A couples to the magnetic dipole of spin-B.
The energy level diagram for the diagonal states of the reference hamiltonian, H 0 , appears in Fig. 3a. States of different multiplicity are split by the large inter-chromophore exchange interaction, J, and for J < 0 , the singlet state is higher in energy than the triplet and quintet manifolds (Eq. 6). The singlet-quintet splitting between 1 TT and 5 TT is three times the singlet-triplet splitting. The much smaller intra-chromophore axial ZFS interaction, D, splits magnetic, M, sublevels within each manifold. The rhombicity parameter, E-the smallest energy scale in the hamiltonian-gives rise to fine structure splittings between M-sublevels of H 0 .
When evaluating the hamiltonian, we express it in terms of spherical tensor operators and apply the Wigner-Eckhart theorem (see Supplementary Information S2) 38,39 . In our evaluation, the "renormalization" of D, for example, to D/3 for quintets (Fig. 3a), is a direct consequence of the Wigner-Eckart theorem. We assume that the ZFS, D and E parameters are independent of nuclear coordinates and that the time-dependent nuclear motion appears in J. A recent paper that numerically simulates spin dynamics in singlet fission in Ref. 28 makes a similar approximation.
The application of the Wigner-Eckhart theorem to the hamiltonian in Eq. (6) also shows that transitions from 1 TT to 5 TT 0 and 5 TT ±2 are allowed while all other transitions from 1 TT are forbidden (see Supplementary  Information S3). Forming linear combinations from the sums and differences of each pair of degenerate states breaks the degeneracy of the 5 TT ±1 states by 2|E| (Fig. 3b). These linear combinations of the M-sublevels are analogous to those of the ℓ = 2 spherical harmonics that are taken to construct the d-orbitals 40 , and so we label the states accordingly (Fig. 3b). Motivated by degenerate perturbation theory, we form symmetric and antisymmetric combinations of the ±M-states of the degenerate hamiltonian (Fig. 3b). In the new basis, the | 1 TT� state only couples to |Q x 2 −y 2 � = | 5 TT 2 � + | 5 TT −2 � / √ 2 and to |Q z 2 � = | 5 TT 0 � . These states are not simply the eigenstates of H ZFS = H A + H B . The state-selectivity follows from the symmetries of how the Cartesian Q-states, In a rate theory for relaxation, the relaxation rate from 1 TT to Q z 2 is proportional to D 2 and the rate from 1 TT to Q x 2 −y 2 is proportional to E 2 (see Supplementary Information S3). Because |D| ≫ |E| , the dominant relaxation channel is from 1 TT to Q z 2 (Fig. 3b). These selection rules are strict for identical chromophores with parallel symmetry, but upon breaking this symmetry, transitions from the quintet 5 TT to the triplet manifold 3 TT become allowed while transitions from the 1 TT to the 3 TT manifold remain forbidden. In the absence of parallel symmetry, transitions also become allowed between 1 TT and all of the quintet sublevels.
Equation (6) can describe how the unpairing of 1 TT to independent triplets T + T takes place (Eq. 4) 28 . By measuring quantum beats in delayed fluorescence spectra, Ref. 11 showed, rather convincingly, that exciton unpairing can occur in crystalline tetracene. In a crystal, the exciton hopping rate is fast. As a result, J, which depends sensitively on the distance between chromophores, will rapidly go to zero. This scenario can be modeled using a quantum quench with the time-dependent hamiltonian H(t) = H A + H B + (1 − θ(t))JS A · S B , where θ(t) is the Heaviside step function, and the quench occurs for t > 0 . After the initial transients associated with decoherence and population relaxation phenomena subside, detailed balance shows that the system's reduced density matrix will factorize into a thermal product state. Because H A and H B commute, at times that are long compared to the spin relaxation times the density matrix becomes ρ ∼ exp (−βH)/Z = exp (−βH A ) exp (−βH B )/Z A Z B = ρ A ⊗ ρ B . The two triplet exciton states become formally independent.
Exciton entanglement diminishes as J becomes smaller. This implies that dimers will preserve entanglement on longer timescales than crystals with mobile excitons. In tetracene crystals, the loss of coherence has been observed to occur on timescales of tens of nanoseconds, which is not much longer than the quantum gate switching times given by the inverse characteristic EPR transition frequencies 11 . Once thermalized, these two triplets offer no quantum advantage over single triplets prepared through more standard intersystem crossing processes. Transitions from 1 TT are only allowed to two Q states (black). The Q z 2 state may be populated through a non-adiabatically fast interconversion process which goes as k fast ∼ D 2 . A much slower process k slow ∼ E 2 allows transition from the 1 TT state to the Q x 2 −y 2 state.

Scientific Reports
| (2020) 10:18480 | https://doi.org/10.1038/s41598-020-75459-x www.nature.com/scientificreports/ For quantum information applications, one should focus attention on dimers where chromophores are covalently bound or packed together as a minority component in a crystal so that the exciton hopping rate to the host is negligible 14 . We refer to the latter category as a "dilute crystal". In these systems, J may fluctuate about a nonzero value, but it cannot go to zero. We impose the condition that |J| ≫ |D| . For organic chromophores, such as polyacenes, it is also often the case that |D| ≫ |E|.
In dimers, transitions can occur from rare fluctuations of the bare energy gaps between states of total |S, M� . These fluctuations are driven by nuclear motions, and we assume that the energy gap embodied in a timedependent J obeys Gaussian statistics. This scenario is valid so long as |J| ≫ |D| and the energy gap obeys linear response with respect to the nuclear motions. The resulting theory is completely analogous to Marcus' theory of nonadiabatic electron transfer and the Förster-Dexter theory of exciton hopping, where the ZFS parameters play the role of the nonadiabatic coupling matrix elements. Transitions between initial and final states take the form k i→f = F|�i|V |f �| 2 , where F is a Franck-Condon weighted density of states. In principle, F incorporates a thermal factor between the various final states that is the result of summing over nuclear fluctuations. Detailed balance, however, gives the condition that k i→f /k f →i = exp(−β E i − E f ) . Given that the various quintet states are split by about 0.05 cm −1 , ignoring the temperature dependence of the prefactor is a safe approximation for temperatures above about 1 K. It is a straightforward matter to extend this analysis to a case where |J| < |D| , by first diagonalizing the JDE hamiltonian and then applying second-order perturbation theory in the exciton-heat bath coupling. This approach would resonate closely with the much earlier work in Ref. 41 on delayed fluorescence in molecular crystals that has since been applied to recent experiments 13 . EPR spectroscopy. Recent experiments have employed an optical pump/EPR probe scheme to observe the fate of exciton polarization following singlet fission [15][16][17][18][19][20][21][22]28,42 . Many of these experiments use field-swept EPR as the probe, where the system is subjected to a static magnetic field, B 0 , along the laboratory z-axis. The static field splits the magnetic sublevels while an oscillatory microwave field, B 1 , polarized in the xy-plane, induces transitions between them. In these experiments, one finds resonances as a function of the static field strength, B 0 .
Many experiments use X-band EPR (8-12 GHz) and for small organic chromophores this is in the strong field limit, where states are split by much more than D. To model these experiments, we introduce the Zeeman term, H Zeeman = gµ B B 0 S z , into Eq. (6) and choose the quantization axis along the lab, or Zeeman, z-axis (see Supplementary Information S4). The static Zeeman field splits states of different M but not states of different S; the 1 TT state is unaffected by the Zeeman field.
Because the Zeeman field induces splittings that are large compared to those of H ZFS , the reference hamiltonian, H 0 , changes. To construct it, we first project out the quintet block to find its eigenstates, |α� = M c M,α | 5 TT M � . These states, the adiabats (Fig. 5b), adiabatically follow B 0 . The hamiltonian is then re-expressed in the adiabatic basis, with the reference hamiltonian, H 0 = α |α�ǫ α �α| , and coupling to | 1 TT� defined accordingly. Transitions occur between the adiabatic sublevels and have a spectrum given by the Golden Rule I = α,β |�α|S x |β�| 2 (P α − P β )δ(ǫ α − ǫ β ) 43 , where P α = Tr(ρ|α��α|) is the population in state |α� . More details appear in Supplementary Information S7. While the time-dependence of the populations can be, and has been, measured, we focus attention on the "prompt" EPR spectrum that interrogates the initial population of the exciton magnetic sublevels immediately following singlet fission, where the short-time approximation P α ∼ |� 1 TT|H|α�| 2 is valid. With these provisions in place, there are no adjustable parameters for the calculated EPR spectra.
Changing the orientation of the dimer relative to the Zeeman field results in a nonperturbative change in the hamiltonian (see Supplementary Information S5). The transitions from 1 TT to the 5 TT ±1 sublevels, that were once symmetry forbidden, are now allowed and state selectivity diminishes. For dimers in a powder or frozen solution, the EPR signal is a sum over an ensemble of molecules that have a broad distribution of orientations with respect to the Zeeman field. The resulting spin polarization is scrambled, which leads to decoherence in the ensemble signal. In quantum information applications this is a source of noise, and it is therefore important not only to fix the molecular axes relative to one another, but also with respect to the laboratory axis.
For dimers with their principal axes fixed in space, the prompt EPR spectra as a function of the polar angle, θ , between z and z ′ , exhibit different numbers of peaks with different frequencies and signed relative intensities (Fig. 4). The coupling term, V, is a function of θ , and all quintet states may be directly accessible from the 1 TT state (see Supplementary Information S5). There is intensity-borrowing from the 5 TT 0 and 5 TT ±2 zero-field states into all 5 TT M sublevels as a function of θ . As θ goes from 0 • to 90 • , state selectivity for the relaxation from 1 TT monotonically decreases. The EPR signal is a much weaker function of the azimuthal angle, φ.
Reference 14 published TT EPR spectra for dilute crystals of pentacene molecules doped into a p-terphenyl matrix. The pentacene molecules adopted both parallel and herringbone configurations in the host, as observed in the pentacene crystal structure. This elegant design allowed them to disperse the dimers and fix their orientations in space. The 5 TT EPR spectrum contains contributions from dimers in both parallel and herringbone geometries. Figure 5a shows the field-swept energies of 5 TT M states for parallel chromophores, where the Zeeman field, B 0 , is directed along the principal x ′ -axis, which is the most well-resolved spectrum in their paper. See our full 9 × 9 hamiltonian in Supplementary Information S6. On this diagram, we indicate the diabats which are pure |S, M� states, along with the adiabats described above (Fig. 5b). The transitions occur in the strong field limit, far from the avoided crossings between the adiabats, which occur at much smaller values of the B 0 field. Indeed, in the vicinity in which the EPR transitions are recorded, the diabats and adiabats very nearly coincide.
We present our calculation for the prompt EPR spectrum in Fig. 5c using the parameters from Ref. 14 . As Fig. 5 shows, our model fits their data extremely well for the quintet part of the spectrum (left panel of Fig. 6 in Ref. 14 ), with one caveat. The polarization pattern (AEAE) we calculate for parallel molecules (Fig. 5c) matches the polarization of the spectrum that Ref. 14  www.nature.com/scientificreports/ included an anisotropic inter-exciton exchange coupling and treated the populations as fit parameters. Ref. 14 did find, however, that setting the populations in quintet M = ±1 levels to zero produced the best fit to the spectrum. This result is completely consistent with the selection rules derived here with B 0 ||x ′ . Additional details appear in Supplementary Information S6.

Discussion
We have provided a derivation and an analysis of a model hamiltonian for singlet fission with an eye toward quantum computing, information, and sensing applications. The model is specified by only three parameters: the inter-chromophore isotropic exchange coupling, J, and the intra-chromophore ZFS axial parameter, D, and rhombicity parameter, E. These parameters can be measured independently or calculated using electronic structure. The model one arrives at under a set of reasonable approximations is something we call the JDE model, named for the J, D and E parameters of that hamiltonian. In particular, we have shown that one can use the magnetic sublevels of the 5 TT space as "qudits" in quantum information applications 44 , where EPR experiments perform the function of quantum gates. The five quintet states offer a quantum advantage over the three states of the spin-polarized triplets, produced either by intersystem crossing or as the final spin unpaired products of a singlet fission process in a crystal. We have shown the conditions under which the 1 TT state transfers to states in the quintet block and have given the conditions for maximal state selectivity, and thereby the most efficient pathway to optical spin polarization for those transitions.
To decrease the transition rates from the 2S+1 TT M manifold into the incoherent unpaired triplets, one needs to keep the value of J large. This implies that molecular dimers that are covalently bound to one another or doped as an impurity component into a host crystal are ideal candidates for generating optically spin-polarized quantum states near room temperature. . Calculated prompt EPR spectra for the rotation of an oriented and parallel sample with respect to an applied field. θ is the polar angle between the lab-fixed Zeeman axis and the principal z ′ -axis. Peak intensities are proportional to differences in state populations and are relative to intensities at θ = 0 • . Differences in the number and sign of the peaks are the result a θ-dependent coupling between the | 1 TT� and the adiabatic |α� states described in the text. Colored lines follow specific transitions, indicated in the inset where states are ordered in energy with respect to the applied field. The parameters used to calculate the spectra are consistent with those reported in Fig. 5, where we discuss the spectrum for θ = 90 • .
Scientific Reports | (2020) 10:18480 | https://doi.org/10.1038/s41598-020-75459-x www.nature.com/scientificreports/ We have identified, for the first time, the selection rules for relaxation between the various doubly excited TT levels in chromophores with parallel symmetry at both zero and large Zeeman fields. At zero field, fluctuations in J transfer population from the 1 TT state into the maximally entangled quintet state, Q z 2 . This transition rate goes as D 2 . There is one and only one other allowed transition, which is into the Q x 2 −y 2 state, but this rate is proportional to E 2 and is much slower. One can make the relaxation even more state selective by synthesizing molecules with large |D|/|E| ratios.
In the strong field conditions, characteristic of both time-resolved field swept EPR experiments and quantum computing applications, we find that relaxation can be kept state-selective provided that the principal axes of the two chromophores are parallel to each other and to the Zeeman field. When a molecule's principal z ′ -axis does not align with the Zeeman axis, several symmetries are broken, and transitions are possible to all sublevels in the quintet block. In samples where the molecules have a broad distribution of orientations relative to the Zeeman axis, the ensemble will exhibit decoherence. This is a different source of decoherence than, for example, inhomogeneous broadening, that can be removed through echo techniques. This means that one needs to also devise a method to immobilize and control the orientation of the singlet fission chromophores. Recent, elegant, experimental work has shown that this is possible 14 .
Finally, using our model and analysis, we calculated the prompt EPR spectra of pentacene dimers doped into a p-terphenyl crystal to compare with the experimental measurements reported in Ref. 14 . With the parameters J, D and E provided, there are no adjustable parameters in this calculation. The results match the major features of the measured quintet spectrum, including the polarization pattern (AEAE), very well.
Singlet fission can create strongly spin-polarized products and thereby generate nearly pure quantum states at room temperature, but there are several design principles that one should follow. First, keep the inter-chromophore exchange, J, large. Second, immobilize the molecules and align their principal axis to the Zeeman field. Both requirements are satisfied in immobilized and oriented, covalently linked dimers, or in dimer pairs that are embedded in a crystal host that inhibits exciton diffusion 14 .
Singlet fission can offer many of the quantum advantages found in color centers, like the NV center in diamond, but with a bottom-up approach to design that is currently unavailable in color centers whose defects are implanted in the material post-synthesis. Singlet fission, in contrast, is able to capitalize on the arsenal of  14 . The experiments performed in Ref. 14 are in the strong field limit, where the transitions are far removed from the avoided crossings between adiabatic states (b). In this limit, the diabatic states (solid lines) of the Zeeman hamiltonian are very close to the adiabatic states (dashed lines). Open circles indicate computed transitions which compare favorably to the field values of transitions reported by Ref. 14 (closed circles). With B 0 ||x ′ , the 1 TT state transfers only to the quintet M = 0, ±2 states. Wavy arrows indicate the direction of transitions for D > 0 . (c) Simulated prompt EPR spectrum, where positive changes in the intensity indicate induced absorption and negative changes indicate stimulated emission. The relative peak intensities were calculated from coupling matrix elements between the singlet and quintet states using a short-time approximation to the Pauli master equation. These populations reproduce the AEAE polarization pattern for the quintet part of the spectrum reported in Ref. 14 .

Scientific Reports
| (2020) 10:18480 | https://doi.org/10.1038/s41598-020-75459-x www.nature.com/scientificreports/ synthetic techniques developed in organic chemistry to design molecules. Our work provides a quantitative model for computing dynamics and fitting spectra, and qualitative design principles for the synthetic design of new organic molecules for quantum information applications. This is an important step in establishing the relationship between molecular structure and function in an emerging class of organic, novel quantum materials.