Integrated 9 Be + multi-qubit gate device for the ion-trap quantum computer

H. Hahn, 2 G. Zarantonello, 2 M. Schulte, A. Bautista-Salvador, 2, 4 K. Hammerer, and C. Ospelkaus 2, 4 Physikalisch-Technische Bundesanstalt, Bundesallee 100, 38116 Braunschweig, Germany Institute of Quantum Optics, Leibniz University Hannover, Welfengarten 1, 30167 Hannover, Germany Institute for Theoretical Physics and Institute for Gravitational Physics (Albert-Einstein-Institute), Leibniz University Hannover, Appelstrasse 2, 30167 Hannover, Germany Laboratory for Nanoand Quantum Engineering, Leibniz University Hannover, Schneiderberg 39, 30167 Hannover, Germany


INTRODUCTION
Following the proposal by Cirac et al., 1 trapped atomic ions have shown to be a promising and pioneering platform for implementing elements of quantum information processing (QIP). 2,3 Qubits are encoded in the internal states of individual ions, and shared motional modes are used as a "quantum bus" for multi-qubit operations. Toward a large-scale universal quantum processor based on trapped-ion qubits, the "Quantum Charge-Coupled Device" (QCCD) 4,5 is considered as a possible scalable hardware implementation. It relies on microfabricated multi-zone ion-trap arrays, in which quantum information is processed in dedicated zones interconnected via ion transport. While some key requirements such as high-fidelity ion transport 6 and fault-tolerant singlequbit gates 7,8 have already been demonstrated in multiple setups, high-fidelity multi-qubit gates 9,10 below the fault-tolerant threshold still remain challenging. In this context, entangling gates driven by microwave fields 11,12 represent a technically less demanding alternative to laser-induced gates, as microwave signals can typically be controlled more easily than optical fields from highly specialized laser systems. The microwave approach avoids spontaneous scattering as a fundamental source of infidelities 13 and experimental fidelities [14][15][16][17] are approaching the fidelities of the best laser-driven gates. 9,10 Here, we focus on the near-field microwave 12 gate approach, where the leading sources of infidelity in implementations so far comprise the spatiotemporal stability of the microwave near-field pattern 14,18 or fluctuating AC Zeeman shifts. 17,19 We note that in the latter work, the error contribution arising from fluctuating AC Zeeman shifts has been reduced to ≲0.1% through the use of a dynamical decoupling scheme. 17 In this letter, we realize a two-qubit gate using a tailored microwave conductor embedded in a surface-electrode trap optimized to produce high oscillating magnetic near-field gradients and low residual fields at the ion position, thus directly addressing the main sources of error in previous near-field gates. The gate is realized on a field-independent hyperfine qubit in 9 Be + ions, a promising ion species for scalable QIP, [20][21][22] following the Mølmer-Sørensen (MS) [23][24][25][26] protocol. The implementation is based on an optimized single-conductor design, which can be thought of as the prototype of a scalable multi-qubit gate module for an ion-trap quantum computer based on surface-electrode trap arrays. The measured gate fidelity of 98.2 ± 1.2% is purely limited by technical imperfections, in agreement with a numerical analysis.

RESULTS AND DISCUSSION
The surface-electrode trap was fabricated at the PTB cleanroom facility employing the single-layer method as detailed in ref. 27 on an AlN substrate (for the present trap, we chose the single-layer process in order to quickly test improvements that were made to the setup and trap orientation compared with ref. 28 ). Gold electrodes are about 10 μm thick and separated by 5 μm gaps. Aiming to remove potential organic residuals on top of the electrode surfaces, the trap was cleaned in an ex situ dry-etching process before being installed in a UHV vacuum chamber at room temperature. Electrical connectivity is provided by wire bonding to a printed-circuit board for DC signal filtering and signal routing. Figure 1 shows a top view of the trap center. Besides the microwave conductor MWM, which produces a magnetic field quadrupole 29 suitable for driving multi-qubit gates, the electrode layout includes two microwave conductors labeled MWC to generate an oscillating magnetic field for global spin-state manipulation. The ions are radially and axially confined at an ion-to-electrode distance of about 70 μm using one radio frequency (RF) and 10 DC control electrodes, respectively. With an applied RF voltage of 150V pp at Ω RF ≃ 2π × 88.2 MHz, the motional-mode frequencies of the radial rocking modes of a two-ion 9 Be + crystal were measured to be (ω r1 , ω r2 ) ≃ 2π × (6.275, 6.318) MHz.
The static magnetic field B 0 defining the quantization axis at an angle of 45°with respect to the trap axis is produced by a hybrid setup consisting of two permanent magnet assemblies and a pair of compensation coils. 30 At the ion position, this setup generates a magnetic field of |B 0 | = 22.3 mT forming a first-order magnetic field-insensitive qubit 31 on the hyperfine levels 2 S 1/2 |F = 1, m F = 1〉 ≡ |↑〉 and 2 S 1/2 |F = 2, m F = 1〉 ≡ |↓〉 with an unperturbed transition frequency of ω 0 ≃ 2π × 1082.55 MHz, cf. Fig. 2. Here, F and m F represent the quantum numbers for the ion's total angular momentum and its projection on the quantization axis, respectively.
Ions are loaded using laser ablation combined with photoionization. 29 State initialization is done via optical pumping to the | 2, 2〉 state (also referred to as the bright state) and subsequent Doppler cooling on the closed-cycle transition 2 S 1/2 |2, 2〉 ↔ 2 P 3/2 m J ¼ 3 2 ; m I ¼ 3 2 (where m J and m I are the projections of the total electronic and nuclear angular momenta onto the quantization axis). Resolved sideband cooling is performed by a pair of counterpropagating Raman beams aligned along theẑ direction. Each sideband cooling cycle consists of a global π rotation on the hyperfine transition labeled "I" in Fig. 2, followed by an optical red sideband pulse on the qubit transition and a repumping sequence to transfer all population back to the initial bright state. The repumping sequence comprises multiple microwave induced π rotations on the qubit transition and laser pulses on the 2 S 1/2 |1, transition similar to ref. 18 Employing interleaved sideband cooling sequences on both radial rocking modes of a two-ion crystal, we measure a mean phonon number n of the modes n r1 ' 0:27 and n r2 ' 0:11 assuming a thermal distribution. The heating rate was determined to be _ n r2 ' 28s À1 . In an earlier measurement with similar radial-mode frequencies, the heating rates of a single-ion's low-frequency (LF) and high-frequency (HF) radial mode were found to be about _ n LF ' 116s À1 and _ n HF ' 122s À1 , respectively.
To perform a MS two-qubit entangling gate on the qubit transition using near-field microwaves, 12 we subsequently initialize the qubits in |↑↑〉 and apply a bichromatic microwave current to MWM at the frequencies ω RSB = ω 0 + Δ − (ω r2 + δ) and ω BSB = ω 0 + Δ + (ω r2 + δ). Here, δ is the gate detuning from the highfrequency rocking mode of N = 2 ions at ω r2 , and Δ is the differential AC Zeeman shift of the unperturbed qubit transition induced by the bichromatic field. In the ideal case, the implemented dynamics can be described by the Hamiltonian where Ω is the gate Rabi frequency, a r2 (a y r2 ) represents the mode's annihilation (creation) operator and we define σ ± = 1/2(σ x ± iσ y ) with σ x and σ y being the Pauli matrices. Following, 32,33 we apply the interaction on |↑↑〉 in order to produce the maximally entangled state Ψ where K is an integer number) and calculate the resulting state preparation fidelity F Ψ h jρ Ψ j i ¼ 1=2ðP "" þ P ## Þ þ ρ "";## by determining the far off-diagonal element ρ ↑↑,↓↓ of the system's density matrix ρ as well as the population probabilities in |↓↓〉, |↑↓〉, and |↓↑〉, and |↑↑〉 given by P ↓↓ , P ↑↓,↓↑ , and P ↑↑ , respectively.
Experimentally, we infer these quantities by adding a π/2 analysis pulse after the gate operation and observe the global ion fluorescence as a function of the analysis pulse's phase ϕ a employing fluorescence detection on the closed-cycle transition for 400 μs. Prior to detection, we transfer the population in |↑〉 back to the bright state and shelve the population in |↓〉 to |1, −1〉 (also referred to as the dark state) using a sequence of π rotations on the transitions "I" and "II"-"IV", respectively. The populations are determined by repeating the experiment 200 times for each value of ϕ a , and fitting the resulting histograms to a sum of weighted Poisson distributions representing 0, 1, and 2 ions bright. We calibrate the mean of the distributions to a reference two-ion Ramsey experiment which we assume to produce an ideal outcome following the same procedure of 14 (see Supplementary Methods section). To account for off-resonant optical pumping effects, we modify the three Poissonians to include all depumping processes from the shelved dark state to the bright state during the detection interval. 34 By consequence, we expect the resulting populations to compensate imperfect state preparation and detection within the present level of accuracy. Finally, ρ ↑↑,↓↓ can be deduced by calculating the parity Π(ϕ a ) = P ↓↓ (ϕ a ) + P ↑↑ (ϕ a ) − P ↑↓+↓↑ (ϕ a ), while ϕ a is varied and extracting the magnitude |A Π | of the parity oscillation equal to |2ρ ↑↑,↓↓ |. 33 Figure 3 shows the population and parity signal after performing the two-qubit gate operation using the MWM conductor. From sinusoidal fits to the extracted populations (solid lines), we calculate a corresponding gate fidelity of 98.2 ± 1.2%. The error on the fidelity is derived from In more detail, we apply a power of~5.5 W to each sideband tone, resulting in gate Rabi frequency of Ω/2π = 1.071 kHz and an inferred near-field gradient of~19 Tm −1 . For a single sideband, the corresponding residual magnetic field at the ion position is highly suppressed by the optimized conductor geometry, resulting in an on-resonance Rabi frequency of about Ω c /2π~450 kHz. In case the bichromatic gate drive is applied, the accompanied differential AC Zeeman shift is measured to be Δ/2π = 4.37 kHz and is predominantly caused by σ-components of the residual magnetic fields as the induced shifts of the π-components mainly cancel each other assuming an equal power in both sideband tones. In order to suppress off-resonant spin excitations, we adiabatically shape the envelope of the gate pulse at its beginning and end with a 2 -μs long cumulative error function using a microwave envelope stabilization circuit based on a fast analog multiplier and a digital PI controller 35 with the regulator setpoint generated by an arbitrary waveform generator. 36 We optimize the gate by fixing the pulse duration to the expected value based on the gate Rabi frequency and subsequently scanning the gate detuning resulting in the highest fidelity. Following this procedure, we find an optimal gate time and detuning of τ = 808 μs and δ/2π = 3.4 kHz, respectively, corresponding to K = 3 loops in motional-phase space. Here, the mismatch to the theoretically predicted detuning at δ theory /2π = 3.71 kHz was tracked back to a systematic frequency offset from an independent radial-mode frequency measurement as well as a radial mode frequency "chirp" of 0.3 Hz μs −1 during the gate pulse. Qualitatively, a similar effect is also observed in other experiments using near-field gradients and appears to be inherent to warm-up processes in the microwave-generating structures. 19 This hypothesis is supported by our observation of a saturation behavior of the "chirp" at~1 ms. Consequently, we reduce the impact of the "chirp" by preceding the gate with a 400 -μs long warm-up pulse with the duration chosen conservatively to avoid excessive heating of trap structures.
In order to identify current infidelity contributions in producing the maximally entangled state, we simulate the dynamics of the system using a master equation considering experimentally determined input parameters. This becomes necessary, as the exact propagator may no longer be obtained analytically in the presence of additional error sources. The master equation is given by where the Hamiltonian is H ¼H MS þ H m þ H z þ H spec , and ðΩ B ðtÞσ þ j a y r2 e Àiδt þ Ω R ðtÞσ þ j a r2 e iδt Þ þ H:c: is an extension to the ideal case presented in Eq. (1). Here, we have assumed equal Rabi frequencies and phases for both ions, which is true in the experiment to the best of our knowledge. Further, H m ¼ δ ϵ ðtÞa y r2 a r2 describes the instability of the rocking mode frequency, H z ¼ Δ ϵ ðtÞ=2 P j σ z j gives the uncompensated AC Zeeman shift resulting from shot-to-shot microwave power fluctuations, general Ω B and Ω R allow an imbalance in the two sideband Rabi frequencies, and time dependencies of the pulse shape are taken into account by the time-dependent parameters. Couplings via additional, off-resonant motional modes are included by the term H spec ¼ Ω r1 =2 P N j¼1 ðσ þ j þ σ À j Þða r1 e iðΔνþδÞt þ a y r1 e ÀiðΔνþδÞt Þ whereby we limit ourselves to the nearest mode only (with Ω r1 ≃ Ω), which contributes the largest error of this kind. In addition to the unitary dynamics, motional heating to a thermal state with n th ) 1 and qubit decoherence are considered by the Lindblad terms 32 L h ρ ¼ γ h ðD½a r2 ρ þ D½a y r2 ρÞ with the heating rate γ h in phonons per second and L d ρ ¼ γ d =2 P j D½σ z j ρ with the decoherence rate γ d , respectively, where D½Ôρ ¼ÔρÔ y ÀÔ yÔ ρ=2 À ρÔ yÔ =2. Table 1 lists contributions of the different error sources to the infidelity 1 À F . These values result from numerical simulations of the quantum dynamics according to Eq. (2) considering the ideal gate dynamics with addition of the corresponding noise in the form we described above. All simulations were done with QuTiP 37 and used a truncated Hilbert space for the motional mode. For our analysis, including the first 25 Fock states was sufficient to reach convergence given the low initial thermal distribution and the small motional displacements during the gate.
We examined the following effects, which we considered to be the most relevant, in more detail: the largest error according to our investigation results from the frequency instability of the rocking mode, which establishes the gate dynamics. This effect consists of two parts. On the one hand, normally distributed variations of the frequency with a standard deviation of ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi hðδ ϵ =δÞ 2 i q ¼ 1:1 10 À2 , inferred from a measured instantaneous linewidth of 2π × 101 Hz in a calibration scan directly before the gate measurement. On the other hand, a frequency "chirp" within each gate that we model by a linear increase of 0.3 Hz μs −1 within the first 600 μs, and subsequent constant frequency leading to in total 1.3% infidelity. While the mode fluctuations can be reduced by actively stabilizing the amplitude and frequency of the trap RF signal, 38,39 the "chirp" can be reduced by e.g., longer warm-up pulses. Simulations (see Fig. 4) also us to identify the individual contribution of each effect in view of further improvements.
The second largest contribution is caused by driving the identical spin-spin interaction via the additional low-frequency Fig. 3 Population and parity oscillation as a function of the phase ϕ a of a π/2 analysis pulse applied after preparing the maximally entangled state Þutilizing the near-field pattern of the MWM conductor. The solid lines are sinusoidal fits to the observed data, while each data point represents the average of 200 experiments. The error bars result from the fit of the weighted Poissonians to the experimental data. a Population in |↓↓〉 (P ↓↓ ; red circles), |↑↓〉 and |↓↑〉 (P ↑↓,↓↑ ; green squares) and |↑↑〉 (P ↑↑ ; blue diamonds). The fits yield P ↑↑ + P ↓↓ = 0.990 ± 0.021. b Parity Π(ϕ a ) = P ↓↓ (ϕ a ) + P ↑↑ (ϕ a ) − P ↑↓,↓↑ (ϕ a ) oscillation with a fitted amplitude of A Π = 0.975 ± 0.012 rocking mode. The measured frequency spacing was Δν = ω r2 −ω r1 = 2π × 42.5 kHz and ground-state cooling to n r1 ¼ 0:27 was applied, resulting in an infidelity of 5.2 × 10 −3 from the competing gate dynamics. As this effect scales inversely with the radial-mode splitting, it can be suppressed by engineering the trap potentials or suitable pulse sequences. 40,41 Heating of the motion and decoherence of the qubits contribute~0.1% and <0.1% infidelity, respectively. Again, both effects do not represent a fundamental limit to the gate performance and can be improved experimentally. 31,42 Off-resonant scattering on carrier transitions can lead to undesired excitations inside and outside the qubitmanifold and thus contribute a gate error. Here, an excitation other than on the qubit transition is much less probable due to the higher-frequency difference of the driving field, which is >150 MHz detuned from the next spectator transition. Direct simulation of this effect was not performed due to the vastly different timescales of the gate dynamics (~kHz) and the carrier processes (~GHz) which would have considerably increased the runtime of the numerical simulations. We instead performed direct measurements on a single qubit 19,43 to evaluate the extent of this error, which is then quantified to be <2.3 × 10 −3 . Infidelities below 6.3 × 10 −4 resulted from distortion of the pulse shape, whereby we combine here the influence of adiabatic switching on and off as well as small changes of the Rabi frequency and AC Zeeman shift during the pulses which result from power transients on the ideally rectangular signal. Stabilization of the microwave power allowed to reduce the shot-to-shot fluctuations of the power, and accordingly of the AC Zeeman shift, to an extent that the simulated infidelity of 1.1 × 10 −4 contributes only insignificantly. The same applies to the imbalance of Rabi frequencies, cf. Table 1.
In conclusion, we have demonstrated a microwave-driven twoqubit gate with 9 Be + ions using a single microwave conductor with an optimized design embedded in a surface-electrode ion trap. The design of the MWM conductor has been developed to generate a high magnetic near-field gradient with low residual field at the ion position, thus suppressing AC Zeeman shift fluctuations, an inherent error source of the near-field approach, to a simulated infidelity contribution of~10 −4 . In contrast, according to the presented error budget, the main contributions can all be decreased upon technical improvements; the by far biggest of these (1.3%) is consistent with the measured two-qubit gate infidelity of 1.8 ± 1.2%. In addition to technical modifications to the apparatus, more elaborate gate schemes employing Walsh modulation 44 or continuous dynamic decoupling 17 can be applied in order to increase the gate fidelity as required for fault-tolerant quantum computation. In the future, the MWM conductor design can be used as an entangling gate unit of a "QCCD" architecture purely employing microwave-driven quantum gates. Moreover, the conductor design can also be integrated into a scalable multilayer trap. 27,28

DATA AVAILABILITY
The datasets generated during and/or analyzed during the current study are available from C. Ospelkaus (christian.ospelkaus@iqo.uni-hannover.de) on reasonable request. Fig. 4 Simulated infidelity assuming different variations of the relative motional mode frequency and lengths of the frequency "chirp" while considering otherwise ideal gate dynamics. In all cases, the gate parameters are as specified in the text. The experimental condition of the presented two-qubit gate is given by the intersection of the vertical line with the solid blue line. The inset illustrates the observed linear frequency increase of the selected motional mode of about 0.3 Hz μs −1 at the beginning of the gate pulse until time T chirp The infidelity values result from numerical simulations of the quantum dynamics according to Eq. (2), including the respective noise effect with a strength given by the measured parameter specified in the second column H. Hahn et al.