Anomalous Fano Resonance in Double Quantum Dot System Coupled to Superconductor

We analyze the influence of a local pairing on the quantum interference in nanoscopic systems. As a model system we choose the double quantum dot coupled to one metallic and one superconducting electrode in the T-shape geometry. The analysis is particularly valuable for systems containing coupled objects with considerably different broadening of energy levels. In such systems, the scattering of itinerant electrons on a discrete (or narrow) energy level gives rise to the Fano-type interference. Systems with induced superconducting order, along well understood Fano resonances, exhibit also another features on the opposite side of the Fermi level. The lineshape of these resonances differs significantly from their reflection on the opposite side of the Fermi level, and their origin was not fully understood. Here, considering the spin-polarized tunneling model, we explain a microscopic mechanism of a formation of these resonances and discuss the nature of their uncommon lineshapes. We show that the anomalous Fano profiles originate solely from the pairing of nonscattered electrons with scattered ones. We investigate also the interplay of each type of resonances with the Kondo physics and discuss the resonant features in differential conductivity.

Impurities or nanoobjects like quantum dots (QDs) hybridized to superconductors (SC) adopt some SC properties via proximity effects. As a consequence, the ground state of a QD is represented by either single particle state |↑〉, |↓〉, or a superposition of empty and doubly occupied states u v 0 d d + ↑↓ [1][2][3] . The fingerprints of this local pairing can be observed in the Andreev spectroscopy as two quasiparticle peaks 1,4,5 . Currently, dynamic development in fabrication of complex nanodevices on the top of SC substrate allows to construct SC-based systems built of multiple QD's 6 , quantum rings 7,8 , monoatomic chains 9 , gate-controlled carbon nanotubes (CNT) 10,11 , multiwall CNT quantum dots 12 , modified Aharonov-Bohm rings with a QD embedded within one of the ring's arms 13 , SQUID interferometers with gate-controlled CNT quantum dots 14 or quantum dots connected to Rashba chains 15 . In such systems, the various paths for electron propagation give rise to quantum interference effects. Therefore, deep understanding of mutual relations between the proximity induced pairing and the quantum interference is highly demanded. A classic model to analyze such relations consists of a QD (QD 1 ) coupled directly to (i) one metallic electrode and (ii) one superconducting electrode as well as side-coupled to second QD (QD 2 ) [a schema of the system is shown in Fig. 1(a)]. In such system, the main charge transport between electrodes leads directly through the central quantum dot (i.e., QD 1 ), Fig. 1(b). Additional path includes the electron hopping between the central dot (QD 1 ) and the side dot (QD 2 ). Different paths for electron transport overlap giving rise to quantum interference effects. As the interfacial quantum dot (i.e., QD 1 ) is connected to superconducting reservoir, scattering on a side level is accompanied by the local pairing. In metal-hybrid structures interference patterns can be observed in the spectral function and transport characteristics as asymmetric Fano features emerging at energies equal to the energy level of the side dot(s) [16][17][18][19][20][21][22][23] . In the presence of local pairing, two resonant structures emerge simultaneously on both sides of the Fermi level [24][25][26] . A shape of the feature located at the energy of the side dot resembles the ordinary Fano resonance. However, a structure on the opposite side of the Fermi level seems to diverge from the ordinary Fano profile 24,25 . Although a particle-hole mixing of states rationalizes an appearance of two resonances instead of just one, astonishing difference in their profiles is intriguing. One could even argue whether shape of additional resonance should be referred to as the Fano-like.
The Fano-like profiles have been reported in numerous works in various fields of physics including atomic [27][28][29] , molecular 30 physics, photonics 31,32 , plasmonics 33,34 , electron-phonon interaction [35][36][37][38][39][40] , microwave physics 41,42 , metamaterials and nonlinear optics 34,43,44 , ultra cold gases 45 , or nuclear physics [46][47][48][49] . The Fano resonances turn out to be particularly relevant also for nanoscale physics. In various systems, in which nanoobjects with different broadenings of energy levels are tunnel-coupled, similar resonances appear on the background of the Breit-Wigner resonance (or, equivalently, the Lorentz distribution). Such Fano-like resonances were predicted and observed, e.g., in double [50][51][52][53] and triple 26,54 quantum dots systems in various configurations 50 . The asymmetric resonances were also predicted in "bridge" realization, where two electrodes were tunnel coupled to a single QD and, additionally, to each other directly 21,55 . In such a realization, the Fano effect arises as a result of interference of waves traveling directly between electrodes with a localized state. It was predicted that the Fano resonances appearance in a similar configuration can enhance the effectiveness of a Cooper spliter device 56 .
Recently, it has been noticed that the Fano resonances can be useful in indicating the existence of the Majorana bound states [57][58][59] . In systems, in which the quantum dot is weakly connected to the Rashba chain, the scattering on the Majorana zero mode (MZM) suppresses the local density of states (LDOS) of QD only by one half 59,60 . This is because electrons scattered by the Majorana quasiparticle change their phase only by the fraction of π, while in the ordinary Fano effect, we observe π − 0 phase shift features. Taking this into account, one can distinguish scattering of electrons on the MZM from scattering on the topologically trivial zero energy states. Rich interplay of the Fano resonances with strong correlations effects was analyzed by number of authors both in metallic 21,61-63 and superconducting environment 24,64-69 . Among others, it was found that suppression of the Kondo state by its coexistence with the Fano antiresonance reveals a novel Fano-Kondo resonance 66,67 . A pedagogical review of the Fano resonances in nanoscale physics was done by A. E. Miroshmishenko in ref. 70 . Interplay of the Fano resonance itself with local pairing was less widely explored. P. Orellana and coworkers 71,72 analyzed configuration with one quantum dot placed between two metallic electrodes and side dot coupled only to the SC electrode. In such realization, scattering on narrow quasiparticle states gives rise to two Fano-like features on background of single particle broad level.
In this work, we present the analysis of the local pairing for electrons scattered on the side structure. In the considered model, the SC electrode is connected directly to the interfacial dot, thus quasiparticle states are considered as broad continuum while scattering occur on the dot decoupled from the SC environment (cf. Fig. 1). We discuss origin of appearing resonances and reveal the microscopic mechanism of their formation. We analyze the shape of obtained resonances by comparing them with the Fano profiles and calculate characteristic Fano factors such as asymmetry parameter. We also discuss the interplay of each resonant feature with the Kondo resonance and inspect the appearance of resonant features in differential conductivity.

Formulation of the Problem
A heterojunction depicted in Fig. 1(a) can be modeled by the Anderson impurity Hamiltonian in the following formˆˆˆˆˆ∑ refers to s-wave superconducting electrode. Electron energies k ξ β are measured with respect to chemical potentials µ β β = N S ( , ). Two quantum dots connected with spin-dependent interdot hoppings t σ σ = ↑ ↓ ( , ) are represented by the following terms: Figure 1. (a) The schematic illustration of the analyzed system. It consists of two quantum dots (QD 1 and QD 2 ). QD 1 (interfacial one) is directly coupled with superconducting (SC) and normal metal (N) electrodes with spinindependent couplings Γ S and Γ N , respectively. Coupling t σ between QD 1 and QD 2 (side one) is spin-dependent [cf., Eq. (1)]. (b) The schematic illustration of the scattering processes occurring in the strongly spin-polarized tunneling model (t 0 ≠ ↑ and t 0 = ↓ ). Orange (green) arrows indicate the propagation paths for σ = ↑ (σ =↓, respectively) electrons. Only spin-↑ electrons can directly scatter on side dot QD 2 (vertical orange arrows), whereas spin-↓ electrons are paired with them (a blue spring represents the local pairing S Γ ) and scattered indirectly.
where i ε is the energy level of i-th quantum dot, U i stands for intra-dot Coulomb interactions = i ( 1, 2), and σ denotes the spin opposite to σ (e.g., ↓ = ↑). The hybridization of the interfacial (i.e., i 1 = ) quantum dot to the external reservoirs N S ( , ) β = is given by It is useful to introduce the wide band limit constant coupling strength between the interfacial dot and both reservoirs: . In the deep superconducting atomic limit ∆ Γ ( ) S  , the influence of superconducting electrode on interfacial quantum dot QD 1 is reduced to the induced local pairing. A problem of the "proximized" quantum dot was widely explored by many authors 2,73-77 including ourselfs 3,[78][79][80] . In such conditions, the Hamiltonian of the interfacial dot coupled to the SC reservoir (i.e., To gain a clear picture of the interplay between interference effects and the local pairing we will mostly focus on noncorrelated regime, i.e., U 1 = U 2 = 0 (excluding the section, where correlations are studied explicitly). Information on spectral properties and Andreev transmittance is encoded in particular Green's functions In the present work, we assume spin-dependent interdot hopping t σ , therefore, Green's functions for each spin component are not identical and for each index σ = ↑ ↓ , these functions need to be calculated separately. In the equilibrium conditions the equation of motion technique 81 yields the following expression for a Fourier transform of the retarded Green function matrix for the interfacial quantum dot (QD 1 ): pectral function (local density of states) ρ ω σ ( ) 1 of QD 1 for each spin component σ is given by standard formula 11 . The position of the Fermi level is ω µ = S , which is located in the middle of the superconducting energy gap. In the following, for a sake of simplicity, we also take that 0 Fano-like resonances in nanoscopic systems. If low dimensional structures with discrete energy spectrum (such as, e.g., quantum dots) are coupled to reservoirs characterized by continuum of states, the energy levels of nanoobjects are broadened to form the Breit-Wigner (Lorentz) distribution with half-width controlled by QD-bath coupling strength (i.e., Γ N ). Consequently, if subparts of a device are coupled to environment with different coupling strengths, broadening of particular energy levels differ significantly. Quantum interference of electron waves resonantly transmitted through narrow (quasidiscrete levels) and those transmitted through broad levels give rise to asymmetric Fano-like profiles observed in density of states and differential conductivity. For electrons whose energy tends to resonant energy from one side (e.g., res ω ω → + ), scattering does not change their phase, while for electrons reaching the resonant level from the other side ω ω → − ( ) res phase is shifted by a factor of π 67,70 . Therefore, in the Fano-like profiles a constructive enhancement (i.e., the 0 phase shift) is accompanied by an antiresonant deep (i.e., the π phase shift).
In the original work 82 and later review 83 , U. Fano described the asymmetric lineshapes in the absorption spectra of noble gases observed previously by H. Beutler 84 . He introduced the formula for such profiles on a flat background. This function can be expressed as where ω res and R Γ stand, respectively, for the energy of the resonant level and the effective broadening of the resonant level, and q represents a phenomenological factor called an asymmetry parameter.
Because it will be very useful in the further discussion included in the next section (Profile analysis), it is worth noticing a few properties of this function. Function (6) has one local minimum and one maximum. Arguments of the extremes are located at ω ω , respectively. Due to the ideal anti-resonance, the minimum value of this function is equal to 0, while the maximum depends only on the asymmetry parameter and it is equal to 1 + q 2 . Away from the resonant energy (i.e., for  ω ω res or ω ω  res ) the function reaches value equals to 1. For asymmetry parameter close to unity and small broadening of the resonant energy level ( ) R Γ the close proximity of the minimum and the maximum forms a well-pronounced asymmetric profile. Note that, for the asymmetry parameter approaching to zero, the Fano function corresponds to a symmetric deep, whereas for the parameter going to the infinity, the resonant feature resembles the Lorentz function.
In the nanoscopic systems, the broadening of the resonant level is dependent on the coupling to the continuum of states. Namely, in the case of the double quantum dot, it is proportional to the square of the interdot coupling. For the strong interdot coupling, the broadening of the side level becomes comparable with the broadening of the central dot (cf., e.g., refs. 85,86 and references therein). In such a case, the interferometric structures evolve into molecular states. Nevertheless, this issue is out of the scope of the present work, where we focus on the Fano-like features.
Fano-like resonances in a presence of superconducting electrode (symmetric case: t ↑ = t ↓ ≠ 0). In a hybrid system, where single QD is coupled to SC reservoir, due to proximity effects single particle QD's level evolves into two quasiparticle peaks representing so-called Andreev bound states (ABS). These states in noncorrelated regime emerge at ω = ±E 1 , where quasiparticle energy levels are represented by E S 1 1 2 2 ε = +Γ , and they are weighted by a corresponding BCS coefficients . In a particular case of 0 1 ε = , the Andreev states are symmetric Lorentzians separated by Γ S . Density of states of the QD in such conditions (for ε = 0 1 ) can be expressed as If one QD (QD 1 ) is coupled to both metallic and superconducting electrodes and, additionally, side-coupled to the second quantum dot (QD 2 ) with the spin-independent coupling (i.e., the system shown in Fig 24,25 . First one appears for energies close to energy level of the side dot ω ε ≈ ( ) 2 . Asymmetric lineshape of this resonance resembles the characteristic Fano-like shape. Second feature emerges on opposite side of the Fermi level (ω ≈ −ε 2 ). This resonance, however, differs significantly from the former one. First notable observation is a sharp spike apparent on one side of the resonance near ω ≈ −ε 2 . In terms of the Fano function such imbalance emerges for very large asymmetry parameters q. Second peculiar observation is that local minima near this resonance (in particular this for ω =−ε 2 ) have a finite value while the ordinary Fano function vanishes for ω ω Strongly asymmetric spin-polarized tunneling model (t ↑ ≠ 0 and t ↓ = 0). To understand the origin of both resonances, we consider a ("toy") model in which we completely suppress the interdot coupling for one spin component (e.g., for spin σ = ↓). In the metal-hybrid structures two spin channels are independent. In a presence of superconducting electrode spin-↑ and spin-↓ electrons are bound into the local pairs. Thus, every physical process that involves one spin component affects also the other one. By applying the strongly asymmetric condition (i.e., t 0 = ↓ ), we effectively decompose the effect of the direct interference, which occurs only for electrons coupled to the side dot (i.e., these with σ = ↑) from the effect of bounding them into the local pairs [cf. Fig. 1(b)]. The latter effect can be observed in the spectral function of electrons decoupled from the side dot (σ = ↓). As spin-↓ electrons are not directly scattered the resonant characteristics appearing in their spectral function originate solely from pairing with scattered electrons.
The spectral function of QD 1 for both directions of electron spin are shown in Fig. 2(b). In such conditions, for directly scattered electrons (σ = ↑), we obtain only the resonant feature near 2 ω ε ≈ , while the second resonance disappears. Counter-wise, for opposite spin electrons (with σ = ↓) only the feature located near ω ≈ −ε 2 remains. It is worth noting that the shape of the resonant features remains unchanged, i.e., the shape of the resonance near ω ε ≈ 2 (ω ≈ −ε 2 ) for the symmetric case (t ↑ = t ↓ ≠ 0) is identical as the resonant feature that remains in ( ) , respectively) for the perfectly-polarized interdot coupling (t ↑ ≠ t ↓ = 0), cf. Fig. 2(a,b). Thus, this indicate (solid blue line) and ρ ω ↓ ( ) 1 (dashed red line) of QD 1 for the spin-dependent coupling (the strongly asymmetric condition) with QD 2 . The model parameters are: t 1 for the spin-dependent coupling (the strongly asymmetric condition) with QD 2 . In the inset the region near ω ≈ 0 is shown. The model parameters Scientific RepoRtS | (2020) 10:2881 | https://doi.org/10.1038/s41598-020-59498-y www.nature.com/scientificreports www.nature.com/scientificreports/ that the resonance near ω ε ≈ 2 originates purely from the direct scattering of electrons on a side level, while the resonant characteristic near ω ≈ −ε 2 is merely a response to the pairing of a given electron with its scattered partner.
In other words, one spin component is directly scattered on side structure while the other one "feels" the scattering only by bounding into a local pair with directly scattered one. Such the conclusion is also visible if we compare the analytic formula for single particle Green's functions for each spin component at QD 1 . These functions for perfectly spin-polarized tunneling (t are parts responsible for the scattering [cf. also Eq. (5) Therefore, for the nonpolarized case, we observe both the ordinary Fano feature near 2 ω ε ≈ and the "anomalous" Fano resonance near ω ≈ −ε 2 for electrons with spin-↑ as well as for electrons with spin-↓. For t 0 = ↓ , the above equations for σ = ↑ ↓ , reduces to Eq. (8) and each resonant feature occurs in different spin channel. The local density of states at QD 1 for each spin component is given by the imaginary part of an adequate Green In the case of t 0 = From the above equations, it is difficult to see if the shape around the resonant energies ω ≈ ±ε 2 can be described as the Fano-like shape. Moreover, for ε = 0 2 the sharp resonant peak at ω ≈ −ε 2 evolves into a symmetric Lorenzian [cf. Fig. 2(c)]. In the next section, we present analysis of resonant features in both spin channels for = ↓ t 0. For a sake of simplicity we focus on the case of ε = 0 1 , which is studied further in this work.
Profiles analysis (asymmetric case: t ↓ = 0). The Fano resonances were successfully used as a probe for electron phase coherence in quantum dots 87 . It was shown that dephasing time can be determined from the asymmetry parameter (q) of measured profiles. This issue was particularly relevant to take the meaning of Fano profiles appearing in single-electron transistor 88 ). The influence of such dephasing on Fano resonances was also analyzed by one of us 25 . Thus, a proper evaluation of the asymmetry parameter for a given profile turn out to be a relevant issue. In the case of resonances that appear on non-flat backgrounds, straight fitting of the regular Fano function may produce highly inaccurate values. The problem becomes even more complicated if a given shape deviates from the regular Fano profile. In here analyzed system the resonances near ω ≈ −ε 2 exhibit features that do not match the ordinary Fano shape.
In this section, we will analyze the obtained resonant lines to compare them with the Fano profiles and indicate to what extent a given profile can be approximated by the Fano function [cf. Eq. (6)]. We develop a feasible procedure of fitting the Fano parameters to the assumed form. In the case of resonances that deviate from the ordinary Fano shape, we take into account and estimate the correction factor 0 φ . In general, the Fano-like resonances can be represented as a function www.nature.com/scientificreports www.nature.com/scientificreports/ resonance appears on a flat singular background). In our case, the resonant features emerge on the background of quasiparticle Andreev states described by Eq. (7). Therefore, we assume that one can approximate the density of states ( ) 1 ρ ω σ by a product of the ordinary Fano curve and a background composed of the Andreev states, i.e., by where ω S( ) is given by (7). If so, it should be possible to find the relations between the model parameters and the parameters used in the function F ( ) ω σ : q σ , R, Γ σ , and res , ω σ in such way that the constructed function ω ω σ F S ( ) ( ) will reproduce the density of states ρ ω σ ( ) 1 with high accuracy.
The resonance for directly scattered electrons (near ω ≈ ε 2 ). We will start with an analysis of the LDOS ( ) The sign of asymmetry parameter ↑ q is governed by a position of ε 2 (i.e., for ε > 0 2 one gets q 0 < ↑ and for ε < 0 2 one has > ↑ q 0). The other two parameters ω ↑ res , and Γ ↑ R, can be found by comparison of positions of local minimum and maximum of ρ ω ↑ ( ) 1 and ω ↑ F ( ). Arguments ω for which the Fano function takes the minimum is given by , while for the maximum ω ω = +Γ . This yields where q ↑ is an asymmetry parameter estimated previously from Eq. (10). In Fig. 3(a) we examine the convergence of the obtained function with the exact prototype. One can note that for spin-↑ electrons the product function with Γ ↑ ↑ q , R, and ω ↑ res , estimated by the above procedure reproduces original ( ) 1 ρ ω ↑ with very high accuracy. Asymmetry parameter q ↑ for this fit of the Fano resonant feature in ( ) 1 ρ ω ↑ as a function of ε 2 is shown in Fig. 4(b) as dotted green line.

The resonance for indirectly scattered electrons (near ω ≈ −ε 2 ). Approximation with the regular
Fano function. Situation is more complicated for opposite spin electrons, i.e., for spin-↓ electrons. The minimum value of the exact function is finite, while minimal value for the ordinary Fano resonance is equal to 0. This means that product of the ordinary Fano resonant curve and arbitrary background [in particular, also function ] will never reproduce the original function ( ) 1 ρ ω ↓ for spin-↓ electrons accurately. However, we can still assume that spectral function can be approximated by the Fano function with small correction 0 φ . Let us assume that the exact function can be approximated by a product of ω S( ) and where 0 φ represents a (small) deviation (correction) dependent only on model parameters (and it is independent of ω) and ω ↓ F ( ) has a form of Eq. (6) with q ↓ , R, Γ ↓ , and ω ↓ res , parameters. One should note that if correction param- www.nature.com/scientificreports www.nature.com/scientificreports/ eter 0 φ will be small enough (say less than 1) then one can state that the resonant shape can be approximated by the Fano function with high accuracy. Thus, auxiliary parameter φ 0 can be considered as a measure to what extent one can fit the original Fano shape into a given resonant feature. Therefore, a high value of correction 0 φ indicate that the Fano function may not be adequate for fitting to the exact function. Introduction of 0 φ slightly rearranges the background of the ordinary Fano function, i.e., away from the resonant energy F ( ) ω ↓ ⁎ reaches 1 0 φ + instead of just 1. This would cause problems in achieving acceptable convergence of assumed form and the exact function away from the resonant feature. To neutralize this obstacle we need to normalize our assumption. Therefore, instead of ω ω . Extreme arguments ω − and ω + can be calculated numerically (in experimental realizations these values can be directly read from the data). Therefore, we should add just one step in our procedure. First, we get 0 φ comparing local minima of ( ) 1 ρ ω ↓ and assumed form (13) Fig. 3(b,c)]. However, good convergence with a small correction is achieved only for resonant energies close to the center of Andreev states ( ) S 2 ε ≈ ±Γ (see solid blue line in Fig. 4(a)). For resonant energies close to the Fermi level, spectral function ( ) 1 ρ ω ↓ is represented by almost a symmetric peak [cf. Fig. 2(c)]. For the original Fano curve (6) such a situation is met when the asymmetry parameter approach infinity q ( ) → ∞ . On the other hand, very large asymmetry parameter implies huge values for the maximum of the Fano function F ( ) ω ↓ . For the product to remain finite, the correction must grow together with the asymmetry parameter. Thus, in assumed , "correction" φ 0 (as well as asymmetry parameter q ↓ ) become enormously large (cf. Fig. 4(a); they tend to infinity if 0 2 ε → ). A correction value is also inadequate if the resonant energy is located far outside the Andreev states (i.e., . These make the statement about such resonances as the Fano-like one somewhat exaggerated. Asymmetry parameter q ↓ of F ( ) ω ↓ for this fit of the resonant feature in ( ) 1 ρ ω ↓ as a function of ε 2 is shown in Fig. 4(b) as solid blue line.

Approximation with the inverse of the Fano function.
To underline the fact that the ordinary Fano function is not the best way to approximate the resonant feature near ω ≈ −ε 2 (at least at same range of the model parameters), we will try to fit another asymmetric function with a well-defined asymmetry parameter and compare the result with fitting of the ordinary Fano function. The problem of an inadequate correction for a wide spectrum of ε 2 can be reduced by a slight rearrangement of the assumption. The resonant feature in ρ ω ↓ ( ) 1 is composed of a finite deep accompanied by an over-sized peak. We noticed that if we add a small parameter to the regular Fano function, and then we take the inverse of that structure, the resulting function should have similar features. Thus, we www.nature.com/scientificreports www.nature.com/scientificreports/ propose to approximate the original function by a product of the Andreev states [i.e., S( ) ω ] and an inversion of the Fano function. Taking into account the normalization as previously, our assumption should be in the following form has a form of Eq. (6) with ↓  q , R,  Γ ↓ , and ω ↓  res , . Using the procedure similar to the previous one, we determine the corresponding parameters and compared the obtained result with the exact function ρ ω ↓ ( ) 1 . Here, maximum (minimum) of  F ( ) ω ↓ at ω + (ω − , respectively) corresponds to minimum (maximum) of ( ) 1 ρ ω ↓ . Thus, using the properties of the Fano curve ω ↓ F ( )  in the similar manner as previously, one gets that . We found that the new assumption reproduces the original function as accurately as the previous one [cf. Fig. 3(b,c)]. The advantage of such a fit is that correction φ 0  for such assumption is considerably smaller for a wide range of 2 ε excluding S 2 ε ≈ ±Γ , where it expands to infinity (Fig. 4(a)). If the resonant energy is close to Γ S , the spectral function can be approximated by the regular Fano shape with a small correction using Eq. (13) again, as described in previous section. The inverse of asymmetry parameter for this fit of the resonant feature in ρ ω ↓ ( ) 1 as a function of ε 2 is shown in Fig. 4(b) as dotted red line. Note, that in this case, ↓  q 1/ (rather than q  ↓ ) is a measure of the asymmetry comparable to the ordinary Fano asymmetry parameter (as for  → ↓ q 1/ 0 resonant feature is represented by a symmetric deep and for → ∞ ↓ q 1/  by the Lorentz distribution). One should note that, in contrast to the Fano function, which originates from a rigorous examination of transmission rates in noble gases 28 , function (14) is a hypothetical (semiempirical) function that can be fitted into "anomalous Fano" curves more accurately at wide range of parameters. Nevertheless, using such a function one can estimate the parameter  ↓ q , which is responsible for a measure of an asymmetry of the resonant feature (and thus, it is fragile for decoherence).
Concluding, the LDOS function for indirectly scattered electrons can be treated as a normalized product of the Andreev states and (i) Fano resonance or (ii) inverse Fano, both with a small correction. The first approach reproduces well the original function only for resonant energies close to Γ S , i.e., it describes a case when the resonant energy coincide with the Andreev states. The second approach works well for resonant energies much smaller and much larger than Γ S . To have a full insight into the behavior of LDOS function, it would be useful to combine these two approaches or determine the ratio / S 2 ε Γ and then use appropriate product.

Effects of correlations between electrons on the quantum dots. In nanoscopic systems, the
Coulomb repulsion between electrons often plays an important role, therefore, in this section, we briefly discuss the interplay of correlations effects with the analyzed features. In the model, QD 2 is not directly connected to any external reservoir and on-site interactions on QD 2 (i.e., U 2 ) lead only to appearance of an additional narrow state in the spectrum of QD 2 located at U  24 ). To account for the correlations on QD 1 , we adopt procedure used previously in ref. 24 . In the presence of correlations, the matrix of Green's functions can be represented by is approximated using the decoupling scheme, which approximates higher order Green's functions and reduces them to lower order once; details are given in ref. 24 . It yields: N , where Σ η=1,3 are given by k kN kN kN being the Fermi distribution at temperature T , 〈 〉 σ n denotes an average occupancy of QD 1 with spin-σ electrons (calculated self-consistently), and i /2 N 0 Σ = Γ . In this section, we will investigate the spin-dependent energy spectrum in the perfectly polarized case, i.e., = ↓ t 0. The symmetric spin interdot coupling case was described in ref. 24 . We inspect two cases: (i) the strongly proximized case, where the In the correlated case U ( 0) 1 ≠ each Andreev state splits into low and high energy branches separated by energy U 1 , cf., e.g., ref. 3 for a detailed discussion of this issue. In realistic systems, the Coulomb interactions are usually larger than the energy gap of the superconductor (i.e., ∆  U 1 ). Consequently, the high energy states coincide with the continuum outside the gap. Therefore, high energy branches of Andreev states are beyond considered energy scale and they do not introduce any meaningful physics to the scope of this work. The detailed analysis of the competition between the local pairing and correlations in the absence of the scattering is conducted in refs. 2,3 . In Fig. 5, obtained for strongly correlated regime U N 1 = Γ , we present the features of the scattering on the background of two low energy Andreev states in the case of spin-polarized tunneling [ Fig. 5(a)] and, for a comparison, in the nonpolarized case [ Fig. 5(b)]. One can note that resonant features, described earlier for noncorrelated case, remain qualitatively unchanged despite of strong correlations. Additionally, one can note a small step near ω ≈ 0. This can be described as underdeveloped Abrikosov-Suhl state also known as Kondo peak. These resonances appear as a result of screening of electron spin located on the quantum dot (or impurity) by opposite spin itinerant electrons from the metallic electrode.
To inspect the interplay between scattering features and Kondo state we will analyze the conditions more suitable for full development of the Kondo state. Namely, we put energy of QD 1 dot slightly below the Fermi level ( 15 ) N 1 ε = − . Γ , comparable hybridizations Γ = Γ ( ) S N and low temperature = . Γ k T 0 01 B N . In such conditions, two Andreev states overlap on each other forming a structure resembling Lorentzian. One should remember that this structure is still built of two quasiparticle states, which become well-separated if the hybridization to the SC electrode Γ ( ) S is considerably larger than broadening Γ N (cf. Fig. 7 Fig. 6(e) as well as Fig. 5(b). Consequently, the ordinary Fano feature originating from direct scattering and the anomalous Fano feature coincide if the energy of the side dot 2 ε is equal to 0. As it can be seen in Fig. 6(f), in such a case, both these resonances also coincide with the Kondo spike. We note that, for such parameters, the destructive interference plays a dominant role as zero energy state is strongly suppressed (forming a structure resembling the Fano-Kondo feature). Contribution of the anomalous Fano resonance in this case is visible as a small spike slightly below the Fermi level and slight enhancement of the Fano-Kondo feature just above the Fermi level [cf. maximal values in the inset of Fig. 6(f) and in Fig. 6(d) for the fully spin-polarized case]. One should note that a similar Fano-Kondo structure was predicted, e.g., for the double-quantum-dot system coupled to ferromagnetic electrodes [c.f., Fig. 4(f) in ref. 89 ]. However, in that work, the enhancement [and the features presented in the inset of Fig. 6(f)] associated with the anomalous Fano resonance are not present there due the absence of the local pairing in the system considered in ref. 89 . www.nature.com/scientificreports www.nature.com/scientificreports/ Resonant features in differential conductivity. The spectral function is not a directly measurable quantity. Therefore, the resonant features described in this paper can be investigated experimentally only by inspection of differential conductivity = G V dI dV ( ) / . For junctions with one metallic and one superconducting electrode low energy charge transport is supported solely by so-called Andreev reflections. In such processes, single electron of a given spin from the metallic lead is converted into a Cooper pair propagating in the superconductor with simultaneous reflection of a hole (with the opposite spin) back to the metal. This process, however, involves electrons of both spins equally. If an electron of a given spin and of energy ω is supposed to be converted into a Cooper pair propagating in the superconductor, it needs to "pick" additional electron of the opposite spin and of energy −ω. Particularly, for energies close to −ε 2 (where the "anomalous" Fano resonance emerges in spectral function of σ = ↓ electrons) electrons are paired with the opposite spin electrons of energy ε 2 , for which the ordinary Fano resonance emerges. Consequently, in the picture of the Andreev conductivity G dI dV / A A = , even for perfect spin-polarized case, resonant features near ω ε ± 2 become a mixture of the ordinary Fano and the "anomalous" Fano resonances as seen in Fig. 7(a). Indeed, the total Andreev current can be expressed by  To detect the resonant features separately one should rather inspect single particle transport. Therefore, we assume that interfacial quantum dot QD 1 is connected to an additional metallic electrode. Assuming that chemical potential of the superconducting electrode is tuned such, that no current is contributed on average from it (so-called floating lead), the only charge transport left is a single-particle current between two metallic electrodes. We calculate single-particle differential conductivity assuming that both normal electrodes are coupled to QD 1 . We assume that the energy unit is equal to a sum of the hybridizations of both metallic electrodes, i.e., Γ + Γ = Γ N N N 1 2 . Single-particle current I V ( ) σ is calculated using the Landauer formula 90 11 2 is a single particle transmittance. In Fig. 7, we also present differential conductivity G dI dV / = σ σ as a function of the applied voltage considering two cases (i) the toy model with the perfect spin-polarized scattering [t 0 Fig. 7(b)] and (ii) the more realistic case where both spin components can be tunneled between the dots [ Fig. 7(c)]. In the first case, the ordinary Fano resonance emerges in the conductivity of directly scattered electrons σ = ↑ [blue line in Fig. 7(b)], while the feature related to pairing with scattered electron is visible as sharp spike in conductance of the opposite spin electron near ω ≈ −ε 2 . For the nonpolarized case, electrons of each spin σ = ↑ ↓ , are both directly scattered and bound into a pair with the scattered electron of spin σ [ Fig. 7(c)]. Therefore, in the picture of conductivity, we can detect the regular Fano shape near 2 ω ε ≈ and the anomalous resonant feature near ω ≈ −ε 2 .

conclusions
In the present work, we inspect the energy spectrum of the double-quantum-dot system coupled to a superconducting reservoir in the T-shape geometry. In the analyzed system, combined effect of the electron scattering and the local pairing gives rise to two resonant features on the opposite sides of the Fermi level. Considering the perfectly spin-polarized interdot tunneling regime, we show that one of the resonances emerges as a result of the direct scattering. The other one emerges as a result of pairing of a given electron with a scattered electron of the opposite spin. Therefore, an existence of a pair of characteristic Fano and anomalous Fano resonances can be considered as a fingerprints for an occurrence of a bound state in the given system. We also obtained characteristics for differential conductivity and identified the features associated with the discussed resonances. These results of the work, derived for both strongly asymmetric = ↓ t ( 0) and symmetric = ≠ ↓ ↑ t t ( 0 ) cases, suggest that the regular Fano (near ω ε ≈ 2 ) and the anomalous Fano (near ω ≈ −ε 2 ) features could be detected in real nanoscopic systems. Such resonances can be observed in a variety of complex nano-systems (coupled to a superconductor) where the broadening of energy levels for each subsystem is considerably different. Although the spin-polarized model is hardly achievable experimentally without using the magnetic field, it allows to uncover the mechanism behind the formation of resonant features on both sides of the Fermi level appearing also in the realistic nonpolarized model (i.e., the symmetric case of = ↑ ↓ t t ) [24][25][26] . Note also that systems in such a configuration (i.e., T-shape one) can be investigated experimentally (cf. refs. [16][17][18][19][20][21][22][23]91,92 and references therein).
In this work, we showed that the resonant feature that originates from direct scattering can be described in terms of the Fano-like function with great details. Particularly, for the double-quantum-dot system coupled to a (normal) metal and a superconductor (Fig. 1), the spectral function of directly scattered electrons was approximated by product of the Fano line-shape and the Andreev states. A convergence of such approximation and the exact spectral function turned out to be very accurate. To achieve satisfactory convergence for the resonant feature on the opposite side of the Fermi level, one needed to impose additional correction φ ( ) 0 to the Fano function [Eq. (13)]. Using such an assumption, we managed to achieve a good convergence, but "correction" parameter φ 0 becomes enormously large, when the resonant level approaches the Fermi surface [e.g., for ε | | < . Γ 0 1 N 2 , parameter 0 φ becomes two order of magnitude higher than assumed energy unit (i.e., φ ≈ Γ 300 N 0 ), also the asymmetry parameter in such cases becomes as large as q 80 ≈ ↓ ]. Therefore, we proposed to approximate such resonances by the inversion of the Fano function [Eq. (14)] rather than the direct Fano one. Using this assumption, we achieved a high convergence with keeping correction  φ 0 small for a wide range of the model parameters. We also discussed the interplay of both Fano-like features with the Kondo resonance in the presence of correlations.