arxiv: v2 [cond-mat.mes-hall] 12 Nov 2012

Transcription

1 Tailoring population inversion in Landau-Zener-Stückelberg interferometry of flux qubits Alejandro Ferrón Instituto de Modelado e Innovación Tecnológica (CONICET-UNNE), 34 Corrientes and Centro Atómico Bariloche and Instituto Balseiro, 84 San Carlos de Bariloche, Río Negro, Argentina. arxiv: v2 [cond-mat.mes-hall] 2 Nov 22 Daniel Domínguez and María José Sánchez Centro Atómico Bariloche and Instituto Balseiro, 84 San Carlos de Bariloche, Río Negro, Argentina. We distinguish different mechanisms for population inversion in flux qubits driven by dc+ac magnetic fields. We show that for driving amplitudes such that there are Landau-Zener-Stückelberg interferences, it is possible to have population inversion solely mediated by the environmental bath. Furthermore, we find that the degree of population inversion can be controlled by tailoring a resonant frequency Ω p in the environmental bath. To observe these effects experiments should be performed for long driving times after full relaxation. PACS numbers: r,85.25.Cp,3.67.Lx,42.5.Hz Population inversion, where the most highly populated state is an excited state, is among the most interesting phenomena in maser and laser physics []. The usual mechanism to get population inversion (PI) requires driven quantum systems with three or more energy levels [ 4]. Interesting systems to study PI are artificial atoms made with mesoscopic Josephson devices [5, 6]. Among them, a well known circuit is the flux qubit (FQ) [7, 8]. When driven by a dc+ac magnetic flux, it has transitions between energy levels at avoided crossings. A rich structure of Landau-Zener-Stückelberg (LZS) interferences [9] combined with multi-photon resonances [ 2] is observed. LZS interference patterns have also been measured in charge qubits [3], Rydberg atoms [4], ultracold molecular gases [5], optical lattices [6] and single electron spins[7]. In FQ, diamond-like patterns were found for large ac amplitudes, from which the energy level spectrum has been reconstructed []. Population inversion was observed in the second diamond, i. e. for ac amplitudes that excite to the third and four energy levels through Landau-Zener transitions. PI was also measured in another flux qubit device driven at large frequencies [8]. Here we will show that an even more notable effect is awaiting to be observed in these systems if the ac driving pulses are applied for longer time scales: PI could be observable even for ac amplitudes where only the two lowest levels of the FQ participate. In the following we will explain, through realistic time-dependent numerical simulations, (i) how this type of PI can occur, (ii) why it has not been observed experimentally yet and (iii) how the occurrence and the degree of PI can be tailored by changing the structure of the environmental bath. The FQ consists of a superconducting ring with three Josephson junctions [7] enclosing a magnetic flux Φ = fφ (Φ = h/2e) with phase differences ϕ, ϕ 2 and ϕ 3 = ϕ + ϕ 2 2πf. Two of the junctions have coupling energy, E J, andcapacitance, C, and the otherhase J,3 = αe J and C 3 = αc. The hamiltonian is: [7] ( 2 H = 2E C + 2 ) +E J V(ϕ t,ϕ l ), () +2α ϕ 2 t ϕ 2 l with ϕ t = (ϕ +ϕ 2 )/2, ϕ l = (ϕ ϕ 2 )/2, E C = e 2 /2C and V(ϕ t,ϕ l ) = 2+α 2cosϕ t cosϕ l αcos(2πf+2ϕ l ). The FQ has several discrete levels with eigenenergies E i and eigenstates Ψ i which depend on f, α and η = 8E C /E J. Typical experiments have α.6.9 and η..6 [8,, ]. For α /2 and f /2, the potential V(ϕ t,ϕ l ) has the shape of a double-well with two minima along the ϕ l direction. In this regime the system can be operated as a quantum bit [7, 8] and approximated by a two-level system (TLS) [7, 9]. When driven by a magnetic flux f(t) = f dc + f ac sin(ωt), the hamiltonian is time periodic H(t) = H(t + τ), with τ = 2π/ω. In the Floquet formalism, which allows to treat periodic forces of arbitrary strength and frequency [2, 2], the solutions of the Schrödinger equation are of the form Ψ β (t) = e iε βt/ Φ β (t), where the Floquet states Φ β (t) satisfy Φ β (t) = Φ β (t+τ) = k Φk β e ikωt, andareeigenstates of the equation [H(t) i / t] Φ β (t) = ε β Φ β (t), with ε β the associated quasi-energy. Experimentally, the system is affected by the electromagnetic environment that introduces decoherence and relaxation processes. A standard theoretical model is to linearly couple the system to a bath of harmonic oscillators [22 25] with a spectral density J(ω) and equilibrated at temperature T. For the FQ, the dominating source of decoherence is flux noise, in which case the bath degrees of freedom couple with the system variable ϕ l (see [25]). For weak coupling (Born approximation) and fast bath relaxation (Markov approximation), a Floquet-Born-Markov master equation for the

2 2 reduced density matrix ˆρ in the Floquet basis, ρ αβ (t) = Φ α (t) ˆρ(t) Φ β (t), can be obtained [22]: dρ αβ (t) dt = i (ε α ε β )ρ αβ (t)+ α β L α β αβ(t) ρ α β (t). (2) The L α β αβ(t) are approximated by their average over one period τ [22] since the time scale t r for full relaxation satisfies t r τ, obtaining L α β αβ = R α β αβ + ( δ ββ R α αηη +δ αα R β βηη ). The rates Rβ α βα η R α β αβ = n gn αα Γn αα Γ n β β are sums of n-photon exchange terms, with gαα n = J(εn αα ( ) coth( εn αα 2k ) BT ), ε n αα = ε α ε α + n ω, and Γ n αα = k Φk α ϕ l Φ k+n α. This formalism has been extensively employed to study relaxation and decoherence for time dependent periodic evolutions in double-well potentials and in TLS [22 24]. Here we use it to model the ac driven FQ for the full hamiltonian of Eq. (), beyond the TLS approximation. We calculate the Floquet states Φ k β and quasienergies ε β [2], and then we integrate numerically Eq. (2), obtaining ρ αβ (t) as a function of t. P R A B C D t= 2π/ω Asymptotic P R.5 f =.245 ac 2 T [mk] f ac FIG. : Population P R as a function of the driving amplitude f ac. Continuous line: P R(t = τ t exp). Dashed line: asymptotic average population P R. Inset: P R vs. T for a driving amplitude f ac =.245. The calculations were performed for f dc =.55, ω =.3E J/ and an ohmic bath with γ =. at T = 2mK (for E J/h 3GHz). Vertical lines separate regimes A, B, C, D described in the text. Black circle: f ac corresponding to the inset, Fig.2(a) and Fig.3. Black square: f ac corresponding to Fig.2(b). WestartbyconsideringaFQcoupledtoabathwithan ohmic spectral density J(ω) = γω. We take γ =., corresponding to weak dissipation as in [], and the bath at T =.4E J /k B ( 2mK for E J /h 3GHz). TheFQdeviceparametersareα =.8andη =.25. The static field is taken as f dc =.55, and for the driving microwave field we choose ω =.3E J / 9MHz and different amplitudes f ac. The initial condition is the ground state Ψ of the static hamiltonian H H(f dc ). Experimentally, the probability of having a state of positive or negative persistent current in the FQ is measured [8]. The probability of a positive current measurement ( right side of the double-well potential) can be calculated as P R (t) = Tr(ˆΠ Rˆρ(t)), with ˆΠ R the operator that projects wave functions on the ϕ l > subspace [2]. For a static field f dc /2, the ground state has P R. In Fig. we show, as a function of f ac, P R (t) after a driving time t = τ, which is similar to the time scale t exp of the experiments []. We also plot the asymptotic P R lim t P R (t) τ, averaged over one period τ. We find different regimes: (A) For small f ac, P R, since the system is slightly perturbed from the ground state. (B) When further increasing f ac, new regimes are found whenever the extreme driving amplitudes f dc ±f ac reach an avoided crossing in the energy level spectrum {E i (f))}. Since the slopes dei df are proportional to the average loop current of the FQ, at avoided crossings there is interference between positive and negative loop current states, which results in LZS oscillations [9] in the dependence of P R with f ac. The first case is found when the avoidedcrossingat f = /2between the groundstate level E and the first excited level E is reached, and the transfer of population to the E level lowers P R. The FQ of [, ] has a decoherence time t φ 2ns > τ and a large interwell relaxation time t r µs. Due to this time scaleseparation(t φ < t exp t r ), a TLS model with classical noise, valid for t < t r, can explain the experimentally observed LZS oscillations of P R vs. f ac, where the minima correlate with the zeros of Bessel functions J n (af ac /ω) (a a normalization constant) [9, ]. The results of Fig. for t = τ t exp are in agreement with this finding. However, the finite time P R (t exp ) is very different from the asymptotic P R, that even shows population inversion, P R < /2. (C) At higher values of f ac the avoided level crossing between E and E 2 is reached and the ground level E is repopulated due to fast E 2 E intrawell transitions (the corresponding states have the same sign of the average loop current). Thus, P R increases close to ( cooling effect, see []). When more than two levels are involved, intrawell relaxation in a time scale t i is possible. In the experiments, t i 5ns t r []. This explains that P R (t exp ) P R, since intrawell processes dominate relaxation to the asymptotic state and t i < t exp. (D) For f ac such that the symmetrically located (with respect to f = /2) avoided crossing between E and E 2 is reached, there are new LZS oscillations. Furthermore the levels E and E 3 are also involved in the dynamics (since the avoided crossings between E 2 E 3 and E E are traversed by the driving) and a more complex dependence of P R with f ac emerges. In this case, we find population inversion, which is also observed experimentally in []. A multilevel extension of the semiclassical model of [] can describe this behavior as well [26]. InFig.2weshowtheexplicittimedependenceofρ αβ (t)

3 3 ρ ab ρ ab..5 ρ ab 2e-7 γ (n) e-7 ρ ab f.25.3 ac P - P 3 ωt/2π P P 2 P P (a) (b).5 P i.5 P i FIG. 2: (color online) Off-diagonal (in Floquet basis) density matrix element ρ ab (left axis) and level population P i (right axis) as a function of time (log) for f dc =.55, ω =.3E J/, γ =.. (a) f ac =.245 and (b) f ac =.9. Inset: n-photon relaxation terms γ (n) (the corresponding n value is indicated on the top of each curve) in the neighborhood of f ac =.245 (indicated by the arrow). for two values of f ac in regimes (B) and (D) respectively. The off-diagonal elements ρ αβ decay exponentially to zero in a time scale t φ [see Fig.2(a) and (b), left axis]. In the Floquet basis, after full decoherence, the relaxation process is determined by the evolution of the diagonal ρ αα (t) for t > t φ [27]. Thus the asymptotic regime can be obtained from the non-trivial solution of β L ββ,ααρ ββ = after imposing dραα dt = in Eq.(2). In this way, we calculated the P R shown in Fig. as a function of f ac. In orderto understand the emergence of population inversion in the asymptotic limit, we evaluate the population of the eigenstates of H, P i (t) = Ψ i ˆρ(t) Ψ i. We distinguish two mechanisms of PI: (i) Third-level-mediated population inversion: For large f ac [regime (D) in Fig.], the populations P i relax to the asymptotic regime in a short time scale similar to t i t φ. In Fig.2(b) we see that the populations P 2,P 3 start increasing, and later they decrease transferring population to P. This leads to population inversion, P > P, mediated by E 2 and E 3, which is the usual mechanism for this effect. In our case the higher levels are reached through LZS transitions[, 26] instead of resonant transitions. (ii) Bath-mediated population inversion: In Fig.2(a) we show results for a value of f ac =.245 such that only LZS interference between the two lowest energy levels occurs [regime (B) in Fig.]. Thus one can restrict to two Floquet states α = a, b. We find that the decay time t φ of the off-diagonal ρ ab is similar to the previous case but the level populations relax to their asymptotic value in a very large time scale t r t φ. Only for t t r population inversion can arise. The large t r explains the difference between P R and the finite time P R (t t exp ) in Fig.. To understand the underlying mechanism we decompose the relaxation rate γ r = t r as γ r = 2(R aabb + R bbaa ) = γ () + n γ(n), with γ (n) = 2(g n ab Γn ab 2 + g n ba Γn ba 2 ) the terms that describe virtual n-photon transitions to bath oscillator states [24]. In the inset of Fig.2(a) we show that when there is PI the term γ (n= ) is the largest while γ (n=). This indicates that the dominant mechanisms leading to PI is a transition to a virtual level at energy E + ω > E (one photon absorption, n = ), followed by a relaxation to E. Previous works have found PI, under various approximations, but for the asymptotic regime of TLS [28 3]. However the time-dependent dynamics with different time scales has been overlooked. In fact, the relevant point from our findings is that the asymptotic regime is difficult to reach in the experiment, since PI needs long times (t t r t φ ) to emerge when mediated by the bath. Moreover, the difference between P R (t exp ) and the asymptotic P R, is enhanced when decreasing temperature (see inset of Fig.) P R Stationary t= τ Ω p FIG. 3: Population P R as a function of the frequency Ω p (in units of E J/ ) for T = 2mK, f dc =.55, f ac =.245, ω =.3E J/, γ =. and κ =.. Continuous line shows the asymptotic P R and dashed line corresponds to P R(t = τ t exp). The circles indicate the values of + and. The inset shows P R as a function of Ω p around the n = 9 resonance. See text for details. An interesting consequence of the bath-mediated mechanism is that it enables - by changing the spectral structure of the bath - to modify the asymptotic state of the FQ. In particular, a realistic modeling of the environmental bath for FQ includes a read-out dc SQUID inductively coupled to the FQ [7, 25, 32]. For this case, the bath spectral density becomes J sb (ω) = γωω 4 p/[(ω 2 p ω 2 ) 2 + (2πκωΩ p ) 2 ], with Ω p the resonant frequency of the SQUID detector, and 2πκΩ p the width of the resonance at ω = Ω p [25]. In Fig.3 we show P R as a function of Ω p for the same f ac of Fig.2(a), but after solving Eq.(2) with J sb (ω). We identify three regimes in

4 4 the overall behavior of P R vs Ω p that we describe below. Defining (f) = E (f) E (f), + = (f dc +f ac ) and = (f dc f ac ) we get: () Multiphoton resonances with the bath. For Ω p < the population P R has a strong dependence on Ω p, displaying resonances for Ω p nω + ε b ε a (with ε b,ε a the quasienergies of the Floquet states mostly weighted in the two lowest H eigenstates). These resonances correspond to the maxima at Ω p of the J sb (ε a ε b +n ω), and have been previously discussed for n = within a perturbative approach [23]. Near these resonances it is possible to fully control the population P R, going from P R to complete population inversion P R, with small changes of Ω p. In the inset of Fig.3 we show in detail the case near the n = 9 resonance. (2) Bath mediated population inversion. In the opposite limit, Ω p > +, the behavior is similar to the one discussed previously in Fig.2(a), since for large Ω p, the spectral density J sb (ω) γω (3) Bath mediated cooling. In the intermediate regime, Ω p +, transitions to an effective energy level [23, 32] at E + Ω p < E predominate. From this effective level it is possible to decay and repopulate the ground state, obtaining P R. The resonances at Ω p nω + ε b ε a are still observed. A full picture of the effect of a structured bath can be observed in Fig.4, which shows a map of the asymptotic population P R as a function of Ω p and f ac. The onset of third level mediated PI at higher f ac (corresponding to the second diamond of []) is also observed. Notice that this type of PI is almost independent of Ω p, while the bath-mediated mechanism is active for Ω p > +. In conclusion, by performing a realistic modeling of the FQ we are able to analyse different scenarios for population inversion in strongly driven quantum systems, understanding the parameter regimes for their occurrence. One puzzling situation is that the bath-mediated PI discussed here has not been observed in []. The large time scale separation t r t φ in the device of [] suggests that the PI could be beyond the experimental time window [33]. Indeed, our results explain that to observe this effect the experiments should be carried out for long driving times, after full relaxation with the bath degrees of freedom (t t r ). A remarkable point we found is the dependence of the LZS oscillations on the frequency Ω p ofthe SQUIDdetector. The FQfabricatedwith Nb junctions as in [, ] have typically a small gap and thus they are in the regime Ω p > +, where bath-mediated PI can take place. The FQ fabricated with Al junctions, as in [8, 32], have a large gap and thus they are expected to be in the regime Ω p <. Amplitude spectroscopy measurements carried out in these later systems could give different diamond patterns as a function of f dc,f ac with resonances and LZS interference effects from the oscillator levels at E i +n Ω p. In principle, the frequency Ω p can be controlled by varying in the SQUID detector the driving current or the shunt capacitor [32], or by adding a tunable resonator to the circuit [34]. In these cases there is room to fully tailor the population P R. Furthermore, the discussed population inversion scenarios could also apply to other quantum systems in which LZS interference has been measured [3 7]. We acknowledge discussions with Sergio Valenzuela and financial support from CNEA, CONICET (PIP22882 and PIP2295) and AN- PCyT (PICT and PICT27-824). FIG. 4: (color online) Contour map of the asymptotic population P R as a function of the driving amplitude f ac and the resonance frequency Ω p (in units of E J/ ) for T = 2mK, f dc =.55, ω =.3E J/, γ =. and κ =.. The black and green lines represent the positions of + and. Black circle: f ac corresponding to Fig.2(a) and Fig.3; black square: f ac corresponding to Fig.2(b). [] J. Weber, Rev. Mod. Phys. 3, 68 (959) [2] S. Y. Han, R. Rouse, and J. E. Lukens, Phys. Rev. Lett. 76, 344 (996). [3] J. Q. You et al, Phys. Rev. B 75, 456 (27). [4] O. Astafiev et al, Nature 449, 588 (27). [5] Y. Makhlin, G. Schön, and A. Shnirman, Rev. Mod. Phys. 73, 357 (2). [6] J. Q. You and F. Nori, Nature 474, 589 (2). [7] T.P.Orlando et al, Phys. Rev. B 6, 5398 (999). [8] I. Chiorescu et al, Science 299, 869 (23). [9] S.N. Shevchenko, S. Ashhab and Franco Nori, Phys. Rep. 492, (2). [] W. D.Oliver et al, Science 3, 653 (25); D. M. Berns et al, Phys. Rev. Lett. 97, 552 (26). [] D. M. Berns et al, Nature 455, 5 (28); W D. Oliver and S. O. Valenzuela, Quantum Inf. Process 8, 26 (29). [2] A. Izmalkov et al, Phys. Rev. Lett., 73 (28). [3] Y. Nakamura, Yu. A. Pashkin, and J. S. Tsai, Phys. Rev. Lett. 87, 2466 (2); M. Sillanpaa et al, Phys. Rev. Lett. 96, 872 (26). C. M. Wilson et al, Phys. Rev. Lett. 98, 2573 (27). [4] S. Yoakum et al, Phys. Rev. Lett. 69, 99 (992).

5 5 [5] M. Mark et al, Phys. Rev. Lett. 99, 32 (27). [6] S. Kling et al, Phys. Rev. Lett. 5, 253 (2) [7] P.Huang et al, Phys. Rev. X, 3 (2). [8] G. Z. Sun et al, Appl. Phys. Lett. 94, 252 (29); [9] A. Ferrón and D. Dominguez, Phys. Rev. B 8, 455 (2). [2] J. H. Shirley, Phys. Rev. 38, B979 (965). [2] A. Ferrón, D. Domínguez and M. J. Sánchez, Phys. Rev. B. 82,34522 (2). [22] S. Kohler et al, Phys. Rev. E. 58, 729 (998), H.-P. Breuer, W. Huber and F. Petruccione, Phys. Rev. E 6, 4883 (2); D. W. Hone, R. Ketzmerick and W. Kohn, Phys. Rev. E 79, 529 (29). [23] M. C. Goorden, M. Thorwart and M. Grifoni, Phys. Rev. Lett. 93, 2675 (24); Eur. Phys. J. B 45, 45 (25). [24] J. Hausinger and M. Grifoni, Phys. Rev. A 8, 227 (2). [25] C. H. van der Wal et al, Eur. Phys. J. B 3, (23); G. Burkard, R. H. Koch, and D. P. DiVincenzo, Phys. Rev. B 69, 6453 (24). [26] X. Wen and Y. Yu, Phys. Rev. B. 79, (29). [27] This is expected for ε α ε β L α β αβ, which is fulfilled for very weak coupling with the environment, see [22]. [28] Yu. Dakhnovskii and R.D. Coalson, J. Chem. Phys. 3, 298 (995); I.A. Goychuk, E.G. Petrov, and V. May, Chem. Phys. Lett. 253, 428 (996). [29] L. Hartmann et al, Phys. Rev. E 6, R4687 (2). [3] T.M. Stace, A. C. Doherty, and S. D. Barrett, Phys. Rev. Lett. 95, 68 (25). [3] L. Du, M. Wang, Y.Yu., Phys.Rev. B82, 4528 (2). [32] I. Chiorescu et al., Nature (London) 43, 59 (24); J. Johansson et al, Phys. Rev. Lett. 96, 276 (26) [33] Moreover, t r/t φ in [] is much larger than for the parameters used here, where for numerical convenience we have chosen a larger gap and a larger γ. [34] F. Deppe et al, Nature Physics 4, 686(28); A. Fedorov et al, Phys. Rev. Lett. 5, 653 (2).