Observation of coherent oscillation in single-passage Landau-Zener transitions.

Landau-Zener transition (LZT) has been explored in a variety of physical systems for coherent population transfer between different quantum states. In recent years, there have been various proposals for applying LZT to quantum information processing because when compared to the methods using ac pulse for coherent population transfer, protocols based on LZT are less sensitive to timing errors. However, the effect of finite range of qubit energy available to LZT based state control operations has not been thoroughly examined. In this work, we show that using the well-known Landau-Zener formula in the vicinity of an avoided energy-level crossing will cause considerable errors due to coherent oscillation of the transition probability in a single-passage LZT experiment. The data agree well with the numerical simulations which take the transient dynamics of LZT into account. These results not only provide a closer view on the issue of finite-time LZT but also shed light on its effects on the quantum state manipulation.

Landau-Zener transition (LZT) has broad applications in atomic and molecular physics, quantum optics, condensed matter physics, chemical physics, and quantum information science. For example, LZT has been applied to investigating the jump time and quantum Zeno and anti-Zeno effects of cold atoms in accelerated optical lattices12, the behavior of molecular magnets at low temperature34, nonequilibrium phase transitions5, and it is also exploited as a tunable beam splitter of wave functions to generate entangled multipartite states67. LZT also plays a key role in determining whether random optimization problems can be solved using the quantum adiabatic algorithm8. Recently, LZT's potential for robust manipulation of coherent quantum states has attracted much attention in the context of quantum information processing679101112131415161718 because LZT may provide a simple and effective solution to the realization of high fidelity quantum state control without the need for precise timing.

The time-dependent Hamiltonian describing LZT in quantum two-level systems can be written in the generic form aswhere σ, are Pauli matrices, ε(t) = vt is the energy difference between the two diabatic (crossing) basis states (i.e., the eigenstates | ↑〉 and | ↓〉 of the σ operator) controlled by an external parameter which depends linearly on time t, and Δ is the constant gap between the two instantaneous eigenenergy states |+〉 and |−〉 at the center of the avoided crossing ε = 0, as depicted in Fig. 1(a). In such systems, when ε(t) is swept through the avoided crossing, transitions between |±〉 with energies can occur and the transition probability is given by the well-known Landau-Zener (LZ) formulawhere v = |dε/dt| is the Landau-Zener speed and we have set the reduced Planck constant . Equation (2) gives the probability of finding the system in the excited (ground) state at ε = ε(t → +∞) when it is started in the ground (excited) state at ε = ε(t → −∞). By defining α = Δ2/4v as the adiabaticity parameter the LZ formula can be simplified to P = exp(−2πα). Although analytical solution to the problem cannot be obtained when ε and/or ε are finite, it is well known that for , the LZ formula provides an excellent approximation to the actual transition probability and P↓ ≈ P. However, when is not satisfied, the LZ formula may become quantitatively inaccurate or even qualitatively incorrect. In spite of some theoretical studies on the effects of finite |ε,|/Δ on P↓, there is an acute lack of adequate experimental evidence.

Figure 1

Circuit and experimental procedure.

(a) A general avoided energy-level crossing with diabatic basis states (the dashed lines) and adiabatic basis states (the solid lines). The constant gap Δ between the two instantaneous eigenenergy states at ε = 0 is the tunneling amplitude, i.e., the coupling strength. (b) Optical micrograph of the sample with Al/AlO/Al Josephson junctions on the silicon substrate. MW, Φ and RO represent the microwave pulse to prepare the initial state, the sweeping flux bias to induce LZT and the 5-ns flux bias pulse to tilt the potential well in order to readout the qubit state, respectively. The planar coils and dc-SQUID magnetometer are coupled inductively to the qubit. (c) A typical time profile of the manipulation and measurement waveforms employed to perform the single-passage Landau-Zener experiment.

On the other hand, understanding LZT with moderate values of |ε,|/Δ is in urgent need because this region of parameter space is important to quantum information processing. For instance, in superconducting qubits the tuning range of energy level spacing is usually limited to a couple of GHz or even as narrow as a few hundreds of MHz while Δ/2π could be as large as 102 MHz671920212223242526. For quantum state control based on sweeping ε through avoided crossings, understanding LZT probability's dependence on |ε,|/Δ and the sweeping time is essential to high fidelity operation. The fidelity of various techniques based on LZT relies critically on the accuracy of the LZ formula which predicts a simple exponential dependence of P on the adiabaticity parameter α only. Therefore, for a qubit starting from the ground (excited) state the probability of finding it in the excited (ground) state after a single passage through the avoided crossing is assumed to be determined entirely by α (i.e., Δ2/4v) but not the detail of the process such as |ε,|/Δ. Here, using an artificial atom — a superconducting phase qubit — coupled to a microscopic two-level system (TLS), we test the accuracy of the LZ formula in the region of |ε,|/Δ < 4.3. We show that in contrast to conventional wisdom, in the region of parameter space most relevant to superconducting qubits, P↓ could deviate significantly from P determined by the LZ formula. Our experiment and numerical simulation demonstrate P↓ can oscillate coherently as a function of ε for constant α when |ε,| is comparable to Δ, which is named as coherent Landau-Zener oscillation (LZO).

Results

In our experiment we use a superconducting phase qubit. However, since a single phase qubit does not have an intrinsic avoided energy-level crossing, we utilize an avoided level crossing arising from interaction between the qubit and a microscopic TLS2728. As discussed below in more detail, when the transition frequency of the qubit ω10 is close to that of the TLS ω, which is fixed, the first and second excited states of the coupled qubit-TLS system form an effective quantum two-level system described by the LZ Hamiltonian (1). Note that to make quantitative comparisons between the theory/numerical simulation and the experiment without free parameters, all relevant system parameters, including the energy relaxation and dephasing time of the qubit and the energy gap Δ, are obtained from direct measurements.

A microscopic picture of the superconducting phase qubit is shown in Fig. 1(b). The qubit consists of an L ≈ 770 pH superconducting loop intersected by a Al/AlO/Al Josephson tunnel junction with a critical current I0 ≈ 1.4 μA and a junction capacitance C ≈ 240 fF. By varying the magnetic flux applied to the superconducting loop the potential energy of the qubit becomes asymmetrical. The ground state and the first excited state in the upper potential well, represented by |0〉 and |1〉 respectively, can be used as the computational basis states of the qubit. For an isolated qubit, the transition frequency between |0〉 and |1〉, ω10, is a single-valued function of the external flux bias Φ which is coupled inductively to the superconducting loop through an on-chip flux bias line.

As shown in Fig. 2(a), however, the microwave spectrum of the qubit ω10(Φ) has a rather large avoided energy-level crossing at Φ ≈ −2.8 mΦ0 (with respect to the flux bias point at which ω10/2π ≈ 16.348 GHz) indicating significant interaction between the qubit and a microscopic TLS30. The transition frequency between the TLS' ground state |g〉 and excited state |e〉 and the qubit-TLS coupling strength are ω/2π = 16.450 ± 0.002 GHz and Δ/2π = 70.0 ± 0.5 MHz from the spectrum and vacuum Rabi oscillation, respectively. Note that in this coupled qubit-TLS system the time-dependent energy difference between the two diabatic states involved in LZT is ε(t) = ω10(t) − ω which depends linearly on the flux bias to a good approximation. The relationship between the flux bias and ε can be found from Fig. 2(a).

Figure 2

Spectroscopy and Rabi oscillation.

(a) Microwave spectroscopy measurement of the coupled qubit-TLS system. The splitting is Δ/2π = 70.0 ± 0.5 MHz centered at ω/2π = 16.450 ± 0.002 GHz. The beginning point of the flux bias for the single passage LZ sweeping is denoted as 0 mΦ0, corresponding to ε/2π = (ω10 − ω)/2π ≈ −100 MHz in the upper abscissa. (b) Rabi oscillation and T1 at ε, respectively. The experimental data (the red circles and blue triangles) agree well with the theoretical fits (the black solid lines).

Fig. 1(c) illustrates the experimental procedure used to observe coherent LZO. We begin by setting the initial diabatic energy of the effective quantum two-level system ε at about 100 MHz below ω with a static flux bias. The qubit is prepared in its ground state by waiting for much longer than the energy relaxation time T1 ≈ 70 ns of the qubit. A microwave pulse is then applied to the qubit when it is biased at a fixed value ε/2π ≈ −100 MHz. The microwave pulse coherently transfers the population of the qubit-TLS from |0g〉 to one of the system's eigenstates |−〉 through a process that is discussed in detail in Methods. The lack of oscillation in T1 measurement taken at ε as shown in Fig. 2(b) confirms that the initial state of the qubit-TLS system at t = 0 is indeed the eigenstate |−〉. As illustrated in Fig. 1(c), a time-dependent flux Φ(t) = Φ/t is then superimposed between t = 0 and t = t onto the static flux bias to sweep ε linearly from ~ −100 MHz to its maximum value ε. The corresponding LZ speed v is thus (ε − ε)/t. This is followed immediately by a 5-ns readout pulse which performs a projective measurement of the probability P↓ of finding the qubit in state |1〉 (i.e., the coupled system is in state |1g〉 corresponding to | ↓〉 in Fig. 1(a)).

We first measure P↓ vs. t at a constant value of ε by keeping Φ fixed while increasing t from almost 0 ns to 45 ns. The maximum t is selected to avoid too much influence of the qubit's energy relaxation. By stepping Φ from 0 to −11 mΦ0 the value of corresponding ε is then varied from about −1.4Δ to ~ 4Δ, in which the condition is no longer satisfied. This procedure is repeated at each ε to obtain P↓(ε, t). Fig. 3(a) shows the dependence of P↓ on ε and t. It can be seen that P↓ vs. t decays exponentially for ε < 0 (Φ ∈ [−2.8 mΦ0, 0]) with a characteristic time T1 due to energy relaxation. As ε becomes positive, P↓ vs. t becomes oscillatory. Since the avoided crossing is traversed only once, the observed oscillation in P↓ vs. t with constant ε must be a consequence of the moderate value of ε/Δ and is not caused by the Landau-Zener-Stückelberg interference which requires multiple passages through the avoided crossing.

Figure 3

Coherent LZO.

(a) Experimentally measured P↓ vs. ε and t. (b) Numerically calculated P↓ vs. ε and t with all input parameters obtained from the experiment. The white dashed lines correspond to the value of modified adiabaticity parameter α′ = 10. Notice that in the region below (above) the lower (upper) white dashed line, one has α′ ≫ 1 thus the system evolves adiabatically and no oscillation in P↓ is expected as confirmed experimentally. The LZ speed v equals the slope of any straight lines originated from the lower-left corner of the t-ε plane. For example, the yellow dashed line in (a) has v = 400/19.5 ≈ 20.5 MHz/ns. For the sake of clarity, the temporal evolution of the system along the yellow dashed line (constant v) is presented separately in Fig. 4(b).

It is worth noting that the observed oscillation is not a consequence of ill-prepared initial states with non-negligible probability amplitude in the excited state |+〉 of the effective Hamiltonian (1) because the microwave pulse used to initialize the system resonantly couples |0g〉 to |−〉 and has negligible coupling to |+〉 due to large frequency detuning. Furthermore, this process of transferring the system to the desired initial state |−〉 via a resonant microwave pulse is robust in the sense that it does not depend sensitively on the accuracy of the pulse duration t. Deviation in pulse duration simply leaves some probability amplitude in |0g〉 which has no effect on LZT other than reducing the visibility of the oscillation (see Methods for detail on the initial state preparation). Therefore, we are confident that the oscillation observed in the non-adiabatic region of the parameter space arises neither from Landau-Zener-Stückelberg interference nor unwanted probability amplitude of |+〉 in the initial state. This is also supported by the good agreement between the results of experiment and numerical simulation shown in Fig. 3(b), which uses |−〉 as the initial state at the start of the single passage sweep.

Discussion

By replacing ε with ε + vt, solving the problem of sweeping ε in a finite range is transformed to finding P↓ at finite time. Previous studies have discovered that LZ transition probability reaches the asymptotic value given by equation (2) at (assume ε = 0 at t = 0), where is called the Landau-Zener time, and oscillates in the vicinity of avoided crossing (corresponds to t ≤ t) due to the transient dynamics123132. Since corrections to the standard LZ formula are significant only if the adiabaticity parameter α ≤ 1, the region of ε within which transient dynamics plays an important role is given by , where “≤” means less than or comparable to. By examining the experimental data we find that the region of most noticeable coherent LZO coincides with which agrees well with the result of numerical simulation. These results unambiguously show that coherent LZO is originated from the transient dynamics of the LZT.

Such coherent LZO has little effect on the adiabatic evolution. Because in the true adiabatic regime, by definition the system always stays in the instantaneous ground state and no LZT could occur. In order to find the region of approximate adiabatic evolution in our experiment, it is necessary to modify the definition of the adiabaticity parameter to . For the adiabatic theorem to hold is required. As shown in Fig. 3(a), the white dashed lines represent α′ = 10. It is clear that there is no coherent LZO observed in the region .

In a popular analogy to optics an avoided crossing acts as an effective beam splitter, with a transmission coefficient corresponding to P in the LZ formula, for quantum wave functions. This beam splitter analogy has been applied successfully to the visualization and explanation of the behavior of superconducting and semiconductor qubits69192033343536. In this analogy, a single sweep through the avoided crossing is equivalent to passing a beam of light through the beam splitter only once. When , P↓ ≈ P and thus a greater LZ speed corresponds to a higher transmission coefficient of the beam splitter according to the LZ formula. But when is not satisfied, P↓ differs greatly from P. As an example, P↓ vs. t, and thus the LZ speed v, with ε/2π = 200 MHz is shown in Fig. 4(a). The maximum in the difference δP↓ between the experimental P↓ and those obtained from the LZ formula (2), shown in the inset of Fig. 4(a), can reach 0.21. The observation of coherent LZO strongly suggests that when is not met corrections to the LZ formula should be considered to avoid conceptual difficulties.

Figure 4

P↓ oscillation at constant LZ speed.

(a) Measured P↓ as a function of t (the red squares) with ε/2π = 200 MHz which clearly shows oscillation in the region of α < 1 which compares well with the result of numerical calculation (the solid line). This is in stark contrast to the smooth exponential decay expected from the LZ formula (the dashed line). The inset is the difference δP↓ = P↓ − P. (b) The measured (the red squares) and numerically calculated (the solid line) P↓ vs. t with constant LZ speed v ≈ 20.5 MHz/ns corresponding to evolving along the yellow dashed line in Fig. 3(a). Again, P↓ oscillates in the region where the adiabatic condition α > 1 is not satisfied which is not expected from the LZ formula (the dashed line).

Coherent LZO also has significant consequences on the coherent manipulation of quantum states of single qubits and coupled two-qubit systems based on LZT1012. For this approach of quantum state control, the LZ transition probability P plays a central role since each single passage through the avoided crossing results in a unitary operation U given by Ref. 12where φ is the Stokes phase37, which has no effect on the single-passage LZT process discussed here and thus can be set to zero for the sake of convenience. As mentioned above, the transition frequency ω10 of most artificial atoms, in particular the superconducting qubits, is limited to a couple of GHz. Because the speed of two-qubit operations is proportional to the inter-qubit coupling strength Δ, increasing |ε,|/Δ by reducing Δ is undesirable. Hence, evaluating U according to the LZ formula (2) could result in significant errors when is not satisfied. In order to conduct a quantitative analysis, the experimental data along the yellow dashed line corresponding to constant v, in the Fig. 3(a) are extracted and shown in Fig. 4(b). Oscillation in P↓ is clearly observed and it is qualitatively different from the exponential decay predicted by equation (2), when decoherence is taken into account. Suppose the initial quantum state is | ↓〉. Then after a single passage through the avoided crossing, if one replaces P with P↓ as in the asymptotic situation, the deviation δP↓ = P↓ − P would be quite large. For example, when t = 12.8 ns the deviation δP↓ = 0.229, which is unacceptably large for coherent quantum state transformation.

In conclusion, we have investigated the effect of finite energy (ε) sweep (or equivalently finite time) on LZT probability P↓ experimentally. Single-passage technique is used to isolate the effect of finite ε on P↓ from that of interference caused by passing the avoided crossing multiple times. We find that P↓(ε/Δ, α = const) oscillates when ε is comparable to Δ and α < 1. The good agreement between the experiment and numerical calculation strongly supports the notion that coherent LZO is caused by the underlying transient dynamics of the finite time LZT which cannot be described by the LZ asymptotic formula. In this region of the LZT parameter space, corrections to the LZ formula must be taken into account, otherwise it will lead to substantial errors in quantum state operations based on LZT. The result also shows that when applying the simple beam splitter analogy one should not automatically assume that greater α (i.e., faster sweep) corresponds to larger transmission coefficient (i.e., greater P) as implied by the asymptotic LZ formula.

Methods

Initial state preparation

We first derive an analytical result explaining the lack of oscillation at the very beginning of LZT. The Hamiltonian of the qubit-TLS system coupled to a microwave field is given as: (in the basis {|0g〉, |1g〉, |0e〉, |1e〉})where Ω is the Rabi frequency, ω is the microwave frequency, ω is the energy difference between the ground state |g〉 and the excited state |e〉 of TLS, Δ is the coupling strength between the qubit and TLS, δ = ω10 − ω and δ = ω − ω10 are detunings. By rotating the frame, Hamiltonian (3) can be transformed to the following time-independent form2930Next, we rewrite H1 in which the subspace spanned by {|1g〉, |0e〉} is diagonalized:where . The basis for H2 is now {|0g〉, |−〉, |+〉, |1e〉}. Note that in our experiment, before turning on the microwave the state is at Ψ(t = 0) = |0g〉. By turning on the microwave (Ω ≠ 0), |0g〉 is coupled to both |−〉 and |+〉. The resonance between |0g〉 and |±〉 occurs when , from which we obtain the resonant condition:In the limit of , we have ω = ω10, which corresponds to the usual two-state Rabi oscillation. Note that in this limit there is also a solution ω = ω. However, in this case the coupling strength is . The reason is that although the microwave frequency could match that of TLS, coupling between the microwave and TLS is negligible which is confirmed by the absence of Rabi oscillation between the two states of the TLS in a separate experiment. In the other limit of ω10 − ω = 0, we have , and the dynamics have been thoroughly studied in Ref. 30.

In our experiment, we have (ω − ω10)/2π ≈ 100 MHz and Δ/2π ≈ 70 MHz, which means LZT occurs in the region where (ω − ω10) ~ Δ. Because the frequency of the applied microwave is ω = λ−, which can be determined from the measured energy spectrum shown in Fig. 2(a), |0g〉 is resonantly coupled to |−〉, which is the eigenstate of . Although there is in principle also a coupling between |0g〉 and |+〉, the effective coupling is much smaller because of the large detuning, as discussed below.

For Δ/2π ≈ 70 MHz, (ω − ω10)/2π ≈ 100 MHz, the resonance between |0g〉 and |−〉 occurs at , we obtain and . In our experiments, the coupling strength between |0g〉 and |−〉 is about 20 MHz and that between |0g〉 and |+〉 would be . Because 6.3 MHz is comparable with 20 MHz, one may think that coupling between |0g〉 and |+〉 cannot be neglected. However, there is also a large detuning of about 122 MHz between |0g〉 and |+〉. Therefore, the effective coupling between |0g〉 and |+〉 is reduced to (6.32/122) ≈ 0.33 MHz and thus can be safely neglected. To be more precise, we calculated the population P± (where P± is the population of state |±〉) after the application of a π pulse numerically, and it is found that . Based on this analysis, when ω = λ−, the dynamics can be described by the Hamiltonian in the subspace {|0g〉, |−〉}:At the resonance , H3 becomeswhere I is a 2 × 2 identity matrix.

For initial state Ψ(t = 0) = |0g〉, the amplitude of |−〉 isConsidering , the amplitude of |1〉 is

Therefore, with a microwave pulse of duration t which is used to prepare the initial state, we haveThis is the reason why in the experiment we observe a usual Rabi oscillation instead of Rabi beating30 which is indicated by the red circles, as shown in Fig. 2(b). In addition, when microwave is turned off, the subspace {|+〉, |−〉} is isolated from |0g〉 and |1e〉. Projected into the subspace {|+〉, |−〉}, the system stays in the eigenstate |−〉. This explains why we observe a monotone decay of P↓ with no oscillations, as indicated by the blue triangles in Fig. 2(b).

Effect of t and ε/Δ on LZT

In this section, we discuss two factors that may affect the LZT probability, i.e., the width of the microwave pulse t used to prepare the initial state at ε and the end of the normalized diabatic energy sweeping ε/Δ, respectively.

After a microwave pulse, by projecting into the subspace {|1g〉, |0e〉}, the system is in the eigenstate |−〉. Then the dynamics of LZT can be described by H with a time-dependent ω10(t), i.e.,To investigate the Landau-Zener diffraction effect, we sweep ω10(t) across ω, i.e.,When , we expect that the Landau-Zener asymptotic formula holds and P↓ ≈ P. Thus no oscillations should occur in P↓. This is confirmed by the result of numerical simulation shown in Fig. 5(c), where ε/2π = 1450 MHz ≈ 20.7Δ, reproducing the exponential decay behavior described by the asymptotic LZ formula independent of the initial state of the qubit-TLS system. In this case, the population in the qubit state |1〉 can be expressed aswhere the first term reflects the effect of microwave duration t in preparing the initial state, the second term corresponds to the LZT probability, and the third term represents the relaxation effect. However, as ω10(t = t) moves towards ω, the situation does not hold any more, and we observe oscillation features in the t direction, as shown in Fig. 5(a) (experiment) and 5(b) (numerical simulation). Notice that Fig. 5(a) and 5(b) also confirm that the effect of imprecise π pulse is an incomplete transfer of system from |0g〉 to |−〉, which reduces the probability amplitude of |−〉 from the maximum value, instead of resulting in non-negligible probability amplitude in the unwanted |+〉 state.

Figure 5

The effect of t on finite time LZT (a) Experimentally measured and (b) numerically calculated P↓ vs. t and t for ε/2π = 200 MHz showing that the effect of imprecise π pulse is to reduce the visibility of the oscillation of P↓ vs. t by reducing the probability amplitude of the desired |−〉 state. Because the microwave pulse is resonant with the |0g〉 ↔ |−〉 transition while largely detuned from the |0g〉 ↔ |+〉 transition even a significant deviation from a π pulse would only result in negligible transfer of population to |+〉. Furthermore, since |0g〉 does not participate in the single passage LZ process, the observed oscillation could neither be due to LZS interference nor non-negligible population in |+〉 at the beginning of each ε sweep. For comparison, we also present the numerically calculated P↓(t, t) for ε/2π = 1450 ≈ 20.7Δ in (c). The result shows the exponential decay behavior described by the asymptotic LZ formula as expected for .

Author Contributions

G.S. and S.H. conceived the experiments; G.S. carried out the measurements and analyzed the data with the help of X.W., M.G., D.Z., Y.Y., S.Z., J.C., P.W. and S.H.; X.W. performed the numerical calculations; G.S. and S.H. wrote the paper.