arxiv: v1 [quant-ph] 13 Jul 2018

Transcription

1 Fast Amplification and Rephasing of Entangled Cat States in a Qubit-Oscillator System Z. Xiao, 1,, T. Fuse,, S. Ashhab, 3 F. Yoshihara, K. Semba, M. Sasaki, M. Takeoka, 1,,, 5 and J. P. Dowling 1 Hearne Institute for Theoretical Physics, and Department of Physics & Astronomy, Louisiana State University, Baton Rouge, Louisiana 783, USA National Institute of Information and Communications Technology, --1, Nukui-Kitamachi, Koganei, Tokyo , Japan 3 Qatar Environment and Energy Research Institute, Hamad Bin Khalifa University, Qatar Foundation, Doha, Qatar NYU-ECNU Institute of Physics at NYU Shanghai, Shanghai 6, China 5 CAS-Alibaba Quantum Computing Laboratory, USTC, Shanghai 1315, China arxiv:187.97v1 quant-ph 13 Jul 18 We study a qubit-oscillator system, with a time-dependent coupling coefficient, and present a scheme for generating entangled Schrödinger-cat states with large mean photon numbers and also a scheme that protects the cat states against dephasing caused by the nonlinearity in the system. We focus on the case where the qubit frequency is small compared to the oscillator frequency. We first present the exact quantum state evolution in the limit of infinitesimal qubit frequency. We then analyze the first-order effect of the nonzero qubit frequency. Our scheme works for a wide range of coupling strength values, including the recently achieved deep-strong-coupling regime. Introduction. The interaction of a two-level atom qubit with a quantized field oscillator has been widely studied over the past few decades. There have been numerous experimental realizations of such systems, including superconducting circuits 1 9, and systems of atoms coupled to superconducting microcavities 1. Mathematically, such qubit-oscillator systems, described by the quantum Rabi model, can be simplified to the Jaynes- Cummings model, with the rotating-wave approximation RWA, in the case of a near-resonant field 11, 1. We focus on the case where the qubit frequency is small compared to that of the oscillator, where the RWA cannot be used. Previous studies have mainly focused on time-independent coupling coefficients 13, 1. In this work, we examine the qubit-oscillator system with a time-dependent coupling coefficient, and we demonstrate a scheme to generate Schrödinger-cat states of a very large size, and we develop a second scheme that also protects the cat states from dephasing. Our amplification scheme offers a potentially simple and fast alternative compared to previous methods 15, 16. We present the state evolution in the exact analytical form in the limit of an infinitesimal qubit frequency, and then we examine the first-order effect of a small nonzero qubit frequency. State evolution under time-dependent coupling with infinitesimal qubit frequency. The Hamiltonian of the combined system of the qubit and the oscillator is Ĥt= ˆ + â â gtˆσ x â + â, 1 where and are the frequencies of the oscillator and the qubit, respectively, and gt is the time-dependent coupling constant, â and â are, respectively, the annihilation and creation operators of the oscillator, and ˆσ x,z are the Pauli operator of the qubit 17. We focus on the situation where is small compared to. We will first examine the zeroth order effect of the small by taking the limit of /, which gives Ĥ t= â â gtˆσ x â + â. The initial eigenstates of the Hamiltonian at t= are the entangled states 13, 1: E = + ˆD x g ± x ˆD + g, 3 ignoring a factor of 1/, and where are the two qubit eigenstates of the Pauli matrix ˆσ x with eigenvalues ±1, the is the N-photon Fock state in the oscillator, and ˆD ± g are displacement operators. Note that, in the limit of /, the energy eigenstates E N+ and E N are almost degenerate, with the energy eigenvalues of E = N + 1/ g /. We now consider the state evolution under an arbitrary time-dependent coefficient gt. The energy eigenstates E t, determined by the instantaneous value of gt, do not reflect the evolution of quantum states. In other words, an initial state E generally does not, for a time-dependent gt, evolve into e iφt E t, with a phase factor φ t, at a later time t. On the other hand, a quantum state in the form of Ẽ evolves into the quantum state e iφ N t Ẽ Ẽ t = + x ˆD t, where gt ± x ˆD + gt, and the phase factor φ N t is given in the Supplemental Material. Here the complex variable gt obeys the

2 a gt/g c gt/g, gt/g i ii iii time π/ gt/g gt/g time π/ b Im gt/g d Imgt/g, Im gt/g Re gt/g Regt/g, Re gt/g FIG. 1. Evolution of gt under time-dependent gt with infinitesimal. a Three different time dependences of gt. b The traectories of gt in the complex plane corresponding to a, which, as seen in Eq., determine the evolution of the dynamical evolution eigenstates. The initial time is t= and the traectories start from g/g=1. c The time dependences of gt and gt, in the case of sinusoidal driving force gt=g1 + cos t/. d The traectories of gt and gt corresponding to c, in the complex plane. The amplitude of gt keeps expanding, showing cat-state amplification on resonance. The blue and red dots in c and d indicates where the modulation stops t=π/ in the scheme shown in Fig.. equation gt=igt gt. 5 Note that the initial g can be set to any value. The proof is given in the Supplemental Material. We refer to Ẽ t as the dynamical evolution eigenstates, which take the similar form with the energy eigenstates in Eq. 3, but with gt replaced with gt. Now the dynamics of the system is governed by Eq. 5, which shows that gt does not respond instantly to changes in gt. Note that, apart from a few special cases, the dynamical evolution eigenstates are generally not energy eigenstates. To understand the evolution of quantum states, we give the following three scenarios, in which we set the initial g=g so that Ẽ = E. i Suppose gt is adiabatically changed over time red curve in Fig. 1a. Any energy eigenstate has ample time adust to the adiabatically changing Ĥ t and also remains an energy eigenstate red curve in Fig. 1b. ii Now suppose gt is constant at a certain value at t<t, and then set to zero instantaneously at t=t blue curve in Fig. 1a. Since neither nor is infinitely large, the states Ẽ t cannot adust instantaneously, and they remain the same at t=t +. However Ẽ t + are no longer the energy eigenstates, and the states begin to evolve. Taking the ground dynamical evolution eigenstate Ẽ ± t = + x gt/ ± x + gt/, 6 as an example, the amplitude of the coherent state component of the state, ± gt before the adustment, should begin to revolve around the origin after the adustment, consistent with the evolution of a regular coherent state in a free oscillator blue circle in Fig. 1b. When gt is instantaneously set to a nonzero value, gt revolves around this value in the complex plane. iii Now we consider the intermediate scenario. We assume that gt is adusted over a finite period of time to zero and then kept stabilized as shown in the green curve in Fig. 1a. In this scenario, gt will start changing as gt starts changing. Its traectory is less intuitive than the extreme scenarios, but can be understood from Eq. 5. After gt becomes constant again, gt evolves in circular motion around the new constant g. A very useful application of varying g through time is to amplify the absolute value of the amplitude of the entangled-cat-state components of Eq. 6. We can achieve this without having to make the absolute value of g large, which is difficult in real experiments, but rather by modulating g periodically on resonance with the frequency. A somewhat similar modulation has been proposed to control the coupling of two flux qubits via a quantum bus 18. Generally, any modulation of gt on resonance with the oscillator frequency leads to a steady increase in the coherent-state amplitude in Eq. 6. As a specific example, the case of a sinusoidal driving function gt=g1 + cos t/ is shown in Figs. 1c and d. The magnitude gt will grow linearly with time. This behavior is easy to understand from Eq. 5, whose solution corresponds to a simple harmonic oscillator being driven by an external force. Since the photon number in the coherent state is proportional to g/, it will grow quadratically as a function of time. Simply growing for two oscillator periods increases the absolute amplitude of the coherent state component by a factor of 3.3. So we will focus on two oscillator periods of growth in later discussions. First-order effect of nonzero-qubit frequency. So far we have ignored the effect of the small by taking the limit of /. We now examine the first-order effect in in the full Hamiltonian in Eq. 1. Note that for a qubit to avoid effects from surrounding noise and unwanted excitations by the pulse operations, has to be reasonably large. For example, in recent experiments 6, 7, /.1. At any time t we can express a general state of in-

3 3 terest as a superposition of dynamical evolution eigenstates: ϕt = N,± C int te Ẽ t. As the state ϕt evolves, so does every dynamical evolution eigenstate Ẽ t and its corresponding coefficient C t. Under the full Hamiltonian Ĥt, the C t generally change over time. If we consider up to the first order in /, C t can be expressed as: C+ C C C 1... { i t exp tm tdt}, where M σz t is the matrix int of the operator ˆ in the basis e Ẽ t at time t. For generality, we have made the parameter timedependent = t. A detailed derivation of this result is provided in the Supplemental Material. An intuitive way to understand the effect of nonzero is the following. If was to be considered infinitesimal, the quantized oscillator has equally spaced energy levels. Nonzero disrupts such equally spaced energy levels. We have shown above that nonzero leads to the change in C t causing any general quantum state to dephase, which we will illustrate below. To minimize this dephasing, we need to make t tm tdt as close to zero as possible. This type of dephasing is well known for coherent state solutions to the harmonic oscillator and was first described by Schrödinger in 196 when he showed how Bohr s correspondence principle applies to the hydrogen atom 19. For large principle quantum number n, the electron should resemble a classical particle in a Kepler orbit around the proton. This never happens with the stationary eigenstates of the hydrogen atom, regardless of how big n is, but in the limit of large n, the energy levels are nearly equally spaced and can be approximated by a harmonic oscillator. Schrödinger then showed that the electron in a coherent state, a particular superposition of different n, moves as a wave packet in a Kepler orbit. However since the levels are not truly equally spaced, dephasing occurs and the electron packet spreads in the orbit until the front of the wave packet wraps around and interferes with the back. The effect has been seen in experiments with single-electron Rydberg atoms. We see exactly the same effect here. With zero the energy levels are exactly equally spaced and the coherent state of the oscillator does not dephase. With small but nonzero the energy levels are not exactly equally spaced and dephasing now occurs, as we show below. Let us examine a special situation where we can actually reduce the term t tm tdt completely to zero. We can take advantage of the fact that as long as g is a constant, M σz t is periodic with a period π. Keeping g a constant, first we keep at a certain nonzero value for a k π period of time, with k being an integer; then we flip to the opposite sign for another k π period of time. As a result, the two parts k π tm σz tdt q q a c p b d p FIG.. Cat-state amplification and rephasing in the Wigner representation of the oscillator state proected onto the qubit state + x + x. The parameters are: =.1; g=.833. a The initial state, which is taken to be the ground state, reasonably resembles the cat state g g + +, with a fidelity of After the initial state goes through two oscillator periods of the sinusoidally driven cat-stateamplification process, the resulting state is shown in b, which is the cat state gπ/ + + gπ/, with g/ = 3.3 and a fidelity of.986. We now let the state in b freely evolve for 1 oscillator periods, and the resulting state of the oscillator is shown in c. The fidelity between the state gπ/ + + gπ/ and the states in c is only.85, since the cat state is dephased by the nonzero. However, if we insert one ˆσ x-π pulse in the middle of the 1 oscillator periods, the evolved state, shown in d, has a.9817 fidelity with the state in b, analogous to a Hahn spin-echo rephasing effect. and k π tm k π σz tdt cancel each other, eliminating the first-order dephasing effect of for the duration of, k π. π-pulses and their application for cat state rephasing. Above we have discussed the method of eliminating the first-order effect of finite by flipping the sign of. In real experiments, manipulating g or directly is not always easy. Usually g and cannot be changed too fast, and the range of their change can be limited. For example, in the recent experiments 6, the qubit and the oscillator are coupled via a Josephson inductances which effectively limits g to being higher than a certain bound. Here we will show a scheme to manipulate g or indirectly. By applying π-pulses to the qubit alone, which is a commonly used technique in dynamical decoupling 1,, we can achieve the effect of flipping the signs of g or. There are three basic types of π-pulses, each of which amount to applying a Pauli operator ˆ, ˆσ x or ˆσ y to the qubit. Let us first examine the ˆ -π-pulse. Since ˆσ x and ˆ anti-commute, we have Ĥgˆ = ˆ Ĥ g in Eq. 1. Therefore, without directly altering the qubit-

4 oscillator coupling coefficient g, in the Hamiltonian the sign of g is flipped by applying ˆ -π pulses to the qubit alone. For ˆσ x and ˆσ y pulses, similar arguments apply. By applying ˆσ x -π pulses to the qubit, the sign of is flipped. By applying ˆσ y -π pulses to the qubit, the signs of both g and are simultaneously flipped. With infinitesimal and the simplified Hamiltonian Ĥ t, when the initial state is the ground state E + t = + gt x + gt x +, modulating gt at the oscillator frequency can increase the amplitude of the coherent state component, as shown in Figs. 1c and d. With nonzero but still small compared to and full Hamiltonian Ĥt, the ground state, which can be writen as E + t = + x ψ + x ψ + is very close to E + t. Although E +t is not strictly a cat state, for convenience we will still refer to the process for increasing the amplitude in the oscillator components in Eq. 6 as entangled cat-state amplification. This can be achieved either by directly modulating gt, or by indirectly flipping the sign of gt by applying ˆ -π pulses to the qubit alone. Here, we modulate gt directly sinusoidally to amplify a cat state. Cat-state amplification. We begin with g constant and start with the ground state E + as the initial state, with the Wigner representation of ψ + ψ + shown in Fig. a. We apply a short sinusoidal driving force by modulating gt= cost after t=, until reaching t= π, when two oscillator periods of catstate amplification is completed, as shown with the red and blue dots in Fig. 1c and d. The Wigner representation of the resulting state, which has a high fidelity with gπ/ + + gπ/, is shown in Fig. b. We can repeat this process to further increase the absolute value of amplitudes of the oscillator states, but we will focus on the case where the cat state amplification lasts for only two oscillator periods. Note that the amplified cat in Fig. b is tilted by the angle about π/. This is due to the fact that gπ/ lags gπ/ by about a quarter of a cycle, as indicated by the dots in Fig. 1c and d. Rephasing the Amplified Cat. Compared to the catstate amplification, which increases the amplitude by a factor of 3.3 within ust the first two oscillator periods, the dephasing effect takes a longer time due to being small compared to. Therefore the dephasing effect of the small is not very concerning if the cat state amplification lasts only for a short time. But if we now let the state shown in Fig. b evolve freely for π/, then the state is eventually dephased, as shown in Fig. c. To counter this, we instead in Fig. d allow the cat to evolve freely for five oscillator periods, apply one σ x π pulse to the qubit, and then allow free evolution for another five periods. We can see from the Wigner representation of the evolved state, shown in Fig. d, that the amplified cat state is recovered in analogy to the rephasing of the Hahn spin-echo method 3. If needed, more ˆσ x -π pulses can be applied in the same manner to preserve a state for a longer time. Also we can apply ˆσ x -π pulses more frequently, to counter the dephasing by a larger. Conclusion. We have analytically solved the equations that describe the evolution of the quantum state in a qubit-oscillator system, under a Hamiltonian with timedependent coupling coefficients and infinitesimal qubit frequency. We have introduced dynamical evolution eigenstates analogous to the energy eigenstates under a time-independent coupling coefficient, and have shown that the dynamical evolution eigenstate tremendously simplifies the expressions for the quantum state evolution. We have also calculated the effect of the nonzero qubit frequency and obtained its first order form. Making use of a sinusoidal driving field and π pulses, we have designed schemes to generate cat states with large amplitudes and rephase a cat state from the effect of nonzero qubit frequency. We point out that our formalism is quite general and can be used for general entangled-cat-state engineering. Acknowledgments. T.F., F.Y., K.S., M.S., and M.T. would like to acknowledge support from Japan Science and Technology Agency Core Research for Evolutionary Science and Technology Grant No. JPMJCR1775. Z.X and J.P.D. would like to acknowledge AFOSR, ARO, DARPA, NSF, and NGAS. xiaozhihao@hotmail.com tfuse@nict.go.p 1 I. Chiorescu, P. Bertet, K. Semba, Y. Nakamura, C. Harmans, and J. Mooi, Nature 31, 159. A. Wallraff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. Huang, J. Maer, S. Kumar, S. M. Girvin, and R. J. Schoelkopf, Nature 31, M. Devoret, S. Girvin, and R. Schoelkopf, Annalen der Physik 16, T. Niemczyk, F. Deppe, H. Huebl, E. Menzel, F. Hocke, M. Schwarz, J. Garcia-Ripoll, D. Zueco, T. Hümmer, E. Solano, et al., Nature Phys. 6, P. Forn-Díaz, J. Lisenfeld, D. Marcos, J. J. García-Ripoll, E. Solano, C. Harmans, and J. Mooi, Phys. Rev. Lett. 15, F. Yoshihara, T. Fuse, S. Ashhab, K. Kakuyanagi, S. Saito, and K. Semba, Nature Phys. 13, F. Yoshihara, T. Fuse, S. Ashhab, K. Kakuyanagi, S. Saito, and K. Semba, Phys. Rev. A 95, N. Didier, J. Bourassa, and A. Blais, Phys. Rev. Lett. 115, Y. Yin, H. Wang, M. Mariantoni, R. C. Bialczak, R. Barends, Y. Chen, M. Lenander, E. Lucero, M. Neeley, A. O Connell, et al., Phys. Rev. A 85, S. Haroche and J.-M. Raimond, Exploring the Quantum: Atoms, Cavities, and Photons. Oxford, 13 Chap C. Gerry and P. Knight, Introductory Quantum Optics.

5 5 Cambridge University Press, 5 p M. O. Scully and M. S. Zubairy, Quantum Optics. Cambridge University Press, 1997 p S. Ashhab and F. Nori, Phys. Rev. A 81, J. Casanova, G. Romero, I. Lizuain, J. García-Ripoll, and E. Solano, Phys. Rev. Lett. 15, M. Hofheinz, H. Wang, M. Ansmann, R. C. Bialczak, E. Lucero, M. Neeley, A. O connell, D. Sank, J. Wenner, J. M. Martinis, et al., Nature 59, B. Vlastakis, G. Kirchmair, Z. Leghtas, S. E. Nigg, L. Frunzio, S. M. Girvin, M. Mirrahimi, M. H. Devoret, and R. J. Schoelkopf, Science 3, D. Braak, Q.-H. Chen, M. T. Batchelor, and E. Solano, Journal of Physics A: Mathematical and Theoretical 9, Y.-D. Wang, A. Kemp, and K. Semba, Phys. Rev. B 79, E. Schrödinger, Naturwissenschaften 1, V. M. Kodach, J. Kalkman, D. J. Faber, and T. G. van Leeuwen, Biomed. Opt. Express 1, L. Viola, E. Knill, and S. Lloyd, Phys. Rev. Lett. 8, G. de Lange, Z. H. Wang, D. Ristè, V. V. Dobrovitski, and R. Hanson, Science 33, E. L. Hahn, Phys. Rev. 8, SUPPLEMENTAL MATERIAL: FAST AMPLIFICATION AND REPHASING OF ENTANGLED CAT STATES IN A COUPLED-QUBIT OSCILLATOR Derivation of the state evolution with infinitesimal qubit frequency Under the limit of /, the Hamiltonian of the combined system of the qubit and the oscillator is Ĥ t= â â gtˆσ x â + â. 7 Let us examine the evolution of an initial state in the form of: Φ = ˆD g. 8 We can choose particular ± and N, but the following calculation applies to all ± and N. We can also see that with various choices of N and ±, Φ represents a complete orthonormal basis regardless of the value of g. By examining the evolution of initial state Φ with all choices of ± and N, we can understand the evolution of any general initial state, which itself can be expressed as a superposition of Φ with different ± and N. Also note that we can set g to any value. To calculate the evolution of the initial state Φ to the final time t F, we divide the time period of interest into K + 1 small segments, so that the entire period is divided by time points, t 1, t,, t i, t i+1,, t K, t F. Each segment is considered to be small enough so that, in a single segment gt and Ĥ t do not change much and are treated as constants. Therefore the final state at t=t F can be expressed as Φ t F = exp iĥ t K t F t K / exp iĥ t t +1 t / exp iĥ t 1 / ˆD g. 9 Now let us solve for Φ t F. At each time point t=t, starting with t=, we carry out the following procedures. 1 We make the ansatz that, at any time point t=t, the evolved state is in the form of Φ t =e iφ N t ˆD gt, 1 where ± and are the same as the initial state Φ. This is obviously true at t= and we will show, by mathematical induction, that indeed at each time point the state can be expressed in such a form. The complex number gt generally changes at different time points, and a phase e iφ t can accumulate as well. To make the expression more compact, let τ = t +1 t, g =gt, g = gt. The state at t +1, evolving from the state at the previous time point t, can be expressed as follows:

6 6 Φ t +1 =exp iĥ t t +1 t / Φ t { =exp i â â + 1 } ± g â + â τ =expiφ N t exp{ ii g g g / } exp ˆD g g ˆD g ± x =expiφ N t exp{ ii g g g / } exp { exp g g â ± g expiφ N t ˆD { i g â â + 1 ± g â + â τ } { i â â + 1 } ± g â + â τ g g g =expiφ N t exp{ ii g g g / } { exp g g â ± g e iτ g { exp i â â + 1 } ± g â + â τ =expiφ N t exp{ ii g g g / } ˆD exp g g g { exp i { =expiφ N t exp ii { exp i N + 1 g i â ± g g g ˆD g g g e iτ e iτ + g g g e +iτ + g g g N + 1 } g i τ ˆD g ± x } g g g 1 e iτ } τ ˆD g g g g } ˆD â ± g g g g e iτ, g e +iτ g } where the symbol I means an imaginary part. In the second to the last step of Eq. 11, we have used â â + 1 ± g â + â ˆD g = ˆD g â g â g ± g â g + â g + 1 = ˆD = ˆD g g = N + 1 â â + 1 g i N + 1 g i g i ˆD g, 11 1 and in the third to the last step we have used exp { iτ â â + 1 ± g â + â } â ± g = â ± g e iτ exp { iτ â â + 1 ± g â + â }. 13 Now recall that Φ t =expiφ N t ˆD gt, 1 Therefore we have proved, by mathematical induction,

7 7 that the quantum state at t=t +1 also has the form Φ t +1 =expiφ N t +1 ˆD gt +1, 15 with and gt +1 =gt + gt gt exp it +1 t 16 φ N t +1 =φ N t + I{ g g g 1 exp iτ } N + 1 τ + g i τ. 17 Since t +1 t =τ, we have gt +1 = gt 1 it +1 t gt it +1 t 18 and Therefore, which leads to φ N t +1 =φ N t + Ii g g g τ N + 1 τ + g i τ. 19 gt +1 gt t +1 t =igt gt, for which the solution is d gt=igt gt, 1 dt t gt= ge it + e it igt e it dt. As the state evolves under changing g, the phase φ N t F will accumulate and the final state is tf Φ t F =e iφ N t F ˆD gt F ± x. 3 where g = gt is given by Eq. and φ N t F = Ii g gg + g N + 1 dt. Notice in the phase φ N t F, the part of t F Ii g gg + g dt is the same for every N and ±, therefore we can simplify it as φ N t F = Nt F. Combining Φ N+ t F and Φ N t F, we arrive at Eqs. and 5 in the main text. Derivation of the state evolution with non-zero qubit frequency We now consider the case where is not infinitesimal and use the full Hamiltonian Ĥt=Ĥ t ˆ. As explained above, the dynamical energy eigenstates form a complete basis. Therefore at any time t, we can express any quantum state of interest as a superposition of dynamical energy eigenstates: ϕt = N,± C t Ẽ t. As the state ϕt evolves, so does every dynamical energy eigenstate Ẽ t and its corresponding coefficient C t. Now we need to solve the evolution of the coefficients C t to completely determine the evolution of state ϕt. Starting from initial time zero, we divide the time period of interest, t into K + 1 segments:, t 1, t,..., t K, t, with each segment being infinitesimal. For convenience we will adopt the matrix form, therefore the state at time t, which evolves from the initial state ϕ under Hamiltonian Ĥt, can be written as: ϕt = int C te Ẽ t N,± = C + t C 1+ t... C t C 1 t... Ẽ + t e it Ẽ 1+ T t... Ẽ t e it Ẽ 1 t... = exp iĥ t K t t K / exp iĥ t 1 t t 1 / exp iĥ t 1/ C + C C C 1... Ẽ + Ẽ 1+ T... Ẽ Ẽ 1....

8 8 From t to t +1, so long as t +1 t, the state evolves in the following way: { ϕt +1 = exp i Ĥ t } ˆ t +1 t / ϕt { = exp i Ĥ t } ˆσ int z t +1 t / C t e Ẽ t = C t exp i 1 ˆ t +1 t N,± N,± exp iĥ int t t +1 t / e Ẽ t = C t exp i 1 ˆ t +1 t e int+1 t int e Ẽ t +1 N,± = N,± N,± C t e in t +1 Ẽ N ± t +1 Ẽ N ± t +1 e +in t +1 exp i 1 ˆσ int+1 z t +1 t e Ẽ t +1 = C + t C 1+ t... C t C 1 t... 5 exp i 1 M t +1 t +1 t Ẽ + t it+1 +1 e Ẽ 1+ t Ẽ t +1 e where we have used the relation { exp i Ĥ t ˆ t +1 t / } it+1 Ẽ = exp i ˆ t +1 t / 1 t T, exp iĥ t t +1 t / expot +1 t t +1 t exp i ˆ t +1 t / exp iĥ t t +1 t /. The term M σz t in Eq. 5 is the matrix expansion of operator ˆ in the basis Ẽ t at time t: M σz t= M +,+ M +, t M,+ t M, t t 6, 7 where and M ±,± M ±, Ẽ ± t ˆ Ẽ ± t Ẽ 1± t ˆ Ẽ ± t e+it... t= Ẽ ± t ˆ Ẽ 1± t e it Ẽ 1± t ˆ Ẽ 1± t..., Ẽ ± t ˆ Ẽ t Ẽ 1± t ˆ Ẽ t e+it... t= Ẽ ± t ˆ Ẽ 1 t e it Ẽ 1± t ˆ Ẽ 1 t On the other hand, ϕt +1 = N,± int+1 C t +1 e Ẽ t +1 = C + t +1 C 1+ t C t +1 C 1 t Ẽ + t it+1 +1 e Ẽ 1+ t Ẽ t +1 e it+1 Ẽ 3 T 1 t

9 9 Comparing Eqs. 5, 3 and using mathematical induction, we have: C+ t C 1+ t... C t C 1 t... = C + C C C 1... J, t, 31 where J, t= exp 1 M σz t 1 t 1 i exp i 1 exp i M σz t t t 1 1 M σz tt t K. 3 Note that generally at different time points t and t, M σz t and M σz t do not commute, making the analytical calculation of the evolution matrix J, t very complicated. But in the case when we only consider up to the first-order effect of small, we can ignore the second order term, which means we consider the commutation M σz t, M σz t O, and therefore M σz t and M σz t approximately commute. This enables us to calculate the evolution matrix J, t and state ϕt up to the first order: and J 1, t=exp i 1 t tm σz tdt, 33 ϕ 1 t = C + C C C 1... { exp i 1 t Ẽ + t e it Ẽ 1+ t... Ẽ t e it Ẽ where for generality we can consider the parameter to be time dependent = t. } tm σz tdt 3 T 1 t...,