# Non-Equilibrium Josephson and Andreev Current through Interacting Quantum Dots

###### Abstract

We present a theory of transport through interacting quantum dots coupled to normal and superconducting leads in the limit of weak tunnel coupling. A Josephson current between two superconducting leads, carried by first-order tunnel processes, can be established by non-equilibrium proximity effect. Both Andreev and Josephson current is suppressed for bias voltages below a threshold set by the Coulomb charging energy. A -transition of the supercurrent can be driven by tuning gate or bias voltages.

###### pacs:

74.45.+c,73.23.Hk,73.63.Kv,73.21.La^{†}

^{†}: New J. Phys.

, ,

## 1 Introduction

Non-equilibrium transport through superconducting systems attracted much interest since the demonstration of a Superconductor-Normal-Superconductor (SNS) transistor [1]. In such a device, supercurrent suppression and its sign reversal (-transition) are achieved by driving the quasi-particle distribution out of equilibrium by means of applied voltages [2, 3, 4, 5]. Another interesting issue in mesoscopic physics is transport through quantum dots attached to superconducting leads. For DC transport through quantum dots coupled to a normal and a superconducting lead, subgap transport is due to Andreev reflection [6, 7, 8, 9, 10, 11]. Also transport between two superconductors through a quantum dot has been studied extensively. The limit of a non-interacting dot has been investigated in [12]. Several authors considered the regime of weak tunnel coupling where the electrons forming a Cooper pair tunnel one by one via virtual states [13, 14, 15]. The Kondo regime was also addressed [13, 16, 17, 18, 19]. Multiple Andreev reflection through localized levels was investigated in [20, 21]. Numerical approaches based on the non-crossing approximation [22], the numerical renormalization group [23] and Monte Carlo [24] have also been used. The authors of [25] compare different approximation schemes, such as mean field and second-order perturbation in the Coulomb interaction. In double-dot systems the Josephson current has been shown to depend on the spin state of the double dot [26]. Experimentally, the supercurrent through a quantum dot has been measured through dots realized in carbon nanotubes [27] and in indium arsenide nanowires [28].

In this Letter we study the transport properties of a system composed of an interacting single-level quantum dot between two equilibrium superconductors where a third, normal lead is used to drive the dot out of equilibrium. A Josephson coupling in SNS heterostructures can be mediated by proximity-induced superconducting correlations in the normal region. In case of a single-level quantum dot, superconducting correlations are indicated by the correlator , where is the annihilation operator of the dot level with spin . To obtain a large pair amplitude, i.e. the equal-time correlator , at least two conditions need to be fulfilled: (i) the states of an empty and a doubly-occupied dot should be nearly energetically degenerate and (ii) the overall probability of occupying the dot with an even number of electrons should be finite. For a non-interacting quantum dot, i.e. vanishing charging energy for double occupancy, this can be achieved by tuning the level position in resonance with the Fermi energy of the leads, [12]. In this case, the Josephson current can be viewed as transfers of Cooper pairs between dot and leads and the expression of the current starts in first order in the tunnel coupling strength .

The presence of a large charging energy destroys this mechanism since the degeneracy condition is incompatible with a finite equilibrium probability to occupy the dot with an even number of electrons. Nevertheless, a Josephson current can be established by higher-order tunnelling processes (see, for example, [13, 14, 15]), associated with a finite superconducting correlator at different times. The amplitude of the Josephson coupling is, however, reduced by a factor , i.e., the current starts only in second order in , and the virtual generation of quasiparticles in the leads suppresses the Josephson current for large superconducting gaps . In particular, it vanishes for .

The main purpose of the present paper is to propose a new mechanism that circumvents the above-stated hindrance to achieve a finite pair amplitude in an interacting quantum dot, and, thus, restores a Josephson current carried by first-order tunnel processes that survives in the limit . For this aim, we attach a third, normal, lead to the dot that drives the latter out of equilibrium by applying a bias voltage, so that condition of occupying the dot with an even number of electrons is fulfilled even for .

We relate the current flowing into the superconductors to the nonequilibrium Green’s functions of the dot. In the limit of a large superconducting gap, , the current is only related to the pair amplitude. The latter is calculated by means of a kinetic equation derived from a systematic perturbation expansion within real-time diagrammatic technique that is suitable for dealing with both strong Coulomb interaction and nonequilibrium at the same time.

## 2 Model

We consider a single-level quantum dot connected to two superconducting and one normal lead via tunnel junctions, see figure 1.

The total Hamiltonian is given by . The quantum dot is described by the Anderson model , where is the number operator for spin , is the energy level, and is the charging energy for double occupation. The leads, labeled by , are modeled by , where is the superconducting order parameter (). The tunnelling Hamiltonians are . Here, are the spin- and wavevector-independent tunnel matrix elements, and and represent the annihilation (creation) operators for the leads and dot, respectively. The tunnel-coupling strengths are characterized by .

## 3 Current formula

We start with deriving a general formula for the charge current in
lead by using the approach of Ref. [29]
generalized to superconducting leads.
Similar formulae that relate the charge current to the Green’s function of the
dot in the presence of superconducting leads have been derived in previous
works, in particular for equilibrium situations, see e.g.
Refs. [16, 22].
The formula derived below is quite general, as it allows for arbitrary bias
and gate voltages, temperatures, and superconducting order parameters for a
quantum dot coupled to an arbitrary number of normal and superconducting leads.
For this, it is convenient to use the operators
and
in Nambu formalism.
The current from lead is expressed as
^{1}^{1}1Note that but
for .,
with ,
where indicate the Pauli matrices in Nambu space
and the electron charge.
Evaluating the commutator leads to

(1) |

with and the lead–dot lesser Green’s functions that are the Fourier transforms of . In the following, we assume the tunnelling matrix elements to be real (any phase of can be gauged away by substituting ). The Green’s function is related to the full dot Green’s functions and the lead Green’s functions by a Dyson equation in Keldysh formalism: , where is the retarded (lesser) dot Green’s function, and and the lead advanced (lesser) Green’s function. Using this relation and assuming energy-independent tunnel rates , we obtain for the current with

(2) | |||||

(3) |

where , and is the Fermi function, with being the temperature and the Boltzmann constant. The dot Green’s functions and are defined as the Fourier transforms of and , respectively. The two weighting functions and are given by

The terms and involve excitation energies above and below the superconducting gap, respectively. For , only the part of that involves normal (diagonal) components of the Green’s functions contributes, and the current reduces to the result presented in [29]. For superconducting leads, this part describes quasiparticle transport that is independent of the superconducting phase difference. The other part of involves anomalous (off-diagonal) components of the Green’s functions and is, in general, phase dependent. The contribution to the Josephson current stemming from this term is the dominant one in the regime considered in [13, 14, 15]. The excitation energies above the gap are only accessible either for transport voltages exceeding the gap or by including higher-order tunnelling, involving virtual states with quasiparticles in the leads, and, therefore, vanishes for large . In this case , that involves only anomalous Green’s functions with excitation energies below the gap, dominates transport. It is, in general, phase dependent, and describes both Josephson as well as Andreev tunnelling.

In the following we consider the limit , where the current simplifies to

(4) |

with being the phase of and the pair amplitude of the dot that has to be determined in the presence of Coulomb interaction, coupling to all (normal and superconducting) leads and in non-equilibrium due to finite bias voltage.

We now consider a symmetric three-terminal setup with , and , and . The quantities of interest are the the current that flows between the two superconductors (Josephson current) and the current in the normal lead (Andreev current) .

Furthermore, we focus on the limit of weak tunnel coupling, . In this regime, an Josephson current through the dot in equilibrium would be suppressed even in the absence of Coulomb interaction, , since the influence of the superconductors on the quantum-dot spectrum could not be resolved for the resonance condition . This can, e.g., be seen in the exactly-solvable limit of together with , where the Josephson current is with . This provides an additional motivation to look for a non-equilibrium mechanism to proximize the quantum dot.

## 4 Kinetic equations for quantum-dot degrees of freedom

The Hilbert space of the dot is four dimensional: the dot can be empty, singly occupied with spin up or down, or doubly occupied, denoted by , , , with energies , , . For convenience we define the detuning as . The dot dynamics is fully described by its reduced density matrix , with matrix elements . The dot pair amplitude is given by the off-diagonal matrix element . The time evolution of the reduced density matrix is described by the kinetic equations

(5) |

We define the generalized transition rates by , which are the only quantities to be evaluated in the stationary limit. Together with the normalization condition , (5) determines the matrix elements of . Furthermore, in (5) we retain only linear terms in the tunnel strengths and the detuning . Hence, we calculate the rates to the lowest (first) order in for . This is justified in the transport regime .

The rates are evaluated by means of a real-time diagrammatic technique [30], that we generalize to include superconducting leads. This technique provides a convenient tool to perform a systematic perturbation expansion of the transport properties in powers of the tunnel-coupling strength. In the following, we concentrate on transport processes to first order in tunnelling (a generalization to higher orders is straightforward). This includes the transfer of charges through the tunnelling barriers as well as energy-renormalization terms that give rise to nontrivial dynamics of the quantum-dot degrees of freedom.

We find for the (first-order) diagonal rates the expressions . The N lead also contributes to the rates where with being the chemical potential of the normal lead and the Digamma function. Notice that vanishes when or . The superconducting leads do not enter here due to the gap in the quasi-particle density of states. These leads, though, contribute to the off-diagonal rates .

For an intuitive representation of the system dynamics we define, in analogy to [31], a dot isospin by

(6) |

From (5), we find that in the stationary limit the isospin dynamics can be separated into three parts, , with

(7) | |||||

(8) | |||||

(9) |

where is the -direction and is an effective magnetic field in the isospin space. The accumulation term (7) builds up a finite isospin, while the relaxation term (8) decreases it. Finally, (9) describes a rotation of the isospin direction.

## 5 Non-equilibrium Josephson current

In the isospin language the current in the superconducting leads is

(10) |

where the upper(lower) sign refers to the left(right) superconducting lead. The component contributes to the Andreev current, while is responsible for the Josephson current. To obtain subgap transport, we first need to build up a finite isospin component along the -direction, i.e. we need a population imbalance between the empty and doubly occupied dot [(this is generated by the accumulation term in (7)]; second, we need a finite which rotates the isospin so that it acquires an inplane component. In order to have a finite Josephson current (), we need the -component, , of the effective magnetic field producing the rotation to be non zero.

The Josephson current and the Andreev current read

(11) | |||||

(12) | |||||

These results take into account only first-order tunnel processes, i.e. the rates are computed to first order in . The factor ensures that no finite dot-pair amplitude can be established if the chemical potential of the normal lead, , is inside the interval by at least . In this situation both the Josephson and the Andreev currents vanish. On the other hand, this factor takes the value if and the value if . Hence, the sign of the Josephson current can be reversed by the applied voltage (voltage driven -transition). The considerations above establish the importance of the non-equilibrium voltage to induce and control proximity effect in the interacting quantum dot. In figure 2 we show in a density plot (a) and (b) for as a function of the voltage and the level position . Both the control of proximity effect by the chemical potential and the voltage driven -transition are clearly visible.

If the detuning is too large, , it becomes difficult to build a superposition of the states and , which is necessary to establish proximity. As a consequence, the Josephson and the Andreev current are algebraically suppressed by and , respectively. Figure 3 shows the Josephson current as a function of . The fact that the Josephson current is non zero for is due to the term , i.e. of the interaction induced contribution to the -component of the effective field acting on the isospin. The term has a maximum at , which causes this effect to be more pronounced at the onset of transport. The fact that the value of the Josephson current varies on a scale smaller than temperature indicates its nonequilibrium nature.

A -transition of the Josephson current can also be achieved by changing the sign of , as shown in figure 4 where is plotted as a function of the phase difference for different values of the level position. Notice that the current for () is different from zero only due to the presence of the term acting on the isospin.

## 6 Conclusion

In conclusion, we have studied non-equilibrium proximity effect in an interacting single-level quantum dot weakly coupled to two superconducting and one normal lead. We propose a new mechanism for a Josephson coupling between the leads that is qualitatively different from earlier proposals based on higher-order tunnelling processes via virtual states. Our proposal relies on generating a finite non-equilibrium pair amplitude on the dot by applying a bias voltage between normal and superconducting leads. The charging energy of the quantum dot defines a threshold bias voltage above which the non-equilibrium proximity effect allows for a Josephson current carried by first-order tunnelling processes, that is not suppressed in the limit of a large superconducting gap. Both the magnitude and the sign of the Josephson current are sensitive to the energy difference between empty and doubly-occupied dot. A -transition can be driven by either bias or gate voltage. In addition to defining a threshold bias voltage, the charging energy induces many-body correlations that affect the dot’s pair amplitude, visible in a bias-voltage-dependent shift of the -transition as a function of the gate voltage.

## References

## References

- [1] Baselmans J J A, Morpurgo A F, van Wees B J and Klapwijk T M 1999 Nature 397 43
- [2] Volkov A F Phys. Rev. Lett.1995 74 4730
- [3] Wilhelm F K, Schön G and Zaikin A D 1998 Phys. Rev. Lett.81 1682
- [4] Yip S-K 1998 Phys. Rev.B 58 5803
- [5] Giazotto F, Heikkilä T T, Taddei F, Fazio R, Pekola J P and Beltram F 2004 Phys. Rev. Lett.92 137001
- [6] Fazio R and Raimondi R 1999 Phys. Rev. Lett.80 2913; Fazio R and Raimondi R 1999 Phys. Rev. Lett.82 4950
- [7] Kang K 1998 Phys. Rev.B 58 9641
- [8] Schwab P and Raimondi R 1999 Phys. Rev.B 59 1637
- [9] Clerk A A, Ambegaokar V and Hershfield S 2000 Phys. Rev.B 61 3555
- [10] Shapira S, Linfield E H, Lambert C J, Seviour R, Volkov A F and Zaitsev A V 2000 Phys. Rev. Lett.84 159
- [11] Cuevas J C, Levy Yeyati A and Martín-Rodero A 2001 Phys. Rev.B 63 094515
- [12] Beenakker C W J and van Houten H 1992 Single-Electron Tunneling and Mesoscopic Devices(Berlin: edited by H. Koch and H. Lübbig, Springer) p 175–179
- [13] Glazman L I and Matveev K A 1989 JETP Lett. 49 659
- [14] Spivak B I and Kivelson S A 1991 Phys. Rev.B 43 3740
- [15] Rozhkov A V, Arovas D P and Guinea F 2001 Phys. Rev.B 64 233301
- [16] Clerk A A and Ambegaokar V 2000 Phys. Rev.B 61 9109
- [17] Avishai Y, Golub A and Zaikin A D 2003 Phys. Rev.B 67 041301
- [18] Sellier G, Kopp T, Kroha J and Barash Y S 2005 Phys. Rev.B 72 174502
- [19] López R, Choi M-S and Aguado R 2007 Phys. Rev.B 75 045132
- [20] Levy Yeyati A , Cuevas J C, López-Dávalos A and Martín-Rodero A 1997 Phys. Rev.B 55 R6137
- [21] Johansson G, Bratus E N, Shumeiko V S and Wendin G 1999 Phys. Rev.B 60 1382
- [22] Ishizaka S, Sone J and Ando T 1995 Phys. Rev.B 52 8358
- [23] Choi M-S, Lee M and Belzig W 2004 Phys. Rev.B 70 020502(R)
- [24] Siano F and Egger R 2004 Phys. Rev. Lett.93 047002
- [25] Vecino E, Martín-Rodero A and Levy Yeyati A 2003 Phys. Rev.B 68 035105
- [26] Choi M-S, Bruder C and Loss D 2000 Phys. Rev.B 62 13569
- [27] Buitelaar M R, Nussbaumer T and Schönenberger C 2002 Phys. Rev. Lett.89 256801; Cleuziou J-P, Wernsdorfer W, Bouchiat V, Ondarçuhu T, and Monthioux M 2006 Nature Nanotechnology 1 53; Jarillo-Herrero P, van Dam J A and Kouwenhoven L P 2006 Nature 439 953; Jørgensen H I, Grove-Rasmussen K, Novotný T, Flensberg K and Lindelof P E 2006 Phys. Rev. Lett.96 207003
- [28] van Dam J A, Nazarov Y V, Bakkers E P A M, De Franceschi S and Kouwenhoven L P 2006 Nature 442 667; Sand-Jespersen T, Paaske J, Andersen B M, Grove-Rasmussen K, Jørgensen H I, Aagesen M, Sørensen C, Lindelof P E, Flensberg K and Nygård J Preprint cond-mat/0703264
- [29] Meir Y and Wingreen N S 1992 Phys. Rev. Lett.68 2512
- [30] König J, Schoeller H and Schön G 1996 Phys. Rev. Lett.76 1715; König J, Schmid J, Schoeller H and Schön G 1996 Phys. Rev.B 54 16820
- [31] Braun M, König J and Martinek J 2004 Phys. Rev.B 70 195345