# Superconducting electron and hole lenses

###### Abstract

We show how a superconducting region (S) sandwiched between two normal leads (N), in the presence of barriers, can act as a lens for propagating electron and hole waves by virtue of the so-called crossed Andreev reflection (CAR). The CAR process which is equivalent to the Cooper pair splitting into the two N electrodes provides a unique possibility of constructing entangled electrons in solid state systems. When electrons are locally injected from an N lead, due to the CAR and normal reflection of quasiparticles by the insulating barriers at the interfaces, sequences of electron and hole focuses are established inside another N electrode. This behavior originates from the change of momentum during electron-hole conversion beside the successive normal reflections of electrons and holes due to the barriers. The focusing phenomena studied here is fundamentally different from the electron focusing in other systems like graphene pn junctions. In particular due to the electron-hole symmetry of superconducting state, the focusing of electrons and holes are robust against thermal excitations. Furthermore the effect of superconducting layer width, the injection point position, and barriers strength is investigated on the focusing behavior of the junction. Very intriguingly, it is shown that by varying the barriers strength, one can separately control the density of electrons or holes at the focuses.

###### pacs:

74.45.+c, 73.63.-b, 74.25.F-## I Introduction

One of the fundamental aspects of quantum mechanics is non-locality which manifests as non-classical correlations between spatially separated particles which is usually called as quantum entanglement Horodecki et al. (2009). During the last decades a tremendous amount of efforts has concentrated to exploit the entanglement as the main resource for quantum information processing and its technological applications Nielsen and Chuang (2011). In solid state systems among variety of proposed methods to create entanglement between electrons, one of the most accessible sources of entangled electrons is Cooper pair (CP) splitting Lesovik et al. (2001); Recher et al. (2001). In conventional superconductors CPs are composed of two electrons in a spin singlet state which can be used to produce two entangled electrons. The process of spatially separating two electrons of a CP without losing their correlations is commonly called as crossed Andreev reflection (CAR) which takes place in the system of a superconductor (S) coupled to two normal leads (N) Byers and Flatté (1995); Deutscher and Feinberg (2000); Deutscher (2002); Falci et al. (2001); Feinberg (2003); Kalenkov and Zaikin (2007). In this process an incoming electron from one N lead is able to enter to the S region via attracting another electron from the other N lead which leaves a hole inside it.

Unlike Andreev reflection (AR) which takes place locally at any NS interface, CAR has been proved to be hardly detectable in experiments, due to the cancellation of its contribution with that of the elastic co-tunneling (EC) in the cross conductance measurements Russo et al. (2005); Cadden-Zimansky and Chandrasekhar (2006); Cadden-Zimansky et al. (2009); Kleine et al. (2009a). Among proposed setups to detect CAR, one method which has been experimentally verified, uses two quantum dots between each N lead and the superconductor. Hofstetter et al. (2009); Herrmann et al. (2010); Hofstetter et al. (2011); Schindele et al. (2012); Das et al. (2012). In other proposals, the replacement of normal leads with strong ferromagnets (F) or even spin half metal has been suggested Yamashita et al. (2003); Beckmann et al. (2004); Cadden-Zimansky et al. (2007). Alternatively, some other works have focused on the diffusive regime of transport in NSN structures Morten et al. (2006); Brinkman and Golubov (2006); Golubev et al. (2009). Other groups have revealed that the Coulomb interactions, in particular in the quantum dot systems, can lead to the nonlocal Andreev processes Futterer et al. (2009); Michałek et al. (2013). Moreover, very recently it has been suggested to use devices based on graphene or topological insulators and generally the two dimensional materials in which the low energy excitations behave as Dirac fermions Cayssol (2008); Nilsson et al. (2008); Linder et al. (2009); Haugen et al. (2010); Chen et al. (2011); Linder and Yokoyama (2014); Schroer et al. (2015).

In this paper we investigate the focusing effect for both electrons and holes injected locally from an electron source and passing through a planar normal-insulator-superconductor-insulator-normal (NISIN) junction. Transport properties of such ballistic systems and particularly the effect of barriers on CAR and EC signals have been already investigated Falci et al. (2001); Feinberg (2003); Kalenkov and Zaikin (2007). Here we will show any planar ballistic NISIN junction under certain circumstances can lead to the focusing of both electrons and holes. The underlying mechanism for the focusing property here is the combination of CAR or equivalently non-local electron-hole conversion (EHC) induced by superconducting correlations with large momentum transfer due to the potential barriers at the interfaces. The main characteristic of NISIN junctions is depicted schematically in Fig 1 where some semiclassical trajectories of electrons and holes originating from a point source are shown. It clearly shows that electrons and holes trajectories are completely separated since the momentum and velocity directions for electron and hole excitations are parallel and antiparallel, respectively. When electronic waves injected from a point source inside N arrive at an NS interface, they are partially reflected and transmitted as both electrons and holes where the last scattering process occurs due to EHC. In other words there exist four different scattering processes in NISIN junctions including normal and Andreev reflections which take place locally, as well as EC and CAR which correspond to the transmission as an electron and a hole, respectively. In particular for the transmitted waves each EHC process results in the negative refraction of electronic wave function and then successive scattering processes can give rise to the sequence of electron and hole focal points inside N (see Fig. 1).

It must be mentioned that prior to the present work, Gómez et al. have shown that an NSN junction based on graphene can act as a lens for electron and hole excitations Gómez et al. (2012). However in their investigation the lensing property is mostly accompanied by the chiral relativistic characteristics of the charge carriers and the so-called Klein tunneling which occurs at the pn type interfaces in graphene. In fact even without the superconducting correlations, the Dirac systems like graphene has been proven to show negative refraction at the pn interfaces and they can behave as Veselago-like lenses Cheianov et al. (2007); Moghaddam and Zareyan (2010). On the other hand, here we show that in NISIN junctions, the combination of EHC with the effect of barriers can lead to the electron and hole focusing. This phenomenon does not require a special Dirac band structure of graphene or topological insulators and it can be in principle observed in the conventional superconducting devices as well as those based on the Dirac materials.

The paper is organized as follows. After the introduction, in Sec. II, the model consisting of superconductor between normal metals and at the presence of barriers is introduced. The formulation which is based on the Gorkov’s equations and Green’s functions are presented afterwards. In Sec. III we show the numerical results for local particle and current densities. These results reveal the electron and hole focusing phenomena in the NISIN structure. Our findings are followed by the discussion over the basic underlying mechanisms, effect of various parameters and irregularities on the results, and possible experimental realizations. Finally, Sec. IV is devoted to the conclusions.

## Ii Theoretical model and formalism

In order to have a quantitative description of the anticipated focusing effects in NISIN junctions, we exploit the Green’s function formalism. In superconducting structures and for the time-independent situations, the full Green’s function is defined with the following matrix form in the Nambu space and the energy representation,

(1) |

in which

(2) | |||

(3) |

are the so-called normal and anomalous parts of the Green’s function. Here field operators and represent the creation and annihilation of excitations with spin at point and time . The statistical averaging is denoted by . In ballistic regime, the full Green’s function can be obtained from the matrix form of Gorkov equations Fetter and Walecka (2003),

(4) |

with Pauli matrices and operating over the Nambu (electron-hole) space. The single particle Hamiltonian describes a typical metal with quadratic dispersion relation as,

(5) |

in which indicates the Fermi energy. For the superconducting gap we use the sharp interface model in which the gap has a constant finite value inside the superconductor () and vanishes outside it. In fact the more precise spatial form of the pair potential should be obtained using the self-consistency relation of the gap . However due to the presence of barriers at the interfaces, the superconducting gap inside S is only weakly affected by the adjacent N regions and therefore using self-consistency equation to obtain the spatial variations of is not crucial. The potential function accounts for the two barriers at the interfaces induced by atomically thin insulating layers. In the absence of external potential which represents the electron injection, the Gorkov equation can be solved to obtain the unperturbed Green’s function . Then in the presence of external perturbation , the full Green’s function can be obtained using the Dyson equation,

(6) |

From the Green’s function one can calculate all of the local properties of the system, including local densities and currents in the framework of linear response theory. In particular, the electron and hole local density of states (LDOS) can be extracted from the diagonal components of matrix , namely using the relations and . Then assuming weak electron injection from the source at and expanding over , the lowest order corrections of the electrons and holes LDOS read,

(7) | |||

(8) |

where and indicate the normal and anomalous parts of the unperturbed Green’s function. On the same ground one can immediately see that the magnitudes of local current densities, up to linear order in the injected current at the source point, are given by,

(9) | |||

(10) |

## Iii Numerical results and discussion

In the following we will present the numerical results obtained by discretized Green’s function calculation for LDOS and current densities. Using the translational symmetry along vertical direction, the unperturbed Green’s function can be written in the Fourier transformed form,

(11) |

Hence, Eq. (4) reduces to a one dimensional form for which can be solved in general by discretization methods. The only important point is to set the increment small enough in comparison with the Fermi wavelength which is the smallest length scale of the problem. Finally by numerical implementation of the two integrations in Eqs. (11) and (6) the perturbed Green’s function at the presence of the electron source can be obtained.

In order to have observable CAR signal, the superconductor must have width comparable with the superconducting coherence length in which is the Fermi velocity Kalenkov and Zaikin (2007). In fact the nonlocal EHC is efficient only when which increases the probability of CAR. For thinner or very wide superconducting layers the CAR is weakened in favor of EC and local AR, respectively. In typical superconductors with meV, the superconducting correlations extend over m, which gives the typical size of the proposed superconducting lenses. In addition, the presence of potential barriers at the interfaces leads to the further suppression of the local AR and EC through the junction and subsequently increases the efficiency of electron or hole waves focusing mediated by the combination of CAR and normal reflections from the barriers.

In the following the electron and hole focusing property of a superconducting NISIN heterostructure is examined by studying the spatial variations of the induced local density and local currents inside the second N region. We further investigate the influence of width of superconducting region, the distance of the electron injection from the interfaces, and the strength of the barriers at the interfaces on the focusing properties of the superconducting lens. Finally we comment on the experimental feasibility of the theoretical model and discuss over the obtained results and their limitations.

### iii.1 Electron and hole focusing

In order to see the focusing signature of the superconducting layer, the spatial variations of the LDOS for the electrons and holes are shown in Fig. 2. We define the dimensionless excess density of states denoted by in which is the constant density of states of the two dimensional electron gas (2DEG). The width of S layer is assumed to be and the source is placed at . Now the sequences of the focal points where electrons and holes are focused, can be clearly seen in Fig. 2. These focal points are accompanied by oscillations around them which originate from the quantum interference of refracted electronic waves. In the case of electron density due to the EC, there is an extra overall push function which shows oscillatory decreasing dependence with the distance from the second interface. On the other hand, the focuses followed by oscillations around them can be seen in the hole density. As we will discuss in the following, the focuses of either electrons or holes originate from the combination of EHC at the N/S interfaces and normal reflection of bogoliubov quasiparticles inside the superconductor. It must be emphasized that the second process is only possible when the barriers are present. Moreover without any barrier at the interfaces the CAR itself will be suppressed in favor of the AR and EC processes.

It is instructive to describe the appearance of electron and hole focuses based on a semiclassical argument. As illustrated in Fig. 1, the trajectories of propagating electronic waves, originated from the injection point inside N, are refracted at the N/S interfaces when an EHC takes place. It is well-known that the velocity of a hole as an excitations below the Fermi level have opposite direction with respect to its momentum. However the velocity and momentum of the electrons have the same direction. Therefore when an EHC occurs, because of momentum conservation along direction, the electronic waves undergo a negative refraction which is revealed in the classical trajectories as well. The negative refraction here is a reminiscent of the famous Veselago lenses in optics which are made at the interfaces between a conventional dielectric medium and a metamaterial with negative refraction index. Now the refractions caused by EHC and successive normal reflections of the Bogoliubov quasiparticles from the barriers, can result in a sequence of focuses inside the lead . It should be noticed that the presence of barriers are not only crucial for the normal reflections but also they provide the source of the large momentum transfer in the EHC processes. In Fig. 1, for the sake of clarity, only those trajectories which lead to a focal point are shown. However there are other trajectories corresponding to the normal and Andreev reflections back into N as well as the diverging waves of electrons and holes inside N which are not shown here.

The semiclassical trajectories of electrons and holes can be helpful in finding the positions of the focal points. It must be mentioned that although the trajectories are defined by both the source point and the angle of incidence, however the positions of the focuses do not vary by the angle of trajectories, warranted by the simple trigonometric relations. Depending on the number of EHCs during a process (being even or odd), an electron or a hole focus shows up which can be seen from Fig. 1(a) and (b), respectively. When the number of normal reflections back into the superconductor increases, new focuses far from the second interface can be established. Then, in principle, we end up with two infinite sequences of separate electron and hole focuses. Using the trigonometric relations the positions of the electron and hole focuses are obtained as below,

(12) |

One should note that the first electron focus corresponds to while the hole focuses start from or depending on the width being larger or smaller than , respectively.

### iii.2 Electron and hole local current densities

Very similar to the LDOS, the local current densities of electrons and holes denoted by and can indicate the focusing signature of the superconducting layer between two normal electrodes. The intensities of electron and hole currents with the corresponding focal points are shown in Fig. 3. There are also small oscillations around the focuses but not as clear as the Friedel type oscillations revealed in the LDOS spatial variations. Here a single geometry but for different configurations of the barriers is considered, in order to see the effect of barriers more clearly. Figures 3 (a)-(c) correspond to the cases with two barriers, single barrier at the first interface, and finally a single barrier at the second interface with corresponding sketches for the electron and hole trajectories. We see that in the presence of both barriers the first two focuses of the holes are very strong. These two focal points corresponds respectively, to the EHC at the first and second interfaces. Other focuses for the holes as well as those of electrons are much weaker since they originate from higher order processes when successive EHC and normal reflections take place. When one of the barriers are absent then the normal reflections cannot lead to further focuses. Subsequently the EHC is sufficiently strong only at one of the interfaces which leads to a single hole focus without any electron focusing. This suggest a way to control the intensity of different focuses by varying the strength of the two barriers, separately.

As mentioned earlier, the CAR is very effective for S layers with certain range of the widths being on the order of the coherence length (). Hence, we expect that varying the width of S layer, can strongly influence the focusing properties. This can be seen clearly in Fig. 4 where the spatial variations of the electron and hole local currents for various values of are shown. Starting from a thin S layer with , EC processes is dominant and therefore the overall push of the electron current is larger as indicated in Fig. 4(a). By increasing the width the EC signal weakens and due to the dominance of CAR the electron and hole focuses are more intense for the junctions with intermediate widths e.g. . Then further increase in the width, weakens the CAR in favor of local AR and subsequently the intensity of the currents and excess densities around the focuses will be suppressed. For very long junctions only normal reflection and AR could take place and CAR as well as the focusing signature diminish. As another interesting point it must be mentioned that when the width of the junction is smaller than the distance of the injection point from the first interface (), two first hole focuses are very intense. These focuses originate from the EHC at the first and second interfaces, respectively and their positions are and according to Eq. (12). On the other hand, when the width is larger than the distance , the EHC at the first interface cannot give rise to a hole focus inside N, as can be understood from Fig. 1(b). Therefore only one single hole focus located at becomes intense when . Other hole focuses with as well as those of electrons which are related to the higher order processes have much weaker intensities as explained before.

To see the effect of the injection point distance on the EC and CAR signals, Fig. 5 shows the spatial variations of and for some different values of when the width is . The results show that both electron and hole currents with main contributions from EC and CAR decreases if the injection point becomes farther away from the S layer. Moreover as it has been mentioned above, depending on the distance of the injection point being smaller or larger than , there will be a single or two intense hole focuses, respectively. On the same ground, when the injection point is close to the S layer (), the first electron focus which originates from two successive EHC at the two interfaces is very strong. However in the opposite limit of far injection points (), beside EHC processes, two or more reflections of Bogoliubov excitations from the interfaces back into S layer take place to establish the first electron focusing. Therefore the intensity of resulting electron focuses will be much weaker as it is clear from Fig. 5(a).

Before closing this part, we should comment on the effect of the excitation energy on the focusing properties of the NISIN junctions. Throughout the paper the results are presented for the certain value of energy inside the superconducting gap, . However we have checked the dependence of the results on and it has been found that the change in the excess densities and currents due to the excitation energy varied from to is very small. This is a simple consequence of the fact that the superconducting gap and therefore is very smaller than the Fermi energy and all of the wave vectors can be approximated with . So the scattering and focusing properties of the junction are determined only by the geometrical parameters and as well as the Fermi wave vector, almost irrespective of the excitation energy. Above the superconducting gap () the EHC is suppressed and subsequently the focusing patterns disappear when .

### iii.3 The effect of barriers strength

In what follows, we examine the influence of barriers strength in more details. As it has been revealed by our results in the previous part, the presence of at least one potential barrier at the interfaces is crucial for the emergence of hole focusing. On the other hand the electron focusing will be absent unless we have effective barriers at both interfaces. In the absence of potential barriers, the dominant scattering processes are AR and EC where both of them satisfy the momentum conservation. The presence of barriers at the interfaces favors the normal reflection and CAR which involve a finite momentum transfer and simultaneously suppresses the other scattering processes i.e. local AR and EC.

Now as one can see from Fig. 6 when potential barriers are very weak () which leads to very transparent interfaces, CAR and the intensity of focusing is negligibly small. By increasing the potential strength the EHC and resulting CAR becomes more effective at the intermediate values . Nevertheless for very strong barriers the competition between normal reflection and CAR is in favor of the first process. Therefore the intensity of the local hole current around the focuses declines by further increase in the barriers strength and becomes fully suppressed for . Unlike the hole current which reveals a maximum value around , the electron local current, originated from EC, monotonically decreases with in particular in the regime of strong barriers ().

### iii.4 On the experimental feasibility

As it has been already mentioned, we show that a generic superconducting NISIN heterostructure, can focus both electrons and holes. Moreover, the focusing property of the junction is irrespective of the band structure and dispersion relation of the normal leads and the superconductor. In fact there is no need for special band structure or other features like the Dirac spectrum considered in previous studies Gómez et al. (2012). Then devices based on conventional systems like 2DEG as well as graphene and other 2D materials can be used as the normal leads contacted to a superconducting region at the middle to see the lensing phenomena for electrons and holes.

Although the specific band structure is not an important requirement for the superconducting lens studied here, it should be noticed that this electronic lens could work properly under certain circumstances. Nevertheless, we will discuss that these requirements can be fulfilled in currently available experimental setups and subsequently our predictions can be checked experimentally. First of all we assume that the Fermi velocities inside normal and superconducting region are almost equal to each other. Fortunately this assumption can be essentially satisfied in the common superconducting heterostructures based on metals, semiconductors, and even structures based on graphene and topological insulators hee (56); Xu et al. (2007). Another important requirement is the presence of effective potential barriers or equivalently insulating layers at the interfaces. As it has been shown in the previous part, the maximum feasibility of the focusing property occurs for . However, in a practically wide range of barrier strengths, the focusing at least for the holes can be observed. Very common an insulating layer is established due to the oxidation or mismatches in the electronic characters of the normal and superconducting layers. Therefore it seems that controlling the barrier strength would be challenging. However by changing the materials, different barrier strength can be achieved Kleine et al. (2009b); Fritzsche et al. (2009) and even it is possible to have highly transparent interfaces like Al/Au Cadden-Zimansky and Chandrasekhar (2006). Therefore choosing suitable materials for the two normal and the superconducting leads, different barrier strengths for the two interfaces and particularly the extreme cases of NISN or NSIN structures with single barrier can be achieved.

On top of the above mentioned requirements, there are some other parameters including disorder or finite size in -direction which may affect the focusing character of the junction. Although it has not been considered explicitly here, however similar to the graphene-based Veselago lenses, a finite disorder strength, suppresses the focusing signals and the focal points are washed out and broadened by the impurity scatterings Gómez et al. (2012). Hence it is crucial to have clean and almost impurity-free samples to see clear focusing patterns for holes and electrons. On the other hand in NISIN junctions with finite width in vertical direction, the back reflection of electronic waves from the boundaries could lead to further interferences. But as long as this vertical width is large enough compared to , , and , the intensity of reflected waves will be much smaller than the focused waves coming from CAR. It is worth noting that the focuses positions remain unchanged at the presence of a finite vertical width. In fact as we can understand from the classical trajectories shown in Fig. 1, those with small enough angles will not see the boundaries before being focused by the S layer. All of these trajectories cross at the focal points and give rise to an intense focusing. However the trajectories with larger angles will be reflected by the horizontal boundaries and depending on their angles they will pass through different points. Therefore unless for very narrow junction in which most of the electronic waves are reflected, the finite vertical width cannot strongly influence the focusing pattern.

Throughout the paper we have considered zero temperature which simplifies the Green’s function calculations, while at finite temperatures one needs to calculate finite temperature Green’s function which technically speaking needs to sum over Matsubara frequencies. But even without calculations it is possible to discuss over the effect of temperature on focusing, based on physical intuition. To this end we note that according to our model and the obtained results, the excitation energy can only slightly affect the CAR and EC signals. Moreover due to the electron-hole symmetry as a basic character of superconductors, for all values of the excitation energy inside the superconducting gap (), the positions of focuses do not depend on . This can be simply understood from the semiclassical arguments and the fact that wave vectors dependence on is negligible. Now the main effect of temperature is to excite quasiparticles with energies below the thermal energy (). Thus as long as the S layer is in its superconducting state, most of the thermally excited quasiparticles will have subgap energies and all of them will be focused almost in the same way. So we can say in brief that due to the electron-hole symmetry as a basic character of superconductors, the focusing properties of the junction cannot be influenced by thermal excitations when with indicating the critical temperature of the superconductor. As a result there will not be any observable broadening of the focal points due to the finite temperature. This is in fact an advantage of superconducting lens studied here in comparison with Veselago-like lenses in materials with chiral relativistic dispersion like graphene Gómez et al. (2012); Cheianov et al. (2007); Moghaddam and Zareyan (2010). For instance in the pn junctions based on graphene, even when one considers a symmetric potential profile, only electrons and holes at the Fermi level have the same magnitude momenta. Then in graphene-based Veselago lenses, the electron-hole symmetry is broken for the excitations above the Fermi level. Hence, depending on the excitation energy, the focusing takes place in slightly different points which leads to the broadening of the focuses at finite temperatures Cheianov et al. (2007).

We finally comment on the possible experimental setups for the injection of electrons and detection of lensing effects. In order to inject electrons locally from a point source one can either use a tip of scanning tunneling microscope or alternatively a small submicron size gate with a size comparable with the Fermi wavelength or below it. Even for larger electron sources, the focusing would take place but in a broadened form of the focuses inside N. Although we have not considered this case explicitly, however in order to understand the effect of a broad injection point, Fig. 7 shows the density variations at the presence of two point sources separated by in vertical direction. We see that irrespective of extra interferences, we will have well separated focuses corresponding to each of the injection points. Therefore we can simply expect similar imaging signature for more injection points and even a broad electron source. Similar to the injection itself, the spatial variation of the electron and hole local currents and densities inside N can be detected via the scanning probe microscopes (SPMs) which have been already used for imaging electron motions in 2DEG Topinka et al. (2001); Aidala et al. (2007).

## Iv Conclusions

In summary, we study the focusing of electron and hole excitations induced by CAR and multiple reflections from the interfaces of an NISIN junction when incoming electrons are locally injected inside the left metallic contact. It is shown that the barrier potentials at the interfaces, intensifies the focal points due to the CAR, while suppressing the EC signal. In fact, the barrier potentials provide the required abrupt changes in particle momentum during a CAR. In addition, the intensities and the configuration of electron and hole focuses can be controlled by varying the location of electron injector, width of the superconductor, and the potential barriers. These findings suggest promising applications of the superconducting lenses for remote manipulation of the entanglement in spatially separated electron or hole pairs.

###### Acknowledgements.

We would like to acknowledge the encouragement and instructive comments of late Malek Zareyan in the early stages of the work. H.C. and A.G.M. thank the International Center for Theoretical Physics (ICTP) for their hospitality and support during a visit in which part of this work was done. A.G.M. acknowledges financial support from Iran Science Elites Federation under Grant No. 11/66332## References

- Horodecki et al. (2009) R. Horodecki, P. Horodecki, M. Horodecki, and K. Horodecki, Rev. Mod. Phys. 81, 865 (2009).
- Nielsen and Chuang (2011) M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, 10th ed. (Cambridge University Press, New York, USA, 2011).
- Lesovik et al. (2001) G. Lesovik, T. Martin, and G. Blatter, Euro. Phys. J. B 24, 287 (2001).
- Recher et al. (2001) P. Recher, E. V. Sukhorukov, and D. Loss, Phys. Rev. B 63, 165314 (2001).
- Byers and Flatté (1995) J. M. Byers and M. E. Flatté, Phys. Rev. Lett. 74, 306 (1995).
- Deutscher and Feinberg (2000) G. Deutscher and D. Feinberg, Appl. Phys. Lett. 76, 487 (2000).
- Deutscher (2002) G. Deutscher, J. Supercond. 15, 43 (2002).
- Falci et al. (2001) G. Falci, D. Feinberg, and F. W. J. Hekking, Europhys. Lett. 54, 255 (2001).
- Feinberg (2003) D. Feinberg, Eur. Phys. J. B 36, 419 (2003).
- Kalenkov and Zaikin (2007) M. S. Kalenkov and A. D. Zaikin, Phys. Rev. B 75, 172503 (2007).
- Russo et al. (2005) S. Russo, M. Kroug, T. M. Klapwijk, and A. F. Morpurgo, Phys. Rev. Lett. 95, 027002 (2005).
- Cadden-Zimansky and Chandrasekhar (2006) P. Cadden-Zimansky and V. Chandrasekhar, Phys. Rev. Lett. 97, 237003 (2006).
- Cadden-Zimansky et al. (2009) P. Cadden-Zimansky, J. Wei, and V. Chandrasekhar, Nature Phys. 5, 393 (2009).
- Kleine et al. (2009a) A. Kleine, A. Baumgartner, J. Trbovic, and C. SchÃ¶nenberger, Europhys. Lett. 87, 27011 (2009a).
- Hofstetter et al. (2009) L. Hofstetter, S. Csonka, J. Nygård, and C. Schönenberger, Nature (London) 461, 960 (2009).
- Herrmann et al. (2010) L. G. Herrmann, F. Portier, P. Roche, A. L. Yeyati, T. Kontos, and C. Strunk, Phys. Rev. Lett. 104, 026801 (2010).
- Hofstetter et al. (2011) L. Hofstetter, S. Csonka, A. Baumgartner, G. Fülöp, S. d’Hollosy, J. Nygård, and C. Schönenberger, Phys. Rev. Lett. 107, 136801 (2011).
- Schindele et al. (2012) J. Schindele, A. Baumgartner, and C. Schönenberger, Phys. Rev. Lett. 109, 157002 (2012).
- Das et al. (2012) A. Das, Y. Ronen, M. Heiblum, D. Mahalu, A. V. Kretinin, and H. Shtrikman, Nature Commun. 3, 1165 (2012).
- Yamashita et al. (2003) T. Yamashita, S. Takahashi, and S. Maekawa, Phys. Rev. B 68, 174504 (2003).
- Beckmann et al. (2004) D. Beckmann, H. B. Weber, and H. v. Löhneysen, Phys. Rev. Lett. 93, 197003 (2004).
- Cadden-Zimansky et al. (2007) P. Cadden-Zimansky, Z. Jiang, and V. Chandrasekhar, New Journal of Physics 9, 116 (2007).
- Morten et al. (2006) J. P. Morten, A. Brataas, and W. Belzig, Phys. Rev. B 74, 214510 (2006).
- Brinkman and Golubov (2006) A. Brinkman and A. A. Golubov, Phys. Rev. B 74, 214512 (2006).
- Golubev et al. (2009) D. S. Golubev, M. S. Kalenkov, and A. D. Zaikin, Phys. Rev. Lett. 103, 067006 (2009).
- Futterer et al. (2009) D. Futterer, M. Governale, M. G. Pala, and J. König, Phys. Rev. B 79, 054505 (2009).
- Michałek et al. (2013) G. Michałek, B. R. Bułka, T. Domański, and K. I. Wysokiński, Phys. Rev. B 88, 155425 (2013).
- Cayssol (2008) J. Cayssol, Phys. Rev. Lett. 100, 147001 (2008).
- Nilsson et al. (2008) J. Nilsson, A. R. Akhmerov, and C. W. J. Beenakker, Phys. Rev. Lett. 101, 120403 (2008).
- Linder et al. (2009) J. Linder, M. Zareyan, and A. Sudbø, Phys. Rev. B 80, 014513 (2009).
- Haugen et al. (2010) H. Haugen, D. Huertas-Hernando, A. Brataas, and X. Waintal, Phys. Rev. B 81, 174523 (2010).
- Chen et al. (2011) W. Chen, R. Shen, L. Sheng, B. G. Wang, and D. Y. Xing, Phys. Rev. B 84, 115420 (2011).
- Linder and Yokoyama (2014) J. Linder and T. Yokoyama, Phys. Rev. B 89, 020504 (2014).
- Schroer et al. (2015) A. Schroer, P. G. Silvestrov, and P. Recher, Phys. Rev. B 92, 241404 (2015).
- Gómez et al. (2012) S. Gómez, P. Burset, W. J. Herrera, and A. L. Yeyati, Phys. Rev. B 85, 115411 (2012).
- Cheianov et al. (2007) V. V. Cheianov, V. Fal’ko, and B. L. Altshuler, Science 315, 1252 (2007).
- Moghaddam and Zareyan (2010) A. G. Moghaddam and M. Zareyan, Phys. Rev. Lett. 105, 146803 (2010).
- Fetter and Walecka (2003) A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle Systems (Dover Publications, New York, USA, 2003).
- hee (56) Nature 446, 2007 (56).
- Xu et al. (2007) S.-Y. Xu et al., Nature Phys. 10, 943 (2007).
- Kleine et al. (2009b) A. Kleine, A. Baumgartner, J. Trbovic, and C. SchÃ¶nenberger, EPL (Europhysics Letters) 87, 27011 (2009b).
- Fritzsche et al. (2009) J. Fritzsche, R. B. G. Kramer, and V. V. Moshchalkov, Phys. Rev. B 80, 094514 (2009).
- Topinka et al. (2001) M. A. Topinka, R. M. LeRoy, B. J.and Westervelt, S. E. J. Shaw, R. Fleischmann, E. J. Heller, K. D. Maranowski, and A. C. Gossard, Nature (London) 410, 183 (2001).
- Aidala et al. (2007) K. E. Aidala et al., Nature Phys. 3, 464 (2007).