Experimentally testing quantum critical dynamics beyond the Kibble-Zurek mechanism

We experimentally probe the distribution of kink pairs resulting from driving a one-dimensional quantum Ising chain through the paramagnet-ferromagnet quantum phase transition, using a single trapped ion as a quantum simulator in momentum space. The number of kink pairs after the transition follows a Poisson binomial distribution, in which all cumulants scale with a universal power-law as a function of the quench time in which the transition is crossed. We experimentally verified this scaling for the first cumulants and report deviations due to noise-induced dephasing of the trapped ion. Our results establish that the universal character of the critical dynamics can be extended beyond the paradigmatic Kibble-Zurek mechanism, which accounts for the mean kink number, to characterize the full probability distribution of topological defects.


Introduction
The understanding of nonequilibrium quantum matter stands out as a fascinating open problem at the frontiers of physics. Few theoretical tools account for the behavior away from equilibrium in terms of equilibrium properties. One such paradigm is the so-called Kibble-Zurek mechanism (KZM) that describes the nonadiabatic dynamics across a phase transition and predicts the formation of topological defects, such as vortices in superfluids and domain walls in spin systems [1]. Pioneering insights on the KZM were conceived in a cosmological setting [2] and applied to describe thermal continuous phase transitions [3,4,5]. The resulting KZM was latter extended to quantum phase transitions [6,7,8,9]. Its central prediction is that the average total number of topological defects n T , formed when a system is driven through a critical point in a time scale τ Q , is given by a universal power-law n T ∼ τ −β Q . The power-law exponent β = Dν/(1 + zν) is fixed by the dimensionality of the system D and a combination of the equilibrium correlation-length and dynamic critical exponents, denoted by ν and z, respectively. Essentially, the KZM is a statement about the breakdown of the adiabatic dynamics across a critical point. As such, it provides useful heuristics for the preparation of ground-state phases of matter in the laboratory, e.g., in quantum simulation and adiabatic quantum computation [10]. It has spurred a wide variety of experimental efforts in superfluid Helium [11,12,13], liquid crystals [14,15], convective fluids [16,17], superconducting rings [18,19,20], trapped-ions [21,22,23], colloids [24], and ultracold atoms [25,26,27,28,29], to name some relevant examples. This activity has advanced our understanding of critical dynamics, e.g., by extending the KZM to inhomogeneous systems [30,31]. In the quantum domain, experimental progress has been more limited and led by the use of quantum simulators in a variety of platforms [32,33,34]. Beyond the focus of the KZM, the full counting statistics encoded in the probability distribution can be expected to shed further light. The number distribution of topological defects has become available in recent experiments [35]. In addition, the distribution of kinks has recently been explored theoretically in the one-dimensional transverse-field quantum Ising model (TFQIM) [36], a paradigmatic testbed of quantum phase transitions [37]. Indeed, it has been argued that the distribution of topological defects is generally determined by the scaling theory of phase transitions and thus exhibits signatures of universality beyond the KZM [36]. We aim at validating this prediction experimentally. The TFQIM is described by the Hamiltonian that we shall consider with periodic boundary conditionsσ z 1 =σ z N +1 and N even. Additionally,σ z m andσ x m are Pauli matrices at the site m and g plays the role of a (dimensionless) magnetic field that favors the alignment of the spins along the x-axis. This system exhibits a quantum phase transition between a paramagnetic phase (|g| 1) and a doubly degenerated phase with ferromagnetic order (|g| 1). There are two critical points at g c = ±1. We shall consider the nonadiabatic crossing of the critical point g c = −1 after initializing the system in the paramagnetic phase and ending in a ferromagnetic phase. Our quantum simulation approach relies on the equivalence of in the TFQIM in one spatial dimension and an ensemble of independent two-level systems. This seminal result forms the basis of much of the progress on the study of quantum phase transitions and can be established via the Jordan-Wigner transformation [38], Fourier transform and Bogoliubov transformation [39]. We detail these steps in the Supplementary Note 1, where it is shown that the TFQIM Hamiltonian can be alternatively written in terms of independent modes as [37,40] whereγ k are quasiparticle operators, with k labeling each mode and taking values k = π N (2m − 1) with m = − N 2 + 1, . . . , N 2 . The energy k of the k-th mode is k (g) = 2J (g − cos k) 2 + sin 2 k. Other physical observables can as well be expressed in both the spin and momentum representations. We are interested in characterizing the number of kinks. As conservation of momentum dictates that kinks are formed in pairs, we focus on the kink-pair number operator N ≡ N m=1 1 −σ z mσ z m+1 /4. The latter can be equivalently written asN = k>0γ † kγ k , wherê γ † kγ k is a Fermion number operator, with eigenvalues {0, 1}. As different k-modes are decoupled, this representation paves the way to simulate the dynamics of the phase transition in the TFQIM in "momentum space": the dynamics of each mode can be simulated with an ion-trap qubit, in which the expectation value ofγ † kγ k can be measured. To this end, we consider the quantum critical dynamics induced by a ramp of the magnetic field in a time scale τ Q that we shall refer to as the quench time. We further choose g(0) < −1 in the paramagnetic phase. In momentum space, driving the phase transition is equivalently described by an ensemble of Landau-Zener crossings. This observation proved key in establishing the validity of the KZM in the quantum domain [8,32,33]: the average number of kink pairs N = n after the quench scales as in agreement with the universal power law n KZM ∝ τ − ν 1+zν Q with critical exponents ν = z = 1 and one spatial dimension (D = 1).

Results
The full counting statistics of topological defects is encoded in the probability P (n) that a given number of kink pairs n is obtained as a measurement outcome of the observableN . To explore the implications of the scaling theory of phase transitions, we focus on the characterization of the probability distribution P (n) of the kink-pair number in the final nonequilibrium state upon completion of the crossing of the critical point induced by (3). Exploiting the equivalence between the spin and momentum representation, the dynamics in each mode leads to two possible outcomes, corresponding to the mode being found in the excited state (e) or the ground state (g), with probabilities p e = p k and p g = 1 − p k , respectively. Thus, one can associate with each mode k > 0 a discrete random variable, with excitation probability p k = γ † kγ k . The excitation probability of each mode is that of Bernoulli type. As the kink-pair number n is associated with the number of modes excited, P (n) is a Poisson binomial distribution, with the characteristic function [36] associated with the sum of N/2 independent Bernoulli variables. The mean and variance of P (n) are thus set by the respective sums of the mean and variance of the N/2 Bernoulli distributions characterizing each mode, n = k>0 p k and Var(n) = k>0 p k (1 − p k ). The full distribution P (n) can be obtained via an inverse Fourier transform. The theoretical prediction of the kink distribution P (n) relies on knowledge of the excitation probabilities {p k }, which can be estimated using the Landau-Zener formula p k = exp − π Jτ Q k 2 [8] (details are given in Supplementary Note 2). Experimental design. Experimentally, the excitation probabilities {p k } can be measured by probing the dynamics in each mode, that is simulated with the ion-trap qubit. The dynamics of a single mode k is described by a Landau-Zener crossing, that is implemented with a 171 Yb + ion confined in a Paul trap consisting of six needles placed on two perpendicular planes, as shown in Fig. 1a. The hyperfine clock transition in the ground state S 1/2 manifold is chosen to realize the qubit, with energy levels denoted by |0 ≡ |F = 0, m F = 0 and |1 ≡ |F = 1, m F = 0 , as shown in Fig. 1b.  Figure 1: Experiment setup. a A single 171 Yb + is trapped in a needle trap, which consists of six needles on two perpendicular planes. b The qubit energy levels are denoted by |0 and |1 , which is the hyperfine clock transition of the trapped ion. c The qubit is driven by a microwave field, generated by a mixing wave scheme in the high pass filter (HPF). Operations on the qubit are implemented by programming the arbitrary wave generator (AWG). The quantum critical dynamics of the one-dimensional TFQIM is detected by measuring corresponding Landau-Zener crossings governing the dynamics in each mode.
State preparation and manipulation: At zero static magnetic field, the splitting between |0 and |1 is 12.642812 GHz. We applied a static magnetic field of 4.66 G to define the quantization axis, which changes the |0 to |1 resonance frequency to 12.642819 GHz, and creates a 6.5 MHz Zeeman splittings for 2 S 1/2 , F = 1. In order to manipulate the hyperfine qubit with high control, coherent driving is implemented by a wave mixing method, see the scheme in Fig. 1c. First, an arbitrary wave generator (AWG) is programmed to generate signals around 200 MHz. Then, the waveform is mixed with a 12.442819 GHz microwave (generated by Agilent, E8257D) by a frequency mixer. After the mixing process, there will be two waves around 12.242 GHz and 12.642 GHz, so a high pass filter is used to remove the 12.224 GHz wave. Finally, the wave around 12.642 GHz is amplified to 2 W and used to irradiate the trapped ion with a horn antenna. Measurement: For a typical experimental measurement in a single mode, Doppler cooling is first applied to cool down the ion [41]. The ion qubit is then initialized in the |0 state, by applying a resonant light at 369nm to excite S 1/2 , F = 1 to P 1/2 , F = 1. Subsequently, the programmed microwave is started to drive the ion qubit. Finally, the population of the bright state |1 is detected by fluorescence detection with another resonant light at 369nm, exciting S 1/2 , F = 1 to P 1/2 , F = 0. Fluorescence of the ion is collected by an objective with 0.4 numerical aperture (NA). A 935nm laser is used to prevent the state of the ion to jump to metastable states [41]. The initialization process can prepare the |0 state with fidelity > 99.9%. The total error associated with the state preparation and measurement is measured as 0.5% [42]. In the Supplementary Note 1, we rewrite the explicit Hamiltonian of the k-th mode as a linear combination of a single two-levels system [8]. In this way, experimentally, the Hamiltonian of a single k-mode can be explicitly mapped to a qubit Hamiltonian describing a two-level system driven by a chirped microwave pulse,Ĥ where Ω R = 4J/ and ∆ k (t) = 4J[g(t) + cos k]/( sin k) are the Rabi frequency and the detuning of the chirped pulse, respectively. To measure the excitations probability after a finite ramp with normalized quench parameter A = Jτ Q / , we can vary g(t) from -5 to 0. For this purpose, a chirped pulse with length of T p = 5AT R sin k and g(t) = 5t/T p − 5 is used in the experiment, with T R = 1/Ω R denoting the Rabi time. A typical operation process driven by the programmed microwave is illustrated on the Bloch sphere, see Supplementary Notes 3 and 4 for details. First, the qubit is prepared in the ground state ofĤ TLS k at g(0) = −5. Then, g(t) is ramped from −5 to 0 with quenching rate 1/τ Q . Finally, the ground state of the finalĤ TLS k is rotated to qubit |0 , so that the excited state is mapped to |1 , which can be detected as a bright state. The measured excitation probability for different quench times is shown in Fig. 2. Figure 2: Excitation probability in the ensemble of Landau-Zener crossings. The quantum critical dynamics of the one-dimensional TFQIM is accessible in an ion-trap quantum simulator, which implements the corresponding Landau-Zener crossings governing the dynamics in each mode. The Rabi frequency Ω R = 4J/ was set to 2π × 20 kHz in experiment. For each value of k, the excitation probability is estimated from 10000 measurements. The shaded region describes the excitation probability predicted by the Landau-Zener formula. Deviations from the latter become apparent for fast quench times, especially for large values of k.
Full probability distribution of topological defects. From the experimental data, the distribution P (n) of the number of kink pairs is obtained using the characteristic function (S11) of a Poisson binomial distribution in which the probability p k of each Bernoulli trial is set by the experimental value of the excitation probability upon completion of the corresponding Landau-Zener sweep. This result is compared with the theoretical prediction valid in the scaling limit -the regime of validity of the KZM-in which P (n) approaches the normal distribution, away from the adiabatic limit [36] P (n) 1 A comparison between P (n) and Equation (7) is shown in Fig. 1. The matching between theory and experiment is optimal in the scaling limit far away from the onset of the adiabatic dynamics or fast quenches, for which non-universal corrections are expected. We further note that the onset of adiabatic dynamics enhances non-normal features of the experimental P (n) that cannot be simply accounted for by a truncated Gaussian distribution, that takes into account the fact that the number of kinks n = 0, 1, 2, 3, . . . The experimental kink-pair number distribution for a TFQIM with N = 100 spins is compared with the Gaussian approximation derived in the scaling limit, Eq. (7), ignoring high-order cumulants. The mean and width of the distribution are reduced as the quench time is increased. The experimental P (n) is always broader and shifted to higher kinkpair numbers than the theoretical prediction. Non-normal features of P (n) are enhanced near the sudden-quench limit and at the onset of adiabatic dynamics.
To explore universal features in P (n), we shall be concerned with the scaling of its cumulants κ q (q = 1, 2, 3 . . . ) as a function of the quench time. We focus on the first three, theoretically derived in the scaling limit in [43]. The first one is set by the KZM estimate for the mean number of kink pairs κ 1 (τ Q ) = n KZM . The second one equals the variance and as shown in the SM is set by Figure 4: Universal scaling of the cumulants κ q (τ Q ) of the kink-pair number distribution P (n) as a function of the quench time. The experimental data (symbols) is compared with the scaling prediction (dashed lines) and the numerical data (solid lines) for the first three cumulants with q = 1, 2, 3. The universal scaling of the first cumulant κ 1 (τ Q ) = n is predicted by the KZM according to Eq.(4). Higher-order cumulants are also predicted to exhibit a universal scaling with the quench time. All cumulants of the experimental P (n) exhibit deviations from the universal scaling at long quench times consistent with a dephasing-induced anti-KZM behavior. Further deviations from the scaling behavior are observed at fast quench times. The range of quench times characterized by universal behavior is reduced for high-order cumulants as q increases. The error bars are much smaller than the size of the symbols used to depict the measured points. being dictated by the KZM. We compared this theoretical prediction with the cumulants of the experimental P (n) in Fig. 4, represented by dashed lines and symbols, respectively. As discussed in the Supplementary Information, deviations from the power law occur at very fast quench times, satisfying τ Q < /(π 3 J). However, this regime is not probed in the experiment, as all data points are taken for larger values of τ Q . Within the parameters explored, deviations at fast quenches with Jτ Q / ∼ 1 are due to the fact that the excitation probability p k in each mode is not accurately described by the Landau-Zener formula, as shown in Fig. 2. For moderate ramps, the power-law scaling of the cumulants is verified for q = 1, 2. The power-law scaling for κ q (τ Q ) with q > 1 establishes the universal character of critical dynamics beyond the KZM. The latter is explicitly verified for κ 2 (τ Q ), as shown in Fig. 4. Beyond the fast-quench deviations shared by the first cumulant, the power law scaling of the variance of the number of kink pairs κ 2 (τ Q ) extends to all the larger values of the quench time explored, with barely noticeable deviations. By contrast, for the third cumulant κ 3 (τ Q ) the experimental data is already dominated by non-universal contributions both at short quench times, away from the scaling limit. For slow ramps, the third cumulant exhibits an onset of adiabatic dynamics due to finite size effects around Jτ Q / = 10 2 . Finitesize effects are predicted to lead to a sharp suppression of the third cumulant. However, the experimental data shows that κ 3 starts to increase with τ Q . This is reminiscent of the anti-Kibble-Zurek behavior that has been reported in the literature for the first cumulant in the presence of heating sources [44,45] and we attribute it to noise-induced dephasing of the trapped-ion qubit. The nature of these deviations is significantly more pronounced in high-order cumulants, reducing the regime of applicability of the universal scaling.

Discussion
In summary, using a trapped-ion quantum simulator we have probed the full counting statistics of topological defects in the quantum Ising chain, the paradigmatic model of quantum phase transitions. The statistics of the number of kink pairs has been shown to be described by the Poisson binomial distribution, with cumulants obeying a universal power-law with the quench time in which the phase transition is crossed. Our findings demonstrated that the scaling theory associated with a critical point rules the formation of topological defects beyond the scope of the Kibble-Zurek mechanism, which is restricted to the average number. Our work could be extended to probe systems with topological order in which defect formation has been predicted to be anomalous [46]. We anticipate that the universal features of the full counting statistics of topological defects may be used in the error analysis of adiabatic quantum annealers, where the Kibble-Zurek mechanism already provides useful heuristics [10].

Methods
A detailed description of analytical methods and experimental techniques employed to obtain the results can be found in the Supplementary Information accompanying this work.

Data availability
The data that support the plots and other findings within this paper are available from the corresponding author upon reasonable request.

-Supplementary Notes-Experimentally testing quantum critical dynamics beyond the Kibble-Zurek mechanism
Inserting (S2) into the Ising Hamiltonian Eq. (S1) we find that the Hamiltonian is quadratic in the Fermi operatorŝ This Hamiltonian includes terms such asĉ mĉm+1 orĉ † m+1ĉ † m that violate the fermion-number conservation. However, the parity operatorΠ = (−1) N m=1ĉ † mĉm is a constant of motion having the value +1 and −1. We expandĉ m into the Fourier modeŝ where the pseudomomenta k take values Using the choice of even N and considering anti-periodic boundary conditionĉ N +1 = −ĉ N , the Hamiltonian maps toĤ whereĤ k is defined byĤ k ≡ 2J (g − cos (ak))σ z + 2J sin (ak)σ x . (S7) Next, we use the Bogoliubov transformation to rewrite the Hamiltonian in terms of a new set of fermionic operatorsγ k whose number is conserved. These new operators are defined via a unitary transformation in terms of the pairĉ k ,ĉ † −k : where u k , v k are numbers satisfying u 2 k + v 2 k = 1, u −k = u k , and v −k = −v k . It can be checked that canonical fermion anti-commutation relation for theĉ k imply that the same relations are also satisfied byγ k , that is, The Hamiltonian is diagonalized in terms ofγ k , where k (g) = 2J (g − cos (ak)) 2 + sin 2 (ak), with the ground state definedγ k |0 = 0. The ground state free energy reads Supplementary Note 2 Cumulants of the kink-pair number distribution P (n) The kink-pair number distribution P (n) describes the probability of observing n pairs of kinks. The number of kink pairs is measured by the observablê with eigenvalues n = 0, 1, 2, . . . The kink-pair number distribution is thus defined as the quantum expectation value of the projector operator δ[N − n] on the subspace with n pairs of kinks where the expectation value is taken over the final stateρ upon completion of the quench. The Fourier transform of P (n) is the characteristic function The cumulants κ q of the distribution P (n) are defined via the expansion of the cumulant-generating function, which is the logarithm of the characteristic function, We next focus on the characterization of P (n) for the state that results from the nonadiabatic crossing of the phase transition in the TFQIM and derive the exact expression for the first few cumulants {κ q }. To this end, we note that the nonadiabatic crossing of the phase transition results in the production of quasi-particles (and kinks) in pairs. Specifically, due to the conservation of momentum, quasiparticles with wavectors +k and −k are excited jointly. The distribution of the number of pairs of kinks is a Poisson binomial distribution associated with N/2 Bernoulli trials each with probability p k for k > 0. Its characteristic function reads with k ∈ {π/N, 3π/N, . . . , (N − 1)π/N } for an Ising chain with period boundary conditions and p k denoting the excitation probability in the k-mode. The Poisson binomial distribution is well studied and its cumulants are known to be of the form Their explicit evaluation can be performed using the Landau-Zener formula for p k . In the continuum limit, this yields In terms of the error function, it is found that an expression that is valid even for moderately fast quenches, as long as described by the Landau-Zener formula but possibly away from the scaling limit. However, for large τ Q > /(2π 3 J), erf √ π 2d approaches unity and the average number of kink pair simply reads in agreement with the KZM [2]. A similar derivation shows that the variance is given by This results in the exact expression that in the scaling limit reduces to To characterize the non-normal features of the distribution of kinks, we also provide the exact expression for the third cumulant κ 3 , that equals the third centered moment µ 3 = (n − n ) 3 . In particular, we note that κ 3 is related to the skewness γ 1 of the distribution, as γ 1 = κ 2 /κ 3/2 2 . In the continuum, it reads Supplementary Figure. 1: Universal scaling of the cumulants κ q (τ Q ) of the kink-pair number distribution P (n) as a function of the quench time. The solid lines show the exact result for the first three cumulant κ q (τ Q ) with q ∈ 1, 2, 3, as predicted by equations (S17), (S19), and (S21). The dashed lines show the universal scaling of the first cumulant according to Eqs. (S18), (S20), and (S22). In the experimental regime, marked by the shaded area, the scaling behavior holds.
This expression is simplified in the scaling limit, performing a power series expansion with τ Q → ∞.
To leading order, we obtain the expression quoted in the main text Note that the scaling expressions Eqs. (S18), (S20) and (S22) are obtained from the more general ones, Eqs. (S17), (S19) and (S21), whenever the error functions can be replaced by unity. The error function increases monotonically and swiftly from zero value to unity as the argument is increased from zero value. Indeed, it saturates at unity for fairly moderate values of the argument, e.g., erf(2) = 0.995. By imposing 2qπ 3 Jτ Q > 2, we predict the onset of the scaling limit for κ q at This condition is satisfied in all the experimental realizations we have performed. As a result, it is justified to consider the scaling limit, in which it is actually possible to derived a closed form expression of the characteristic function [3] in terms of the polylogarithmic function Li 3/2 (x) = ∞ p=1 x p /p 3/2 . From this expression, it follows that all cumulants are actually nonzero and proportional to n KZM . The kink-pair number distribution is manifestly non-normal. The Gaussian (normal) approximation discussed in the main text follows from neglecting all κ q with q > 2. One can also approximate 1 − 1/ √ 2 ≈ 3/π 2 , so that and thus P (n) 1 The experimental data shows deviations from the scaling regime at short and long quench times that become more and more pronounced with increasing order q of the cumulant κ q . Deviations at long times are attributed to dephasing-induced heating and result in an anti-KZM behavior whereby the value of κ q increases with τ Q . Deviations at short times are of a different nature: For Jτ Q ∼ 1 deviations from the Landau-Zener formula occur, as shown in Fig. 2 of the main text.
For the sake of clarity we elaborate on the connection between the distribution of kink pairs and the distribution of the total number of kinks. To account for the latter, we note that its characteristic function is simply obtained as where the factor of 2θ (instead of the θ) accounts for the creation of kinks in pairs. Cumulants κ T q of the total kink number distribution are simply related to cumulants κ q of the distribution of pairs by the identity Figure. 2: Scheme to measure the excitation probability. The quantum critical dynamics of the one-dimensional TFQIM is detected by measuring corresponding Landau-Zener crossings governing the dynamics in each mode. For each mode, a typical process to measure the excitation probability in three stages is shown in a, b and c.
frequency of the qubit simulator is set around 20 kHz, which depends on the driven power of our microwave amplifier. Higher microwave power can shorten the Rabi time T R = 1/Ω R , and reduce the total operation time. To simulate the quench process, one would like to vary g from −∞ to 0, an idealized evolution that needs infinite time, which can not be realized in experiment. However, to explore the universality associated with the crossing of the critical point one can initialize the system sufficiently deep in the paramagnetic phase. According to the KZM, it is sufficient to choose an initial value of the magnetic field such that the corresponding equilibrium relaxation time is much smaller than the time left until crossing the critical point. The system is then prepared out of the "frozen region" in the language of the adiabatic-impulse approximation [8]. For the quench time τ Q ≥ 1, the initial value of g = −5 is far out of the frozen time. We can simulate the TFQIM with an initial g = −5 and no excitations, by preparing the initial state of the qubit in the ground state ofĤ s k,i =Ĥ s k (g = −5) before the quench process. The initial state can be derived by solving the eigenvectors of Eq. (S30), which is where θ k,i = − arctan sin(ka) 5+cos(ka) . The scheme to detected the quantum critical dynamics of the one-dimensional TFQIM in each mode by using a single qubit is shown in Fig. 2. The whole process can be divided to three steps. Before the process, the qubit has been pumped to |0 state by using a 369nm laser to excite transition S 1/2 , F = 1 → P 1/2 , F = 1. In the first stage, the ion-trap qubit is prepared into the state |ψ − k,i by a resonant microwave pulse. The second stage is the quench process; the Hamiltonian is time dependent and varies fromĤ s k,i ≡Ĥ s k (g = −5) toĤ s k,f ≡Ĥ s k (g = 0). This is implemented by driving a chirped pulse into the qubit. The chirped pulse is in the form of Eq. (S31) with g(t) = 5(t/T p − 1), where T p = 5Jτ Q sin(ka)/ ( Ω R ) is the pulse length. The third stage is to measure the excitation probability after the quench, which is to measure the occupation probability on the excited eigenstate ofĤ s k,f . The excited eigenstate is |ψ − k,f = sin θ k,f |0 +cos θ k,f |1 , where θ k,f = −ka. As the fluorescence detection on 171 Yb + can only discriminate |0 and |1 sates, we need to rotate the state |ψ − k,f to |1 and then detect the bright state probability. We thus use a qubit rotation and a fluorescence detection in this stage. The calculated waveform to set the AWG to perform the pre-rotation, quench and post-rotation in the three stages will be discussed in the next section.
Driving Waveform with an Arbitrary Wave Generator: We consider the driving microwave from AWG is A cos(ω c t + Φ(t)), where A is the amplitude, ω c is the carrier frequency and Φ(t) is an arbitrary phase function. The carrier frequency can be mixed up with some frequency ω 0 from a local oscillator (LO) by using a frequency mixer. The waveform at the output of the mixer is A cos(ω c t + Φ(t)) cos(ω o t) = A/2 (cos((ω o + ω c )t + Φ(t)) + cos((ω o − ω c )t + Φ(t))) .
When the wave passes through the high pass filter, the wave form is filtered as A/2 cos((ω o + ω c )t + Φ(t)). In the experiment, we choose the frequency ω o + ω c resonant with the qubit transition ω Q = ω o + ω c , so the magnetic field of the final driving waveform is B(t) = B 0 cos(ω Q t + Φ(t)). The interaction between the magnetic field and spin is H I = −µ · B, where µ is the magnetic dipole of the spin qubit. The Hamiltonian can be further simplified after the rotating wave approximation in the interaction frameĤ where the Rabi frequency Ω R = − 0| µ · B |1 / . The phase function Φ(t) corresponds to the azimuthal angle in the Bloch sphere. Equation (S30) can be also expressed in interaction frame H s I,k = 2 Ω R |1 0| exp i ∆(t)dt + Ω * R |0 1| exp −i ∆(t)dt , so in the quench stage, the phase function is Φ(t) = T P 0 ∆(t)dt = Ω R sin ka Ω R τ Q sin ka t 2 − (5 + cos ka) t .
We rotate the state along a vector in the equatorial plane to prepare a state from |0 , or measure state to |1 , which means Φ(t) is constant in these two stages. The pulse length for the preparation and measurement stages is determined by the polar angle of the ground state at the beginning and end of the quench process, t = θ/Ω R , respectively. The whole expression of Φ(t) in the three stages can be derived as where t 1 = (2π − θ k,i )/Ω R , t 2 = t 1 + T p , t 3 = t 2 + (2π + θ k,f )/Ω R and φ f = −5τ Q (2.5 + cos ka).
Qubit preparation and detection error: There are some limitations in the preparation and measurement of the qubit. We measured the error through a preparation and detection experiment. First, we prepare the qubit in the |0 state by optical pumping method, and detect the ion fluorescence. Ideally, no photon can be detected as the ion is in the the dark state. However,