Distinguishing quantum from classical oscillations in a driven phase qubit

Rabi oscillations are coherent transitions in a quantum two-level system under the influence of a resonant perturbation, with a much lower frequency dependent on the perturbation amplitude. These serve as one of the signatures of quantum coherent evolution in mesoscopic systems. It was shown recently [N. Gronbech-Jensen and M. Cirillo, Phys. Rev. Lett. 95, 067001 (2005)] that in phase qubits (current-biased Josephson junctions) this effect can be mimicked by classical oscillations arising due to the anharmonicity of the effective potential. Nevertheless, we find qualitative differences between the classical and quantum effect. First, while the quantum Rabi oscillations can be produced by the subharmonics of the resonant frequency (multiphoton processes), the classical effect also exists when the system is excited at the overtones. Second, the shape of the resonance is, in the classical case, characteristically asymmetric; while quantum resonances are described by symmetric Lorentzians. Third, the anharmonicity of the potential results in the negative shift of the resonant frequency in the classical case, in contrast to the positive Bloch-Siegert shift in the quantum case. We show that in the relevant range of parameters these features allow to confidently distinguish the bona fide Rabi oscillations from their classical Doppelganger.


I. INTRODUCTION
Superconducting phase qubits [1,2] provide a clear demonstration of quantum coherent behaviour in macroscopic systems. They also have a very simple design: a phase qubit is a current-biased Josephson junction (see Fig. 1(a)), and its working states |0 , |1 are the two lowest metastable energy levels E 0,1 in a local minimum of the washboard potential. The transitions between these levels are produced by applying an RF signal at a resonant frequency ω 10 = (E 1 − E 0 )/ ≡ ∆E/ . The readout utilizes the fact that the decay of a metastable state of the system produces an observable reaction: a voltage spike in the junction or a flux change in a coupled dc SQUID. In the three-level readout scheme ( Fig. 1(b)) both |0 and |1 have negligible decay rates. A pulse at a frequency ω 21 = (E 2 − E 1 )/ transfers the probability amplitude from the state |1 to the fast-decaying state |2 . Its decay corresponds to a singleshot measurement of the qubit in state |1 . Alternatively, instead of an RF readout pulse one can apply a dc pulse, which increases the decay rate of |1 . One of the effects observed in driven phase qubits is Rabi oscillations [2,3]: coherent transitions in a quantum two-level system under the influence of a resonant perturbation, with a much lower frequency dependence on the perturbation amplitude A via Ω R = √ A 2 + δ 2 , where δ = ω − ω 10 is the detuning of the driving frequency ω from the resonance frequency ω 10 . In resonance, Ω R = A. Multiphoton Rabi oscillations, at ω 10 = nω (n stands for an integer), observed in such qubits [4,5], were also interpreted as a signature of coherence.
The quantum coherent dynamics in phase qubits was tested by several complementary methods (see, e.g., Refs. [3,5,6,7]). In particular, in Ref. [6] Rabi oscillations between the ground and excited state of the phase qubit were measured by applying a 25 ns pulse at ω 10 followed by a measurement pulse at ω 21 . The probability to find the system in the upper state oscillated with the amplitude of the resonant signal, as it should in case of Rabi oscillations. In Ref. [7] Rabi oscillations were seen instead as the resonant pulse duration varied; a dc readout was used.
However, it is known that the quantum mechanical behaviour of a quantum two-level system can still be similar to the dynamical behaviour of the classical non-linear oscillator [8,9,10,11,12,13]. In particular, it was recently pointed out [14,15] that due to the nonlinear behaviour of current-biased Josephson junctions, a similar effect can arise in a purely classical way. Though direct tests -à la Bell -of whether a given system is quantum or not are possible, even their simplified versions (e.g., Ref. [13]) are demanding. This motivated us to further investigate the classical behaviour of a phase qubit and find that there is a possibility to distinguish the quantum Rabi oscillations from its classical double, by the shape of the resonance, by the fact that the classical effect can be also produced by the overtones, nω 10 , of the resonance frequency, and by the sign of the resonant frequency shift. (An observation that a non-Lorentzian shape of the resonance should exclude a classical explanation was made already in Ref. [16], where the spectroscopy of two coupled qubits was performed. Also, a symmetric versus asymmetric Stark shift in a qubit playing the role of a detector was proposed to distinguish the classical and quantum behaviour of a nanomechanical oscillator [11].) Classical and quantum resonances, as a function of applied drives are also studied in Ref. [17].

II. MODEL
The phase qubit [1,2] is a current-biased Josephson junction. The Josephson potential, as a function of the phase difference φ, forms local minima, in which the quantized metastable levels serve as the working states of the qubit. Here I c is the critical current of the junction, and I dc is the static bias current. A perturbation can be produced by applying a time-dependent bias current, I ac sin ωt. In the quantum case, the system can be reduced to a two-level model, described by the Hamiltonian where σ z,x are Pauli matrices, C is the capacitance of the junction, and χ ≈ 1/4 in the relevant range of parameters.
One can obtain from here the Rabi oscillations (coherent oscillations of the probability to find the system in the upper/lower state with the frequency Ω R ) when excited near the resonance or at its subharmonics; the shape of the resonance is a symmetric Lorentzian, as determined by, e.g., the average energy of the system versus the driving frequency ω.
Unlike the flux qubit [18,19], where the interlevel distance ∆E is determined by the tunneling, here ∆E is close to the "plasma" frequency ω p = [2eI c / C] 1/2 of small oscillations near the local minima of the potential, Eq. (1): characteristically ∆E/ ≈ 0.95ω p [2]. The same frequency determines the resonance in the system in the classical regime. This means that, in principle, the same AC signal could cause either the Rabi oscillations or their classical double [14].

III. CLASSICAL REGIME
In the classical regime, the phase qubit can be described by the RSJJ (resistively shunted Josephson junction) model [20,21], in which the equation of motion for the superconducting phase difference across the junction, characterized by the normal (quasiparticle) resistance R, reads C 2e Introducing the dimensionless variables, we obtainφ where and the dot stands for the derivative with respect to τ . The solution is sought in the phase-locked Ansatz, We substitute Eq. (7) in Eq. (5) and expand sin φ to third order, which yields Therefore sin φ 0 = η, and introducing we obtainψ which describes the anharmonic driven oscillator [22].
Here we briefly point out several features of the solution of Eq. (10) (for more details see Chap. 5 in Ref. [22]).
(ii) The amplitude b of the small driven oscillations at the main resonance is b max = ǫ/(γ 0 α). When this amplitude is not small, then the phase-locked Ansatz becomes invalid. Then the solution of Eq. (5) describes the escape from the phase-locked state, which means the appearance of a non-zero average voltage on the contact. This voltage is proportional to the average derivative of the phase,φ.
(iii) The position of the resonances for small oscillations is shifted due to the anharmonicity of the potential. This, for the main resonance γ res , is given by: Note that the resonance shift is negative (e.g. as in [11]).
(iv) The shape of the resonances, as a function of the driven frequency γ, is essentially non-symmetrical. The asymmetry of the main resonance becomes pronounced at at small enough damping.
(v) The parametric-type resonance at γ ≈ 2γ 0 appears when the damping is sufficiently low, i.e. at α < ηǫ At the relevant parameters, η ≃ 0.9, γ 0 ≃ 0.6, this means the following: which is fulfilled (see (ii)) when the solution close to the main resonance corresponds to the escape from the phase-locked state.
We note that for the anharmonic driven oscillator, described by Eq. (10), both the resonances at γ 0 /2 and 2γ 0 appear due to the anharmonicity of the potential energy and are of the same order. For the driven flux qubit in the classical regime [19] the equation for the phase variable θ(t) = θ 0 + ψ(t) can also be expanded for small oscillations ψ(t) about the value θ 0 ; restricting ourselves here to the linear in ψ terms, the equation can be rewritten in the form: In this case the genuine parametric resonance at γ ≈ 2γ 0 appears due to the term containing sin γτ · ψ [22]. This explains the prevailing of this resonance over the resonance at γ ≈ γ 0 /2, due to the small anharmonicity of the potential in Ref. [19]. Now we proceed to numerically solve the equation of motion (5) for the relevant set of parameters close to the experimental case. We also investigate the behaviour of the energy of the system [20,21] which determines the thermally activated escape probability from the local minimum of the potential [23] in Eq. (1). The classical Rabi-like oscillations are displayed in Fig. 2. In Fig. 2(a) the modulated transient oscillations of the phase difference φ are plotted. These oscillations result in the oscillating behaviour of the energy of the system as shown in Fig. 2(b). Averaging over fast oscillations, we plot in Fig. 2(c) with green solid curve the damped oscillations of the energy, analogous to the quantum Rabi oscillations [14]. These curves are plotted for the following set of the parameters: η = 0.95, α = 10 −3 , ǫ = 10 −3 , and γ = γ 0 . For comparison we also plotted the energy averaged over the fast oscillations for different parameters, changing one of these parameters and leaving the others the same. The dashed black curve in Fig. 2(c) is for the smaller damping, α = 10 −4 ; the solid (violet) line and the dash-dotted (black) line in Fig. 2(d) demonstrate the change in the frequency and the amplitude of the oscillations respectively for ǫ = 2 · 10 −3 and η = 0.9. We notice that the effect analogous to the classical Rabi oscillations exist in a wide range of parameters. In Fig. 3 the effect of the driving current on the time-averaged energy of the system is shown for different driving amplitudes: Fig. 3(a) for weaker amplitudes, close to the main resonance, to show the asymmetry and negative shift of the resonance; and Fig. 3(b) for stronger amplitudes, to show the resonances at γ 0 /2 and 2γ 0 (which are also shown closer in the insets). We note that the parametric-type resonance at 2γ 0 originates from the third-order terms when the solution of the equation for ψ is sought by iterations [22]; when there are two or more terms responsible for this resonance, the respective resonance may become splitted, which is visible in Fig. 3(b) for the lowest curve. An analogous tiny splitting of the resonance was obtained for the driven flux qubit in Fig. 4 of Ref. [19].

IV. QUANTUM REGIME
In the quantum regime, the phase qubit can be described by the Bloch equations for the density matrix components. In order to take into account the relaxation and dephasing processes, the corresponding rates Γ relax ≡ ω p · λ relax and Γ φ ≡ ω p · λ φ are included in the Liouville equation phenomenologically [24]. Then the evolution of the reduced density matrix ρ, taken in the form is described by the Bloch equations [24,25]: where A and C stand for the off-diagonal and diagonal parts of the dimensionless Hamiltonian: From these equations we obtain Z(τ ), which defines the occupation probability of the upper level, We choose the initial condition to be X(0) = Y (0) = 0, Z(0) = 1, which corresponds to the system being in the ground state; we also consider the zero-temperature limit in which the equilibrium value of Z is Z 0 = 1. When the system is driven close to resonance, ω ≈ ∆E, the upper level occupation probability P (τ ) exhibits Rabi oscillations. The damped Rabi oscillations are demonstrated in the upper inset in Fig. 4, which is analogous to the classical oscillations presented in Fig. 2. After averaging the time dependent probability, we plot it versus frequency in Fig. 4 for two values of the amplitude, demonstrating the multiphoton resonances. Figure 4 demonstrates the following features of the multiphoton resonances in the quantum case: (a) in contrast to the classical case, the resonances appear only at the subharmonics, at ω ≈ ∆E/n; (b) the resonances have Lorentzian shapes (as opposed to the classical asymmetric resonances); (c) with increasing the driving amplitude the resonances shift to the higher frequencies -the Bloch-Siegert shift, which has the opposite sign from its classical counterpart. The Bloch-Siegert shift (the shift of the principal resonance at ω ≈ ∆E) is plotted numerically in the right inset in Fig. 4. Analogous shifts of the positions of the resonances were recently observed experimentally [5].

V. CONCLUSIONS
In conclusion, the following criteria can be proposed to distinguish classical Rabi-type oscillations in current-biased Josephson junctions: (1) the appearance of resonances both at the overtones of the main resonant frequency (e.g., 2γ 0 ) and its subharmonic (γ 0 /2); (2) the asymmetric shape of the resonances; (3) a negative shift of the resonant frequency when increasing the driving amplitude. In recent publications these features were not reported: the multiphoton resonances were observed [4,5]; the resonances observed have a Lorentzian shape [3,16,26] and have a positive frequency shift [5], which means that the observed resonant excitations in the system were in the quantum regime.