Fundamental limits of electron and nuclear spin qubit lifetimes in an isolated self-assembled quantum dot

Combining external control with long spin lifetime and coherence is a key challenge for solid state spin qubits. Tunnel coupling with electron Fermi reservoir provides robust charge state control in semiconductor quantum dots, but results in undesired relaxation of electron and nuclear spins through mechanisms that lack complete understanding. Here, we unravel the contributions of tunnelling-assisted and phonon-assisted spin relaxation mechanisms by systematically adjusting the tunnelling coupling in a wide range, including the limit of an isolated quantum dot. These experiments reveal fundamental limits and trade-offs of quantum dot spin dynamics: while reduced tunnelling can be used to achieve electron spin qubit lifetimes exceeding 1 s, the optical spin initialisation fidelity is reduced below 80%, limited by Auger recombination. Comprehensive understanding of electron-nuclear spin relaxation attained here provides a roadmap for design of the optimal operating conditions in quantum dot spin qubits.


INTRODUCTION
Semiconductor quantum dots (QDs) offer excellent quantum optical properties and well-defined quantum states of individual spins-an attractive combination for quantum information processing devices 1 . Recent proof-of-concept demonstrations with QDs include heralded entanglement of two remote spins 2 , generation of photonic cluster states 3 , spin-controlled photon switching 4 as well as implementation of electron-nuclear quantum spin interfaces 5 and nuclear spin quantum computing 6 . The stability of the spin states, measured by their lifetimes, is crucial in all these applications. Quantum dot is described by a central spin of a single charge (electron or hole) coupled to N ≈ 10 3 −10 5 nuclear spins via hyperfine interaction 7 . The lack of translational motion combined with the mismatch in electron and nuclear spin energies suppresses relaxation 8 , providing long spin lifetimes required for spin qubits. However, thorough understanding of spin relaxation is complicated by the multitude and complexity of the residual environment couplings, which include electron-phonon interactions [9][10][11][12] , quadrupolar coupling of nuclear spins to strain 13 , nuclear spin diffusion, and electron cotunneling 14,15 arising from proximity of the Fermi reservoir. Moreover, impurity charge traps 16,17 adjacent to QDs degrade spin qubit lifetimes. Thus, it remains an open question to establish the maximum (intrinsic) spin lifetimes that can be achieved at any given magnetic field and temperature, as opposed to spin relaxation arising from QD device design and imperfections.
Phonon-assisted electron spin relaxation enabled by spin-orbit interaction is a dominant mechanism 9,10,18 at high magnetic field B z ≳ 2 T, but the limit to electron spin lifetime T 1,e at low fields remains unexplored. In case of nuclear spins, cotunnelingmediated relaxation was identified as dominant mechanism 15 , while direct verification is lacking, since bias control of cotunneling is restricted to a narrow range compatible with single-electron QD state.
Here we study a series of structures where electron cotunneling is controlled directly by the thickness t B of the tunnel barrier separating the dot from the Fermi reservoir. We find that at B z ≳ 2 T and temperatures θ > 4.2 K nuclear spin relaxation is dominated by a higher-order process assisted by phonons 19,20 and noncollinear hyperfine interaction 13 , rather than by cotunneling, which is dominant only at low fields B z ≲ 2 T. Electron spin lifetimes exceeding T 1,e > 1 s are found at B z ≈ 0.4 T, with a fundamental maximum T 1,e ≈ 20 s estimated for an isolated dot at θ = 4.2 K, bounded by phonon relaxation and direct hyperfine interaction at high and low magnetic fields, respectively. While coupling to Fermi reservoir degrades T 1,e , it is shown to play a crucial role in counteracting Auger recombination 21 and enabling electron spin initialisation with near-unity fidelity 22 . Figure 1a sketches conduction band energy profile which is controlled with external bias V S to tune an InAs QD into Coulomb blockade regime, where it is charged deterministically 23 with a single electron (1e). The trion state with two electrons and one hole can be accessed through resonant optical excitation, and observed in resonance fluorescence (ResFl). Magnetic field B z along the sample growth axis splits the electron spin-up " j i and spin-down # j i energies (Fig. 1b), enabling selective excitation of the optically allowed transition between " j i electron and the trion "#* j i with a spin-up hole * j i. Weak recombination enabled by hyperfine and heavy-light hole mixing β ≪ 1 can 'shelve' the dot 14 into # j i state, quenching ResFl intensity I ResFl until electron returns to " j i through a spin flip with rate ξ ↑↓ . Such shelving provides an efficient way both for initialisation and readout of the electron spin 14,22 . Furthermore, hyperfine interactionĤ hf / ðÎŝÞ of electron spin s with nuclear spins I enables electron-nuclear flip-flops, so that repeated electron spin initialisation creates a net nuclear spin polarisation P N , which can be monitored through optically measured hyperfine shifts E hf in the splitting of the " j i and # j i states 7,17 . Electron and nuclear spin dynamics are measured using a pump-delay-probe protocol (Fig. 1c). The decay of optically pumped electron [nuclear] spin polarisation over dark period T Dark is probed by measuring I ResFl [E hf ]. Figure 1d shows an example of time-resolved ResFl, which is used to derive the residual electron spin polarisation P e after a delay T Dark (see the "Methods" section). Measurements of P e at different T Dark reveal electron spin relaxation (symbols in Fig. 1e), while examples of nuclear spin relaxation E hf (T Dark ) are shown in Fig. 1f. Fitting (lines in Fig. 1e, f) is used to derive the intrinsic spin-relaxation rates of electron Γ e = 1/T 1,e = 2ξ ↑↓ and nuclei Γ N = 1/T 1,N .

Effect of Auger recombination on spin initialisation
We make a systematic comparison of spin dynamics in a thinbarrier sample (t B = 37 nm, Fig. 2a-c), similar to structures used previously 10,14,15 , and a thick-barrier structure (t B = 52 nm, Fig.  2d-f), approximating a QD isolated from the Fermi reservoir. Examining the bias dependence of continuous excitation resonance flourescence intensity I ResFl in a thin-barrier sample, we observe a 1e plateau at B z = 0 T (triangles in Fig. 2a), while at high B z (circles and squares in Fig. 2a) ResFl is strongly suppressed, indicating spin 'shelving' 10,14 . A striking difference is observed in a thick-barrier sample (Fig. 2d), where fluorescence intensity and spin 'shelving' contrast are reduced (Fig. 2a), which may at first suggest the lack of electron spin initialisation. However, this is ruled out by time-resolved ResFl (e.g. Fig. 1d), which reveals spinpumping fluorescence pulses of similar intensity for all t B . We ascribe the reduction in I ResFl and the peculiar two-stage electron spin decay in a thick-barrier sample (squares in Fig. 1e) to the Auger process 21,24 , where electron-hole recombination ejects the second electron with a rate γ A (Fig. 1b). Following the ejection, an empty QD does not contribute to ResFl, hence P e ≈ 1 observed initially. During T Dark an unpolarised electron can return from the Fermi reservoir with recharging rate r, giving rise to the fast component of the P e (T Dark ) decay (squares in Fig. 1e at T Dark < 100 μs), whereas the slow component corresponds to resident electron spin relaxation with rate Γ e .
Using rate equation modelling (see details in the "Methods" section) of the four-level system shown in Fig. 1b, we find good description of the experiments (dashed line Fig. 1e) and derive r ≈ 1.26 × 10 5 s −1 , Γ e ¼ T À1 1;e % 3:3 s −1 . Importantly, the level P e ≈ 0.77 reached after dot recharging (T Dark ≈ 100 μs) gives a direct measure of the electron spin initialisation fidelity in a thickbarrier sample, revealing the fundamental limitations arising in an isolated (t B → ∞) QD. In a thin-barrier sample, Auger recombination is counteracted by fast recharging: the resulting spin initialisation fidelity is higher, but can never reach unity. The maximum spin initialisation fidelity is an algebraic function of QD properties such as relaxation rates and heavy-light hole mixing. (The exact expression can be found in Supplementary Note 2.) Analysis shows that fidelity is improved for faster recharging r, larger trion mixing β and slower spin flip ξ ↑↓ . Conversely, in the limit of infinitely slow recharging r → 0 spin initialisation becomes Electron and nuclear spins in a quantum dot. a Schematic of an InAs quantum dot embedded in a n-i-Schottky diode structure. Electron spin is coupled to nuclei via hyperfine interaction (Ĥ hf ) and to phonons (Ĥ eÀph ). Tunnel barrier thickness t B controls the cotunneling interaction (Ĥ cotun ). Quantum dot charge state is controlled with bias V S , which tunes the energies of one-electron (1e) and two-electron (2e) states with respect to Fermi energy E F . b Energy levels of an empty dot 0 j i, electron with spin up ( " j i, s z = +1/2) or down ( # j i, s z = −1/2), and a trion "#* j iwith hole moment j z = +3/2, which has a small admixture β ≪ 1 of a trion "#+ j iwith opposite hole moment j z = −3/2. Arrow labels show the rates of radiative recombination γ R , Auger recombination γ A , recharging r, and electron spin flip ξ ↑↓ . Electron spin splitting is due to Zeeman effect (μ B g e B z , where g e is electron g-factor and μ B is Bohr magneton) and nuclear hyperfine shift E hf ∝ A hf P N . c Timing diagram of a pump-delay-probe experiment where optical excitation is used to initialise and probe the spins of either the electron or nuclei. impossible (P e → 0), imposing a practical lower limit on the tunnel coupling with the Fermi reservoir. Figure 2 shows that electron (Fig. 2b, e) and nuclear (Fig. 2c, f) spin relaxation rates are reduced at the centre of the 1e plateau 15,25 . This Coulomb blockade regime is of most interest, as it corresponds to a stable electron spin qubit, and is examined in more detail in Fig. 3. The dependence of Γ e on magnetic field is well described (solid lines in Fig. 3a) by

Fundamental limits of electron spin lifetimes
where for the field-dependent mechanism assisted by spin-orbit interaction and phonons we find Γ e;ph % 2:27 ± 0:48 s À1 T Àk ph and k ph ≈ 4.1 ± 0.13 in both samples. The exponent is in good agreement with k ph = 4 predicted 26 and observed experimentally 9 for this mechanism in high-temperature regime μ B g e B z ≲ k B θ, where phonon thermal occupation factor gives rise to an additional / B À1 z factor in Γ e . This condition is well satisfied for our experiments at B z ≤ 8 T, θ ≥ 4.5 K and typical g-factor values |g e | ≈ 0.4. By contrast, in previous studies at B z ≤ 12 T, θ ≈ 1 K the onset of low-temperature regime was observed 9 , where phonon thermal occupation factor is ≈1, resulting in k ph = 5.
Cotunneling involves virtual injection of a second electron into the dot, followed by return of a spin-flipped electron to the Fermi reservoir 14,27 . The fitted cotunneling-induced relaxation rate in a thin-barrier (t B = 37 nm) sample Γ e,cotun ≈ 532 ± 65 s −1 = (0.0019 ± 0.0002 s) −1 is larger than Γ e,cotun ≈ (1.65 ± 0.21 s) −1 found for a thicker t B = 52 nm. Since Γ e,cotun is field independent, the increase in Γ e at very low B z ≲ 0.2 T (squares in Fig. 3a) is likely due to energy-conserving electron-nuclear flip-flops, which become allowed when electron Zeeman energy is comparable to nuclear quadrupolar energy. By extrapolating the phonon (dotted line) and hyperfine (dash-dotted line) mechanisms we roughly estimate the fundamental minimum of the electron spin relaxation rate in an isolated (t B → ∞) QD as Γ e;min \ð20 sÞ À1 , expected to occur at B z ≈ 0.4 T for θ ≈ 4.5 K. Similarly slow electron spin relaxation rates were reported in strain-free III-V QDs 28 , although at lower temperatures θ < 0.1 K. Figure 3b shows Γ e measured in cotunneling-dominated lowfield regime in samples with different barriers. A considerable variation between individual dots for the thin barrier t B = 37 nm can be due to random atomic-scale positioning of the individual Si dopants 29 and Si segregation 30 at the interface between tunnel barrier and Fermi reservoir. By contrast, in a thick-barrier sample (t B = 52 nm) the dot is coupled to a large number of dopants, smoothing out atomic-scale variations and leading to consistent Γ e .
In low-field regime the barrier thickness t B controls both the electron spin relaxation rate Γ e and recharging rate r, but we find that r exceeds Γ e by approximately five orders of magnitude, as exemplified in Fig. 1e. The recharging of an empty QD is a firstorder tunnelling process, whereas cotunneling in a charged QD is a second-order process 31,32 , which qualititatively explains the difference in rates. Moreover, the charge state of the dot itself may affect the conduction band energy profile, altering the tunnel ResFl intensity, 1e 2e . At non-zero B z = 1 T (circles) and B z = 8 T (squares) the centre of the 1e plateau is suppressed due to spin 'shelving' , leaving two peaks arising from resonant tunnel coupling to the Fermi reservoir. b Bias dependence of the electron spin relaxation rate measured at B z = 1 T (circles) and B z = 8 T (squares). c Bias dependence of the nuclear spin relaxation rate measured for the same QD as in a and b. Note the offset in the horizontal axis, which is likely due to additional electric fields arising from the optically generated trapped charges near the QD. Vertical dashed lines separate different charge states 0e, 1e, 2e. d-f Same as a-c but for a QD from a thick-barrier sample (t B = 52 nm). All error bars are 95% confidence intervals.
coupling. An accurate first-principle quantitative description of the rates would require development of a detailed theoretical model.

Nuclear spin relaxation mechanisms
The marked difference in nuclear spin relaxation rate Γ N (solid symbols in Fig. 3a) of the two samples at low magnetic fields suggest cotunneling as the dominant channel, whereas at B z ≳ 2 T cotunneling is negligible. This is in contrast to previous studies under similar conditions (t B = 35 nm, B z = 5 T), which identified cotunneling and nuclear spin diffusion 15 as dominant mechanisms. We examine diffusion by varying the spin pumping time (squares in Fig. 3c): Taking the difference in Γ N at short and long pumping, the diffusion rate is estimated to be as small as . Slow diffusion is due to quadrupolar freezing of nuclear spin flip-flops 15,17 . Moreover, optical pumping through wetting layer states used in this work (as opposed to resonant optical pumping in ref. 15 ), is likely to polarise nuclear spins not only in the dot but also in its vicinity, further suppressing the diffusion. Thus cotunneling and spin diffusion alone do not account for all the relevant nuclear spin relaxation mechanisms.
To explain the entire Γ N (B z ) dependence, we treat the average spin of the QD electron as a random process. Uniquely for selfassembled QDs, noncollinear hyperfine interaction permits nuclear spin relaxation without electron spin flip 13 -this mechanism is expected to be more efficient than direct electron-nuclear spin flips 15 . The transition rate Γ N,ij between states i j i and j j i of a single nuclear spin is proportional to spectral power density 33 of the fluctuating electron spin s z (t) at the nuclear spin transition frequency ν ij . Using first-order perturbation theory 13 we have where A hf is the hyperfine constant, M ij ¼ hijÎ z jji is the matrix element of the nuclear spin operatorÎ z , and electron spin correlation time is approximated by τ e ≈ 1/Γ e . This model describes a higher-order nuclear spin relaxation process, mediating by electron spin relaxation, which in turn is dominated by phonons at high fields or cotunneling at low fields. In the high field limit jM ij j 2 / B À2 z , Γ e / B 4 z and ν 2 ij / B 2 z (see the "Methods" section) leading to Γ N / const, which agrees with the weak field dependence of Γ N observed for both samples at B z ≳ 4 T (solid symbols in Fig. 3a). Moreover, temperature dependence at high field (Fig. 3d) is close to linear Γ N ∝ θ at θ ≲ 15 K, matching the Γ e ∝ θ dependence 34 of the underlying phonon-mediated electron spin relaxation process. Superlinear growth of Γ N at θ ≳ 15 K is likely due to two-phonon processes 26,35 , with scaling predicted to range between ∝θ 2 and ∝θ 9 .
For quantitative description (see the "Methods" section) we use experimentally measured Γ e and estimate ν ij and M ij from magnetic resonance spectra 36,37 . The results (dashed lines in Fig.  3a) are in good agreement for the thin barrier (t B = 37 nm), where electron correlation time is short. The discrepancy for the thickbarrier sample (t B = 52 nm) is most prominent at B z ≲ 2.5 T, revealing the limitations of the electron-spin fluctuation model (Eq. (2)) in the previously unexplored regime of a nearly isolated long-lived electron spin. To examine the cause, we note that τ 2 e ν 2 ij ) 1 except for possible quadrupolar anti-crossings of the nuclear spin levels 37,38 , so that Eq. (2) can be rewritten as The right side of this equation is a function of magnetic field and quantum dot structural properties, such as chemical composition and strain inhomogeneity, but it does not depend on tunnel coupling to Fermi reservoir. If τ e = 1/Γ e , the equation predicts independence of Γ N /Γ e on tunnel barrier. This is seen to be the case in Fig. 4 for samples with t B = 37 and 42 nm, whereas the thick barrier sample (t B = 52 nm) shows excessive Γ N /Γ e . The exact reason for the increased Γ N /Γ e in the thick barrier sample is not clear. One likely possibility is additional nuclear spin relaxation mechanisms where hyperfine interaction fluctuates without electron spin flips, resulting in τ e < 1/Γ e . For example, modulation of the hyperfine interaction can occur through electron wavefunction density shifts, arising from fluctuating electric fields of the itinerant carriers in the Fermi reservoir 33,39 , or charge traps 16 . Charge noise in the studied structures is indeed present and evidenced, e.g. by fluctuating electron spin relaxation rates at the edges of 1e plateau (circles in Fig. 3e). Future experiments using, e.g. bias modulation spectroscopy may elucidate the roles of different nuclear spin relaxation mechanisms and lead to more accurate theoretical models. Further improvements to nuclear spin relaxation description can be sought through a microscopic model that takes into account quadrupolar anti-crossings of the individual nuclear spin levels 37,38 , which may accelerate relaxation and reenable frozen spin diffusion. A contribution of direct nuclear-phonon interaction 33,40 is also possible, as its rate Γ N ≈ 10 −4 −10 −3 s −1 is comparable to the lowest Γ N observed here in electron-charged QDs.

DISCUSSION
Experiments presented here establish a comprehensive picture of electron-nuclear spin relaxation in self-assembled QDs in a wide range of practically accessible conditions. Present experiments require B z ≳ 0.15 T to initialise the spins 22 and resolve the Zeemansplit optical transitions for spin probing. Extension to lower fields could shed light on the less explored regime where electron spin relaxation abruptly slows down from Γ e ≈ 10 9 s −1 at zero field 8,41 to Γ e ≈ 1 s −1 observed here at 0.15 T. For the practically interesting range B z ≳ 0.15 T, electron spin relaxation is fundamentally limited by phonon coupling, which is similar in other types of QDs. Hence, electron spin lifetimes exceeding 1 s should be achievable in GaAs electrostatic 28 and epitaxial 42 QDs, as well as in II-VI QDs 43,44 . By contrast, nuclear spin relaxation studied here is specific to selfassembled III-V QDs, and is governed by noncollinear hyperfine interaction. All experiments here were conducted in Faraday geometry, whereas noncollinear interaction is expected to be even stronger for magnetic field tilted away from the sample growth axis 37 , which may lead to faster nuclear spin relaxation in Voigt geometry. The techniques employed here, can also be applied to establish the less explored fundamental limits of nuclear spin dynamics in electron-charged strain-free QDs 28,42 , where noncollinear interaction will be small, but nuclear spin diffusion might be more prominent.

Samples and experimental techniques
The samples are low-density InAs self-assembled QDs (≲1 QD per μm 2 ) grown on a GaAs substrate. The dot layer is positioned at the centre of a λ/ 2 optical cavity formed by a bottom Bragg mirror consisting of 15 GaAs/ AlAs pairs and a top reflector with 2 pairs (estimated quality factor Q ≈ 60). Cavity mode is centred at 950 nm, which matches the long-wavelength tail of the QD wavelength distribution. The Fermi reservoir is formed by a doped GaAs layer (Si concentration of 1.1 × 10 18 cm −3 , thickness ≈ 80 nm). The doped layer is located beneath QDs and is separated by a GaAs layer of thickness t B = 37-52 nm, depending on the structure. Each sample is processed into a Schottky diode structure with an Au/(In-Ge) ohmic back contact 45 annealed from the top surface, and a 5 nm-thick semitransparent Ti Schottky top contact. External bias is applied to the top contract and controls the charge states of QDs. In order to form an electron spin qubit, the dot is charged deterministically with one electron (1e). This is achieved by tuning the energy of the 1e state to ≈10 meV below the Fermi energy E F (Fig. 1a), while the two-electron (2e) state remains depopulated, since its energy exceeds E F by ≈10 meV, which is ≳20 times the thermal energy k B θ at liquid helium temperature θ ≈ 4.2 K (Boltzmann constant k B ≈ 86.17 μeV K −1 ). The dot is then charged by an electron tunnelling from the Fermi reservoir.
The sample is mounted in a bath cryostat equipped with a superconducting coil producing magnetic field up to 8 T in Faraday geometry (field parallel to sample growth direction and optical axis z). An aspheric lens mounted near the sample is used for optical excitation of the QD and for light collection. Photoluminescence (PL) spectroscopy (see Supplementary Fig. 1) is used for initial QD characterisation. In nuclear spin dynamics experiments the dot is excited using diode lasers operating at 850 nm (resonant with InGaAs wetting layer). Nuclear spin polarisation (cooling) is achieved with a circularly polarised high power (≳100 times the power of ground state exciton saturation) laser, with typical pump pulse duration of T Pump ≈ 8 s. A short (T Probe ≈ 10 ms) low power (approximately corresponding to ground state exciton saturation) probe pulse is used to excite PL, which is then analysed on a double grating spectrometer to derive the hyperfine shifts E hf in the splitting of a QD Zeeman doublet. The relaxation of the nuclear spin polarisation is derived by measuring E hf in the probe as a function of delay T Dark between the pump and the probe. The resulting E hf (T Dark ) dependencies (e.g. Fig. 1f) are fitted with stretched or compressed exponentials / e ÀðTDark=T1;NÞ η , where η is the parameter describing stretching (η < 1) or compression (η > 1).
In ResFl experiments the dot is excited with a linearly polarised singlemode tunable diode laser. The scattered laser is rejected using crosspolarised detection 46 , and the collected fluorescence is directed to an avalanche photodiode detector, whose photon-counting pulses are measured with a pulse counter and a digital oscilloscope. Typical linewidths measured in continuous excitation ResFl spectra on a negatively charged trion are ≈0.5 GHz at low power (non-saturating excitation). Electron spin initialisation at finite magnetic field is witnessed through ResFl intensity I ResFl , which is significantly reduced when electron is initialised into the # j i state, taking the dot out of resonance with optical driving of the " j i $ "#* j i transition (compare squares and circles with triangles in Fig. 2a). Electron spin relaxation is accelerated and spin shelving is destroyed when the bias is tuned to the level where Fermi reservoir is resonantly tunnel-coupled with 1e (V S ≈ 0.1 V in Fig. 2a) or 2e (V S ≈ 0.33 V in Fig. 2a) quantum dot state, resulting in two peaks in the I ResFl (V S ) dependence (squares and circles in Fig. 2a).
In electron spin dynamics measurements pulsing of the resonant laser is achieved with acousto-optical modulators providing on/off ratio better than 10 7 . The power of the pump and probe pulses is close to ResFl saturation conditions and typical duration is T Pump ≈ T Probe ≈ 5-10 μs, which is significantly longer than the spin pumping time. As a result, timeresolved ResFl exhibits short pulses (Fig. 1d) with amplitudes I ResFl,Pump and I ResFl,Probe . At the start of each measurement cycle, and prior to optical pump pulse, the bias is adjusted for resonant electron tunnelling in order to counteract optical nuclear spin pumping and depolarise the electron (see further details in Supplementary Note 1). The rising edge of the fluorescence pulse corresponds to the rise time of the laser intensity in a pulse that pumps the initially depolarised electron. The exponentially decaying falling edge of the fluorescence pulse traces the gradual shelving (initialisation) of the dot into the # j i electron spin state. The loss of electron spin polarisation during the delay T Dark results in partial recovery of the fluorescence intensity measured in the probe pulse. The residual polarisation P e (i.e. electron spin polarisation at the start of the probe G. Gillard et al. normalised by polarisation at the end of the pump pulse) is then derived as P e = (I ResFl,Pump − I ResFl,Probe )/I ResFl,Pump . This way complete loss of electron polarisation (P e = 0) is observed as I ResFl,Probe = I ResFl,Pump , while I ResFl,Probe = 0 implies no loss (P e = 1) or Auger recombination that empties the dot. By measuring P e at different T Dark the decay of electron spin polarisation is obtained as shown in Fig. 1e by the symbols.
Unless stated otherwise, all error estimates in the text and error bars in figures are 95% confidence intervals.

Modelling of the electron spin relaxation dynamics
We simulate the dynamics of the four-level system shown in Fig. 1b using a simplified noncoherent rate equation model. The relaxation rates of all possible transitions are shown in Fig. 1b, and when resonant optical pumping is present we add a transition " j i ! "#* j i with rate P. We assume symmetric rates in electron spin flips " j i $ # j i, which is justified when electron Zeeman energy is smaller than the thermal energy k B θ. The system of first-order differential rate equations is ∂ ∂t In simulations we set γ R = 10 9 s −1 , which is typical for InAs/GaAs QDs 47 . In order to simulate the pump-delay-probe experiment we use initial population probabilities ðp 0 j i ; p " j i ; p # j i ; p "#* j i Þ ¼ ð0; 1=2; 1=2; 0Þ, and propagate the equations numerically over pump pulse (with P ≠ 0), dark delay (P = 0), and probe (P ≠ 0). The values of γ A , P, β, r, ξ ↑↓ are used as fitting parameters. The evolution of the trion state population p "#* j i ðtÞ reproduces the experimentally measured time-resolved ResFl (e.g. Fig.  1d). The simulated p "#* j i ðtÞ traces are integrated over pump and probe intervals and are used to calculate the residual electron polarisation P e in the same way I ResFl,Pump and I ResFl,Probe are used to calculate P e from experimental data (e.g. Fig. 1e). The fitting parameters are adjusted to achieve two-objective optimisation: one objective is to minimise the root mean square difference between simulated and experimental P e (T Dark ) traces, the other objective is to match the characteristic exponential time in the falling edge of the ResFl intensity produced by the pump pulse. An example of the Pareto-optimal fitted P e (T Dark ) is shown by the dashed line in Fig. 1e, in good agreement with experiment (squares). The calculated falling edge time (87 ns) is also in good agreement with the experimental value ≈91 ns.
From fitting we find r ≈ 1.26 × 10 5 s −1 (95% confidence interval [0.47 × 10 5 s −1 , 3.34 × 10 5 s −1 ]) and ξ ↑↓ ≈ 1.65 s −1 [1.16 s −1 , 2.34 s −1 ], which correspond to the characteristic timescales of the fast and slow components, respectively, in the two-stage decay (squares in Fig. 1e). For Auger rate we obtain γ A ≈ 1.09 × 10 7 s −1 [0.66 × 10 7 s −1 , 3.11 × 10 7 s −1 ]. This is approximately five times higher than γ A ≈ 0.23 × 10 7 s −1 reported previously from time-resolved ResFl experiments 21 . The discrepancy could be due to the difference in QD structures, high optical pump power used in our experiments, limitations of a non-coherent rate equation model, and uncertainty in the fitted parameters. The uncertainty is increased by the interdependence of the fitting parameters γ A , P, β, r, ξ ↑↓ , which is inevitable since five parameters are used to fit essentially four degrees of freedom (fast and slow rates of the two stage P e (T Dark ) decay, P e following fast decay and the characteristic time of the falling edges in fluorescence pulses). This uncertainty can also be understood to arise from the limited information provided by the ResFl measurement, which does not distinguish between spin shelving into the # j i state and Auger recombination into the 0 j i state, restricting the ability to monitor the full dynamics of the four-level system. The remaining best fit values are P ≈ 2.5 × 10 9 s −1 [0.34 × 10 9 s −1 , 3.16 × 10 9 s −1 ] and β ≈ 5.2 × 10 −3 [0.75 × 10 −3 , 20.1 × 10 −3 ].

Modelling of the nuclear spin relaxation rate
We start by noting that at all magnetic fields used in this study (B z = 0.15-8 T) nuclear spin relaxation in an empty dot (0e) is at least an order of magnitude slower than at the centre of the 1e charging plateau. (In 0e regime we measure Γ N ≈ 6.0 × 10 −4 s −1 at B z = 0.15 T, which reduces at higher fields below the minimum measurable level of Γ N < 10 −4 s −1 .) This suggests that electron is the dominant mediator of nuclear spin relaxation in a Coulomb blockade regime (1e). Electron-nuclear coupling is given by the hyperfine Hamiltonian H hf ¼ X k A hf;k ðŝ xÎx;k þŝ yÎy;k þŝ zÎz;k Þ; (6) where the summation goes over all nuclei, A hf,k is the hyperfine constant of the kth nucleus, ðŝ x ;ŝ y ;ŝ z Þ are the components of the electron spin-1/2 operator and ðÎ x;k ;Î y;k ;Î z;k Þ are the components of the spin operator of the kth nucleus. The flip-flop term / ðŝ xÎx þŝ yÎy Þ of this interaction permits spin exchange between nuclear spin I and electron spin s, but at sufficiently large magnetic field, where electron Zeeman energy significantly exceeds nuclear Zeeman and quadrupolar spin splitting, such flip-flops are strongly suppressed. For a depolarised nuclear spin ensemble in an InAs QD, this threshold field would be on the order of 0.03 T, corresponding to statistical fluctuation / A hf;k ffiffiffi ffi N p of the nuclear hyperfine field acting on the electron. However, in self-assembled QDs the principal strain axis is generally misaligned 36,37 from the growth axis z, resulting in nuclear eigenstates which are superpositions ofÎ z eigenstates. Under these conditions nuclear spin states are mixed by the non-flip-flop part of the hyperfine interaction A hf;kÎz;kŝz . This noncollinear interaction enables transitions between nuclear spin states i j i and j j i without transfer of spin to the electron. Using first-order perturbation and Weisskopf-Wigner approximation one arrives to Eq. (2) for nuclear spin relaxation rate between the pair of states i j i and j j i, where we have assumed the same hyperfine constant A hf,k = A hf /N for all nuclei.
In order to calculate matrix element M ij ofÎ z we consider spin I = 3/2 and assume that the principal component of electric field gradient is characterised by quadrupolar shift frequency ν Q , and is tilted by angle α from the z-axis. Using first-order perturbation approach we calculate M ij in the opposite limits of small magnetic field (γ N B z ≪ ν Q , where γ N is the nuclear gyromagnetic ratio) and high magnetic field (γ N B z ≫ ν Q ). For each individual nucleus the pair ( i j i, j j i) with the largest off-diagonal matrix element is then selected leading to: For intermediate fields we interpolate the matrix element with a monotonic function: jMðB z Þj 2 ¼ 1=ðjMj À2 Bz !0 þ jMj À2 Bz !1 Þ, and the nuclear spin transition frequency is taken to be ν 2 ij ¼ ðγ N =2πÞ 2 ðB 2 z þ B 2 z;min Þ, where non-zero B z;min reflects the fact that at B z = 0 the spin states are split by nuclear quadrupolar effects. This model is a simplification since nuclear spin levels (anti)cross 38 at B z ≈ 2πν Q /γ N , where nucleus experiences a nearly zero effective magnetic field. The perturbative approach breaks down as M ij is enhanced and ν ij ≈ 0 at these anticrossing points (typically occurring at B z ≈ 1 T) 37 , hence the introduction of B z;min which softens the singularities in Eq. (2). While spin relaxation would be enhanced in such resonant nuclei, it would also prevent their optical polarisation, thus we effectively neglect their contribution to the overall measured nuclear spin decay in our simplified model.
In order to calculate the nuclear spin decay rate according to Eq. (2), we take τ e ¼ Γ À1 e , equivalent to assuming that electron spin flips are the only source of noise acting on the nuclear spins. For Γ e we use Eq. (1) taking best fit parameters for each QD sample. Based on NMR spectroscopy of similar QDs 36,37 , we model strain inhomogeneity within the quantum dot by considering a uniformly distributed quadrupolar shift ν Q ∈ [0, 16] MHz with principal axis uniformly distributed on a part of a sphere with α ∈ [0°, 76°]. The large values of α ≈ 76°account for As nuclei, whose quadrupolar shifts are dominated by atomic scale disorder, arising from Ga and In alloying 36 . The gyromagnetic ratio is also varied uniformly γ N /2π ∈ [7.4, 9.2] MHz T −1 to account for five different isotopes 48 present in the dot ( 113 In, 115 In, 69 Ga, 71 Ga, 75 As). We take B z;min ¼ 0:06 T and use an average value A hf = 50 μeV for all the isotopes 49 . The number of nuclei is taken to be N = 4 × 10 4 . The nuclear spin relaxation rate at a given magnetic field B z is then calculated by averaging over all parameter distributions to take into account the contributions of the individual nuclei in a QD.
Despite the simplifications, the model is in good agreement with experimental dependence Γ N (B z ) in the t B = 37 nm sample (thin dashed line in Fig. 3a) and t B = 42 nm sample (Supplementary Fig. 2). Qualitative