Stability of Majorana bound states in the presence of spin-flip scattering
Subhajit Pal Colin Benjamin

  
  
  colin.nano@gmail.comSchool of Physical Sciences, National Institute of Science Education & Research, HBNI, Jatni-752050, India 
74.50.+r,74.78.Na,85.25.-j,85.25.Cp,85.65.+h,75.50.Xx,85.80.Fi

A popular evidence of the existence of Majorana bound states(MBS) is a quantized zero-bias conductance peak(ZBCP) which is robust to scattering by impurities, a consequence of its topological protection. In this work we examine the stability of this MBS induced ZBCP in a metal-superconductor junction in the vicinity of a spin flipper. We analytically calculate the differential charge conductance for metal-spin flipper-superconductor junctions with two distinct types of superconductors: (a) spin less p-wave superconductor and (b) spin-orbit-coupled s-wave superconducting wire in presence of a Zeeman field(SOCSW). We see that the quantized ZBCP remains stable in presence of spin flip scattering for metal-p-wave superconductor junction, while it loses its stability when the p-wave superconductor is replaced by a SOCSW. Further, the scattering matrix of the normal metal-p-wave superconductor junction satisfies BDI symmetry class regardless of the presence or absence of spin flip scattering. In the BDI symmetry class both Hamiltonian as well as scattering matrix satisfies particle-hole, time reversal and chiral symmetries. However, in case of normal metal-SOCSW junction the Hamiltonian as well as scattering matrix belongs to symmetry class D in absence of spin flip scattering. In the symmetry class D both Hamiltonian and the scattering matrix satisfy particle-hole symmetry relation, but do not satisfy time reversal and chiral symmetry relations. In presence of spin flip scattering the scattering matrix for SOCSW belongs to symmetry class A for which the scattering matrix does not satisfy either particle-hole or time reversal or chiral symmetry relations. The reason for ZBCP at a metal-p wave superconductor junction is perfect Andreev reflection regardless of spin flip scattering, while for a metal-SOCSW junction it is the exact cancellation between normal and Andreev reflection probabilities at zero bias and not perfect Andreev reflection, in absence of spin flip scattering. In presence of spin flip scattering, on the other hand, there is no exact cancellation at zero bias which leads to absence of quantized ZBCP for a metal-SOCSW junction. 


Stability of Majorana bound states in the presence of spin-flip scattering
    Colin Benjamin
    December 30, 2023
==========================================================================



§ INTRODUCTION

In recent years Majorana bound states<cit.> (MBS) have been studied both theoretically and experimentally in many setups. MBS are quasiparticle excitations within the superconductor at zero energy which are their own anti-particles<cit.>. These Majorana zero modes can be observed at the boundary with a topological superconductor<cit.>. Possible signatures of MBS have been reportedly seen in nanowire-superconductor hybrid device<cit.>, topological insulator-superconductor heterostructure<cit.>, atomic chains on a superconductor<cit.> and superconducting nanowire between two ferromagnetic leads<cit.>. Further there have been reports of the generation of Majorana states at interfaces of normal metals and spin orbit coupled superconducting wire (SOCSW)<cit.>. Our main motivation in this work is to propose a reliable check on the stability of Majorana bound states at metal-superconductor interfaces. To show this we consider a spin flipper at a metal-superconductor interface, and see that the quantized zero bias conductance peak (ZBCP) remains stable in presence of spin flip scattering for p-wave superconductor(pSc) but unfortunately it loses its stability due to spin flip scattering in case of a  SOCSW. As a consequence, ZBCP does not take universal values when a spin flipper is placed at interface of a metal-SOCSW. 

To better understand the stability of MBS to spin-flip scattering, we study the symmetry of the scattering matrix(S-Matrix) of a Metal-pSc as well as a Metal-SOCSW junction. In 1D, the Hamiltonian of a spinless pSc preserves all three symmetries- particle-hole, time reversal and chiral, implying that the Hamiltonian for a spinless pSc is in topological class BDI<cit.>. An one dimensional normal metal-spin less pSc junction is considered in Ref. [setiawan], and by analytically calculating the different reflection probabilities it is shown that the zero bias conductance is robustly quantized at 2e^2/h in the topological regime. The Hamiltonian in Ref. [setiawan] for the normal metal-pSc junction belongs to symmetry class BDI. In Ref. [fulga], another way to determine the symmetry class is proposed, which is using the symmetries of the scattering matrix. The symmetries of the S-Matrix can also be classified into ten symmetry classes<cit.>. In Ref. [setiawan] S-Matrix for normal metal-pSc junction satisfies particle hole, time reversal and chiral symmetries and thus again belongs to the symmetry class BDI which  has a topological nature in 1D. Both Hamiltonian and S-Matrix for the normal metal-pSc junction in Ref. [setiawan] belong to the BDI symmetry class. In our work, there is a spin flipper at metal-pSc interface. We find that irrespective of spin flip scattering at the interface between metal and pSc, the zero bias S-Matrix always lies in the symmetry class BDI. The topological quantum number for class BDI is the number of negative eigenvalues of the S-Matrix, which does not change with change of parameters regardless of spin flip scattering. Thus, the topological properties of the predicted MBS at a Metal-pSc junction are robust to spin-flip scattering.

Further, in Ref. [setiawan], similar to metal-pSc junction the metal-SOCSW junction is also studied analytically and it is seen that the zero bias conductance peak (ZBCP) is quantized. The normal Metal-SOCSW junction Hamiltonian belongs to the symmetry class D, which in 1D is topological. The S-Matrix thus also belongs to the same symmetry class D. In Ref. [kli], Majorana Fermion wave functions occurring in a Metal-SOCSW junction are derived. The SOCSW Hamiltonian used in Ref. [kli] is similar to that of Ref. [setiawan] and the SOCSW Hamiltonian in Ref. [kli] is also in the symmetry  class D. There have been two previous works which deal with a 1D linearized SOCSW Hamiltonian, Ref. [setiawan] and Ref. [kli]. In these two papers, the linearized SOCSW Hamiltonian has two branches: interior and exterior. The interior branch Hamiltonian as well as the exterior branch Hamiltonian for Ref. [setiawan] is in D symmetry class, while that in Ref. [kli] is in symmetry class DIII. On the other hand, the S-Matrix of Normal metal-SOCSW junction of Ref. [setiawan], satisfies only particle-hole symmetry relation. Therefore, both Hamiltonian and S-Matrix for normal metal-SOCSW  junction in Ref. [setiawan] belongs to the symmetry class  D for which the topological quantum number is the determinant of the S-Matrix. In our case, when spin flipper does not cause any spin flip, we see similar behavior as in Ref. [setiawan], i.e., the S-Matrix for normal metal-SOCSW junction satisfies only particle-hole symmetry relation and there is a ZBCP which is brought out not because of perfect Andreev reflection as in normal metal-pSc junction but because of exact cancellation of normal and Andreev reflection, a fact already noticed in Ref. [setiawan]. For a normal metal-SOCSW junction,
we find the topological quantum number for class D (Determinant of the S-Matrix) does not change with change of parameters. But, in presence of spin flip scattering, the S-Matrix for normal metal-SOCSW junction has no symmetry and is in class A and the topological quantum number for class D, i.e., determinant of the S-Matrix changes with change of parameters. In presence of spin flip scattering the determinant of the S-Matrix for SOCSW is complex with absolute value 1. Thus, a 1D normal metal-SOCSW junction which purportedly hosts a MBS loses whatever topological nature it possesses when subject to spin flip scattering and as a result the ZBCP is no longer stable at a metal-SOCSW interface.

The paper is organized as follows: in the next section we first present our model for spin-flip scattering at a normal metal-pSc interface  and then at a normal metal-SOCSW interface by writing the Hamiltonian, wavefunctions and boundary conditions to determine the different scattering probabilities. Following this we discuss our results by comparing the zero-bias conductance in metal-SOCSW junction with that in metal-pSc junction in presence of spin flip scattering. We next discuss the symmetry class of the S-Matrix for both Metal-pSc junction and Metal-SOCSW junction in presence of spin flip scattering. Finally, we conclude with a brief summary of our work.



§ THEORY

In this work we contrast the Majorana states arising at metal-pSc interfaces with those arising at metal-SOCSW interfaces. The exact setting we will use is shown in Figs. 1(a) and 2(a), it represents a spin flipper at x=-a while at x=0 a δ-like potential barrier which separates the normal metal from pSc or SOCSW. The regions I (x<-a) and II (-a<x<0) are metallic while for x>0 there is a pSc or SOCSW in the topological regime.



 §.§ Spin flipper in the vicinity of Metal-pSc junction
 

  §.§.§  Hamiltonian
 We consider a one-dimensional normal metal (NM)-normal metal (NM)-pSc junction wherein a spin flipper is embedded between two metallic regions at x=-a. The interface at x=0 (see Fig. 1) is modeled by a δ-like potential barrier (strength Z) and the problem is solved using the BTK<cit.> approach. The normal metal is spinful and spinflip processes are confined to the normal metal.
When a spin up/down electron with energy E is incident from the metallic region, at the x=-a interface, it interacts with the spin flipper through an exchange interaction which may cause a mutual spin flip. The incident electron can be reflected back to metallic region I, or transmitted to metallic region II, with spin up or down. When this transmitted electron with spin-up is incident at x=0 interface it can be normally reflected back as an electron with spin-up from the interface or could be Andreev reflected, only as spin down hole back to region II. In the p-wave superconductor the spin up electron and spin down hole form one channel of transport, while the spin down electron and spin-up hole form another channel of transport for energies above the superconducting gap. There is no spin mixing in the p-wave superconductor. These are two separate channels. So, in our case the p-wave superconductor is still spinless regardless of the presence of spin flipper in its vicinity. Similar studies on spinful semiconductor coupled with spinless s-wave superconductor junction where spin flip processes are confined to the semiconductor region are dime a dozen. For e.g., in Ref. [cheng], spin-orbit coupling is only present in the semiconductor in a superconductor-semiconductor-superconductor junction, while the superconductor is spinless.

In our work spin flipper is a point like magnetic impurity similar to Ref. [AJP]. The effect of magnetic impurities on the quantum transport is vast<cit.>. This includes novel Kondo effect in which the conduction electrons scatter in a metal due to magnetic impurities. This results in a characteristic change in electrical resistivity with temperature<cit.>. In a diffusive normal metal-insulator-s-wave superconductor junction, normalized zero-bias conductance is suppressed by magnetic impurity scattering in diffusive normal metal for low transparent junction, while for high transparent junction zero-bias conductance is enhanced by magnetic impurity scattering<cit.>. Further, the effect of magnetic impurity on Majorana bound state is a topic of current interest<cit.>.
The Hamiltonian for spin flip interaction, from Refs. <cit.> is

    H_Spin flipper   =   -Js⃗.S⃗'⃗, 
    s⃗.S⃗'   =   s_zS'_z+1/2(s^-S'^++s^+S'^-),


J denotes the strength of exchange coupling between electron's spin s⃗ and  spin S⃗'⃗ (of spin-flipper).  S'^± = S'_x± iS'_y are the spin raising and lowering operators for spin flipper while s^± = s_x± is_y are the spin raising and lowering operators for electron/hole with s_k=ħ/2[ σ_k   0;   0 σ_k ], k=x,y,z, with σ_k being the usual Pauli matrices. 
When the electron interacts with the magnetic impurity there are two possible processes, one is ‘Both no flip and mutual spin flip’ and another is ‘Only no flip’. In first process, an electron (with spin s and magnetic moment m) can scatter from the magnetic impurity (with spin S' and magnetic moment m') while flipping its own spin as well as flipping the spin of magnetic impurity with finite probability although there is no certainty for mutual spin flip, there can be another process without any flip. In second process, an electron (with spin s and magnetic moment m) scatters from the magnetic impurity (with spin S' and magnetic moment m') without flipping either its own spin or the spin of magnetic impurity and there is no possibility of spin-flip scattering.

The Bogoliubov-de Gennes (BdG) Hamiltonians for normal metal (NM) and p-wave superconductor (pSc), from Ref. [setiawan] are-


    H_NM=(-ħ^2∂_x^2/2m^*-μ_NM)τ_z,
          H_pSc=(-ħ^2∂_x^2/2m^*-μ_pSc)τ_z-iΔ_pSc∂_xτ_x,


 μ_NM and μ_pSc are the respective chemical potentials, m^* the mass of electron, Δ_pSc≥0 is the p-wave pairing potential, τ_μ=σ_μ⊗ I, with I being 2× 2 identity matrix and σ_μ (μ=x,y,z) are Pauli matrices. For simplicity, we consider ħ=μ_NM=2m^*=1. The energy spectrum are then ε_NM,±(k)=±(k^2-1) and ε_pSc,±(k)=±√((k^2-μ_pSc)^2+(Δ_pSck)^2), respectively. In Fig. 1 (Ref. <cit.>) the energy spectra of the pSc for different values of μ_pSc is plotted. It is seen that the energy spectrum becomes gapless at μ_pSc=0. But for nontopological regime, there is a energy gap |μ_pSc| in the energy spectrum. In this work we only concentrate on the topological regime<cit.>, i.e., μ_pSc>0. For μ_pSc>Δ_pSc^2/2, the positive energy spectrum for pSc shows the characteristic “double-well” BCS structure with minima at ϵ_1=Δ_pSc√(μ_pSc-Δ_pSc^2/4) for k=±√(μ_pSc-Δ_pSc^2/2) and a local maximum at ϵ_2=μ_pSc for k=0, as seen in Fig. 1 of Ref. <cit.>.



  §.§.§ Wavefunctions-

The wave functions for the different regions of our system as shown in Fig. 1 can be written  for a spin-up electron incident at x=-a interface as-


    ψ_NM^I(x)=[ 1; 0; 0; 0 ]e^i(x+a)ϕ_m'^S'+r_ee^↑↑[ 1; 0; 0; 0 ]e^-i(x+a)ϕ_m'^S'+r_ee^↑↓[ 0; 1; 0; 0 ]e^-i(x+a)ϕ_m'+1^S'+r_eh^↑↑[ 0; 0; 0; 1 ]e^i(x+a)ϕ_m'+1^S'+r_eh^↑↓[ 0; 0; 1; 0 ]e^i(x+a)ϕ_m'^S',
    ,


    ψ_NM^II(x)=t_ee^'↑↑[ 1; 0; 0; 0 ]e^i(x+a)ϕ_m'^S'+t_ee^'↑↓[ 0; 1; 0; 0 ]e^i(x+a)ϕ_m'+1^S'+b_ee^↑↑[ 1; 0; 0; 0 ]e^-ixϕ_m'^S'+b_ee^↑↓[ 0; 1; 0; 0 ]e^-ixϕ_m'+1^S'
    +c_eh^↑↑[ 0; 0; 0; 1 ]e^ixϕ_m'+1^S'+c_eh^↑↓[ 0; 0; 1; 0 ]e^ixϕ_m'^S'+a_eh^↑↑[ 0; 0; 0; 1 ]e^-i(x+a)ϕ_m'+1^S'+a_eh^↑↓[ 0; 0; 1; 0 ]e^-i(x+a)ϕ_m'^S',,
 

    ψ_pSc(x)=t_ee^↑↑[ η_-;   0;   0;   1 ]e^ik_-xϕ_m'^S'+t_ee^↑↓[   0; η_-;   1;   0 ]e^ik_-xϕ_m'+1^S'+t_eh^↑↑[ η_+;   0;   0;   1 ]e^ik_+xϕ_m'+1^S'+t_eh^↑↓[   0; η_+;   1;   0 ]e^ik_+xϕ_m'^S',,


where η_±=E+k_±^2-μ_pSc/Δ_pSck_± and ϕ_m'^S' is the eigenspinor of spin flipper: with the S'_z operator acting as- S'_zϕ_m'^S'=m'ϕ_m'^S', m' is the magnetic moment of spin flipper. r_ee^↑↑ and r_ee^↑↓ are the normal reflection amplitudes of an incoming electron with spin up (↑) reflected as an electron with spin up (↑) and spin down (↓) respectively, while r_eh^↑↑ and r_eh^↑↓ are the Andreev reflection amplitudes of an incoming electron with spin up (↑) reflected as a hole with spin up (↑) and spin down (↓) respectively. Similarly, t_ee^↑↑, t_ee^↑↓, t_eh^↑↑, t_eh^↑↓ are the transmission amplitudes into pSc. In Eqs. (<ref>,<ref>) we approximate the wave vector in normal metal by the Fermi wave vector k_F=√(2m^*μ_NM)/ħ=1  (since ħ=μ_NM=2m^*=1) with E<<E_F. Wave vector's k_± in p-wave superconductor are solutions of-

    E^2=(k^2-μ_pSc)^2+(Δ_pSck)^2.

Solutions of Eq. (<ref>) for various values of chemical potential μ_pSc>0 (in the topological regime) with energy E are mentioned in Table I of Ref. <cit.>.



  §.§.§ Boundary conditions
 Boundary conditions at x=-a are- ψ_NM^I(x)=ψ_NM^II(x) (continuity of wavefunction) and, 2i∂_xτ_zψ_NM^I(x)-2i∂_xτ_zψ_NM^II(x)=2iJs⃗.S⃗'⃗τ_zψ_NM^I(x) (discontinuity in first derivative). Boundary conditions at x=0 are- ψ_NM^II(x)=ψ_pSC(x) (continuity of wavefunction) and, (-2i∂_xτ_z+Δ_pSCτ_x)ψ_pSC(x)+2i∂_xτ_zψ_NM^II(x)=-2iZτ_zψ_NM^II(x) (discontinuity in first derivative).

The action of spin-flipper Hamiltonian s⃗.S⃗'⃗ on the wavefunction ψ_NM^I via the boundary condition at x=-a gives for spin up electron-

    s⃗.S⃗'τ_z[ 1; 0; 0; 0 ]ϕ_m'^S'=m'/2[ 1; 0; 0; 0 ]ϕ_m'^S'+F/2[ 0; 1; 0; 0 ]ϕ_m'+1^S',

where F=√((S'-m')(S'+m'+1)) is the spin flip probability of spin flipper. Similarly, the action of exchange operator for spin down electron gives-

    s⃗.S⃗'τ_z[ 0; 1; 0; 0 ]ϕ_m'+1^S'=-m'+1/2[ 0; 1; 0; 0 ]ϕ_m'+1^S'+F/2[ 1; 0; 0; 0 ]ϕ_m'^S'.

Further, the action of exchange operator on spin up holes gives-

    s⃗.S⃗'τ_z[ 0; 0; 1; 0 ]ϕ_m'+1^S'=m'+1/2[ 0; 0; 1; 0 ]ϕ_m'+1^S'-F/2[ 0; 0; 0; 1 ]ϕ_m'^S',

and finally the action of exchange operator on spin down holes gives-

    s⃗.S⃗'τ_z[ 0; 0; 0; 1 ]ϕ_m'^S'=-m'/2[ 0; 0; 0; 1 ]ϕ_m'^S'-F/2[ 0; 0; 1; 0 ]ϕ_m'+1^S'.

Using the above equations and solving the boundary conditions we obtain 16 equations. From these 16 equations we can compute the different scattering probabilities: R_ee^↑↑=|r_ee^↑↑|^2, R_ee^↑↓=|r_ee^↑↓|^2, R_eh^↑↑=|r_eh^↑↑|^2, R_eh^↑↓=|r_eh^↑↓|^2. Similarly, if we consider a spin down electron incident from metallic region I, we get the different scattering probabilities as follows- R_ee^↓↑, R_ee^↓↓, R_eh^↓↑, R_eh^↓↓. 



 §.§ Spin flipper in the vicinity of metal-SOCSW junction


Next, we consider a spinful normal metal (NM)-spin flipper-normal metal (NM)-Insulator (I)-SOCSW junction as shown in Fig. 2. As done before, we model the metal superconductor interface as a δ-like potential barrier with strength Z.



  §.§.§ Hamiltonian

Using the same convention as before, i.e., ħ=2m^*=μ_NM=1, the BdG Hamiltonians for normal metal (NM) and spin orbit coupled superconducting wire (SOCSW) can be written as<cit.>-


    H_NM   =   (-∂_x^2-1)τ_z,
    
    H_SOCSW   =   -∂_x^2τ_z-iβ∂_xτ_zσ_z+B_Zσ_x+Δ_0τ_x,


 with β being strength of spin orbit interaction, B_Z is Zeeman field, and Δ_0≥0- the proximity induced s-wave pairing potential. We also consider uniform electron masses throughout the system and fix chemical potential of the SOCSW and Zeeman interaction in the lead to be zero<cit.>. The positive branches of the energy spectrum<cit.> of SOCSW are given as ε_SOCSW,±(k)=(k^4+β^2k^2+Δ_0^2+B_Z^2±2√(k^4(β^2k^2+B_Z^2)+Δ_0^2 B_Z^2))^1/2. The energy spectrum for SOCSW is shown in Fig. 4(b,c,d) of Ref. <cit.> for nontopological (B_Z<Δ_0), topological (B_Z>Δ_0) and transition regimes (B_Z=Δ_0). There is a gap in the energy spectrum except for B_Z=Δ_0. We focus only on the topological regime B_Z>Δ_0<cit.>. In the limit of strong spin orbit interaction (SOC)  (β≫ B_Z,Δ_0), the energy spectrum of the SOCSW has two branches<cit.>-(i) interior branch, (ii) exterior branch. In Fig. 5(b) of Ref. <cit.>  the energy spectrum for SOCSW is plotted in the limit of strong SOC. The energy spectrum about the minima at k=0 form the interior branches, while the energy spectrum about the minima at k=±β constitute the exterior branches.  For small energies the SOCSW Hamiltonian, Eq. (<ref>) can be linearized about these minima (k=0 and k=±β) by introducing the ansatz for the wavefunction-


    ψ_SOCSW↑(x)≈ R_↑(x)+L_↑(x)e^-iβ x,
          ψ_SOCSW↓(x)≈ R_↓(x)e^iβ x+L_↓(x),


where R_σ(x) and L_σ(x) (with σ=↑,↓) denote the slowly varying right moving and left moving waves respectively. By inserting this ansatz for the superconducting wavefunction in the SOCSW Hamiltonian Eq. (<ref>) and neglecting all the terms involving e^± iβ x (fast oscillating terms), we get the effective Hamiltonian corresponding to the interior and exterior branches, which is given as<cit.>-

    ℋ_SOCSW^(L)=1/2∫ dx  Ψ_SOCSW^(L)(x)^†H_SOCSW^(L)Ψ_SOCSW^(L)(x)

where L=i, e represents the interior and exterior branches respectively. The BdG Hamiltonians for the interior and exterior branches are written as-


    H_SOCSW^(i)=-iβτ_zσ_z∂_x+B_Zσ_x+Δ_0τ_x,
          H_SOCSW^(e)=-iβτ_zσ_z∂_x+Δ_0τ_x.



In addition to Ref. [setiawan], Ref. [kli] too linearizes the SOCSW Hamiltonian because in the limit of strong spin orbit coupling (SOC) the magnetic field and proximity-induced superconductivity are considered as small perturbations. The BdG Hamiltonian used for the superconducting nanowire in Ref. [kli] is given by-

    H_SOCSW^l=-∂_x^2τ_z-iβσ_z∂_x+B_Zσ_xτ_z+Δ_0σ_yτ_y

The linearized BdG Hamiltonians from Eq. (<ref>) for the interior and exterior branches of the superconducting nanowire in the regime of strong SOC as given in Ref. <cit.> is-


    H_SOCSW^l(i)=-iβσ_z∂_x+B_Zσ_xτ_z+Δ_0σ_yτ_y,
          H_SOCSW^l(e)=iβσ_z∂_x+Δ_0σ_yτ_y,


If a Hamiltonian possesses particle-hole symmetry then τ_yσ_yℋ^*(-k)τ_yσ_y=-ℋ(k) ( or, ℋ(k)=-ℋ^*(-k), or ℋ(k)=-τ_xℋ^*(-k)τ_x, or ℋ(k)=-τ_zσ_yℋ^*(-k)σ_yτ_z, or ℋ(k)=-σ_yℋ^*(-k)σ_y, or ℋ(k)=-τ_yℋ^*(-k)τ_y), where k is the wavevector and τ_x=[ 0 0 1 0; 0 0 0 1; 1 0 0 0; 0 1 0 0 ], τ_y=[  0  0 -i  0;  0  0  0 -i;  i  0  0  0;  0  i  0  0 ] and τ_z=[  1  0  0  0;  0  1  0  0;  0  0 -1  0;  0  0  0 -1 ] are the Pauli matrices in particle-hole space, while σ_x=[ 0 1 0 0; 1 0 0 0; 0 0 0 1; 0 0 1 0 ], σ_y=[  0 -i  0  0;  i  0  0  0;  0  0  0 -i;  0  0  i  0 ] and σ_z=[  1  0  0  0;  0 -1  0  0;  0  0  1  0;  0  0  0 -1 ] are Pauli matrices in spin space.

Similarly, if Hamiltonian possesses time reversal symmetry then ℋ(k)=τ_zℋ^*(-k)τ_z ( or, ℋ(k)=τ_yℋ^*(-k)τ_y, or ℋ(k)=τ_xℋ^*(-k)τ_x, or ℋ(k)=σ_yℋ^*(-k)σ_y, or ℋ(k)=ℋ^*(-k)) and finally if Hamiltonian possesses chiral symmetry then ℋ(k)=-τ_zℋ(k)τ_z. 


For 1D SOCSW Hamiltonian only particle-hole symmetry holds. The 1D SOCSW Hamiltonian in Eq. (<ref>) can be written as- 

    H_SOCSW(k)=[  k^2+β k      B_Z      Δ_0        0;      B_Z  k^2-β k        0      Δ_0;      Δ_0        0 -k^2-β k      B_Z;        0      Δ_0      B_Z -k^2+β k ]

 Using Eq. (<ref>) we get-

    τ_yσ_yH_SOCSW^*(-k)τ_yσ_y=      [  0  0  0 -1;  0  0  1  0;  0  1  0  0; -1  0  0  0 ][  k^2-β k      B_Z      Δ_0        0;      B_Z  k^2+β k        0      Δ_0;      Δ_0        0 -k^2+β k      B_Z;        0      Δ_0      B_Z -k^2-β k ][  0  0  0 -1;  0  0  1  0;  0  1  0  0; -1  0  0  0 ]
    
                                               =      [ -k^2-β k     -B_Z     -Δ_0        0;     -B_Z -k^2+β k        0     -Δ_0;     -Δ_0        0  k^2+β k     -B_Z;        0     -Δ_0     -B_Z  k^2-β k ]

From Eqs. (<ref>,<ref>), we see  τ_yσ_yH_SOCSW^*(-k)τ_yσ_y=-H_SOCSW(k). Similarly, we can show that time reversal symmetry does not hold as: H_SOCSW(k)≠τ_zH_SOCSW^*(-k)τ_z, or H_SOCSW(k)≠τ_yH_SOCSW^*(-k)τ_y, or H_SOCSW(k)≠τ_xH_SOCSW^*(-k)τ_x, or H_SOCSW(k)≠σ_yH_SOCSW^*(-k)σ_y, or H_SOCSW(k)≠ H_SOCSW^*(-k). Finally chiral symmetry also does not hold as: H_SOCSW(k)≠-τ_zH_SOCSW(k)τ_z. Thus, the SOCSW Hamiltonian in Eq. (<ref>) only satisfies particle-hole symmetry relation and therefore belongs to symmetry class D.
The SOCSW Hamiltonian Eq. (<ref>) from Ref. [kli] also only satisfies particle-hole symmetry relation discussed above.
Thus, the SOCSW Hamiltonians Eq. (<ref>) used in Ref. [setiawan] and Eq. (<ref>) from Ref. kli are ofcourse Hermitian, as is wont, but unlike the p-wave case (section I) which is time reversal symmetric as well as satisfies particle-hole and chiral symmetries, the SOCSW Hamiltonians satisfy only particle-hole symmetry.
Therefore, according to the symmetry properties of the Hamiltonian, all the Hamiltonians in  Ref. [setiawan] and Ref. [kli] belong to the same symmetry class D<cit.>. The exterior branches of Hamiltonian Eq. (<ref>) in Ref. [setiawan] also satisfies particle-hole symmetry relation τ_yσ_yH_SOCSW^(e)*(-k)τ_yσ_y=-H_SOCSW^(e)(k), and therefore belongs to D class<cit.>. But the exterior branch Hamiltonian Eq. (<ref>) of Ref. [kli] satisfies the time reversal symmetry relation H_SOCSW^l(e)(k)=τ_yH_SOCSW^l(e)*(-k)τ_y as also the particle-hole symmetry relation H_SOCSW^l(e)(k)=-τ_xH_SOCSW^l(e)*(-k)τ_x and thus according to the classification of symmetries of the Hamiltonian it belongs to DIII class<cit.>. The interior branch Hamiltonians Eqs. (<ref>,<ref>) in both Ref. [setiawan] and Ref. [kli] satisfy only particle-hole symmetry relation, therefore belong to the same symmetry class D. Thus, as a whole the linearized Hamiltonians for all cases belong to symmetry class D.



  §.§.§ Wavefunctions
 The wavefunction in the normal metal regions I and II for spin up electron incident with energy E is as mentioned in Eqs. (<ref>,<ref>).  For SOCSW the spinors for the interior and exterior branches are described in terms of the slowly varying left and right moving waves as Ψ_SOCSW^(i)=[    R_↑(x);    L_↓(x);  L_↓^†(x); -R_↑^†(x) ] and Ψ_SOCSW^(e)=[    L_↑(x);    R_↓(x);  R_↓^†(x); -L_↑^†(x) ]. After diagonalizing the Hamiltonians in Eqs. (<ref>,<ref>), we will get the eigenfunctions for both interior and exterior branches.
The wavefunction in the SOCSW is sum of solutions for exterior as well as an interior branch (Eq. (<ref>)) and in the topological regime can be written as<cit.>-


    ψ_SOCSW(x)=t_1^(i)[   -u_-ϕ_m'^S'; -v_-ϕ_m'+1^S';  v_-ϕ_m'+1^S';    u_-ϕ_m'^S' ]e^ik_-^(i)x+t_2^(i)[   u_+ϕ_m'^S'; v_+ϕ_m'+1^S'; v_+ϕ_m'+1^S';   u_+ϕ_m'^S' ]e^ik_+^(i)x+t_1^(e)[   v_0ϕ_m'^S';            0; u_0ϕ_m'+1^S';            0 ]e^i(k_0^(e)-β)x+t_2^(e)[            0; u_0ϕ_m'+1^S';            0;   v_0ϕ_m'^S' ]e^i(k_0^(e)+β)x

                                                                      
where the first two terms on the right-hand side denote contributions from the interior branches, while the other two terms give contributions from the exterior branches. t_(1,2)^(i) and t_(1,2)^(e) are the transmission amplitudes into SOCSW. 
In Eq. (<ref>), for E≥|Δ_λ|, u_λ^2=E+(E^2-Δ_λ^2)^1/2/2E and u_λ^2+v_λ^2=1, while for 0≤ E<|Δ_λ|, u_λ^2=E+i(Δ_λ^2-E^2)^1/2/2|Δ_λ| and u_λ^2+v_λ^2=E/|Δ_λ|, where λ=±,0, and Δ_±=Δ_0± B_Z. The wave vectors in Eq. (<ref>) are k_±^(i)=(E^2-Δ_±^2)^1/2/β for the interior branch, and k_0^(e)=(E^2-Δ_0^2)^1/2/β for the exterior branch. In the limit of strong SOC (β≫ B_Z, Δ_0), we neglect terms proportional to k_+^(i), k_-^(i), k_0^(e)≪1 in our calculation.



  §.§.§ Boundary conditions
 The boundary conditions at x=-a are- ψ_NM^I(x)=ψ_NM^II(x) (continuity of wavefunction) and, 2i∂_xτ_zψ_NM^I(x)-2i∂_xτ_zψ_NM^II(x)=2iJs⃗.S⃗'⃗τ_zψ_NM^I(x) (discontinuity in first derivative). The boundary conditions at x=0 are- ψ_NM^II(x)=ψ_SOCSW(x) (continuity of wavefunction) and, (-2i∂_xτ_z+βτ_zσ_z)ψ_SOCSW(x)+2i∂_xτ_zψ_NM^II(x)=-2iZτ_zψ_NM^II(x) (discontinuity in first derivative). Substituting the wavefunctions in the above boundary conditions we get 16 equations. Solving the 16 equations we get different normal and Andreev reflection probabilities: R_ee^↑↑=|r_ee^↑↑|^2, R_ee^↑↓=|r_ee^↑↓|^2, R_eh^↑↑=|r_eh^↑↑|^2, R_eh^↑↓=|r_eh^↑↓|^2. 
Similarly, if we consider a spin down electron incident from normal metal (region I), we can easily calculate the different reflection probabilities as: R_ee^↓↑=|r_ee^↓↑|^2, R_ee^↓↓=|r_ee^↓↓|^2, R_eh^↓↑=|r_eh^↓↑|^2, R_eh^↓↓=|r_eh^↓↓|^2. We do not repeat them here but in analogy to spin up case the wavefunctions can be easily written and solved for the afore mentioned boundary conditions.



 §.§ Differential charge conductance

Using the well established definitions as in Refs. [jin, kashiwaya], we calculate the net differential charge conductance as-

    G_c=G_c^↑+G_c^↓,  G_c^↑=G_0(1+R_eh^↑↑+R_eh^↑↓-R_ee^↑↑-R_ee^↑↓)
           G_c^↓=G_0(1+R_eh^↓↑+R_eh^↓↓-R_ee^↓↑-R_ee^↓↓)

with, G_0=e^2/h and G_c^↑- the differential charge conductance when spin up electron is incident from region I, while G_c^↓- the differential charge conductance when spin down electron is incident from region I.



§ RESULTS & ANALYSIS
 In Table I, we compare the zero-bias conductance of a NM-spin flipper-NM-I-SOCSW junction and NM-spin flipper-NM-I-pSc junction for both transparent (Z=0) and tunnel (Z=3) regimes. For spin flip case, F≠ 0, i.e., S'≠ m', see Eqs. (<ref>-<ref>), implying there is finite possibility for spin-flipper to flip its own spin while interacting with an electron. On the other hand, for no flip case, F=0, i.e., S'=m'. We take two different values of spin orbit coupling strength β (β=1 and β=2) for SOCSW in second and third column of Table I. We also take two different values of μ_pSc=0.01, 0.001, for pSc in fourth and fifth column of Table I. 
For no flip (F=0) case we see that the normalized zero bias conductance G_c/G_0 is quantized at 2 for NM-spin flipper-NM-I-SOCSW junction, while for NM-spin flipper-NM-I-pSc junction it is quantized at 4 regardless of other parameters like S', m', J, Z, β, μ_pSc, etc. The reason for this is that the Andreev and normal reflection probabilities exactly cancel at zero bias in Eq. (<ref>) for NM-spin flipper-NM-I-SOCSW junction while for NM-spin flipper-NM-I-pSc junction there is perfect Andreev reflection (i.e., normal reflection probabilities vanish) at zero bias and thus from Eq. (<ref>), G_c/G_0 is quantized at 4. In Table I, we have  three non-zero values of spin flip probability (F=1,√(3),3) for transparent junction (Z=0) while two non-zero values (F=1,3) for tunnel junction (Z=3), since the case of F=√(3) is plotted in Fig. 3. 








We notice that for NM-spin flipper-NM-I-SOCSW junction, G_c/G_0 is no longer quantized and depends on the parameters S', m', J, Z, β, etc. However, this is not the case for NM-spin flipper-NM-I-pSc junction, wherein G_c/G_0 is robustly quantized and is independent of spin flip scattering. For NM-spin flipper-NM-I-pSc junction there is perfect Andreev reflection at zero bias even in presence of spin flip scattering (F≠0). Thus, from Eq. (<ref>), zero bias conductance is quantized at 4e^2/h and does not depend on the parameters S', m', J, Z, μ_pSc, etc. For no flip(F=0) case G_c/G_0 is robustly quantized in the topological regime for both cases, see Table I. 
Thus, in a metal-pSc junction the “Majorana states”are not affected by presence of spin flipper, while for metal-SOCSW junction the “Majorana states”are affected by presence of spin flipper. In Fig. 3 we plot the differential charge conductance as a function of energy E for different values of interface barrier strength Z in the topological regime. In Fig. 3(a) for a NM-spin flipper-NM-I-pSc junction we see that G_c/G_0 at E=0 is quantized and independent of Z and spin flip probability (F) of spin-flipper. Thus the topological character of the zero-bias conductance peak is evident, implying the stability of Majorana state in NM-spin flipper-NM-I-pSc junction. The differential conductance for NM-spin flipper-NM-I-pSc junction decreases with increase of Z (transparency of junction). Furthermore, the width of the zero bias conductance peak (ZBCP) also decreases with increasing Z. A discontinuity appears in the differential conductance curve at the gap edge ϵ_1=Δ_pSc√(μ_pSc-Δ_pSc^2/4)=0.007, where Δ_pSc=0.07, μ_pSc=0.01 (topological regime) and μ_pSc>Δ_pSc^2/2 regardless of Z. The reason for this discontinuity is that the wave vectors k_± in p-wave superconductor are different for below (E<ϵ_1) and above (E>ϵ_1) the gap (see Table I of Ref. <cit.>).
In Fig. 3(b) we do the same for NM-spin flipper-NM-I-SOCSW junction. We see that G_c/G_0 at E=0 is not quantized and depends on interface transparency Z. The zero-bias differential conductance (G_c/G_0 at E=0) increases with increase of Z (junction transparency). The differential conductance in NM-SOCSW case exhibits much more complex structure than that for a junction with p wave superconductor. Further, we observe a discontinuity in the differential conductance for SOCSW case at the energy corresponding to energy gaps |Δ_-|(=0.5Δ_0) and Δ_0. The reason for this behavior is that u_λ and v_λ are different for below (0≤ E<|Δ_λ|) and above (E≥|Δ_λ|) the energy gap, where λ=±,0.
Thus, in presence of spin flip scattering the topological character of zero-bias conductance peak seen in case of NM-spin flipper-NM-I-SOCSW junction is affected, suggesting the absence of Majorana states in such junctions.

The reason for quantized zero bias conductance for NM-spin flipper-NM-I-SOCSW junction in absence of spin-flip scattering  is exact cancellation between normal and Andreev reflection probabilities in the conductance formula (Eq. (<ref>)). Differential charge conductance is given as- G_c=e^2/h(2+(R_ee^↑↑+R_ee^↑↓+R_ee^↓↑+R_ee^↓↓)-(R_eh^↑↑+R_eh^↑↓+R_eh^↓↑+R_eh^↓↓)). In no flip process and at zero energy (E=0), 
 R_ee^↑↑+R_ee^↑↓+R_ee^↓↑+R_ee^↓↓=R_eh^↑↑+R_eh^↑↓+R_eh^↓↑+R_eh^↓↓=1. Thus, normal reflection probabilities exactly cancel the Andreev reflection probabilities in conductance formula and this leads to the quantized conductance (G_c=2e^2/h) at zero bias. But, in presence of spin-flip scattering (S'=1/2, m'=-1/2, F=1) and for parameters E=0, β=1 and B_Z=1.5Δ_0 from Eq. (<ref>) we get-

    R_ee^↑↑+R_ee^↑↓+R_ee^↓↑+R_ee^↓↓=   400+264J^2+33J^4+320JZ+144 J^3Z+640Z^2+352J^2Z^2+256JZ^3+256Z^4/400+200J^2+29J^4+320JZ+112J^3Z+640Z^2+288J^2Z^2+256JZ^3+256Z^4,
    
     
    R_eh^↑↑+R_eh^↑↓+R_eh^↓↑+R_eh^↓↓=   400+136J^2+25J^4+320JZ+80 J^3Z+640Z^2+224J^2Z^2+256JZ^3+256Z^4/400+200J^2+29J^4+320JZ+112J^3Z+640Z^2+288J^2Z^2+256JZ^3+256Z^4.

 From Eqs. (<ref>,<ref>), it is evident that- R_ee^↑↑+R_ee^↑↓+R_ee^↓↑+R_ee^↓↓≠ R_eh^↑↑+R_eh^↑↓+R_eh^↓↑+R_eh^↓↓. Thus, there is no exact cancellation at zero bias between normal and Andreev reflection probabilities in presence of spin flip scattering and as a consequence zero bias conductance isn't quantized at 2e^2/h and depends on parameters like J, Z, etc. Remarkably, by tuning parameters we can get a far greater deviation from charge conduction quantization at 2 in case of Normal metal-SOCSW junction with spin flipper. For S=5/2, m'=-5/2, J=3, Z=3, β=12, a=π, Δ_0=0.001, B_Z=1.5Δ_0 we get G_c/G_0=0.25 at E=0, which is much smaller than the quantized value of 2.
We have plotted the normalized conductance (G_c/G_0) as a function of energy E for the above mentioned parameters in Fig. 4 below.

 
Next we try to find out possible reasons for this behavior by looking at the symmetries of scattering matrix and the topological quantum number. 


 §.§ Spin flip scattering and symmetry class

To understand the deviation from topological character due to spin flip scattering in case of NM-SOCSW junction, we analyze the scattering matrix for our system, which is a 4×4 matrix below the gap, given by-

    S_NS=[ s_ee s_eh; s_he s_hh ]

where s_ee=[ r_ee^↑↑ r_ee^↑↓; r_ee^↓↑ r_ee^↓↓ ], s_eh=[ r_eh^↑↓ r_eh^↑↑; r_eh^↓↓ r_eh^↓↑ ], s_he=[ r_he^↓↑ r_he^↓↓; r_he^↑↑ r_he^↑↓ ], s_hh=[ r_hh^↓↓ r_hh^↓↑; r_hh^↑↓ r_hh^↑↑ ] are all 2×2 block matrix. r_e(h)e(h)^σσ' is the normal reflection amplitude of an incoming electron (hole) with spin σ reflected as a 
electron (hole) with spin σ'. Similarly, r_e(h)h(e)^σσ' is the Andreev reflection amplitude of an incoming electron (hole) with spin σ reflected as a hole (electron) with spin σ', σ=↑,↓ and σ'=↑,↓. Following the similar procedure as mentioned before in section II, if we consider a spin up (down) hole incident from normal metal region I, we will get different reflection amplitudes as: r_hh^↑↑ (r_hh^↓↑), r_hh^↑↓ (r_hh^↓↓), r_he^↑↑ (r_he^↓↑), r_he^↑↓ (r_he^↓↓).


  §.§.§ Normal metal-spin flipper-normal metal-p-wave superconductor junction

For no flip process (F=0), and parameters J=1, a=π, Z=1, E=0, μ_pSc=0.01 and Δ_pSc=0.07, the scattering matrix for p-wave superconductor is given as-

    S_NS=[  0  0  0 -i;  0  0 -i  0;  0  i  0  0;  i  0  0  0 ]

S_NS is an unitary matrix with S_NS.S_NS^†=1 and Det S_NS=1. The scattering matrix S_NS and its determinant Det S_NS (Det S_NS=1) do not change with change of parameters like J, Z, a. In Eq. (<ref>) block reflection matrices satisfy the particle-hole symmetry- s_ee=s_hh^* and s_eh=s_he^*. Using this particle-hole symmetry we can choose a basis where all the block matrices have purely real elements. This is called Majorana basis<cit.> in which the scattering matrix is a real orthogonal matrix. The Majorana character of Bogoliubov quasiparticle is hidden in particle-hole basis, but becomes evident in Majorana basis. Thus, to transform the particle-hole basis into Majorana basis of the scattering matrix S_NS, we do a unitary transformation on the scattering matrix S_NS such that the transformed scattering matrix is-

    S=Ω S_NSΩ^†=[  0  0  0 -1;  0  0 -1  0;  0 -1  0  0; -1  0  0  0 ]

where Ω=1/√(2)[  1  0  1  0;  0  1  0  1; -i  0  i  0;  0 -i  0  i ], see Ref. <cit.> for details of Ω.
The matrix S is also an unitary matrix with Det S=Det S_NS=1. The transformed scattering matrix S and its determinant Det S remain unchanged with change of different junction parameters like J, Z, a. Since, the matrix S satisfies S=S^* (particle-hole symmetry), S=S^T (time-reversal symmetry) and S^2=1 (chiral symmetry), thus according to the classification of symmetries of the S-Matrix for the 1D NM-p-wave superconductor junction is in symmetry class BDI<cit.>. In 1D BDI class is topological with the topological quantum number Q for class BDI defined as the number of negative eigenvalues of scattering matrix<cit.>-

    Q=ν(r)


where ν(r) is the number of negative eigenvalues of scattering matrix S. The eigenvalues of S are- -1, -1, 1, 1. Thus, in no flip process topological quantum number for class BDI is 2. In Figs. 5(a) and 5(b) we plot topological quantum number Q as a function of Z and J respectively. Figs. 5 (a) and (b) show that Q is constant at 2 independent of J and Z. Thus, topological quantum number Q for  BDI remains robust against any change of parameters in absence of spin-flip scattering. 

Now in presence of spin-flip scattering(F≠0), the S-Matrix S_NS (Eq. <ref>) for Normal Metal-p-wave superconductor junction remains identical. Similar to the preceding case of no spin-flip scattering, the S-Matrix and its determinant do not change with change of different parameters for spin-flip scattering too. Thus, in presence of spin flip scattering the 1D normal metal-p-wave superconductor junction also belongs to the same symmetry class BDI with topological quantum number Q=2. In Figs. 5(c) and 5(d) we plot topological quantum number Q as a function of Z and J respectively in presence of spin flip scattering. We see that Q remains constant at 2 and does not change with change of Z and J similar to no flip process. The number of Majorana fermion at the edge of the topological superconductor is determined by the topological quantum number. In our work for the case of p-wave junction, topological quantum number for class BDI is the number of negative eigenvalues of the scattering matrix for p-wave superconductor, which is 2 irrespective of spin flip scattering. Thus, the number of Majorana fermion is two at the edge of the superconductor.


  §.§.§ Normal metal-spin flipper-normal metal-SOCSW junction

We do a similar analysis as was done for normal metal-p-wave superconductor junction for the normal metal-SOCSW junction. Our purpose is to understand the loss of ZBCP in presence of spin flip scattering in a NM-spin flipper-NM-I-SOCSW junction. We first analyze the case where there is no spin flip scattering. The scattering matrix, for no flip process (S'=1/2, m'=1/2, F=0), and for parameters J=1, Z=1, E=0, β=0.5, a=π, and B_Z=1.5Δ_0 is-

    S_NS=[  0.065-0.261i -0.543-0.043i -0.130+0.326i  0.696+0.152i;  0.326-0.130i -0.369-0.609i  0.261-0.065i   -0.217-0.5i;    0.217-0.5i  0.261+0.065i -0.369+0.609i -0.326-0.130i;  0.696-0.152i  0.130+0.326i  0.543-0.043i  0.065+0.261i ]

where S_NS is a unitary matrix with S_NS.S_NS^†=1 and Det S_NS=-1. The S-Matrix S_NS however changes with change of parameters unlike the S-Matrix for Normal metal-p-wave superconductor junction. For a=π/2 it  is-

    S_NS=[ -0.133+0.668i  0.594-0.096i -0.214+0.228i  0.255+0.107i; -0.214+0.228i  0.255+0.107i  0.659-0.173i -0.346-0.492i;  0.346-0.492i  0.659+0.173i  0.255-0.107i  0.214+0.228i;  0.255-0.107i  0.214+0.228i -0.594-0.096i -0.133-0.668i ]

where the other parameters are same as for Eq. (<ref>). As shown in Eqs. (<ref>,<ref>), the S-Matrix for a 1D  normal metal-SOCSW junction in absence of spin-flip scattering changes with change in junction length a, while its determinant (Det S_NS=-1) does not depend on a, in both cases a=π and a=π/2, Det S_NS=-1. We have checked that the S-Matrix S_NS changes with other parameters like J, Z, although its determinant remains constant at -1 (Det S_NS=-1). Similarly, as mentioned before, when we transform the particle-hole basis of the S-Matrix into Majorana basis by doing an unitary transformation on S_NS, the transformed S-Matrix is-

    S=[ -0.109+0.087i  0.043+0.022i  0.848+0.391i       -0.326i;  0.913-0.196i -0.196-0.261i  0.087+0.109i  0.022-0.043i; -0.022-0.043i  0.087-0.109i -0.196+0.261i -0.913-0.196i;        0.326i -0.848+0.391i -0.043+0.022i -0.109-0.087i ]
    
    
       
    S=[  0.127+0.148i  0.861+0.206i -0.028+0.087i  0.129+0.392i;  0.053-0.074i -0.005-0.412i -0.195-0.012i -0.747+0.474i;  0.747+0.474i -0.195+0.012i -0.005+0.412i -0.053-0.074i;  0.129-0.392i  0.028+0.087i -0.861+0.206i  0.127-0.148i ]

S is again an unitary matrix with Det S=Det S_NS=-1. Similar to S_NS, the transformed S-Matrix S changes with change of parameters like J, Z, a, but its determinant remains unchanged. We verify that the matrix S only satisfies particle-hole symmetry relation τ_yσ_yS^*τ_yσ_y=S, thus the S-Matrix for no spin-flip scattering in case of a 1D normal metal-SOCSW junction belongs to symmetry class D.

In 1D, symmetry class D is topological with the topological quantum number<cit.> for class D defined as the determinant of the scattering matrix<cit.>-

    Q'=(S)

In Figs. 6(a) and (b), topological quantum number Q' is plotted as a function of Z and J respectively. We see that the Q' is constant at -1 and does not change with change in Z and J. Thus, topological quantum number Q' for class D remains robust against any change of parameters in absence of spin-flip scattering.


In presence of spin-flip scattering (S'=1/2, m'=-1/2, F=1), and parameters J=1, Z=1, E=0, β=0.5 and B_Z=1.5Δ_0, the S-Matrix for a 1D normal metal-SOCSW junction is given as-

    S_NS=[ -0.055+0.435i  0.428+0.199i -0.109+0.207i  0.698+0.207i; -0.436+0.492i -0.055+0.435i  0.184-0.039i -0.232-0.536i; -0.232-0.536i  0.698+0.207i  0.151-0.301i -0.022-0.123i;  0.184-0.039i -0.109+0.207i -0.767-0.456i  0.151-0.301i ]
    
    
       
    S_NS=[ -0.166-0.525i -0.313+0.121i  0.249+0.299i  0.656+0.037i;  0.546+0.028i -0.166-0.525i  0.289-0.140i  0.161-0.518i;  0.161-0.518i  0.656+0.037i -0.158+0.343i -0.145-0.332i;  0.289-0.140i  0.249+0.299i  0.733-0.254i -0.158+0.343i ]

where S_NS is a unitary matrix with S_NS.S_NS^†=1. The determinant of the S-Matrix for a=π/2 and a=π are e^-i2.76667 and e^i2.92505 respectively with absolute value 1. In contrast to no flip case, both the S-Matrix (S_NS) and its determinant (Det S_NS) change with change of parameters like J, Z. Similarly, as discussed before, we do a unitary transformation on S_NS to transform its basis from particle-hole to Majorana. In Majorana basis the transformed scattering matrix is given as-

    S=[ -0.123-0.098i  0.901+0.245i  0.003-0.164i -0.161+0.225i; -0.418-0.021i -0.123-0.098i -0.474+0.166i  -0.73-0.041i;  0.739+0.041i  0.161-0.225i  0.219+0.231i -0.495-0.168i;  0.474-0.166i -0.003+0.164i -0.785+0.057i  0.219+0.231i ]
    
    
       
    S=[  0.043-0.200i  0.427-0.068i  0.843-0.048i -0.227-0.084i;  0.929-0.253i  0.043-0.200i -0.141-0.093i  0.026+0.041i; -0.026-0.041i  0.227+0.084i -0.368+0.018i -0.885-0.143i;  0.141+0.093i -0.843+0.048i  0.351+0.027i -0.368+0.018i ]

where S is an unitary matrix with Det S=Det S_NS. Similar to S_NS, its determinant Det S_NS depends on the various junction parameters. We have checked that the S-Matrix S does not satisfy any symmetry relation and therefore belongs to class A in presence of spin flip scattering.

In Figs. 6(c) and 6(d) we plot both absolute value and argument of the complex Q' (determinant of S) as a function of Z and J respectively for a 1D normal metal-SOCSW junction(Fig. 2). We notice that Abs(Q') remains constant at 1, while the argument of the complex Q' (Arg(Det S)) changes with J and Z.



§ CONCLUSION

To conclude, we have shown that zero bias quantized conductance Majorana peaks remain unaffected, in presence of spin flipper, at metal-p-wave superconductor interface while zero bias conductance peak at metal-SOCSW interface loses its quantization in presence of spin flip scattering. Further, we find that in presence of spin flip scattering a normal metal-p-wave superconductor junction belongs to the symmetry class BDI and topological quantum number for class BDI does not change with change of parameters. But, for normal-metal-SOCSW junction the topological quantum number for class D, i.e., the determinant of the scattering matrix is complex and Det S changes with change of parameters due to spin flip scattering. The S-Matrix for a  1D normal metal-SOCSW junction satisfies only particle-hole symmetry relation and belongs to topological symmetry class D in absence of spin flip scattering, while in presence of spin flip scattering the S-Matrix does not satisfy any symmetry relation and thus belongs to the non-topological class A. The reason for ZBCP in a 1D normal metal-SOCSW junction in absence of spin-flip scattering is the exact cancellation at zero bias of the normal and Andreev reflection probabilities. However, in presence of spin flip scattering in  a 1D normal metal-SOCSW junction the exact cancellation of normal and Andreev reflection probabilities is no longer valid and as a consequence the ZBCP loses it quantization.


§ AUTHORS' CONTRIBUTIONS

C.B. conceived the proposal, S.P. did the calculations on the advice of C.B., C.B. and S.P. analyzed the results and wrote the paper. Both authors reviewed the manuscript.


§ FUNDING

This work was supported by the grant “Non-local correlations in nanoscale systems: Role of decoherence, interactions, disorder and pairing symmetry” from SCIENCE & ENGINEERING RESEARCH BOARD, New Delhi, Government of India, Grant No.  EMR/20l5/001836,
Principal Investigator: Dr. Colin Benjamin, National Institute of Science Education and Research, Bhubaneswar, India.


§ DATA AVAILABILITY

Data available on request from the
authors. The data that support the findings of this study are available from the corresponding author upon reasonable request.
99
BenBeenakker C. W. J 2013 Search for Majorana Fermions in Superconductors Annu. Rev. Con. Mat. Phys. 4 113.
kitKitaev A. Y 2001 Unpaired Majorana fermions in quantum wires Phys. Usp. 44 131.
colinBenjamin C, Pachos J. K 2010 Detecting Majorana bound states Phys. Rev. B 81 085101.
mazMazziotti M. V, Scopigno N, Grilli M and Caprara S 2018 Majorana Fermions in One-Dimensional Structures at LaAlO_3/SrTiO_3 Oxide Interfaces, Condens. Matter. 3(4) 37.
DengDeng M. T, Yu C. L, Huang G. Y, Larsson M, Caroff P and Xu H. Q 2012 Anomalous Zero-Bias Conductance Peak in a Nb-InSb Nanowire-Nb Hybrid Device Nano Lett. 12 6414.
DasDas A, Ronen Y, Most Y, Oreg Y, Heiblum M and Shtrikman H 2012 Zero-bias peaks and splitting in an Al-InAs nanowire topological superconductor as a signature of Majorana fermions Nat. Phys. 8 887.
MouMourik V, Zuo K, Frolov S. M, Plissard S. R, Bakkers E. P. A. M and Kouwenhoven L. P 2012 Signatures of Majorana Fermions in Hybrid Superconductor-Semiconductor Nanowire Devices Science 336 1003.
xuXu Jin-Peng et al 2015 Experimental Detection of a Majorana Mode in the core of a Magnetic Vortex inside a Topological Insulator-Superconductor Bi_2Te_3/NbSe_2 Heterostructure Phys. Rev. Lett. 114 017001.
perPerge S. N et al 2014 Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor Science 346 602.
jeonJeon S et al 2017 Distinguishing a Majorana zero mode using spin-resolved measurements Science 358 772-776.
valVal'kov V. V and Aksenov S. V 2018 Spin-polarized-current switching mediated by Majorana bound states J. Magn. Magn. Mater 465 88-92.
setiawanSetiawan F, Brydon P. M. R, Sau J. D and DasSarma S 2015 Conductance spectroscopy of topological superconductor wire junctions Phys. Rev. B 91 214513.
hasFulga I. C, Hassler F and Akhmerov A. R 2012 Scattering theory of topological insulators and superconductors Phys. Rev. B 85 165409; Beenakker C. W. J 2015 Random-matrix theory of Majorana fermions and topological superconductors Rev. Mod. Phys. 87 1037.
sauTewari S and Sau J. D 2012 Topological Invariants for Spin-Orbit Coupled Superconductor Nanowires Phys. Rev. Lett. 109 150408.
fulgaFulga I. C, Hassler F and Akhmerov A. R and Beenakker C. W. J 2011 Scattering formula for the topological quantum number of a disordered multimode wire Phys. Rev. B 83 155429.
kliKlinovaja J and Loss D 2012 Composite Majorana fermion wave functions in nanowires Phys. Rev. B 86 085408.
BTKBlonder G. E, Tinkham M and Klapwijk T. M 1982 Transition from metallic to tunneling regimes in superconducting microconstrictions: Excess current, charge imbalance, and supercurrent conversion Phys. Rev. B 25 4515.
chengCheng M and Lutchyn R. M 2012 Josephson current through a superconductor/semiconductor-nanowire/superconductor junction: Effects of strong spin-orbit coupling and Zeeman splitting Phys. Rev. B 86 134522.
AJPde Menezes O. L. T and Helman J. S 1985 Spin flip enhancement at resonant transmission Am. J. Phys 53 1100.
hurHur K. L 2015 Quantum dots and the Kondo effect Nature 526 203-204.
tanaTanaka Y and Kashiwaya S 2004 Anomalous charge transport in triplet superconductor junctions Phys. Rev. B 70 012507.
samuLara-Avila Samuel et al 2015 Influence of Impurity Spin Dynamics on Quantum Transport in Epitaxial Graphene Phys. Rev. Lett. 115 106602.
lnolnoshita T 1998 Kondo Effect in Quantum Dots Science 281 526-527.
yokoyYokoyama T et al 2005 Influence of magnetic impurities on charge transport in diffusive-normal-metal/superconductor junctions Phys. Rev. B 71 094506.
gueGuerci D and Nava A 2019 Probing Majorana edge states by measuring transport through an interacting magnetic impurity arXiv:1907.06444.
LiuLiu H. D and Yi X. X 2011 Geometric phases in a scattering process Phys. Rev. A 84 022114.
MaruCordourier-Maruri G, Omar Y, de Coss R and Bose S 2014 Graphene-enabled low-control quantum gates between static and mobile spins Phys. Rev. B 89 075426.
FCCiccarello F, Palma G. M and Zarcone M 2007 Entanglement-induced electron coherence in a mesoscopic ring with two magnetic impurities Phys. Rev. B 75 205415.
ysrPal S and Benjamin C 2018 Yu-Shiba-Rusinov bound states induced by a spin flipper in the vicinity of a s-wave superconductor Sci Rep 8 11949.
tewSau J. D, Tewari S, Lutchyn R. M, Stanescu T. D and DasSarma S 2010 Non-Abelian quantum order in spin-orbit-coupled semiconductors: Search for topological Majorana particles in solid-state systems Phys. Rev. B 82 214509.
jinCheng Q and Jin B 2013 Quantum transport in normal-metal/ferromagnet/spin-triplet superconductor junctions Physica B 426 42.
kashiwayaKashiwaya S, Tanaka Y, Yoshida N and Beasley M. R 1999 Spin current in ferromagnet-insulator-superconductor junctions Phys. Rev. B 60 3572.
ARAkhmerov A. R, Dahlhaus J. P, Hassler F, Wimmer M and Bennakker C. W. J 2011 Quantized Conductance at the Majorana Phase Transition in a Disordered Superconducting Wire Phys. Rev. Lett. 106 057001.
ryuRyu S Periodic table of topological insulators and superconductors Online: https://www.youtube.com/watch?v=cKzUuQyZjFo. 

