Next Article in Journal
Effects of Structure and Composition of Adsorbents on Competitive Adsorption of Gaseous Emissions: Experiment and Modeling
Next Article in Special Issue
Phase-Selective Epitaxy of Trigonal and Orthorhombic Bismuth Thin Films on Si (111)
Previous Article in Journal
Editorial for Special Issue “Multifunctional Nanomaterials and Hybrid Structures for Sensors, Actuators and Smart Technologies”
Previous Article in Special Issue
Effects of the Spatial Extension of the Edge Channels on the Interference Pattern of a Helical Josephson Junction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust and Fragile Majorana Bound States in Proximitized Topological Insulator Nanoribbons

1
Peter Grünberg Institute 9, Forschungszentrum Jülich & JARA Jülich-Aachen Research Alliance, 52425 Jülich, Germany
2
JARA-Institute for Green IT, RWTH Aachen University, 52056 Aachen, Germany
3
Blackett Laboratory, Imperial College London, South Kensington Campus, London SW7 2AZ, UK
*
Authors to whom correspondence should be addressed.
Nanomaterials 2023, 13(4), 723; https://doi.org/10.3390/nano13040723
Submission received: 30 December 2022 / Revised: 2 February 2023 / Accepted: 8 February 2023 / Published: 14 February 2023
(This article belongs to the Special Issue Topological Materials in Low Dimensions)

Abstract

:
Topological insulator (TI) nanoribbons with proximity-induced superconductivity are a promising platform for Majorana bound states (MBSs). In this work, we consider a detailed modeling approach for a TI nanoribbon in contact with a superconductor via its top surface, which induces a superconducting gap in its surface-state spectrum. The system displays a rich phase diagram with different numbers of end-localized MBSs as a function of chemical potential and magnetic flux piercing the cross section of the ribbon. These MBSs can be robust or fragile upon consideration of electrostatic disorder. We simulate a tunneling spectroscopy setup to probe the different topological phases of top-proximitized TI nanoribbons. Our simulation results indicate that a top-proximitized TI nanoribbon is ideally suited for realizing fully gapped topological superconductivity, in particular when the Fermi level is pinned near the Dirac point. In this regime, the setup yields a single pair of MBSs, well separated at opposite ends of the proximitized ribbon, which gives rise to a robust quantized zero-bias conductance peak.

1. Introduction

Three-dimensional (3D) topological insulators (TIs) have received a lot of attention in the last decade due to their interesting electronic properties, in particular due to their topologically protected surface states with a spin-momentum-locked Dirac-cone energy spectrum [1]. This interest has only increased when the possibility emerged to realize exotic forms of superconductivity by combining TIs and ordinary s-wave superconductors in heterostructures [2,3,4,5,6,7,8], exploiting the superconducting proximity effect [9]. Due to strong spin-orbit coupling in the TI, the induced superconductivity transforms into p-wave pairing for the TI surface states [10].
A promising application of p-wave superconductivity is to realize Majorana bound states (MBSs) in a spinless fermionic channel [11], forming a so-called Majorana wire. These MBSs come in pairs of states at zero energy, which are localized at opposite ends of the wire. A pair of MBSs can also be understood as a pair of quasiparticles forming equal-weight superpositions of a particle and hole state. As such, the MBS is a quasiparticle that is its own anti(quasi)particle, with an associated creation/annihilation operator that is self-adjoint [12]. Furthermore, these MBSs are anyons with non-Abelian exchange statistics [13]. By combining these exotic properties, MBSs can be exploited for quantum information processing with the promise of being immune against the most common sources of decoherence [14,15,16,17,18].
With the abovementioned properties, a superconductor-TI nanowire or nanoribbon heterostructure appears to be a natural Majorana wire candidate. Unfortunately, the surface-state spectrum suffers from a spin degeneracy, which naturally arises due to confinement quantization of the spin-momentum-locked Dirac cone spectrum with antiperiodic boundary conditions for the Dirac spinor solutions [19]. This prevents the realization of topological p-wave superconductivity with a single spinless channel that is underlying the realization of a Majorana wire. The unwanted degeneracy can be lifted, however, by applying an external magnetic field along the wire [19,20,21,22]. The magnetic flux through the cross section of the TI nanowire modifies the boundary condition for the surface states and thereby lifts the degeneracy that prevents the formation of a topologically nontrivial regime [23,24].
Even when the proximitized TI nanowire is brought into the topological regime with an external magnetic field, there is no guarantee that well-separated MBSs form at opposite ends of the wire. For this, a sizeable proximity-induced superconducting gap must be induced in the surface-state spectrum of the TI nanowire. In recent works, the realization of such a fully gapped nontrivial regime was brought into question, as it was found that the particle and hole states fail to couple due to a mismatch in transverse momentum and, hence, the proximity effect fails to induce a superconducting gap [25]. To overcome this mismatch, it has been proposed to consider a vortex in the superconducting condensate that envelopes the TI wire [25], or to break the transverse symmetry of the wire with an electric field by introducing an electrostatic gate in the device layout [26].
In this article, we scrutinize the conditions for fully gapped topological superconductivity in a proximitized TI nanoribbon (i.e., a nanowire with rectangular cross section) structure, considering a realistic device layout that is compatible with state-of-the-art nanofabrication processes [6]. In particular, we consider a selectively grown TI nanoribbon that is covered from the top by a conventional superconductor (see Figure 1a), e.g., Nb [27]. Through careful consideration of the proximity effect, we find that this setup naturally yields optimal conditions for fully gapped topological superconductivity when the ribbon cross section is pierced by (close to) half a magnetic flux quantum, with a single pair of robust MBSs appearing at the ends of the nanoribbon. For other flux values, we identify different gapless and gapped phases with different numbers of end-localized MBSs, depending on the position of the Fermi level with respect to the Dirac point. Some of these MBSs appear to be fragile when disorder is introduced in the system and suffer from hybridization with other MBS solutions on the same end of the ribbon. We also consider a tunneling spectroscopy setup for distinguishing robust and fragile MBSs and identifying these different phases in experiments.
Starting with this introduction, the article is divided into six sections. In Section 2, we cover the details of our simulation approach, including the continuum model Hamiltonian (Section 2.1), the treatment of the superconducting proximity effect (Section 2.2), and the tight-binding modeling approach (Section 2.3). In Section 3, we discuss the (proximity-induced) spectral gap and the topologically trivial and nontrivial regimes of a top-proximitized TI nanoribbon. In Section 4, we discuss the different phases of this system, with different numbers of robust and fragile MBSs, which can be probed with tunneling spectroscopy (Section 4.1). We proceed with a discussion of the simulation results in Section 5 and conclude in Section 6.

2. Model

2.1. Topological Insulator

For obtaining the TI nanoribbon spectrum, we consider a tight-binding model (see Section 2.3) that is derived from the following 4-band continuum model Hamiltonian for the Bi2Se3 family of TI materials [28,29]:
H 0 ( k ) = [ ϵ ( k ) μ ] + M ( k ) σ z + A ( k y s x k x s y ) σ x + A z k z σ y ϵ ( k ) C 0 C ( k x 2 + k y 2 ) C z k z 2 , M ( k ) M 0 M ( k x 2 + k y 2 ) M z k z 2 .
H 0 ( k ) is a 4 × 4 matrix that is written as a linear combination of a tensor product of Pauli matrices s a and σ b ( a , b { x , y , z } ) acting on the spin and orbital (pseudospin) subspaces, respectively. Note that the identity matrices are not written explicitly here. The parameters C 0 , C , C z , M 0 , M , M z , A , A z can be obtained for different TI materials, such as Bi2Se3 or Bi2Te3 [29]. Here, we will neglect in-plane (⊥) versus out-of-plane (z) anisotropy (in terms of model parameters, A = A z A , C = C z C , M = M z M ) and asymmetry between valence and conduction bands ( C = 0 ) for simplicity (these simplifications will not significantly affect the surface-state spectrum near the Dirac point [30], which is the focus of this work). We consider remaining model parameters A = 3 eV · Å , M = 15 eV · Å 2 , and M 0 eV = 0.3 to represent a 3D TI material with inverted (direct) bulk band gap at the Γ point equal to 0.6 eV and a Dirac velocity for the 3D TI surface states equal to v D = A / 4.6 × 10 5 m / s , which are both comparable to Bi2Se3 [28] (note that 0.6 eV reflects the band separation at the Γ point and not the overall bulk band gap, which is closer to 0.3 eV ).
For describing a charge carrier density fluctuations in the ribbon (i.e., electrostatic disorder), we can add a disorder term S dis ϕ ( r ) to the parameter C 0 ( μ = C 0 corresponds to the chemical potential being pinned at the Dirac point in a pristine system without disorder). This disorder is characterized by a fluctuation amplitude S dis and the function ϕ ( r ) , which we consider to be a unit-normalized white-noise profile [ ϕ ( r ) ϕ ( r ) = δ ( r r ) ] or a Gaussian random field [ ϕ ( r ) ϕ ( r ) = e ( r r ) 2 / ( 2 λ 2 ) ] with spatial correlation length λ . Here, we consider an electrostatic disorder strength of the order of the TI surface-state subband spacing (∼ 2 π A / P , see below, which is in the 10 meV range) and a spatial correlation length in the few-nm range [31].

2.2. Proximity-Induced Superconductivity

In our setup, we consider a 3D TI nanoribbon that is proximitized via its top surface by an s-wave superconductor. We make use of the Bogoliubov–de Gennes (BdG) formalism to treat the proximity-induced superconducting pairing [32], which yields the following model Hamiltonian:
H = 1 2 k Ψ k H 0 ( k ) Δ ̲ ̲ Δ ̲ ̲ H 0 * ( k ) = H BdG Ψ k ,
with the Nambu spinor Ψ k = ( c k , c k T ) with c k = ( c k A , c k B , c k A , c k B ) . Here, c k α ( c k α ) are creation (annihilation) operators that form a basis for H 0 ( H 0 * ) with degrees of freedom α , β { A , B , A , B } .
We consider conventional s-wave pairing (induced by the superconductor on top of the TI nanoribbon) in the BdG formalism, which is given by a momentum-independent pairing
Δ ̲ ̲ = i { Δ } s y { Δ } s y ,
with complex-valued pairing potential Δ = { Δ } + i { Δ } . Combining Equations (1)–(3), we obtain the following 8-by-8 BdG-Hamiltonian matrix:
H BdG ( k ) = [ ϵ ( k ) μ ] τ z + M ( k ) τ z σ z + A ( k y s x k x τ z s y ) σ x + A z k z τ z σ y = { Δ } τ y s y i { Δ } τ x s y ,
with τ c ( c { x , y , z } ) Pauli matrices acting on the particle–hole subspace.
Note that the superconducting proximity effect encompasses several aspects, such as an induced superconducting gap, an induced pairing potential, and induced particle–hole correlations [33]. For our simulation approach, based on the BdG Hamiltonian in Equation (4), we do not explicitly include the superconductor on top and only consider the pairing potential that is induced at the interface (see Ref. [34] for an explicit derivation of such a pairing term at the interface) as input. With this input, we can calculate the surface-state quasiparticle spectrum at low energies, with induced spectral (superconducting) gap, as well as particle–hole correlations throughout the complete TI nanoribbon, for example [35].

2.3. Tight-Binding Model

For simulating a proximitized TI nanoribbon, we discretize the BdG Hamiltonian of Equation (4) via the standard procedure onto a regular cubic grid with lattice constant a = 1 nm , resulting in a tight-binding model Hamiltonian with on-site and nearest-neighbor hopping matrices, H onsite and H hop x , y , z , respectively.
For the BdG-Hamiltonian H BdG ( k ) described in Equation (4), this results in the following matrices:
H onsite = C 0 4 C + 2 C z a 2 μ τ z + M 0 4 M + 2 M z a 2 τ z σ z , = { Δ } τ y s y i { Δ } τ x s y H hop x = C a 2 τ z + M a 2 τ z σ z + i A 2 a τ z s y σ x ,
with comparable hopping matrices along the y and z directions.
All the simulation results presented in this work are obtained with this tight-binding model (see Appendix A for details on the implementation of the simulation approach and retrieval of the spectral gap). For our setup, we consider a nanoribbon with infinite length L or L = 1 μ m along the x-direction, and a square cross section (width W = 10 nm along y, height H = 10 nm along z, perimeter P = 2 W + 2 H = 40 nm ) that is proximitized by an s-wave superconductor covering its top surface (see Figure 1a). Because the TI is not an intrinsic superconductor, the pairing potential decays quickly away from the TI-superconductor interface, i.e., over atomic distances [35]. It is therefore reasonable to assume a nonzero Δ (here, we consider Δ to be real, without loss of generality, and equal to 5 meV ) only on the topmost layer of the TI nanoribbon lattice model (∼ 1 nm thick) that is considered to be in direct contact with the superconductor, while Δ = 0 elsewhere in the lattice.
For the orbital effect of an external magnetic field, we consider the Peierls substitution method for the hopping terms:
t i j t i j exp i q r i r j d r · A ( r ) ,
Here, t i j represents the hopping matrix from site i with position r i to j with position r j , which is modified by Peierls substitution with a phase depending on the vector potential A with corresponding external magnetic field B = × A , reduced Planck’s constant , and the charge q = e (with e the elementary charge) of the charge carrier (either a particle or a hole). We consider a constant external magnetic field oriented along the nanoribbon (in the frame of reference considered here, along the x-direction) with A = ( 0 , | B | ( z H ) , 0 ) . We consider this vector potential such that it vanishes on the topmost layer of the nanoribbon ( z = H ) and is compatible with A = B = 0 for z > H , avoiding any supercurrent | Δ | A in our description [32]. These assumptions correspond to a simplified treatment of the experimental setup: a vector magnetic field that vanishes completely inside the superconductor on top of the TI nanoribbon while neglecting any shielding current.

3. Spectral Gap

In Figure 1c, we present the gap in the quasiparticle spectrum of a top-proximitized TI nanoribbon as a function of magnetic flux piercing the cross section of the ribbon and of chemical potential (note that μ = 0 corresponds to the position of the Dirac point of the TI surface-state Dirac cone). The parameter space is divided into different gapped regions that are separated by gapless phase boundaries or regions. In general, the gap lies somewhere between zero and | Δ | , with Δ the superconducting pairing potential considered on the top surface of the TI nanoribbon (in direct contact with the proximitizing superconductor), which provides a natural upper bound for the proximity-induced superconducting gap. Only near μ = 0 and integer multiples of the flux quantum does the gap exceed | Δ | . Here, the evaluated spectral gap is a trivial insulating gap due to confinement quantization ∼ π A / P , rather than a proximity-induced superconducting gap.
The spectral gap is either topologically trivial or nontrivial, and the top-proximitized TI nanoribbon only forms a quantum wire with unpaired MBSs in the nontrivial regime. The nature of the gap can be determined with the following Z 2 topological invariant [11],
M = sign { Pf [ H TINR ( k = 0 ) ] Pf [ H TINR ( k = π / a ) ] } ,
with Pf short for the Pfaffian and H TINR ( k ) the tight-binding model Hamiltonian over the cross section of the TI nanoribbon with wave number k along the direction of the ribbon. The trivial and nontrivial regions are indicated by color and are in good qualitative agreement with the diamond-tiled phase diagram that can be obtained analytically for a cylindrical TI nanowire model with Δ = 0 [23]. A (nonproximitized) cylindrical TI nanowire has the following surface-state (particle) spectrum:
E l ( k ) = ± v D k 2 + ( 2 π l + 1 / 2 η ) 2 / P 2 μ ,
with k the wave number (for propagation along the nanowire), l = 0 , ± 1 , the quantum number for quantized angular momentum, v D the Dirac velocity of the Dirac cone, η = Φ / Φ 0 the total magnetic flux piercing the nanoribbon cross section in multiples of flux quanta ( Φ 0 h / e ), and P = 2 π R the diameter of the cylindrical nanowire with radius R. From this expression, it can be seen that the spectrum is a subband-quantized Dirac cone that is flux quantum-periodic. By evaluating the number of Fermi points ν with k > 0 (or k < 0 ), i.e., the number of forward (or backward)-propagating surface states at zero energy from the different subbands that cross the chemical potential, the topological invariant above can be obtained in an alternative way by evaluating M = ( 1 ) ν . In other words, the system becomes topologically nontrivial when there is an odd number of such Fermi points and corresponding propagating modes. Each diamond in the phase diagram represents a bounded region with a given number of Fermi points, which is even for the diamonds in grayscale, and odd for the diamonds in redscale. Hence, for a piercing magnetic flux that is a half-integer multiple of one flux quantum, the system always has an odd number of Fermi points and remains in the nontrivial regime for all values of the chemical potential μ (within the TI nanoribbon bulk gap, as we are only considering the topological surface states inside the bulk gap). Conversely, the system is always in the trivial regime without an external magnetic field, or when the piercing magnetic flux is an integer multiple of Φ 0 .
Note that Fermi points, i.e., zero-energy surface states, can only be considered in general for Δ = 0 , as the states can otherwise gap out around zero energy. Therefore, we resort to the more general topological invariant in Equation (7) and the perfect diamond tiling gets slightly deformed, with different diamonds in the phase diagram becoming connected. Further note that the diamond height (its extent as a function of chemical potential) is equal to the subband energy spacing and smaller than 2 π A / P (with A = v D ), which is the expected spacing, based on the cylindrical wire model when substituting the diameter with the perimeter of the square cross section. This can also be seen in the surface-state spectrum presented in Figure 2a and has been reported before [30]. It may originate from the pile-up of wave function density near the corners, which appears to increase the effective perimeter of the cross section.
An interesting finding is that the spectral gap in the nontrivial region is maximal for the diamond closest to the Dirac point ( μ = 0 ), which corresponds to the region with a single Fermi point, and remains large for the diamonds at higher μ for a piercing flux close to a half-integer flux quantum (see Figure 1d). The size of this spectral gap agrees well with the perturbative estimate E gap ψ p Δ ̲ ̲ ψ h ( W / P ) Δ Δ / 4 for our TI nanoribbon with square cross section, with ψ p and ψ h particle and hole surface states at the Fermi point with Δ = 0 . The reduction factor 1 / 4 originates from the ratio of the section of the perimeter with nonzero Δ (only the top surface) to the complete perimeter, which gets enveloped by the particle and hole surface states.
In addition to the separation into trivial (grayscale) and nontrivial (redscale) regions, the spectral gap also reveals different gapped phases within the trivial and nontrivial regions themselves, separated by gapless phase boundaries (note that gap closings without a change of topological invariant were also reported in Ref. [34]). This suggests that the regions with and without unpaired MBSs subdivide further into different phases with additional distinct properties. As some of the regions stretch out over multiple diamonds, we can already rule out that the properties are strictly related to the number of Fermi points when Δ = 0 .
To reveal the properties of the different phases in the phase diagram, based on the spectral gap, we take a closer look at four points (indicated by ■, ★, •, and ▲), which lie in different diamonds or regions separated by a gapless boundary in Figure 1c. Their quasiparticle spectrum is presented in Figure 2b. In general, we see the expected number of Fermi points of the different subbands, based on the diamond to which the point belongs (a single Fermi point in the diamond of ■, two Fermi points in the diamond of •, and three Fermi points in the diamond of ★ and ▲), and local minima of the proximity-induced superconducting quasiparticle gap forming near them. Interestingly, ★ belongs to the same phase that stretches out over the complete single-Fermi point diamond below, in which ■ also lies, while ▲ sits in the same diamond and represents a different phase, separated from ★ by a gapless boundary. The qualitative difference between the two phases in the quasiparticle can be narrowed down to the number of local minima that can be attributed to a single spinless channel. In the diamond with three Fermi points where ★ and ▲ lie, there is either a single such local minimum or three of them, depending on the relative positioning of the Fermi points in reciprocal space. For ★, two Fermi points overlap in momentum, which effectively turns these channels into a trivial spinful channel. Hence, the number of spinless channels is reduced to one, which is the same numbers as in the diamond below. This spectral property has consequences for the formation of MBSs, as will be discussed in the following section.

4. Robust and Fragile Majorana Bound States

In Figure 3a, we present the low-energy quasiparticle spectrum as a function of flux of a top-proximitized TI nanoribbon with finite length. In this way, we also reveal the states that are localized at the ends of the ribbon. We fix the chemical potential to two different values to explore the different trivial and nontrivial phases, as discussed in the section above. Near Φ = n Φ 0 ( n Z ), there is a completely trivial insulating phase without any subgap states. For other values of the flux, however, subgap states appear in the spectrum, even in regions that are trivial according to the topological invariant in Equation (7), and they are localized at both ends of the nanoribbon (see Figure 4a).
Close to Φ = ( n + 1 / 2 ) Φ 0 and for the complete nontrivial diamonds nearest to the Dirac point (the connected region containing ■ and ★, for example), there are only two subgap states. These can be identified as a pair of MBSs forming at opposite ends of the top-proximitized TI nanoribbon. In other words, this phase corresponds to the conventional Majorana quantum wire system. This phase and its MBSs are robust against local disorder, as can be seen in Figure 3b,c, where it is presented how the low-energy spectra of Figure 3a are affected by increasing electrostatic disorder throughout the nanoribbon.
In the region with •, four subgap states form, two on each end of the ribbon. In this case, the TI nanoribbon effectively has two independent spinless channels (see previous section and Figure 2), which both give rise to the formation of a pair of MBSs forming at opposite ends of the nanoribbon. However, as there are two MBSs on each wire end, these MBSs can couple to the other MBS on the same end when they are exposed to local electrostatic disorder. Hence, in the presence of disorder, these MBSs hybridize into non-self charge-conjugate bound states with finite energy. We refer to such MBSs as fragile MBSs. In contrast, a single unpaired MBS can only suffer from hybridization with the MBS on the opposite end of the ribbon, which can be suppressed by making the proximitized section significantly longer than the MBS localization length v D / E gap . Therefore, we refer to it as a robust MBS pair. The low-energy spectrum of ▲ shows six subgap states, which can be identified as three pairs of MBSs that form at opposite ends of the TI nanoribbon. When electrostatic disorder is present, two out of the three MBSs can hybridize locally and are thus fragile, while a single robust pair should remain protected as long as the TI nanoribbon remains in the nontrivial regime (similar to what happens in a tri-junction of Majorana wires [16]). When the electrostatic disorder strength becomes of the order of the diamond height 2 π A / P , the TI nanoribbon is not guaranteed to remain in in the same phase when the flux is not an integer ( Φ = n Φ 0 ) or half-integer [ Φ = ( n + 1 / 2 ) Φ 0 ] multiple of the flux quantum. In this strongly disordered regime, the spectral gap will fluctuate strongly along the disordered nanoribbon and locally cross gapless phase boundaries. Because of this, many states with near-zero energies can appear, which renders it difficult to interpret the spectrum in terms of end-localized MBS pairs.
Now, we can put all these findings together with the phase diagram and spectral properties obtained in the previous section. It becomes clear that the different gapped phases can be identified by the number of MBS pairs forming at opposite ends of the proximitized TI nanoribbon, with the number always being even (odd) in the trivial (nontrivial) regime. When there is some amount of electrostatic disorder in the proximitized TI nanoribbon, however, it is expected that any even number of MBSs will hybridize locally with the other MBSs at the same end, while a single pair of MBSs should remain immune from hybridization in a nontrivial regime with an odd number of pairs in total. This single pair will then survive near zero energy at opposite ends of the proximitized TI nanoribbon. We can thus identify and label the phases in Figure 1c (also see Figure 1b) by their number of robust MBS pairs (zero or one) and fragile MBS pairs (an even number). As the phase with a single robust MBS pair and zero fragile pairs stretches out over a large chemical potential window near Φ = ( n + 1 / 2 ) Φ 0 , it is the most robust phase with respect to electrostatic disorder. It does not suffer from fragile MBS, nor from strong fluctuations of the spectral gap. In the subsection below, we discuss characteristic signatures of these different phases and their MBSs in a tunneling spectroscopy setup.

4.1. Tunneling Spectroscopy

In this subsection, we consider the tunneling conductance of a metallic tunneling probe that is attached to an uncovered end of a proximitized TI nanoribbon (see Appendix A for details), as a function of piercing magnetic flux and energy of the injected carriers from the tunneling probe (corresponding to bias voltage across the tunneling junction). Note that, in the experimental setup, the uncovered end should be shorter than the induced coherence length v D / E gap [9] for obtaining a hard proximity-induced gap and clearly revealing the subgap states. The results of these simulations are shown in Figure 3d–f for the same TI nanoribbon and parameters as in Figure 3a–c, and in Figure 4b for fixed values of the piercing magnetic flux. We consider conductance normalized to G 0 e 2 / h .
Overall, we find that the tunneling conductance reveals the subgap spectrum of the top-proximitized TI nanoribbon without disorder well, with a quantized conductance peak of 2 e 2 / h near zero bias (corresponding to perfect Andreev reflection) when a single unpaired MBS is localized on the side of the tunneling probe. Due to the finite length of the ribbon, there is hybridization of MBS pairs across the length of the ribbon, resulting in a splitting of the conductance peak away from zero bias. This splitting is modulated by the piercing magnetic flux.
When electrostatic disorder is introduced, the tunneling conductance of the phase with only a single (robust) MBS pair remains qualitatively the same. The zero-bias conductance peak is wider in the case of strong disorder, but it remains quantized and pinned at zero energy (see Figure 4b). The tunneling conductance of the phases with (fragile) MBS pairs gets heavily affected by disorder, however, and a quantized conductance peak cannot be easily identified. Instead, the conductance near zero bias displays irregular signatures that are very sensitive to the bias and piercing flux. This can be expected as the subgap spectrum itself is heavily affected by the disorder. Hence, it will be harder to identify phases with fragile MBSs via tunneling spectroscopy, and to determine whether the phase lies in the trivial or nontrivial regime, in particular when the electrostatic disorder strength is of the order of the subband spacing or larger.

5. Discussion

It is important to note here that the realization of fully gapped topological superconductivity in proximitized TI nanoribbons has been considered before with comparable simulation approaches [23,24,25,26,34,36]. It was pointed out by de Juan et al. in Ref. [25] that some form of transverse asymmetry is required to open up a superconducting gap in the TI nanoribbon surface-state quasiparticle spectrum. Without asymmetry, the induced gap is expected to vanish ( E gap ψ p Δ ̲ ̲ ψ h 0 ) because of a mismatch of quantized transverse momentum. Transverse asymmetry can be induced by a superconducting vortex enveloping the nanoribbon [25], by electrostatic gating [26], or by considering more sophisticated hybrid structures with multiple superconductors and gates, which may also enhance the proximity-induced gap away from the Dirac point [34], for example. Our results, however, suggest that the strong decay of the pairing potential away from the intrinsic superconductor is already sufficient to realize the required transverse asymmetry when only bringing one of the side surfaces (e.g., the top surface) in direct contact with the superconductor. In this way, a sizeable proximity-induced gap, ( W / P ) | Δ | or ( H / P ) | Δ | , can be naturally achieved near a half-flux quantum of piercing magnetic flux, especially when the Fermi level is close to the Dirac point [34]. Hence, our findings suggest that a rather straightforward device layout (a TI nanoribbon onto which superconducting material is deposited) should already be ideally suited for realizing the topologically nontrivial regime with well-separated unpaired MBSs.
Furthermore, we note that multiple (fragile) MBSs are not unique to top-proximitized TI nanoribbons. They also appear in proximitized semiconductor nanowires with a multi-subband treatment [37], for example. In that case, however, the nontrivial regions with a different number of Fermi points are disconnected, such that electrostatic disorder can more easily push the system into a trivial regime [38].
For future work, the consideration of material and sample-specific parameters for the, e.g., TI model Hamiltonian, ribbon dimensions, and disorder strength would be interesting to explore. Equally relevant is the consideration of more complicated MBS architectures that allow for braiding, with multiple proximitized nanoribbons with different orientations as building blocks (e.g., a Y-junction [24]). For such structures, the orbital effect from an external magnetic field that is misaligned with one of the ribbons should also be considered. This has been shown to induce a steering effect in TI nanoribbon structures [30,39], and also affects the topological phase in semiconductor nanowires, for example [40,41,42,43].
Regarding the experimental feasibility, we expect that the phase diagram presented here can be resolved in state-of-the-art TI nanoribbon samples. Quasi-ballistic transport of topological surface states [44,45] (also in combination with tunability of the Fermi level with respect to the Dirac point via electrostatic gating [46,47,48,49,50]), surface-state subband quantization [51], as well as proximity-induced superconductivity [8], have all been demonstrated in TI nanowires or ribbons. Aside from electrostatic gating, heterostructure engineering can be considered to tune the intrinsic Fermi level within a few meV from the Dirac point [52,53].
Finally, we comment on the consideration of tunneling spectroscopy to probe the different topological phases of top-proximitized TI nanoribbons. From alternative MBS platforms (in particular, semiconductor nanowires), we know that a (quantized) zero-bias conductance peak as MBS signature must be considered with care, as such a signature can also have a trivial origin [54]. Nonetheless, the conditions for such false positives are less likely to appear in top-proximitized TI nanoribbons because of two important differences. First, topological TI surface states have more intrinsic robustness against disorder than the low-energy modes of semiconductor nanowires due to spin-momentum locking [55]. Second, a quantum dot at the end of the ribbon, which is one of the important mechanisms for retrieving a zero-bias conductance peak with trivial origin in semiconductor nanowires [54], is less likely to form because of the linear Dirac-cone spectrum.

6. Conclusions

With a detailed three-dimensional tight-binding model, we investigate numerically the spectral gap of three-dimensional topological insulator nanoribbons with magnetic flux-piercing and proximity-induced superconductivity, induced by a superconductor on the top surface. The spectral gap reveals a rich phase diagram as a function of flux and chemical potential, with different gapped phases with paired and unpaired Majorana bound states appearing at opposite ends of the nanoribbon. These Majorana bound states can be robust or fragile with respect to local hybridization due to electrostatic disorder. When the Fermi level in the topological insulator nanoribbon is close to the Dirac point and the piercing magnetic flux is close to half a flux quantum, we retrieve the optimal conditions for realizing fully gapped topological supperconductivity. With these conditions, there is a single pair of robust MBSs on opposite sides of the nanoribbon over an extended range of chemical potential and flux values. This phase gives rise to a robust quantized zero-bias conductance peak in tunneling spectroscopy.

Author Contributions

Conceptualization, K.M.; Software, D.H.; Investigation, D.H. and K.M.; Writing—original draft, D.H.; Writing—review & editing, D.B., M.R.C., P.S. and K.M.; Visualization, D.H. and K.M.; Supervision, P.S., D.G. and K.M.; Project administration, P.S., D.G. and K.M.; Funding acquisition, P.S. and D.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been supported financially by the German Federal Ministry of Education and Research (BMBF) via the Quantum Future project ‘MajoranaChips’ (Grant No. 13N15264) within the funding program Photonic Research Germany, by Germany’s Excellence Strategy Cluster of Excellence ‘Matter and Light for Quantum Computing’ (ML4Q) EXC 2004/1–390534769, and by the Bavarian Ministry of Economic Affairs, Regional Development and Energy within Bavaria’s High-Tech Agenda Project “Bausteine für das Quantencomputing auf Basis topologischer Materialien mit experimentellen und theoretischen Ansätzen” (Grant allocation No. 07 02/686 58/1/21 1/22 2/23).

Data Availability Statement

Data available in a publicly accessible repository.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Simulation Approach

For the tight-binding simulations, we use the Python package Kwant [56] with parallel sparse direct solver MUMPS [57] and Adaptive [58] for parameter sampling (e.g., flux and chemical potential) on a nonuniform adaptive grid. The Pfaffian is calculated numerically with the algorithm of Ref. [59]. The source code and raw data of our simulations are archived in a public repository [60].
For the calculation of the spectral gap E gap , we use an algorithm that is based on that of Ref. [41], but slightly modified to speed up the convergence. The modification consists of checking whether modes exist for Δ = 0 at E = 0 (i.e., a Fermi point). If not, the algorithm from Nijholt et al. is used [41]. If modes exist at E = 0 , we look for a local minimum of the spectral gap as a function of k around the Fermi point with an adaptive search algorithm [58]. For Δ small compared to the subband spacing (as is the case in our setup), the difference between the k value of the minimum and the subband crossing is small (see Figure A1) so the search of this local minimum quickly yields a precise result.
Figure A1. The quasiparticle spectrum of a top-proximitized TI nanoribbon near a local minimum of the gap (indicated by horizontal pink dashed lines), which develops near the Fermi point of the subband crossing of the spectrum with Δ = 0 . A narrow interval around each Fermi point (indicated by vertical red dashed lines) is considered to obtain fast convergence towards the value of the spectral gap with an adaptive search algorithm over k values in that interval.
Figure A1. The quasiparticle spectrum of a top-proximitized TI nanoribbon near a local minimum of the gap (indicated by horizontal pink dashed lines), which develops near the Fermi point of the subband crossing of the spectrum with Δ = 0 . A narrow interval around each Fermi point (indicated by vertical red dashed lines) is considered to obtain fast convergence towards the value of the spectral gap with an adaptive search algorithm over k values in that interval.
Nanomaterials 13 00723 g0a1
For the tunneling spectroscopy simulations, we consider a metal lead with the following continuum model Hamiltonian:
H metal ( k ) = 2 | k | 2 2 m * μ ,
representing a free electron gas with effective mass m * = 0.01 m e and Fermi level μ = 20 eV . These artificial model parameters allow us to consider a lattice-matched few-channel metallic tunneling probe with which the tunneling conductance can be resolved, also when exceeding 2 e 2 / h (due to Andreev reflection signatures of multiple MBSs adding up, for example). We discretize this Hamiltonian on the same cubic grid as for the TI nanoribbon model, with lattice constant a = 1 nm , and consider a metal lead with translational invariance along x and square cross section of 2-by-2 lattice sites. We attach this lead at one end of a proximitized TI nanoribbon, directly in the center of its cross section. We consider a 10-nm-long section of the TI nanoribbon that is not covered by a superconductor, separating the metallic tunneling probe from the part of the ribbon that is in direct contact with the superconductor (with Δ 0 on the top surface).

References

  1. Hasan, M.Z.; Kane, C.L. Colloquium: Topological insulators. Rev. Mod. Phys. 2010, 82, 3045–3067. [Google Scholar] [CrossRef]
  2. Fu, L.; Kane, C.L. Superconducting Proximity Effect and Majorana Fermions at the Surface of a Topological Insulator. Phys. Rev. Lett. 2008, 100, 096407. [Google Scholar] [CrossRef]
  3. Xu, J.P.; Liu, C.; Wang, M.X.; Ge, J.; Liu, Z.L.; Yang, X.; Chen, Y.; Liu, Y.; Xu, Z.A.; Gao, C.L.; et al. Artificial Topological Superconductor by the Proximity Effect. Phys. Rev. Lett. 2014, 112, 217001. [Google Scholar] [CrossRef]
  4. Wiedenmann, J.; Liebhaber, E.; Kübert, J.; Bocquillon, E.; Burset, P.; Ames, C.; Buhmann, H.; Klapwijk, T.M.; Molenkamp, L.W. Transport spectroscopy of induced superconductivity in the three-dimensional topological insulator HgTe. Phys. Rev. B 2017, 96, 165302. [Google Scholar] [CrossRef]
  5. Flötotto, D.; Ota, Y.; Bai, Y.; Zhang, C.; Okazaki, K.; Tsuzuki, A.; Hashimoto, T.; Eckstein, J.N.; Shin, S.; Chiang, T.C. Superconducting pairing of topological surface states in bismuth selenide films on niobium. Sci. Adv. 2018, 4, eaar7214. [Google Scholar] [CrossRef] [PubMed]
  6. Schüffelgen, P.; Rosenbach, D.; Li, C.; Schmitt, T.W.; Schleenvoigt, M.; Jalil, A.R.; Schmitt, S.; Kölzer, J.; Wang, M.; Bennemann, B.; et al. Selective area growth and stencil lithography for in situ fabricated quantum devices. Nat. Nanotechnol. 2019, 14, 825–831. [Google Scholar] [CrossRef]
  7. Fischer, R.; Picó-Cortés, J.; Himmler, W.; Platero, G.; Grifoni, M.; Kozlov, D.A.; Mikhailov, N.N.; Dvoretsky, S.A.; Strunk, C.; Weiss, D. 4π-periodic supercurrent tuned by an axial magnetic flux in topological insulator nanowires. Phys. Rev. Res. 2022, 4, 013087. [Google Scholar] [CrossRef]
  8. Bai, M.; Wei, X.K.; Feng, J.; Luysberg, M.; Bliesener, A.; Lippertz, G.; Uday, A.; Taskin, A.A.; Mayer, J.; Ando, Y. Proximity-induced superconductivity in (Bi1-xSbx)2Te3 topological-insulator nanowires. Commun. Mater. 2022, 3, 20. [Google Scholar] [CrossRef]
  9. McMillan, W.L. Theory of Superconductor—Normal-Metal Interfaces. Phys. Rev. 1968, 175, 559–568. [Google Scholar] [CrossRef]
  10. Potter, A.C.; Lee, P.A. Engineering a p+ip superconductor: Comparison of topological insulator and Rashba spin-orbit-coupled materials. Phys. Rev. B 2011, 83, 184520. [Google Scholar] [CrossRef] [Green Version]
  11. Kitaev, A.Y. Unpaired Majorana fermions in quantum wires. Phys.Uspekhi 2001, 44, 131–136. [Google Scholar] [CrossRef]
  12. Wilczek, F. Majorana returns. Nat. Phys. 2009, 5, 614–618. [Google Scholar] [CrossRef]
  13. Stern, A. Non-Abelian states of matter. Nature 2010, 464, 187–193. [Google Scholar] [CrossRef] [PubMed]
  14. Freedman, M.; Kitaev, A.; Larsen, M.; Wang, Z. Topological quantum computation. Bull. Am. Math. Soc. 2003, 40, 31–38. [Google Scholar] [CrossRef]
  15. Alicea, J.; Oreg, Y.; Refael, G.; von Oppen, F.; Fisher, M.P.A. Non-Abelian statistics and topological quantum information processing in 1D wire networks. Nat. Phys. 2011, 7, 412–417. [Google Scholar] [CrossRef]
  16. van Heck, B.; Akhmerov, A.R.; Hassler, F.; Burrello, M.; Beenakker, C.W.J. Coulomb-assisted braiding of Majorana fermions in a Josephson junction array. New J. Phys. 2012, 14, 035019. [Google Scholar] [CrossRef]
  17. Hyart, T.; van Heck, B.; Fulga, I.C.; Burrello, M.; Akhmerov, A.R.; Beenakker, C.W.J. Flux-controlled quantum computation with Majorana fermions. Phys. Rev. B 2013, 88, 035121. [Google Scholar] [CrossRef]
  18. Aasen, D.; Hell, M.; Mishmash, R.V.; Higginbotham, A.; Danon, J.; Leijnse, M.; Jespersen, T.S.; Folk, J.A.; Marcus, C.M.; Flensberg, K.; et al. Milestones Toward Majorana-Based Quantum Computing. Phys. Rev. X 2016, 6, 031016. [Google Scholar] [CrossRef]
  19. Zhang, Y.; Vishwanath, A. Anomalous Aharonov-Bohm Conductance Oscillations from Topological Insulator Surface States. Phys. Rev. Lett. 2010, 105, 206601. [Google Scholar] [CrossRef]
  20. Ostrovsky, P.M.; Gornyi, I.V.; Mirlin, A.D. Interaction-Induced Criticality in Z2 Topological Insulators. Phys. Rev. Lett. 2010, 105, 036803. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Bardarson, J.H.; Brouwer, P.W.; Moore, J.E. Aharonov-Bohm Oscillations in Disordered Topological Insulator Nanowires. Phys. Rev. Lett. 2010, 105, 156803. [Google Scholar] [CrossRef] [PubMed]
  22. Rosenberg, G.; Guo, H.M.; Franz, M. Wormhole effect in a strong topological insulator. Phys. Rev. B 2010, 82, 041104. [Google Scholar] [CrossRef]
  23. Cook, A.; Franz, M. Majorana fermions in a topological-insulator nanowire proximity-coupled to an s-wave superconductor. Phys. Rev. B 2011, 84, 201105. [Google Scholar] [CrossRef]
  24. Cook, A.M.; Vazifeh, M.M.; Franz, M. Stability of Majorana fermions in proximity-coupled topological insulator nanowires. Phys. Rev. B 2012, 86, 155431. [Google Scholar] [CrossRef]
  25. de Juan, F.; Bardarson, J.H.; Ilan, R. Conditions for fully gapped topological superconductivity in topological insulator nanowires. SciPost Phys. 2019, 6, 60. [Google Scholar] [CrossRef]
  26. Legg, H.F.; Loss, D.; Klinovaja, J. Majorana bound states in topological insulators without a vortex. Phys. Rev. B 2021, 104, 165405. [Google Scholar] [CrossRef]
  27. Schüffelgen, P.; Rosenbach, D.; Neumann, E.; Stehno, M.P.; Lanius, M.; Zhao, J.; Wang, M.; Sheehan, B.; Schmidt, M.; Gao, B.; et al. Stencil lithography of superconducting contacts on MBE-grown topological insulator thin films. J. Cryst. Growth 2017, 477, 183–187. [Google Scholar] [CrossRef]
  28. Zhang, H.; Liu, C.X.; Qi, X.L.; Dai, X.; Fang, Z.; Zhang, S.C. Topological insulators in Bi2Se3, Bi2Te3 and Sb2Te3 with a single Dirac cone on the surface. Nat. Phys. 2009, 5, 438–442. [Google Scholar] [CrossRef]
  29. Liu, C.X.; Qi, X.L.; Zhang, H.; Dai, X.; Fang, Z.; Zhang, S.C. Model Hamiltonian for topological insulators. Phys. Rev. B 2010, 82, 045122. [Google Scholar] [CrossRef]
  30. Moors, K.; Schüffelgen, P.; Rosenbach, D.; Schmitt, T.; Schäpers, T.; Schmidt, T.L. Magnetotransport signatures of three-dimensional topological insulator nanostructures. Phys. Rev. B 2018, 97, 245429. [Google Scholar] [CrossRef] [Green Version]
  31. Chong, Y.X.; Liu, X.; Sharma, R.; Kostin, A.; Gu, G.; Fujita, K.; Davis, J.C.S.; Sprau, P.O. Severe Dirac Mass Gap Suppression in Sb2Te3-Based Quantum Anomalous Hall Materials. Nano Lett. 2020, 20, 8001–8007. [Google Scholar] [CrossRef] [PubMed]
  32. Tinkham, M. Introduction to Superconductivity; Dover Publications: New York, NY, USA, 2004. [Google Scholar]
  33. Zhu, J.X. Bogoliubov-de Gennes Method and Its Applications; Springer International Publishing: Berlin/Heidelberg, Germany, 2016. [Google Scholar] [CrossRef]
  34. Sitthison, P.; Stanescu, T.D. Robustness of topological superconductivity in proximity-coupled topological insulator nanoribbons. Phys. Rev. B 2014, 90, 035313. [Google Scholar] [CrossRef]
  35. Chiu, C.K.; Cole, W.S.; Sarma, S.D. Induced spectral gap and pairing correlations from superconducting proximity effect. Phys. Rev. B 2016, 94, 125304. [Google Scholar] [CrossRef]
  36. Legg, H.F.; Loss, D.; Klinovaja, J. Metallization and proximity superconductivity in topological insulator nanowires. Phys. Rev. B 2022, 105, 155413. [Google Scholar] [CrossRef]
  37. Stanescu, T.D.; Lutchyn, R.M.; Sarma, S.D. Majorana fermions in semiconductor nanowires. Phys. Rev. B 2011, 84, 144522. [Google Scholar] [CrossRef]
  38. Aghaee, M.; Akkala, A.; Alam, Z.; Ali, R.; Ramirez, A.A.; Andrzejczuk, M.; Antipov, A.E.; Aseev, P.; Astafev, M.; Bauer, B.; et al. InAs-Al Hybrid Devices Passing the Topological Gap Protocol. arXiv 2022, arXiv:2207.02472. [Google Scholar] [CrossRef]
  39. Kölzer, J.; Moors, K.; Jalil, A.R.; Zimmermann, E.; Rosenbach, D.; Kibkalo, L.; Schüffelgen, P.; Mussler, G.; Grützmacher, D.; Schmidt, T.L.; et al. In-plane magnetic field-driven symmetry breaking in topological insulator-based three-terminal junctions. Commun. Mater. 2021, 2, 116. [Google Scholar] [CrossRef]
  40. Lim, J.S.; Serra, L.; López, R.; Aguado, R. Magnetic-field instability of Majorana modes in multiband semiconductor wires. Phys. Rev. B 2012, 86, 121103. [Google Scholar] [CrossRef]
  41. Nijholt, B.; Akhmerov, A.R. Orbital effect of magnetic field on the Majorana phase diagram. Phys. Rev. B 2016, 93, 235434. [Google Scholar] [CrossRef] [Green Version]
  42. Kiczek, B.; Ptok, A. Influence of the orbital effects on the Majorana quasi-particles in a nanowire. J. Phys. Condens. Matter 2017, 29, 495301. [Google Scholar] [CrossRef]
  43. Winkler, G.W.; Antipov, A.E.; van Heck, B.; Soluyanov, A.A.; Glazman, L.I.; Wimmer, M.; Lutchyn, R.M. Unified numerical approach to topological semiconductor-superconductor heterostructures. Phys. Rev. B 2019, 99, 245408. [Google Scholar] [CrossRef]
  44. Dufouleur, J.; Veyrat, L.; Teichgräber, A.; Neuhaus, S.; Nowka, C.; Hampel, S.; Cayssol, J.; Schumann, J.; Eichler, B.; Schmidt, O.G.; et al. Quasiballistic Transport of Dirac Fermions in a Bi2Se3 Nanowire. Phys. Rev. Lett. 2013, 110, 186806. [Google Scholar] [CrossRef] [PubMed]
  45. Kim, H.S.; Shin, H.S.; Lee, J.S.; Ahn, C.W.; Song, J.Y.; Doh, Y.J. Quantum electrical transport properties of topological insulator Bi2Te3 nanowires. Curr. Appl. Phys. 2016, 16, 51–56. [Google Scholar] [CrossRef]
  46. Cho, S.; Dellabetta, B.; Zhong, R.; Schneeloch, J.; Liu, T.; Gu, G.; Gilbert, M.J.; Mason, N. Aharonov–Bohm oscillations in a quasi-ballistic three-dimensional topological insulator nanowire. Nat. Commun. 2015, 6, 7634. [Google Scholar] [CrossRef]
  47. Jauregui, L.A.; Pettes, M.T.; Rokhinson, L.P.; Shi, L.; Chen, Y.P. Magnetic field-induced helical mode and topological transitions in a topological insulator nanoribbon. Nat. Nanotechnol. 2016, 11, 345–351. [Google Scholar] [CrossRef]
  48. Ziegler, J.; Kozlovsky, R.; Gorini, C.; Liu, M.H.; Weishäupl, S.; Maier, H.; Fischer, R.; Kozlov, D.A.; Kvon, Z.D.; Mikhailov, N.; et al. Probing spin helical surface states in topological HgTe nanowires. Phys. Rev. B 2018, 97, 035157. [Google Scholar] [CrossRef]
  49. Kim, H.S.; Hwang, T.H.; Kim, N.H.; Hou, Y.; Yu, D.; Sim, H.S.; Doh, Y.J. Adjustable Quantum Interference Oscillations in Sb-Doped Bi2Se3 Topological Insulator Nanoribbons. ACS Nano 2020, 14, 14118–14125. [Google Scholar] [CrossRef]
  50. Rosenbach, D.; Moors, K.; Jalil, A.R.; Kölzer, J.; Zimmermann, E.; Schubert, J.; Karimzadah, S.; Mussler, G.; Schüffelgen, P.; Grützmacher, D.; et al. Gate-induced decoupling of surface and bulk state properties in selectively-deposited Bi2Te3 nanoribbons. SciPost Phys. Core 2022, 5, 017. [Google Scholar] [CrossRef]
  51. Münning, F.; Breunig, O.; Legg, H.F.; Roitsch, S.; Fan, D.; Rößler, M.; Rosch, A.; Ando, Y. Quantum confinement of the Dirac surface states in topological-insulator nanowires. Nat. Commun. 2021, 12, 1038. [Google Scholar] [CrossRef]
  52. Kellner, J.; Eschbach, M.; Kampmeier, J.; Lanius, M.; Młyńczak, E.; Mussler, G.; Holländer, B.; Plucinski, L.; Liebmann, M.; Grützmacher, D.; et al. Tuning the Dirac point to the Fermi level in the ternary topological insulator (Bi1-xSbx)2Te3. Appl. Phys. Lett. 2015, 107, 251603. [Google Scholar] [CrossRef]
  53. Eschbach, M.; Młyńczak, E.; Kellner, J.; Kampmeier, J.; Lanius, M.; Neumann, E.; Weyrich, C.; Gehlmann, M.; Gospodarič, P.; Döring, S.; et al. Realization of a vertical topological p–n junction in epitaxial Sb2Te3/Bi2Te3 heterostructures. Nat. Commun. 2015, 6, 8816. [Google Scholar] [CrossRef]
  54. Pan, H.; Das Sarma, S. Physical mechanisms for zero-bias conductance peaks in Majorana nanowires. Phys. Rev. Res. 2020, 2, 013377. [Google Scholar] [CrossRef]
  55. Dufouleur, J.; Veyrat, L.; Dassonneville, B.; Xypakis, E.; Bardarson, J.H.; Nowka, C.; Hampel, S.; Schumann, J.; Eichler, B.; Schmidt, O.G.; et al. Weakly-coupled quasi-1D helical modes in disordered 3D topological insulator quantum wires. Sci. Rep. 2017, 7, 45276. [Google Scholar] [CrossRef]
  56. Groth, C.W.; Wimmer, M.; Akhmerov, A.R.; Waintal, X. Kwant: A software package for quantum transport. New J. Phys. 2014, 16, 063065. [Google Scholar] [CrossRef]
  57. Amestoy, P.R.; Duff, I.S.; L’Excellent, J.Y.; Koster, J. A fully asynchronous multifrontal solver using distributed dynamic scheduling. SIAM J. Matrix Anal. Appl. 2001, 23, 15–41. [Google Scholar] [CrossRef]
  58. Nijholt, B.; Weston, J.; Hoofwijk, J.; Akhmerov, A. Adaptive: Parallel active learning of mathematical functions. Zenodo 2019, 1182437. [Google Scholar] [CrossRef]
  59. Wimmer, M. Efficient Numerical Computation of the Pfaffian for Dense and Banded Skew-Symmetric Matrices. ACM Trans. Math. Softw. 2012, 38, 1–17. [Google Scholar] [CrossRef] [Green Version]
  60. Heffels, D.; Burke, D.; Connolly, M.R.; Schüffelgen, P.; Grützmacher, D.; Moors, K. Replication Data for: Robust and Fragile Majorana Bound States in Proximitized Topological Insulator Nanoribbons. Available online: https://data.fz-juelich.de/dataset.xhtml?persistentId=doi:10.26165/JUELICH-DATA/LZOYYD (accessed on 1 December 2022).
Figure 1. (a) A 3D TI nanoribbon that is proximitized by an s-wave superconductor via its top surface and pierced by a magnetic flux originating from an external magnetic field oriented along the ribbon. With this setup, MBSs form at opposite ends of the nanoribbon. The setup has different phases with different numbers of robust and fragile MBSs, which are shown schematically in (b); (c) the spectral gap in the proximitized TI nanoribbon as a function of magnetic flux and chemical potential. The grayscale (redscale) colormap indicates whether the gap is topologically trivial (nontrivial). The phase boundaries between the trivial and nontrivial regions are indicated with green dashed lines. The number of robust plus the number of fragile MBS pairs is indicated for some representative phases; (d) a line cut of the spectral gap shown in (c) as a function of chemical potential with magnetic flux equal to half a flux quantum (indicated by a vertical black dotted line in (c)). The gap of each surface-state subband that has a Fermi point (FP) when Δ = 0 is presented separately, as well as the gap of the lowest subband without Fermi point, evaluated at k = 0 .
Figure 1. (a) A 3D TI nanoribbon that is proximitized by an s-wave superconductor via its top surface and pierced by a magnetic flux originating from an external magnetic field oriented along the ribbon. With this setup, MBSs form at opposite ends of the nanoribbon. The setup has different phases with different numbers of robust and fragile MBSs, which are shown schematically in (b); (c) the spectral gap in the proximitized TI nanoribbon as a function of magnetic flux and chemical potential. The grayscale (redscale) colormap indicates whether the gap is topologically trivial (nontrivial). The phase boundaries between the trivial and nontrivial regions are indicated with green dashed lines. The number of robust plus the number of fragile MBS pairs is indicated for some representative phases; (d) a line cut of the spectral gap shown in (c) as a function of chemical potential with magnetic flux equal to half a flux quantum (indicated by a vertical black dotted line in (c)). The gap of each surface-state subband that has a Fermi point (FP) when Δ = 0 is presented separately, as well as the gap of the lowest subband without Fermi point, evaluated at k = 0 .
Nanomaterials 13 00723 g001
Figure 2. (a) The subband-quantized surface-state spectrum of a TI nanoribbon pierced by half a magnetic flux quantum. The two chemical potential values indicated with line cuts in Figure 1c and considered in (b) are indicated with horizontal lines; (b) the low-energy quasiparticle spectrum of a top-proximitized TI nanoribbon with infinite length for different combinations of piercing magnetic flux and chemical potential, as indicated in Figure 1c, with the color indicating the particle–hole (p-h) mixing (blue for hole, red for particle). The spectrum for Δ = 0 is presented with thin black lines and the extent of the spectral gap is indicated with horizontal pink dashed lines.
Figure 2. (a) The subband-quantized surface-state spectrum of a TI nanoribbon pierced by half a magnetic flux quantum. The two chemical potential values indicated with line cuts in Figure 1c and considered in (b) are indicated with horizontal lines; (b) the low-energy quasiparticle spectrum of a top-proximitized TI nanoribbon with infinite length for different combinations of piercing magnetic flux and chemical potential, as indicated in Figure 1c, with the color indicating the particle–hole (p-h) mixing (blue for hole, red for particle). The spectrum for Δ = 0 is presented with thin black lines and the extent of the spectral gap is indicated with horizontal pink dashed lines.
Nanomaterials 13 00723 g002
Figure 3. The (ac) low-energy quasiparticle spectrum of a top-proximitized TI nanoribbon with finite length as a function of piercing magnetic flux Φ , and the (df) tunneling conductance as a function of piercing magnetic flux and energy E (or bias over the tunneling junction). The pink lines correspond to the spectral gap in the TI nanoribbon with infinite length and without disorder. The top (bottom) panel of each subfigure corresponds to the upper (lower) horizontal line cut shown in Figure 1c at μ 1.3 × 2 π A / P ( μ 0.75 × 2 π A / P ). The results are obtained (a,d) without disorder, and (b,c,e,f) with disorder (considering S dis = 0.1 × 2 π A / P in (b,e), and S dis = 2 π A / P in (c,f), with the different disorder strengths also indicated in Figure 1c).
Figure 3. The (ac) low-energy quasiparticle spectrum of a top-proximitized TI nanoribbon with finite length as a function of piercing magnetic flux Φ , and the (df) tunneling conductance as a function of piercing magnetic flux and energy E (or bias over the tunneling junction). The pink lines correspond to the spectral gap in the TI nanoribbon with infinite length and without disorder. The top (bottom) panel of each subfigure corresponds to the upper (lower) horizontal line cut shown in Figure 1c at μ 1.3 × 2 π A / P ( μ 0.75 × 2 π A / P ). The results are obtained (a,d) without disorder, and (b,c,e,f) with disorder (considering S dis = 0.1 × 2 π A / P in (b,e), and S dis = 2 π A / P in (c,f), with the different disorder strengths also indicated in Figure 1c).
Nanomaterials 13 00723 g003
Figure 4. The (a) one-dimensional wave function density summed over the n lowest-energy (in absolute value) solutions and (b) tunneling conductance as a function of energy of a 1 μ m-long top-proximitized TI nanoribbon with and without disorder with different disorder strengths (indicated in Figure 1c) for four different combinations of chemical potential and piercing magnetic flux (also indicated in Figure 1c). The tunneling conductance in (b) corresponds to the vertical gray dashed line cuts in Figure 3d–f. The spectral gap of the TI nanoribbon considering infinite length and no disorder is indicated by vertical pink dashed lines.
Figure 4. The (a) one-dimensional wave function density summed over the n lowest-energy (in absolute value) solutions and (b) tunneling conductance as a function of energy of a 1 μ m-long top-proximitized TI nanoribbon with and without disorder with different disorder strengths (indicated in Figure 1c) for four different combinations of chemical potential and piercing magnetic flux (also indicated in Figure 1c). The tunneling conductance in (b) corresponds to the vertical gray dashed line cuts in Figure 3d–f. The spectral gap of the TI nanoribbon considering infinite length and no disorder is indicated by vertical pink dashed lines.
Nanomaterials 13 00723 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Heffels, D.; Burke, D.; Connolly, M.R.; Schüffelgen, P.; Grützmacher, D.; Moors, K. Robust and Fragile Majorana Bound States in Proximitized Topological Insulator Nanoribbons. Nanomaterials 2023, 13, 723. https://doi.org/10.3390/nano13040723

AMA Style

Heffels D, Burke D, Connolly MR, Schüffelgen P, Grützmacher D, Moors K. Robust and Fragile Majorana Bound States in Proximitized Topological Insulator Nanoribbons. Nanomaterials. 2023; 13(4):723. https://doi.org/10.3390/nano13040723

Chicago/Turabian Style

Heffels, Dennis, Declan Burke, Malcolm R. Connolly, Peter Schüffelgen, Detlev Grützmacher, and Kristof Moors. 2023. "Robust and Fragile Majorana Bound States in Proximitized Topological Insulator Nanoribbons" Nanomaterials 13, no. 4: 723. https://doi.org/10.3390/nano13040723

APA Style

Heffels, D., Burke, D., Connolly, M. R., Schüffelgen, P., Grützmacher, D., & Moors, K. (2023). Robust and Fragile Majorana Bound States in Proximitized Topological Insulator Nanoribbons. Nanomaterials, 13(4), 723. https://doi.org/10.3390/nano13040723

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop