Radical pairs may play a role in xenon-induced general anesthesia

Understanding the mechanisms underlying general anesthesia would be a key step towards understanding consciousness. The process of xenon-induced general anesthesia has been shown to involve electron transfer, and the potency of xenon as a general anesthetic exhibits isotopic dependence. We propose that these observations can be explained by a mechanism in which the xenon nuclear spin influences the recombination dynamics of a naturally occurring radical pair of electrons. We develop a simple model inspired by the body of work on the radical-pair mechanism in cryptochrome in the context of avian magnetoreception, and we show that our model can reproduce the observed isotopic dependence of the general anesthetic potency of xenon in mice. Our results are consistent with the idea that radical pairs of electrons with entangled spins could be important for consciousness.

Understanding consciousness remains one of the big open questions in neuroscience 1 , and in science in general. The study of anesthesia is one of the key approaches to elucidating the processes underlying consciousness 2,3 , but there are still significant open questions regarding the physical mechanisms of anesthesia itself 4 .
One anesthetic agent that has been studied extensively is xenon. Xenon has been shown experimentally to produce a state of general anesthesia in several species, including Drosophila 5 , mice 6 , and humans 7 . While the anesthetic properties of xenon were discovered in 1939 8 , the exact underlying mechanism by which it produces anesthetic effects remains unclear even after decades of research 9 . Our focus here is on this underlying physical mechanism, and there are important hints provided by two recent publications.
First, Turin et al. 5 showed that when xenon acts anesthetically on Drosophila specific electron spin resonance (ESR) signals, consistent with free electrons, can be observed. Based on this observation, Turin et al. 5 proposed that the anesthetic action of xenon may involve some form of electron transfer. They supported their proposal by density-functional theory (DFT) calculations showing the effect of xenon on nearby molecular orbitals. Second, Li et al. 6 showed experimentally that isotopes of xenon with non-zero nuclear spin had reduced anesthetic potency in mice compared with isotopes with no nuclear spin.
If the process by which xenon produces anesthetic effects includes free-electron transfer as well as nuclear-spin dependence, a mechanistic framework proposed to explain xenon-induced general anesthesia should possess these characteristics. Here we show that a model involving a radical pair of electrons (RP) and the subsequent modulation of the RP spin dynamics by hyperfine interactions (HFIs) is consistent with these assumptions.
The radical pair mechanism (RPM) was first proposed more than 50 years ago 10 . The rupture of a chemical bond can create a pair of electrons, where the electrons are localized on two different molecular entities, whose spins are entangled in a singlet state 11 . The magnetic dipole moment associated with electron spin can interact and couple with other magnetic dipoles and external magnetic fields, including HFIs and Zeeman interactions 12 . As a consequence of such interactions, the initial singlet state can evolve into a more complex state that has both singlet and triplet components. Eventually, the coherent oscillation of the RP between singlet and triplet states ceases and the electrons may recombine (for the singlet component of the state) or diffuse apart to form various triplet products 13 . The coherent spin dynamics and spin-dependent reactivity of RPs allow magnetic interactions which are six orders of magnitude smaller than the thermal energy, k B T , to have predictable and reproducible effects on chemical reaction yields 14 The RPM has become a prominent concept in quantum biology 16 . In particular, it has been studied in detail for the cryptochrome protein as a potential explanation for avian magnetoreception 14,[17][18][19][20][21][22] . In the present work we apply the principles and methods used to investigate cryptochrome 15 to xenon-induced general anesthesia.
General anesthetics produce widespread neurodepression in the central nervous system by enhancing inhibitory neurotransmission and reducing excitatory neurotransmission across synapses 23,24 . Three ligand-gated ionchannels in particular have emerged as likely molecular targets for a range of anesthetic agents 25 : the inhibitory glycine receptor, the inhibitory γ-aminobutyric acid type-A (GABA A ) receptor, and the excitatory N-methyl-Daspartate (NMDA) receptor 26 . NMDA receptors require both glycine and glutamate for excitatory activation 27 , and it has been suggested that xenon's anesthetic action is related to xenon atoms participating in competitive inhibition of the NMDA receptor in central-nervous-system neurons by binding at the glycine binding site in the NMDA receptor 25,26,28,29 . Armstrong et al. 25,28 propose that the glycine binding site of the NMDA receptor contains aromatic residues such as tryptophan (Trp), similar to those found in cryptochrome, which is consistent with the possibility that the anesthetic action of xenon involves a spin-dependent process similar to the RPM thought to occur in cryptochrome. Let us note that other targets have also been proposed for many anesthetics, e.g. tubulin 30 . In the following we focus on the NMDA receptor to be specific, but our model is quite general and could well apply to other targets.
We suggest that the (possibly partial) electron transfer related to xenon's anesthetic action that is evidenced by Turin et al. 5 plays a role in the recombination dynamics of a naturally occurring RP and that, for isotopes of xenon with a non-zero nuclear spin, this nuclear spin couples with (at least one of) the electron spins of such a RP, affecting the reaction yields of the RP and hence xenon's anesthetic action. Such a mechanism is consistent with the experimental results of Li et al. 6 that xenon isotopes with non-zero nuclear spin have reduced anesthetic potency compared to isotopes with zero nuclear spin.
Aromatic residues are important for the binding of xenon and glycine at the glycine binding site of the NMDA receptor. While Armstrong et al. 28 have shown that phenylalanine residues within the NMDA receptor are necessary for the binding of xenon in the active site, the redox inactivity of phenylalanine 31,32 makes it unlikely to be involved in RP formation. Trp, however, is redox active, as evidenced by its participation in RP formation in the context of cryptochrome, and could feasibly be involved in the formation of a RP here. It has been suggested that water (a source of oxygen) may be present in the NMDA receptor 25,27 , and Aizenman et al. 33 as well as Girouard et al. 34 suggest that reactive oxygen species (ROS) may be found within the NMDA receptor. Further, Turin and Skoulakis 35 found that when a sample of xenon gas was administered to Drosophila without oxygen gas present in the sample, no spin changes were observed in the flies. Motivated by these observations, we propose that a Trp residue located in the glycine binding site of the NMDA receptor could be oxidized by a nearby ROS, forming a [O ·− 2 TrpH ·+ ] RP. We consider a scenario in which a xenon atom interacts with the O ·− 2 radical electron, shifting some of the electron density from the oxygen onto the xenon atom (see also "Methods", Fig. 10), which results in HFI between the xenon nuclear spin and the O ·− 2 radical electron spin. The xenon nuclear spin is consequently able to influence the spin-dynamics and ultimate product yields of the [O ·− 2 TrpH ·+ ] RP, provided that the hyperfine interaction is sufficiently strong compared to the spin relaxation and reaction rates. Similarly to our above remark on the target receptor, here we focus on one specific potential radical pair for concreteness, but other implementations of the RPM may also be viable explanations for the mentioned experimental observations.
Given that cryptochrome is one case in which the RPM has been studied extensively, we propose that by recognizing commonalities in the biological and chemical environments in which magnetoreception and xenoninduced general anesthesia are thought to take place, analytical and numerical techniques that have been used to study cryptochrome may be adapted and applied to the case of xenon-induced anesthesia, potentially providing insight into general anesthetic mechanisms. We explore the feasibility of such a mechanism by determining and analyzing the necessary parameters and conditions under which the spin-dependent RP product yields can explain the experimental isotope-dependent anesthetic effects reported by Li et al. 6 .

Results
Predicting experimental xenon anesthesia results using the RPM model. Quantifying experimental anesthetic potency. In the work of Li et al. 6 , a metric referred to as the "loss of righting reflex ED50" (LRR-ED50) was defined using the concentration of xenon administered to mice, in which the mice were no longer able to right themselves within 10 s of being flipped onto their backs. The LRR-ED50 metric was reported to be correlated with consciousness in mice, and was measured experimentally for 132 Xe, 134 Xe, 131 Xe, and 129 Xe to be 70(4)%, 72(5)%, 99(5)%, and 105(7)%, respectively 6 .
Here we define the anesthetic potency as the inverse of the LRR-ED50 metric. In order to quantify the relative anesthetic potency of the various xenon isotopes, the potency of 132 Xe (with I = 0 ) is normalized to 1. The inverse of the LRR-ED50 value of isotopes 131 Xe and 129 Xe is then divided by that of 132 Xe. The relative isotopic anesthetic potencies are quantified as Pot 0 = 1 , Pot 3/2 = 0.71 (8) , and Pot 1/2 = 0.67(8) for 132 Xe, 131 Xe, and 129 Xe, respectively, as shown in Table 1, where Pot 0 is the relative potency of xenon with nuclear spin I = 0 , and likewise for Pot 3/2 and Pot 1/2 .
RPM model. We have developed a RPM model to predict anesthetic potency by making a connection with the relative singlet yield for different isotopes, as described in more detail below. The model that we have used here was developed using a hypothetical xenon-NMDA receptor RP system, based on the information about xenon action sites mentioned previously and shown in Fig. 1. Our model was also developed by considering the cryptochrome case as related to magnetoreception. Here we have proposed that the xenon-NMDA receptor RP system may involve xenon, ROS, and Trp, located in the glycine-binding site of the NMDA receptor. A spin- 2 radical electron (electron A) spin couples with non-zero xenon nuclear spin. The number of xenon atoms located in the active site when anesthetic action takes place is not yet completely clear, and the work of Dickinson et al. 25 suggests that the number of xenon atoms simultaneously present in the active site ranges probabilistically between zero and three, with single-xenon occupation proposed to be most probable. Here we show that the simplest and most probable case of a single xenon atom occupying the active site allows us to explain the anesthetic potency ratios derived from the experimental results of Li et al. 6 The additional degrees of freedom implicit in more complex Hilbert spaces, such as the cases of two and three-xenon occupation states, only aid the model in explaining the experimental results of Li et al. 6 We therefore focus our analysis on the single-xenon case.
The Hamiltonian of the RPM in the case of xenon-induced anesthesia depends not only upon Zeeman interactions and the number of xenon atoms present in the glycine binding site of the NMDA receptor, but also upon the assumed HFIs. Due to the likely randomized orientation of the relevant receptors in the brain, we consider only the isotropic Fermi contact contribution to the HFIs. In the simulation involving a single xenon atom occupying the active site, it is assumed that the O ·− 2 radical electron couples only to the xenon nucleus. Most of our calculations below are based on the simplifying assumption that the TrpH ·+ electron (electron B) couples only to the nuclear spin of the indole nitrogen in TrpH ·+ , following the treatment of Hore 15 , but we also briefly consider a separate scenario in which electron B couples only to the nuclear spin of the TrpH ·+ β-hydrogen, which we select because it has the largest isotropic hyperfine coupling constant of all 14 nuclear spins in tryptophan according to Maeda et al. 36 . The Hamiltonian describing these interactions is given as where Ŝ A and Ŝ B are the spin operators of radical electrons A and B, respectively, Î 1 is the nuclear spin operator of the xenon nucleus, Î 2 is the nuclear spin operator of the TrpH ·+ residue (where I = 1 for the indole nitrogen, and I = 1/2 for the β-hydrogen), a 2 = γ e a ′ 2 is the hyperfine coupling constant describing the strength of the HFI  www.nature.com/scientificreports/ between radical electron B and the nuclear spin in TrpH ·+ ( a ′ 2 = 347 µ T for the indole N, and a ′ 2 = 1600 µ T for the β-hydrogen 36 ), a 1 = (γ n /|γ129 Xe |)γ e a ′ 1 is the hyperfine coupling constant describing the HFI strength between radical electron A and the xenon nucleus, where γ n is the nuclear gyromagnetic ratio of a given xenon isotope, γ129 Xe is the nuclear gyromagnetic ratio of 129 Xe, γ e is the gyromagnetic ratio of an electron, and ω is the Larmor precession frequency of the electrons about an external magnetic field 15 . The Larmor precession frequency is defined as ω = γ e B , where B is the external magnetic field strength (assumed to be along the z-axis). In the calculation of a 1 , the correction factor of (γ n /|γ129 Xe |) is included to account for the unique nuclear gyromagnetic ratios of the various xenon isotopes, where the magnitude of this correction factor is normalized to 1 for 129 Xe. The isotopic xenon gyromagnetic ratios are taken to be 37 γ129 Xe = −7.441 × 10 7 rad(sT) −1 , γ131 Xe = 2.206 × 10 7 rad(sT) −1 , and γ132 Xe = 0 . Consequently, the hyperfine coupling constant a ′ 1 is multiplied by factors [−1, (2.206/7.441), 0] when calculated for each respective isotope. It should be pointed out that we focus here on the HFIs previously described, and that interactions between the two electron spins, as well as potential interactions between the electron spins and other nuclei are neglected.
Using DFT to calculate hyperfine coupling constants. As a means of checking if the prediction results of our xenon-induced anesthesia model are reasonable, we use DFT to study the interaction between O ·− 2 and xenon, and to calculate the value of the expected hyperfine coupling constant a ′ 1,exp . Depending on the functional and basis-set used, we find a range of values a ′ 1,exp ∈ [955, 2644] µ T. We also use DFT analysis to calculate a value Figure 2. The dependence of the agreement between relative anesthetic potency and singlet yield ratio on changes in the hyperfine coupling constant a ′ 1 and the RP spin-coherence relaxation rate r for a ′ 1 ∈ [0, 5000] µ T and r ∈ [5.0 × 10 5 , 1.0 × 10 9 ] s −1 , using a ′ 2 = 347 µ T, B = 50 µ T, and k = 5.5 × 10 7 s −1 . The singlet yield ratio SR is calculated using a 1 = (γ n /|γ129 Xe |)γ e a ′ 1 , where the magnitude of the isotopic nuclear gyromagnetic correction factor, (γ n /|γ129 Xe |) , is normalized for 129 Xe. The model can explain the experimentally derived 6 relative anesthetic potency of xenon for values of a ′ 1 and r where |Pot 1/2 − SR 1/2 |, |Pot 3/2 − SR 3/2 | ≤ 0.08 . (a) The absolute difference between Pot 1/2 and SR 1/2 . (b) The absolute difference between Pot 3/2 and SR 3/2 . www.nature.com/scientificreports/ of a ′ 2,exp = 347 µ T for the hyperfine coupling between the radical electron on TrpH ·+ and the nitrogen nuclear spin. Our value for the hyperfine coupling between the radical electron on TrpH ·+ and the β-hydgrogen nuclear spin is taken from Maeda et al. 36 . More details about our DFT analysis can be found in "Methods".
Determination of singlet yield ratios. The eigenvalues and eigenvectors of the Hamiltonian can be used to determine the ultimate singlet yield ( S ) for all times much greater than the radical-pair lifetime ( t ≫ τ) 15 : where, following the methodology used by Hore 15 in the context of cryptochrome, M is the total number of nuclear spin configurations, P S is the singlet projection operator, |m� and |n� are eigenstates of Ĥ with corresponding energies of ω m = �m|Ĥ|m� and ω n = �n|Ĥ|n� , respectively, k = τ −1 is inverse of the RP lifetime, and r = τ −1 c is the inverse of the RP spin-coherence lifetime.
The spin-dependent RP singlet yield is calculated for each xenon isotope under consideration. The singlet yield ratio (SR) for each xenon nuclear spin value is then calculated by dividing the singlet yield obtained using the given spin value by the singlet yield obtained using spin I = 0 , resulting in the ratio of spin-1/2 singlet yield to spin-0 singlet yield being expressed as SR 1/2 , and the ratio of spin-3/2 singlet yield to spin-0 singlet yield expressed as SR 3/2 , with the singlet yield ratio of spin-0 being normalized to SR 0 = 1 . The calculated singlet yield ratios are compared with the xenon potency ratios (Pot) derived using the data reported by Li et al. 6 , as described above.
Here we analyze the sensitivity of the singlet yield ratios to changes in the HFI between radical electron A and the xenon nucleus ( a ′ 1 ), RP reaction rate (k), RP spin-coherence relaxation rate (r), and external magnetic field strength (B). The dependence of the quantity |Pot − SR| on a ′ 1 and r is shown in Figs. 2 and 3, the dependence of |Pot − SR| on a ′ 1 and k is shown in Figs. 4 and 5, and the dependence of |Pot − SR| on r and k is shown in Figs. 6 and 7. The relationship between SR and B can be seen in Fig. 8. We also look at the dependence of the singlet yield ratios on the RP spin-coherence relaxation rate and the RP reaction rate in the case where electron B couples with the TrpH ·+ β-hydrogen. The dependence of |Pot − SR| on r and k in the electron B-β-hydrogen case can be seen in Fig. 9. Our goal is to find regions in parameter space such that the spin-dependent singlet yield ratios match the anesthetic isotopic potency ratios, i.e., the quantities |Pot 1/2 − SR 1/2 | and |Pot 3/2 − SR 3/2 | should be smaller than the experimental uncertainties in the anesthetic potency.

Discussion
The main question that we set out to answer in this study was whether the RPM can provide a viable explanation for the isotope effects observed by Li et al. 6 . Our above results answer this question in the affirmative. A simple radical pair model with reasonable parameter values can reproduce the observed anesthetic potencies.
To be concrete, we proposed a specific scenario involving a [O ·− 2 TrpH ·+ ] RP, where the oxygen molecule interacts with a single xenon atom. This proposal is motivated by the suggestions that Trp is present in the relevant Figure 4. The dependence of the agreement between relative anesthetic potency and singlet yield ratio on changes in the hyperfine coupling constant a ′ 1 and the RP reaction rate k for a ′ 1 ∈ [0, 5000] µ T and k ∈ [5.0 × 10 5 , 1.0 × 10 9 ] s −1 , using a ′ 2 = 347 µ T, B = 50 µ T, and r = 5.0 × 10 6 s −1 . The singlet yield ratio SR is calculated using a 1 = (γ n /|γ129 Xe |)γ e a ′ 1 , where the magnitude of the isotopic nuclear gyromagnetic correction factor, (γ n /|γ129 Xe |) , is normalized for 129 Xe. The model can explain the experimentally derived 6 relative anesthetic potency of xenon for values of a ′ 1 and k where |Pot 1/2 − SR 1/2 |, |Pot 3/2 − SR 3/2 | ≤ 0.08 . (a) The absolute difference between Pot 1/2 and SR 1/2 . (b) The absolute difference between Pot 3/2 and SR 3/2 .  28 and that ROS may be found within the NMDA receptor 33,34 , as well as the observation that oxygen should be present when xenon is administered to Drosophila in order to observe the ESR signals that correlate with xenon-induced general anesthesia 35 . The RP spin-coherence relaxation rate, RP reaction rate, and hyperfine coupling parameters were varied to reproduce the experimental results of Li et al. 6 , and the resulting range of values of the hyperfine coupling constant a ′ 1 was found to agree with the range of values calculated for xenon in the proximity of O ·− 2 using DFT analysis. One potential challenge for the specific scenario involving O ·− 2 are the requirements on the spin relaxation rate r. Hogben et al. 38 and Player and Hore 39 have discussed a radical pair involving O ·− 2 in the context of magnetoreception. These authors argue that for free O ·− 2 in solution the spin relaxation lifetime is expected to be below 1 ns due to molecular rotation. The relaxation rate requirements that we find from our calculations come within about an order of magnitude of this value, as shown in Figs. 2, 3, 6 and 7. The gap widens to about two orders of magnitude if one considers β-hydrogen as the dominant nuclear spin for the electron on Trp instead of nitrogen, as shown in Fig. 9. However, the same authors have also pointed out that the oxygen spin relaxation rate could be reduced if the biological environment reduces the molecular symmetry and inhibits molecular rotation. Moreover Kattnig 40 has shown that the involvement of a paramagnetic scavenger can significantly reduce the constraints on the spin relaxation rate in the context of magnetoreception, which suggests the possibility that the same principle may apply to isotope effects in the present context. Alternatively, other RPs not including O ·− 2 might also be able to explain the isotope effects.
To gain a deeper understanding of the RP formation and recombination dynamics within our proposed scenario, it would be of interest to perform DFT modelling of the electron transfer between ROS and Trp with and without the presence of xenon. Such modelling could be expanded on by exploring the molecular dynamics of the binding site using quantum mechanics/molecular mechanics (QM/MM) simulation techniques, similar to those used in the case of cryptochrome 41 , to more clearly understand the stability of the radicals and other characteristics of the currently proposed RP system. In particular it would be of interest to attempt to theoretically reproduce the experimental ESR results of Turin et al. 5 .
When considering the sensitivity of the model to changes in external magnetic field, as shown in Fig. 8, the range of B values that produce agreement between SR and Pot is given as B ∈ [0, 2525] µ T, and includes the geomagnetic field at different geographic locations (25 to 65 µT) 42 . This result indicates that for external field values vastly larger (by approximately 2 orders of magnitude) than the geomagnetic field, the anesthetic potency of xenon may be significantly different than that observed by Li et al. 6 ; specifically, the potency of 131 Xe may be greater than experimentally observed. It would be of interest to investigate the experimental effects of the external magnetic field strength on xenon-induced general anesthesia in vivo. For example, such an experiment could In conclusion, our results suggest that xenon-induced general anesthesia may fall within the realm of quantum biology, and be similar in nature to the proposed mechanism of magnetoreception involving the cryptochrome protein 17 .
This also raises the question whether the action of other anesthetics involves similar mechanisms to the one proposed here. It would be interesting to explore isotopic nuclear-spin effects as well as magnetic field effects in experiments with other general anesthetic agents that are thought to function similarly to xenon, such as nitrous oxide and ketamine 43,44 .
General anesthesia is clearly related to consciousness, and it has been proposed that consciousness (and other aspects of cognition) could be related to large-scale entanglement [45][46][47][48] . RPs are entangled and could be a key Figure 6. The dependence of the agreement between relative anesthetic potency and singlet yield ratio on the relationship between r and k for r ∈ [1.0 × 10 6 , 1.0 × 10 8 ] s −1 and k ∈ [5.0 × 10 6 , 5.0 × 10 8 ] s −1 , using a ′ 1 = 1896 µ T, a ′ 2 = 347 µ T, and B = 50 µ T. The singlet yield ratio SR is calculated using a 1 = (γ n /|γ129 Xe |)γ e a ′ 1 , where the magnitude of the isotopic nuclear gyromagnetic correction factor, (γ n /|γ129 Xe |) , is normalized for 129 Xe. The model can explain the experimentally derived 6 relative anesthetic potency of xenon for values of r and k where |Pot 1/2 − SR 1/2 |, |Pot 3/2 − SR 3/2 | ≤ 0.08 . (a) The absolute difference between Pot 1/2 and SR 1/2 . (b) The absolute difference between Pot 3/2 and SR 3/2 .  49 . Viewed in this-admittedly highly speculative-context, the results of the present study are consistent with the idea that general anesthetic agents, such as xenon, could interfere with this large-scale entanglement process, and thus with consciousness.

Methods
DFT analysis. We used the Gaussian package 50 with B3LYP functional and TZVP basis set for our DFT calculation of a ′ 2,exp . The range of values calculated for a ′ 1,exp using DFT analysis varied widely ( a ′ 1,exp ∈ [955, 2644] µT), depending on the functional and basis set used. The calculation is also very sensitive to the separation between O ·− 2 and Xe, and the potential energy is minimized for a separation of ∼ 3 Å. Ultimately, the ORCA package 51 was used for our Xe-O ·− 2 DFT calculations, and the molecular structure was optimized using the dispersion-corrected PBE0 functional and def2-TZVP basis set. The resulting calculation gives a ′ 1,exp = 1767 µ T, which is close to the value of a ′ 1 for which our model minimizes the quantity CE = (|Pot 3/2 − SR 3/2 | + |Pot 1/2 − SR 1/2 |). The orbitals obtained from the molecular structure optimization were used to calculate orbital energies as well as the hyperfine coupling constant a ′ 1,exp , using the B3LYP functional and def2-TZVP basis set. Relativistic effects were treated by a scalar relativistic Hamiltonian using the zeroth-order regular approximation (ZORA). Solvent effects were included with the conductor-like polarizable continuum model (CPCM), using a dielectric constant of 2.
Our DFT computations show that the unpaired Xe-O ·− 2 electron is bound, and that the highest occupied molecular orbital (HOMO) resides primarily on the O 2 molecule but is extended slightly onto the xenon atom, as shown in Fig. 10.   Figure 7. The dependence of the agreement between relative anesthetic potency and singlet yield ratio on the relationship between r and k for r ∈ [1.0 × 10 6 , 1.0 × 10 8 ] s −1 and k ∈ [5.0 × 10 6 , 5.0 × 10 8 ] s −1 , using a ′ 1 = 1896 µ T, a ′ 2 = 347 µ T, and B = 50 µ T. The singlet yield ratio SR is calculated using a 1 = (γ n /|γ129 Xe |)γ e a ′ 1 , where the magnitude of the isotopic nuclear gyromagnetic correction factor, (γ n /|γ129 Xe |) , is normalized for 129 Xe. The model can explain the experimentally derived 6 relative anesthetic potency of xenon where |Pot 1/2 − SR 1/2 | ≤ 0.08 and |Pot 3/2 − SR 3/2 | ≤ 0.08 intersect, shaded in yellow. www.nature.com/scientificreports/ Isotropic hyperfine coupling is a result of the Fermi contact interaction between a nucleus and electron 52 , and here this is a measure of the contribution of the valence s-orbital to the molecular orbital occupied by the unpaired radical electron. If the unpaired electron resides mainly in a p-or a d-orbital, the Fermi contact interaction can also be significant where the spin density at the nucleus is induced by configuration interaction or spin polarization 53 . In our case, the unpaired electron mainly occupies the p-orbital which results in a significant isotropic HFI. The isotropic hyperfine coupling constant can be calculated as 54,55 where g e is the electron g-factor, µ 0 is the permeability of free space, and |�(0)| 2 is the electron probability density evaluated at the nucleus.

Data availability
The datasets generated and analysed during the current study are available from the corresponding author on reasonable request.