Reveal the non-local coherent nature from a dissipative Majorana teleportation

TL;DR

This paper proposes a transport scheme to observe the non-local coherent nature of Majorana devices, showing how dissipation influences conductance and provides a benchmark for validating Majorana qubits.

Contribution

It introduces a dissipative environment-based transport scheme that reveals the non-local coherence of Majorana devices through distinctive conductance behavior.

Findings

Conductance increases with decreasing temperature in Majorana teleportation.

Dissipation suppresses incoherent signals, highlighting coherent transport.

A clear conductance scaling behavior serves as a benchmark for Majorana coherence.

Abstract

The non-local coherent nature of the Majorana devices is one of the key factors for realizing decoherence-free topological qubits. Direct observation of this coherent nature could provide a first-step benchmarking scheme to validate Majorana qubit quality. We propose a simple transport scheme with a Majorana island device along with a dissipative environment in the electrodes. We found that the dissipative environment renormalizes the quantum transport in significant different ways: As reducing temperature, while the conductance for Majorana coherent teleportation increases, all other incoherent signals are strongly suppressed due to dissipation. This special conductance scaling behavior is a clear benchmark to reveal the non-local coherent nature of Majorana devices.

Figures

Equations

$$H_{\text{wire}}=\sum _{\alpha }{ϵ}_{\alpha }{\gamma }_{\alpha }^{†}{\gamma }_{\alpha }+U_N.$$

$$H_{\text{tot}}=H_{\text{leads}}+H_{\text{wire}}+H_{\text{T}}+H_{\text{env}},$$

$$H_{\text{T}}=\sum _{j,\sigma }{t}_j{d}_{\sigma }^{†}\left(r_j\right)c_{\sigma }\left(r_j\right)e^{-i{\phi }_j}+\text{H.c.},$$

$$G=$$

$$\times f\left({\xi }_p\right)\delta \left({ϵ}_{\alpha }-{\xi }_p+U_{N+1}-U_N+E\right)P\left(E\right),$$

$$P\left(E\right)$$

$$\displaystyle\begin{array}[]{cc}\propto T^{\frac{r}{2}-1}e^{\frac{E}{2T}}\left|\mathit{\Gamma}\left(\frac{r}{4}+i\frac{E}{2\pi T}\right)\right|^{2}&,\text{for}\;T\ll\omega_{R}\end{array},$$

$$G_{2e}=$$

$$\times f\left({\xi }_{p_2}\right)\delta \left(E+U_{N+2}-U_N-{\xi }_{p_1}-{\xi }_{p_2}\right)P\left(E\right),$$

$$I$$

$$\frac{\partial W\left(\left\{n_{dot}\right\}\right)}{\partial t}=$$

$$+\sum _{\alpha }W\left(n_1,\dots ,n_{\alpha -1},1,n_{\alpha +1},\dots \right){\delta }_{n_{\alpha },0}\left({\Gamma }_{\alpha \to 0}^l+{\Gamma }_{\alpha \to 0}^r\right)$$

$$+\sum _{\alpha }W\left(n_1,\dots ,n_{\alpha -1},0,n_{\alpha +1},\dots \right){\delta }_{n_{\alpha },1}\left({\Gamma }_{0\to \alpha }^l+{\Gamma }_{0\to \alpha }^r\right).$$

$${\Gamma }_{0\to \alpha }^j$$

$${\Gamma }_{\alpha \to 0}^j$$

$$P\left(E\right)=\frac{{\left(2\pi \right)}^{\frac{r}{2}-2}}{{\omega }_R\Gamma (\frac{r}{2})}{\left(\frac{T}{{\omega }_R}\right)}^{\frac{r}{2}-1}{e}^{\frac{E}{2T}}{\left|\Gamma \left(\frac{r}{4}+i\frac{E}{2\pi T}\right)\right|}^2,$$

$$W\left(\left\{n_{dot}\right\}\right)=W_{eq}\left(\left\{n_{dot}\right\}\right)\left(1+\frac{eV}{T}\Psi \left(\left\{n_{dot}\right\}\right)\right),$$

$$\Psi \left(n_1,\dots ,n_{\alpha -1},1,n_{\alpha +1},\dots \right)=\Psi \left(n_1,\dots ,n_{\alpha -1},0,n_{\alpha +1},\dots \right)+\frac{g_r}{g_r+g_l}.$$

$$G=\frac{e^2}{h}\frac{g_l{g}_r}{g_r+g_l}\frac{\Lambda \delta }{T}\sum _N\sum _{\alpha }{W}_{eq}\left(N\right)\left[1-F_{\mathrm{eq}}\left({ϵ}_{\alpha }|N\right)\right]\left[\int dEP\left(E\right)f\left(E+{ϵ}_{\alpha }+U_{N+1}-U_N\right)\right],$$

$$W_{\mathrm{eq}}\left(N\right)=\frac{1}{1+e^{-\left(U_{N+1}-U_N\right)/T}}.$$

$$G=\frac{\Lambda e^2}{h}\frac{g_l{g}_r}{g_r+g_l}\frac{1}{1+e^{-\left(U_{N+1}-U_N\right)/T}}\left[\int dEP\left(E\right)\frac{\left(E+U_{N+1}-U_N\right)/T}{e^{\left(E+U_{N+1}-U_N\right)/T}-1}\right].$$

$$G^{\text{peak}}=\frac{\Lambda e^2}{h}\frac{g_l{g}_r}{g_r+g_l}\frac{{\left(2\pi \right)}^{\frac{r}{2}-2}}{\Gamma (\frac{r}{2})}{\left(\frac{T}{{\omega }_R}\right)}^{\frac{r}{2}}\left[\int ds{\left|\Gamma \left(\frac{r}{4}+i\frac{s}{\pi }\right)\right|}^2\frac{s}{\sinh s}\right]\propto T^{\frac{r}{2}}.$$

$$W_{\mathrm{eq}}\left(N\right)=\frac{1}{1+{\sum }_{\alpha }{e}^{-\left(U_{N+1}+{ϵ}_{\alpha }-U_N\right)/T}},$$

$$G=\frac{\Lambda e^2}{h}\frac{\delta }{T}\frac{g_l{g}_r}{g_r+g_l}\frac{1}{1+{\sum }_{\alpha }{e}^{-\left({ϵ}_{\alpha }+U_{N+1}-U_N\right)/T}}\sum _{\alpha }\int dEP\left(E\right)\frac{1}{1+e^{\left(E+{ϵ}_{\alpha }+U_{N+1}-U_N\right)/T}}.$$

$$W_{\mathrm{eq}}(N)=\frac{{\sum }_{\left\{n_{dot}\right\}}{\delta }_{even,\sum n_{\alpha }}\exp \left[-\frac{1}{T}\left({\sum }_{\alpha }{ϵ}_{\alpha }{n}_{\alpha }+U_N\right)\right]}{{\sum }_{\left\{n_{dot}\right\}}{\delta }_{even,\sum n_{\alpha }}\exp \left[-\frac{1}{T}\left({\sum }_{\alpha }{ϵ}_{\alpha }{n}_{\alpha }+U_N\right)\right]+{\sum }_{\left\{n_{dot}\right\}}{\delta }_{odd,\sum n_{\alpha }}\exp \left[-\frac{1}{T}\left({\sum }_{\alpha }{ϵ}_{\alpha }{n}_{\alpha }+U_{N+1}\right)\right]},$$

$$1-F_{\mathrm{eq}}\left({ϵ}_{\alpha }|N\right)=\frac{{\sum }_{\left\{n_{i\ne \alpha }\right\}}{\delta }_{even,\sum n_{\alpha }}\exp \left[-\frac{1}{T}\left({\sum }_{\alpha }{ϵ}_{\alpha }{n}_{\alpha }\right)\right]}{{\sum }_{\left\{n_{dot}\right\}}{\delta }_{even,\sum n_{\alpha }}\exp \left[-\frac{1}{T}\left({\sum }_{\alpha }{ϵ}_{\alpha }{n}_{\alpha }\right)\right]}.$$

$$G=\frac{\Lambda e^2}{h}\frac{\delta }{T}\frac{g_l{g}_r}{g_r+g_l}\frac{1}{F_{even}+F_{odd}{e}^{-\left(U_{N+1}-U_N\right)/T}}\sum _{\alpha }{D}_{even}\left({ϵ}_{\alpha }\right)\int dEP\left(E\right)\frac{1}{1+e^{\left(E+{ϵ}_{\alpha }+U_{N+1}-U_N\right)/T}},$$

$$F_{even/odd}=$$

$$D_{even/odd}\left({ϵ}_{\alpha }\right)=$$

$$W_{\mathrm{eq}}(N)=\frac{{\sum }_{\left\{n_{dot}\right\}}{\delta }_{odd,\sum n_{\alpha }}\exp \left[-\frac{1}{T}\left({\sum }_{\alpha }{ϵ}_{\alpha }{n}_{\alpha }+U_N\right)\right]}{{\sum }_{\left\{n_{dot}\right\}}{\delta }_{odd,\sum n_{\alpha }}\exp \left[-\frac{1}{T}\left({\sum }_{\alpha }{ϵ}_{\alpha }{n}_{\alpha }+U_N\right)\right]+{\sum }_{\left\{n_{dot}\right\}}{\delta }_{even,\sum n_{\alpha }}\exp \left[-\frac{1}{T}\left({\sum }_{\alpha }{ϵ}_{\alpha }{n}_{\alpha }+U_{N+1}\right)\right]},$$

$$1-F_{\mathrm{eq}}\left({ϵ}_{\alpha }|N\right)=\frac{{\sum }_{\left\{n_{i\ne \alpha }\right\}}{\delta }_{odd,\sum n_{\alpha }}\exp \left[-\frac{1}{T}\left({\sum }_{\alpha }{ϵ}_{\alpha }{n}_{\alpha }\right)\right]}{{\sum }_{\left\{n_{dot}\right\}}{\delta }_{odd,\sum n_{\alpha }}\exp \left[-\frac{1}{T}\left({\sum }_{\alpha }{ϵ}_{\alpha }{n}_{\alpha }\right)\right]}.$$

$$G=\frac{\Lambda e^2}{h}\frac{\delta }{T}\frac{g_l{g}_r}{g_r+g_l}\frac{1}{F_{odd}+F_{even}{e}^{-\left(U_{N+1}-U_N\right)/T}}\sum _{\alpha }{D}_{odd}\left({ϵ}_{\alpha }\right)\int dEP\left(E\right)\frac{1}{1+e^{\left(E+{ϵ}_{\alpha }+U_{N+1}-U_N\right)/T}}.$$

$$I_{2e}=-2e\left(W_0{\Gamma }_{0\to 2}^l-W_2{\Gamma }_{2\to 0}^l\right),$$

$$\frac{\partial W_0}{\partial t}=$$

Full text

Introduction–. The realization of Majorana zero modes (MZMs) Read and Green (2000); Kitaev (2001) provides a promising platform to study novel fundamental physics, e.g. non-Abelian braiding statistics M. Leinaas and Myrheim (1977); Fredenhagen et al. (1989); Ivanov (2001); Nayak et al. (2008), and has potential applications in quantum information processing and topological quantum computation Kitaev (2002); Nayak et al. (2008). Many proposals for realizing MZMs in topological superconductors (SCs) have been put forward  Fu and Kane (2008); Sato et al. (2009); Sau et al. (2010); Lutchyn et al. (2010); Oreg et al. (2010); Sau et al. (2010); Alicea (2010); Cook and Franz (2011); Alicea (2012), and lead to recent experimental progress in the realization and detection of MZM in both one-dimensional Mourik et al. (2012); Rokhinson et al. (2012); Deng et al. (2012); Das et al. (2012); Churchill et al. (2013); Finck et al. (2013); Nadj-Perge et al. (2014); Albrecht et al. (2016); Deng et al. (2016); Zhang et al. (2017); Gül et al. (2018); Zhang et al. (2018) and two-dimensional platforms Sun et al. (2016); Wang et al. (2018); Machida et al. (2019); Liu et al. (2018a); Fornieri et al. (2019); Ren et al. (2019); Yang et al. (2019). The quantized Majorana conductance at $2e^{2}/h$, observed in nanowire devices Zhang et al. (2018), closes one chapter in tunneling spectroscopy based on the simplest device set-ups. Additionally, some clues of Majorana conductance plateau were also shown in the vortex core of topological SC Zhu et al. (2019); Chen et al. (2019). Those experimental observations are gearing up for next-step discoveries Zhang et al. (2019) towards the realization of non-Abelian braiding experiments and Majorana qubits.

So far, experimental activities of Majorana research mainly focus on the Majorana resonance behavior or its oscillating splittings Albrecht et al. (2016). These experiments, aiming at verifying the existence of Majorana bound states, reveal little information on the non-local coherence nature between two Majorana pairs. It is this non-local coherent property which enables the non-local information storage and makes Majorana qubit robust against decoherence. Therefore experimentally demonstrating the Majorana non-local coherence would be a milestone towards the final braiding experiment and topological qubits. Among the earlier proposed schemes, a well-studied one is the fractional Josephson effect Kitaev (2001); H.-J. Kwon et al. (2004), which has undergone a lot of experimental activities Rokhinson et al. (2012); Laroche et al. (2019). However, the contamination mechanisms of the $4\pi$ periodicity mix both incoherent quasi-particle events Fu and Kane (2009); Lutchyn et al. (2010) and dynamical Landau-Zener processes Houzet et al. (2013); Badiane et al. (2013). Therefore, the deviation of $4\pi$ periodicity may not be a valuable benchmark for Majorana device coherence. Another important scheme is the non-local coherent teleportation proposed by Fu Fu (2010) for a floating Majorana superconducting island with finite charging energy. The electron transport through the island is mediated by the teleportation process: an electron injected into the MZM at one end can tunnel out from the other MZM at the other end while still maintaining its phase coherence even for a long wire-distance. In order to benchmark the coherent nature, an Aharonov-Bohm (AB) interference experiment Fu (2010) is thought to be necessary, where the electron coherence in a parallel normal channel needs also to be maintained in a loop-type device. Those requirements greatly increase the complexity and difficulty in experiments Gazibegovic et al. (2017). Therefore, a realistic scheme with more experimentally feasible conditions and sharp transport signals is preferable to reveal the Majorana non-local coherence.

In this letter, we propose such an experimentally feasible scheme based on a Majorana island device coupled to a dissipative bath. The dissipative bath can be realized by making part of the leads highly resistive (comparable to von-Klitzing resistance $h/e^{2}$) Ingold and Nazarov (1992); Mebrahtu et al. (2012, 2013). In our scheme, we only require a simple conductance measurement to reveal the coherence nature without the AB loop structure. The intuitive trust of this scheme is from the fact that the dissipative environments suppress all low-energy transport signals except for the symmetric coherent resonant level system Mebrahtu et al. (2012); Liu et al. (2014). We map our system to a dissipative resonant level model Mebrahtu et al. (2012); Liu et al. (2014), and show that, for symmetric lead-island couplings, the conductance from coherent Majorana teleportation increases as the temperature is reduced. While, for all other non-Majorana regimes where incoherent transports dominate, the conductance with the same parameter conditions show a power-law decay as the temperature is reduced. In addition, we want to point out that the dissipative electrode lead (which could help us to reveal the coherence or decoherence) in our scheme is quite different from the disorder Brouwer et al. (2011); Sau et al. (2012); Liu et al. (2012); Takei et al. (2013); Hui et al. (2015); Cole et al. (2016); Liu et al. (2018b); Chang et al. (2015); Gül et al. (2017) in the Majorana island, which can cause in-gap soft modes and contaminate their coherence. Therefore, the conductance enhancement behavior in the dissipative environment can serve as a clear benchmark for the coherent nature of Majorana devices.

Model– The system set-up is shown in Fig. 1(a). A semiconductor nanowire with Rashba spin-orbit coupling is in proximity with an s-wave superconductor. A magnetic field B is applied in parallel with the nanowire to realize a tunable topological SC phase Lutchyn et al. (2010); Oreg et al. (2010). The central gates control the nanowire chemical potential, while the tunnel-gates control the coupling between the (left and right) leads and the proximitized nanowire. The system also couples to a dissipative bath realized by adding on-chip resistors in source and drain leads Mebrahtu et al. (2012). This can be achieved by replacing part of the electrode (typically Au, 100 nm thick: brown in Fig. 1(a)) with a thin (around 10 nm thick: red in Fig. 1(a)) long resistive metal strip (e.g. Cr). For the dissipation effect to be effective, the resistive part needs to be on the device chip (at low temperature) and close the electrodes which form good Ohmic contacts with the nanowire. An equivalent circuit diagram is shown in Fig. 1(d) including a familiar quantum (lead-island) part coupled to a classical R-C circuit.

The central nanowire island is electrically floated such that electrons in the island feel an electrostatic energy: $U_{N}=E_{c}\left(N-V_{g}C_{g}/e\right)^{2}$, where $E_{c}=e^{2}/2C_{\Sigma}$ is the island charging energy and $C_{\Sigma}=C_{l}+C_{r}+C_{g}$ is the total capacitance. $N$ indicates the total excess electron number of the island. Here, $C_{g}$ ($C_{l}$, $C_{r}$) is the central gate (left, right junction) capacitance, and $V_{g}$ is the voltage applied on the central gate. The nanowire shows induced superconductivity with an energy gap $\Delta_{\textit{B}}$ at the magnetic field $B$ as shown in Fig. 1(e). The gap closes at critical $B_{c}$ indicating a topological phase transition from a trivial phase to a topological phase with MZMs at the island ends. The Hamiltonian of the island is

[TABLE]

Here, $\gamma_{\alpha}$ is the Bogoliubov quasiparticle operator and $\epsilon_{\alpha}$ is the quasiparticle energy. In trivial phases ($B<B_{c}$), the energy levels above the gap can be approximately expressed as $\epsilon_{\alpha}=\sqrt{\Delta_{\textit{B}}^{2}+\xi_{\alpha}^{2}}$ with $\xi_{\alpha}$ being the normal state energy. In topological SC phase ($B>B_{c}$), there are two-fold degenerate ground states separated from the continuum by a topological gap $\Delta_{\textit{Topo}}$.

Including the left and right leads along with dissipative environment (resistor $R_{l}$ and $R_{r}$ in two leads), the total Hamiltonian of the system can be written as

[TABLE]

where $H_{\textit{leads}}=\sum_{j,p\sigma}\xi_{j,p}c_{j,p\sigma}^{\dagger}c_{j,p\sigma}$ describes the free electrons in the left ($j=l$) and right ($j=r$) leads with electron energy $\xi_{p}$ at orbital $p$.

[TABLE]

which describes the tunneling between the nanowire and the leads. $d_{\sigma}\left(\mathbb{r}_{j}\right)$ and $c_{\sigma}\left(\mathbb{r}_{j}\right)$ are the electron field operator in the nanowire and lead separately, with spin $\sigma$ and evaluated at the specific position of the junction-$j$ sup . $t_{j}$ is the average tunneling amplitude at that junction. The tunneling across that single junction gives rise to the conductance $G_{j}=g_{j}(2e^{2}/h)$ with $g_{j}\propto\left|t_{j}\right|^{2}\nu_{j}\nu$ where $\nu_{j}$ ($\nu$) is the density of state (DOS) at the leads (nanowire). Later, we will require symmetric island-lead couplings $g_{l}=g_{r}\equiv g$ by fine tuning a single parameter. $H_{\textit{env}}$ describes the dissipative environment, and can be modeled by the classical circuit Ingold and Nazarov (1992). To incorporate the dissipation effect in our system, we add a phase operator $e^{-i\varphi_{j}}$ in the tunneling Hamiltonian Ingold and Nazarov (1992); Liu et al. (2014); and this operator couples the quantum tunneling with the classical circuit. This phase operator change the total charge of the junction capacitor by one electron (note the charge-phase conjugation $\left[\varphi_{j},Q_{j^{\prime}}\right]=ie\delta_{jj^{\prime}}$ where $Q_{j}$ is the charge of the junction capacitor). For convenience of calculations, we consider small gate capacitance $C_{g}\ll C_{l,r}$, then the only relevant phase is $\varphi=\varphi_{l}-\varphi_{r}$ which corresponds to the charge transfer between the left and right junction; and for general cases, the physics will not change qualitatively Liu et al. (2014). The correlation function of this phase operator can be obtained from the corresponding R-C circuit, and results in $\left\langle e^{i\varphi(t)}e^{-i\varphi(0)}\right\rangle\sim(\omega_{R}t)^{-2r}$ where $r=R/R_{K}$ with von-Klitzing resistance $R_{K}=h/e^{2}$, and $\omega_{R}=(RC)^{-1}$ with $R=R_{l}+R_{r}$ and $C=C_{l}C_{r}/(C_{l}+C_{r})$.

Coherent teleportation vs incoherent transport–. The electron transports through a confined island with charging energy show Coulomb blockade (CB) oscillations Aleiner et al. (2002); Bruus and Flensberg (2004). For a trivial SC island, the low-energy events correspond to Cooper-pair transport with $2e-$ periodicity in CB oscillation Hekking et al. (1993); van Heck et al. (2016). However, the presence of MZMs in a topological SC island causes a regular single electron tunneling with $1e-$ periodicity, which is the same as the periodicity for a normal metal island. Why is the Majorana teleportation special? In fact, the key underlying physics is quite unique. For Majorana teleportation Fu (2010) as shown in Fig. 1 (a), the Majorana wave-function is localized at the wire’s two ends, and overlaps strongly with electrode leads. The lead-Majorana coupling strength is represented by $\Gamma^{\textit{M}}\sim g\Delta_{\textit{Topo}}$ which is the reverse electron lifetime on the island electron. Then, for $T<\Gamma^{\textit{M}}$ and gate-voltage tuned close to CB peaks, the electron can tunnel into one Majorana state and out of the other one and maintain its phase coherence. On the other hand, the normal state wavefunction widely spreads out over the wire as shown in Fig. 1(c). In order to reach the other side, the electrons have to propagate through the whole wire with a typical time $\hbar/\Gamma^{D}$ (note $\Gamma^{D}\sim g_{j}\delta$ with a tiny level spacing $\delta$). Therefore, in the realistic regime, $g\delta\ll T\lesssim g\Delta_{\textit{Topo}}$ ( $\delta\ll\Delta_{\textit{Topo}}$ which requires a long island), only Majorana teleportation maintains phase coherence; and on the contrary, all other normal states generate incoherent transport near CB peak. Although the regular conductance measurement cannot tell the difference: both show 1e-periodic CBs, we will see the dissipative environment could help us to reveal the coherent signature without a more sophisticated interferometer set-up.

Dissipative Majorana teleportation–. As the nanowire enters into the topological region (@slowromancapi@) shown in Fig. 1(e), the transport is dominated by the resonant tunneling through the zero energy Majorana state. With dissipation, this system can be mapped to a dissipative resonant level model Mebrahtu et al. (2012); Liu et al. (2014). For symmetric island-lead couplings $g_{l}=g_{r}$, the dissipation effect cannot qualitatively change the low-energy transport behavior for an originally coherent resonant tunneling. For our model, as long as $r=R/R_{K}<2$ for arbitrary $\Gamma^{\textit{M}}$ or $2\leq r<3$ for $\Gamma^{\textit{M}}$ above a critical value, the CB peak height still increases as the temperature is reduced Mebrahtu et al. (2012); Liu et al. (2014) as shown in Fig. 2 (a). Dissipation only slightly changes the scaling function: At low temperature $T\ll\Gamma^{\textit{M}},\omega_{R}$, the peak height $G_{\textit{M}}^{\textit{peak}}$ reaches $e^{2}/h$ as $e^{2}/h-G_{\textit{M}}^{\textit{peak}}\propto T^{{2}/{(1+r)}}$ Zhang et al. (2016); while at higher $T$, the peak height increases as $T^{-2r}$ when $T$ decreases Liu et al. (2014). On the other hand, we will see dissipation significantly renormalizes the low-energy transport, and results in an opposite temperature dependence for other processes, which are incoherent, in the next few sections.

Dissipative sequential tunneling–. As mentioned earlier, in the non-topological regime, the electron transport near CB peaks through the nanowire loses the phase information. Those tunneling processes can be treated sequentially through the two junctions; and therefore, the conductance can be obtained from a rate equation method Beenakker (1991); Glazman and Pustilnik (2005). Including the dissipation effect, we derive the general form of the zero-bias conductance near the CB peaks (see SI sup for more details):

[TABLE]

which is from the processes that a lead-electron with energy $\xi_{p}$ jumps into an unoccupied level $\alpha$ with energy $\epsilon$α in the island. Note that $\Lambda\sim 1$ is a dimensionless parameter associated to the SC wavefunction sup . The energy conservation requires $\epsilon_{\alpha}+U_{N+1}+E=U_{N}+\xi_{p}$ where the energy exchange $E$ with environment is allowed. The probability function $P\left(E\right)$ of such an exchange can be written as the Fourier transform of the phase-phase correlation function Ingold and Nazarov (1992); sup

[TABLE]

where $\mathit{\Gamma}\left(z\right)$ is the gamma function. $W_{\mathrm{eq}}\left(N\right)$ is the probability that the isolated island contains $N$ electrons in equilibrium, $F_{\mathrm{eq}}\left(\epsilon_{\alpha}|N\right)$ is the equilibrium conditional probability that level $\alpha$ is occupied given the island electron number $N$, and $f(x)$ is the standard Fermi-Dirac distribution sup . Next, we will show how the dissipation affects the conductance in different non-topological regimes ($B<B_{c}$), respectively.

Near the SC gap closing point in region (@slowromancapii@) of Fig.1(e), the nanowire corresponds to a metallic phase and behaves like a normal metal island. For a long nanowire, the level spacing is small such that $\delta\ll T$, the discrete energy spectrum of the island can be treated as a continuum. In this limit, $W_{\mathrm{eq}}\left(N\right)\approx 1/\left[1+e^{-\left(U_{N+1}-U_{N}\right)/T}\right]$ and $F_{\mathrm{eq}}\left(\epsilon_{\alpha}|N\right)$ can be approximated as Fermi-Dirac distribution $f\left(\epsilon_{\alpha}\right)$ sup . Applying the integration in Eq. (4), we obtain the conductance near CB peaks with peak height proportional to $T^{r/2}$ sup . The CB conductance oscillations with $1e$-periodicity for different temperatures is shown in Fig. 2(b), where the dissipation parameter is $r=0.5$. Different from the coherent Majorana teleportation, the dissipation effect significantly suppresses the conductance at low energy. Note that this calculation is valid only near the CB peaks; and beyond the peak width (labeled by horizontal bars in the figures), the co-tunneling events need to be included in the calculation.

Secondly, we consider the case that the system has a finite but smaller gap $\Delta_{\textit{B}}<E_{c}$ (region (@slowromancapiii@)). In this case, adding an unpaired electron is still energetically preferred, so the single electron tunneling contributes to the conductance Averin and Nazarov (1992); van Heck et al. (2016). The degeneracy happens at $U_{N}\approx U_{N+1}\pm\Delta_{\textit{B}}$ ("$+$" for even $N$ and "$-$" for odd $N$); therefore, the position of the peaks will shift to the right (left) by about $N_{\Delta}$ for even (odd) $N$ as shown in Fig. 2(c). In the conductance calculation, we need to include the tunneling processes where electrons enter into the levels above the SC gap in the island, and thus we have $W_{\mathrm{eq}}\left(N\right)\approx 1/\left[1+\sum_{\alpha}e^{-\left(\epsilon_{\alpha}+U_{N+1}-U_{N}\right)/T}\right]$ near the charge degenerate point. In region (@slowromancapiii@), the SC gap is still big enough ($\Delta_{\textit{B}}\gg T$) to prohibit thermal excitations of quasi-particle states, so $\left[1-F_{\mathrm{eq}}\left(\epsilon_{\alpha}|N\right)\right]\approx 1$. Then following Eq.(4) we can obtain the conductance near the CB peaks sup . The conductance peak height shows a special temperature scaling $T^{r/2-\eta_{B}}$, where an extra scaling factor $\eta_{B}$ tends to one for $\Delta_{\textit{B}}\gg T$ and zero at the gap closing point (refer to SI sup for more details). The dissipative conductance at $r=1.9$ and $\Delta_{B}=0.86K$ is shown in Fig. 2(c). We can see that the dissipation still renormalizes the low energy conductance in quite a different way compared with Majorana teleportation.

Finally, we consider the region (@slowromancapiv@) corresponding to a trivial SC phase with $\Delta_{\textit{B}}>E_{c}$. In this region, the island at any $V_{g}$ is a pool of Cooper pairs at BCS ground state. Then the conductance peak arises only at the degenerate point of adding two electrons (e.g. electron number from $N$ to $N+2$) to the nanowire, because single electron tunneling process is blockaded by the SC gap. In this case, the conductance near the CB peaks due to the sequential tunneling of Cooper pairs is (see SI sup for more details)

[TABLE]

which describes the Andreev reflection process that a lead-electron in the level $p_{1}$ jumps into the nanowire, and reflects a hole to the level $p_{2}$ in the lead. The energy conservation requires $\xi_{p_{1}}+\xi_{p_{2}}+U(N)=U(N+2)+E$ where $E$ is the energy exchange with the environment. Near the peak, $W_{\mathrm{eq}}\left(N\right)\approx\left(1+e^{-\left(U_{N+2}-U_{N}\right)/T}\right)^{-1}$. Because the Andreev reflection is a second-order process, the lead-island coupling becomes: $\Gamma^{2e}=Ag^{2}T/2\pi$ ($A\sim 1$ is a dimensionless number representing the amplitude for the Andreev reflection sup ). The transport is incoherent under the condition $\Gamma^{2e}\ll T$ which is is equivalent to saying $g^{2}\ll 1$. The dissipative CB conductance of Cooper pairs at $r=0.5$ is shown in Fig. 3(c). The peak height is also suppressed as $T$ decreases to zero.

Conclusion and discussion–. So far, we have shown that under the condition $g\delta\ll T\lesssim g\Delta_{\textit{Topo}}$, $g^{2}\ll 1$, all the states in trivial phases carry incoherent transport, which will be suppressed as reducing the temperature in the presence of a dissipative environment. While the Majorana teleportation will maintain coherence, and the conductance will increase when the temperature goes down. This behavior could be a clear benchmark to determine the coherent nature of Majorana devices.

In addition, we emphasize that this two-terminal transport experiment for certain conditions can also tell the MZM from other contaminations like Andreev bound states (ABS), especially the zero-energy ABS (or equivalently two decoupled quasi-Majorana state) which can arise from smooth potential at wire end Kells et al. (2012); Liu et al. (2017); Moore et al. (2018a); Vuik et al. (2018); Moore et al. (2018b). The single-terminal tunneling conductance of such devices can also have a quantized zero-bias peak Moore et al. (2018a) which could mimic the MZM. By adjusting the potential of one end sup , the two quasi-Majorana components can have finite coupling regardless of the situation (either MZM or ABS) of the other end, which will generate the splitting in the single-terminal tunneling conductance as shown in Fig. 3 (d), similar to the splitting of Majorana case (Fig. 3 (b)). Denoting the splitting of MZM and ABS as $\epsilon_{0}$ and $\epsilon_{A}$, when $T\ll\epsilon_{0/A}$, the splittings in these two cases will lead to completely different two-terminal CB conductance signals. In the Majorana case, the splitting comes from the coupling between the two Majorana modes at two ends, which will only shift the position of the CB peak by $\epsilon_{0}$ as shown in Fig. 3 (a). In the ABS case, the splitting comes from the coupling between two quasi-Majorana components of ABS located at the same end; therefore, the two quasi-Majorana modes form a localized fermionic state sup . Then the non-local teleportation disappears, while only the CB peaks caused by the Cooper-pair tunneling (2-e) exist as shown in Fig. 3 (c). In this case, the single electron tunneling needs the electron to go through the wire via the localized ABS. Because the coupling of the ABS to the near lead is big and to the far lead is exponentially small, this is a highly asymmetric resonant tunneling process. The resulting conductance is quite small and according to the results in Mebrahtu et al. (2012); Liu et al. (2014), dissipation can significantly suppress the conductance peak, i.e. showing opposite temperature dependence compared to the true Majorana case.

Supplementary Material for

“Reveal the non-local coherent nature from a dissipative Majorana teleportation”

In this supplementary material, we will provide some details about:  A.) Estimations of experimental parameters, B.) the dissipative single electron tunneling conductance near the CB peak, C.) the dissipative Cooper pair tunneling conductance near the CB peak, D.) Distinguish between Majorana and ABS from two-terminal transport.

Appendix A Estimations of experimental parameters.

The equivalent circuit diagram of the set-up is shown in Fig. 1(d) of the main text. Two junctions on both sides of the nanowire are characterized by capacitance $C_{j}$ and conductance $g_{j}$. The total capacitance between the nanowire and the ground is $C_{\Sigma}=C_{l}+C_{r}+C_{g}$ which determines the charging energy $E_{c}=e^{2}/2C_{\Sigma}$. While the capacitance of the RC circuit is $C=C_{l}C_{r}/\left(C_{l}+C_{r}\right)$. The total resistance of the RC circuit is $R=R_{l}+R_{r}$, and the dissipation strength can be described by a dimensionless parameter $r=R/R_{K}$ where $R_{K}=h/e^{2}$ is the von-Klitzing resistance. For the convenience of the discussion, the capacitances of two junctions are taken the same with $C_{l}=C_{r}$ (e.g. choose $2\times 10^{-16}\text{F}$ for the energy scale comparison) and $C_{g}$ can be small $C_{g}\ll C_{l,r}$; for the general capacitance cases, our main conclusion is still valid Liu et al. (2014). Then the charging energy are of the value $E_{c}\approx e^{2}/\left[\left(C_{l}+C_{r}\right)\right]=2.32K$ and $\omega_{R}=\left(RC\right)^{-1}=3r^{-1}K$. With $r$ ranging from $0.5$ to $1.9$ and temperature $T$ ranging from $0.02K$ to $0.12K$ in the main text. The induced SC gap $\Delta_{B}$ is chosen such that: $\Delta_{B}>E_{c}\gg T$ in region (@slowromancapiv@), $E_{c}>\Delta_{B}\gg T$ in region (@slowromancapiii@), and gap-closing ($\Delta_{B}$ becomes the smallest energy scale) in region (@slowromancapii@). In principle, all the parameters are experimentally achievable and meet the conditions of the energy scales we study.

Appendix B The dissipative single electron tunneling conductance near the CB peak.

The nanowire can be treated as a quantum dot shown in Fig. 4. The quantum dot and two leads are connected through the left and right barrier potentials. The energy levels of the quantum dot includes discrete levels with level spacing $\delta$ and an induced SC gap. The Fermi seas for the leads are shown in the red region, where the chemical potential of the left (right) lead is zero ($eV$ by the bias voltage). When the induced SC gap is small ($\Delta_{B}<E_{c}$), the tunneling current also includes the contribution from the single electron sequential tunneling processes. In this case, the current through the left barrier (which is the same as that through the right) is given by

[TABLE]

where $\Gamma_{0\rightarrow\alpha}^{j}$ is the tunneling rate for the cases where an electron tunnels from the lead$-j$ to the state $\left|\alpha\right\rangle$ of the dot, and $\Gamma_{\alpha\rightarrow 0}^{j}$ is the tunneling rate of the reverse process. $n_{\alpha}$ (0 or 1) is the occupation number of the dot state $\alpha$ , and $W\left(\left\{n_{dot}\right\}\right)$ is the probability of the occupation $\left\{n_{dot}\right\}=\left\{n_{1},n_{2},\ldots\right\}$. The electron tunneling from leads to a state $\left|\alpha\right\rangle$ in the dot and from $\left|\alpha\right\rangle$ to the leads change the probability $W\left(\left\{n_{dot}\right\}\right)$. The changing rate should satisfy the master equation:

[TABLE]

Subsequently, we look at how to compute tunneling rate $\Gamma_{\alpha 0}^{j}$ for our model, solve the probability $W\left(\left\{n_{dot}\right\}\right)$, and obtain the conductance formula.

The total Hamiltonian of the system is shown in Eq. (2) of the main text. The states of the leads are labeled by $p$ and the states of the quantum dot are labeled by $\alpha$. The tunneling Hamiltonian is $H_{\textit{T}}=\sum_{j,\sigma}t_{j}d_{\sigma}^{\dagger}\left(\mathbb{r}_{j}\right)c_{\sigma}\left(\mathbb{r}_{j}\right)e^{-i\varphi_{j}}+\textit{H.c.}$, where $d_{\sigma}\left(\mathbb{r}_{j}\right)=L^{\frac{1}{2}}\sum_{\alpha}u_{\alpha}\left(\mathbb{r}_{j},\sigma\right)\gamma_{\alpha}+v_{\alpha}^{*}\left(\mathbb{r}_{j},\sigma\right)\gamma_{\alpha}^{\dagger}$ and $c_{\sigma}\left(\mathbb{r}_{j}\right)=\Omega_{j}^{\frac{1}{2}}\sum_{p,\sigma}\phi_{p}\left(\mathbb{r}_{j}\right)c_{p\sigma}$ are the electron field operators in the nanowire and the lead respectively, with $u_{\alpha}/v_{\alpha}\left(\mathbb{r}_{j},\sigma\right)$ and $\phi_{p}\left(\mathbb{r}\right)$ being the wavefunctions of BDG Hamiltonian and the wavefunction of the lead-electrons respectively, where $\sigma=\uparrow,\downarrow$ is the spin index and $j=l,r$ is the junction index. $L$ is the length of the nanowire and $\Omega_{j}$ is the volume of lead$-j$. Treating $H_{\textit{T}}$ as perturbation we can obtain the tunneling rate

[TABLE]

where $f\left(x\right)=\left(1+e^{x/T}\right)^{-1}$ is the Fermi-Dirac distribution. $\mu_{j}$ is the chemical potential of the lead-$j$ and we have set $\mu_{l}=0$, $\mu_{r}=eV$ as show in Fig. 4. The barrier transmission is proportional to the dimensionless conductance $g_{j}=4\pi^{2}\left|t_{j}\right|^{2}\nu_{j}\nu$, where $\nu_{j}$ is the density of states (DOS) of the lead-$j$ and $\nu=\delta^{-1}$ is the DOS of the nanowire when the SC gap closes ($\Delta_{B}=0$). $\Lambda$ is related to the SC-nanowire wavefunction at the tunneling location: $\Lambda=L\sum_{\sigma}\left|u_{\alpha}\left(\mathbb{r}_{j},\sigma\right)\right|^{2}$. Because the temperature is much smaller than the gap of the s-wave SC shell, almost the entire wavefunction $u_{\alpha}\left(\mathbb{r},\sigma\right)$ that contribute to the transport is located inside the nanowire and thus $\Lambda\sim 1$ which is treated as a constant value in region (@slowromancapii@) and (@slowromancapiii@) of Fig. 1(e) in the main text van Heck et al. (2016).

The function $P\left(E\right)=\left(2\pi\hbar\right)^{-1}\int_{-\infty}^{+\infty}\mathrm{d}t\left\langle e^{i\varphi_{j}(t)}e^{-i\varphi_{j}(0)}\right\rangle e^{iEt/\hbar}$ which can be interpreted as the probability of exchanging energy $E$ with the environment. They follow the standard relation $P\left(-E\right)=P\left(E\right)e^{-E/T}$, which means the ratio between the probability to emit energy into the environment and the probability to absorb energy from the environment is a Boltzmann factor. In the absence of dissipation, we have $P\left(E\right)=\delta\left(E\right)$ and Eq. (10) becomes the familiar tunneling rate Glazman and Pustilnik (2005); van Heck et al. (2016). With $C_{l}=C_{r}$, $\left\langle e^{i\varphi_{j}(t)}e^{-i\varphi_{j}(0)}\right\rangle=\left\langle e^{i\varphi(t)/2}e^{-i\varphi(0)/2}\right\rangle$ where the phase $\varphi=\varphi_{l}-\varphi_{r}$ corresponds to the charge transfer between the left and right junction. For ohmic dissipation and $\omega_{R}/T\gg 1$, we can evaluate the integral

[TABLE]

where $\mathit{\Gamma}\left(z\right)$ is the gamma function.

To solve the master equation for probability function, we assume an expansion form for $W\left(\left\{n_{dot}\right\}\right)$ in order of bias voltage $V$, and only keep the linear order terms

[TABLE]

where $W_{eq}\left(\left\{n_{dot}\right\}\right)$ is the probability distribution in equilibrium when no bias voltage is applied and has the form $W_{eq}\left(\left\{n_{dot}\right\}\right)=Z^{-1}\exp\left[-\left(\sum_{\alpha}\epsilon_{\alpha}n_{\alpha}+U_{N}\right)/T\right]$ where $Z=\sum_{\left\{n_{dot}\right\}}\exp\left[-\left(\sum_{\alpha}\epsilon_{\alpha}n_{\alpha}+U_{N}\right)/T\right]$. The steady state occupation probabilities can be determined by setting $\partial W\left(\left\{n_{dot}\right\}\right)/\partial t=0$. Then for each $\alpha$ in Eq. (9), we can obtain

[TABLE]

Combining Eq. (8), (13) and (14) and keeping only the leading terms in $V$, we obtain the zero-bias conductance:

[TABLE]

where $W_{eq}\left(N\right)=\sum_{\left\{n_{dot}\right\}}W_{eq}\left(\left\{n_{dot}\right\}\right)\delta_{N,\sum_{\alpha}n_{\alpha}}$ is the probability that the isolated dot contains $N$ electrons in equilibrium, and $F_{\mathrm{eq}}\left(\epsilon_{\alpha}|N\right)=W^{-1}_{\mathrm{eq}}(N)\sum_{\left\{n_{dot}\right\}}P_{\mathrm{eq}}\left(\left\{n_{dot}\right\}\right)\delta_{n_{\alpha},1}\delta_{N,\Sigma_{\alpha}n_{\alpha}}$ is the equilibrium conditional probability that level $\alpha$ is occupied given the dot electron number $N$. In the main text, we only consider the limit $T\ll E_{c}$ and the region near the CB peak. In this limit, we can neglect the summation over $N$ in Eq. (15).

B.1 The normal metal regime.

Near the SC gap closing point in region (@slowromancapii@) of Fig. 1(e) of the main text, the nanowire corresponds to a metallic phase and behaves like a normal metal island. In the case $\delta\ll T$, the discrete energy spectrum in the island can be treated as a continuum. In this limit, $F_{\mathrm{eq}}\left(\epsilon_{\alpha}|N\right)$ can be approximated as Fermi-Dirac distribution $f\left(\epsilon_{\alpha}-\mu\left(N\right)\right)$ and $\mu\left(N\right)$ can be determined by $\sum_{\alpha}f\left(\epsilon_{\alpha}-\mu\left(N\right)\right)=N$, and choose $\mu\left(N\right)\approx\mu=0$. Near the CB peak, a lead-electron can tunnel into a dot state near the Fermi surface; then the energy change of the dot is $U_{N+1}-U_{N}$. Therefore, including both $N$ and $N+1$ electron state, the probability of $N$ electrons is

[TABLE]

The terms beyond $U_{N}$ and $U_{N+1}$ are neglected for $T\ll E_{c}$. Performing the integration (the summation over $\alpha$) we can obtain the conductance in the metal phase:

[TABLE]

The results shown in Fig. 2(b) of the main text are obtained numerically from evaluating Eq. (17). The CB peaks appear at $U_{N+1}-U_{N}=0$, from which we can obtain the peak conductance:

[TABLE]

B.2 The single electron tunneling regime with $T\ll\Delta_{B}$.

In the region (@slowromancapiii@) with the limit $T\ll\Delta_{B}$, the SC gap is big enough to prohibit thermal excitations of quasiparticle states, and therefore the equilibrium conditional probability for level $\alpha$ simply becomes $F_{\mathrm{eq}}\left(\epsilon_{\alpha}|N\right)=0$. Near the CB peak, a lead-electron can tunnel into a dot quasiparticle state with energy $\epsilon_{\alpha}$; then the energy change of the dot is $U_{N+1}+\epsilon_{\alpha}-U_{N}$. Therefore, including both $N$ and $N+1$ electron state, the probability of $N$ electrons is

[TABLE]

Feeding Eq. (19) into Eq. (15), we can obtain the conductance:

[TABLE]

The results shown in Fig. 5(a) below and Fig. 3(c) of the main text are obtained by numerically evaluating Eq. (20). Because the states with energy $\epsilon_{\alpha}$ close to $\Delta_{B}$ contribute the most to the conductance, the peak positions will be shifted by about $\Delta_{B}/\left(2E_{c}\right)$. Besides, the conductance shows a special temperature scaling $T^{r/2-\eta_{B}}$, where an extra scaling factor $\eta_{B}$ tends to one for $\Delta_{\textit{B}}\gg T$ and zero at the gap closing point.

B.3 The regime $T\lesssim\Delta_{B}$.

To find the behavior of the scaling factor $\eta_{B}$ we consider the crossover between single electron (still SC) regime (@slowromancapiii@) and a metallic (non-SC) regime (@slowromancapii@). Between (@slowromancapii@) and (@slowromancapiii@), $T<\Delta_{B}$ but not in the limit $T\ll\Delta_{B}$, the quasiparticle states could be thermally excited which can modify the transport signatures. In this phase, the quasiparticle number is not conserved, only its parity is conserved, the total number of excited states will have the same parity with $N$. When $N$ is even, we have

[TABLE]

[TABLE]

Substituting these equations into Eq. (15), we obtain

[TABLE]

where

[TABLE]

When $N$ is odd,

[TABLE]

[TABLE]

Substituting the upper two formulas into Eq. (15)

[TABLE]

By numerically evaluating Eq. (23) and (28), we obtain the CB conductance shown in Fig. 5(b). Together with Eq. (20), we numerically plot the $\eta_{B}$ as a function of SC gap $\Delta_{B}$, and find that $\eta_{B}$ decreases as the gap is reduced as shown in Fig. 5(c). Note that we have $\eta_{B}=0$ at the gap closing point.

For this regime, we in principle need to consider contributions of all the levels above the SC gap $\Delta_{B}$. The level spacing $\delta$ above the gap is very small $\delta\ll T\lesssim\Delta_{B}$, which results in a large number of levels above the gap, and the exact numerical evaluation becomes intractable. In fact, we only want to qualitatively understand the behavior of the scaling factor $\eta_{B}$ in conductance scaling $G\propto T^{r/2-\eta_{B}}$ as shown in the main text, and look at how does $\eta_{B}$ change as we reduce the SC gap (or vary magnetic field). Therefore, we include a truncation and only consider a few energy levels above the gap, which needs a relatively big $\delta$. In the numerical evaluation of the scaling factor $\eta_{B}$ of the main text, we take the level spacing $\delta=0.065K$ and $T\sim 0.1K$. Since $\Delta\sim T$, it seems that the situation is quite different from the real case in which $T\gg\delta$. Here we emphasize that as we decrease $\delta$ toward more realistic situation, the scaling factor $\eta_{B}$ will decrease as shown in Fig. 5(d), which means that with the same $r$, the CB peak conductance will decrease faster as temperature is reduced, when $\delta$ becomes smaller (toward the value of a more realistic case). So, more realistic situations with smaller level spacing can even strengthen our results (peak conductance for non-Majorana incoherent regimes is strongly suppressed as the temperature is reduced in the presence of dissipation).

Appendix C The dissipative Cooper pair tunneling conductance near the CB peak.

In region (@slowromancapiv@), $\Delta_{B}>E_{c}$ and the CB conductance exhibits oscillations with $2e-$ periodicity. Near the CB peak, the transport is dominated by the Cooper pair tunneling or the Andreev reflection process, while the single electron tunneling is suppressed by the SC gap. In this case, the current through the left barrier is given by

[TABLE]

where $\Gamma_{0\rightarrow 2}^{j}$ is the tunneling rate of an Andreev reflection process, which changes the number of electrons in the dot by two by adding a Cooper pair from lead-$j$ and $\Gamma_{2\rightarrow 0}^{j}$ is the tunneling rate of the reverse process which removes a Cooper pair from the dot to lead-$j$. $W_{2}$ and $W_{0}$ are the probabilities that the dot has and doesn’t have an additional Cooper pair, respectively. The tunneling of Cooper pairs into and from the dot will change the probabilities, and give us the standard master equation

[TABLE]

From the second-order perturbation theory, we can obtain the Cooper pair tunneling rates:

[TABLE]

Note here $P\left(E\right)$ becomes

[TABLE]

In Eq. (33) and (32), $A$ is a dimensionless number representing the amplitude of the Andreev reflection:

[TABLE]

For the steady state, we can set $\partial W_{0}/\partial t=0$ and $\partial W_{2}/\partial t=0$ in Eq. (30) and (31). We also have the probability conservation condition $W_{0}+W_{2}=1$. After solving $W_{0}$ and $W_{2}$, we obtain the CB conductance due to the Cooper pair tunneling:

[TABLE]

which is the Eq. (6) in the main text. The CB peaks appear at $U_{N+2}-U_{N}=0$, from which we can obtain the peak conductance:

[TABLE]

Appendix D Distinguish between Majorana and ABS from two-terminal transport.

Here we show that how the two-terminal transport experiment can tell the MZM from the near-zero-energy Andreev bound states (ABS). The zero-energy ABS can arise from smooth potential at wire end Kells et al. (2012); Liu et al. (2017); Moore et al. (2018a); Vuik et al. (2018); Moore et al. (2018b). Fig. 6(f) shows the energy spectrum when there are smooth Gaussian-shape potentials at two wire ends (as shown in Fig. 6(d)). When the Zeeman field is big enough but smaller than the topological transition point $V_{Z}=4.52meV$, zero-energy ABSs show up at two wire ends. The single-terminal tunneling conductance taken at one end of such devices can also have a zero-bias peak which could mimic the MZM. The four decoupled quasi-Majorana components (whose definition can be found in Moore et al. (2018a)) which make up two ABSs are plotted in Fig. 6(e) (two quasi-Majorana components make up one ABS). By making the potential sharper at one end (see Fig. 6(g)), one can couple the two quasi-Majorana states at that end and split the zero-energy spectral line (the green line in Fig. 6(i)) of the corresponding ABS. In the single-terminal tunneling experiment, this will only split the conductance peak (see Fig. 3(d) of the main text), which could also mimic the Majorana splitting due to the finite size (Fig. 3(b) of the main text). However since the two coupled quasi-Majorana states (green solid and dashed lines in Fig. 6(h)) are coupled and make up a localized fermionic ABS at the left end, which can no longer pair up with the states at the other end, the non-local teleportation via two quasi-Majorana states localized at different ends disappears. For two-terminal transport, a single electron can only tunnel through the whole wire via this fermionic ABS, which only strongly couples to one lead and exponentially decouples to the other. For such a highly asymmetric resonant tunneling, the conductance should be quite small; and in addition, the conductance will also be suppressed when decreasing the temperature in the case with dissipation Mebrahtu et al. (2012); Liu et al. (2014). In the Majorana case as shown in Fig. 6(c), even if there is splitting, the two Majorana modes shown in Fig. 6(b) are still mostly localized at each end. The coherent Majorana teleportation still exists and the resonant peak height still increases approaching the quantized value when lowering the temperature in the case with dissipation. Only the resonant peak position will be shifted. In conclusion, by tuning the sharpness of the potential, the dissipative transport can distinguish the coherent Majorana teleportation with the incoherent ABS transport.