arxiv: v2 [cond-mat.mes-hall] 23 Jun 2011

Transcription

1 Charge Qubit Purification by an Electronic Feedback Loop Gerold Kießlich, Gernot Schaller, Clive Emary, and Tobias Brandes Institut für Theoretische Physik, Hardenbergstraße 36, Technische Universität Berlin, D-63 Berlin, Germany arxiv:.377v [cond-mat.mes-hall] 3 Jun We propose the manipulation of an isolated qubit by a simple instantaneous closed-loop feedback scheme in which a time-dependent electronic detector current is directly back-coupled into qubit parameters. As specific detector model we employ a capacitively coupled single-electron transistor. We demonstrate the stabilization of pure delocalized qubit states above a critical detector-qubit coupling. This electronic purification is independent of the initial qubit state and is accomplished after few electron jumps through the detector. Our simple scheme can be used for the efficient and robust initialization of solid-state qubits in quantum computational algorithms at arbitrary temperatures. PACS numbers: 3.65.Ta, 5.4.-a, b, p The ability to measure and control individual quantum systems (e.g. qubits) is a major requirement for contemporary and future information technologies[]. The measurement process the contact between classical measurement device and nano-scale object plays a key role in controlling quantum systems and mostly introduces a source of decoherence. State-of-the-art measurement technologies in solid-state environments mainly rely on the change of current flow through an on-chip detector due to Coulomb repulsion with the charges in the measured system (see e.g. Refs. []). The continuous monitoring of charge qubits by such charge detectors yields relaxation and decoherence as studied experimentally e.g. in Ref. [3] and theoretically in Refs. [4]. However, the detector signal contains detailed information on the qubit dynamics which can be used for its control. An elegant control strategy, which has been successfully applied in a wide variety of physical and engineering systems [5], is stabilization of favored system states due to feedback. For quantum systems such as qubits, several theoretical proposals to enhance quantum coherence by feedback have been made. These are mainly based upon sophisticated techniques, e.g., the on-chip conversion of detector outputs into qubit density matrix evolutions [6, 7] or monitoring the phase of qubit oscillations by detector current quadratures [8]. However, a much simpler scheme was recently used in the first experimental use of feedback to increase quantum coherence of a spin qubit [9]. Here, we study a charge qubit consisting of a single electron confined in a two-state system which is continuously monitored by the current through a single-electron transistor (SET) [see Fig. ]. The SET current will be directly fed back into qubit parameters which allow the stabilization of arbitrary pure qubit states. Model. The SET electronic level is shifted depending on where the electron is located in the two-level system due to Coulomb interaction. Moreover, the charge qubit electron is assumed to be coupled to thermal phonons of the surrounding matrix material. The full Hamiltonian reads H = H SET +H q + d d [ Uˆσ z +(U t +U b ) ], () H q = εˆσ z +T cˆσ x + [ ] ˆσ zg Q A Q +Ω Q a Q a Q, Q with ε ε t ε b and A Q a Q + a Q. The SET with Hamiltonian H SET = ε d d d+ ka=l,r [ε kac ka c ka+ (t ka c ka d + h.c.)] contains a single electronic level (with creation/annihilation operators d /d) which is electrostatically coupled to the charge qubit levels by U U t U b. The charge qubit Hamiltonian is given in the Bloch representation with Pauli matrices ˆσ z = t t b b and ˆσ x = t b + b t. Thephononic modes with wave vector Q and energy Ω Q are diagonally coupled to the qubit assuming weak coupling T c. Equations of Motion. The standard Born-Markovsecular approximation yields a quantum master equation ρ(t) = Lρ(t) for the system (SET + qubit) density matrix ρ(t) with the Liouvillian superoperator L where the SET leads and the phonon environment are treated in thermal equilibrium. Due to the secular approximation the 4 populations and coherences contained in ρ(t) decouple in the energy eigenbasis. In the Markovian limit, tunneling rates and bath expectation values are all evaluated at the transition frequencies of the system. Since the SET Coulomb shift depends on the position of the charge qubit electron, this introduces for finite SET bias effective tunneling rates Γ a and Γ a for lower and upper electron position, respectively, such that the (time and ensemble averaged) SET current comprises damped oscillations for Γ a Γ a, see Fig. (d). These are caused by the qubit evolution [Fig. (e)]. With the SET occupation n(t) = Tr[ˆnρ(t)] and the Bloch vector components conditioned on whether the SET is occupied, S α(t) Tr[ (ˆn ˆσ α )ρ(t) ] and S α(t) Tr{ [( ˆn) ˆσ α ]ρ(t) } (α = x,y,z) respectively, such that S α (t) = S α(t) + Sα (t), the quantum master equation has the following convenient vector form (omit-

2 QPC SET L +V/ (a) phonons Γ L Γ R ε d +U Γ L ε d ε t U ε b T C qubit Γ R R -V/ IL IR Current I S z n IQPC (b) (c) (d) (e) t FIG. : (color online). (a) Sketch of the device: charge qubit and single-electron transistor (SET) as capacitively coupled detector; a quantum point contact (QPC) is attached to the SET to resolve I L(t) and I R(t) for the closed-loop feedback. (b) Single realization of the SET occupation n(t) measured by the QPC current I QPC(t). (c) Single realizations of SET currents I L(t) (blue, solid) and I R(t) (dashed, green). (d) SET current Ī (t), ensemble averaged. and time averaged Ī = t+τ I(t )dt with I(t) = [IL(t) + IR(t)], exhibits τ t damped oscillations caused by the coupling to the qubit. (e) Bloch vector component S z (t) (qubit occupation): ensemble average (red, dashed) shows oscillations damped by the presence of the SET, single realization (blue, solid) exhibits free evolution interrupted whenever an electron jumps at the SET. Parameters: Γ /Γ = 4, T c/γ =, ε =, U/Γ =.5, measurement time Γτ =., initial state: S z () =, n() =. ting the time dependence): ṅ = γ + L γ+ n γ L Sz +γ R Sz, Ṡ = A S +D L S +B R S +n n+(η r,,γ L )T, Ṡ = A S +D R S +B L S +n n (,,γ L )T,() with S i (t) (Si x(t),sy i (t),sz i (t))t, γ a ± (Γ a ±Γ a ), γ ± γ ± L +γ± R, n = n (η r,,γ ) T, and D a γ+ a η p η + +η γ a + η p, γ a + B a Diag( Γa Γ a, ) Γ a Γ a,γ a +, (3) with the phonon rates η p and η r η + η explicitly defined in Ref. []. For zero qubit bias (ε = ) they are η r = gπt c e Tc/ωc and η p = η r coth(βt c ) with the dimensionless coupling strength g, the cut-off frequency ω c, and the inverse temperature β. The free dynamics of S j (t) in Eqs. () can be considered as a precession around vector A j ( T c,,ε+ju ) T with frequency ε +4Tc (j = ) or (ε+u) +4Tc (j = ), respectively. The effect of dissipation is contained in D a : damping of all Bloch vector components by the presence of the SET with γ + L or γ+ R and by the electronphonon coupling which leads to an additional decay and Γt 4 3 (a) (b) (c) Θ Rε Θ RΤ S z (t) FIG. : (color online). Bloch vector component S z (t) vs. time Γt and feedback parameter Θ Rε (a) and Θ RT (b); (c) S z (t) cross-sections at Θ Rε = (blue, dash-dotted), Θ RT =.5 (red, full),.5 (green, dashed). Parameters: ε =, U/Γ =.5, Γ /Γ = 3, T c/γ =. Initial state is the complete statistical mixture with empty SET. With feedback damped qubit oscillations can be induced which would not occur for vanishing SET-qubit coupling Γ = Γ. scattering of the S x/y i (t). The conditioned Bloch vectors S j (t) are mutually coupled via the diagonal matrices B L/R, i.e., by electron hopping onto the SET and out of it, respectively. We provide an analytical solution in the online supplement []. The SET occupation couples to the qubit and is not affected by its dynamics for Γ a = Γ a. We note that the Coulomb interaction U enters Eqs. () implicitly in Γ and explicitly in A. Whenever we set U = in the following we do not exclude Γ Γ in order to maintain the coupling but rather neglect the physical back-action on the qubit. Detector counting statistics and active measurement. The information about the number of charges n a transferred through a multi-lead conductor (our SET detector) at the a-th lead can be gained by the help of n a - resolved quantum master equations (see e.g. Ref. [] for detailed discussions). The unraveled Fourier-transformed master equation for unidirectional transport is ρ = [L + a J eiχa a ]ρ L({χ a })ρ, where the J a (L ) denote jump (non-jump) operators and χ a are the counting fields and the system state ρ is not actively controlled (passive measurement). In contrast, when active control operations on the system depending on the measurement outcomes are performed after single jump events, the Fourier-transformed master equation generalizes to d dt ρ( {χ a },t ) [ = L + e iχa eˆk a J a ]ρ ( {χ a },t ) (4) a where the ˆK a describe the control operation []. In contrast to passive measurements, in an active (control) measurement the (scalar) counting fields χ a are thus upgraded to superoperators, iχ a iχ a + ˆK a. Feedback scheme. We choose the feedback of the SET current acting on T c or/and on ε which are accessible in present-days setups by appropriate gate voltages. The corresponding feedback Hamiltonian is defined by

3 3 H fb = a=l/r I a(t)h a with h a Θ aεˆσ z + Θ atˆσ x and the feedback parameters Θ at and Θ aε. Experimentally, single jump events I L/R (t) may neither be resolved nor measured independently in the SET circuit. Therefore, we propose an additional quantum-point contact (QPC) weakly attached to the SET as sketched in Fig. (a). Its telegraph-like current signal I QPC (t) provides the occupation of the SET [Fig. (b)]. For Γ a > Γ a the upwards/downwards steps are related to single jumps into/out of the SET and, consequently, I L (t)/i R (t) are trains of δ-kicks [Fig. (c)]. The superoperators ˆK L/R in the feedback master equation (4) are then given by ˆK a ρ = i [ h a,ρ ]. With feedback into ε/t c, the diagonal matrices B L/R (3) become rotation matrices around the qubit x/z-axis [], so that whenever an electron jumps into/out of the SET the qubit state is rotated by the corresponding angle. Feedback-induced qubit oscillations. As a remarkable effect of this Markovian feedback, we find damped qubit oscillations even when the qubit is initially prepared in a completely mixed state (i.e., in the origin of the Bloch sphere: S = ). For low electron-phonon coupling without feedback [] or for vanishing SET-qubit coupling (Γ = Γ, i.e., the qubit kicked by a random sequence of δ-pulses) the qubit would remain completely mixed. Fig. shows the oscillations in S z (t) for a feedback from I R (t) [I L (t) provides similar results] into (a) ε and (b) T c. The feedback into ε [Fig. (a)] induces oscillations with the frequency ω Q = 4 (8Tc ) γs (ε = ) for 8T c > γ S where γ S = ( Γ Γ ) is the damping due to SET-qubit coupling [] and additional feedback damping which increases with Θ Rε. At large times, the qubit state is again completely mixed. In contrast, for the T c feedback [Fig. (b)] a frequency shift towards higher(lower) frequencies for Θ RT > (<) occurs, and in the long-term limit a localization [S z (t ) ] takes place. The transient oscillations of the qubit occupation [Fig. (c)] can be directly measured as damped current oscillations [compare Figs. (d) and (e)]. Noise spectrum. The total noise spectrum of the SET current is computed by (Ramo-Shockley theorem) S(ω) = 4 a,b=l,r S ab(ω) for a symmetric SET with S ab (ω) = Tr { [J ab +J a Ω(ω)J b +J a Ω( ω)j b ] ρ } where Ω(ω) = [ iω L(,)], steadystate ρwithl(,) ρ =, and the jump operators J a iχa L(χ L,χ R ) (,), J a,b iχa iχb L(χ L,χ R ) (,). The spectrum provides detailed information on the qubit dynamics such as oscillation frequencies and relaxation rates. Without feedback, S(ω) shows peaks at the bare qubit oscillation frequency (for γ S T c ) and at the shifted frequency due to Coulomb interaction for nonvanishing and weak SET-qubit coupling Γ /Γ, see inset of Fig. 3(a). We can further use it to monitor the effect of feedback on the qubit dynamics, as shown in Fig. 3(a). The most interesting feedback scenario is when both I L (t) and I R (t) (a).5 S Ω..5. S Ω Θ T = Ω (b) Sy n, Sz S z S y t FIG. 3: (color online). (a) SET noise spectrum S(ω) vs. ω for various feedback parameter Θ T = Θ LT = Θ RT =.4 (dashed, red),. (dotted, green), (thick solid, black),. (dash-dotted, blue),.4 (thin solid, violet); (Inset) S(ω) for Θ T =, Γ /Γ =.5, U/Γ = 5. (b) Single realization of Bloch vector components S y,z (t) and SET occupation n(t) vs. time Γt with initially mixed state S = ; ensemble averaged quantities (dashed lines); (Inset) Single realizations in the Bloch sphere; measurement time Γτ =., feedback parameter: Θ T =.45. Parameters: ε =, U/Γ =, Γ /Γ = 9, T c/γ =. are simultaneously coupled back into T c. This yields a frequency shift linear in Θ T Θ LT = Θ RT [compare Fig. (b)] of the -peak which is clearly revealed in S(ω) [Fig. 3(a)]. Moreover, at low frequencies the noise exhibits a pronounced local minimum w.r.t. Θ T at small negative values (not shown) which turns out to correspond to the stabilization of pure qubit states. Generation of pure states The simulation of the quantum trajectory in Fig. 3(b) indeed confirms that by starting in a complete statistical mixture, the qubit will finally approach the pure state ( t i b ) (S y = ) after a few electron jumps through the SET. In the following we discuss in more detail the stabilization of such pure qubit states. The steady state yields n =, S y = γ 4cT c 4T c (4T c +sg)+cγ + (G+η p ),(5) with s sin(θ T ), c cos(θ T ), and G ( Γ+ Γ ). The S x -component is solely governed by electron-phonon relaxation [] (we omit the S z -component for brevity). We note that the SET occupation is not affected by the feedback due to the high SET bias limit. The steady-state purity [P Tr( ρ ) = α ( S α ) ], i.e., the length of the Bloch vector for small electron-phonon coupling η p Γ reads P = (γ ) c cg +8T c (4T c +sg) cγ + G+4T c (4T c +sg). (6) First, for zero feedback (Θ T =, c = ) the purity vanishes as expected. Second, above a critical coupling Γ /Γ + 8T c /Γ following from the condition

4 4 T (a) T c. T c ϕ S z π/6 π/4 π/3 π/ π ϕ T c. 5 5 Sy ( ϕ ( ϕ cos t isin b ) ) ϕ=π/ (b) (c) 3 (d) π/ π/ Θ T FIG. 4: (color online). (a) Location of purity P maxima in feedback parameter plane (Θ T, Γ /Γ), Eq. (7). Above a critical coupling Γ /Γ +8T c/γ (bifurcation) the qubit state is pure with P =, otherwise P <. At the bifurcation the qubit state is maximally delocalized with S y = (ϕ = π/) and above S z with < ϕ < π, ϕ π/ holds. The bifurcation becomes blurred and the complete delocalization is diminished for U/Γ =.5 (dashed line) and/or η r/γ =., βt c = (dotted line). (b)-(d) Normalized steady-state Current I/I with I Γ+Γ (red, solid curves), Fano factor F = S()/(I) (blue, dashed curves), purity P = Tr(ρ ) 8 (black, dotted curves) vs. Θ T. Parameters: ε = U =, T c/γ =, Γ /Γ = 4 (b), 9 (c), 5 (d). B (γ ) 6T c cos ( Θ ± T, two maxima of P exist: ) = G[(γ + ) +D+T c ] γ+ D ±8T c B G(ΓΓ +6T c ),(7) with D G + 4Tc where the stationary qubit state has purity P =. Below the bifurcation (B < ) the state is partially mixed (P < ) and cos(θ T ) = G(G +6Tc ) / holds. Fig. 4(a) shows this bifurcation behavior of the maximal purity for various T c. The bifurcation point (B = ) corresponds to the completely delocalized state with S y = (ϕ = π/) The sign of S y is determined by γ, compare Eq. (5). In the upper(lower) branch the pure qubit state becomes increasingly localized with S z > (S z < ). The effect of finite U and the phononic environment onto the steady-state purity is essentially a blurring of the bifurcation in Fig. 4(a) (dashed and dotted curves). A complete delocalization (S y = ) as for the noninteracting case hence cannot be obtained here anymore. Nevertheless, asymptotic purification with increasing localization can be realized in the limit of large SET-qubit coupling. Steady-state SET transport. Without feedback we obtain the steady-state current I /e = (Γ +Γ )/8 and the Fano factor F S() I = ) +(γs+ηp)(γ Γ (8T c) (Γ+Γ ). The SET current is not affected by the phonons, whereas the noise carries information about the interaction of the qubit with the SET detector γ S and the phononic environment η p. The Fano factor is one-half for vanishing SET-qubit coupling (Γ = Γ ) and diverges for large SET-qubit coupling (Γ /Γ ) which results from the bistable current behavior (quantum Zeno effect). It also diverges for T c since the two subspaces decouple in this limit. Figs.4(b)-(d) presenttheθ T -dependenciesofthesteadystate current and the Fano factor at various detectorqubit couplings Γ /Γ. The current, in principle, follows the electron localization in the qubit. At the critical coupling [Fig. 4(c)] when P = the qubit state is maximally delocalized so that the current approaches its no-feedback value I and the Fano factor displays a local maximum. For larger couplings [Fig. 4(d)] the purity maximum in the upper branch is accompanied with a current maximum and a Fano factor minimum. Hence, the degree of the qubit state purity can be very accurately monitored by varying the feedback parameter Θ T. Conclusions. We have presented a simple way to stabilize exact pure states and generate charge oscillations in an individual charge qubit by applying an electronic feedback loop. Even for moderate detector back-action and coupling to a dissipative environment (e.g. phonons) an asymptotic purification can be achieved. We note that in contrast to ground state relaxation just by cooling and strong electron-phonon coupling our method also works at high temperatures and, in principle, arbitrary pure states can be produced. The detector current and its noise spectrum are shown to be sensitive monitors for the effect of feedback. Financial support by the DFG (BR 58/5, BR 58/8, SCHA 646/-, and SFB 9) and helpful discussions with F. Haupt, K. Mosshammer and C. Pöltl are gratefully acknowledged. Electronic address: gerold.kiesslich@tu-berlin.de [] H. M. Wiseman and G. J. Milburn, Quantum Measurement and Control (Cambridge University Press, Cambridge, ). [] T. Fujisawa et al., Science, 634 (6); S. Gustavsson et al., Phys. Rev. Lett. 96, 7665 (6). [3] K. D. Petersson et al., Phys. Rev. Lett. 5, 4684 (). [4] S. A. Gurvitz, Phys. Rev. B 56, 55 (997); A. Shnirman and G. Schön, Phys. Rev. B 57, 54 (998); H.- S. Goan et al., Phys. Rev. B 63, 536 (); A. N. Korotkov, Phys. Rev. B 63, 543 (); S. Pilgram and M. Büttiker, Phys. Rev. Lett. 89, 4 (); N. P. Oxtoby et al., Phys. Rev. B 74, 4538 (6); I. Neder and F. Marquardt, New J. Phys. 9, (7); S. A. Gurvitz and D. Mozyrsky, Phys. Rev. B 77, 7535 (8). [5] Handbook of chaos control, edited by E. Schöll and H. G. Schuster, (Wiley-VCH, Weinheim, 8). [6] R. Ruskov and A. N. Korotkov, Phys. Rev. B 66, 44 (). [7] S.-K. Wang et al., Phys. Rev. B 75, 5534 (7). [8] A. N. Korotkov, Phys. Rev. B 7, 35 (5). [9] H. Bluhm et al., Phys. Rev. Lett. 5, 683 (). [] G. Kießlich et al., Phys. Rev. Lett. 99, 66 (7). [] See supplementary material. [] G. Schaller et al., Phys. Rev. B 8, 457 (9).

5 Supplementary Material: Charge Qubit Purification by an Electronic Feedback Loop G. Kießlich, G. Schaller, C. Emary, and T. Brandes Institut für Theoretische Physik, Hardenbergstraße 36, Technische Universität Berlin, D-63 Berlin, Germany arxiv:.377v [cond-mat.mes-hall] 3 Jun Equation numbers in the main text are referred to by an asterisk. QUBIT DYNAMICS WITHOUT FEEDBACK We rewrite Eqns. (*) as vector equation ẋ = C x+v with C γ+ γ L et z γ R et z n A +D L B R, () n B L A +D R x (n,s,s ) T, v (γ + L,γ+ r e x +γ L e z, γ L e z) T. For its solution we use the Laplace transform x(z) = [ z C ] ( x +v/z ) with the initial state x (n,s x,sy,sz,sx,sy,sz )T. For symmetric coupling (Γ/ Γ L = Γ R, Γ / Γ L = Γ R ) and ε = U = we find ñ(z) = A (z,n )+za (z)+(γ + +n z)a 3 (z) z[a (z,)+(γ +, +z)a 3 (z)] S x (z) = g(z) Sx, Sy (z) = T cs z +zs y 8T c +zg(z), S z (z) = 4T cs y +g(z)s z 8T c +zg(z), () with A (z,n ) ΓΓ [ΓΓ +n z(γ + +z)] (γ ) (γ + +z), A (z) ΓΓ γ + γ [g(z)s z 4T c S y ], A 3 (z) 8Tc +(γ+ +z)(γ + +z), S y/z S y/z S y/z, and g(z) γ + ΓΓ +z. The inverse Laplace transform for ε = U = yields S x (t) = η ( r e γdt) γ S +η r coth(βt c ) +e γdt S x, [ S y/z (t) = e 4 γdt S y/z cosh(ω Q t) γ DS y/z 8T c S z/y ] sinh(ω Q t), (3) 4ω Q with ω Q 4 γ D (8T c ) (for γ D 8T c damped qubit oscillations occur.) and γ D γ S + η p where γ S ( Γ Γ ) is the decoherence rate due to the SET-qubit coupling which was also found for a qubit coupled to a QPC in Ref. [] where the respective rates are given by the QPC transmissions T and T. For U > the expressions become rather cumbersome such that the decay can be non-exponential (see e.g. Ref. [] and discussions therein). The steady state can be readily obtained by lim z + z x(z) and yields (also for ε,u ) n = ( ) T, η r S = γ S +η r coth(βt c ),,, (4) with η r = gπt c e Tc/ωc being the phonon relaxation rate where g is the dimensionless electron-phonon coupling strength and ω c is a cut-off frequency [3]. For zero phonon temperature (β ) the x component becomes y/(+y) (y η r /γ S ) also discussed in Ref. [4]; the qubit relaxes into a pure state for y (qubit ground state). For large temperatures (β ) the qubit will behave as statistical mixture.

6 FEEDBACK-MODIFIED EQUATIONS OF MOTION The feedback scheme modifies Eqns. (*) such that B a B f a, n n f, n n f, and v vf. For Θ at = we find n f = n, n f = n, v f = v, and B f a = cos(θ aε ) Γ a Γ a sin(θ aε ) Γ a Γ a sin(θ aε ) Γ a Γ a cos(θ aε ) Γ a Γ a, (5) γ + R which represents a rotation of the Bloch vector around the z axis with angle Θ aε. Respectively for Θ aε = the matrix Ba f rotates the Bloch vector around the x axis with angle Θ at : Γa Γ Ba f = a cos(θ at ) Γ a Γ a sin(θ at )γ a + sin(θ at ), Γ a Γ a cos(θ at)γ a + ( ) T n f = γ R,sin(Θ RT ),cos(θ RT )+ γ L, γ R ( ) T n f = γ L,sin(Θ LT ),cos(θ LT )+ γ R γ, L ( T v f = γ + L,γ L e z, γ L [sin(θ LT)e y +cos(θ LT )e z ]). (6) For symmetric coupling and ε = U = the x components decouple, and we can truncate the dynamical system to 5 equations: γ + γ γ γ C f s R γ+ η p T ΓΓ γ c c + R s R = γ (+c R) T c γ+ ΓΓ γ s + R c R γ s, (7) ΓΓ γ L c + L s L γ+ η p T c γ (+c L) ΓΓ γ s + L c L T c γ+ with s a sin(θ at ) and c a cos(θ at ). Electronic address: gerold.kiesslich@tu-berlin.de [] S. A. Gurvitz, Phys. Rev. B 56, 55 (997). [] C. Emary, Phys. Rev. A 78, 35 (8). [3] T. Brandes and T. Vorrath, Phys. Rev. B 66, 7534 (). [4] S. A. Gurvitz, L. Fedichkin, D. Mozyrsky, and G. P. Berman, Phys. Rev. Lett. 9, 668 (3).