Detection of Majorana Bound States by Sign Change of the Tunnel Magnetoresistance in a Quantum Dot Coupled to Ferromagnetic Electrodes

We study quantum transport and tunneling magnetoresistance (TMR) through an interacting quantum dot (QD) attached to ferromagnetic electrodes in the Coulomb blockade regime. The QD is also side-coupled to a superconductor nanowire hosting Majarana bound states (MBSs). It is found that when the electrodes' magnetic moments are arranged in antiparallel, the current's intensity will be enhanced to be larger than that of the parallel configuration by the hybridization between the QD and the MBSs. This change in the current induces anomalous negative TMR unique to the existence of MBSs, providing an efficient detection way of the MBSs. The negative TMR is weakened by the overlap between the two bound states but obviously enhanced by the left-right asymmetry between the QD and the electrodes. We also find that the TMR value changes non-monotonously with the spin polarization of the electrodes. Our results may find real use in energy saving spintronic devices and quantum information processing.


INTRODUCTION
Tunneling of spin-polarized electrons in nano-and micro-scale devices are at the core of quantum information data processings [1,2] and future spintronics [3]. Spin-dependent current is conveniently generated in a system with ferromagnetic electrodes, and its intensity changes when the magnetic configuration of the device changes from parallel to antiparallel alignment [4]. Generally, the electronic current is larger in parallel configuration when transport occurs between the majority-majority and minority-minority spin bands, than in the antiparallel configuration, where electrons tunnel between majority and minority spin bands. The change in the currents is measured by the quantity of tunnel magnetoresistance (TMR) defined as TMR = (J P − J AP )/(J P + J AP ), where J P/AP is the electronic current in parallel (antiparallel) configuration. Usually, the TMR value is positive as J P > J AP [4]. The TMR effect has been extensively studied in various systems, including the relatively simple planar junctions [5], and in more complex mesoscopic doublebarrier junctions [6], junctions with nano-scale granular systems [7,8], as well as single-electron transistors [9], etc. When the central region between the electrodes is sufficiently small that the charging energy is larger than the thermal energy k B T, some interesting effects due to quantization of energy levels becomes visible, such as the Coulomb blockade of electric current characterized by the Coulomb steps in the current-voltage curves [10][11][12][13]. Especially, behaviors of the TMR have been investigated both theoretically and experimentally in systems composing of ferromagnetic electrodes and quantum dots (QDs) [14][15][16][17][18][19], which is a nanoscale zero-dimensional structure [13]. This is because such a kind of device is ideal candidate to study the fundamental interactions between spins and charges [3,13], and is being considered for applications in energy saving devices as well as for quantum computing [1,2].
Recently, there is much interest in preparation and detection of Majorana fermions, which is a kind of elementary particle being of its own antiparticle, due to its promising applications in quantum computation free from decoherence [20][21][22]. It has been be prepared in the edges of topological superconductors as a pair of zero-energy Majorana bound states (MBSs) that obeys non-Abelian statistics enabling topologically protected quantum computation and high-efficiency spintronic devices [20,22]. The MBSs has also be prepared in some other systems, including the non-centrosymmetric superconductors [23], three-or two-dimensional topological insulators coupled to superconductors [24], electrostatic defects in topological superconductors [25], p-wave superconductors [26], the semiconducting [27], or ferromagnetic [28] nanowires with native strong spin-orbit interaction proximitization to a conventional s-wave superconductors, and Josephson junctions [29], etc. The detection of MBSs is quite challenging because the unique chargeless and zero-energy properties, and then most of conventional detection techniques for elementary particles failed. For example, the MBSs formed in the nanowire induces a zerobias anomaly in the differential conductance when a bias voltage is applied [30], which has been demonstrated in experiments. However, it can also arise from some other mechanisms, such as the quantum interference and Kondo effects [31]. Anomalous change of the thermoelectric effect in systems composing of a QD side-coupled to topological superconductor hosting a Majorana edge state has also been found and proposed as an alternative detection means in recent years [32][33][34]. The hybridization between QD and MBSs breaks the particle-hole symmetry of the system and leads to large enhancement and sign reversion of the thermopower. These, unfortunately, can also not been solely attributed to the existence of the MBSs.
In recent experiments, semiconductor spacers of InAs QD has been inserted in between nickel or cobalt leads [15,35]. The dot's size and discretized energy levels can be fully adjusted by the length and thickness of the spacers or the gate voltages. Moreover, the spin polarization of the current injected from the ferromagnetic electrodes and TMR can also be electrically adjusted by a gate nearby the QD. Some new characteristics, such as the anomalies of the TMR caused by the intradot Coulomb repulsion energy in the QD, were explained in the subsequent theoretical work [36]. In the present manuscript, we study the properties of current-voltage and TMR in a QD sandwiched between the left and right ferromagnetic electrodes and side-coupled to a semiconducting nanowire hosting a pair of MBSs at the opposite ends of the wire (see the schematic plot of the system in Figure 1). We assume that the MBSs are coupled only to spin-up electrons in the QD under a strong magnetic field in experiments. Our calculation results show that due to the presence of the MBSs, the current of the antiparallel configuration can be larger than that of the parallel one, inducing a negative TMR. Such an effect is a good evidence of the existence of the Majorana fermions. The negative TMR will be weakened by the coupling between the two MBSs but obviously enhanced by the left-right asymmetry between the QD and the electrodes. Moreover, we find that the TMR strength depends non-monotonously on the spin polarization of the electrodes.

MODEL AND METHODS
The Hamiltonian of the QD coupled to MBSs and ferromagnetic electrodes takes the following form [18,32,37,38] where c † kβσ (c kβσ ) creates (annihilates) an electron of momentum k, energy ε kβσ and spin σ =↑, ↓ in the ferromagnetic electrode β = L, R. For the QD, d † σ (d σ ) is the creation (annihilation) operator of an electron having energy level ε d , spin-σ and intradot Coulomb interaction U, which has usually been neglected in previous work [32]. In experiments, ε d and U are tunable by external gate voltages. The coupling strength between the QD and the electrodes is described by T kβσ , which is spin-dependent due to the ferromagnetism on the electrodes. The last term H MBSs in Equation (1) stands for the zero-energy MBSs located on the opposite ends of the semiconducting nanowire and their coupling to the QD [39,40]: in which δ M is the overlap strength between the two MBSs with operator satisfying both η j = η † j (j = 1, 2) and {η i , η j } = δ i,j . The hopping amplitude between MBSs and spin-↑ electrons in the QD is accounted by λ j . Following previous work [39,40], we write η j in terms of the regular fermionic operators f as

and then H MBSs becomes
Within the standard Keldysh Green's function technique, the spin-dependent electric current is obtained as [37,38] where e is the electron charge,h the reduced Planck's constant, can be expressed with the help of the retarded Green's function G r σ (ε) as [37,38] T where is the spin-dependent line-width function. Ferromagnetism on the electrodes can be expressed by defining a spin-polarization parameter P β as In the present paper, we consider the parallel (P L = P R ) and antiparallel (P L = −P R ) configurations of the two electrodes. In both of the two cases, we have Ŵ β By applying the equation of motion method, the retarded Green's function in Equation (5) is obtained as (we have truncated the higher-order Green's functions by following reference [39,40], in which the simultaneous tunneling of the electron of opposite spin has been neglected): where the self-energies considering the MBSs are given by 0 andε ± = ε ± ε d + i(Ŵ L σ + Ŵ R σ )/2. The occupation number in Equation (6) is calculated self-consistently from To our knowledge, there are three kinds of schemes that are mainly used for studying the transport phenomena [41][42][43]. The first kind is the relatively rough Hartree-Fock-like scheme used in reference [41,42]; the second kind is the Hubbard-I scheme adopted in the present work as well as other reference [39,40] that can deal with relatively higher-order transport processes; the third kind of scheme is the one that can account for multiple transitions from contacts to the central site, such as the subtle Kondo effect. Detailed introduction of these schemes can be found in of reference [43]. In the Coulomb blockade regime, these three schemes give essentially the same results, see for example reference [43]. From theoretical point of view, here we adopt the Hubbard-I scheme for the sake of sophistication, and the procedure to evaluate the QD Green's function will be relatively simplified.

RESULTS AND DISCUSSION
In this section, we present our numerical results for the spindependent current, TMR and differential conductance varying as functions of the bias voltage. The intradot Coulomb interaction U = 1 is chosen as the energy unit with symmetric bias voltage V as µ L = −µ R = eV/2. Figure 2 shows the spinpolarized current, differential conductance and the current's spin-polarization p = (J ↑ − J ↓ )/(J ↑ + J ↓ ) for parallel and antiparallel configurations, respectively. The currents exhibit typical Coulomb blockade effect as shown in Figures 2A,B,D,E [3,10,13]. The reason is that by applying a bias voltage between the electrodes, the QD can be charged by single electrons when the energy levels ε d are within the transport window. The current then is increased. With an increase from n to n + 1 electrons on the QD, an additional potential of e/C is induced with C denoting the capacitance of the device [10,13]. The transfer of an electron through the QD is inhibited when the potential of the charged dot is larger than the applied bias voltage. Now the electrons are blockaded on the QD and the current will not increase with increasing bias voltage. With further increased bias voltage, electrons may tunnel through the dot level of ε d + U again when the potential is overcome. In this way, the electrons are added on the dot one by one and the current shows Coulomb staircases. The Coulomb blockade effect can be clearly seen from the differential conductance, which develops two sharp peaks around ε d and ε d + U as given in the figure.
In the parallel configuration, the current intensity of J ↑ is obviously larger than that of J ↓ in the absence of QD-MBSs. The reason is that the spin-up electrons will enter and leave the QD faster than the spin-down ones because of the ferromagnetism on the electrons, i.e., Ŵ L ↑ = Ŵ R ↑ > Ŵ L ↓ = Ŵ R ↓ . In the antiparallel configuration, however, the current intensities of both the two spin components are the same because the transmission coefficient is spin-independent in the absence of hybridization between the QD and MBSs (λ 1 = λ 2 = 0), which is not shown in the figure. Turning on the coupling between the QD and MBSs, the intensity of spin-up currents in both parallel and antiparallel configuration are enhanced as shown in Figures 2A,D, whereas that of the spin-down electrons are unchanged as the MBSs is coupled only to the spin-up electrons. Moreover, the enhancement of the intensity of J ↑ in the antiparallel is more obvious than that of the parallel one, which results in negative TMR as will be shown in the following. For the differential conductance, an obvious zero bias anomaly (ZBA) emerges in addition to the increase in the intensity [39,40]. The ZBA, which has been theoretically predicted and experimentally observed, originates from the QD-MBSs coupling and serves as a signal of the existence of the Majornana fermions [39,40]. In the absence of bias voltage and Coulomb interaction, the conductance is 0.5 e 2 /h by the coupling between the QD and the MBSs [44][45][46]. Figures 2C,F shows that the current's spin polarization is also enhanced in the presence of the MBSs in both parallel and antiparallel configurations. This is because the intensity of the spin-up current is increased by the MBSs whereas that of the spin-down one remains unchanged. Moreover, the spin polarization of the current changes more drastically near the dot level of ε d as compared to ε d + U due to the spin blockade effect.
As demonstrated by previous work, the influence of the MBSs on the current is rather weak since a pair of Majarana fermions are charge neutral [20][21][22]. This can be seen from the total current J = J ↑ + J ↓ in Figures 3A,B. The ZBA in the differential conductance, which is measured by the electrical tunnel spectroscopy, then is experimentally interesting (indicated by the arrows in the figure). To enlarge the signal of the MBSs in the current, we present the TMR in Figure 3C. For weak QD-MBSs hybridization strength λ = 0.01U, the TMR is positive as indicated by the solid line. It has two plateaus in the bias regimes of eV < ε d and eV > ε d + U, in which the current remains almost unchanged [3]. In the Coulomb blockade regime (ε d < eV < ε d + U), the intensity of TMR increases with FIGURE 3 | Total current, non-linear differential conductance for the magnetic moments' configurations of parallel (A) and antiparallel (B), and the corresponding TMR (C). Other parameters are as in Figure 2.
increasing bias voltage, which is quite different from the behavior of the current. This is because the current from the parallel configuration originates from transport between the majoritymajority and minority-minority spin bands, and that in the antiparallel configuration between majority and minority spin bands. In the absence of the QD-MBSs coupling, the former transport process is rather easier than the latter case. When the electrons on the QD is coupled to the MBSs, however, there is more transport channels in the antiparallel configuration and then the current intensity is increased accordingly. With increasing bias voltage, the leakage current in the parallel configuration then will be larger than that of the antiparallel one, resulting in increased TMR. Around the two transport channel ε d and ε d + U, the TMR develops a dip and peak, respectively. The reason is that the current of the antiparallel configuration originating from transport through ε d (ε d + U) is enhanced (suppressed). With increased QD-MBSs coupling strength, the strength of the TMR is weakened as shown by the dashed and dotted lines in Figure 3C. For sufficiently large λ, the TMR around transport channel ε d is negative (see the dotted line), and then serves as a detection method for the existence of the Majorana fermions in the superconductor nanowire side-coupled to the QD. We emphasize that this phenomenon originates from the fact that the current intensity of the antiparallel configuration is more sensitive to the MBSs as compared to the parallel one.
In Figure 4 we present the TMR varying with respective to the bias voltage for different values of the spin-polarization of the electrodes (p = P L = P R ) and overlap amplitude between the two MBSs δ M . For p = 0, the TMR is zero as indicated by the solid line in Figure 4A. With increasing spin polarization, the line-shape of the TMR essentially resembles those in Figure 3C. The strength of TMR changes with respective to the spin polarization in a non-monotonic way. Figure 4A shows that the TMR is negative around dot level ε d for p <= 0.6, and is positive for P = 0.8. The TMR in the inset of Figure 4A has a minimal value in the regime of 0.5 < p < 0.6 for the indicated dot level and bias voltage. This is because that with increasing ferromagnetism on the electrodes p, the transport processes will be more and more dominated by only one spin-component electrons, and then the difference between the two kinds of the current becomes less prominent. In the extreme case of p ∼ 1 (half-metallic), electron transport in the antiparallel is totally blockaded and the TMR equals to 1. Under this condition, the signal of the MBSs disappears. Figure 4B presents the function of δ M on the TMR for fixed spin polarization of the electrodes p = 0.6. When the two states are decoupled from each other (δ M = 0), an obvious negative TMR emerges around the dot level ε d . With increasing δ M , we find that the strength of the TMR is enhanced, especially in the Coulomb blockade regime. For strong enough overlap between the two states δ M = 0.05U, the TMR is positive in the whole bias voltage regime. The reason is that when the two MBS are coupled, their functions on the TMR are canceled out and the abnormal increase of the current in the antiparallel configuration induced by the MBSs disappears. In Figure 5, we present the total currents in both of the two configurations and the TMR varying with the bias voltage for different values of left-right asymmetry Ŵ R / Ŵ L . First of all, the currents' strength is monotonously enhanced by increasing Ŵ R / Ŵ L . This is because that for fixed ingoing tunneling rate Ŵ L , the electrons will tunnel through the dot faster with increased outgoing tunneling rate Ŵ R . Obviously, this holds true for both of the two spin component electrons regardless of the magnetic configuration. The TMR in Figure 5C, however, changes in a non-monotonic way with varying Ŵ R / Ŵ L . In the regimes of eV < ε d and eV > ε d +U, the TMR is positive and its strength increases with increasing Ŵ R / Ŵ L . This is the usual case in the absence of the QD-MBSs. In the Coulomb blockade regime ε d < eV < ε d + U, however, the value of TMR decreases with increasing left-right asymmetry. As a consequence of it, the negative TMR becomes more prominent and becomes a good means to detect the existence of the MBSs. As mentioned above, this is induced by that fact that the QD-MBSs coupling strengthes the current intensity of the antiparallel configuration as compared to the parallel one, which is an unique function of the MBSs on the electron transport. Finally, we briefly discuss the experimental realization of the present devices. The nanowire hosting the MBSs can be fabricated with InAs grown by molecular beam epitaxy with several nanometers of epitaxial Al layer. It has been experimentally proven that an hard superconducting gap can be induced on such a kind of nanowires by applying a critical magnetic field exceeding 2T along the wire axis. A QD is formed in the bare InAs segment at the end of the wire due to density of state gradients at the edges of the Al shell.

SUMMARY
In summary, spin-dependent current and TMR in a QD sandwiched between two ferromagnetic leads and side-coupled to a pair of MBSs formed at the opposite ends of a superconductor nanowire is investigated within the nonequilibrium Green's technique. An unique negative TMR induced by the hybridization between the QD and the MBSs is found, which serves as a detection means of the Majorana fermions. This negative TMR is more likely to emerge in longer nanowire in which the two MBSs are well-separated from each other and the overlap between them is weak. By increasing the leftright asymmetry of the coupling strength between the QD and electrodes, the negative TMR becomes more obvious. It is also found that the intensity of the TMR depends on the spinpolarization of the electrodes in a non-monotonic way and is positive for large spin-polarization regardless of the existence of the MBSs.

DATA AVAILABILITY STATEMENT
All datasets presented in this study are included in the article/supplementary material.

AUTHOR CONTRIBUTIONS
W-GM and L-WT contributed the ideas equally and performed the numerical calculations. L-WT derived the formulae in the paper and wrote the original manuscript.