arxiv:cond-mat/ v1 [cond-mat.mes-hall] 2 Feb 1998

Transcription

1 Transport through an Interacting Quantum Dot Coupled to Two Superconducting Leads arxiv:cond-mat/9828v [cond-mat.mes-hall] 2 Feb 998 Kicheon Kang Department of Physics, Korea University, Seoul 36-7, Korea (October 2, 28) We derive a formula for the current through an interacting quantum dot coupled to two supercouducting leads, using the non-equilibrium Green s function formalism. It is shown that the formula takes an especially simple form, when the Andreev reflections in the junctions are negligible because of large Coulomb repulsion in the dot. Our formalism provides a new framework to investigate many-body effects of mesoscopic systems in the presence of the superconducting leads. For an example, we describe the resonant tunneling of quasiparticles in the presence of strong Coulomb repulsion in the dot. Further, we discuss the boson assisted transport. Multiple sharp peaks are found in the current-voltage curve, which originate from the singularity of BCS density of states and the boson assisted tunneling. PACS numbers: 72..Bg, Hk, 73.4.Gk, r Recent advancement of modern nanofabrication technique has stimulated a lot of interest in the study of electron transport in quantum dot devices,2. In addition to the Coulomb blockade and the resonant tunneling, a novel non-equilibrium Kondo effect is expected to occur in the low temperature regime for an ultra small quantumdotorartificialatomwhichiscoupledtotwometallic leads 3. More recently, transport through a nanoparticle is studied experimentally, where the particle is connected to superconducting (Al) leads rather than normal metallic leads. The level spacing of a superconducting nanoparticle can be either larger or smaller than its bulk superconducting energy gap, depending on the size of the particle and other parameters such as temperature or external magnetic field. When the level spacing of the nanoparticle is sufficiently large, the superconducting orderinthedotdisappearsandtheparticlecanberegarded as an ordinary atom 2. In this paper, we present a formulation of calculating the current through an interacting quantum dot with single level coupled to two superconducting leads (Fig. ). The superconducting leads are characterized by the energy gaps L, R, and the chemical potentials µ L, µ R, where µ L µ R = ev, V being the bias voltage across the leads. The subscript L(R) represents the left (right) lead. Using the non-equilibrium Green s function technique, we derive an expression of the steady state current due to quasiparticle processes. It is shown that the current can be expressed in terms of the single-particle spectral density of states of the dot electron, if the Andreev reflection (AR) processes are negligible. In fact, the AR contributions to the current can be neglected when the Coulomb interaction energy (U) in the dot exceeds both the superconducting gaps L, R and the coupling between the dot and the leads Γ L,Γ R. It is because, on this condition, the 2e-charge fluctuation of the dot is energetically unfavourable. Note that this will be a typical condition of double barrier tunneling structures with sufficiently small quantum dot (see e.g. ). In this case the formulation becomes quite simple, because there is no AC component of the current. To describe the system, we start with the following model Hamiltonian H = H L +H R +H D {d ;d } + ) V k (c k d +d c k, () k L,R where H L, H R and H D represent left, right BCS superconducting leads and the quatum dot, respectively. The last term in the Eq. () corresponds to the coupling between the dot and each of the leads. d (d ) creates (destroys) an electron with spin in the dot, and c k (c k) creates (destroys) an electron with the momentum k and spin in either the left (L) or the right (R) lead. In general, H D includes intra-dot electron-electron interactions, and coupling to an environment like a bosonic bath through the electron-boson interactions 9. As noted above, when U L,R,Γ L,R the Andreev reflection processes are prohibited, and DC current flows by the resonant tunneling of quasiparticles. The current flowing out of the lead α (α = L,R) is written by I α = ē h V k k α dω 2π (G < d,k (ω) G< k,d (ω) ), (2) where G < d,k (ω) is the Fourier transform of the Keldysh Green s function G < d,k c (t) i k ()d (t). By using the Dyson s equation in the Keldysh s formalism, the current due to quasiparticles is given by I α = ē Vk 2 dω [( g < h 2π k (ω) g k > (ω)) G < (ω) k α + g < k (ω)(gr (ω) Ga (ω))], (3) where G < (t) i d ()d (t), and G r and G a are the corresponding retarded and advanced Green s functions,

2 respectively. The Green s functions denoted by g are the unperturbed Green s functions of the superconducting lead. In deriving the Eq.(3) from the Eq.(2), we have neglected terms containing the off-diagonal Green s functions of dot electron in the Nambu representation because those terms should involve Andreev reflection processes, which is negligible in our case. It should be noted that the full Green s functions of the dot represented by G may include the Andreev reflections, but contribution of those can be neglected when the charging energy in the dot is sufficiently large. Therefore the Eq. (3) can be adopted as the formula for the current in the system. Using the relations g < k (ω) = 2πi[ u 2 k f α(e k )δ(ω E k ) +v 2 k ( f α(e k ))δ(ω +E k ) ], (4a) g > k (ω) = 2πi[ u 2 k ( f α(e k ))δ(ω E k ) ] +v 2 k f α(e k )δ(ω +E k ) ], (4b) ( ( where u 2 k = 2 + ǫ k µ α E k µ α ), vk 2 = 2 ǫ k µ α E k µ α ), ǫ k represents the (normal) single particle energy of the lead, E k = µ α ± (ǫ k µ α ) α, and f α is the Fermi distribution function of the lead α, one find I α = 2ie h dωγ S α (ω) [ f α (ω)(g r (ω) G a (ω))+g < (ω) ]. (5) Here Γ S α (ω) = Γ αζ α (ω µ α ), with Γ α = π k α V k 2δ(µ α ǫ k ) being the constant characterizing coupling between the dot and the lead α which is assumed to be independent of the energy, and ζ α (ω) being the dimensionless BCS density of states of the lead α: ζ α (ω) = ω / ω 2 2 α if ω > α, and otherwise. It is straightforward to generallize the Eq. (5) to the system with a quantum dot containing multi-levels. In that case, Γ S α and the Green s functions are replaced by M M matrices, M being the number of the levels, and then I α is given by the trace of the M M matrix in the expression (5). In the steady state, the current is conserved, that is I L +I R =, and thus one can symmetrize the formula for the current as I = I L = I R = 2 (I L I R ) = ie dω {( Γ S L h (ω)f L(ω) Γ R (ω)f R (ω) ) (6) (G r (ω) Ga (ω))+( Γ S L (ω) ΓS R (ω)) G < (ω)}, with the condition of the current conservation (I L = I R ), dω ( Γ S L (ω)+γ S R (ω)) G < (ω) = dω ( Γ S L (ω)f L(ω)+Γ S R (ω)f R(ω) ) (7) (G r (ω) Ga (ω)). The current formula (6) is reduced to the simple form I = 2e dω h Γ S (ω){f L (ω) f R (ω)}ρ (ω), (8) where Γ S (ω) = ΓS L (ω)γs R (ω) Γ S L (ω)+γs R (ω) and ρ (ω) = π ImGr (ω), by using the approximation dω ( Γ S L(ω) Γ S R(ω) ) G < (ω) = dω ΓS L (ω) ΓS R (ω) Γ S L (ω)+γs R (ω) (9) ( Γ S L (ω)f L(ω)+Γ S R (ω)f R(ω) ) (G r (ω) Ga (ω)). This approximation is made by multiplying the factor Γ S L (ω) ΓS R (ω) ontheintegrandsofbothsidesoftheeq.(7), Γ S L (ω)+γs (ω) R andseems tobe reasonablesince ΓS L (ω) ΓS R (ω) isasmooth Γ S L (ω)+γs (ω) R function of ω except for the points ω = µ α ± α and its absolute value is not larger than one. Further, this approximation becomes exact in the small coupling limit, Γ k B T, because both G < and G r G a have sharp δ-function-like peaks at resonant level of the dot. The Eq. (8) is the main result of this paper. Note that the spectral function ρ (ω) should be calculated in the presence of the leads, and it includes resonant tunneling, spin flips, inelastic scattering, etc. For the normal metallic leads, Γ S α(ω) = Γ α, and thus the Eq. (8) is reduced to Meir and Wingreen s formula 3. That is, the Eq. (8) is a generallization of Meir and Wingreen s formula to the case of the superconducting leads. As an example of the quantum dot, let us consider an Anderson impurity. The Hamiltonian of the dot is H D = ε d d +Un n, () where n = d d is the number operator of the dot electron with spin. If one neglects higher order effects such as the Kondo-like correlations, the Green s function can besimplyevaluatedbythemethod ofthe equationofmotionandtruncationofhighordergreen sfunctions 4,4,5. To begin with, let us introduce 2 2 Nambu representation of the (retarded) Green s function where γ = yields G (ω) = γ, γ, () ). The equation of motion for G ( d d ( ) ω ε G ω +ε (ω) = + k L,R V k z ψ k, γ (2) +U z γ (2), γ, 2

3 ( ) ( ) ck where z =, ψk = c, γ (2) = ( ) k n d n d and is the unit matrix. This equation cannot be solved exactly due to the last term of the R.H.S. In order to get an approximate solution of G, we consider the equation of motion for the Green s function γ (2), γ and truncate higher order Green s functions which appear on the equation of motion. Adopting this procedure and taking U limit one can find that G (ω) = [G (ω)] n ω ε +i ( Γ S L (ω)+γs R (ω)), (3) which leads to the resonant-tunneling-like expression of the current I = 4e ( n ) dω[f L (ω) f R (ω)] h Γ S L (ω)γs R (ω) (ω ε ) 2 + ( Γ S L (ω)+γs R (ω)) 2. (4) Thefactor( n )arisesfromstrongcoulombrepulsion. Similar expression to the Eq.(4) has been obtained by Yeyati et al. 6, where the dot Hamiltonian is simply simulated by a single non-interacting level, instead of the Eq. (). Resonance broadening of the spectral function isincludedin theeq.(4). Theeffectofresonancebroadening in the I V characteristic has been discussed in the Ref. 6. It is notable that for a small resonance broadening Γ α α, the current via single level directly reflects the BCS singularity in the density of states of the superconducting lead via the function Γ α (ω) (see Fig.2(a)). While the Eq. (4) is valid for α Γ α, Kondolike correlation becomes important on the other limit, α Γ α. InordertodescribetheKondoeffectproperly, more elaborate calculation, such as non-crossing approximation (NCA) 7, is required in calculating the spectral function of the dot electron. For normal metallic leads, NCA has been used to investigate anomalies of the I V characteristic arising from the Kondo resonances in finite bias voltages 5. The formalism given in this paper can address the problem of the non-equilibrium Kondo effect in the presence of the superconducting leads. Another interesting example is the boson-assisted transport. Consider a quantum dot in which the electrons interact with bosonic modes ω q with the coupling strength g q. The Hamiltonian of such a quantum dot is given by H D = + q ε d d +Un n hω q b qb q + ˆn q g q (b q +b q), (5) where ˆn = n and b q (b q ) creates (destroys) a boson. For the full Hamiltonian H of the Eq. () with the quantum dot of the Eq. (5), a canonical transformation with φ = i q (g q/ hω q )(b q b q) gives the new effective Hamiltonian 9,8 where and H D = H T = H = e iˆnφ He iˆnφ = H L +H R + H D + H T, (6a) k L,R ε n +Ūn n + q hω q b q b q, (6b) V k ( c k d e iφ +e iφ d c k ). (6c) The level position and the Coulomb energy in the dot are renormalized by the electron-boson interaction as ε = ε q g2 q/ hω q, Ū = U 2 q g2 q/ hω q. The phase factors e ±iφ which appear in the tunneling term are related to the boson assisted tunneling. With the transformed Hamiltonian H, the current formula (8) is modified as I = 2e h dω Γ S (ω){f L (ω) f R (ω)} ρ (ω), (7) where ρ is the spectral function of the composite particle denoted by d d e iφ : ρ (ω) = π Im d, d. In general, exact calculation of ρ (ω) is not possible. However, in the limit of weak coupling between the dot and leads, that is for Γ = Γ L +Γ R k B T, the leads can be neglected and the spectral function can be evaluated exactly for that case. We consider the Einstein model, in which all bosons have a same energy hω and are coupled to the electrons in the quantum dot with strength g. In this model the spectral function is given by (see e.g. 8 ) ρ (ω) = e α(2n+) l= I l (2α[N (N +)] /2 ) e βl hω/2 {( n )δ(ω ε l hω ) (8) + n δ(ω ε Ū l hω ) }, where N = /(e β hω ) and α = g 2 /( hω ) 2 are the boson occupation number and the dimensionless electronboson coupling constant, respectively. The renormalized quantum dot level and the Coulomb repulsion can be written as ε = ε α hω, Ū = U 2α hω. I l (x) is the modified Bessel function. Applying the current formula (7) with (8), one can get I V characteristics in the presence of the electron-boson coupling. 3

4 Fig.2 shows the I V characteristics for different values of α. Here we consider the spin degenerate case (ε = ε = ε). In the absence of the electron-boson coupling (Fig.2(a)), the I V curve consists of two sharp peaks at ev = 2( ε + ) and at ev = 2(ε + U + ). These peaks reflect the singularity in the BCS density of states of the leads. This feature, which makes better spectral resolution than the case of the normal metallic leads, have been verified experimentally. Note that there is no thermal broadening in contrast to the case of the normal metallic leads. For normal metallic leads, there are two thermally broadend steps at ev = 2 ε,2(ε + U). For α, the main peaks are located at ev = 2( ε + ),2( ε + Ū + ). Further, we obtain additional side peaks at ev = 2( ε + l hω + ) and at ev = 2( ε+ū +l hω + ) with positive integer l related to boson emission. Boson absorption processes are negligible because k B T hω is taken here. For weak electron-boson coupling, only single boson emission process contributes to the I V curve (Fig2.(b)). As the coupling becomes stronger, multiple side peaks appear (Fig2.(c)). For normal metallic leads, side steps arise from the boson assisted tunneling, but those are not clearly shown because of thermal broadening. In conclusion, we have derived a formula for the current through an interacting quantum dot coupled to two suerconducting leads, when the Andreev reflecction processes are negligible because of the strong Coulomb repulsion in the dot. This formula provides a new framework to study interaction effects in the transport through mesoscopic systems coupled to superconductors. For example, quasiparticle resonant tunneling has been described by an equation of motion method in the limit of large Coulomb repulsion in the dot. Further, current-voltage charactersistics of boson assisted tunneling have been studied. Splitting of BCS singularities, which arises from the boson assisted tunneling, has been shown in the I V characteristics. The author thanks B.I. Min for useful discussions. This work has been supported by the Korea Science and Engineering Foundation. Electronic address : kckang@silicon.korea.ac.kr M. A. Kastner, Rev. Mod. Phys. 64, 849 (992); Phys. Today 46, 24 (993). 2 D. C. Ralph and R. A. Buhrman, Phys. Rev. Lett. 72, 34 (994). 3 S. Hershfeld, J. H. Davies and J. W. Wilkins, Phys. Rev. Lett. 67, 372 (99); Phys. Rev. B 46 (992). 4 Y. Meir, N. S. Wingreen and P. A. Lee, Phys. Rev. Lett. 66, 348 (99). 5 Y. Meir, N. S. Wingreen and P. A. Lee, 7, 26 (993); N. S. Wingreen and Y. Meir, Phys. Rev. B 49, 4 (994). 6 A. L. Yeyati, A. Martin-Rodeo and F. Flores, Phys. Rev. Lett. 7, 299 (993). 7 T. K. Ng, Phys. Rev. Lett. 76, 487 (996). 8 C. A. Stafford, Phys. Rev. Lett. 77, 277 (996). 9 J. König, H. Schoeller and G. Schön, Phys. Rev. Lett. 76, 75 (996); J. König, J. Schmid, H. Schoeller and G. Schön, Phys. Rev. B 54, 682 (996). K. Kang and B. I. Min, Phys. Rev. B 55, 542 (997). D. C. Ralph, C. T. Black and M. Tinkham, Phys. Rev. Lett. 74, 324 (995); C. T. Black, D. C. Ralph and M. Tinkham, Phys. Rev. Lett. 76, 688 (996); D. C. Ralph, C. T. Black and M. Tinkham, Phys. Rev. Lett. 78, 487 (997). 2 P. W. Anderson, J. Phys. Chem. Solids, 28 (959); J. von Delft, A. D. Jaikin, D. S. Golubev and W. Tichy, Phys. Rev. Lett. 77, 389 (996); R. A. Smith and V. Ambegaokar, Phys. Rev. Lett. 77, 4962 (996); K. A. Matveev and A. I. Larkin, Phys. Rev. Lett. 78, 3749 (997). 3 Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 252 (992). 4 G. Czycholl, Phys. Rep. 43, 277 (986). 5 K. Kang and B. I. Min, Phys. Rev. B 52, 689 (995). 6 A. L. Yeyati, J. C. Cuevas, A. López-Dávalos and A. Martín-Rodero, Phys. Rev. B 55, R 637 (997). 7 N. E. Bickers, Rev. Mod. Phys. 59, 845 (987). 8 G. D. Mahan, Many-Particle Physics, 2nd ed. (Plenum, New York 99). FIG.. Schematic diagram of a quantum dot connected to two superconducting leads. The superconducting leads are characterized by the chemical potentials µ L,µ R and the energy gaps L, R. An eletron current I flows from left to right if µ L µ R = ev >. FIG. 2. Current-voltage characteristics (a) without, and (b)(c) with boson assisted tunneling. Here L = R =, Γ L = Γ R = Γ/2 and µ L = ev/2, µ R = ev/2 are considered. Other parameters used here are k BT =., ε = ε = ε = 2, U = 8, hω =.6 in the unit of. Full lines represent the results for the superconducting leads, while dashed lines being corresponding results for the normal metallic leads. 4

5 Quantum Dot I µ L, L µ R, R

6 I =(2e =h) I =(2e =h) I =(2e =h) (a) = (b) = : (c) = : ev = Fig.2