Full Counting Statistics of Superconductor-Quantum Dot-Superconductor Contacts in the Presence of Interactions

In this paper, we discuss the full counting statistics of superconducting quantum dot contacts. We discuss the effects both of phonon and onsite electronic interaction focusing on the experimentally most relevant case of strong onsite electronic interactions. We find that in general, the Josephson effect and multiple Andreev reflections in these systems are strongly suppressed due to the onsite interaction. However, in case resonant phonons are found, the effect of the onsite interaction can be overcome.


Introduction
Superconducting nanoscopic systems represent one of the most interesting classes of systems [1] [2]. This includes specifically the possibility to study and manipulate entangled quantum systems [3] [4] [5] as well as the possibility to generate even higher order correlations (i.e. multiple Andreev reflections) [6] [7] [8]. Andreev reflection refers to the retroreflection of a hole leaving a Cooper pair behind. Such processes may happen multiple times if the two leads contacted are both superconductors.
The properties of superconductor-quantum dot-superconductor (S-QD-S) systems are particularly relevant and rich due to the manifold energy scales being present in the system (superconductor gap ∆ , phonon frequency D ω , Cou-without overcomplifications [11]).
The highly nonlinear transport characteristics of the S-QD-S setup we wish to analyze have been studied in the non-interacting setup [6] and noise properties have been analyzed in non-equilibrium in less complicated systems [12], but a full analysis of the system in question including all aspects of the current flow is missing so far.

Description of the System
The quantity of our primary interest is the full counting statistics (FCS) in terms of the cumulant generating function (CGF) χ. It represents a very convenient tool for the calculation of a variety of transport properties. It is directly related to the probability distribution function P(Q) to transfer Q elementary charges during a fixed very long measurement time τ. By simple derivation with respect to some parameters (counting fields) χ gives all cumulants (irreducible moments) of P(Q). This allows to measure and understand in depth the different transport properties contributing to the CGF. The Hamiltonian for the system under consideration is given by The superconducting electrodes on the left and right can be described by BCS Hamiltonians with the gap ∆ of the superconducting terminal. We assume both gaps to be identical in order to simplify the analysis as no large-scale new features are to be expected.
In equilibrium the respective Hamiltonians can be written as ( . The voltage is applied symmetrically so that 2 The corresponding Green's functions are given in [13] and are abbreviated as ω ω respectively. The quasi-particle density of states is strongly energy dependent The tunnel Hamiltonian has to take into account the time-dependent phase of the SCs.
For the resonant level model applied here [14] it is convenient to work with dressed electronic states by applying a polaron transformation.
which leads to a Hamiltonian where the electron-phonon interaction characterized by the coupling 0 λ and resonance frequency 0 ϖ is completely absorbed in the tunnel part of the Hamiltonian. We are left with The general solution of the system described by Equation (1) seems out of question and approximations to proceed are generally necessary.
We want to first consider 0 0 λ = (without electron phonon interaction). We want to consider the case of strong onsite interaction and an odd number of electrons on the quantum dot. In this case the most prominent effect is the emergence of a spin 1/2 Kondo resonance around the Fermi edge for temperatures below the Kondo temperature K T . In normal conducting systems the Kondo temperature is directly related to the onsite interaction via [15] 2 where n Γ refers to the tunnel rate between the normal conductor and the quantum dot. Due to the additional energy scale Δ in the problem, two scenarios may occur: for large Δ K T the Kondo resonance couples to the quasiparticles in the superconductor leading to a behavior similar to the one for normal conducting systems [16] whereas for small Δ K T the Kondo resonance is weakly coupled to the superconductors due to the absence of mobile electrons at the Fermi edge [17]. In this case we can greatly simplify the problem as an effective model in this case is the resonant level model with small hybridisations K Γ of the superconductor with the dot (which we assume to be equal for simplicity).
We evaluate the cumulant generating function ( ) , χ λ ϕ using the generalized Keldysh technique [18]. ( ) t λ is the counting field and ϕ is the phase difference of the two superconductors. ( ) , χ λ ϕ in the case of charge transport is given by where ( ) q P ϕ is the probability of charge q to be transferred during a given (long) measurement time τ .
Partial derivations of ( ) , χ λ ϕ with respect to λ give direct access to the cumulants (irreducible moments). For our system the connection to the Hamiltonian is given by to find the CGF as the counting field derivative.

H. Soller Open Journal of Applied Sciences
Compared to the case of contacts between normal metals and superconductors the counting field derivative has normal but also anomalous contributions leading to ( ) 2 2 ln , e . . e . . 2 We proceed by defining the exact-in-tunneling λ-dependent dot Green's functions and the free electrode Green's functions in Keldysh space: Equation (10) can now be integrated with respect to λ to access the CGF. The normal contributions give rise to MARs and quasi-particle tunneling whereas the anomalous part gives rise to Josephson tunneling. We discuss both parts separately.

Josephson Tunneling
We first evaluate the second part of the expression in Equation (10). Following [12] and assuming the limit of the Kondo resonance Γ C being weakly coupled the result for the CGF is The inelastic tunnel processes associated with the emission and absorption of phonons are described by the function ( ) P  being the Fourier transform of the phonon-phonon correlation function. For this correlation function we assume the phonons are thermally distributed, which may be due to coupling to a thermal environment given by the substrate of backgate [19]. The effect of coupling to an external bath can be characterized by an additional coupling constant B γ and in the following we assume the bath to be fully Ohmic. In this case the phonon spectral density has Lorentzian shape The phonon correlation function can now be analytically to be:  .
In the case of small c Γ treated here we find the first harmonic for 2 V < ∆ : which corresponds to the dc and ac Josephson current depending on whether a voltage is applied or not.
Typically we would therefore expect a reduced Josephson current given the low transparency of the contact. We would not expect pronounced features of the phonons given that only the phase would depend on voltage.

Multiple Andreev Reflection
We go over to the evaluation of Equation (10). The calculation can be carried forward as in [12]. In the case of low tunneling transparency the electronic transport can be described by a sequential tunneling picture with a product of tunneling rates for each MAR where F n is the Fermi function.
We observe that multiple Andreev reflections are strongly suppressed as 2 , n a n c T ∝ Γ . The suppression can only be overcome in case conductance is boosted by a resonant phonon.
For voltages above the gap we recover the previously discussed situation for normal conducting systems as in [18].

Conclusions
To conclude, we have derived the cumulant generating function for a superconducting quantum dot involving both Josephson tunneling as well as multiple Andreev reflections. Whereas we discovered that the Josephson effect is typically strongly suppressed in these setups and only slightly affected by resonant phonons, multiple Andreev reflections are strongly affected and show clear signs when resonant phonons are involved.
We expect that these results pave the way towards future usage of these junctions as transistors as they show a clearly nonlinear behavior. We would also believe that the corresponding current and noise characteristics can be observed in experimental setups soon. The corresponding analysis of the correspondence between the model developed here and the experimental results would clearly mark the next step also given the approximations in the Kondo limit taken here.