Robust Macroscopic Schr¨odinger’s Cat on a Nucleus

We propose a scheme to generate spin cat states, i

In continuous-variable systems, cat states are generated using nonlinearity induced by light-matter interactions in the Jaynes-Cummings model [21,22], e.g., a harmonic oscillator coupled to a superconducting qubit [23][24][25] or an atom in an optical cavity [26].However, errors occurring in the qubit can propagate to the cavity, corrupting such nonclassical states [25].Alternatively, Greenberger-Horne-Zeilinger (GHZ) states, which are spin cat states (SCSs) for permutationally symmetric systems, can be realized by entangling multiple qubits [27][28][29][30], but such SCSs are fragile-with their decoherence rate proportional to the square of the number of qubits [31].
We propose a new platform to generate cat states deterministically, on-demand and without post selection using a nucleus, which is controlled without entanglement with an ancilla, hence, is robust against error propagation, unlike continuous-variable systems [12,[23][24][25][26].
Here, we consider SCSs (implying macroscopic to differentiate them from "kitten" states [32]) on a nuclear spin, e.g., shown in Fig. 1, but, more generally, a superposition of any two maximally separated quasi-classical states, which are minimum uncertainty states that follow a classical evolution under a linear drive [5,6].We are motivated by recent experiments on high-spin nuclei (spin I > 1) in a solid-state device [33,34] accompanied by an enhanced electric field-driven quadrupole interaction, two orders of magnitude faster than the dephasing timescale, that could achieve ultrafast non-linear dynamics.
A highly coherent nuclear spin could open the possibility of robust generation and long-lived storage of SCSs.So far, quadrupole interaction has only been utilized for realizing spin squeezed states (SQSs), that too, on an ensemble of nuclear spins [35,36].However, individual addressing of high-spin systems is important for applications like quantum information processing and error correction [15].We present two key results: i) a scheme to go beyond squeezed states towards generating SCSs and ii) an experimentally feasible control protocol for storing and manipulating SCSs on an individual nucleus.

II. SETUP
We consider a high-spin donor (Group-V element with a high-spin nucleus) implanted in a nanoelectronic device based on enriched 28 Si (or any other spinless Group-IV atoms), as described in [33].The chip is affixed to a dilution refrigerator maintained at a low temperature (T ≈ 0.01 K) and placed inside a superconducting coil that produces a static magnetic field B 0 ≈ 1 T. As shown in Fig. 2(a), the chip includes electrostatic gates to control the electrochemical potential of the donor, a microwave antenna to supply electromagnetic pulses, and a single electron transistor (SET) for reading out the donor spin state [37].
An implanted Group-V donor contains five electrons in its outermost shell, four of which form bonds with neighboring Si atoms and the remaining electron can be removed by ionizing the donor using electrostatic gates.The quantum state of the nuclear spin is initialized using flip-flop drive [38], where an electron and nuclear spin pair is flipped together in opposite directions.The donor atom is subject to a static potential, and a target nuclearspin eigenstate can be achieved by applying a combination of flip-flop pulses and electron spin-resonance pulses that move population from higher and lower nuclear spin states towards a target state [34].Then, the electrostatic potential is increased to ionize the donor atom by removing the electron and only the nuclear spin is manipulated using RF pulses.[39].The readout process using electron spin resonance can be made dependent on the nuclear spin state, thus achieving nuclear readout via spin dependent tunneling of the electron [33,40].
The differential thermal expansion between the silicon substrate and the electrostatic gates induce lattice strain in the chip [33], thereby distorting the bonds between the donor atom and neighboring Si atoms, as shown in Fig. 2(b).The electric-field gradient (EFG) generated from distortion of the bonds results in a quadrupole interaction of the high-spin nucleus, which has a nonzero quadrupole moment due to its asymmetric nuclear charge distribution.For a spatially varying potential V (x, y, z), the EFG tensor comprises elements V αβ := , α, β ∈ {x, y, z}, which is real, traceless and symmetric.Thus, V αβ can be diagonalized in a principal axis system {x ′ , y ′ , z ′ }, such that and the EFG is characterized by the asymmetry parameter η := for 0 ≤ η ≤ 1.For a nucleus with electric quadrupole moment q n , the quadrupole interaction is where ω q is the quadrupole coupling strength, e the elementary charge, ℏ the reduced Planck's constant, the nuclear spin operator vector Î := [ Îx ′ , Îy ′ , Îz ′ ] and invariant I 2 = I • I.
In addition to the quadrupole interaction, we manipulate the nuclear spin state using a static magnetic dipole interaction due to B 0 and time dependent perturbation driven by magnetic field B 1 with time-dependence ϵ(t) (|ϵ(t)| ≤ 1).For a nucleus with gyromagnetic ratio γ, the nuclear spin dynamics are described by The 2I+1 degenerate levels of a nuclear spin undergo Zeeman splitting due to γ Î • B 0 which creates an equal spacing between the energy levels.The electric quadrupole interaction Ĥq causes an unequal shift of the levels such that the transition frequency ω i,i−1 between two consecutive levels i and i − 1, for i ∈ [I, −I + 1], is

III. ONE-AXIS TWISTING
Quadrupole interaction is equivalent to one-axis twisting (OAT) due to Î2 z ′ if η = 0, two-axis counter twisting (TACT) due to Î2 z ′ and Î2 y ′ if η = 1 [36,42], and partial TACT otherwise.Generally, TACT is preferred over OAT as TACT generates faster squeezing in comparison to OAT and can help produce SQSs for metrology [43].However, TACT cannot produce SCSs, and, in general, quadrupolar non-linear evolution with η > 0 cannot go beyond SQSs on a nuclear spin [35].
We devise a way to convert TACT to OAT by adding a large linear term to the quadrupole interaction, induced here by the Zeeman coupling.The resulting dynamics depend on the orientation of the EFG and strength of magnetic field B 0 .When the Zeeman splitting is much larger than quadrupole coupling, i.e. γB 0 ≫ ω q , the component perpendicular to B 0 does not affect the nuclear spin because quadrupolar transitions (∆m I = ±2) are /2ω q /ω q 3/2ω q 2/ω q FIG. 3. Degree of superposition N r eff under free evolution (ϵ(t) = 0) of a coherent spin state to a SCS at t = π /2ωq.For differing nuclear spins, the maximum value N r eff ∝ 2I, showing the macroscopicity of these SCSs, but the I = 3 /2 case lacks a "hill-on-mesa" and is actually a "spin 'kitten' state".
prevented by the large energy difference.In this regime, only the component of the quadrupole interaction parallel to B 0 generates non-linear rotation, which effectively results in OAT, regardless of the value of η.In contrast to TACT, our scheme allows generation of SCSs for any orientation of the quadrupolar axes of symmetry {x ′ , y ′ , z ′ } with respect to linear interaction along the z axis.For numerical simulations in this work, we set B 0 = B 0 ẑ, B 1 = B 1 ŷ, η = 0, and EFG to be symmetric about z axis, i.e. ẑ′ = ẑ, but, our results can be generalized to other configurations.
Despite the femtometer scale of a nucleus, nuclear-SCSs are macroscopic with their size, quantified by the relative quantum Fisher information [44,45], scaling linearly with I.The size or degree of superposition N r eff of a nuclear-spin state ψ for measurement operator Ô is i.e., proportional to the variance of the observable.For an appropriate Ô, the variance is I 2 for bimodal SCSs and we calculate the size of a nuclear SCS to be 2I, which denotes the maximum quantum Fisher information of the system and is equivalent to a GHZ state [46] made of N = 2I entangled qubits [47].A nuclear SCS with 2I + 1 dimensions, e.g., a spin-7 /2 SCS with 8 dimensions is equivalent to a 2I = 7 qubit GHZ state with permutation symmetry, and not a log 2 (2I + 1) = 3 qubit state, and can have macroscopicity up to 2I = 9 for a 209 Bi donor.The degree of superposition N r eff (Eq.( 3)) increases and decreases periodically during OAT, as shown in Fig. 3. N r eff = 1 for a spin coherent state on the Îx axis of the Bloch sphere at t = 0, increases as the spin state squeezes, then reaches the flat-top of the "mesa" where it interferes with itself and finally increases again ("hill") as two distinct components form resulting in a SCS at t = π /2ωq.Under further evolution, the SCS returns to a coherent state at t = π /ωq, which again squeezes to a SCS at t = 3π /2ωq, and so on.This periodic process of collapse-and-revival of the SCS has a time period of π /ωq and is a signature of coherence [47].From the plots for different values of nuclear spin in Fig. 3, we note the macroscopicity of the SCSs through the linear scaling of their size (N r eff ) max = 2I.

IV. CONTROL SCHEME
External control of a high-spin nucleus is required for initialization, manipulation and detection of nuclear-SCSs.For initialization, the eigenstate |I, I⟩, which is also a coherent state, should first be rotated to reside on the Îx axis of the Bloch sphere.Then, OAT results in a SCS on the Ô = Îy axis for a half-integer nuclear spin, but, for donors in solid-state devices, nuclear-spin initialization and measurement is limited to the eigenstates of Îz .Observing nuclear SCSs could be possible with angular momentum operations, which are trivial when γB 1 ≫ ω q [42,48] or with a tunable nonlinearity [49].However, quantum control is needed as, experimentally, γB 1 = 2π × 800 Hz ≪ ω q = 2π × 40 kHz such that a RF pulse can only selectively drive a particular nuclear spin transition leading to Givens rotations, and the quadrupole coupling cannot be turned on-andoff.
For applying global SU(2) rotations with an alwayson non-linearity and low drive strength, we use a multitone (MT) pulse that is a superposition of 2I tones of equal amplitudes, corresponding to the 2I nearest neighbor sub-spaces with transition frequencies represented by frequency vector ω := {ω i,i−1 } −I−1 i=I .Each tone has a global phase ϕ and a pulse envelope switched on at a time t 0 is described by the amplitude ϵ MT (t; ω, ϕ, t 0 ) = 1 2I ωj ∈ω cos(ω j (t − t 0 ) + ϕ).Here, ϕ tunes the axis of rotation on the Bloch sphere while preserving the relative phase between different energy levels.The frequency Ω of SU(2) rotations is generalized from the Rabi frequency γB1 /2 for a two-level nuclear spin, to a high-spin nucleus by dividing the Rabi frequency by the number of tones 2I, i.e.Ω = γB1 /4I.Applying this pulse for time ∆t rotates the nuclear spin by angle Θ given by Θ = Ω∆t.In contrast to other schemes that employ numerical techniques with multi-tone control [50,51], we derive exact expressions for phase-locked pulses, which can enable novel experimental applications such as fullstate tomography and interferometry with a high-spin nucleus.
Figure 4(a) shows our Ramsey-like control sequence to realize SCSs, where a square-shaped MT RF pulse is applied to the initial state |I, I⟩ between time t 0 = 0 and t 1 = tπ /2 .Then the pulse is switched off for time T , and again applied between t 2 = T + tπ /2 and t 3 = T + 2tπ /2 , with the amplitude ϵ cat (t; ω, ∆ϕ, T ) =ϵ MT (t; ω, 0, t 0 ) ⊓ (t; t 0 , t 1 ) + ϵ MT (t; ω, ∆ϕ, t 2 ) ⊓ (t; t 2 , t 3 ) (4) for ⊓ the top-hat function.Here, the global rotation of Θ = π /2 would take time ∆tπ /2 = π 2 • 4I γB1 and setting the phase difference ∆ϕ = π /2 between the first and the second global rotations accounts for the angle between a coherent state and the resulting SCS on the equator.The resulting dynamics are shown in Fig. 4(d), where we note that SCSs along any axis can be realized using global rotations (Fig. 4(d)(vi)).The quadrupole coupling cannot be turned off and these states squeeze during rotations.A SCS can be stored for long times by aligning it to the fixed points on the z axis (Fig. 4

(d)(vii)).
For a spin-7 /2 123 Sb nucleus, Fig. 4(b) shows the change in N r eff , for Ô = Îz , as T is varied, similar to Fig. 3.When ∆ϕ = π /2, we observe a rapid variation of N r eff with a frequency of γB0 /π, because nuclear-spin states generated by OAT undergo Larmor precession and when the axis of rotation corresponds with the axis of the SCS, the second global rotation is unable to move the state.To remedy this immobility, we assign a T -dependent shift to the phase difference: ∆ϕ = π /2 + ωT in (4), which is equivalent to a rotating-frame transformation, and results in the periodic collapses-and-revivals of a SCS solely due to OAT, which would enable efficient experimental detection.

V. VIRTUAL PHASE CONTROL
In this section, we show that cat-state generation is possible with just phase modulation of multi-tone pulses in the generalized rotating frame of reference.First, we note that, in the limit of large Zeeman interaction γB 0 ≫ ω q , we can neglect non-commuting terms, such as Î2 x ′ of the quadrupole interaction (1) so that the nonlinear evolution is only described in terms of the strength of quadrupole tensor, ω eff q , parallel to the magnetic field B 0 .Under this approximation and using the notation [ℓ, ℓ ′ ] := {ℓ, ℓ + 1, . . ., ℓ ′ }, the effective Hamiltonian is ϵ j cos(ω j t + ϕ j ), (5) where ω j , ϵ j , and ϕ j are respectively the frequency, amplitude and phase of the j th component of the multi-tone control.The eigenstates |I, k⟩ of the nuclear spin have energy for k ∈ [−I, I].Now we transform the Hamiltonian to the generalized rotating frame, defined by the diagonal matrix [50] U rot = diag e −ie I t , e −ie I−1 t , . . .
and obtain where We note that global rotations are applied by setting ϵ j = 1 /2I and ϕ j = ϕ; for example, using ϕ = 0 is equivalent to Ĥrot = γB 1 Îx 2I (10) and ϕ = π /2 is equivalent to Ĥrot = γB1 Îy /2I.One-axis twisting can be realized by virtually updating the phases of a multi-tone pulse.To see this, we note that that the effective quadrupolar operator applied for duration T generates the operation e −iT ω eff q Î2 z , which results in a non-linear phase shift between the different eigenstates of a nuclear-spin Hamiltonian.Using the dummy variable k for integers in the interval [−I, I], a state is transformed as where −T ω eff q k 2 is the resulting phase shift on the energy eigenstate |I, k⟩.Now, we can factor out the global phase −T ω eff q I 2 and instead of freely evolving a state under one-axis twisting, apply equivalent phase updates, to the subsequent multi-tone pulse.Such virtual phase updates would have same effect as free non-linear evolution and would result in cat state formation when T = π /2ω q eff .Thus, in the generalized rotating frame, we can start with the initial state |I, I⟩, apply a global π /2 rotation with ϕ j = ϕ, followed by a second by a global π /2 rotation with realizing cat states by only phase modulation of multitone pulses.

VI. DECOHERENCE
Cat states with a lifetime of several milliseconds can be achieved using silicon-based chips with an implanted spin-7 /2 123 Sb donor atom, that have a nuclear spin coherence time T * 2 ≈ 100 ms for two-level transitions inferred from experiments on nuclear electric resonance [33].SCSs along the z-axis would have a reduced coherence time, lower by a factor of 2I, due to the large separation between the two components, resulting in a coherence time of 100 /7 ≈ 14 ms for 123 Sb.Even for a modest driving strength γB 1 = 2π × 800 Hz, a π /2-global rotation would take ∆tπ /2 = π 2 • 4I γB1 = 4.375 ms.Thus, our approach, using two MT pulses and OAT, can yield highly coherent SCSs in less than 9 ms-within the dephasing time on current hardware.Faster Rabi frequencies are easily achieved for nuclear magnetic resonance (NMR) methods [34,40], that could reduce the time needed for MT pulses and lead to faster generation and detection of nuclear SCSs.
We compare the performance of OAT to an intuitively appealing gate sequence for creating SCSs, achieved by sequential Givens rotations [52].Starting with |I, I⟩, we apply a π /2-pulse of frequency ω I,I−1 to form the superposition |I, I⟩ + |I, I − 1⟩, with unit norm implied.Then, we apply a sequence of 2I − 1 π-pulses to create the SCS |I, I⟩ + |I, −I⟩.For observing a collapse of the SCS, repeating the same sequence evolves the SCS into the coherent state |I, −I⟩, opposite to the initial state |I, I⟩.In total, 4I pulses would be needed for each collapse-and-revival using Givens rotations, that would take ∼ 9 ms for a I = 7 /2 with Rabi frequency 2π×800 Hz, in contrast to OAT which takes just π /ωq = 12.5µs when ω q = 2π × 40 kHz.
We simulate our scheme under decoherence from a magnetic field fluctuating at rate Γ m , modeled using Lindblad operator Lm = Îz , and an electric field fluctuating at rate Γ e , modeled using Le = Î2 z .Figure 4(c) shows N r eff vs evolution time beyond 9 ms (approximate time to form the first SCS).We note that OAT is two orders of magnitude faster than both Givens rotations and the dephasing timescale, opening the possibility of realizing hundreds of SCS collapses-and-revivals by using a strong non-linearity relative to decoherence, which makes our SCS robust, in contrast to alternative platforms [53].Additionally, SCSs generated by Givens rotations lie on the Îz axis, where the dephasing rate of the cat state increases with the nuclear spin increases; in comparison, cat states produced by one-axis twisting lie on the equator, perpendicular to the dephasing axis, and are highly protected against decoherence.

VII. CONCLUSION
In summary, Schrödinger cat states are important for investigating foundations of quantum mechanics as well as for applications such as quantum error correction but are fragile and hard to realize.We propose a new platform for realizing macroscopic spin cat states by converting quadrupolar two-axis counter twisting to one-axis twisting on a high-spin donor in a solid-state device.Our scheme can deterministically generate spin cat states ondemand without post-selection and is two orders of magnitude faster than the dephasing timescale, resulting in highly coherent states.We show that, in the generalized rotating frame, a spin coherent state can be converted to a cat states by only phase modulation of multi-tone pulses.Unlike continuous-variable systems, we generate cat states without entanglement with an ancilla, which renders robustness against error propagation [25] and opens the possibility of storing long-lived nuclear-spin cat states.Furthermore, our cat states can be experimentally realized on a single nucleus, which could lead to novel applications in quantum information processing.FIG. 5. Husimi Q-function plots of SCSs for different values of the spin.The spin-1 case lacks bi-modality as the spin coherent states overlap with each other.The spin-3 /2 case corresponds to a "spin 'kitten' state" with small overlap between its two components.Higher spin (≥ 5 /2) form macroscopic cat states.
to reach a target state in the energy ladder, first, the nuclear spin is pumped down from all higher-lying energy levels using flip-flop transitions (orange).Then, the nuclear spin state is pumped up from all the lower-lying energy levels using a combination electron-spin resonance pulses (green) and flip-flop pulses (orange).
The initialization and readout schemes benefit from the coupling of a nuclear spin to an electron spin, which enables fast initialization of the nuclear spin to any energy eigenstate |I, m I ⟩, of which |I, ±I⟩ are spin coherent states.We use the ground state |I, I⟩ as the initial state in our scheme.After initialization, the electrostatic potential is increased to ionize the donor atom by removing the electron and only the nuclear spin is manipulated using RF pulses [39].When the nuclear spin reaches the desired final state, the donor atom is de-ionized such that the nuclear spin again couples to an electron.As shown in Fig. 6 (b), the nuclear-spin state is read out using single-shot detection with pulsed electron spin resonance and spin-dependent tunneling of the electron to the SET, which acts as a cold-electron reservoir [54].

Control
We manipulate the nucleus using RF magnetic pulses that cause NMR transitions with selection rule ∆m I = ±1.The quadrupole interaction shifts the energy levels of a nuclear spin qudit, as shown in Fig. 2 of the main text, which allows selectively addressed transitions by tuning the frequency of Rabi pulses, applied through the timedependent perturbation ϵ(t) in Eq. (A1).At off-resonant drives, the probability of a Rabi flip is given by, where γB 1 is the drive strength (A1), Ω 2 = (γB 1 ) 2 + ∆ω 2 /4 with ∆ω the difference between the detuning of the drive from the transition frequency.With quadrupole interaction ω q ≈ 2π × 40 kHz achieved in experiments [33,34] and drive strength γB 1 = 2π×0.8kHz, frequency difference between two transitions ∆ω ≥ 2ω q , thus, P flip < 0.0004.Thus, cross talks are highly suppressed.Additionally, the quadrupole interaction generates fast collapse and revivals due to non-linear evolution having time-period of a few microseconds, which can be resolved by pulses applied using an Arbitrary Waveform Generator (AWG) with a resolution of 1 ns.While lower-order non-linearity would slow down SCS generation, higher-order non-linearity, for example in the GHz range, would generate OAT at a rate faster than the resolution of an AWG.
Appendix C: Decoherence

Master equation for a high-spin nucleus
We model decoherence using a Lindblad master equation of the form where ρ is the density matrix of the nuclear spin, Γ m is the decay rate due to magnetic field fluctuations, modeled using Lindblad operator Lm = Îz and Γ e is the decay rate due to electric field fluctuations, modeled using Lindblad operator Le = Î2 z .Here, we take the electric field fluctuation only in the z direction, i.e. parallel to the FIG. 7. Collapse and revivals under decoherence (C1).The line joining the revival points (where N r eff is maximum) can be extrapolated to estimate the number of revivals possible within the coherence time of a nuclear spin.The decoherence rates are specified in units of kHz and time in ms.
Zeeman splitting, because other components of the electric quadrupole interaction are washed out, as described above.
We simulate the cat-state collapse-and-revival control scheme for various values of Γ m and Γ e , shown in Fig. 7. Here, magnetic fluctuations slowly dephase the nuclear spin and the coherence of the SCS decreases during collapse-and-revivals.On the other hand, electric field fluctuations have a more drastic effect on SCSs generated by quadrupole interaction, resulting in an overall decrease of the variance, i.e. macroscopicity of the SCSs.

Coherence time of a qudit SCS
We study the scaling of the coherence of a SCS with respect to the Hilbert space dimension of a qudit [55] and consider dephasing purely due to magnetic field fluctuations, of the state |I, I⟩ + |I, −I⟩ for different values of I.The density matrix of such a state consists only of four elements: two diagonal elements corresponding to the population in |I, ±I⟩ and two off-diagonal components that quantify the coherence of the superposition state.In Fig. 8, we plot the magnitude of the off-diagonal component after a fixed time of 1 ms under an amplified dephasing at rate of Γ m = 1 kHz (Γ m ∝ 1/T * 2 ).We note that coherence decreases with the dimension d = 2I + 1 of the nuclear spin and results in a linear decrease in coherence time as the dimension of a qudit increases.Such a trade-off between lifetime and the macroscopic- ity of a SCS is similarly found in other systems such as harmonic oscillators and maximally entangled GHZ states [56].However, compared to other systems, a nuclear spin, with naturally long coherence times, can help increase the lifetime of practically achievable SCSs.

FIG. 2
FIG. 2. a) High-spin donor implanted in a fabricated siliconbased chip.b) Bond distortion due to lattice strain leads to an electric-field gradient at the donor site.c) Energy-level diagram of a high-spin nucleus.

FIG. 6 .
FIG. 6. Initialization (a) and readout (b) of a nuclear spin by coupling it to an electron spin and using flip-flop pulses (orange arrows) and electron-spin resonance pulses (green arrows).

FIG. 8 .
FIG.8.Decrease in coherence, given by the off diagonal component of the density matrix, of a SCS with increase in qudit dimension d = 2I + 1.

2 FIG. 9 .
FIG. 9. Nuclear spin squeezing under non-linear quadrupole interaction for different values of asymmetry parameter η.The quadrupole interaction can produce spin cat states (SCSs) for η = 0, however, for higher values of η, only spin squeezed states are formed.One can note that the period of non-linear evolution decreases as η increases, for e.g. the coherent state reappears faster for η = 1 (third figure from right), in comparison to η = 0.5 (second figure from right).
FIG. 10.Dynamics when η = 1 for different Zeeman interaction strength γB0 and quadrupole strength ωq = 2π ×40 ≈ 250 kHz.SCSs are formed by converting two-axis counter twisting to one-axis twisting by adding a large linear interaction induced by the Zeeman interaction (γB0 ≫ ωq).