All-electrical generation and control of odd-frequency s-wave Cooper pairs in double quantum dots

All-electrical generation and control of odd-frequency -wave Cooper pairs in double quantum dots

Pablo Burset Department of Applied Physics, Nagoya University, Nagoya 464-8603, Japan    Bo Lu Department of Applied Physics, Nagoya University, Nagoya 464-8603, Japan    Hiromi Ebisu Department of Applied Physics, Nagoya University, Nagoya 464-8603, Japan    Yasuhiro Asano Department of Applied Physics, Hokkaido University, Sapporo 060-8628, Japan    Yukio Tanaka Department of Applied Physics, Nagoya University, Nagoya 464-8603, Japan
July 30, 2019

We propose an all-electrical experimental setup to detect and manipulate the amplitude of odd-frequency pairing in a double quantum dot. Odd-frequency pair amplitude is induced from the breakdown of orbital symmetry when Cooper pairs are injected in the double dot with electrons in different dots. When the dot levels are aligned with the Fermi energy, i.e., on resonance, nonlocal Andreev processes are directly connected to the presence of odd-frequency pairing. Therefore, their amplitude can be manipulated by tuning the level positions. Detection of nonlocal Andreev processes by conductance measurements contributes a direct proof of the existence of odd-frequency pair amplitude and is available using current experimental techniques.


Introduction.— The symmetry analysis of Cooper pairs is a key element in the study of superconductivity. For example, Cooper pairs at conventional BCS superconductors form a spin-singlet even-parity state, where the electrons have opposite spins and are coupled in momentum space by the isotropic -wave channel. A current trend in the study of superconductivity is to engineer unconventional superconductors by breaking down symmetries of a BCS superconductor. Consequently, a new type of pairing emerges which is odd in frequency, i.e., odd under an exchange of the time coordinates Berezinskii (1974a); *Berezinskii_1974R; Balatsky and Abrahams (1992); Tanaka et al. (2012); Eschrig (2015). Plenty of theoretical studies suggest ubiquitous presence of odd-frequency superconductivity at inhomogeneous superconducting systems Tanaka and Golubov (2007); *Tanaka_2007c; Yokoyama (2012); Black-Schaffer and Balatsky (2012, 2013a, 2013b); Sothmann et al. (2014); Ebisu et al. (2015); Crépin et al. (2015); Burset et al. (2015); Asano and Sasaki (2015); Mizushima et al. (2015). Unfortunately, experimental evidence for odd-frequency pair amplitude is very limited. Odd-frequency spin-triplet -wave superconductivity can explain the long-range proximity effect Bergeret et al. (2001); *Bergeret_RMP; Khaire et al. (2010), the intrinsic paramagnetic Meissner effect Tanaka et al. (2005); Yokoyama et al. (2011); Asano et al. (2011); Di Bernardo et al. (2015a), and the subgap structure Di Bernardo et al. (2015b) observed in ferromagnet-superconductor hybrids. However, odd-frequency pairs are mixed with conventional ones and their amplitude is not tunable but accidentally determined by the configuration of magnetic moments realized at the junction. To unambiguously establish the presence of odd-frequency pairing, new proposals that filter odd-frequency pairs and allow to control their amplitude are required.

Figure 1: Proximity-induced superconductivity in a DQD three-terminal device. (a) Schematics of a DQD with level positions and contacted by a superconducting lead S and two normal leads L and R. (b) In a local Andreev process, the electrons of a Cooper pair tunnel through one dot into the normal lead. The amplitude for these processes, (blue dashed lines), is enhanced for symmetric devices, i.e., those where each dot is similarly coupled to the leads. (c) In a nonlocal Andreev process, each electron of the Cooper pair tunnels to a different lead. The nonlocal amplitude, (red solid line), is enhanced in asymmetric devices. On resonance (), is odd in frequency if S is a BCS superconductor.

In this Letter, based on recent experimental progress on double quantum dot (DQD) Cooper pair splitter devices Hofstetter et al. (2009); Herrmann et al. (2010), we propose a setup that allows for the detection and manipulation of odd-frequency pairing without using any magnetic elements. Such pair splitters consist of a DQD independently connected to two normal leads and one superconducting electrode as sketched in Fig. 1(a). We show that the symmetry of the induced pairing from the superconducting lead is broken due to the DQD orbital degree of freedom, thus becoming a superposition of symmetric and antisymmetric orbital states. For a spin-singlet superconductor, induced Cooper pairs that are antisymmetric (symmetric) in DQD space must be odd (even) in frequency according to Fermi-Dirac statistics, since parity and spin rotation symmetries are not broken. Cooper pairs transmitted to the same lead are in an even-frequency state [see Fig. 1(b)]. On the other hand, when the levels are on resonance, only odd-frequency pair amplitude is responsible for the splitting of Cooper pairs into different leads [see Fig. 1(c)]. Antisymmetric odd-frequency singlet pairing is greatly enhanced if each dot is coupled differently to the leads, resulting in a measurable contribution to the conductance. Such a connection between symmetry and microscopic transport processes is a unique feature of our proposal. Additionally, the amplitude of odd-frequency Cooper pairs can be controlled by tunning the DQD level positions on- or off-resonance. Alternatively, using a spin-triplet superconductor, the same geometry can be used for the study of Majorana edge states Nakosai et al. (2013). In such a case, odd-frequency triplet pairs are now transmitted to the same electrode. Our proposal opens a novel direction in the study of exotic Cooper pairing owing to the unique connection between symmetry of the Cooper pair and tunneling processes and due to the tunability of the pair amplitude.

Model.— DQD-based three terminal devices Recher et al. (2001); Chtchelkatchev et al. (2002); Bena et al. (2002); Samuelsson et al. (2003) are an ideal platform for exploring the symmetry of induced pairing. Recent experiments are very well modeled by two-level systems and show an exquisite degree of tunability Herrmann et al. (2010); Hofstetter et al. (2009); Herrmann et al. (2012); Fülöp et al. (2014, 2015). Moreover, strong evidence for splitting of Cooper pairs Hofstetter et al. (2011); Das et al. (2012); Schindele et al. (2012); Deacon et al. (2015), which we shall link to the presence of odd-frequency spin-singlet -wave pairing, has been presented. Here, we consider a system with two quantum levels at positions . In the limit of large level separation at the quantum dots, it describes very well a DQD close to the crossing point of the dot resonances Burset et al. (2011). In the absence of external magnetic fields and spin-orbit coupling terms, we describe the system in the combined Nambu-dot space using spinor fields , where () annihilates (creates) an electron with spin at dot . In the following, (), with , are Pauli matrices in dot (Nambu) space, with identity matrix (). The Hamiltonian of the isolated DQD is given by , with and the inter-dot tunneling rate . Transport properties are characterized by the Green’s function


where denotes or for retarded/advanced or Matsubara Green’s function, respectively, with for temperature , Boltzmann constant and integer . Following the geometry described in Fig. 1(a), we include the normal and superconducting leads as self-energies,


with , , and the tunneling rates between dots and leads. We consider the regime where Kondo and exchange interactions between dots can be neglected. The effect of Coulomb repulsion on each dot is to renormalize the level positions and tunneling rares Cuevas et al. (2001); Burset et al. (2011). We assume that the superconducting region is well described by a constant pair potential and neglect its spatial dependence at the surface of the superconductor. The dimensionless Green’s functions at the edge of the superconducting lead, for a BCS superconductor, are .

class element
Table 1: Symmetry classification of Cooper pairs according to frequency/spin/momentum/dot. From left to right, pairs can be even (E, ) or odd (O, ) under time-reversal (frequency), spin (), momentum (), and dot index (L, R), with S for spin singlet (top rows) and T for triplet (bottom rows). The last column shows the corresponding element of the anomalous Green’s function: , with , for local elements and for the nonlocal ones.

Symmetry of induced pair amplitude.— The uncoupled superconducting lead in Eq. (3) represents an even-frequency spin-singlet -wave superconductor which satisfies , for Matsubara frequency. We analyze the symmetry of proximity induced pair amplitude in the DQD system from the anomalous part of the Green’s function of Eq. (1), , with indexes in dot space . Induced superconductivity in the DQD system thus acquires an extra orbital quantum number. Owing to this symmetry, the elements of are divided into even- and odd-orbital terms. Defining , is the only element with odd parity in dot orbital degree of freedom. To be consistent with Fermi-Dirac statistics, must be odd in frequency. Explicitly, we find 111See Appendix for more details on the definition of differential conductance, analytic expressions for the anomalous Green’s function, a detailed description of the symmetry of induced pair amplitude, and results for the case with only one dot coupled to the superconducting lead.


with and . If is satisfied, we find that and . A complete description of the allowed symmetries in the DQD system is given in Table 1, for both spin-singlet and triplet superconductors.

It is possible to enhance odd-frequency over even-frequency pairing on one of the dots, as it is sketched in Fig. 1(c). To study this effect, we define the local ratios 222It is possible to use real frequencies in Eq. (6) instead of Matsubara frequencies. In order to do so, one must define a Green’s function with the same symmetry in both real and imaginary parts as has on Matsubara frequency (see more details in the Supplementary information).


From Eqs. 4 and 5, we see that nonlocal pair amplitudes are proportional to the inter-dot coupling , which is an essential element of our model [OursymmetryanalysisandconclusionsaretriviallyextendedtoarecentlyproposedmodelfortheDQDwithoutinter-dotcouplingwhichincludesathirdsiterepresentingthesuperconductinglead.][]Dominguez_2015. Moreover, is zero when the dot levels are on resonance (). Therefore, dominant odd-frequency pairing on one of the dots requires left-right asymmetry, which can be achieved setting or . Odd-frequency pairing is suppressed when the DQD levels are out of resonance ( and/or ).

Detection and manipulation of odd-frequency pairing.— We consider two different transport measurements. First, a voltage bias is applied symmetrically to both normal leads [Fig. 2(a)]. This configuration is known as Cooper pair splitter setup and has been used in recent experiments Herrmann et al. (2010, 2012); Das et al. (2012); Schindele et al. (2012). At zero temperature, conductance at lead L is given by


with and quasiparticle tunneling transmission . and are the contributions from local and nonlocal Andreev processes, respectively. For subgap voltages (), conductance is mainly given by Andreev processes while the quasiparticle contribution is almost negligible.

Alternatively, a current can flow through lead L if a voltage is applied to lead R [Fig. 3(a)]. This is the basis for a nonlocal conductance measurement Hofstetter et al. (2009, 2011); Tan et al. (2015) which, at zero temperature, is given by


where represents an electron tunneling process. At zero temperature, transmission probabilities for each process are calculated from the retarded Green’s function as , with and 11footnotemark: 1. Specifically, the transmission probability for Andreev processes reduces to


Consequently, we can connect each microscopic process to a symmetry class. Indeed, local even-frequency pair amplitudes provide the probability amplitude for transmission of the two electrons of a Cooper pair into the same lead, i.e., a local Andreev process sketched in Fig. 1(b). On the other hand, nonlocal components account for the probability amplitude of a process where the electrons of a Cooper pair split into different leads: a nonlocal Andreev process [Fig. 1(c)]. If both dot levels are aligned to the chemical potential, i.e., when the DQD is on resonance with , nonlocal pair amplitudes are odd in frequency. The presence of odd-frequency pairing in the DQD and its connection to a specific microscopic process that has been successfully observed in recent experiments is one of the main conclusions of this work.

Figure 2: Cooper pair splitter configuration. (a) A voltage is applied symmetrically to both normal electrodes which allows us to measure the conductances and (left). Right: In an asymmetric DQD system, odd-frequency nonlocal (red line) can be enhanced on dot L and is mainly due to nonlocal Andreev processes (red arrow). (b,c) Map of the ratio (b) and conductance (c) as a function of applied voltage and level position . inside the white dashed line. (d) Conductance (black solid lines), nonlocal Andreev (red dashed lines), and local processes (blue dot-dashed lines) for (left) and (right). For all plots, , , , , , and .

In an ideal setup, we can choose to uncouple one of the dots from the superconductor setting . Local transmission through that dot, is suppressed and Eq. (6) reduces to . As a result, for subgap energies , odd-frequency pairing becomes dominant at dot R (L) if () is satisfied. Consequently, in the Cooper pair splitter configuration, the conductance at one of the leads, Eq. (7), can be completely dominated by nonlocal Andreev processes, which is a signature of the presence of odd-frequency superconductivity 11footnotemark: 1. Decoupling one of the dots requires careful patterning of the DQD, similarly to recent experiments in graphene Tan et al. (2015). In many other experiments, however, DQD are constructed by electrical confinement from electrodes on quasi-one dimensional materials Hofstetter et al. (2009); Herrmann et al. (2010), as sketched in Fig. 2(a). It is thus challenging to decouple one of the dots from the superconductor. Therefore, we consider in the following. We start with strong left-right asymmetry by setting and at the same time. To exclude double occupancy on the dots, we work in the regime with where a single-particle description of transport at the DQD system is allowed Recher et al. (2001). In Fig. 2(b) we show the ratio on dot L, , as a function of and , with a constant. In agreement with our previous analysis, odd-frequency pairing is dominant on dot L for subgap energies as long as the dot levels are close to the chemical potential, i.e., for . At zero temperature, the applied voltage corresponds to the frequency . The conductance at lead L is enhanced for the same bias voltage regime, as shown in Fig. 2(c). Detailed analysis shows that the conductance is mainly given by nonlocal Andreev processes which stem from induced odd-frequency pairing [red lines in Fig. 2(d)].

A small degree of asymmetry is experimentally inevitable. However, by setting in the previous results, the contribution from local Andreev processes is enhanced and becomes comparable to that of nonlocal processes, making it more difficult to establish a connection between conductance and odd-frequency pairing. For weakly asymmetric setups, with , it is better to perform a nonlocal conductance measurement where odd-frequency induced nonlocal Andreev processes only compete with electron tunneling processes [see Fig. 3(a)]. In principle, the two contributions should cancel each other Deutscher and Feinberg (2000); Chtchelkatchev et al. (2002). In a DQD three terminal device, however, the relative position of the dot levels becomes very important to favor Andreev processes, since they mainly take place when the levels are aligned on resonance. For this condition, the pair amplitude is odd in frequency. Therefore, a positive nonlocal conductance proves the presence of odd-frequency pair amplitude [Nonlocalconductancehasalsobeensuggestedtodetectodd-frequencyspin-tripletpairingattheedgeoftwo-dimensionaltopologicalinsulators.][]Crepin_2015. Setting , we show in Fig. 3(b) a map of as a function of and . Within the black dashed line, conductance is positive, i.e., dominated by nonlocal Andreev processes. As , however, electron tunneling processes become more important and the conductance changes sign. The connection between positive Andreev-dominated conductance and odd-frequency symmetry is explicitly shown in Fig. 3(c). When the dot R is on resonance (red solid lines), the nonlocal conductance is positive for subgap energies (bottom) while odd-frequency pair amplitude is dominant on dot L (top). If the dot R is taken out of resonance, the conductance becomes negative and odd-frequency pair amplitude is suppressed.

Figure 3: Nonlocal conductance measurement. (a) A voltage is applied to lead R allowing to measure the conductance on lead L (left). Right: On resonance, odd-frequency nonlocal Cooper pairs provide a dominant contribution to (red arrow), while electron tunneling dominates out of resonance (blue arrow). (b) Map of as a function of the applied voltage and the position of dot R, , for . inside the black dashed line. (c) Ratio on dot L (top) and nonlocal conductance (bottom) for (red solid line) and (blue dashed line). For all plots, , , , and .

Spin-triplet superconducting lead.— When the central superconducting lead is a one-dimensional spin-triplet -wave superconductor, like the case of a metallic nanowire on SrRuO Maeno et al. (2012), induced pairing amplitude at its edges is odd-frequency triplet -wave represented by , which displays Majorana edge states Asano and Tanaka (2013); Nakosai et al. (2013); Zazunov et al. (2016). On the DQD, , , and are now odd-frequency functions, while is even-frequency pairing (see Table 1). Consequently, conductance measured in the Cooper pair splitter configuration is a very useful tool to study the symmetry of edge states at spin-triplet superconductors. For the perfectly symmetric case, where both and are satisfied, the even-frequency term vanishes and Cooper pairs injected on the same lead maintain the odd-frequency symmetry of the superconducting lead. increases proportionally to the difference between the tunneling rates for left and right dots. For example, if the asymmetry originates from the coupling to the normal leads (superconducting lead), even-frequency component follows [].

Conclusions.— We propose a way to generate odd-frequency spin-singlet -wave Cooper pairs on DQD-based three terminal devices. Due to the orbital degree of freedom in the DQD, symmetry of induced Cooper pairs can be broken, featuring a superposition of even and odd-frequency terms. Each symmetry type, however, is responsible for a different transport process; a feature unique of this setup. For spin-singlet superconductors, nonlocal Andreev processes on resonance are uniquely caused by odd-frequency pairing. Therefore, odd-frequency pairs can be detected from standard conductance measurements in asymmetric devices where the contribution of nonlocal Andreev processes is greatly enhanced. Additionally, it is possible to manipulate the amplitude of odd-frequency Cooper pairs by tuning the position of the dot levels. The situation is reversed if the central electrode is a spin-triplet -wave superconductor. Odd-frequency triplet -wave pairing is now associated to local Andreev processes which are the dominant contribution to the conductance if the dots are symmetrically coupled to the leads.

Acknowledgments.— We thank A. Levy Yeyati, N. Nagaosa, and B. Trauzettel for fruitful discussions and comments. We acknowledge financial support by JSPS International Research Fellowship (P.B.). This work was also supported by a Grant-in Aid for Scientific Research on Innovative Areas “Topological Material Science” (Grant Nos. 15H05853 and 15H05852), a Grant-in-Aid for Scientific Research B (Grant Nos. 15H03686, 26287069, and 15H03525).

Appendix A Appendix A. Transport observables

The Hamiltonian describing an isolated double-dot system reads

with for and , where () is the dot electron annihilation (creation) operator, the charging energy, the energy level, and the hybridization between dots. We model the leads as non-interacting Fermi liquids with Hamiltonian

with () the annihilation (creation) operator for electrons on lead with spin , the pair potential and for and zero otherwise. For a conventional superconductor, pair potential is isotropic in momentum space (we assume it constant) while it depends on the momentum for triplet superconductor. In our effective model, quantum dots are considered almost point-like so electron fields are integrated out and the magnitude of the pair potential for triplet superconductor can also be considered constant. We choose for spin-singlet superconductors and for spin-triplet ones, with constant, and Pauli matrices acting in spin space. Consequently, the system is spin degenerate for both types of superconductors and we work only in Nambu-dot space. Finally, the tunneling Hamiltonian between dots and leads reads as

for and the tunneling amplitudes. We define the tunneling rates as for the normal leads and for the superconducting lead, with and the normal state density of states at the Fermi level of normal and superconducting leads, respectively.

In the following, we consider the regime where Kondo correlations can be neglected and where both dots are close to resonance. In this limit, within Hartree-Fock approximation, Coulomb interactions on each dot give rise to the self-energies . Using the equation of motion technique, the interacting Green’s function reduces to the non-interacting one with renormalized couplings and level positions , where and are the bare coupling and level energies Cuevas et al. (2001); Burset et al. (2011).

Figure 4: Conductance measurements in a DQD three-terminal device. Schematics of DQD geometry where a quasi-one dimensional system is contacted by two normal leads L, R and a central superconducting lead S. Quantum dots are formed by electric confinement between leads and can be independently controlled by side gates . (a) Cooper pair splitter configuration where the normal leads are equally biased. A current flows from S to the normal leads where we can measure the conductances and . (b) Nonlocal conductance measurements require that only lead R is biased. A current can still flow to lead L, generating nonlocal conductance .

The steady state current between dot and lead is defined as Meir and Wingreen (1992)

with and . We have used the Fourier transform of the Keldysh lesser Green’s functions and . All matrices are defined in the Nambu-dot space used in the main text. Following Refs. Cuevas et al., 1996; Sun et al., 1999, the current expression can be written in terms of the retarded Green’s function of the DQD, , defined in the main text. For the sake of completeness, we explicitly write the Green’s function in dot space, namely,

where stands for and we have defined . The Pauli matrices and identity matrix act in Nambu space. For subgap voltages (i.e. when ), the current at normal electrodes is given by

where if then and with . All the microscopic processes that give a contribution to the subgap current are summarized in this important result. Indeed, the current measured at electrode L depends on the local Andreev reflection that takes place at that lead, namely,

where we have used that and , with for and labeling Nambu and dot space, respectively. Also included in the current are nonlocal Andreev processes

and electron tunneling processes

On the other hand, for , the quasiparticle contribution to the current is

for , , and with Nambu matrices .

In the Cooper pair splitter setup, we apply a voltage difference between the superconductor and both normal leads. Differential conductance is thus calculated simultaneously biasing the normal leads with . At zero temperature, it reads

with and .

On the other hand, in a nonlocal conductance measurement, the current flowing through lead L originates from a voltage drop applied only to lead R. Therefore, setting , at zero temperature, we find

is only given by electron tunneling and nonlocal Andreev processes and is positive when the latter are dominant.

Appendix B Appendix B. Anomalous Green’s function

We define the anomalous Green’s function . By choosing Matsubara frequencies, we can analyze each component, namely,


where we have defined and the denominator is given by

It is easy to check that the denominator fulfills . Therefore, the symmetry with respect to frequency is given by the numerator. Local terms and and the symmetric nonlocal term , which is zero if , maintain the symmetry of the superconducting lead determined by . On the other hand, it is clear from Eq. (10)(d) that the symmetry of the anti-symmetric nonlocal pairing is determined by . It represents the presence of induced odd-frequency spin-singlet -wave pairing in the DQD system if is an even function of , as it is the case of BCS spin-singlet -wave superconductors. Alternatively, is an even-frequency function if is odd, which is the symmetry of the edge states of a spin-triplet -wave superconductor. is finite only in the presence of inter-dot coupling (). It also vanishes if and are satisfied simultaneously.

For retarded and advanced Green’s functions, , and the symmetry analysis is not straightforward. Due to the infinitesimal imaginary part, retarded and advanced Green’s functions are complex functions with real and imaginary parts that have different dependence on the energy . It is possible, however, to construct a Green’s function that depends on real energy and has the same symmetry as the Matsubara Green’s function Burset et al. (2015). Namely,


where is a real variable which we associate with the frequency . is also a complex function and its real and imaginary parts have the same symmetry with as has with .

Appendix C Appendix C. Symmetry classification of pair amplitude

Regarding its degrees of freedom, the pair amplitude can be even (E) or odd (O) with respect to the frequency, momentum, and dot orbital component (i.e., L or R). It can also be a spin-singlet (S) or triplet (T) component. Total symmetry is constrained to be antisymmetric under the exchange of all degrees of freedom. Following Ref. Tanaka et al., 2012, we adopt the convention frequency/spin/momentum/orbital. In the effective model, momentum has been integrated out and all allowed symmetries are even in this quantity. We list all possible terms in Table 1. For the superconducting leads, the dot degree of freedom must also be even. We thus classify a conventional BCS superconductor, which is even in frequency and spin-singlet, as ESEE. On the other hand, superconductivity at the edge of a one-dimensional quantum wire on top of an unconventional triplet superconductor is odd in frequency Asano and Tanaka (2013). It is thus classified as OTEE. The pair amplitude induced in the double-dot features two symmetries. The local and symmetric components maintain the symmetry of the superconducting lead and are thus either ESEE or OTEE, depending on the spin state of the Cooper pairs in the superconductor. The antisymmetric component , however, changes the orbital symmetry. Consequently, frequency dependence is also changed to OSEO or ETEO.

Figure 5: Conductance in the Cooper pair splitter configuration for the case with only one dot coupled to the superconductor. (a) Sketch of the setup where the superconducting lead is only coupled to dot R (). (b) Ratio on dot L is greater than 1 within the black dashed line. (c) Conductance on dot L as a function of energy and dot levels . (d,e) Conductance on lead L as a function of (d) or (e) (black lines). Red lines correspond to nonlocal Andreev processes and blue line to local processes. (f) Map of conductance on lead R. (g,h) Conductance (black lines), nonlocal Andreev processes (red dashed lines), and local processes (blue dot-dashed lines) as a function of the (g) or (h). For all plots, , , .

Appendix D Appendix D. Superconductor coupled to one dot only

In this section we consider the ideal case where only one of the dots is coupled to the superconducting lead. In the following, we assume as it is sketched in Fig. 5(a). We only consider the case of an even-frequency spin-singlet -wave superconductor. By canceling the coupling between central lead and dot L, the ratios between odd- and even-frequency components become very simple, namely,

It is interesting to notice that the ratios do not depend on the parameters from dot R, namely, , , and . Additionally, when , we find that . Consequently, odd-frequency becomes dominant on dot L and is suppressed on dot R if . Such behavior is reversed in the opposite regime with .

We show in Fig. 5(b) the ratio on dot L as a function of the energy and level position for . Odd-frequency pairing is dominant for subgap energies in the region . We now study the conductance in the Cooper pair splitter configuration. We set , with to allow for a single-particle description. The conductance at lead L and R are shown in Fig. 5(c) and Fig. 5(f), respectively. Their behavior for subgap energies looks very different, with featuring a strong contribution at , while displays maxima around . We analyze the contribution to the conductance from local and nonlocal Andreev processes in Fig. 5(d,e) for and in Fig. 5(g,h) for . It is clear that nonlocal processes are dominant for (red dashed lines) while the local processes are the main contribution to (blue dot-dashed lines).

As a result, in the extreme case where only one of the dots is connected to the superconducting lead, induced odd-frequency pairing can be maximized on the opposite dot. Consequently, an enhanced conductance that is mainly caused by nonlocal Andreev processes stemming from odd-frequency pair amplitude can be measured on that dot.


  • Berezinskii (1974a) V. L. Berezinskii, “New model of anisotropic phase of superfluid he-3,” JETP Letters 20, 287–289 (1974a).
  • Berezinskii (1974b) V. L. Berezinskii, “Novaya model’ anizotropnoy fazy sverkhtekuchego he3,” Pis’ma Zh. Eksp. Teor. Fiz. 20, 628 (1974b).
  • Balatsky and Abrahams (1992) Alexander Balatsky and Elihu Abrahams, ‘‘New class of singlet superconductors which break the time reversal and parity,” Phys. Rev. B 45, 13125–13128 (1992).
  • Tanaka et al. (2012) Yukio Tanaka, Masatoshi Sato,  and Naoto Nagaosa, “Symmetry and topology in superconductors —odd-frequency pairing and edge states—,” Journal of the Physical Society of Japan 81, 011013 (2012).
  • Eschrig (2015) Matthias Eschrig, “Spin-polarized supercurrents for spintronics: a review of current progress,” Reports on Progress in Physics 78, 104501 (2015).
  • Tanaka and Golubov (2007) Y. Tanaka and A. A. Golubov, “Theory of the proximity effect in junctions with unconventional superconductors,” Phys. Rev. Lett. 98, 037003 (2007).
  • Tanaka et al. (2007) Y. Tanaka, Y. Tanuma,  and A. A. Golubov, “Odd-frequency pairing in normal-metal/superconductor junctions,” Phys. Rev. B 76, 054522 (2007).
  • Yokoyama (2012) Takehito Yokoyama, ‘‘Josephson and proximity effects on the surface of a topological insulator,” Phys. Rev. B 86, 075410 (2012).
  • Black-Schaffer and Balatsky (2012) Annica M. Black-Schaffer and Alexander V. Balatsky, “Odd-frequency superconducting pairing in topological insulators,” Phys. Rev. B 86, 144506 (2012).
  • Black-Schaffer and Balatsky (2013a) Annica M. Black-Schaffer and Alexander V. Balatsky, “Proximity-induced unconventional superconductivity in topological insulators,” Phys. Rev. B 87, 220506 (2013a).
  • Black-Schaffer and Balatsky (2013b) Annica M. Black-Schaffer and Alexander V. Balatsky, “Odd-frequency superconducting pairing in multiband superconductors,” Phys. Rev. B 88, 104514 (2013b).
  • Sothmann et al. (2014) Björn Sothmann, Stephan Weiss, Michele Governale,  and Jürgen König, “Unconventional superconductivity in double quantum dots,” Phys. Rev. B 90, 220501 (2014).
  • Ebisu et al. (2015) Hiromi Ebisu, Keiji Yada, Hideaki Kasai,  and Yukio Tanaka, “Odd-frequency pairing in topological superconductivity in a one-dimensional magnetic chain,” Phys. Rev. B 91, 054518 (2015).
  • Crépin et al. (2015) François Crépin, Pablo Burset,  and Björn Trauzettel, “Odd-frequency triplet superconductivity at the helical edge of a topological insulator,” Phys. Rev. B 92, 100507 (2015).
  • Burset et al. (2015) Pablo Burset, Bo Lu, Grigory Tkachov, Yukio Tanaka, Ewelina M. Hankiewicz,  and Björn Trauzettel, “Superconducting proximity effect in three-dimensional topological insulators in the presence of a magnetic field,” Phys. Rev. B 92, 205424 (2015).
  • Asano and Sasaki (2015) Yasuhiro Asano and Akihiro Sasaki, “Odd-frequency cooper pairs in two-band superconductors and their magnetic response,” Phys. Rev. B 92, 224508 (2015).
  • Mizushima et al. (2015) Takeshi Mizushima, Yasumasa Tsutsumi, Masatoshi Sato,  and Kazushige Machida, “Symmetry protected topological superfluid 3 he-b,” Journal of Physics: Condensed Matter 27, 113203 (2015).
  • Bergeret et al. (2001) F. S. Bergeret, A. F. Volkov,  and K. B. Efetov, ‘‘Long-range proximity effects in superconductor-ferromagnet structures,” Phys. Rev. Lett. 86, 4096–4099 (2001).
  • Bergeret et al. (2005) F. S. Bergeret, A. F. Volkov,  and K. B. Efetov, “Odd triplet superconductivity and related phenomena in superconductor-ferromagnet structures,” Rev. Mod. Phys. 77, 1321–1373 (2005).
  • Khaire et al. (2010) Trupti S. Khaire, Mazin A. Khasawneh, W. P. Pratt,  and Norman O. Birge, “Observation of spin-triplet superconductivity in co-based josephson junctions,” Phys. Rev. Lett. 104, 137002 (2010).
  • Tanaka et al. (2005) Y. Tanaka, Y. Asano, A. A. Golubov,  and S. Kashiwaya, “Anomalous features of the proximity effect in triplet superconductors,” Phys. Rev. B 72, 140503 (2005).
  • Yokoyama et al. (2011) Takehito Yokoyama, Yukio Tanaka,  and Naoto Nagaosa, “Anomalous meissner effect in a normal-metal˘superconductor junction with a spin-active interface,” Phys. Rev. Lett. 106, 246601 (2011).
  • Asano et al. (2011) Yasuhiro Asano, Alexander A. Golubov, Yakov V. Fominov,  and Yukio Tanaka, “Unconventional surface impedance of a normal-metal film covering a spin-triplet superconductor due to odd-frequency cooper pairs,” Phys. Rev. Lett. 107, 087001 (2011).
  • Di Bernardo et al. (2015a) A. Di Bernardo, Z. Salman, X. L. Wang, M. Amado, M. Egilmez, M. G. Flokstra, A. Suter, S. L. Lee, J. H. Zhao, T. Prokscha, E. Morenzoni, M. G. Blamire, J. Linder,  and J. W. A. Robinson, “Intrinsic paramagnetic meissner effect due to -wave odd-frequency superconductivity,” Phys. Rev. X 5, 041021 (2015a).
  • Di Bernardo et al. (2015b) A. Di Bernardo, S. Diesch, Y. Gu, J. Linder, G. Divitini, C. Ducati, E. Scheer, M.G. Blamire,  and J.W.A. Robinson, “Signature of magnetic-dependent gapless odd frequency states at superconductor/ferromagnet interfaces,” Nat. Commun. 6, 8053 (2015b).
  • Hofstetter et al. (2009) L. Hofstetter, S. Csonka, J. Nygard,  and C. Schoenenberger, ‘‘Cooper pair splitter realized in a two-quantum-dot y-junction,” Nature 461, 960–963 (2009).
  • Herrmann et al. (2010) L. G. Herrmann, F. Portier, P. Roche, A. Levy Yeyati, T. Kontos,  and C. Strunk, “Carbon nanotubes as cooper-pair beam splitters,” Phys. Rev. Lett. 104, 026801 (2010).
  • Nakosai et al. (2013) Sho Nakosai, Jan Carl Budich, Yukio Tanaka, Björn Trauzettel,  and Naoto Nagaosa, “Majorana bound states and nonlocal spin correlations in a quantum wire on an unconventional superconductor,” Phys. Rev. Lett. 110, 117002 (2013).
  • Recher et al. (2001) Patrik Recher, Eugene V. Sukhorukov,  and Daniel Loss, “Andreev tunneling, coulomb blockade, and resonant transport of nonlocal spin-entangled electrons,” Phys. Rev. B 63, 165314 (2001).
  • Chtchelkatchev et al. (2002) Nikolai M. Chtchelkatchev, Gianni Blatter, Gordey B. Lesovik,  and Thierry Martin, “Bell inequalities and entanglement in solid-state devices,” Phys. Rev. B 66, 161320 (2002).
  • Bena et al. (2002) Cristina Bena, Smitha Vishveshwara, Leon Balents,  and Matthew P. A. Fisher, “Quantum entanglement in carbon nanotubes,” Phys. Rev. Lett. 89, 037901 (2002).
  • Samuelsson et al. (2003) P. Samuelsson, E. V. Sukhorukov,  and M. Büttiker, “Orbital entanglement and violation of bell inequalities in mesoscopic conductors,” Phys. Rev. Lett. 91, 157002 (2003).
  • Herrmann et al. (2012) L. G. Herrmann, P. Burset, W. J. Herrera, F. Portier, P. Roche, C. Strunk, A. Levy Yeyati,  and T. Kontos, “Spectroscopy of nonlocal superconducting correlations in a double quantum dot,”  (2012), arXiv:1205.1972.
  • Fülöp et al. (2014) G. Fülöp, S. d’Hollosy, A. Baumgartner, P. Makk, V. A. Guzenko, M. H. Madsen, J. Nygård, C. Schönenberger,  and S. Csonka, “Local electrical tuning of the nonlocal signals in a cooper pair splitter,” Phys. Rev. B 90, 235412 (2014).
  • Fülöp et al. (2015) G. Fülöp, F. Domínguez, S. d’Hollosy, A. Baumgartner, P. Makk, M. H. Madsen, V. A. Guzenko, J. Nygård, C. Schönenberger, A. Levy Yeyati,  and S. Csonka, “Magnetic field tuning and quantum interference in a cooper pair splitter,” Phys. Rev. Lett. 115, 227003 (2015).
  • Hofstetter et al. (2011) L. Hofstetter, S. Csonka, A. Baumgartner, G. Fülöp, S. d’Hollosy, J. Nygård,  and C. Schönenberger, “Finite-bias cooper pair splitting,” Phys. Rev. Lett. 107, 136801 (2011).
  • Das et al. (2012) Anindya Das, Yuval Ronen, Moty Heiblum, Diana Mahalu, Andrey V Kretinin,  and Hadas Shtrikman, “High-efficiency cooper pair splitting demonstrated by two-particle conductance resonance and positive noise cross-correlation,” Nature Communications 3, 1165 (2012).
  • Schindele et al. (2012) J. Schindele, A. Baumgartner,  and C. Schönenberger, “Near-unity cooper pair splitting efficiency,” Phys. Rev. Lett. 109, 157002 (2012).
  • Deacon et al. (2015) R S Deacon, A Oiwa, J Sailer, S Baba, Y Kanai, K Shibata, K Hirakawa,  and S Tarucha, ‘‘Cooper pair splitting in parallel quantum dot josephson junctions,” Nature Communications 6, 7446 (2015).
  • Burset et al. (2011) P. Burset, W. J. Herrera,  and A. Levy Yeyati, “Microscopic theory of cooper pair beam splitters based on carbon nanotubes,” Phys. Rev. B 84, 115448 (2011).
  • Cuevas et al. (2001) J. C. Cuevas, A. Levy Yeyati,  and A. Martín-Rodero, “Kondo effect in normal-superconductor quantum dots,” Phys. Rev. B 63, 094515 (2001).
  • (42) See Appendix for more details on the definition of differential conductance, analytic expressions for the anomalous Green’s function, a detailed description of the symmetry of induced pair amplitude, and results for the case with only one dot coupled to the superconducting lead.
  • (43) It is possible to use real frequencies in Eq. (6) instead of Matsubara frequencies. In order to do so, one must define a Green’s function with the same symmetry in both real and imaginary parts as has on Matsubara frequency (see more details in the Supplementary information).
  • Domínguez and Levy Yeyati (2016) Fernando Domínguez and Alfredo Levy Yeyati, “Quantum interference in a cooper pair splitter: The three sites model,” Physica E: Low-dimensional Systems and Nanostructures 75, 322 – 329 (2016).
  • Tan et al. (2015) Z. B. Tan, D. Cox, T. Nieminen, P. Lähteenmäki, D. Golubev, G. B. Lesovik,  and P. J. Hakonen, “Cooper pair splitting by means of graphene quantum dots,” Phys. Rev. Lett. 114, 096602 (2015).
  • Deutscher and Feinberg (2000) G Deutscher and D Feinberg, “Coupling superconducting-ferromagnetic point contacts by andreev reflections,” Applied Physics Letters 76, 487–489 (2000).
  • Maeno et al. (2012) Yoshiteru Maeno, Shunichiro Kittaka, Takuji Nomura, Shingo Yonezawa,  and Kenji Ishida, “Evaluation of spin-triplet superconductivity in srruo,” Journal of the Physical Society of Japan 81, 011009 (2012).
  • Asano and Tanaka (2013) Yasuhiro Asano and Yukio Tanaka, “Majorana fermions and odd-frequency cooper pairs in a normal-metal nanowire proximity-coupled to a topological superconductor,” Phys. Rev. B 87, 104513 (2013).
  • Zazunov et al. (2016) A. Zazunov, R. Egger,  and A. Levy Yeyati, “Low-energy theory of transport in majorana wire junctions,”  (2016), arXiv:1603.02969.
  • Meir and Wingreen (1992) Yigal Meir and Ned S. Wingreen, “Landauer formula for the current through an interacting electron region,” Phys. Rev. Lett. 68, 2512–2515 (1992).
  • Cuevas et al. (1996) J. C. Cuevas, A. Martín-Rodero,  and A. Levy Yeyati, “Hamiltonian approach to the transport properties of superconducting quantum point contacts,” Phys. Rev. B 54, 7366–7379 (1996).
  • Sun et al. (1999) Qing-feng Sun, Jian Wang,  and Tsung-han Lin, “Resonant andreev reflection in a normal-metal–quantum-dot–superconductor system,” Phys. Rev. B 59, 3831–3840 (1999).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description