Polarization of the Majorana quasiparticles in the Rashba chain

We demonstrate that the selective equal–spin Andreev reflection (SESAR) spectroscopy can be used in STM experiments to distinguish the zero–energy Majorana quasiparticles from the ordinary fermionic states of the Rashba chain. Such technique, designed for probing the p–wave superconductivity, could be applied to the intersite pairing of equal–spin electrons in the chain of magnetic Fe atoms deposited on the superconducting Pb substrate. Our calculations of the effective pairing amplitude for individual spin components imply the magnetically polarized Andreev conductance, which can be used to ‘filter’ the Majorana quasiparticles from the ordinary in–gap states, although the pure spin current (i.e., perfect polarization) is impossible.

The topologically nontrivial superconducting state of one-dimensional (1D) chains 1 allows for a unique phenomenon of the selective equal-spin Andreev reflection (SESAR). This polarized Andreev spectroscopy has been proposed by J. J. He et al. 2 as a useful tool for probing the Majorana states. SESAR measurements have indeed provided evidence for the zero-energy modes in vortices of the p-wave superconducting Bi 2 Te 3 /NbSe 2 heterostructures 3,4 . Similar ideas have been also considered for the Josephson-type junctions 5,6 and ferromagnet-superconductor interfaces with the spin-orbit coupling 7,8 . In this work we demonstrate that SESAR spectroscopy can test inherent polarization of the Majorana quasiparticles appearing at the edges of the Rashba chain. The parallel and perpendicular components of magnetically polarized Majorana states has initially been pointed out by D. Sticlet et al. 9 and their signatures have been recently studied by a number of authors [10][11][12][13] . In this paper we show that magnetic polarization is detectable in STM experiments owing to SESAR processes, which in the subgap regime could distinguish the Majorana quasiparticles out the ordinary Shiba in-gap states. We provide microscopic arguments explaining such polarization and confront our predictions with the experimental data obtained for Fe atom chain deposited on the surface of Pb superconductor by the STM technique with use of the magnetically polarized tip 14 .
The underlying idea of SESAR for the aforementioned configuration is displayed in Fig. 1. This STM-type setup has been previously used by several experimental groups [15][16][17] , however, ignoring the magnetic polarization. Recently A. Yazdani and coworkers 14 have measured the spin-resolved tunneling current and revealed substantial polarization of the zero-bias conductance in regions, where the Majorana quasiparticles exist. This fact can be interpreted within the popular microscopic model, taking into account the Rashba and Zeeman interactions in addition to the proximity-induced pairing which can realistically capture a topography of the Majorana fermions 9,18-20 . Using this model we have recently emphasized 21 , that amplitude of the intersite pairing (between identical spin electrons) differs several times for ↑ and ↓ sectors, respectively. Obviously, such effect should give rise to noticeable polarization of the Majorana quasiparticles near the chain edges. In practice, the low-energy features can be detected only by the anomalous Andreev spectroscopy, as discussed in detail in ref. 18 . Since efficiency of the particle to hole conversion for the spin-polarized Andreev spectroscopy depends on the anomalous propagator 1 , 0 , one should expect its non-vanishing value at ω = 0 nearby the chain edges. In what follows we show, that this is really the case. We also argue, that SESAR could distinguish the Majorana from the ordinary fermionic quasiparticles.

Results
Microscopic model. Nanoscopic chain of the magnetic Fe atoms deposited on the s-wave conventional superconductor and probed by the polarized STM tip (relevant to the experimental situation 14 21 ). In what follwos, we impose the constant couplings Γ N and Γ S to the STM tip and superconducting substrate, respectively (see Fig. 1).
The low-energy Hamiltonian is effectively given by 19 annihilates (creates) an electron with spin σ at site i, t ij is the hopping integral and μ is the chemical potential. The proximity effect, responsible for the on-site (trivial) pairing, can be modeled as 20 with the pairing potential Δ = Γ S /2. In this scenario the intersite p-wave pairing is driven by the Rashba and the Zeeman interactions We assume the magnetic field to be aligned along ẑ-axis and impose the spin-orbit vector α = (0, 0, α). Fig. 2 we present spatial dependence of the off-diagonal

Spin-polarized Majorana quasiparticles. In
obtained at zero energy for different spins ↑ and ↓, respectively. This anomalous spectral function is very instructive, because its sign exhibits intrinsic polarization of the Majorana modes (previously emphasized in ref. 9 ) whereas its absolute value can be probed by the SESAR spectroscopy (see the next paragraph). Concerning the magnitude, we clearly notice a quantitative difference (almost 5 times) between the spin ↑ and ↓ inter-site pairings. As regards the intrinsic polarization we observe that F iσ (ω = 0) changes its phase by π between opposite sides of the Rashba chain and furthermore each of the spin sectors is characterized by opposite polarizations. This aspect resembles the results reported for the interface of ferromagnet/superconductor bilayers 22 . Such feature can be regarded as a hallmark of the finite-size systems, because otherwise (i.e., in thermodynamic limit L → ∞) the off-diagonal spectral function would identically vanish at zero energy for both pairing channels. Figure 3 illustrates the spatial profiles of the spin-polarized (diagonal) spectral function ρ iσ (ω). As expected, we notice quantitative differences between the Majorana states appearing in ↑ and ↓ spin sectors, whereas their overall profiles seem to be pretty similar. Different magnitudes of these spin-polarized Majorana quasiparticles would show up in the SESAR measurements. The polarized Andreev transport. By applying a bias voltage V between the STM tip and the superconducting reservoir one would induce the nonequilibrium charge transport. Deep in a subgap regime (i.e., for Δ  V e / ) such current is contributed solely by the Andreev scattering, when electrons from the STM tip are converted into the pairs, reflecting holes back to the STM tip. This process can be treated within the Landauer-Büttiker formalism.
We can express the nonmagnetic (γ = 0) and magnetically polarized (γ = σ) Andreev currents by the following formula stands for the Fermi-Dirac distribution function. These Andreev channels are characterized by various (dimensionless) transmittances, that can be expressed via the local and non-local anomalous Green's functions, respectively    The corresponding conductances are presented in Fig. 5. We notice that the differential conductance of the nonmagnetic Andreev reflections dominates well inside the Rashba chain at energies coinciding with the fermion Andreev/Shiba states. The SESAR, on the other hand, is efficient mainly near the Majorana modes whose spatial extent covers roughly 10 sites near the Rashba chain edges. In distinction to ref. 2 , we observe that the spin-polarized currents are present for both spins (↑ and ↓) but with significantly different magnitudes. Our results are relevant to the recent experimental data reported by the Princeton group 14 . We have checked that the spin-polarized Majorana quasiparticles are robust upon varying the model parameters, although some additional subtle effects may be observed, for instance the quantum oscillations 18 .
The results presented in Fig. 5 correspond to the topological regime. By varying the model parameters so that the system is driven to the topologically trivial phase, the zero-energy Majorana peak vanishes and the total transmittance in the spin-polarized channels is strongly suppressed. Such evolution from the topologically trivial to nontrivial state is presented in Fig. 6. Note that the polarized transmittance T ↑ (ω) vanishes almost completely outside the topological regime. In the topological regime the unpolarized transmittance of the Majorana peak is much smaller than the transmittance of the ordinary in-gap states that develop when the system enters the topological regime. On the other hand, the polarized transmittance of the Majorana peak is much larger than the ordinary in-gap states.
In summary, we emphasize that the net spin current = − ↑ ↓ , attainable from the SESAR spectroscopy, is expected to acquire meaningful values of the spatially-resolved conductance only near the Majorana quasiparticles (what can be inferred by inspecting Fig. 5). SESAR can hence filter the Majorana from the ordinary Andreev/Shiba quasiparticles (which always exist in the Rashba chain). This unique virtue of SESAR would be valuable for spotting the Majorana quasiparticles and investigating their topography.

Discussion
We have studied the selective equal-spin Andreev spectroscopy (SESAR) which can empirically detect the polarized Majorana quasiparticles appearing at the edges of the Rashba chain. We have shown that different amplitudes of the inter-site equal-spin pairing imply the magnetic polarization of the Majorana states and yields the spindependent Andreev transport with substantially distinct probabilities in each spin components. Our theoretical results qualitatively agree with the recent finding by A. Yazdani 14 , who reported the spin-polarized features in the subgap spectroscopy. Even though the pure spin current (discussed in ref. 2 ) is impossible -the spin current conductance G V ( ) i spin could nevertheless filter the Majorana quasiparticles from the ordinary Andreev/Shiba states. Our quantitative estimations clearly show also that the non-polarized and spin-polarized Andreev conductances are much smaller than the unitary limit value 2e 2 /h as has been indeed observed by the STM 15,16,23 and by the tunneling measurements via heterojunctions 24 .

Methods
Our calculations have been performed for the Rashba chain, comprising L = 70 atoms. In most of the numerical calculations (except Fig. 6) we have used the following model parameters: magnitude of the induced pairing Δ = 0.2 t, the spin-orbit coupling α = 0.15 t, the chemical potential μ = −2.1 t, and the external magnetic field gμ B B/2 = 0.27 t. Such a choice of parameters locates the system strictly in a topological regime 21 . The spinresolved spectral functions, presented in Fig. 3, have been calculated using the following definition Let us outline a brief scheme for computing the charge tunneling current induced through i-th site of the chain coupled between the STM tip (N electrode) and the superconducting substrate (S electrode), for simplicity neglecting the inter-site hopping t ij = 0. Using the Heisenberg equation we can express such current as denotes the hybridization of i-th site with itinerant electrons of the tip. Since we are interested in the spin-resolved spectroscopy let us exprees (9) as  with the bare Green's functions The lesser matrix Green's function obeys the Keldysh equation r r a a r a , where for brevity we dropped the temporal arguments. In our case the first term vanishes, so we are left with