## Long-lived metastable-qubit memory

Xiaoyang Shi, 1, ∗ Jasmine Sinanan-Singh, 1 Kyle DeBry, 1,2 Susanna L. Todaro, 1, † Isaac L. Chuang, 1 and John Chiaverini 1,2

1 Center for Ultracold Atoms, Research Laboratory of Electronics,

Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA

2 Lincoln Laboratory, Massachusetts Institute of Technology, Lexington, Massachusetts 02421, USA

Coherent storage of quantum information is crucial to many quantum technologies. Long coherence times have been demonstrated in trapped-ion qubits, typically using the hyperfine levels within the ground state of a single ion. However, recent research suggests qubits encoded in metastable states could provide architectural benefits for quantum information processing, such as the possibility of effective dual-species operation in a single-species system and erasure-error conversion for fault-tolerant quantum computing. Here we demonstrate long-lived encoding of a quantum state in the metastable states of a trapped ion. By sympathetically cooling with another ion of the same species and constantly monitoring for erasure errors, we demonstrate a coherence time of 136(42) seconds with a qubit encoded in the metastable 5 D 5 2 / state of a single 137 Ba + ion. In agreement with a model based on empirical results from dynamical-decoupling-based noise spectroscopy, we find that dephasing of the metastable levels is the dominant source of error once erasure errors are removed.

Qubits encoded in metastable states of atomic systems show great potential for quantum information processing [15]. Such states provide access to additional electronic manifolds with high isolation from the more generally used ground-state manifold-coherent and dissipative operations may be performed almost completely independently in these different subspaces, leading to novel utilization of trapped ions and neutral atoms for quantum logic. For instance, the flexibility that comes with access to separate sets of states suggests new architectures for quantum computing that can circumvent some of the challenges presented by dual-species methods [6, 7], including those stemming from mass differences between the constituent qubit species.

One possible limitation of metastable states for quantum information processing when compared to the often-employed ground states is the set of additional channels for scattering and decay that the former possess. It is therefore of paramount importance to investigate the behavior of metastable-states when storing quantum information, especially in the presence of other atomic systems that are being used for ancillary operations relevant to carrying out quantum algorithms and quantum error correction. Previous investigations have demonstrated few-to-tens-of-millisecond-scale [8, 9] to second-scale [10] coherence in metastable levels in neutral atoms and ions, respectively, but a validation of metastablestate coherence commensurate with the longer coherence times shown for ground states is a prerequisite for effective use of these multiple manifolds in high-fidelity quantum-state storage and manipulation.

the specific ion or atom that decayed will be identifiable, and could subsequently be returned to the qubit subspace [7, 11]. While decay of qubit states is in general undesirable, in light of the potential benefits of stimulated-Raman-based quantum logic [13], using metastable qubits in tandem with erasure conversion of Raman scattering errors can nonetheless be advantageous. Specifically, this capability strengthens the prospect of quantum fault tolerance, as erasure quantum error correcting codes can have very high thresholds, due to the fact that erasure errors are easier to detect than more general Pauli errors [14]. Demonstrating that such erasure-error detection is possible while maintaining coherence is therefore a key development toward realizing the promise of metastable-qubit processing.

Interestingly, the same aspects of metastable states that may lead to the above-mentioned limitations can also provide new enhancements not straightforwardly present in ground states. The finite state lifetimes of metastable levels lead to decay from the qubit space-however, since this decay is typically to the ground-state manifold, which can be simultaneously used for quantum state detection with isolation from the metastable levels, a decay event can be detected. Moreover, this leakage event from the metastable qubit space can be converted to an erasure error [4, 11, 12], an error in a specific location, since

Here we demonstrate the storage of a quantum state in the metastable D 5 2 / level of 137 Ba + while sympathetically cooling the memory ion using the ground-state manifold of another ion of the same species. We also simultaneously monitor the fluorescence from the memory ion to check for leakage errors. Since neither the sympathetic cooling light nor the light used for resonance-fluorescence-based state detection couples strongly to the memory-qubit manifold, this error detection technique can be effected in real-time to identify decay events. By proceeding with memory qubit analysis conditioned on the detection of no decay event, a coherence time of 136(42) seconds is observed, approximately four times the metastable state's lifetime, and comparable to ground-state memory times. The use of mid-algorithm measurement allows us to observe leakage errors explicitly to extrapolate the expected qubit lifetime that an erasure error correction scheme should provide, and hence validates the use of metastable states for high-fidelity quantum logic more generally.

While we use post-selection to evaluate only experimental trials in which the metastable qubit does not decay to the ground state, we point out that the leakage error detection method demonstrated here is generally useful, well beyond post-selected experiments. Crucially, since the process of monitoring for leakage errors does not affect the mem-

FIG. 1. Experiment setup and level structure. (a) Level diagram of 137 Ba + . The highlighted states in the D 5 / 2 and S 1 / 2 states are the memory qubit and cabinet shelving states respectively. The hyperfine structure of other levels is neglected for clarity. (b) Metastable qubit transition frequency vs. ambient magnetic field. The solid line is a fit to the data with free parameters being the curvature, an overall frequency shift, and a conversion factor between coil current and magnetic field strength. (c) A photograph of the microfabricated trap used in this work. (d) Surface electrode trap and laser configuration. The electrodes for applying RF currents to drive transitions between hyperfine states are marked in blue and the trapping RF electrodes are marked in purple. The DC electrodes are marked in grey. The laser propagation directions are labeled by arrows. Lasers propagating in the same direction are combined and labeled with parentheses. All beams address both ions equally.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00975v2-figure-e3cc8239-1f12-4a51-a8b4-1841f873e4e9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00975v2-figure-465268a8-2ecb-4993-9c7c-0d7cc4712a3f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00975v2-figure-a9393366-4699-4f84-985b-ba640ef54a98.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00975v2-figure-02b5946a-af52-457b-8710-6bc9a26ed01a.jpg)

ory qubit, the signal from the erasure detection can be used in real-time to verify that no error has occurred [15] and to continue a quantum algorithm; this is to be contrasted with post-selection in a system in which measurement outcomes can be acted upon only after the quantum information is destroyed. A key consequence is that the long coherence time measured here could be realistically utilized for quantum information processing if leveraged as part of an erasure-based error-correcting code.

perfine sublevels | F = 2 , mF = 2 ⟩ and | F = 3 , mF = 2 ⟩ are the | 0 ⟩ and | 1 ⟩ qubit states, respectively. Based on the measured properties of this species [17], we expect the qubit transition frequency to be approximately 62.7 MHz and first-order insensitive to magnetic field fluctuation at a magnetic field of 3.34 G. To verify the magnetic field sensitivity, we perform spectroscopy at different ambient fields and the result is shown in Fig. 1(b). The second-order field sensitivity is calculated to be 385 Hz/Gauss 2 .

As we utilize OMG [6]/dual-type [7] encoding in these experiments, individual addressing is required to change the qubit type of specific ions starting in the nominally same state after ion initialization. This is achieved here using a simple micromotion-assisted individual addressing technique and projection. Furthermore, metastable-qubit transition energies are subject to environmental perturbations, chief among these being magnetic field fluctuations. We mitigate dephasing from this source by encoding in a field-insensitive qubit, by surrounding the ion with a superconducting shield, and by applying dynamical decoupling to the memory qubit.

Experiment . We use 137 Ba + ions for both the memory and the ancilla ion, with the relevant level structure shown in Fig. 1(a). The memory qubit is encoded in the metastable D5 2 / level, which has a lifetime of TD = 30.14(30) s [16]. Hy-

As shown in Fig. 1(c,d), we use a microfabricated, surfaceelectrode, RF Paul trap in a cryogenic system (see more detail in [18]) for confining the ions. The ion-electrode distance is approximately 50 um and the trap RF drive is at 36 MHz. The trap is cooled to 6 K. For driving RF transitions within the D 5 2 / manifold, we can apply an RF current of a few milliamperes through the middle electrodes (marked in blue), with the other end of the electrodes shorted to ground. Only one of these electrodes is used in this manner for the results presented here. Two coils mounted on the intermediate (50 K) stage of the cryocooler in a near-Helmholtz configuration are used for generating a magnetic field perpendicular to the trap surface. There are four rings of niobium (superconducting critical temperature approximately 9.2 K) shielding placed co-axially to the trap-chip normal above and below

FIG. 2. Experiment sequence and dynamical decoupling scheme. (a) Eight-step operation sequence for sympathetically cooled quantum memory experiment. "MM": micromotion; "g": ground state qubit; "m": metastable qubit. (b) The CPMG sequence. π pulses with 90 ◦ phases relative to the initial π / 2 pulse are applied interspersed with a delay time τ which sets the filter function that suppresses (or samples) the noise spectrum.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00975v2-figure-e898a415-f909-452a-928c-dfda0d27f659.jpg)

the chip, with central circular holes of various sizes. These rings significantly suppress magnetic field variation along the vertical direction (chip-normal). More information about, and characterization of, the niobium shielding can be found in the supplementary material, Sec. A 1.

For loading ions, two photons at 554 nm and 403 nm [19] ionize atoms produced by an effusive oven. For Doppler cooling and state detection, the 493 nm laser is frequencymodulated using an electric-optical modulator (EOM). The laser frequency is set so that the carrier drives the S 1 / 2 | F = 2 ⟩ to P 1 / 2 | F = 2 ⟩ transition, and two microwave tones at approximately 1.47 GHz and 4.0 GHz drive the other three transitions between the hyperfine-split S 1 2 / and P1 2 / levels. The two tones are switchable using RF switches to allow for state preparation [20] in a particular sublevel. For repumping population from the D 3 2 / state, the 650-nm laser is also modulated by an EOM with two tones so that, along with the carrier, these tones drive the F = 2 ← → 1 1 , ← → 1 0 , ← → 1 3 , ← → 2 transitions between the D 3 2 / and P 1 2 / level.

in the S 1 / 2 | F = 1 , mF = 0 ⟩ state using a similar scheme as in [20], except we do not perform a polarization-based step with the 493 nm laser. To have a consistent assignment of cooling and memory ion, a micromotion-assisted individual addressing technique is used (see supplementary material, Sec. A 2) to excite the memory ion to the D 5 2 / state. State detection is conducted to verify that only one ion is shelved to the D 5 2 / state. The process is repeated until a one-ion-dark event is detected.

For the 614 nm laser, which is used to quench the D 5 2 / state, two acousto-optic modulators (AOMs) set up in a double-pass configuration are used for frequency modulation and switching. We find that it is sufficient to have a single frequency tone for effective quenching of all the sublevels in the D 5 2 / state. An optical shutter is added to the beam path for higher extinction of the 614 nm laser beam during memory-qubit coherence measurements, to avoid unwanted AC Stark shifts and quenching of the D 5 2 / state. A 1762 nm laser locked to a high finesse cavity (StableLaser Systems) and amplified by a thulium-doped fiber amplifier (Cybel, Perseus) is used for state preparation and pre-readout shelving of the memory qubit.

For experiments with two ions, the pulse sequence is shown in Fig. 2a. The ions are first Doppler-cooled using the 493 nm and the 650 nm laser beams. We then prepare both ions

To measure the memory-qubit coherence time, an RF current is sent to the through-wire electrode to drive the transition between the | 0 ⟩ and | 1 ⟩ states. The radio-frequency signal is generated via direct digital synthesis (DDS, M-labs, Urukul) synchronized to an ovenized crystal oscillator, and we drive pulses to implement the CPMG dynamical decoupling sequence [21], a modified Ramsey sequence (see Fig. 2b). The Doppler cooling beams are kept on during the Ramsey sequence for sympathetic cooling and erasure error detection. To detect leakage errors, we make a measurement of scattered photons using a photomultiplier tube every 4.05 ms to check if the memory ion decayed from the D 5 2 / state. The sequence is terminated when spontaneous decay is detected via fluorescence from both ions.

After the Ramsey sequence, the qubit | 0 ⟩ state is shelved back to the ground state for state detection. To avoid shelving the cooling ion to the D 5 2 / state during the same operation, the 493 nm laser is first switched on for 1 ms to shelve the cooling-ion population into the D 3 2 / state. Two 'cabinet shelving' pulses [20] are then used to drive the | 0 ⟩ population to the S 1 / 2 | F = 1 , mF = 0 1 , ⟩ states to achieve a higher shelving fidelity. The state of the memory qubit is then extracted by fluorescence detection (the cooling ion remaining dark).

Results. The result of the coherence time measurement with leakage detection T ( LD ) 2 is shown in Fig. 3. A coherence time of 136(42) s is extracted from an exponential fit.

FIG. 3. Experimental measurement of qubit coherence with no erasure detection (green squares), with real-time erasure detection (blue circles), and expected coherence time based on measured noise spectrum (purple dashed line, see Sec. A 6). A coherence time of 136(42) seconds is measured with erasure error detection, while 22(2) seconds is measured without erasure error detection, limited by decay from the qubit states (see text). Error bars correspond to one standard deviation. Solid lines are exponential decay fits to the data points.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00975v2-figure-a1e31d20-0aec-4f05-8bb3-566774356b6a.jpg)

We also perform the experiment without leakage detection, and the coherence time in that case, T 2, is measured to be 22(2) s. These two coherence times taken together are consistent with the lifetime of the D 5 / 2 state TD if we assume that the coherence time of the metastable qubit in the presence of decay ( T 2) is dependent on the combination of decay from the qubit states and pure dephasing, the latter characterized by a time T φ . In that case, we would assume the relationship 1 / T 2 ≥ 1 / T 1 + 1 / T φ (note that the appropriate decoherence rate due to state decay here is 1 / T 1 rather than 1 / ( 2 T 1 ) since both qubit states can decay). For T φ = T ( LD ) 2 = 136 42 ( ) s and T 1 = TD , we obtain T 2 ≤ 24 6 7 6 . ( . ) s for the decoherence time in the presence of leakage, consistent with the measured 22(2) s.

To understand the limits to the measured coherence time with leakage detection, we model the dephasing as follows: fluctuations of the qubit energy splitting are captured by including a dephasing term β ( ) t in addition to the bare qubit splitting ω 0 in the Hamiltonian H = ¯ h 2 ( ω 0 + ( ) β t ) σ z , where σ z is the usual Pauli matrix. We employ a set of CPMG pulse sequences, each consisting of a series of N total π pulses spaced by time τ to empirically estimate the frequency noise spectrum S ( ω ) associated with the noise process β ( ) t via qubitnoise spectroscopy[22].

The effect of applying a CPMG sequence is similar to applying a filter to the noise spectrum with a center frequency of ω c = 2 π × 1 2 τ and width 1 / T , for total delay time T = N τ , as described by a frequency domain filter function | ˜ y ( ω , T ) | . We perform experiments for a set of CPMG sequences with

FIG. 4. Qubit-noise spectroscopy and fitted qubit-splitting frequency-noise spectrum. The green, blue, and pink data points signify noise amplitudes based on coherence measurements of the metastable qubit using CPMG sequences of varying length and number of pulses N . The noise spectrum, modeled as described in the text, is a fit to these data. The effective filter function of the CPMG sequence used in the coherence-time measurement shown in Fig. 3 is depicted in the lower panel.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00975v2-figure-5c74d2e3-19c9-4748-87a8-154ce75009e6.jpg)

varying N and measure the coherence decay χ ( T ) :

$$\chi ( T ) = \int _ { 0 } ^ { \infty } S ( \omega ) | \tilde { y } ( \omega, T ) | ^ { 2 } d \omega \, \quad \quad ( 1 )$$

(see Sec. A 5 for a full derivation). We fit (1) with experimental coherence decay data points { χ ( N ) exp ( T ) } to extract S ( ω ) at measured total delays as shown in Fig. 4. We fit a power-law spectrum to the data (see Sec. A 6 for fitting details) and find that the noise spectrum has a minimum at ω ≈ 2 π × 5 Hz and increases rapidly at both lower and higher frequencies. Using this spectrum and the applied CPMG sequence parameters for the coherence time measurement ( τ = 0 1 s), we predict a co-. herence decay as shown by the dashed line in Fig. 3, consistent with the measured coherence.

Discussion. Environmental noise and imperfections in control fields both contribute to the observed noise spectrum. Potential high-frequency noise may arise from phase noise in the qubit drive, performed via direct digital synthesis of an RF current delivered to the excitation electrodes, though separately measured phase noise is too low to account for the full noise amplitude. Another possible source of high-frequency noise is magnetic-field fluctuations: the attenuation provided by the superconducting shield decreases with increasing frequency [23].

One major source of low-frequency noise is fluctuation in the amplitude of the applied RF voltage used for ion confinement. The asymmetry of the trap-electrode structure in the vertical ( z ) direction leads to uncancelled magnetic fields arising from the current in the trap-RF electrodes. Due to the small frequency difference between the trap RF frequency and

the qubit frequency, this magnetic field off-resonantly couples to the qubit and causes AC Zeeman shifts. Fluctuations of the RF amplitude thus lead to qubit-frequency fluctuations. Variation of the trapping field also induces frequency fluctuation due to the quadrupole moment of the D 5 2 / state [24], but as shown in the supplementary material (see Sec. A 3), this effect does not contribute significantly in our case. Noise caused by the trap-RF voltage could be suppressed by the addition of active stabilization of its amplitude [25]. Additional low-frequency noise may be caused by magnetic-field fluctuations in spatial directions with insufficient shielding; more complete shielding could potentially reduce such noise at the cost of reduced optical access.

Decoherence also arises from spontaneous Raman photon scattering from the Doppler cooling light [26], which is also applied during the memory storage period, but this effect is negligible due to the large detunings from relevant transitions (see Sec. A 4).

The coherence time observed here for clock states in metastable levels of a trapped ion is comparable in scale to those achievable in the ground-state manifolds of atoms. This suggests that such a state encoding could be useful for quantum computing in OMG-style architectures [6], for instance, using multiple 137 Ba + ions as a register.

As shown here, select ions could be addressed to encode metastable-state qubits, while other ions left in the ground-state manifold could be used for sympathetic cooling, allowing simultaneous periodic erasure-error detection and qubit readout. Such a mid-circuit measurement capability, as demonstrated here, means quantum error correction could also be applied, providing significant functionality in a singlespecies system. This demonstration of several key OMGarchitecture elements thus provides concrete steps toward the realization of trapped-ion quantum computing at scale.

## ACKNOWLEDGEMENTS

I. L. C. acknowledges support by the NSF Center for Ultracold Atoms. This research was supported by the U.S. Army Research Office through grant W911NF-20-1-0037. Support is also acknowledged from the U.S. Department of Energy, Office of Science, National Quantum Information Science Research Centers, Quantum Systems Accelerator, under Air Force Contract No. FA8702-15-D-0001. Any opinions, findings, conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the Department of Energy.

∗ shix@mit.edu

† Present address: Oxford Ionics, Oxford, U.K.

[1] K. Toyoda, H. Shiibara, S. Haze, R. Yamazaki, and S. Urabe, Experimental study of the coherence of a terahertz-separated

- metastable-state qubit in 40 Ca + , Phys. Rev. A 79 , 023419 (2009).
- [2] J. A. Sherman, M. J. Curtis, D. J. Szwer, D. T. C. Allcock, G. Imreh, D. M. Lucas, and A. M. Steane, Experimental recovery of a qubit from partial collapse, Phys. Rev. Lett. 111 , 180501 (2013).
- [3] P. Schindler, D. Nigg, T. Monz, J. T. Barreiro, E. Martinez, S. X. Wang, S. Quint, M. F. Brandl, V. Nebendahl, C. F. Roos, M. Chwalla, M. Hennrich, and R. Blatt, A quantum information processor with trapped ions, New Journal of Physics 15 , 123012 (2013).
- [4] S. Ma, G. Liu, P. Peng, B. Zhang, S. Jandura, J. Claes, A. P. Burgers, G. Pupillo, S. Puri, and J. D. Thompson, High-fidelity gates and mid-circuit erasure conversion in an atomic qubit, Nature 622 , 279-284 (2023).
- [5] K. DeBry, J. Sinanan-Singh, C. D. Bruzewicz, D. Reens, M. E. Kim, M. P. Roychowdhury, R. McConnell, I. L. Chuang, and J. Chiaverini, Experimental quantum channel discrimination using metastable states of a trapped ion, Phys. Rev. Lett. 131 , 170602 (2023).
- [6] D. T. C. Allcock, W. C. Campbell, J. Chiaverini, I. L. Chuang, E. R. Hudson, I. D. Moore, A. Ransford, C. Roman, J. M. Sage, and D. J. Wineland, omg blueprint for trapped ion quantum computing with metastable states, Applied Physics Letters 119 , 214002 (2021), https://doi.org/10.1063/5.0069544.
- [7] H.-X. Yang, J.-Y. Ma, Y.-K. Wu, Y. Wang, M.-M. Cao, W.X. Guo, Y.-Y. Huang, L. Feng, Z.-C. Zhou, and L.-M. Duan, Realizing coherently convertible dual-type qubits with the same ion species, Nature Physics 18 , 1058 (2022).
- [8] S. Pucher, V. Klüsener, F. Spriestersbach, J. Geiger, A. Schindewolf, I. Bloch, and S. Blatt, Fine-structure qubit encoded in metastable strontium trapped in an optical lattice, Phys. Rev. Lett. 132 , 150605 (2024).
- [9] G. Unnikrishnan, P. Ilzhöfer, A. Scholz, C. Hölzl, A. Götzelmann, R. K. Gupta, J. Zhao, J. Krauter, S. Weber, N. Makki, H. P. Büchler, T. Pfau, and F. Meinert, Coherent control of the fine-structure qubit in a single alkaline-earth atom, Phys. Rev. Lett. 132 , 150606 (2024).
- [10] L. Feng, Y.-Y. Huang, Y.-K. Wu, W.-X. Guo, J.-Y. Ma, H.-X. Yang, L. Zhang, Y. Wang, C.-X. Huang, C. Zhang, L. Yao, B.-X. Qi, Y.-F. Pu, Z.-C. Zhou, and L.-M. Duan, Realization of a crosstalk-avoided quantum network node using dual-type qubits of the same ion species, Nature Communications 15 , https://doi.org/10.1038/s41467-023-44220-z (2024).
- [11] M. Kang, W. C. Campbell, and K. R. Brown, Quantum error correction with metastable states of trapped ions using erasure conversion, PRX Quantum 4 , 020358 (2023).
- [12] Y. Wu, S. Kolkowitz, S. Puri, and J. D. Thompson, Erasure conversion for fault-tolerant quantum computing in alkaline earth rydberg atom arrays, Nature Communications 13 , 4657 (2022).
- [13] I. D. Moore, W. C. Campbell, E. R. Hudson, M. J. Boguslawski, D. J. Wineland, and D. T. C. Allcock, Photon scattering errors during stimulated raman transitions in trapped-ion qubits, Physical Review A 107 , 10.1103/physreva.107.032413 (2023).
- [14] M. Grassl, T. Beth, and T. Pellizzari, Codes for the quantum erasure channel, Phys. Rev. A 56 , 33 (1997).
- [15] S. R. Vizvary, Z. J. Wall, M. J. Boguslawski, M. Bareian, A. Derevianko, W. C. Campbell, and E. R. Hudson, Eliminating qubit-type cross-talk in the omg protocol, Phys. Rev. Lett. 132 , 263201 (2024).
- [16] Z. Zhang, K. J. Arnold, S. R. Chanu, R. Kaewuam, M. S. Safronova, and M. D. Barrett, Branching fractions for P 3 2 / decays in ba + , Phys. Rev. A 101 , 062515 (2020).

- [17] N. C. Lewty, B. L. Chuah, R. Cazan, M. D. Barrett, and B. K. Sahoo, Experimental determination of the nuclear magnetic octupole moment of 137 Ba + ion, Phys. Rev. A 88 , 012518 (2013).
- [18] J. M. Sage, A. J. Kerman, and J. Chiaverini, Loading of a surface-electrode ion trap from a remote, precooled source, Phys. Rev. A 86 , 013417 (2012).
- [19] N. Greenberg, B. M. White, P. J. Low, and C. Senko, Trapping Ba + with seven-fold enhanced efficiency utilizing an autoionizing resonance (2023), arXiv:2307.07627 [quant-ph].
- [20] F. A. An, A. Ransford, A. Schaffer, L. R. Sletten, J. Gaebler, J. Hostetter, and G. Vittorini, High fidelity state preparation and measurement of ion hyperfine qubits with i ≥ 1 2 , Phys. Rev. Lett. 129 , 130501 (2022).
- [21] S. Meiboom and D. Gill, Modified spin-echo method for measuring nuclear relaxation times, Review of Scientific Instruments 29 , 688 (1958).
- [22] Y. Wang, M. Um, J. Zhang, S. An, M. Lyu, J.-N. Zhang, L.M. Duan, D. Yum, and K. Kim, Single-qubit quantum memory exceeding ten-minute coherence time, Nature Photonics 11 , 646 (2017).
- [23] A. Kamitani, T. Yokono, and T. Yokono, Frequency dependence of magnetic shielding performance of HTS plates in mixed states , Tech. Rep. NIFS-PROC-46 (Inst. of Information Sciences and Electronics, Japan, 2000).
- [24] W. M. Itano, External-field shifts of the 199 Hg + optical frequency standard, J. Res. Natl. Inst. Stand. Technol. 105 , 829 (2000).
- [25] K. G. Johnson, J. D. Wong-Campos, A. Restelli, K. A. Landsman, B. Neyenhuis, J. Mizrahi, and C. Monroe, Active stabilization of ion trap radiofrequency potentials, Review of Scientific Instruments 87 , 053110 (2016), https://pubs.aip.org/aip/rsi/articlepdf/doi/10.1063/1.4948734/15778092/053110\_1\_online.pdf.
- [26] R. Ozeri, C. Langer, J. Jost, B. DeMarco, A. Ben-Kish, B. Blakestad, J. Britton, J. Chiaverini, W. Itano, D. Hume, D. Leibfried, T. Rosenband, P. Schmidt, and D. Wineland, Hyperfine coherence in the presence of spontaneous photon scattering, Physical Review Letters 95 , 10.1103/physrevlett.95.030403 (2005).
- [27] C. J. Myatt, B. E. King, D. Kielpinski, D. Leibfried, Q. A. Turchette, C. S. Wood, W. M. Itano, C. R. Monroe, and D. J. Wineland, Trapped ions, entanglement, and quantum computing, in Methods for Ultrasensitive Detection , Vol. 3270, edited by B. L. Fearey, International Society for Optics and Photonics (SPIE, 1998) pp. 131 - 137.
- [28] N. K. Lysne, J. F. Niedermeyer, A. C. Wilson, D. H. Slichter, and D. Leibfried, Individual addressing and state readout of trapped ions utilizing rf micromotion (2024).
- [29] M. J. Biercuk, A. C. Doherty, and H. Uys, Dynamical decoupling sequence construction as a filter-design problem, Journal of Physics B: Atomic, Molecular and Optical Physics 44 , 154002 (2011).
- [30] R. Kubo, Generalized Cumulant Expansion Method, Journal of the Physical Society of Japan 17 , 1100 (1962).
- [31] G. Fubini, Opere scelte , Vol. 2 (Cremonese, 1958) pp. 243-249.
- [32] G. A. Álvarez and D. Suter, Measuring the spectrum of colored noise by dynamical decoupling, Phys. Rev. Lett. 107 , 230501 (2011).

## Appendix A: Supplementary material

## 1. Superconducting shields for reduction of magnetic-field fluctuations

The niobium shield configuration is shown in Fig. 5. Four pieces of niobium sheet with various shapes are mounted around the ions to suppress magnetic field variation. To measure the amount of suppression, we choose two transitions ( | F = 1 , mF = 0 ⟩ to | F = 1 , mF = ± ⟩ 1 ) between the S 1 2 / and D5 2 / states for spectroscopy. These two transitions are similarly sensitive to magnetic field changes, but with opposite sign. The difference in the two transition frequencies thus provides a measure of the magnetic field at the ion position. We then perform spectroscopy at two different ambient fields by changing the current applied to the magentic-field coils and measuring the field change ( ∆ B) as measured at the ion location. The difference between the two applied magnetic fields is approximately 1.4 G. This measurement is repeated at different temperatures of the trap stage, and the result is shown in Fig. 6. A suppression factor of ∼ 240 is measured for the temperature used in the experiments presented in the main paper (6 K).

## 2. Micromotion-assisted individual addressing

Micromotion-assisted individual addressing has been proposed [27] and demonstrated for a two-dimensional ion array [28]. In this section, we describe a realization in a twoion chain. To address the memory ion, the two-ion chain is rotated about the cooling ion, perpendicularly to the trap surface, via controlled variation of the trap-electrode voltages (see Fig. 7a). When the memory ion is off-axis radially, it experiences a larger micromotion amplitude compared to the cooling ion. We add a 2 ms delay after the trap voltages are updated to account for low-pass filters present on the electrodes. We then drive the micromotion sideband of the S1 2 / | F = 1 , mF = 0 ⟩ to D 5 / 2 | F = 2 , mF = 2 ⟩ transition with a global beam. The π -time of the micromotion sideband for the memory ion is about 50 us (6 us for the carrier). Rabi flopping on the micromotion sideband using a global beam is shown in Fig. 7b. The fitted amplitude of oscillation is 99.9%. Combined with state detection after the micromotion-based state preparation, the error for memory-ion preparation is negligible. After the memory ion is prepared in the | 0 ⟩ state, the ion chain is rotated back to the original position.

## 3. Trap-induced qubit frequency fluctuations

The AC Zeeman shift due to the oscillating current in the RF electrodes scales quadratically with magnetic field, and thus linearly with the applied RF power. On the other hand, the quadrupole shift is caused by the gradient of the trapping field, and so scales with the square root of RF power. In a separate experiment, we measure the qubit transition frequency as a function of the trap-RF drive power. The result is shown in Fig. 8. The linear relation between the qubit frequency and RF power suggests that the AC Zeeman shift is likely the dominant source of trap-induced frequency shift. Thus, slow fluctuations of the trap RF amplitude could contribute to the measured low-frequency noise.

## 4. Spontaneous scattering rate

Spontaneous Raman scattering is a source of errors affecting memory-qubit coherence, but since the predominant scattering channel is D 5 2 / → P3 2 / → S1 2 , and since there is also / decay to the D 3 2 / level, many scattering errors are correctable via our leakage detection protocol [11]. We estimate the scattering rate from the Doppler cooling beams on the D 5 2 / qubit using the Kramers-Heisenberg formula [13, 26]

$$\Gamma _ { i \rightarrow f } \approx \sum _ { i, f } g ^ { 2 } \gamma \left | \frac { a _ { i \rightarrow f } ^ { ( 3 / 2 ) } } { \Delta } \right | ^ { 2 }. \quad \quad ( A 1 )$$

where g = E 0 µ 2¯ h , γ is the decay rate of the P3 2 / manifold, µ is the relevant dipole matrix element ⟨ P3 2 / | ⃗ d · σ l | D5 2 / ⟩ , a i ( 3 / 2 ) → f = ∑ ∑ q k ∈ P 3 / 2 ⟨ f | d · σ q | k ⟩⟨ k | d · σ l | ⟩ i / µ 2 is the effective transition amplitude for laser polarization σ l and scattered photon polarization σ q , and ∆ is the detuning from the D 5 2 / -P3 2 / resonance. The measured optical power in the 493 nm (650 nm) beam is approximately 150 µ W (20 µ W), with a 20 µ m(30 µ m) beam diameter.

From the reduced matrix element, we have that µ ≈ 5 a.u. [16] and ∑ i , f | a i ( 3 / 2 ) → f | 2 ∼ 1, which leads to a scattering rate of ∼ × 2 10 -3 Hz, similar for both of these cooling beams. For the undetectable scattering errors, those proceeding via D5 2 / → P3 2 / → D5 2 , / however, the scattering rate can be bounded at a lower rate of ∼ × 5 10 -4 Hz for each of these beams due to the branching ratio of the P 3 2 / → D5 2 / decay (0.23). We therefore expect this effect to be a very small contribution to the total memory error in this demonstration.

## 5. Dynamical-decoupling based noise filtering

Magnetic field noise arises from power fluctuations in the trapping RF and induces AC Zeeman shifts of the memory levels. Such fluctuations of the energy splitting are captured by including a dephasing term to the bare qubit splitting:

$$H = \frac { \hbar } { 2 } ( \omega _ { 0 } + \beta ( t ) ) \sigma _ { z }. \text{ \quad \ \ } ( A 2 )$$

Frequency and phase noise in the driving RF pulses also lead to effective dephasing during the Ramsey sequence. This

FIG. 5. The niobium shielding configuration. (a) Schematics of the relative positions of the niobium shields. Black solid lines represent the niobium plates and the dashed line represents the position of a hole. The blue rectangle marks the trap chip and the cyan circles represent the ions (not to scale). (b) In this CAD rendering of the niobium shielding, the trap is the black chip. The 6 K thermal shield and trap mount are not shown.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00975v2-figure-ae72a073-6efd-41b0-94a0-1523eaa31bd8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00975v2-figure-cb4cdb7a-e4a0-4651-a46e-9f9a6e4115cc.jpg)

FIG. 6. Ratio of measured magnetic field change to applied magnetic field change at different temperatures. The externally applied magnetic field is varied by ∼ 1 4 Gauss and the local change in the . magnetic field is measured via ion spectroscopy as a function of trapstage temperature. The dashed blue line marks the critical temperature of niobium (approximately 9.2 K). The ambient field change is suppressed by a factor of ∼ 240 at 6 K. Error bars are smaller than the markers.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00975v2-figure-3409097e-c8d9-4eed-8b93-d586aec82386.jpg)

noise occurs due to local oscillator noise affecting the timing and frequency of driving RF pulses. Furthermore, the near field effects of the trapping RF field include quadrupole terms which will off resonantly couple E2 transitions leading to higher order Stark shifts. The effect of such energy level shifts is again dephasing noise, though noise contributions from quadrupolar terms are a few orders of magnitude smaller than dipolar terms.

Let us assume the noise β ( ) t is a stationary process ⟨ β ( t ′ + t ) β ( t ′ ) ⟩ = ⟨ β ( ) t β ( 0 ) ⟩ . Then, we may define the noise spectral density using the Wiener-Khinchin theorem as

$$\langle \beta ( t ^ { \prime } + t ) \beta ( t ^ { \prime } ) \rangle = \frac { 1 } { 2 \pi } \int _ { - \infty } ^ { \infty } S ( \omega ) e ^ { i \omega t } d \omega. \quad ( A 3 )$$

We will now connect the noise spectral density to the qubit coherence measured in the ion experiments to allow for noise spectrum determination, and subsequently for tailoring of the dynamical decoupling sequence for optimal noise rejection.

The advantage of CPMG pulses is that we can filter which components of noise spectra the ion experiences [29]. The CPMG sequence is comprised of π -pulses spaced by time τ . If we assume the π -pulses are instantaneous, a good approximation for the long free-evolution times explored here, the time evolution to first-order perturbation approximation is:

$$U _ { \text{CMPG} } & = \prod _ { j = 0 } ^ { N } e ^ { i \sigma _ { z } \int _ { t _ { j } } ^ { \prime j + 1 } \beta ( t ^ { \prime } ) \, d t ^ { \prime } } e ^ { i \pi \sigma _ { y } / 2 } \\ & = e ^ { i \int _ { - \infty } ^ { \infty } y ( t, T ) \beta ( t ) d t \sigma _ { z } },$$

where t j is the time of the j th π pulse with t 0 = 0 , t 1 = τ / 2 and t N + 1 = N τ = T . Here y t ( ) alternates between + 1 and -1 in successive free evolution periods, and the σ i 's are the usual Pauli matrices.

This periodic flipping of the Bloch vector is expressed using a modulation function ˜ y ( ω ) in the frequency domain:

$$\tilde { y } ( \omega, T ) & = \frac { 1 } { \sqrt { 2 \pi } } \int _ { 0 } ^ { \infty } y ( t, T ) e ^ { i \omega t } d t \\ & = \frac { 1 } { \omega \sqrt { 2 \pi } } \sum _ { j = 0 } ^ { N } ( - 1 ) ^ { j } \left ( e ^ { i \omega t _ { j + 1 } } - e ^ { i \omega t _ { j } } \right ).$$

The center of the filter function is at ω c = 2 π × 1 2 τ , and the width is 1 / T .

We calculate the coherence time of the memory qubit by performing a Ramsey protocol with or without a dynamical decoupling sequence during the Ramsey delay time:

$$\langle \sigma _ { y } \rangle = \frac { 1 } { 2 } \left ( 1 + \frac { \langle e ^ { i \theta } + e ^ { - i \theta } \rangle } { 2 } \right ) \text{ \quad \ \ } ( A 6 )$$

where θ ( T ) = ∫ ∞ -∞ y t ( , T ) β ( ) t dt is the accumulated phase during total delay time T . The Ramsey contrast can be re-written

FIG. 7. Micromotion-assisted individual addressing. (a) Schematics of the procedure. The cooling and memory ions are marked with "g" and "m" respectively. To address the memory ion, the two-ion chain is rotated around the cooling ion so that the Rabi frequency on the micromotion sideband is much higher for the memory ion. The dashed circle marks the position of the memory ion when micromotion is compensated. (b) Rabi flopping on the micromotion sideband with the rotated chain. The green and blue dots represent two-ion and one-ion excitation events respectively. The solid line represents fitting Rabi oscillation assuming a thermal state of motion.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00975v2-figure-0314a6b2-a0bd-47d9-b040-a77b6406fe7c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00975v2-figure-fa4beb91-5ca9-416a-a2b0-1596421a60bd.jpg)

quency domain:

FIG. 8. The qubit frequency responds linearly to the output amplitude of the function generator that supplies the trapping field. The solid line is a linear fit to the data points. This shows that the AC Zeeman shift due to the oscillating trapping field is the dominant effect of the trap RF on the qubit. Slow fluctuations of the drive power could thus help explain the observed low-frequency noise. Error bars are smaller than the markers

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00975v2-figure-77ab1973-698f-4315-84c7-7c8ac0bf8dc9.jpg)

using the cumulant expansion:

$$\frac { \langle e ^ { i \theta } + e ^ { - i \theta } \rangle } { 2 } = e ^ { - \chi ( T ) }, \quad \quad ( A 7 ) \quad \text{Us} \quad \\ \quad \text{we}$$

where the coherence χ ( T ) = ∑ k ( -i ) k k ! Ck is the cumulant expansion of the noise with Ck = ⟨ θ k ( T ) ⟩ [30].

We assume β ( ) t is Gaussian (i.e. fully characterized by its auto-correlation), and so only the second cumulant survives χ ( T ) = ⟨ θ 2 ( T ) ⟩ / 2. Using Fubini's theorem [31], we further simplify the integral for θ to obtain the coherence in the fre-

$$\chi ( T ) & = \frac { 1 } { 2 } \int _ { - \infty } ^ { \infty } d t ^ { \prime } \int _ { - \infty } ^ { \infty } d t ^ { \prime \prime } \langle \beta ( t ^ { \prime } ) \beta ( t ^ { \prime \prime } ) \rangle y ( t ^ { \prime \prime } ) y ( t ^ { \prime \prime } ) \\ & = \frac { 1 } { 2 } \int _ { - \infty } ^ { \infty } d t ^ { \prime } \int _ { - \infty } ^ { \infty } d t ^ { \prime \prime } \int _ { - \infty } ^ { \infty } \frac { d \omega } { 2 \pi } S ( \omega ) e ^ { i \omega ( t ^ { \prime } - t ^ { \prime \prime } ) } y ( t ^ { \prime } ) y ( t ^ { \prime \prime } ) \\ & = \int _ { 0 } ^ { \infty } S ( \omega ) | \tilde { y } ( \omega ) | ^ { 2 } d \omega.$$

This gives a relation between coherence and noise spectral density via the filter function, allowing for translation of qubit measurements under a particular modulation scheme to a noise spectrum.

The typical Ramsey protocol with modulation y t ( ′ ) = 1 for 0 &lt; t ′ &lt; T filters the noise as a sinc function with center about 0 and width decreasing with delay time as ∼ 1 / T .

$$\chi ( T ) = \frac { 1 } { 4 \pi } \int _ { - \infty } ^ { \infty } \frac { \sin ^ { 2 } ( \omega T / 2 ) } { ( \omega / 2 ) ^ { 2 } } S ( \omega ) d \omega. \quad ( A 9 )$$

Thus, the coherence for the typical Ramsey protocol is χ ( T ) = S ( 0 ) T / 2 and is primarily due to the noise spectral component at ω = 0. We may however extract noise strength for other frequencies by using other modulating functions.

Considering only the main peak (first harmonic) of ˜ y ( ω ) for more complex modulation sequences, we may approximate the noise spectrum from the measured Ramsey contrast. Using the definition of the δ function as δ ( x ) = 1 2 π ∫ ∞ -∞ e ikx dk , we have that ∫ ∞ -∞ d ω | ˜ y ( ω ) | 2 = ∫ ∞ -∞ dt y t | ( ) | 2 = T . Since the center of the filter function is at ω c = 2 π × 1 2 τ , we may approximate the noise at ω c from measured Ramsey contrast as

$$\chi ( T ) = S ( \omega _ { c } ) T / 2 ; \text{ \quad \ \ } ( A 1 0 )$$

there are higher harmonics to the filter function, but at lower amplitude, and we ignore them in this approximation. We

therefore can straightforwardly identify noise spectrum components to build up a spectrum by measuring ion coherence for experiments with different center frequencies and delay times.

## 6. Maximum likelihood noise model

The CMPG sequence is used to measure the noise spectrum affecting the memory-qubit energy splitting [22]. For these experiments, we keep the number of π -pulses fixed, varying the inter-pulse time τ , and measure the coherence time. We repeat this procedure with N = 2, 20, and 200 π -pulses. Using (A10), we can convert measured coherence times into noise spectral components [32]. To get an approximate noise spectrum from the noise spectroscopy measurements (see Fig. 4), we use a cost function to fit all the experimental data via maximum likelihood estimation. We assume the noise exhibits power-law behavior and we include higher harmonics of the filter function up to 1 MHz, parameterizing the spectrum as

$$S ( \omega, \{ \lambda _ { i } \} ) = \lambda _ { 0 } + \lambda _ { 1 } \omega + \lambda _ { 2 } / \omega + \lambda _ { 3 } / \omega ^ { 2 }, \quad ( \text{A11} )$$

where the λ i 's represent the relative strength of the various power-law frequency-scaling terms in the fit function. We generate coherences χ num from the fitted spectrum via

$$\chi _ { \text{num} } ( N, T, \{ \lambda _ { i } \} ) = \int _ { 0 } ^ { \infty } S ( \omega, \{ \lambda _ { i } \} ) | \tilde { y } ( \omega, N, T ) | ^ { 2 } d \omega. \ \ ( A 1 2 )$$

The cost function is given by weighted least squares error as

$$\mathcal { C } ( \{ \lambda _ { i } \} ) = \sum _ { \{ N, T \} } \frac { \left ( \chi _ { \exp } ( N, T ) - \chi _ { \text{num} } ( N, T, \{ \lambda _ { i } \} ) \right ) ^ { 2 } } { \sigma ^ { 2 } ( N, T ) } ; \ ( A 1 3 )$$

this compares the experimental coherence data χ exp ( N T , ) and respective variances σ 2 ( N T , ) from the N = 2 20 200 mea-, , surements with varying delay time T = N τ to a given power law description χ num. Optimizing this cost function leads to parameters λ 0 = 0 0085 . , λ 1 = 0 0138 . , λ 2 = 0 , λ 3 = 1 6631. . For these parameters, noise terms proportional to 1 / ω 2 and ω dominate the spectrum. The transition point between these two noise scaling regimes is approximately 5 Hz, and so we use a CPMG sequence centered about this minimum with τ = 0 1 s to most effectively suppress the effect of the noise . measured in the system. From the modeled noise spectrum, we calculate the expected qubit coherence time with erasure detection using A10 to be 115 s, in reasonable agreement with the measured coherence time (cf. Fig. 3).
Xiaoyang Shi Jasmine Sinanan-Singh Kyle DeBry Susanna L. Todaro Isaac L. Chuang John Chiaverini
Long-lived metastable-qubit memory
Coherent storage of quantum information is crucial to many quantum
technologies. Long coherence times have been demonstrated in trapped-ion
qubits, typically using the hyperfine levels within the ground state of a
single ion. However, recent research suggests qubits encoded in metastable
states could provide architectural benefits for quantum information processing,
such as the possibility of effective dual-species operation in a single-species
system and erasure-error conversion for fault-tolerant quantum computing. Here
we demonstrate long-lived encoding of a quantum state in the metastable states
of a trapped ion. By sympathetically cooling with another ion of the same
species and constantly monitoring for erasure errors, we demonstrate a
coherence time of 136(42) seconds with a qubit encoded in the metastable
$5D_{5/2}$ state of a single $^{137}$Ba$^+$ ion. In agreement with a model
based on empirical results from dynamical-decoupling-based noise spectroscopy,
we find that dephasing of the metastable levels is the dominant source of error
once erasure errors are removed.