## Abstract

Quantum adiabatic evolutions find a broad range of applications in quantum physics and quantum technologies. The traditional form of the quantum adiabatic theorem limits the speed of adiabatic evolution by the minimum energy gaps of the system Hamiltonian. Here, we experimentally show using a nitrogen-vacancy center in diamond that, even in the presence of vanishing energy gaps, quantum adiabatic evolution is possible. This verifies a recently derived necessary and sufficient quantum adiabatic theorem and offers paths to overcome the conventionally assumed constraints on adiabatic methods. By fast modulation of dynamic phases, we demonstrate near–unit-fidelity quantum adiabatic processes in finite times. These results challenge traditional views and provide deeper understanding on quantum adiabatic processes, as well as promising strategies for the control of quantum systems.

## INTRODUCTION

Coherent control on quantum systems is a fundamental element of quantum technologies that could revolutionize the fields of information processing, simulation, and sensing. A powerful and universal method to achieve this control is the quantum adiabatic technique, which exhibits intrinsic robustness against control errors ensured by the quantum adiabatic evolution (*1*). Besides important applications in quantum state engineering (*2*, *3*), quantum simulation (*4*–*6*), and quantum computation (*7*–*11*), the quantum adiabatic evolution itself also provides interesting properties such as Abelian (*12*) or non-Abelian geometric phases (*13*), which can be used for the realization of quantum gates. However, the conventional quantum adiabatic theorem (*14*, *15*), which dates back to the idea of extremely slow and reversible change in classical mechanics (*14*, *16*), imposes a speed limit on the quantum adiabatic methods, that is, for a quantum process to remain adiabatic, the changes in the system Hamiltonian at all times must be much smaller than the energy gap of the Hamiltonian. On the other hand, to avoid perturbations from the environment, high rates of change are desirable. This tension can impose severe limitations on the practical use of adiabatic methods. Despite the long history and broad applicability, it was found recently that key aspects of quantum adiabatic evolution remain not fully understood (*17*, *18*) and that the condition in the conventional adiabatic theorem is not necessary for quantum adiabatic evolution (*19*, *20*).

In this work, we experimentally demonstrate adiabatic evolutions with vanishing energy gaps and energy level crossings, which are allowed under a recently proven quantum adiabatic condition (*20*) that is based on dynamical phases instead of energy gaps, by using a nitrogen-vacancy (NV) center (*21*) in diamond. In addition, we reveal that using discrete jumps along the evolution path allows quantum adiabatic processes at unlimited rates, which challenges the view that adiabatic processes must be slow. By jumping along the path, one can even avoid path points where the eigenstates of the Hamiltonian are not feasible in experiments. Furthermore, we theoretically and experimentally demonstrate the elimination of all the nonadiabatic effects on the system evolution of a finite evolution time by driving the system along the geodesic that connects initial and final states, as well as combating system decoherence by incorporating pulse sequences into adiabatic driving.

## RESULTS

### Experimental study of the necessary and sufficient quantum adiabatic condition

To describe the theory for experiments, we consider a quantum system driven by a Hamiltonian *H*(λ) for adiabatic evolution. In terms of its instantaneous orthonormal eigenstates ∣ψ* _{n}*(λ)〉 (

*n*= 1,2, …) and eigenenergies

*E*(λ), the Hamiltonian is written as

_{n}*H*(λ) = ∑

_{n}*E*(λ)∣ψ

_{n}*(λ)〉〈ψ*

_{n}*(λ)∣. For a given continuous finite evolution path, ∣ψ*

_{n}*(λ)〉 changes gradually with the configuration parameter λ. In our experiments, λ corresponds to an angle in some unit and is tuned in time such that λ = λ(*

_{n}*t*) ∈ [0,1]. The system dynamics driven by the Hamiltonian is fully determined by the corresponding evolution propagator

*U*(λ). It is shown that one can decompose the propagator

*U*(λ) =

*U*

_{adia}(λ)

*U*

_{dia}(λ) as the product of a quantum adiabatic evolution propagator

*U*

_{adia}(λ) that describes the ideal quantum evolution in the adiabatic limit and a diabatic propagator

*U*

_{dia}(λ) that includes all the diabatic errors (

*20*). In the adiabatic limit,

*U*

_{dia}(λ) =

*I*becomes an identity matrix and the adiabatic evolution

*U*=

*U*

_{adia}(λ) fully describes the geometric phases (

*12*,

*13*) and dynamic phases accompanying the adiabatic evolution (that is, the deviation from adiabaticity

*U*−

*U*

_{adia}vanishes). This decomposition guarantees that both

*U*

_{adia}(λ) and

*U*

_{dia}(λ) are gauge invariant, i.e., invariant with respect to any chosen state basis.

According to the result of (*20*), the error part satisfies the first-order differential equation (ℏ = 1)*U*_{dia}(0) = *I*. The generator *W*(λ) describes all the nonadiabatic transitions. On the basis of ∣ψ* _{n}*(0)〉, the diagonal matrix elements of

*W*(λ) vanish, i.e., 〈ψ

*(0)∣*

_{n}*W*(λ)∣ψ

*(0)〉 = 0. The off-diagonal matrix elements*

_{n}_{n,m}(λ) ≡ ϕ

*(λ) − ϕ*

_{n}*(λ) is the difference of the accumulated dynamic phases ϕ*

_{m}*(λ) on ∣ψ*

_{n}*(λ)〉, and the geometric part*

_{n}*G*

_{n,m}(λ) =

*e*

^{i[γm(λ) − γn(λ)]}

*g*

_{n,m}(λ) consists of the geometric functions

Equation 2 shows that the differences of dynamic phases ϕ_{n,m} are more fundamental than the energy gaps in suppressing the nonadiabatic effects because the energies *E _{n}* do not explicitly appear in these equations. According to (

*20*), when the dynamic phase factors at different path points add destructively

*n*≠

*m*and any λ ∈ [0,1] of a finite path with bounded

*G*

_{n,m}(λ) and

*G*

_{n,m}(λ) and

*G*

_{tot}= ∑

_{n≠m}max∣

*G*

_{n,m}(λ′)∣ and

*20*). In the limit ϵ → 0, the system evolution is adiabatic along the entire finite path with

*U*

_{dia}(λ) →

*I*. For a zero gap throughout the evolution path, the evolution is not adiabatic because ϵ

_{n,m}(λ) = λ is not negligible because of the constructive interference of the dynamic phase factors at different path points. For a large constant gap, the destructive interference gives a negligible ϵ

_{n,m}and hence an adiabatic evolution.

To experimentally verify the adiabatic condition in Eq. 3 by an NV center, we construct the Hamiltonian for adiabatic evolution in the standard way (*2*, *3*), that is, we apply a microwave (MW) field to drive the NV electron spin states ∣*m*_{s} = 0〉 ≡ ∣ −*z*〉 and ∣*m*_{s} = +1〉 ≡ ∣*z*〉 (see Fig. 1A and Materials and Methods for experimental details). The Hamiltonian *H*(λ) under an on-resonant MW field reads_{1}(λ)〉 = ∣+_{λ}〉 and ∣ψ_{2}(λ)〉 = ∣−_{λ}〉. Here_{g}λ. We define the initial eigenstates ∣±*x*〉 ≡ ∣±_{0}〉 and the superposition states

In the traditional approach where the Hamiltonian varies slowly with a nonvanishing gap, the strength of relative dynamic phase ϕ_{1,2} = ϕ_{1}(λ) − ϕ_{2}(λ) rapidly increases with the change of the path parameter λ, giving the fast oscillating factor *e*^{iϕ1,2} with a zero mean [see Fig. 1C for the case of a constant gap Ω(λ) = Ω_{0}]. Therefore, the right-hand side of Eq. 1 is negligible in solving the differential equation, leading to the solution *U*_{dia}(λ) ≈ *I*. As a consequence of the adiabatic evolution *U* ≈ *U*_{adia}(λ), the state initialized in an initial eigenstate of the Hamiltonian follows the evolution of the instantaneous eigenstate (Fig. 1D).

However, a quantum evolution with a nonvanishing gap and a long evolution time is not necessarily adiabatic. In Fig. 1 (E and F), we show a counterexample that increasing the energy gap in Fig. 1C to Ω(λ) = Ω_{0}[2 + cos (Ω_{0}λ*T*)] ≥ Ω_{0} will not realize adiabatic evolution because, in this case, the ϵ_{1,2}(λ) in Eq. 3 and the *G*_{1,2}(λ) are not negligible. For example, ϵ_{1,2}(λ) = *J*_{2}(1)λ ≈ 0.115λ (*J _{n}* being the Bessel function of the first kind) whenever the difference of dynamic phases is a multiple of 2π. This counterexample is different from the previously proposed counterexamples (

*17*,

*18*,

*20*,

*22*), where the Hamiltonian contains resonant terms that increase

*15*), which is based on the energy gap and diverges at

*E*−

_{n}*E*= 0, does not guarantee quantum adiabatic evolution. On the contrary, the condition in Eq. 3 based on dynamic phases, i.e., integrated energy differences, does not diverge for any energy gaps. We note that fast amplitude fluctuations on the control fields (hence energy gaps) can exist in adiabatic methods [e.g., see (

_{m}*23*)] because of their strong robustness against control errors. By adding errors in the energy gap shown in Fig. 1E, the adiabaticity of the evolution is substantially enhanced (Fig. 2), showing that the situations to have nonadiabatic evolution with a fluctuating energy gap are relatively rare.

We demonstrate that adiabatic evolution can be achieved even when the energy spectrum exhibits vanishing gaps and crossings as long as Eq. 3 is satisfied for a sufficiently small ϵ. As an example, we consider the energy gap of the form *a*∣ > 1 [see Fig. 1G for the case of *a* ≈ 2.34, where *e*^{iϕ1,2} parameterized by the parameter λ = *t*/*T* is fast oscillating (Fig. 1G) with a zero mean and realizes quantum adiabatic evolution for a sufficiently large total time *T* (see Fig. 1H). In fig. S1, we show how the adiabaticity can also be preserved when gradually introducing energy level crossings.

### Unit-fidelity quantum adiabatic evolution within a finite time

Without the restriction to nonzero energy gaps, it is possible to completely eliminate nonadiabatic effects and to drive an arbitrary initial state ∣Ψ_{i}〉 to a target state ∣Ψ_{t}〉 of a general quantum system by the quantum adiabatic evolution of a finite time duration. We demonstrate this by driving the system along the geodesic for maximal speed [see, e.g., (*24*, *25*) for more discussion on the geodesic in quantum mechanics]. The system eigenstate _{1}(0)〉 and ∣ψ_{1}(1)〉 along the geodesic by varying λ = 0 to λ = 1, with its orthonormal eigenstate *24*, *25*). An example of the geodesic path for a single qubit is given in Eq. 5, which intuitively can be illustrated by the shortest path on the Bloch sphere (see Fig. 1B). We find that, along the geodesic, the nonzero elements *g*_{2,1}(λ) and *g*_{1,2}(λ) are constant. We adopt the sequence theoretically proposed in (*20*) that changes the dynamic phases at *N* equally spaced path points λ = λ* _{j}* (

*j*= 1,2, … ,

*N*). By staying at each of the points

*U*

_{dia}(1) =

*I*because

*W*(λ) commutes and

In other words, by jumping on discrete points λ* _{j}*, the system evolution at λ = 1 is exactly the perfect adiabatic evolution

*U*

_{adia}although the evolution time is finite, and an initial state ∣Ψ

_{i}〉 will end up with the adiabatic target state ∣Ψ

_{t}〉 =

*U*

_{adia}∣Ψ

_{i}〉. To realize the jumping protocol, we apply rectangular π pulses at the points λ

*without time delay between the pulses because, between the points λ*

_{j}*, the Hamiltonian has a zero energy gap and its driving can be neglected (see Fig. 3C). The simulation results in Fig. 3 show how the transition from the standard continuous protocol to the jumping one gradually increases the fidelity of adiabatic evolution.*

_{j}We experimentally compare the jumping protocol with the continuous one along the geodesic given in Eq. 5 by measuring the fidelity between the evolved state and the target state ∣Ψ_{t}〉 that follows the ideal adiabatic evolution *U*_{adia}. The continuous protocol has a constant gap and a constant sweeping rate as in Fig. 1C. As shown in Fig. 4 (A and B), for the case of a geodesic half circle (θ_{g} = π), the jumping protocol reaches unit fidelity within the measurement accuracy, while the standard continuous driving has much lower fidelity at short evolution times. The advantage of the jumping protocol is more prominent when we traverse the half-circle path back and forth [see Fig. 4 (C to F) for the results of a total path length of 6θ_{g}]. We observe in Fig. 4 that the constant gap protocol provides unit state transfer fidelity only when the initial state is an eigenstate of the initial Hamiltonian ∣Ψ_{i}〉 = ∣*x*〉 and when the relative dynamic phase accumulated in a single half circle is *k* = 1,2, …) (see Materials and Methods). However, the phase shifts on the system eigenstates accompanying adiabatic evolution cannot be observed when the initial state is prepared in one of the initial eigenstates. Therefore, in Fig. 4, we also compare the fidelity for the initial state ∣Ψ_{i}〉 = ∣*y*〉, which is a superposition of the initial eigenstates ∣± *x*〉. The results confirm that the jumping protocol achieves exactly the adiabatic evolution *U*_{adia} within the experimental uncertainties.

### Robustness of quantum adiabatic evolution via jumping

To demonstrate the intrinsic robustness guaranteed by adiabatic evolutions, in Fig. 5, we consider large random driving amplitude errors in the jumping protocol. We add random Gaussian distributed errors with an SD of 50% to the control amplitude. To simulate white noise, we change the amplitude after every 10 ns in an uncorrelated manner. Despite the large amplitude errors, which can even cause energy level crossings, during the evolution (see Fig. 5A for a random time trace), a change of fidelity is hardly observable in Fig. 5B. Additional simulations in fig. S2 also demonstrate the robustness to amplitude fluctuations with different kinds of noise correlation, i.e., Gaussian white noise, Ornstein-Uhlenbeck process modeled noise, and static random noise. The robustness of the jumping protocol can be further enhanced by using a larger number *N* of points along the path (fig. S3).

While it is different from dynamical decoupling (DD) (*20*, *26*), the jumping protocol can suppress the effect of environmental noise through a mechanism similar to DD. Therefore, the fidelity is still high even when the evolution time is much longer than the coherence time,

### Avoiding unwanted path points in adiabatic evolution

Without going through all the path points, the jumping protocol has advantages to avoid path points (i.e., Hamiltonian with certain eigenstates) that cannot be realized in experiments. As a proof-of-principle experiment, we consider the Landau-Zener (LZ) Hamiltonian (*27*)* _{α}* (α=

*x*,

*y*,

*z*) being the Pauli matrices. Because Δ is nonzero in the LZ Hamiltonian, tuning the system eigenstates to the eigenstates ∣±

*z*〉 of σ

*requires*

_{z}*B*→ ± ∞. Therefore, for a perfect state transfer from ∣Ψ

_{z}_{i}〉 = ∣−

*z*〉 to ∣Ψ

_{t}〉 = ∣

*z*〉, by using the standard continuous protocol, it is required to adiabatically tune

*B*from −∞ to +∞ (see insets of Fig. 6A). The experimental implementation of

_{z}*B*= ± ∞ however requires an infinitely large control field, which is a severe limitation. In our experiment, a large

_{z}*B*field can be simulated by going to the rotating frame of the MW control field with a large frequency detuning. The experimental realization of

_{z}*B*→ ± ∞ can be challenging in other quantum platforms. For example, for superconducting qubits where Δ/(2π) could be as large as 0.1 GHz, the tuning range of

_{z}*B*/(2π) is usually limited to a couple of gigahertz or even of the same order of magnitude as Δ/(2π) (

_{z}*28*). For a two-level quantum system comprising Bose-Einstein condensates in optical lattices, the maximum ratio of

*B*/Δ is determined by the band structure (

_{z}*29*). For singlet-triplet qubits in semiconductor quantum dots, the exchange interaction for the control of

*B*is positively confined (

_{z}*30*). On the contrary, with the jumping approach, one can avoid the unphysical points such as

*B*= ± ∞ as an infinitely slow and continuous process is not required and can achieve high-fidelity state transfer as shown in Fig. 6.

_{z}As a remark, we find that our jumping protocol with *N* = 1 (i.e., a Rabi pulse) specializes to the optimized composite pulse protocol (*29*) but has the advantage of requiring no additional strong π/2 pulses at the beginning and the end of the evolution. Moreover, by applying the jumping protocol with *N* = 1 to the adiabatic passage proposed in (*31*), we obtain the protocol that has been used to experimentally generate Fock states of a trapped atom (*32*). When, instead of a single target point, high-fidelity adiabatic evolution along the path is also desired, we can use the jumping protocol with a larger *N*.

## DISCUSSION

In summary, our experiments demonstrated that energy level crossings and vanishing gaps allow and can even accelerate quantum adiabatic evolutions, challenging the traditional views that adiabatic control must be slow and that unit-fidelity adiabatic processes require an infinite amount of evolution time. By experimentally verifying a recently derived quantum adiabatic condition, we have shown that the quantum dynamic phases are more fundamental than energy gaps in quantum adiabatic processes. Owing to rapid changes of these phases, nonadiabatic transitions can be efficiently suppressed and fast-varying Hamiltonians can still realize quantum adiabatic evolutions. Our results break the limit imposed by the conventional adiabatic methods that originate from the traditional concept of extremely slow change in classical mechanics (*14*, *16*), allowing fast quantum adiabatic protocols with unit fidelity within finite evolution times. In addition, the freedom of using vanishing gaps provides the ability to avoid unphysical points in an adiabatic path and allows incorporating pulse techniques (*26*) into a quantum adiabatic evolution to suppress environmental noise for long-term robust adiabatic control. While it is possible to mimic the infinitely slow quantum adiabatic evolution by using additional counterdiabatic control, i.e., shortcuts to adiabaticity (*29*, *33*–*37*), the implementation of the counterdiabatic control can be exceedingly intricate because it may need interactions absent in the system Hamiltonian (*36*, *37*). Furthermore, the counterdiabatic control unavoidably changes the eigenstates of the initial Hamiltonian and introduces additional control errors (*36*, *37*). However, because our protocol uses the intrinsic adiabatic path that follows the eigenstates of the Hamiltonian, no additional control is required. As a consequence, our methods avoid the use of difficult or unavailable control resources and share the intrinsic robustness of adiabatic methods. With the removal of the prerequisites in the conventional adiabatic conditions, namely nonzero gaps and slow control, our results provide new directions and promising strategies for fast, robust control on quantum systems.

## MATERIALS AND METHODS

### Adiabatic evolution along the geodesics of a general quantum system

For two arbitrary states (e.g., entangled states and product states) of a general quantum system, ∣Ψ_{i}〉 and ∣Ψ_{t}〉, one can write _{i,t} and θ_{g} being real. Here, θ_{g} is the path length connecting ∣Ψ_{i}〉 and ∣Ψ_{t}〉 by the geodesic, and we set ϕ_{i, t} = 0 by a proper gauge transformation (*24*). The geodesic (*24*, *25*) that connects ∣Ψ_{i}〉 and ∣Ψ_{t}〉 by varying λ = 0 to λ = 1 can be written as ∣ψ_{1}(λ)〉 = *c*_{i}(λ)∣Ψ_{i}〉 + *c*_{t}(λ)∣Ψ_{t}〉, where the coefficients _{1}(λ)〉 in terms of the system eigenstates, we chose an orthonormal state ∣ψ_{2}(0)〉 ∝ (*I* −∣Ψ_{i}〉〈Ψ_{i}∣) ∣Ψ_{t}〉 if _{t}〉 is equivalent to ∣Ψ_{i}〉 up to a phase factor [i.e., _{2}(0)〉 can be an arbitrary orthonormal state. Then, the geodesic and its orthonormal state can be written as *g*_{n,m} = 0 for other combinations of *n* and *m*. Along the geodesic, if one changes the dynamic phases with a π phase shift only at each of the *N* equally spaced path points λ* _{j}* = (2

*N*)

^{−1}(2

*j*− 1), with

*j*= 1,2, … ,

*N*, then the operators

*W*(λ) at different λ commute and we have

*U*=

*U*

_{adia}does not have any nonadiabatic effects.

### Hamiltonian of the NV center under MW control

Under a magnetic field *b _{z}* along the NV symmetry axis, the Hamiltonian of the NV center electron spin without MW control reads

*S*is the electron spin operator,

_{z}*D*≈ 2π × 2.87 GHz is the ground-state zero field splitting,

*b*is the magnetic field, and γ

_{z}_{e}= − 2π × 2.8 MHz G

^{−1}is the electron spin gyromagnetic ratio (

*21*). Following the standard methods to achieve a controllable Hamiltonian for quantum adiabatic evolution (

*2*,

*3*), we applied an MW field

*m*

_{s}= 0 and

*m*

_{s}= 1 levels to form a qubit with the qubit states ∣

*z*〉 ≡ ∣

*m*

_{s}= 1〉 and ∣ −

*z*〉 ≡ ∣

*m*

_{s}= 0〉. The MW frequency ω

_{MW}may also be tuned by the parameter λ to realize a controllable frequency detuning δ(λ) with respect to the transition frequency of

*m*

_{s}= 0 and

*m*

_{s}= 1 levels. In the standard rotating frame of the MW control field, we had the general qubit Hamiltonian under the MW control (

*21*)

*∣±*

_{z}*z*〉 = ±∣±

*z*〉 and [cos (θ

_{g}λ)σ

*+ sin (θ*

_{x}_{g}λ)σ

*] ∣±*

_{y}_{λ}〉 = ±∣±

_{λ}〉, where the states ∣±

_{λ}〉 are given in Eq. 5.

By setting the MW detuning to δ(λ) = 0, we achieved the Hamiltonian in Eq. 4, which in terms of the Pauli operators reads

By varying the parameter λ, the system eigenstates follow the geodesics along the equator of the Bloch sphere where the north and south poles are defined by the states ∣ ± *z*〉. Here, the energy gap Ω(λ) is directly controlled by the amplitude of the MW field. On the other hand, by using a constant Rabi frequency Ω(λ) = Δ and a tunable frequency detuning δ(λ) = *B _{z}*(λ), we obtained the LZ Hamiltonian

*H*

_{LZ}(λ) given in Eq. 7.

### Adiabatic evolution by continuous driving with a constant gap

Consider a conventional adiabatic driving where a constant amplitude driving field rotates around the *z* axis, with the Hamiltonian *t*/*T* along a circle of latitude with σ_{θ} = σ* _{z}*cosθ + σ

*sinθ in a total time*

_{x}*T*. The difference of the accumulated dynamic phases at λ = 1 on the two eigenstates is ϕ = Ω

*T*. One can show that the system evolution at λ = 1 reads

The ideal adiabatic evolution is obtained by using Eq. 9 in the adiabatic limit *T* → ∞ (i.e., ϕ → ∞)

Without the part of dynamic phases, *U*_{adia} describes geometric evolution, and for a cyclic evolution (i.e., θ_{g} = 2π), the geometric evolution is described by the Berry’s phases ±π(cosθ − 1). By comparing *U*_{adia} and *U* or by using the results of (*20*), the nonadiabatic correction is given by

In the adiabatic limit *T* → ∞ (i.e., ϕ → ∞), *U*_{dia} = *I* is the identify operator. We note that, when the phase factor of the state is irrelevant, one can perform perfect state transfer by this driving if the initial state is prepared in an initial eigenstate of the driving Hamiltonian *H*(λ) (i.e., an eigenstates of σ_{θ}). From Eq. 10, *U*_{dia} is diagonal on the basis of σ_{θ} when *U*′ ∝ *I*. As a consequence, when *U*′ ∝ *I* and ∣Ψ_{i}〉 is prepared as an eigenstate of σ_{θ} [and hence *H*(λ = 0)], the evolved state *U*∣Ψ_{i}〉 matches the target state *U*_{adia}∣Ψ_{i}〉 up to a phase factor. For the case of the evolution along the geodesic (e.g., θ = π/2) and *k* = 1,2, …), we had *U*′ ∝ *I* and therefore a perfect population transfer for the initial eigenstates of σ* _{x}*.

### Numerical simulations

In the simulations, we modeled dephasing noise and random fluctuations by adding them to the Hamiltonian (Eq. 8) via δ(λ) → δ(λ) + δ_{0} and Ω(λ) → Ω(λ)(1 + δ_{1}). Here, δ_{0} is the dephasing noise from static and time-dependent magnetic field fluctuations with a _{1} is the random static changes in the driving amplitude. δ_{0} follows the Gaussian distribution with the mean value μ = 0 and the SD σ = 2π × 130 kHz. The probability density of δ_{1} has the Lorentz form *f*(δ_{1}, γ) = 1/{πγ[1 + (δ_{1}/γ)^{2}]} with γ = 0.0067. All the parameters in the distribution function were extracted from fitting the free induction decay and the decay of Rabi oscillation.

### Experimental setup

The experiments were performed with a home-built optically detected magnetic resonance platform, which consists of a confocal microscope and an MW synthesizer (fig. S5). A solid-state green laser with a wavelength of 532 nm was used for initializing and reading out the NV spin state. The light beam was focused on the NV center through an oil immersion objective (numerical aperture, 1.4). The emitted fluorescence from the NV center was collected by a single-photon counting module (avalanche photodiode). Here, we used an NV center embedded in a room-temperature bulk diamond grown by chemical vapor deposition with [100] faces. It has ^{13}C isotope of natural abundance and nitrogen impurity less than 5 parts per billion. To lift the degeneracy of the ∣*m*_{s} = ± 1〉 states, a static magnetic field of 510 G was provided by a permanent magnet. The magnetic field was aligned by adjusting the three-dimensional positioning stage on which the magnet was mounted and by simultaneously monitoring the counts of the NV center. The direction of the magnetic field was well aligned when the counts showed no difference between with and without the magnet. Manipulation of the NV center was performed by MW pulses applied through a home-made coplanar waveguide. The MW pulses were generated by the I/Q modulation of the Agilent arbitrary wave generator (AWG) 81180A and the vector signal generator (VSG) E8267D and then amplified by the Mini Circuits ZHL-30W-252+. An atomic clock was used to synchronize the timing of the two. The AWG supplies the I and Q data with a frequency of 400 MHz, and the VSG generates the 3898-MHz carrier. The output frequency is 4298 MHz, which matches the transition frequency between the NV *m*_{s} = 0 and *m*_{s} = + 1 states.

### Experimental sequences

As the magnetic field is 510 G, we first applied the green laser for 3 μs to initialize the NV center electronic spin to the level of *m*_{s} = 0 and to polarize the adjacent ^{14}N nuclear spin simultaneously (*38*). The preparation of the NV electron spin in an equal superposition state of *m*_{s} = 0 and *m*_{s} = 1 was realized by applying an MW π* _{x}*/2 (π

*/2) pulse, i.e., by the rotation around the*

_{y}*x*(

*y*) axis with an angle of π/2. Then, the NV electron spin was driven according to a desired path. To experimentally characterize the evolution path, we sampled the path with several points and measured the spin state through tomography. π

*/2 or π*

_{x}*/2 pulses were applied to read out the off-diagonal terms. Last, the spin state was read out by applying the laser pulse again and measuring the spin-dependent fluorescence. Typically, the whole sequence was repeated 10*

_{y}^{5}times to get a better signal-to-noise ratio. The schematic diagram of the pulse sequence is shown in fig. S6.

In driving the NV electron spin along the path given in Eq. 5, we used an on-resonant MW field and swept the MW phase θ_{g}λ with the path parameter λ. In driving the NV electron spin along the path of the LZ Hamiltonian (Eq. 7), the MW phase was a constant, Δ was set by the Rabi frequency, and *B _{z}* was the MW frequency detuning that varied as

*B*= − Δcot(θ

_{z}_{g}λ). For continuous driving, the path parameter λ varies with a constant rate

*d*λ/

*dt*=

*f*

_{rot}. In the jumping protocol, λ jumps from point to point: λ = λ

*= (2*

_{j}*N*)

^{−1}(2

*j*− 1) with

*j*= 1,2, … ,

*N*. In this work, the jumping protocol had a constant driving Rabi frequency Ω

_{0}and λ = λ

*if (*

_{j}*j*− 1)

*T*/

*N*≤

*t*<

*jT*/

*N*for a path with

*N*pulses applied in a total time

*T*. In the experiments with the back-forward motion along the geodesic, we reversed the order of the parameter λ in the backward path, that is, in the jumping protocol, we repeated the subsequent parameters (λ

_{1}, λ

_{2}, … , λ

_{N − 1}, λ

*, λ*

_{N}*, λ*

_{N}_{N − 1}, … , λ

_{2}, λ

_{1}), while for the standard protocol of continuous driving, we used the rate

*d*λ/

*dt*=

*f*

_{rot}for a forward path and the rate

*d*λ/

*dt*= −

*f*

_{rot}for a backward path and repeated the process.

We removed the irrelevant dynamic phases if the initial state was not prepared in an initial eigenstate to reveal the geometric evolution. At the beginning of state readout, we compensated the dynamic phases by applying an additional driving with an MW π phase shift (i.e., Ω → − Ω) at the point of the target state for a time equalling to the time for adiabatic evolution. This additional driving did not change the geometric phases and state transfer because it was applied at the final path point.

## SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/5/6/eaax3800/DC1

Fig. S1. Transition of adiabatic control from a constant gap to a gap with energy level crossings.

Fig. S2. Robustness of the jumping protocol against different kinds of control noise.

Fig. S3. Enhancing the robustness of jumping protocol in the presence of control noise.

Fig. S4. Coherence protection during adiabatic evolutions.

Fig. S5. Sketch of the experimental setup.

Fig. S6. Experimental pulse sequence.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is **not** for commercial advantage and provided the original work is properly cited.

## REFERENCES AND NOTES

**Acknowledgments:**

**Funding:**The authors at USTC were supported by the National Key Research and Development Program of China (grant nos. 2018YFA0306600, 2013CB921800, 2016YFA0502400, and 2017YFA0305000), the NNSFC (grant nos. 81788101, 11227901, 11722544, 91636217, and 11775209), the CAS (grant nos. GJJSTD20170001 and QYZDY-SSW-SLH004), the Anhui Initiative in Quantum Information Technologies (grant no. AHY050000), the CEBioM, the IPDFHCPST (grant no. 2017FXCX005), the Fundamental Research Funds for the Central Universities, and the Thousand Young Talents Program. The authors of Ulm University were supported by the ERC Synergy grant BioQ, the EU projects HYPERDIAMOND and AsteriQs, the BMBF projects NanoSpin and DiaPol, and the IQST that is financially supported by the Ministry of Science, Research and Arts Baden-Württemberg.

**Author contributions:**J.D. supervised the entire experiment. Z.-Y.W. and M.B.P. formulated the theory. J.D. and F.S. designed the experiments. K.X., X.X., and P.W. prepared the setup. K.X. and T.X. performed the experiments. T.X., Z.-Y.W., and K.X. performed the simulation. Z.-Y.W., K.X., T.X., M.B.P., Y.W., and J.D. wrote the manuscript. All authors discussed the results and commented on the manuscript.

**Competing interests:**All authors declare that they have no competing interests.

**Data and materials availability:**All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

- Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).