The analytic pair product approximation (APPA) now links the chemical structure of molecular spin qubits to their Hahn-echo coherence times, but only in the high-field limit and only for baths of spin-I=12I=\frac{1}{2} protons. The open problem, stated explicitly by the source authors, is to extend this analytic, design-rule-generating framework to nuclei with I>12I>\frac{1}{2} (quadrupolar nuclei) and to heteronuclear baths — which requires going beyond the current APPA equations — and to carry out the systematic investigation of crystalline host matrices needed to create controllable sparse bath-spin distributions, extending quantitative design rules into the 100 μs100\,\mu\mathrm{s}–millisecond coherence regime targeted for metal-organic-framework and defect-mimicking molecular-qubit architectures.

Background

A molecular spin qubit is a magnetically addressable electron spin hosted in a synthetic molecule, typically an organometallic complex or organic radical; its transverse coherence time T2T_2, often measured by the two-pulse Hahn echo sequence in pulsed EPR (electron paramagnetic resonance), bounds its usefulness for quantum sensing and quantum information. In molecular crystals the dominant decoherence mechanism at low temperature is often spectral diffusion: the electron spin couples via hyperfine interactions to the surrounding bath of nuclear spins (protons, deuterons, metal nuclei), and slow nuclear spin pair flip-flops (mutual exchange of two nuclear spin orientations via their dipole–dipole coupling) stochastically shift the electron's resonance frequency and decay the echo. Cluster-correlation expansion (CCE) is a systematic numerical method that decomposes the bath dynamics into contributions from clusters of 1, 2, 3, … nuclei; it handles arbitrary nuclear spins but yields numbers rather than chemical insight. Witzel, de Sousa and Das Sarma (Phys. Rev. B 72, 161306 (2005)) developed the analytic pair theory of spectral diffusion for quantum-dot electron spins, and Jeschke (J. Magn. Reson. Open 14–15, 100094 (2023)) formulated the corresponding nuclear-pair treatment of the electron spin echo (coining the term analytic pair product approximation, APPA). Suchaneck, Tesi, van Slageren and Köhn (JACS Au 6, 3927–3938 (2026)) showed that in the high-field limit (above roughly 1 T), where single-nucleus hyperfine-anisotropy contributions are quenched and only nuclear pair flip-flops decay the echo, APPA reproduces full generalized-CCE simulations for two crystalline benchmark molecular qubit systems and turns into quantitative design rules: which proton substitution patterns, isotope substitutions and local environments lengthen the echo. The authors restrict their theory to spin-I=12I=\frac{1}{2} nuclei and state that generalizing to other nuclear spins, in particular quadrupolar nuclei with I>12I>\frac{1}{2} (nuclei with spin >12>\frac{1}{2} possess an electric quadrupole moment coupling to local electric-field gradients, adding Hamiltonian terms that are nonsecular, i.e. do not commute with the dominant Zeeman interaction and are normally dropped in high-field secular approximations, absent for protons), would require going beyond the APPA equations; they further calculate that the spatial range of bath nuclei contributing to the echo grows strongly with echo delay time (roughly 15 Å at 10 μs10\,\mu\mathrm{s}, far larger at 100 μs100\,\mu\mathrm{s}), so that at long times the average proton density of the wider crystalline environment becomes decisive and dedicated, systematically engineered host matrices with sparse bath-spin distributions are needed. Numerical tools for arbitrary-spin baths already exist (the PyCCE package of Onizhuk and Galli, Adv. Theory Simul. 4, 2100254 (2021)), and host-material screening has been demonstrated for solid-state point defects (Kanai et al., Proc. Natl. Acad. Sci. 119, e2121808119 (2022), over 12,000 hosts), but no analytic pair-level theory or chemically actionable design framework exists for molecular qubits in molecular crystals beyond the proton, spin-12\frac{1}{2} case.

Problem Statement

As stated in the Conclusions of Suchaneck, Tesi, van Slageren and Köhn (JACS Au 6, 3927 (2026)), the framework that links chemical structure to molecular spin qubit coherence times must be extended along two jointly supporting directions: (1) Theory — generalize the analytic pair product approximation for Hahn-echo decoherence beyond its current restriction to spin-I=12I=\frac{1}{2} nuclei in the high-field limit, deriving pair-level echo-decay expressions that include nuclei with I>12I>\frac{1}{2} and their quadrupolar couplings to electric-field gradients — nonsecular terms that the high-field secular approximation discards and that must be retained within the high-field regime — and heteronuclear baths with different gyromagnetic ratios, reduce to the existing APPA equations when all bath nuclei have I=12I=\frac{1}{2}, and are validated against cluster-correlation-expansion simulations on realistic crystal structures; (2) Host-matrix engineering — carry out the systematic investigation of crystalline host matrices needed to create controllable sparse distributions of bath spins, turning the observation that the contributing-nucleus range grows with echo delay time into quantitative design rules (which matrices, isotope substitutions and substitution patterns suppress bath-induced decoherence) that remain valid at 100 μs100\,\mu\mathrm{s} and longer, millisecond-scale coherence times. A minimal resolution requires contributions to both directions: for direction (1), the derivation, reduction to the all-I=12I=\frac{1}{2} APPA equations, and CCE validation stated above; for direction (2), quantitative design rules established for a concrete crystalline host-matrix and isotope-substitution class and shown to suppress bath-induced decoherence at 100 μs100\,\mu\mathrm{s} and beyond. An exhaustive high-throughput survey of host matrices is beyond this minimal standard rather than a condition of completion. Both deliverables serve one overarching objective: a predictive, chemistry-actionable theory of how molecular and crystalline structure set the coherence time of molecular spin qubits, applicable beyond proton baths and to the long-time regime relevant for metal-organic-framework and defect-mimicking architectures.

References

  1. Sarah Suchaneck, Lorenzo Tesi, Joris van Slageren, Andreas Köhn, "Linking Chemical Structure and Coherence Times of Molecular Spin Qubits", JACS Au 6(7), 3927–3938 (2026), doi:10.1021/jacsau.6c00489, https://doi.org/10.1021/jacsau.6c00489 (preprint: doi:10.26434/chemrxiv-2025-lhdf5)
  2. W. M. Witzel, R. de Sousa, S. Das Sarma, "Quantum theory of spectral-diffusion-induced electron spin decoherence", Phys. Rev. B 72, 161306 (2005), doi:10.1103/PhysRevB.72.161306, https://doi.org/10.1103/PhysRevB.72.161306
  3. Gunnar Jeschke, "Nuclear pair electron spin echo envelope modulation", J. Magn. Reson. Open 14–15, 100094 (2023), doi:10.1016/j.jmro.2023.100094, https://doi.org/10.1016/j.jmro.2023.100094
  4. M. Onizhuk, G. Galli, "PyCCE: A Python Package for Cluster Correlation Expansion Simulations of Spin Qubit Dynamics", Adv. Theory Simul. 4, 2100254 (2021), doi:10.1002/adts.202100254, https://doi.org/10.1002/adts.202100254
  5. S. Kanai, F. J. Heremans, H. Seo, G. Wolfowicz, C. P. Anderson, S. E. Sullivan, M. Onizhuk, G. Galli, D. D. Awschalom, H. Ohno, "Generalized scaling of spin qubit coherence in over 12,000 host materials", Proc. Natl. Acad. Sci. 119, e2121808119 (2022), doi:10.1073/pnas.2121808119, https://doi.org/10.1073/pnas.2121808119
  6. Timothy J. Krogmeier, Anthony W. Schlimgen, Kade Head-Marsden, "A Perturbative Non-Markovian Treatment to Low-Temperature Spin Decoherence", J. Phys. Chem. Lett. 17, 4900–4905 (2026), doi:10.1021/acs.jpclett.6c00125, arXiv:2601.09651, https://arxiv.org/abs/2601.09651
  7. Boning Li, Yifan Quan, Xufan Li, Guoqing Wang, Robert G. Griffin, Avetik R. Harutyunyan, Paola Cappellaro, "Exploring the mechanisms of transverse relaxation of copper(II)-phthalocyanine spin qubits", arXiv:2511.03199 (2025), https://arxiv.org/abs/2511.03199
  8. S. L. Bayliss, P. Deb, D. W. Laorenza, M. Onizhuk, G. Galli, D. E. Freedman, D. D. Awschalom, "Enhancing Spin Coherence in Optically Addressable Molecular Qubits through Host-Matrix Control", Phys. Rev. X 12, 031028 (2022), doi:10.1103/PhysRevX.12.031028, https://doi.org/10.1103/PhysRevX.12.031028
  9. C. Ryan, V. Briganti, C. Hogan, M. O'Neill, A. Lunghi, "Spin decoherence in molecular crystals: Nuclear vs electronic spin baths", J. Chem. Phys. 163, 134104 (2025), doi:10.1063/5.0278659, https://doi.org/10.1063/5.0278659