Anomalous Inner-Gap Structure in Transport Characteristics of Superconducting Junctions with Degraded Interfaces.

Quantitative description of charge transport across tunneling and break-junction devices with novel superconductors encounters some problems not present or not as severe for traditional superconducting materials. In this work, we explain unexpected features in related transport characteristics as an effect of a degraded nanoscaled sheath at the superconductor surface. A model capturing the main aspects of the ballistic charge transport across hybrid superconducting structures with normally conducting nanometer-thick interlayers is proposed. The calculations are based on a scattering formalism taking into account Andreev electron-into-hole (and inverse) reflections at normal metal-superconductor interfaces as well as transmission and backscattering events in insulating barriers between the electrodes. Current-voltage characteristics of such devices exhibit a rich diversity of anomalous (from the viewpoint of the standard theory) features, in particular shift of differential-conductance maxima at gap voltages to lower positions and appearance of well-defined dips instead expected coherence peaks. We compare our results with related experimental data.

Background

The Bardeen-Cooper-Schrieffer (BCS) theory of superconductivity (S) is based on the assumption that electrons in the BCS state are bound into Cooper pairs by a sufficiently weak attractive interaction between them. These pairs are correlated, so that, in order to break a pair, one has to change energies of all other pairs. It means that, unlike in the normal-metal (N) state, there is an energy gap Δ for single-particle excitations. The most complete and convincing evidence for the appearance of the energy gap came from the well-known Giaever effect of electron tunneling through N-I-S and S-I-S heterostructures where I stands for a nanometer-thick insulating layer that permits only single-tunneling processes [1]. When an N-I-S trilayer with a conventional s-wave paring-symmetry superconductor is biased with a voltage V, an electron can be injected into the S side merely if its energy exceeds the gap value Δ. It means that at the environment temperature T → 0, the tunneling current I across the junction is vanishing below VΔ = Δ/e. The nonlinearities in I-V curves can be more clearly seen in the quasiparticle differential conductance dI/dV-versus-V which exhibits coherence peaks at V = ± V for an N-I-S structure and at V = ± (Δ + Δ)/e for an S1-I-S2 trilayer with conventional s-wave superconductors. This fact can be used as direct measure of the superconducting-order parameter magnitude and its dependence on external parameters.

Nowadays, electron tunneling spectroscopy is a well-developed technique aimed to provide extensive information about energy spectra of conducting electrodes and an insulating interlayer between them [1]. It has been a disappointment that tunneling spectroscopy has so far given us no such evidence for novel superconductors. Our work was stimulated by systematic deviations of transport characteristics measured for nominally trilayered metal-barrier-metal heterostructures based on novel superconductors from related theoretical predictions for conventional tunnel junctions with a transmission probability D < < 1 [2], as well as for S-c-S devices with D ≤ 1 where c stands for a constriction [3].

The first finding of such kind was revealed for tunnel junctions with some transition metals, especially for Nb-based devices and intensively discussed in the 1970s long before the discovery of high-temperature superconductivity. We mean the so-called knee, a sudden current decrease in low-temperature I-V characteristics of Nb-based S-I-S junctions at voltages slightly just above a strong increase of the current that was expected to occur at the sum of VΔ values for the two superconducting electrodes [4] (see also [5]). Moreover, the current values at very low temperatures and at voltages lower than the sum of VΔ were considerably larger than theoretically anticipated for the BCS densities of states [6]. Later, it was shown experimentally that this anomaly may be caused by the presence of a normal conducting layer on the superconducting niobium film [6, 7]. This statement was confirmed theoretically [8] in the framework of the proximity-effect model which treats incoherent single-particle scatterings from a normal layer to a superconducting one within a tunnel-barrier approximation [9]. Below, we discuss the origin of the “knee” phenomenon analyzing ballistic transport in tunneling S-I-S structures where each superconducting electrode is covered by a nanometer-scale degraded sheath.

The next puzzling experimental feature was the so-called peak-dip-hump structure clearly observed in conductance-versus-voltage curves for break junctions based on high-Tc BSCCO samples [10]. In these experiments, the peak was attributed to twice the gap magnitude although the values obtained were systematically lower than those provided by other experimental techniques, see Table 3.1 in [10]. But the main problem of the measurements was the presence of a pronounced dip that was revealed in data obtained by STM or ARPES [11] as well. Its derivation has not been clear and, in particular, was believed to represent a physical quantity responsible for pairing in the studied compound [12]. It was suggested [13] that the dip appears at the voltage bias (2Δ + Ω)/e where Ω is the energy of some collective excitations that performs a role of glue for paired electrons. The author of ref. [10] studied how the estimated doping level influences the peak and observed a large spread of the measured peak positions. Below, we reproduce numerically the “peak-dip-hump” structure and propose our explanation of each of the three elements relating them to the presence of a degraded layer on the surface of the high-Tc superconductor. It should be noticed that at early stages of single-electron tunneling experiments, four-layered N-I-n-S heterostructures with a normal (n) interlayer between an ordinary superconductor and the barrier were fabricated and studied. Comparing related data with our simulations, we show that, at least qualitatively, the measured conductance spectra of S-I-n-S samples having the n interlayer for sure do follow our predictions.

The last unusual feature which will be analyzed in the paper relates the inner-gap structure in S-c-S junctions where c stands for a constriction with a transmission probability D ≤ 1. The authors of ref. [14] fabricated two types of symmetric S-c-S trilayers, Al-/InAs-nanowire-Al and Nb-InAs-nanowire-Nb junctions and measured the differential resistance dV/dI as a function of the voltage bias at various temperatures. At 0.4 K, the Al-based devices exhibited subharmonic gap features at V = 2ΔAl/(en); n is an integer. Due to the theory [3], such features are the manifestation of multiple Andreev reflections within the constriction between two identical superconductors and they should expose themselves as local minima in the resistance-versus-voltage dependence as it was indeed observed for Al-based samples, Fig. 1a in ref. [14]. Similar dV/dI measurements at T < 4 K for a typical Nb-based junction depicted in Fig. 1b in ref. [14] revealed a subharmonic gap structure corresponding to 2ΔNb/e and ΔNb/e, but they were peaks rather than expected dips. Moreover, it was not the only observation of the difference between the two types of junctions. As was stressed in ref. [14], see related references therein, in the literature, there is a large number of findings with resistance peaks at V = 2ΔNb/(en) in dV/dI curves for Nb-based S-c-S junctions and at the same time dips in dV/dI curves for Al-based structures. Below, we explain the discrepancy as a result of the degradation of Nb-c and c-Nb interfaces.

Our main results which throw some light on unexpected features of transport characteristics of superconducting heterostructures with degraded interfaces are summarized in the “Conclusions” section.

Methods

General Considerations

We start with a model N-I-n-S system which consists of a normal counter-electrode, a nanometer-thick insulating barrier, a superconducting electrode studied, and a nanoscale clean non-superconducting (n) interlayer between I and S films. Confinement of electrons in the n-metal film results in discrete quantum-well states that can be probed directly by single-electron tunneling spectroscopy. Below, we provide an intuitive geometric picture for the localized resonances that play a decisive role in the formation of related transport characteristics. The important impact of the bound states on the charge transport across S-N-S Josephson junctions and ability to probe them by tunneling experiments are well known (see the review [15]). In this work, we show that some unusual features revealed in novel superconductors by single-electron tunneling and break-junction techniques can be understood as an effect of the bound states within a degraded nanoscale sheath at the superconductor surface.

To calculate energies E of the bound states in ballistic junctions, we should take into account the principal difference between backscattering processes at I-n and n-S interfaces. In the first case, an electron (hole) incident on the I-n interface from the n side is retroreflected into an electron (hole) of the same energy E and of the same absolute value of momentum but travelling in the opposite direction to the incoming charge. On the contrary, the reflected from the n-S interface quasiparticle has the same energy E and almost the same momentum is travelling in the opposite direction, and its charge is opposite in sign to that of the incident quasiparticle. Such process is known as Andreev reflection [16], an elastic quasielectron-into-quasihole transformation of Bogoliubov quasiparticles and inverse (with a missing charge of 2e absorbed into the superconducting ground state as a Cooper pair).

Taking it into account, we can easily understand the origin of Andreev bound states within the energy gap that is formed in the n interlayer of the thickness dn and find their energies Ē from the demand of coherent superposition of reflected from the I-n and n-S interface quasiparticle waves. An electron “round-trip” inside the n interlayer consists of four scattering processes, an Andreev electron-into-hole transformation at the n-S interface, a specular hole-into-hole reflection at the I-n boundary, a hole-into-electron transformation at the n-S interface, and an electron-into-electron reflection at the I-n boundary, and four passages across the n interlayer, two of them as a quasielectron and two of them as a quasihole with the phases accumulated and , respectively; ke(h) is an electron (hole) wave vector. Note that a hole is moving in the direction opposite to that of its wave vector. The reflection coefficients can be calculated from the boundary conditions for the wave functions. Each backscattering from the insulating layer provides a phase shift of π whereas the Andreev reflection within the energy gap of an s-wave superconductor contributes an additional phase shift χeh(he)(ε) = − arccos(ε/Δ), where ε = E − EF and EF is the Fermi energy. Adding the phases accumulated along an electron “round-trip” in the n interlayer with two subsequent Andreev reflections, we get the following expression for the phase shiftwhich is valid for an s-wave superconducting electrode, θ is the incident angle, and vF is the Fermi velocity.

The Bohr-Sommerfeld quantization rule requires an electron wave-function phase accumulation Δφ along an enclosed propagating loop inside the n interlayer to be an integer multiple of 2π. Thus, the lowest bound level follows from the relationvalid when dn ≠ 0. For vanishing dn, we get the classical result , i.e., a well-pronounced singularity at the energy gap value. When the thickness dn is finite, . It is clear from Eq. (2) that its effect on current-voltage characteristics is determined by the ratio α = 2dnΔ/(ℏvF) = kFdn(Δ/EF) = dn/ξ where ξ is the relevant length scale which is determined by electronic characteristics of the n-metal and the gap value of the S-electrode in contact. In the most traditional superconductors, like Pb and Sn, due to comparatively large vF and small Δ, ξ ≥ 100 nm and, hence, the effect of the n interlayer can be observed merely for dn of the order of tens of nanometers. Such interlayers can be introduced within the junction only artificially. If however the Fermi velocity is small and at the same time the energy gap is large, as it is in transition metals and novel high-Tc compounds, a completely different situation is expected. In this case, degradation of the superconductor surface on the length scale of the order of several nanometers which is typical for these materials can result in crucial modifications in related transport characteristics. We shall discuss these changes and the information about the S-electrode spectra that can be extracted from such measurements.

Charge Scattering Characteristics

In the following, we approximate the spectrum of electrons in all conducting layers by parabolic bands and limit ourselves to the planar geometry of the studied heterostructures with an x-axis normal to interfaces. Next, we assume that the electron wave function may be factorized into in-plane and tunneling-direction components and will discuss only the latter one considering the in-plane momentum k∥ constant. For a non-superconducting metallic interlayer, it reads as , where is the wave vector of an excitation with the energy ε, which is a solution of the equation , the sign ± corresponds to electron (e) and hole (h) excitations, respectively, EF is the Fermi energy, and m is the quasiparticle mass. To take into account the finite value of the mean free path of a quasiparticle excitation ln in a nanometer-thick n interlayer, we have introduced an additional imaginary term ±i/(2ln) in the wave vector of an electron (a hole) in the n interlayer, see also [17]. Last, for the sake of simplicity, we assume that Fermi energies in the conducting layers, including the superconducting one, are the same.

Notice that the stepwise approximation for the superconducting pair potential in the n-S bilayer used below and known as a rigid-boundary condition is not self-consistent. According to Likharev [18], deviation of the self-consistent solution from the steplike function strongly decreases when the interface resistivity is much bigger than that of metal electrodes. It is very difficult to estimate the thickness of the degraded surface layer as well as the transparency of its interface with the bulk. Thus, our results can only claim a qualitative explanation of the anomalous in-gap features.

To get analytical expressions for elementary scattering amplitudes in the N-I-n-S structure, we model barrier I by a short-range repulsive potential of a rectangular shape U(x) = U0 and the thickness dI. The transparency of such barrier equals to [19] where is the decaying depth of a quasiparticle wave function. When κdI < < 1, it reads DI = (1 + (kF2 + κ2)2dI2/(4kF2))− 1. Introducing a dimensional parameter Z = (kF2 + κ2)dI/(2kF), we get a simple Lorentzian DI = 1/(1 + Z2) (in the limit κ > > kF, it coincides with related results of the paper [2]). Within the same κdI < < 1 approximation, we can derive reflection and transmission probability amplitudes for an electron (hole) to be transmitted through the barrier

The elastic scattering process by which normal current is transferred to a supercurrent at the interface between a superconductor and a non-superconducting metal, the Andreev effect [16], lies in the fact that an electron (hole) incident on the interface from the normal side is retroreflected into a hole (electron). In the general case, the related reh (an electron is scattered into a hole) and rhe (a hole-electron transformation) scattering characteristics look as in [20]where for |ε| > Δ(θ) and for |ε| < Δ(θ). Here, Δ(θ) is an order parameter which is constant Δ for s-wave pairing realized in all conventional superconductors, and Φ(θ) is the order parameter phase. In high-Tc compounds, it is widely accepted that such superconductors have a -wave pairing symmetry. The energy gap of such a pairing state manifests a sign change at some directions of the Fermi wave vector, and its angle dependence is Δ(θ) = Δ cos[2(θ − γ)] with γ, the misorientation angle between the surface normal and the crystalline axis along which the order parameter reaches maximum.

Results and Discussion

Transport Characteristics of N-I-n-S Junctions

We start with a simplest case of superconducting heterostructures with normal interlayers. To calculate the current I-versus-voltage V curves for N-I-n-S junctions with an s-wave superconductor, we use the Landauer-Büttiker formalism applied to superconductor-based structures [21]where RN is the normal-state resistance, θ is the injection angle between an electron wave vector and the x-axis, the reference potential of a superconducting side is put to zero, f(ε) is the Fermi-Dirac distribution function, and D(ε, θ) is the electron penetration probability. It is easier to calculate the latter quantity in an N part of the structure whereand Ree(ε, θ) and Reh(ε, θ) are angle-dependent probability amplitudes for an electron entering the N-electrode to be scattered back as an electron or as a hole, respectively. To obtain them, we interpret the charge transmission across a heterostructure as a sequence of an infinite number of interface scattering events including Andreev electron-hole and vice versa transformations at N/S boundaries [22, 23]. Then, we getwhere is the complex-valued phase shift acquired during an electron (hole) path from one edge of the n interlayer to the other. Equation (6) provides an insight into an effect of the non-superconducting interlayer on transport characteristics of N-I-n-S junctions. Note that formally, it is reduced to multiplication of the standard formula (3) for the Andreev scattering amplitude by a factor exp(iφe + iφh). Such procedure permits to use all previously developed expressions for N-I-S and S-I-S junctions by only modifying related Andreev scattering coefficients.

Figure 1 exhibits differential-conductance curves for a planar three-dimensional N-I-n-S structure with an s-wave superconductor and a normal n interlayer where the probability to transfer barrier I is much less than unity (the tunneling regime). The impact of the n layer thickness dn and that of the mean free path ln is controlled by the parameters α = 2dnΔ/(ℏvF) and β = dn/ln, respectively. In the limit dn → 0, we get a well-known coherence peak in the differential-conductance curve at V = V (Fig. 1, solid curve) that reflects the presence of a surface bound state in the n interlayer of a vanishing thickness, i.e., at the IS interface. Its asymmetry arises due to the non-analytical behavior of the scattering amplitudes reh(he)(ε) at ε = eV. With increasing dn, the bound state formed due to the interference of electron and hole waves is shifted to lower voltages and a pronounced dip appears slightly below V = V (Fig. 1, dashed curve). For larger dn, we can see the appearance of an additional hump structure (Fig. 1, dotted curve) above V = V which reminds about the presence of a non-analytical square-root dependence of Andreev-reflection amplitudes (3). All three main elements of the peak-dip-hump structure are shown by arrows in the inset in Fig. 1 that illustrates the effect of the electron mean free path on the conductance spectra of the tunneling N-I-n-S(s-wave) structure.

Fig. 1

Main panel: Differential conductance-versus-voltage characteristics for a planar three-dimensional N-I-n-S(s-wave) junction with various thicknesses of the normal n interlayer in the tunneling regime (D I << 1); parameters α = 2d n Δ /(ℏv F) and β = d n/l n. Inset: Effect of the electron mean free path on the conductance spectra of a tunneling N-I-n-S(s-wave) four-layered device; the three arrows show the main elements of the peak-dip-hump structure discussed in the text

Figure 2 shows that the presence of a peak-dip-hump structure does not strongly depend on the barrier transparency. It is more pronounced in the tunneling regime (DI << 1) but is well reproduced even for a high-transparency device with DI = 0.5. As is expected, the temperature and charge scattering processes significantly destroy the anomalous features (see the insets in Figs. 1 and 2).

Fig. 2

Main panel: Effect of the barrier transparency on the peak-dip-hump structure in conductance spectra of a planar three-dimensional N-I-n-S(s-wave) junction; parameters α = 2d n Δ /(ℏv F) and β = d n/l n. Inset: Temperature effect on differential conductance-versus-voltage characteristics for a planar three-dimensional N-I-n-S structure with a high-transparency tunneling barrier (D I = 0.5)

Before turning to novel high-temperature superconductors, let us discuss old experiments on four-layered N-I-n-S structures having the n interlayer for sure. Here, we refer to three papers with n = Cu and S = Pb [24, 25] as well as with n = Al and S = NbZr alloy [26]. In the first paper [24], the authors studied Cu-Pb sandwiches with a 500-nm-thick lead film and copper interlayers of different thicknesses. At 0.06 K, they observed a sharp dip at about the lead energy gap for dCu between 40 and 120 nm, and for higher thicknesses, the amplitude of the singularity was vanishing, probably due to the impurity scattering inside the Cu interlayer [24]. These results were confirmed in ref. [25] for dCu = 35 and 75 nm and dPb = 700 nm. The authors [24, 25] tried to interpret their data within the McMillan tunneling model of the superconducting proximity effect [9] but experienced large discrepancies, in particular in the dip region. To explain them, the authors of ref. [24] supposed that the dip is the result of quasiparticle interference in the copper interlayer. Our simulations, Figs. 1 and 2, are based just on this assumption and qualitatively agree with the data shown in [24-26]. Measurements of In-I-Al-NbZr alloy junctions [26] also revealed the peak-dip structure that was interpreted as a bound state within the Al interlayer. Notice that additional peaks and dips in the differential conductance of tunnel junctions in ref. [26] can arise due to resonant tunneling processes across localized states in the interlayer [27].

Further confirmation of the bound-state assertion for the conductance peak came from an experiment [28] with a ferromagnetic (F) nanometer-thick film which replaced the normal interlayer in the four-layered junction. In this case, the charge reflected at the F-S interface is created in the electron density of states with a spin opposite to that of the incident quasiparticle and Eq. (1) is modified as follows

Since, for example, the difference in wave vectors in Ni is of the order of several nm−1, effects very similar to those predicted for N-I-n-S junctions with dn of the order of tens of nanometers can be now observed in N-I-F-S ones with a few nanometer-thick Ni interlayer. Indeed, an anomalous “double-peak” conductance feature revealed in ref. [28] at V near ΔNb/e, in our opinion, is a peak-dip-hump structure with a maximum strongly smeared by local inhomogeneities within the F interlayer.

We should note that the calculations shown in Fig. 1 relate to a conventional s-wave superconductor while the peak-dip-hump structure was revealed in high-Tc cuprates with d-wave pairing. Thus, we should now discuss the latter case where the order parameter is an angle function Δ(θ) = Δ cos[2(θ − γ)]. If the misorientation angle γ = 0, we have the same phase shift (1) as in the conventional s-wave superconductor with the only exception, the angle dependence of the gap parameter . The conductance spectrum for such an orientation is shown in Fig. 3. As in the case of an s-wave superconductor, it exhibits a peak-dip-hump structure for sufficiently thick normal interlayers. On the contrary, if γ = π/4, we get the following expression for the angle-dependent phase shift

Fig. 3

Main panel: Differential conductance-versus-voltage characteristics for a planar three-dimensional N-I-n-S(d-wave) structure with different thicknesses of the normal n interlayer in the tunneling regime (D I << 1); the angle γ = 0; parameters α = 2d n Δ /(ℏv F) and β = d n/l n. Inset: Effect of the electron mean free path on the conductance spectra of the tunneling N-I-n-S(d-wave) structure

It is important that for such film orientation, independently on the n layer thickness dn, the phase shift Δφ vanishes at ε = 0. Hence, in this case, we should always observe a zero-bias peak which is known to be a d-wave superconductivity indicator [29]. For comparatively thick interlayers, in addition to the zero-bias anomaly, we register, at first, a small peak at voltage bias VΔ and with the further increase of dn, the emergence of a peak-dip-hump structure, see Fig. 4.

Fig. 4

Main panel: Differential conductance-versus-voltage characteristics for a planar three-dimensional N-I-n-S(d-wave) structure with different thicknesses of the normal n interlayer in the tunneling regime (D I << 1); the angle γ = 45°; parameters α = 2d n Δ /(ℏv F) and β = d n/l n Inset: Effect of the electron mean free path on the conductance spectra of a tunneling N-I-n-S(d-wave) structure

Since high-Tc samples are often composites of randomly oriented superconducting grains, we predict that, in the presence of an additional n interlayer, their conductance spectra should exhibit a pronounced peak-dip-hump structure which is present in all orientations and a weak zero-bias anomaly specific for a certain tunneling direction.

Transport Characteristics of S-n-I-n-S Junctions

Physical mechanisms of the charge transmission across heterostructures with two superconducting electrodes are more complicated. In this section, we limit ourselves to s-wave pairing taking into account the fact that, after averaging over the realizations of the disorder, the order parameter in composites of randomly oriented superconducting grains with d-wave paring has global s-wave symmetry [30].

Let us start with the simplest tunneling transmission processes in S-n-I-n-S junctions. The presence of bound states at the electrode surfaces gives rise to unusual features in transport characteristics when the energies of the discrete levels align by applying the external bias. When properly aligned, we get the tunneling current peak and a negative differential resistance just above the resonant bias. Our calculations were done using a standard formula for tunneling current in a symmetric S-n-I-n-S configuration [1] where the tunneling density of states should be replaced by the electron penetration probability (5). The results are shown in Fig. 5. Calculated current-voltage characteristic for a three-dimensional S-n-I-n-S junction well reproduces the “knee” feature, see the “Background” section. Moreover, the first derivative shown in the inset in Fig. 5 demonstrates a peak, a region of the negative differential conductance which can reveal itself as a dip, and a hump at V = 2Δ/e. Note that just the position of the hump but not that of the peak determines the doubled value of the energy gap Δ in the studied s-wave superconductor.

Fig. 5

Main panel: Current-voltage characteristics for a planar three-dimensional S-n-I-n-S junction with identical s-wave superconductors and various thicknesses of the normal n interlayer in the tunneling regime (D I << 1). Inset: Conductance spectra of the S-n-I-n-S device, α = 1

Now, let us discuss a high-transparency junction formed by two superconductors with degraded surfaces. When the length of the transition region between the two superconducting electrodes is less than elastic and inelastic lengths, we can again describe the transport across the device in terms of Andreev-reflection amplitudes and barrier scattering characteristics. But in this case, comparing to the tunneling regime, a new specific feature largely complicates the calculations. Let us look at an electron-like quasiparticle injected, for example, from the left electrode. Its energy after transferring the barrier is increased by eV. The Andreev-reflected hole-like excitation is moving in the opposite direction to that of the injected electron, and because of the opposite charge sign, its energy also increases by eV after the transmission across the barrier. These scattering events will continue back and forth in the normal interspace between the superconductors, and, as a result, each round trip of an electron-like quasiparticle will increase its energy by the value of 2eV. Finally, from each side of the barrier, we get an infinite set of scattering states with different energies shifted by 2eV. It results in recurrence relations for amplitudes of electron-like and hole-like wave functions which have been solved numerically. Our numerical simulations at zero temperature repeated in general terms similar calculations developed earlier [20, 31–33] with the only exception of the presence of an additional phase shift originated from charge passages across the n interlayer. Due to the standard theory of multiple Andreev reflections, the peaks positions in the conductance spectra are determined by the formula V = 2Δ/n where n is an integer. Account of the additional phase shifts in the degraded near-surface regions replaces the peaks by dips (see Fig. 6). It is just what was found in ref. [14] and some previous publications for related Nb-based samples. The difference between Al-based junctions that follow the conventional theory [3] and Nb-based structures which contradict it arises from the existence of a proximity layer, either normal, or superconducting with reduced critical temperature, at the surface of a superconducting Nb film whereas the Al layer is usually not spoiled.

Fig. 6

Differential conductance-versus-voltage spectra for a planar S-n-I-n-S junction with identical s-wave superconductors and high-transparency transition region (D I = 0.7) without n interlayer (α = 0) and for a finite thickness of the normal n interlayer. Shaded regions show positions of the 2Δ/n features (peaks in the first case and dips in the second case) in the conductance spectra

Conclusions

We have presented a scattering-like approach for studying Andreev bound states in a normal interlayer within ballistic superconducting heterostructures. Calculated transport characteristics of such devices exhibit a rich diversity of anomalous (from the viewpoint of the standard theory) transport characteristics, in particular shift of differential-conductance maxima at gap voltages to lower positions and appearance of well-defined dips instead expected coherence peaks. We have compared our results with related experimental data and explained unexpected features observed in transport characteristics of heterostructures with transition metals and novel superconductors as an effect of the additional non-superconducting interlayer at their surfaces. Note that our results are qualitatively similar to those obtained by other authors in the framework of physical assumptions that significantly differ from ours [8, 34, 35].