## Leakage of Majorana mode into correlated quantum dot nearby its singlet-doublet crossover

T. Zienkiewicz, 1 J. Bara´ nski, 1, ∗ G. G´ orski, 2 and T. Doma´ nski 3, †

1 Polish Air Force University, ul. Dywizjonu 303 no. 35, 08-521 Deblin, Poland 2 Faculty of Mathematics and Natural Sciences, University of Rzesz´w, 35-310 Rzesz´w, Poland o o 3 Institute of Physics, M. Curie-Sk/suppress lodowska University, 20-031 Lublin, Poland (Dated: August 2, 2019)

We study quasiparticle spectrum of the correlated quantum dot deposited on superconducting substrate which is side-coupled to the Rashba nanochain, hosting Majorana end modes. Ground state of an isolated quantum dot proximitized to superconducting reservoirs is represented either by the singly occupied site or BCS-type superposition of the empty and doubly occupied configurations. Quantum phase transition between these distinct ground states is spectroscopically manifested by the in-gap Andreev states which cross each other at the Fermi level. This qualitatively affects leakage of the Majorana mode from the side-attached nanowire. We inspect the spin-selective relationship between the trivial Andreev states and the leaking Majorana mode, considering (i) perfectly polarized case, when tunneling of one spin component is completely prohibited, and (ii) another one when both spins are hybridized with the nanowire but with different couplings.

## I. INTRODUCTION

Recent development of the hybrid structures, comprising quantum dots (QDs) coupled the topologically superconducting nanowires [1-4] provide new challenges going beyond mere observation of the Majorana bound states (MBS). Since energy levels of QDs in such hybrid systems are experimentally tunable one can inspect interplay of the topological states with the correlation effects and proximity-induced electron pairing. Due to natural tendency of the Majorana modes to exist at boundaries of finite size systems, one may expect their leakage into any side-attached quantum dot (QD) [5] or more complex magnetic nanoislands [6].

It turns out, however, that efficiency of such process substantially depends on various parameters. For instance, in a weak coupling limit the Majorana modes show up either by the constructive or destructive interferometric lineshapes [7]. True leakage of Majorana mode into nontopological region is possible only in strongly hybridized structures [8-11] as indeed reported experimentally [1-4]. But even under such conditions, leakage of the Majorana mode into the side-coupled QD region could be affected by various effects, e.g. by correlations. We address such issue here, studying an interplay between electron pairing and Coulomb repulsion of the correlated QD proximitized to superconducting substrate where the quantum phase transition from the spinless (BCS-like) to the spinfull (singly occupied) configurations may occur. We analyze whether the Majorana mode would be able to leak into such configurations of the quantum dot.

The paper is organized as follows. We present the physical situation of our interest (Sec. II) and introduce the relevant microscopic model (Sec. III). Next, we con-

∗ e-mail: j.baranski@law.mil.pl

† e-mail: doman@kft.umcs.lublin.pl

FIG. 1. Schematic view of the correlated quantum dot (QD) side-coupled to the Rashba nanowire and both deposited on superconducting (SC) substrate, where topological superconducting nanowire hosts two Majorana end-modes η 1 and η 2 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00349v1-figure-b8175110-8032-4425-856c-55cf56b6017f.jpg)

sider leakage of the Majorana mode into the quantum dot for the fully (Sec. IV) and partly polarized (Sec. V) cases. Finally, we summarize the results (Sec. VI) and give a brief reminder about the singlet-doubled transition in absence of Majorana mode (Appendix A).

## II. PHYSICAL SETUP

We consider the strongly correlated quantum dot attached to the Rashba nanowire, both proximitized to the s -wave superconductor (SC). This setup (Fig. 1) could be realized in the scanning tunneling microscopy (STM) measurements, using Fe or Co nanoscopic chains on the superconducting Pb [12], Al or Re [4] substrates.

In realistic situations the spin-orbit coupling along with the Zeeman effect break spin-rotational symmetry of the system. Such nanowire brought in contact with superconductor develops the intersite pairing of parallel spins. They are 'tilted' with respect to z -axis, but one can project such pairing onto ↑ and ↓ components. For each of these sectors the intersite-pairing is characterized by effective amplitudes, mainly dependent on the applied magnetic field [13].

Upon entering the topological phase, there appear two Majorana quasiparticles simultaneously in both spin channels but with different spectral weights, what has been indeed reported in the STM measurements using a ferromagnetic tip [12]. Such issue has been recently addressed by several groups [11, 13, 14]. J. Klinovaja with co-workers [8] have shown that spin up and down tunneling amplitudes between MBSs and quantum dot depend on the spin-orbit interaction length. Distance between the quantum dot and topological nanowire determines oscillations of the tunneling amplitudes. By changing such dot-nanowire distance or varying the magnetic field one can thus tune the hybrid system either to the fully or to the partly polarized tunneling regimes.

/negationslash

For perfect spin polarization of the QD-chain tunneling (e.g. t ↑ = 0, t ↓ = 0) one would expect signatures of the zero-energy mode to appear only in one spin channel (for ↑ electrons). However, the proximity induced on-dot pairing between opposite spins mixes both these channels. Some aspects of this situation have been addressed in Refs [7, 15, 16]. Since leakage of the Majorana (zero-energy) mode is sensitive to electronic states near the Fermi level, we would like to focus on crossing of the subg-gap (Andreev) states caused by competition between on-dot pairing and Coulomb repulsion (which is spectroscopically manifested by the singlet-doublet quantum phase transition). We show that signatures of the Majorana mode look completely different for both spin channels. In our considerations we take into account the case of (i) perfectly polarized tunneling of only one spin component while the other one is completely prohibited and (ii) another case where both spin electrons can be tunnel transferred, but with different amplitudes.

## III. MICROSCOPIC MODEL

To inspect the mutual relationship between the Majorana mode and in-gap features of the correlated QD we restrict our attention to the limit of large energy gap ∆ of the superconducting reservoir. Under such conditions the single level quantum dot is affected by the proximityinduced pairing. The proximitized QD is described by the Hamiltonian ∑ σ /epsilon1d † σ d σ + Un n ↓ ↑ + Γ S 2 ( d d ↑ ↓ + d d † ↓ † ↑ ), where local pairing is represented by the pair source/sink terms. We next consider coupling to the Majorana zero mode (MZM) represented by ∑ σ λ ( d † σ η 1 + η d 1 σ ) + i/epsilon1 m η η 1 2 . In our approach we recast the self-hermitian operators η † 1 2 , = η 1 2 , by the standard fermionic ones [17] η 1 = 1 √ 2 ( f + f † ) and η 2 = -i √ 2 ( f -f † ). Summarizing, the effective low energy Hamiltonian can be expressed as

$$\underset { \sigma } { \text{a-} } \underset { \text{tes} } { \text{if} } \, H _ { Q D } ^ { \text{eff} } & \simeq \sum _ { \sigma } \epsilon d _ { \sigma } ^ { \dagger } d _ { \sigma } + U n _ { \downarrow } n _ { \uparrow } + \frac { \Gamma _ { S } } { 2 } ( d _ { \uparrow } d _ { \downarrow } + d _ { \downarrow } ^ { \dagger } d _ { \uparrow } ^ { \dagger } ) \quad ( 1 ) \\ \text{in} \\ \text{h is} \quad \, + \sum _ { \sigma } t _ { m \sigma } ( d _ { \sigma } ^ { \dagger } - d _ { \sigma } ) ( f + f ^ { \dagger } ) + \epsilon _ { m } \left ( f ^ { \dagger } f + \frac { 1 } { 2 } \right ). \\ \text{an-} \quad \, \text{We assume that both spin components } ( \sigma = \uparrow, \downarrow ) \, \text{can be}$$

We assume that both spin components ( σ = ↑ ↓ , ) can be transferred between dot and Majorana host and t mσ = λ / σ √ 2 represent the hopping integrals of such processes.

Our objective is to determine the Green's functions G ( ω ) = 〈〈 Ψ;Ψ † 〉〉 defined in the Nambu matrix notation Ψ = ( † d , d ↑ † ↑ , d ↓ , d † ↓ , f, f † ). The equation of motion technique applied to noncorrelated problem yields

$$\lim _ { U = 0 } \mathcal { G } ^ { - 1 } ( \omega ) = \begin{pmatrix} \omega - \epsilon + i \Gamma _ { N } / 2 & 0 & ( \\ 0 & \omega + \epsilon + i \Gamma _ { N } / 2 & - \Gamma \\ 0 & - \Gamma _ { S } / 2 & \omega - \epsilon + \\ \Gamma _ { S } / 2 & 0 & ( \\ - t _ { m \uparrow } & t _ { m \uparrow } & - t \\ - t _ { m \uparrow } & t _ { m \uparrow } & - t \end{pmatrix} \\ \text{To account for the correlation effects one can numerically dia},$$

$$\begin{smallmatrix} 0 & 0 & \Gamma _ { S } / 2 & - t _ { m \uparrow } & - t _ { m \uparrow } \\ \colon + i \Gamma _ { N } / 2 & - \Gamma _ { S } / 2 & 0 & t _ { m \uparrow } & t _ { m \uparrow } \\ \cdot \Gamma _ { S } / 2 & \omega - \epsilon + i \Gamma _ { N } / 2 & 0 & - t _ { m \downarrow } & - t _ { m \downarrow } \\ 0 & 0 & \omega + \epsilon + i \Gamma _ { N } / 2 & t _ { m \downarrow } & t _ { m \downarrow } \\ t _ { m \uparrow } & - t _ { m \downarrow } & t _ { m \downarrow } & \omega - \epsilon _ { m } & 0 \\ t _ { m \uparrow } & - t _ { m \downarrow } & t _ { m \downarrow } & 0 & \omega + \epsilon _ { m } \end{smallmatrix} \,. \quad \, ( 2 ) \\ \intertext { numerically diagonalize } 8 \times 8 \text{ Hamiltonian matrix, determining the }$$

To account for the correlation effects one can numerically diagonalize 8 × 8 Hamiltonian matrix, determining the eigenenergies and transition elements between them. Another equivalent route can rely on the superconducting atomic limit solution [18], extending it to the present model

$$\mathcal { G } ^ { - 1 } ( \omega ) = \left ( \begin{array} { c c c c c c } G ( \omega ) & 0 & 0 & F ( \omega ) & 0 & 0 \\ 0 & - G ^ { * } ( - \omega ) & F ^ { * } ( - \omega ) & 0 & 0 & 0 \\ 0 & F ^ { * } ( - \omega ) & G ( \omega ) & 0 & 0 & 0 \\ F ( \omega ) & 0 & 0 & - G ^ { * } ( - \omega ) & 0 & 0 \\ 0 & 0 & 0 & 0 & \frac { \omega - c _ { m } } { \omega + c _ { m } } & 0 \\ 0 & 0 & 0 & 0 & 0 & \frac { 1 } { \omega + t _ { m } } \end{array} \right ) ^ { - 1 } - \left ( \begin{array} { c c c c c c } 0 & 0 & 0 & 0 & t _ { m \uparrow } & t _ { m \uparrow } \\ 0 & 0 & 0 & 0 & - t _ { m \uparrow } & - t _ { m \uparrow } \\ 0 & 0 & 0 & 0 & t _ { m \downarrow } & t _ { m \downarrow } \\ 0 & 0 & 0 & 0 & - t _ { m \downarrow } & - t _ { m \downarrow } \\ t _ { m \uparrow } & - t _ { m \uparrow } & t _ { m \downarrow } & - t _ { m \downarrow } & 0 & 0 \\ \end{array} \right ) ( 3 ) \\ \text{where}$$

where

$$G ( \omega ) & = \frac { \alpha _ { s } \, u _ { d } ^ { 2 } } { \omega - \left ( \frac { U _ { d } } { 2 } + E _ { d } \right ) } + \frac { \beta _ { s } \, v _ { d } ^ { 2 } } { \omega - \left ( \frac { U _ { d } } { 2 } - E _ { d } \right ) } + \frac { \alpha _ { s } \, v _ { d } ^ { 2 } } { \omega + \left ( \frac { U _ { d } } { 2 } + E _ { d } \right ) } + \frac { \beta _ { s } \, u _ { d } ^ { 2 } } { \omega + \left ( \frac { U _ { d } } { 2 } - E _ { d } \right ) }, \\ F ( \omega ) & = \frac { \alpha _ { s } \, u _ { d } v _ { d } } { \omega + \left ( \frac { U _ { d } } { 2 } - E _ { d } \right ) }, - \frac { \beta _ { s } \, u _ { d } v _ { d } } { \omega + \left ( \frac { U _ { d } } { 2 } - E _ { d } \right ) }, + \frac { \beta _ { s } \, u _ { d } v _ { d } } { \omega + \left ( \frac { U _ { d } } { 2 } - E _ { d } \right ) },$$

$$F ( \omega ) = \frac { \stackrel { \searrow } { \alpha _ { s } } u _ { d } v _ { d } } { \omega - \left ( \frac { U _ { d } } { 2 } + E _ { d } \right ) } - \frac { \stackrel { \searrow } { \beta _ { s } } u _ { d } v _ { d } } { \omega - \left ( \frac { U _ { d } } { 2 } - E _ { d } \right ) } - \frac { \stackrel { \searrow } { \lambda _ { s } } u _ { d } v _ { d } } { \omega + \left ( \frac { U _ { d } } { 2 } + E _ { d } \right ) } + \frac { \stackrel { \searrow } { \beta _ { s } } u _ { d } v _ { d } } { \omega + \left ( \frac { U _ { d } } { 2 } - E _ { d } \right ) },$$

with the energy E d = √ ( /epsilon1 + U/ 2) 2 +(Γ S / 2) , the usual 2 BCS-type coefficients u 2 d = 1 2 [ 1 + /epsilon1 + U/ 2 E d ] = 1 -v 2 d and the relative spectral weights

$$\alpha _ { s } = \frac { e ^ { \beta U / 2 } + e ^ { - \beta E _ { d } } } { 2 e ^ { \beta U / 2 } + e ^ { - \beta E _ { d } } + e ^ { \beta E _ { d } } } = 1 - \beta _ { s }. \quad ( 6 ) \quad \epsilon / \Gamma _ { S } 0 ^ { + } \\ - 1 \cdot$$

In what follows we will quasiparticle spectrum of the correlated quantum dot, assuming that the topological nanowire is long enough so that any overlap between the end-modes is negligible ( /epsilon1 m /similarequal 0).

## IV. FULLY POLARIZED CASE

Let us first consider the fully spin-polarized case when tunneling of one spin component, say ↓ , is totally prohibited. This situation could occur for very strong magnetic fields applied along the topological nanowire. We show that combined effect of MZM leakage and proximity induced pairing give rise to zero modes apparent in energy spectrum of both spin components even if tunneling rate of one spin is turned off. We inspect the influence of these zero states on the characteristic features of a phase transition from the BCS-like singlet state ( S = 0) to the correlation-dominated doublet configuration. We show that even though zero modes appear in both spin channels, each of them have completely different character.

/negationslash

When a quantum dot or other nanoobject is tunnel coupled to the Rashba chain hosting Majorana particles some part of these zero-energy modes can be transferred to QD region. In consequence, the spectral function of QD reveals additional feature pinned to the Fermi level. For the perfectly polarized tunneling amplitude ( t m ↓ = 0 and t m ↑ = 0) one would expect such effect to modify only the spectral function of the directly coupled spin component [19-22]. However, due to the superconducting proximity effect the electrons (say ↑ ) leaking from nanochain into QD region are bound into local pairs with their opposite spin ( ↓ ) partners [7, 16]. Thereby, the Majorana mode becomes apparent in the energy spectrum of electrons for which the direct tunneling to nanochain is prohibited. This feature in ↓ spin sector, however, should be considered merely as a superconducting response for the leaking Majorana mode of the directly coupled spin ↑ . Different nature of these zero-energy states is evident, both for the uncorrelated case (Fig. 2) and in a behavior of the Andreev states of the correlated quantum dot near a quantum phase transition from the BCS-type singlet to the correlation driven doublet state (Fig. 3).

## A. Energy spectrum of uncorrelated QD

Let us focus first on the noncorrelated case, U = 0. By inspecting the matrix Green's function (2) in the limit Γ N → 0 we notice that it is characterized by five poles: zero-energy state, two Andreev bound states

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00349v1-figure-1cb56842-b1cd-42f0-9f0b-831990615f94.jpg)

FIG. 2. Density of states for spin ↑ (left) and ↓ (right) electrons obtained for t m ↑ = 0 6Γ . S , t m ↓ = 0 in absence of correlations ( U = 0).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00349v1-figure-cf55c2ad-6611-49ce-9439-bfa17228c35b.jpg)

at ± √ /epsilon1 2 +(Γ S / 2) 2 and another two 'molecular' states ± √ /epsilon1 2 +(Γ S / 2) 2 +(2 t m ↑ ) 2 , resulting from hybridization of the Andreev bound states with the Rashba chain [7]. In Fig. 2 we visualize qualitative differences between the spectral weights of ↑ and ↓ spin sectors obtained in the strong coupling limit t ↑ = 0 6Γ . S .

In the spectrum of ↑ electrons we observe the well pronounced zero-energy mode directly leaking from the proximitized Rashba chain in much the same fashion as initially predicted for the quantum dot attached to the Kitaev wire [5]. Contrary to that, the opposite spin electrons are not directly coupled to the chain therefore zeroenergy mode emerges as a consequence of the local pairing. Some tiny Majorana mode shows up in this ↓ spin sector only when the QD energy /epsilon1 is close to the Fermi level, otherwise it quickly fades away.

Furthermore, one can notice a small dip in the zeroenergy signature of spin ↑ electrons appearing near /epsilon1 = 0. This comes from destructive feedback effect when the original QD level coincides with Majorana mode. Such destructive interference pattern has been described by our group [7, 16]. We have pointed out that even when spectral weight of the zero-energy mode in spectrum of directly coupled spins is suppressed, the zero-energy mode in opposite spin channel would be enhanced. As the QD-chain tunneling amplitude of spin ↓ electrons is turned off, spin ↓ electrons do not take part in destructive interference. Consequently in the left panel of Fig. 3 we observe well shaped zero mode. It is worth noticing that spectral weights of original Andreev states ( ± √ /epsilon1 2 +(Γ S / 2) ) in spectrum of 2 ↑ electrons are reduced in comparison with the molecular ones. Contrary to that, for electrons with prohibited dot-chain tunneling we observe that pure Andreev states are dominant.

## B. Majorana near singlet-doublet crossover

We shall now discuss the correlation effects, driven

FIG. 3. Density of states for spin ↑ (upper panel) and ↓ (bottom panel) electrons obtained for t m ↑ = 0 3Γ . S , t m ↑ = 0 and U/ Γ S = 0 (left panel), U/ Γ S = 1 (middle panel), U/ Γ S = 3 5 (right panel). . Following Ref. [23] we denote by ξ d = /epsilon1 -U/ 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00349v1-figure-7df3ac46-62d3-43ae-835e-6b1fd6fa1729.jpg)

by the Coulomb repulsion U . In absence of the Rashba nanowire upon varying the ratio U/ Γ S there appear the regular Andreev bound states at ± [ U/ 2 -√ ( /epsilon1 + U/ 2) 2 +(Γ S / 2) ] 2 which eventually cross each other at the singlet-doublet (see Appendix A for an explanation of this quantum phase transition). Proximity induced zero-energy feature in the spin ↓ sector does not affect significantly this characteristic crossing, because it does not originate directly from the leaking Majorana mode. The strong Coulomb repulsion disfavors on-dot pairing, therefore upon increasing U/ Γ S also the side effects of the local pairing are gradually suppressed. One of them is the zero-energy mode appearing in the spectrum of spin ↓ electrons, originating solely from electron pairing. For this reason, upon traversing the singlet-doublet quantum phase transition the zero-energy mode of spin ↓ electrons is gradually washed out from the strongly correlated regime (see right bottom lower panel in Fig. 3).

In the case of spin ↑ electrons the spectrum looks completely different. Upon increasing the correlations strength the ABS states are split and form the molecular branches. These new states no longer cross each other. Approaching the singlet-doublet transition we observe emergence of the zero-energy state (directly leaking from the Rashba nanochain) while the molecular Andreev states repel one another instead of the regular crossing (shown in Fig. 5). Such avoided-crossing behavior (see the upper right panel in Fig. 3) can be regarded as additional signature, of the proximity induced Majorana mode that repels the 'trivial' (finite-energy) Andreev states. This result generalizes the previously discussed topological/nontopological different nature of the bound states in hybrid structures comprising the uncorrelated quantum dots attached to the topological superconducting wires hosting the Majorana modes [3, 8, 9, 11, 24]. We hope that empirical detection would be feasible using the spin-polarized STM measurements and varying the level /epsilon1 of correlated QD by the gate potential.

## V. PARTIAL SPIN-POLARIZATION

Theoretical studies indicate that for specific conditions, it would be possible to achieve nearly perfectly spin polarized tunneling between the quantum dot and chain [8]. However in realistic experimental setups both spin electrons can be transferred, although with substantially different tunneling amplitudes. In this subsection we briefly address such partially polarized case. Motivated by the spin polarized STM measurements [12] we assume in our approach the tunneling rate of spin ↓ electrons to be t m ↓ = 0 1Γ . S which is three times weaker than the coupling amplitude for spin ↑ electrons ( t m ↑ = 0 3Γ . S ). Fig. 4 compares the density of states of each spin electrons in the non-correlated (left panels) and strongly correlated regime (right panels), respectively.

In the strongly correlated limit we can practically observe a convolution of the features typical for each spin

FIG. 4. Density of states for spin up (upper panel) and spin down (bottom panel) electrons obtained for t m ↑ = 0 3Γ . S , t m ↓ = 0 1Γ . S and U/ Γ S = 0 (left panel), U/ Γ S = 3 5 . (right panel).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00349v1-figure-3c48da6c-9291-4586-948c-bb79813137af.jpg)

sectors of the fully polarized case. Yet one can clearly resolve the avoided-crossing behavior in the dominant ↑ spin coupling channel with the well pronounced Majorana mode separating them. Spectrum of the opposite spin sector, on the other hand, is predominantly reminiscent of the continuous Andreev band state branches with only very residual signatures of the zero-energy modes. Correlations could thus be very useful for distinguishing the qualitatively different character of the spin-polarized spectra of the quantum dot.

## VI. SUMMARY

We have studied the energy and spin-dependent spectra of the proximitized correlated quantum dot attached to the topologically superconducting nanochain, hosting the Majorana end-modes. We have analyzed influence of the correlations (responsible for a quantum phase transition from the spinless BCS-type to the spinful configuration) on efficiency of the Majorana mode leakage. Our study predicts that the zero-energy states might appear in both spin sectors (due to the superconducting proximity effect), however their spectroscopic signatures are going to be qualitatively different. The spin channel which is directly coupled to the Majorana mode is characterized by (a) the well separated Andreev branches (of an avoided-crossing behavior) coexisting with (b) the zeroenergy feature of a sizable spectral weight. The other spin sector, which is not directly coupled to the Majorana mode, in the strongly correlated limit is characterized by the continuous Andreev branches traversing the zero energy without any trace of the leaking Majorana mode. Such spin sector can eventually allow the Majorana mode to appear but only in the weak correlation regime, when the proximitized quantum dot stays in the BCS-type configuration.

Correlations can hence be very beneficial for spinselective leakage of the Majorona modes on the quantum dots. Such effects would be easily detectable, using the polarized STM measurement analogous to those already reported in Ref. [12].

FIG. 5. Spectral function ρ ω ( ) of the correlated quantum dot as a function of the energy ξ d = /epsilon1 + U/ 2 obtained at low temperature T = 0 1Γ . S in the 'superconducting atomic limit' Γ N → 0 + for several values of the Coulomb potential U = 0 (upper panel) U = Γ S (middle panel), U = 2Γ S (bottom panel).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00349v1-figure-5061141a-cb89-4183-9b19-e4b8705f4e47.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00349v1-figure-2e5190d8-9230-4a0e-a73e-18b0fc2d50cb.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00349v1-figure-20124b87-c172-42ad-803d-871533f1a25e.jpg)

## ACKNOWLEDGMENTS

This research has been conducted in the framework of the project 'Analysis of nanoscopic systems coupled with superconductors in the context of quantum information processing' No. GB/5/2018/209/2018/DA funded in period 2018-2021 by the Ministry of National Defense Republic of Poland (J.B.,T.Z.). This project has been also supported by National Science Centre (NCN, Poland) under the grants UMO-2018/29/B/ST3/00937 (G.G.) and UMO-2017/27/B/ST3/01911 (T.D.).

## Appendix A: Quantum phase transition in absence of Majorana mode

For illustration of the singlet-doublet phase transition we recall briefly the results obtained for the correlated quantum dot in absence of the Rashba chain. For this purpose we analyze the subgap Andreev states in the superconducting atomic following several earlier works

- [1] M. T. Deng, S. Vaitiekenas, E. B. Hansen, J. Danon, M. Leijnse, K. Flensberg, J. Nyg˚rd, a P. Krogstrup, and C. M. Marcus, 'Majorana bound state in a coupled quantum-dot hybrid-nanowire system,' Science 354 , 1557 (2016).
- [2] F. Nichele, A. C. C. Drachmann, A. M. Whiticar, E. C. T. O'Farrell, H. J. Suominen, A. Fornieri, T. Wang, G. C. Gardner, C. Thomas, A. T. Hatke, P. Krogstrup, M. J. Manfra, K. Flensberg, and C. M. Marcus, 'Scaling of Majorana zero-bias conductance peaks,' Phys. Rev. Lett. 119 , 136803 (2017).
- [3] M.-T. Deng, S. Vaitiek˙ enas, E. Prada, P. San-Jose, J. Nyg˚rd, P. Krogstrup, R. Aguado, a and C. M. Marcus, 'Nonlocality of Majorana modes in hybrid nanowires,' Phys. Rev. B 98 , 085125 (2018).
- [4] H. Kim, A. Palacio-Morales, T. Posske, L. R´ ozsa, K. Palot´s, L. Szunyogh, M. Thorwart, a and R Wiesendanger, 'Toward tailoring Majorana bound states in artificially constructed magnetic atom chains on elemental superconductors,' Sci. Adv. 4 , eaar5251 (2018).
- [5] E. Vernek, P. H. Penteado, A. C. Seridonio, and J. C. Egues, 'Subtle leakage of a Majorana mode into a quantum dot,' Phys. Rev. B 89 , 165314 (2014).
- [6] E. Mascot, S. Cocklin, S. Rachel, and D. K. Morr, 'Quantum engineering of Majorana fermions,' ArXiv eprints (2018), arXiv:1811.06664.
- [7] G. G´ orski, J. Bara´ski, I. Weymann, n and T. Doma´ski, n 'Interplay between correlations and Majorana mode in proximitized quantum dot,' Sci. Rep. 8 , 15717 (2018).
- [8] S. Hoffman, D. Chevallier, D. Loss, and J. Klinovaja, 'Spin-dependent coupling between quantum dots and topological quantum wires,' Phys. Rev. B 96 , 045440 (2017).
- [9] E. Prada, R. Aguado, and P. San-Jose, 'Measuring Majorana nonlocality and spin structure with a quantum

[23, 25, 26]. The true eigenstates are represented by the doublet configurations |↑〉 , |↓〉 (corresponding to spin S = 1 2 ) and the BCS-type singlet states u d | 0 〉 + v d |↓↑〉 or -v d | 0 〉 + u d |↓↑〉 ( S = 0). It can be shown that when correlations are weak U &lt; Γ , the ground state is repS resented by singlet state in full range of initial dot level ( /epsilon1 ). In such case, the density of states consists of two Andreev bound states branches separated by the pairing gap (Fig. 5). For stronger correlations, these subgap branches approach each other and can eventually cross at the singlet-doublet phase boundaries. The critical interaction for ground state transition at the half filling is U = Γ S (middle panel in Fig. 5). In the doublet regime ( U &gt; Γ ), the subgap resonances cross each other, givS ing rise to a loop structure (bottom panel in Fig. 5). In section IV, we noticed that the states appearing in spin ↓ electron spectrum exhibit similar behavior (c.f. bottom panels in Fig. 3). Presence of the zero-energy state in this spin channel did not affect significantly the characteristics of the singlet-doublet transition. In contrast, for spins directly coupled to nanochain, we note the qualitatively different avoided-crossing behavior.

dot,' Phys. Rev. B 96 , 085418 (2017).

- [10] A. Ptok, A. Kobia/suppress lka, and T. Doma´ nski, 'Controlling the bound states in a quantum-dot hybrid nanowire,' Phys. Rev. B 96 , 195430 (2017).
- [11] P. Szumniak, D. Chevallier, D. Loss, and J. Klinovaja, 'Spin and charge signatures of topological superconductivity in Rashba nanowires,' Phys. Rev. B 96 , 041401 (2017).
- [12] S. Jeon, Y. Xie, J. Li, Z. Wang, B. A. Bernevig, and A. Yazdani, 'Distinguishing a Majorana zero mode using spin-resolved measurements,' Science (2017), 10.1126/science.aan3670.
- [13] M.M. Ma´ ska and T. Doma´ nski, 'Polarization of the Majorana quasiparticles in the Rashba chain,' Sci. Rep. 7 , 16193 (2017).
- [14] J. Li, S. Jeon, Y. Xie, A. Yazdani, and B. A. Bernevig, 'Majorana spin in magnetic atomic chain systems,' Phys. Rev. B 97 , 125119 (2018).
- [15] R. Chirla and C. P. Moca, 'Fingerprints of Majorana fermions in spin-resolved subgap spectroscopy,' Phys. Rev. B 94 , 045405 (2016).
- [16] J. Bara´ski, n A. Kobia/suppresska, l and T. Doma´ nski, 'Spinsensitive interference due to Majorana state on the interface between normal and superconducting leads,' J. Phys.: Condens. Matter 29 , 075603 (2017).
- [17] S. R. Elliott and M. Franz, 'Colloquium: Majorana fermions in nuclear, particle, and solid-state physics,' Rev. Mod. Phys. 87 , 137 (2015).
- [18] E. Vecino, A. Mart´ ın-Rodero, and A. Levy Yeyati, 'Josephson current through a correlated quantum level: Andreev states and π junction behavior,' Phys. Rev. B 68 , 035105 (2003).
- [19] M. Lee, J. S. Lim, and R. L´ opez, 'Kondo effect in a quantum dot side-coupled to a topological superconductor,' Phys. Rev. B 87 , 241402 (2013).

- [20] R. L´ opez, M. Lee, L. Serra, and J. S. Lim, 'Thermoelectrical detection of Majorana states,' Phys. Rev. B 89 , 205418 (2014).
- [21] I. Weymann, 'Spin Seebeck effect in quantum dot side-coupled to topological superconductor,' J. Phys.: Condens. Matter 29 , 095301 (2017).
- [22] I. Weymann and K. P. W´ ojcik, 'Transport properties of a hybrid Majorana wire-quantum dot system with ferromagnetic contacts,' Phys. Rev. B 95 , 155427 (2017).
- [23] J. Bauer, A. Oguri, and A. C. Hewson, 'Spectral properties of locally correlated electrons in a Bardeen-Cooper-Schrieffer superconductor,'
- J. Phys.: Condens. Matter 19 , 486211 (2007).
- [24] A. Ptok, A. Kobia/suppress lka, and T. Doma´ nski, 'On controlling the bound states in quantum-dot hybrid-nanowire,' (2017), arXiv:1710.06387.
- [25] A. Mart´ ın-Rodero and A. Levy Yeyati, 'The Andreev states of a superconducting quantum dot: mean field versus exact numerical results,' J. Phys.: Condens. Matter 24 , 385303 (2012).
- [26] J. Bara´ski and T. Doma´ski, 'In-gap states of a quann n tum dot coupled between a normal and a superconducting lead,' J. Phys.: Condens. Matter 25 , 435305 (2013).
Tomasz Zienkiewicz Jan Barański Grzegorz Górski Tadeusz Domański
Leakage of Majorana mode into correlated quantum dot nearby its singlet-doublet crossover
We study quasiparticle spectrum of the correlated quantum dot deposited on
superconducting substrate which is side-coupled to the Rashba nanochain,
hosting Majorana end modes. Ground state of an isolated quantum dot
proximitized to superconducting reservoirs is represented either by the singly
occupied site or BCS-type superposition of the empty and doubly occupied
configurations. Quantum phase transition between these distinct ground states
is spectroscopically manifested by the in-gap Andreev states which cross each
other at the Fermi level. This qualitatively affects leakage of the Majorana
mode from the side-attached nanowire. We inspect the spin-selective
relationship between the trivial Andreev states and the leaking Majorana mode,
considering (i) perfectly polarized case, when tunneling of one spin component
is completely prohibited, and (ii) another one when both spins are hybridized
with the nanowire but with different couplings.