Interplay between pairing and correlations in spin-polarized bound states

  1. Szczepan Głodzik1,
  2. Aksel Kobiałka1,
  3. Anna Gorczyca-Goraj2,
  4. Andrzej Ptok3ORCID Logo,
  5. Grzegorz Górski4,
  6. Maciej M. Maśka2ORCID Logo and
  7. Tadeusz Domański1

1Institute of Physics, M. Curie-Skłodowska University, 20-031 Lublin, Poland
2Institute of Physics, University of Silesia, 41-500 Chorzów, Poland
3Institute of Nuclear Physics, Polish Academy of Sciences, 31-342 Kraków, Poland

4Faculty of Mathematics and Natural Sciences, University of Rzeszów, 35-310 Rzeszów, Poland

  1. Corresponding author email

This article is part of the Thematic Series "Topological materials".

Guest Editor: J. J. Palacios
Beilstein J. Nanotechnol. 2018, 9, 1370–1380. doi:10.3762/bjnano.9.129
Received 15 Dec 2017, Accepted 09 Apr 2018, Published 07 May 2018


We investigate single and multiple defects embedded in a superconducting host, studying the interplay between the proximity-induced pairing and interactions. We explore the influence of the spin–orbit coupling on energies, polarization and spatial patterns of the bound (Yu–Shiba–Rusinov) states of magnetic impurities in a two-dimensional square lattice. We also address the peculiar bound states in the proximitized Rashba chain, resembling the Majorana quasiparticles, focusing on their magnetic polarization that has been recently reported by S. Jeon et al. (Science 2017, 358, 772). Finally, we study leakage of these polarized Majorana quasiparticles into side-attached nanoscopic regions and confront them with the subgap Kondo effect near to the singlet–doublet phase transition.

Keywords: bound states in superconductors; Majorana quasiparticles; subgap Kondo effect


Magnetism is usually detrimental to superconductivity because it breaks the Cooper pairs (at the critical field strength Hc2). There are, however, a few exceptions in which these phenomena coexist, e.g., in iron pnictides [1], CeCoIn5 [2]. Also, sometimes magnetic fields induce superconductivity [3]. Plenty of other interesting examples can be found in nanoscopic systems, where magnetic impurities (dots) exhibit a more subtle relationship with the electron pairing driven by the proximity effect [4,5]. Cooper pairs easily penetrate the nanoscopic impurities, inducing the bound (Yu–Shiba–Rusinov) states that manifest the local pairing in coexistence with magnetic polarization. Such bound states have been observed in various systems [6-14]. In-gap states (appearing in pairs symmetrically around the Fermi level) can be nowadays controlled electrostatically or magnetically [12] whereas their topography, spatial extent and polarization can be precisely inspected by the state-of-art tunneling measurements [15,16].

It has been reported that adatoms deposited on a two-dimensional (2D) superconducting surface develop Yu–Shiba–Rusinov (YSR) states, extending to a dozen of intersite distances and they reveal particle–hole oscillations [11]. Bound states of these magnetic impurities in superconducting NbSe2 are characterized by the star shape [17] typical for the rotational symmetry of its triangular lattice. More complex objects, such as dimers, reveal other spatial features, showing the bonding and antibonding states [18]. In a somewhat different context it has been pointed out [19] that exchange coupling between numerous quantum defects involving their intrinsic spins can couple them ferromagnetically. This can be used (e.g., in metallic carbon nanotubes) for a robust transmission of magnetic information over large distances.

In all cases the bound YSR states are also sensitive to interactions. One of them is the spin–orbit coupling (usually meaningful at boundaries, e.g., surfaces) [20-22]. Such interaction in one-dimensional magnetic nanowires can induce the topologically nontrivial superconducting phase, in which the YSR states undergo mutation to Majorana (zero-energy) quasiparticles. Coulomb repulsion between the opposite spin electrons can bring additional important effects. In the proximitized quantum dots it can lead to a parity change (quantum phase transition) with further influence on the subgap Kondo effect (driven by effective spin-exchange coupling with mobile electrons). Furthermore, such spin exchange can be amplified by the induced electron pairing, and can have constructive influence on the Kondo effect [23,24].

We study here the polarized bound states, taking into account the spin–orbit and/or Coulomb interactions. In particular, we consider: (i) a single magnetic impurity in a 2D square lattice of a superconducting host, (ii) a nanoscopic chain of magnetic impurities on the classical superconductor (i.e., proximitized Rashba nanowire) in its topologically trivial/nontrivial superconducting phase, and (iii) a strongly correlated quantum dot side-attached to the Rashba chain, where the Kondo and the leaking Majorana quasiparticle can be confronted with each other. These magnetically polarized YSR and Majorana quasiparticles as well as the subgap Kondo effect can be experimentally verified using tunneling heterostructures with ferromagnetic lead (STM tip).

Results and Discussion

Single magnetic impurity

Let us start by considering a single magnetic impurity on the surface of an s-wave superconductor in presence of spin–orbit interactions. This situation can be modeled by the Anderson-type Hamiltonian


We describe the superconducting substrate by


where [Graphic 1] ([Graphic 2]) denotes creation (annihilation) of an electron with spin σ at the i-th site, t is a hopping integral between the nearest neighbors, μ is the chemical potential, and [Graphic 3] is the number operator. For simplicity, we assume a weak attractive potential U < 0 between itinerant electrons and treat it within the mean-field decoupling

[Graphic 4]

where [Graphic 5] is the local superconducting order parameter and n = [Graphic 6]. The Hartree term can be incorporated into the local (spin-dependent) chemical potential μ → [Graphic 7] ≡ μ − [Graphic 8]. The second term in Equation 1 refers to the local impurity


which affects the order parameter χi near the impurity site i = 0, inducing the YSR states [25,26]. In this work we focus on the magnetic term J [4,27], disregarding the potential scattering K.

The spin–orbit coupling (SOC) can be expressed by


where the vector [Graphic 9] refers to positions of the nearest neighbors of the i-th site, and [Graphic 10] = (σx, σy, σz) stands for the Pauli matrices. The unit vector [Graphic 11] shows the direction of the spin–orbit field, which can be arbitrary. Here we restrict our considerations to the in-plane [Graphic 12] = (1, 0, 0) polarization, which will be important for nontrivial superconductivity in nanowires discussed in the subsection ’Magnetically polarized Majorana quasiparticles’. The other (out-of-plane) component could eventually mix ↑ and ↓ spins [22].

Impurities break the translational invariance, therefore the pairing amplitude χi and occupancy niσ have to be determined for each lattice site individually. We can diagonalize the Hamiltonian (Equation 1) by the unitary transformation


where [Graphic 13] are quasiparticle fermionic operators with eigenvectors uinσ and vinσ. This leads to the Bogoliubov–de Gennes (BdG) equations


where Dij = δijUχi, and the single-particle term is given by

[Graphic 14]

with the spin–orbit coupling term

[Graphic 15]

Here, [Graphic 16] and [Graphic 17] (where [Graphic 18] is opposite to σ) correspond to in-plane and out-of-plane spin–orbit field, respectively, and satisfy [Graphic 19].

Solving numerically the BdG equations (Equation 6) we can determine the local order parameter χi and occupancy niσ


where f(ω) = [1 + exp(ω/kBT)]−1. In what follows, we shall inspect the spin-resolved local density of states

[Graphic 20]

For its numerical computation we replace the Dirac delta function with the Lorentzian function δ(ω) = ζ/[π(ω2 + ζ2)] with a small broadening ζ = 0.01 t. We have solved the BdG equations, considering a single magnetic impurity in a square lattice, comprising Na × Nb = 41 × 41 sites. We assumed U/t = −3, μ/t = 0, and determined the bound states for two representative values of the spin–orbit coupling λ upon varying J.

The magnetic potential has substantial influence on the local order parameter χ0. In particular, at some critical value Jc this quantity discontinuously changes its magnitude and sign (see the upper panel in Figure 1), signaling a first-order phase transition [28-30]. This quantum phase transition at Jc is an artifact of the classical spin approximation. When spin fluctuations are allowed, a Kondo-like crossover is obtained instead of a first-order phase transition [31,32]. In general, the quasiparticle spectrum at the impurity site is characterized by two bound states ±EYSR inside the gap Δ of the superconducting host (displayed in the bottom panel of Figure 1). These energies ±EYSR and the related spectral weights depend on J. At J = Jc the YSR bound states cross each other EYSR(Jc) = 0 and their crossing signifies the ground-state parity change [33] from BCS-type (spinless) to the singly occupied (spinful) configurations [8,15,21,34]. Let us remark that this quantum phase transition is also accompanied with a reversal of the YSR polarization (see bottom panel in Figure 1). A similar behavior can be observed also for multiple impurities, at several critical values of J [35].


Figure 1: The local order parameter obtained at zero temperature for weak λ/t = 0.1 (red line) and strong spin–orbit coupling λ/t = 1 (blue line). The bottom panel shows the energies and magnetic polarization ρ0↑(ω) − ρ0↓(ω) of YSR states obtained in the weak-coupling limit λ/t = 0.1.

Within the BdG approach we can inspect spatial profiles of the YSR states by integrating the spectral weights

[Graphic 21]

in the interval [Graphic 22] capturing the quasiparticles at negative/positive energies ±EYSR [36]. Figure 2 illustrates the results obtained for λ = 0 (left panel) and λ = t (right panel). We clearly notice a fourfold rotational symmetry (typical for the square lattice) and the spatial extent of YSR states reaching several sites away from the magnetic impurity. The non-vanishing difference of the spectral weight |uin|2 −|uin|2 at the positive energy ω = +EYSR and of |vin|2 −|vin|2 at the negative energy ω = −EYSR implies the effective spin-polarization of the bound states (their polarization is illustrated in the bottom panel of Figure 1).


Figure 2: Spatial profiles of the YSR states [Graphic 23] obtained for |J| < Jc in the absence of spin–orbit coupling (left panel) and for strong in-plane coupling λ = t (right panel). The spin–orbit field is chosen along the x-axis and leads to an additional imaginary hopping term along the y-axis, which elongates the YSR states in the y-direction. The impurity spin is oriented along the (0, 0, 1) direction.

For a quantitative estimation of the spatially varying magnetization (driven by the particle–hole asymmetry) we have computed the displaced moving average [Graphic 24], which corresponds to an averaged spectral weight contained in a ring of the radius r and a small half-width δr. This quantity is sensitive only to the radial distance r from the magnetic impurity, averaging the angular anisotropy. Our results, presented in Figure 3, clearly indicate the spatial particle–hole oscillations [Graphic 25] of the YSR states (compare the blue and red lines). Such particle–hole oscillations decay exponentially with r in agreement with previous studies [11,37,38]. The dominant (particle or hole) contributions to the YSR bound states are displayed by the alternating color of the background in Figure 3. We notice that the spin–orbit coupling seems to suppress these particle–hole oscillations.


Figure 3: Hole-like (blue line) and electron-like (red line) displaced moving average [Graphic 26] as a function of the radial distance r from the impurity site obtained for |J| < Jc using δr = 0.5a. The blue and red background color indicates the dominant type (hole or particle) of the YSR states at a given distance r. The left and right panels correspond to λ = 0 and λ = t, respectively.

Summarizing this section, we point out that the quantum phase transition at Jc depends on the spin–orbit coupling λ and it has experimentally observable consequences in the magnetization induced near the impurity site. For weak magnetic scattering |J| < Jc the impurity is partly screened, whereas for stronger couplings |J| > Jc the impurity polarizes its neighborhood in the direction of its own magnetic moment. Similar effects have been previously discussed in [21], but here we additionally consider the role of spin–orbit coupling. First of all, such interaction shifts the quantum phase transition (to larger values of J) and secondly it enhances the spatial extent of YSR states and gradually smoothes the particle–hole oscillations.

Magnetically polarized Majorana quasiparticles

In this section we increase the number of impurities. Let us now imagine a nanoscopic chain of magnetic impurities (for instance Fe atoms) deposited on the surface of a conventional s-wave superconductor. We study the magnetically polarized bound states, focusing on the proximity-induced nontrivial superconducting phase. In practice, the quasiparticle spectrum can be probed within STM-type setups, by attaching a conducting [39,40], superconducting [41], or a magnetically polarized tip [42]. We assume the spin–orbit interaction aligned perpendicularly to the wire and the magnetic field parallel to it, leading to the effective intersite pairing of identical spins and (under specific conditions) inducing zero-energy end modes resembling Majorana quasiparticles. This issue has been recently studied very intensively but here we simply focus on the spin-polarized aspects of this problem.

Due to the spin–orbit interaction, momentum and spin are no longer “good” quantum numbers. By solving the problem numerically, however, we can estimate the percentage with which the true quasiparticles are represented by the initial spin. We have recently emphasized [43], that the amplitude of intersite pairing (between identical spin electrons) differs several times for ↑ and ↓ sectors. This leads to an obvious polarization of the YSR and Majorana quasiparticles (the latter appearing near the nanochain edges).

Let us consider the STM-type geometry relevant to the recent experimental situation addressed by A. Yazdani and co-workers [42], which can be described by the following Hamiltonian


We assume here that the STM tip describes a polarized fermion gas

[Graphic 27]

where the energy [Graphic 28] can be controlled by some finite detuning of the chemical potentials μN↑ − μN↓. Individual atoms of the nanochain are coupled with such STM tip through

[Graphic 29]

For simplicity, we assume constant couplings

[Graphic 30]

The low-energy physics of such proximitized Rashba nanowire can be described by [44]


where [Graphic 31] annihilates (creates) an electron of spin σ at site i with energy εi, and tij is the hopping integral. The effective intersite (p-wave) pairing is induced through a combined effect of the Rashba and the Zeeman terms


The proximity effect, which induces the on-site (trivial) pairing, can be modelled as [45]


with the local pairing potential Δi = ΓS/2.

Figure 4 shows evolution of the spin-dependent spectrum ρiσ(ω) as a function of a varying magnetic field. At a critical value (B ≈ 0.2) we observe the emergence of zero-energy quasiparticles, whose spectral weights strongly depend on the spin σ.


Figure 4: The effective quasiparticle spectrum ρiσ(ω) as a function of a magnetic field B aligned along the nanochain obtained for σ = ↑ (upper panel) and σ = ↓ (bottom panel). The magnetic field B is expressed in units of t/(gμB/2).

For a better understanding of the polarized zero-energy quasiparticles, we present in Figure 5 the spatial profiles of the zero-energy (Majorana) quasiparticles. As usually such quasiparticles emerge near the edges of a nanoscopic chain, practically over 10 to 15 sites (see inset). Note the substantial quantitative difference between these zero-energy quasiparticles appearing in ↑ and ↓ spin sectors. This “intrinsic polarization” of the Majorana modes has been previously suggested in [46], and recently we have proposed [47] their empirical detection by means of selective equal-spin Andreev reflection (SESAR) spectroscopy.


Figure 5: Magnetically polarized spectrum ρi,↑(ω) − ρi,↓(ω) obtained at ω = 0 for peripheral sites of the Rashba chain.

The main idea is to apply a bias voltage V between the STM tip and the superconducting substrate, inducing a charge transport that, in a subgap regime ([Graphic 32]) originates from the Andreev (particle to hole) scattering mechanism. The polarized Andreev current can be expressed by the Landauer–Büttiker formula


where transmittance is defined as

[Graphic 33]


[Graphic 34]

The anomalous Green’s functions can be computed numerically from the solution of the Bogoliubov–de Gennes equations of this model (Equation 10). The net spin current [Graphic 35] turns out to be predominantly sensitive to the Majorana end-modes. Its differential conductance [Graphic 36] can thus distinguish the polarized Majorana quasiparticle (near V = 0) from the YSR states (appearing at finite voltage).

Bound states can leak to other side-attached nanoscopic objects. This proximity effect has been also predicted for the Majorana quasiparticles by E. Vernek et al. [48] and it has been indeed observed experimentally by M. T. Deng and co-authors [49]. Inspired by this achievement, extensive studies have been carried out regarding the YSR states coalescencing into the zero-energy Majorana states in side-coupled quantum dots driven by electrostatic or magnetic fields [50-52]. This issue would be particularly important when attempting to braid the Majorana end modes, e.g., in T-shape nanowires upon turning on and off the topological superconducting phase in its segments. We briefly analyse here the polarized zero-energy Majorana modes leaking into the multi-site quantum dot (comprising ten lattice sites) side-attached to the proximitized Rashba chain discussed above.

Figure 6 displays the spatial profile of the polarized spectrum obtained at ω = 0 as a function of the gate voltage Vg, which detunes the energies Vg = εi − μ of the multi-site (1 ≤ i ≤ 10) quantum dot. For numerical calculations we used the model parameters λ = 0.15t, μ = −2t, Δi = 0.2t and B > Bc, which guarantee the Rashba chain to be in its topologically nontrivial superconducting phase, hosting the zero-energy Majorana quasiparticles (intensive black or red regions). We clearly observe that for some values of Vg these Majorana modes spread over the entire quantum dot region. By inspecting Figure 6 we furthermore notice the pronounced spatial oscillations of these zero-energy modes. In our opinion, this is a signature of a partial delocalization of the polarized Majorana quasiparticles. Surprisingly, this process seems to be less efficient in the minor spin (σ = ↓) section. This effect has to be taken into account, when designing nanostructures for a controllable spatial displacement of the Majorana modes (criticial for the realization of quantum computations with use of the Majorana-based qubits) either by electrostatic or magnetic means. Some proposals for such nanodevices have been recently discussed by several authors [52,53].


Figure 6: Leakage of the spin-polarized Majorana quasiparticles from the topological superconducting phase of the Rashba chain (i ≥ 10) onto the side-attached multi-site ([Graphic 37]) quantum dot. The upper and bottom panel show ρ(ω) at ω = 0 for ↑ and ↓ spin, respectively.

In summary of this section, we emphasize that the Majorana modes coalescing from the YSR states in the proximitized Rashba nanowire are characterized by their magnetic polarization. Indeed, such a feature has been recently observed by STM spectroscopy with use of a polarized tip [42]. We have studied here the evolution of the polarized quasiparticle states with respect to the magnetic field (Figure 4) and investigated the spatial oscillations of the Majorana zero-energy modes near the chain edges (Figure 5). Finally, we analyzed leakage of the polarized Majorana modes on the multi-site quantum dots, revealing their partial delocalization (Figure 6).

Majorana vs Kondo effect

In previous section we have discussed the polarized Majorana modes leaking into side-attached objects, such as single impurities or segments of normal nanowires. In this section we shall focus on the correlation effects [54-56], confronting the Majorana quasiparticle with the Kondo effect (both manifested at zero energy). This can be practically achieved using STM-type configurations sketched in Figure 7. In particular, we consider the subgap Kondo effect, effectively driven by the Coulomb repulsion U and coupling of the quantum dot (QD) with the normal lead ΓN in presence of electron pairing (induced via ΓS), which has a significant influence on the spin-polarized bound states of the QD. The basic mechanism of this subgap Kondo effect showing up near the quantum phase transition has been earlier considered by us in absence of the Rashba nanowire [24,57]. Our considerations can be practically verified within STM geometry [39,40] using magnetic atoms (e.g., Fe) and side-coupled nonmagnetic atoms (for instance Ag or Au) deposited on the superconducting substrate (such as Pb or Al) probed with a conducting STM tip [42].


Figure 7: Schematic illustration of the quantum dot (QD) coupled between the metallic (N) and superconducting (S) leads and hybridized with the Rashba nanowire, hosting the Majorana quasiparticles η1 and η2 at its edges.

The topological superconducting phase, hosting the Majorana modes, can be driven in semiconducting wires [58,59] or in nanochains of magnetic atoms [39-42] through nearest-neighbor equal-spin pairing. The efficiency of such p-wave pairing differs for each spin [47], giving rise to polarization of the Majorana quasiparticles, with noticeable preference for the ↑ sector (see Figure 4). In order to study the correlation effects we shall assume here a complete polarization of the Majorana quasiparticles. We thus focus, for simplicity, on the topological state originating from intersite pairing of only ↑ electrons and consider its interplay with the correlations. Let us remark, however, that the superconducting lead mixes both the QD spins with the side-attached Majorana quasiparticle [60]. In consequence we shall observe an interesting and spin-dependent relationship between the Majorana and Kondo states that could be probed by the polarized Andreev (particle-to-hole conversion) mechanism.

Our setup (Figure 7) can be described by the following Anderson-type Hamiltonian


where [Graphic 38] corresponds to the metallic electrode, [Graphic 39] refers to the s-wave superconducting substrate and the correlated QD is modeled by [Graphic 40], where ε denotes the energy level and U stands for the repulsive interaction between opposite spin electrons. The QD is coupled to both β = N,S reservoirs through [Graphic 41] and we assume a wide bandwidth limit, using the constant couplings Γβ. It can be shown [61-64] that for energies [Graphic 42] the super-conducting electrode induces the static on-dot pairing

[Graphic 43]

Taking into account the finite magnitude of superconducting gap [50] does not affect the main conclusions of our study.

The effective Majorana modes of the nanowire can be modeled by [65]

[Graphic 44]

where [Graphic 45] are Hermitian operators and εm corresponds to an overlap between Majoranas. We recast these operators by the standard fermionic ones [66] [Graphic 46] and [Graphic 47]. Finally, the Hamiltonian of Equation 15 simplifies to


with the auxiliary coupling [Graphic 48]. The subgap Kondo physics originates in this model from the Coulomb term [Graphic 49] and the effective spin-exchange interactions due to [Graphic 50]. It has been shown [23,24] that under specific conditions the on-dot pairing can cooperate with the subgap Kondo effect. This particular situation occurs only near the quantum phase transition.

Let us examine how the subgap Kondo effect gets along with the Majorana mode. Earlier studies of the correlated quantum dot coupled to both normal (conducting) electrodes indicated that the side-attached Rashba chain leads to a competition between the Kondo and Majorana states [67-72]. For a sufficiently long wire (εm = 0) the Kondo effect persists only in the spin-channel ↓, whereas for ↑ electrons there appears a dip in the spectral density at ω = 0. The resulting tunneling conductance is then partly reduced (from the perfect value 2e2/h) to the fractional value 3e2/2h [67,68,71-73]. In contrast, for the short Rashba wires (with εm ≠ 0) the Kondo physics persists in both spin channels.

In our present setup (Figure 7) the correlated quantum dot is between the metallic and superconducting reservoirs, therefore the Kondo effect is additionally affected by on-dot pairing. Its influence is mainly controlled by the ratio US and partly by the level ε, determining whether the QD ground state is in the spinful or spinless configuration [23,24,62,64,74]. Obviously the latter one cannot be screened. For instance, for the half-filled QD (ε = −U/2) the spinful (doublet) configuration occurs in the regime U ≥ ΓS.

For studying the correlations we adopt perturbative treatment of the Coulomb potential, treating it self-consistently to the second order in the normal and anomalous channels [62,75]. Specific expressions have been provided by us in [24]. Figure 8 shows the spectral function ρσ(ω) for both spins obtained at zero temperature for the Coulomb potential U, covering the (spinless) singlet and (spinful) doublet configurations. In the weak interaction regime we observe appearance of two YSR states. For U ≈ ΓS these peaks merge, signaling the quantum phase transition. The Kondo effect shows up only in the correlated limit (U > ΓS), but its spectroscopic signatures are qualitatively different for each of the spins. Leakage of the Majorana quasiparticle suppresses the low-energy states of ↑ electrons. We notice that the initial density (for tm = 0) is reduced by half, whereas we observe a constructive influence of the Majorana quasiparticle on opposite-spin ↓ electrons.


Figure 8: The polarized spectral function ρσ(ω) obtained at zero temperature for the half-filled QD (ε = −U/2), ΓS = 2ΓN, tm = 0.1ΓN and several values of the Coulomb potential U (as indicated). Energies are expressed in units of ΓN.

Figure 9 shows evolution of the spectral function ρ(ω) for various couplings tm. In the weak-coupling limit we clearly observe a reduction (by half) of the initial density of states. With increasing tm the spectrum develops the three-peak structure that is typical for the “molecular” limit. This behavior indicates that the Majorana and Kondo states have rather a complicated relation, which is neither competitive nor cooperative. In fact, some novel scaling laws have been recently reported by several authors [69,70,76-79] also considering the correlation effects directly in the Rashba nanowire.


Figure 9: The spectral function ρ(ω) of the half-filled quantum dot (ε = −U/2) obtained at T = 0 for ΓSN = 2, UN = 4 and several values of tm (as indicated).


We have studied the polarized bound states of magnetic impurities embedded in an s-wave superconducting material, taking into account the spin–orbit and/or Coulomb interactions. We have shown that spin–orbit coupling strongly affects the subgap states, both of the single impurities and their conglomerate arranged into a nanoscopic chain. For the case of single magnetic impurity the spin–orbit interaction (i) shifts the quantum phase transition towards higher magnetic coupling Jc, (ii) enhances the spatial size of the YSR states, and (iii) smoothes the particle–hole oscillations. For the magnetic chain spin–orbit coupling combined with the Zeeman term induce the topologically nontrivial superconducting state and indirectly give rise to substantial polarization of the Majorana modes (Figure 4), the oscillations of which show up near the chain edges (Figure 5). The polarized Majorana quasiparticles can also leak into other side-coupled objects, such as single or multiple quantum impurities (Figure 6). These polarized Majorana quasiparticles can be controlled by a magnetic field or by an electrostatic potential. This would be important for future quantum computers using qubits based on topologically protected Majorana states. Finally, we have also confronted the Majorana quasiparticles with the subgap Kondo effect, revealing their complex relationship that can be hardly regarded as competitive or collaborative in some analogy to the Kondo effect originating from multiple degrees of freedom [80]. The aforementioned spin-polarized effects can be experimentally verified by polarized ballistic tunneling or by using STM spectroscopy, relying on the selective equal-spin Andreev reflections.


We thank for instructive remarks from R. Aguado, J. Klinovaja, R. Lutchyn, P. Simon, and R. Žitko on different parts of our study. This work was supported by the National Science Centre (Poland) under grants DEC-2014/13/B/ST3/04451 (AK, SG, TD) and DEC-2013/11/B/ST3/00824 (MMM) UMO-2017/25/B/ST3/02586 (AP) and by the Faculty of Mathematics and Natural Sciences of the University of Rzeszów through the project WMP/GD-06/2017 (GG).


  1. Choi, S.; Choi, H. J.; Ok, J. M.; Lee, Y.; Jang, W.-J.; Lee, A. T.; Kuk, Y.; Lee, S.; Heinrich, A. J.; Cheong, S.-W.; Bang, Y.; Johnston, S.; Kim, J. S.; Lee, J. Phys. Rev. Lett. 2017, 119, 227001. doi:10.1103/PhysRevLett.119.227001
    Return to citation in text: [1]
  2. Kenzelmann, M.; Strässle, T.; Niedermayer, C.; Sigrist, M.; Padmanabhan, B.; Zolliker, M.; Bianchi, A. D.; Movshovich, R.; Bauer, E. D.; Sarrao, J. L.; Thompson, J. D. Science 2008, 321, 1652. doi:10.1126/science.1161818
    Return to citation in text: [1]
  3. Meul, H. W.; Rossel, C.; Decroux, M.; Fischer, Ø.; Remenyi, G.; Briggs, A. Phys. Rev. Lett. 1984, 53, 497. doi:10.1103/PhysRevLett.53.497
    Return to citation in text: [1]
  4. Balatsky, A. V.; Vekhter, I.; Zhu, J.-X. Rev. Mod. Phys. 2006, 78, 373. doi:10.1103/RevModPhys.78.373
    Return to citation in text: [1] [2]
  5. Heinrich, B. W.; Pascual, J. I.; Franke, K. J. Prog. Surf. Sci. 2018, 93, 1. doi:10.1016/j.progsurf.2018.01.001
    Return to citation in text: [1]
  6. Yazdani, A.; Jones, B. A.; Lutz, C. P.; Crommie, M. F.; Eigler, D. M. Science 1997, 275, 1767. doi:10.1126/science.275.5307.1767
    Return to citation in text: [1]
  7. Ji, S.-H.; Zhang, T.; Fu, Y.-S.; Chen, X.; Ma, X.-C.; Li, J.; Duan, W.-H.; Jia, J.-F.; Xue, Q.-K. Phys. Rev. Lett. 2008, 100, 226801. doi:10.1103/PhysRevLett.100.226801
    Return to citation in text: [1]
  8. Franke, K. J.; Schulze, G.; Pascual, J. I. Science 2011, 332, 940. doi:10.1126/science.1202204
    Return to citation in text: [1] [2]
  9. Ruby, M.; Pientka, F.; Peng, Y.; von Oppen, F.; Heinrich, B. W.; Franke, K. J. Phys. Rev. Lett. 2015, 115, 087001. doi:10.1103/PhysRevLett.115.087001
    Return to citation in text: [1]
  10. Hatter, N.; Heinrich, B. W.; Ruby, M.; Pascual, J. I.; Franke, K. J. Nat. Commun. 2015, 6, 8988. doi:10.1038/ncomms9988
    Return to citation in text: [1]
  11. Ménard, G. C.; Guissart, S.; Brun, C.; Pons, S.; Stolyarov, V. S.; Debontridder, F.; Leclerc, M. V.; Janod, E.; Cario, L.; Roditchev, D.; Simon, P.; Cren, T. Nat. Phys. 2015, 11, 1013. doi:10.1038/nphys3508
    Return to citation in text: [1] [2] [3]
  12. Jellinggaard, A.; Grove-Rasmussen, K.; Madsen, M. H.; Nygård, J. Phys. Rev. B 2016, 94, 064520. doi:10.1103/PhysRevB.94.064520
    Return to citation in text: [1] [2]
  13. Choi, D.-J.; Rubio-Verdú, C.; de Bruijckere, J.; Ugeda, M. M.; Lorente, N.; Pascual, J. I. Nat. Commun. 2017, 8, 15175. doi:10.1038/ncomms15175
    Return to citation in text: [1]
  14. Assouline, A.; Feuillet-Palma, C.; Zimmers, A.; Aubin, H.; Aprili, M.; Harmand, J.-C. Phys. Rev. Lett. 2017, 119, 097701. doi:10.1103/PhysRevLett.119.097701
    Return to citation in text: [1]
  15. Salkola, M. I.; Balatsky, A. V.; Schrieffer, J. R. Phys. Rev. B 1997, 55, 12648. doi:10.1103/PhysRevB.55.12648
    Return to citation in text: [1] [2]
  16. Flatté, M. E.; Byers, J. M. Phys. Rev. Lett. 1997, 78, 3761. doi:10.1103/PhysRevLett.78.3761
    Return to citation in text: [1]
  17. Ugeda, M. M.; Bradley, A. J.; Zhang, Y.; Onishi, S.; Chen, Y.; Ruan, W.; Ojeda-Aristizabal, C.; Ryu, H.; Edmonds, M. T.; Tsai, H.-Z.; Riss, A.; Mo, S.-K.; Lee, D.; Zettl, A.; Hussain, Z.; Shen, Z.-X.; Crommie, M. F. Nat. Phys. 2016, 12, 92. doi:10.1038/nphys3527
    Return to citation in text: [1]
  18. Kezilebieke, S.; Dvorak, M.; Ojanen, T.; Liljeroth, P. Nano Lett. 2018, 18, 2311. doi:10.1021/acs.nanolett.7b05050
    Return to citation in text: [1]
  19. Santos, H.; Soriano, D.; Palacios, J. J. Phys. Rev. B 2014, 89, 195416. doi:10.1103/PhysRevB.89.195416
    Return to citation in text: [1]
  20. Kim, Y.; Zhang, J.; Rossi, E.; Lutchyn, R. M. Phys. Rev. Lett. 2015, 114, 236804. doi:10.1103/PhysRevLett.114.236804
    Return to citation in text: [1]
  21. Kaladzhyan, V.; Bena, C.; Simon, P. Phys. Rev. B 2016, 93, 214514. doi:10.1103/PhysRevB.93.214514
    Return to citation in text: [1] [2] [3]
  22. Ptok, A.; Głodzik, S.; Domański, T. Phys. Rev. B 2017, 96, 184425. doi:10.1103/PhysRevB.96.184425
    Return to citation in text: [1] [2]
  23. Žitko, R.; Lim, J. S.; López, R.; Aguado, R. Phys. Rev. B 2015, 91, 045441. doi:10.1103/PhysRevB.91.045441
    Return to citation in text: [1] [2] [3]
  24. Domański, T.; Weymann, I.; Barańska, M.; Górski, G. Sci. Rep. 2016, 6, 23336. doi:10.1038/srep23336
    Return to citation in text: [1] [2] [3] [4] [5]
  25. Smith, E. D. B.; Tanaka, K.; Nagai, Y. Phys. Rev. B 2016, 94, 064515. doi:10.1103/PhysRevB.94.064515
    Return to citation in text: [1]
  26. Goertzen, S. L.; Tanaka, K.; Nagai, Y. Phys. Rev. B 2017, 95, 064509. doi:10.1103/PhysRevB.95.064509
    Return to citation in text: [1]
  27. Koerting, V.; Andersen, B. M.; Flensberg, K.; Paaske, J. Phys. Rev. B 2010, 82, 245108. doi:10.1103/PhysRevB.82.245108
    Return to citation in text: [1]
  28. Pershoguba, S. S.; Björnson, K.; Black-Schaffer, A. M.; Balatsky, A. V. Phys. Rev. Lett. 2015, 115, 116602. doi:10.1103/PhysRevLett.115.116602
    Return to citation in text: [1]
  29. Głodzik, S.; Ptok, A. J. Supercond. Novel Magn. 2018, 31, 647. doi:10.1007/s10948-017-4360-6
    Return to citation in text: [1]
  30. Mashkoori, M.; Björnson, K.; Black-Schaffer, A. M. Sci. Rep. 2017, 7, 44107. doi:10.1038/srep44107
    Return to citation in text: [1]
  31. Satori, K.; Shiba, H.; Sakai, O.; Shimizu, Y. J. Phys. Soc. Jpn. 1992, 61, 3239–3254. doi:10.1143/JPSJ.61.3239
    Return to citation in text: [1]
  32. Sakai, O.; Shimizu, Y.; Shiba, H.; Satori, K. J. Phys. Soc. Jpn. 1993, 62, 3181–3197. doi:10.1143/JPSJ.62.3181
    Return to citation in text: [1]
  33. Sakurai, A. Prog. Theor. Phys. 1970, 44, 1472. doi:10.1143/PTP.44.1472
    Return to citation in text: [1]
  34. van Gerven Oei, W.-V.; Tanasković, D.; Žitko, R. Phys. Rev. B 2017, 95, 085115. doi:10.1103/PhysRevB.95.085115
    Return to citation in text: [1]
  35. Morr, D. K.; Yoon, J. Phys. Rev. B 2006, 73, 224511. doi:10.1103/PhysRevB.73.224511
    Return to citation in text: [1]
  36. Röntynen, J.; Ojanen, T. Phys. Rev. Lett. 2015, 114, 236803. doi:10.1103/PhysRevLett.114.236803
    Return to citation in text: [1]
  37. Morr, D. K.; Stavropoulos, N. A. Phys. Rev. B 2003, 67, 020502. doi:10.1103/PhysRevB.67.020502
    Return to citation in text: [1]
  38. Kawakami, T.; Hu, X. Phys. Rev. Lett. 2015, 115, 177001. doi:10.1103/PhysRevLett.115.177001
    Return to citation in text: [1]
  39. Nadj-Perge, S.; Drozdov, I. K.; Li, J.; Chen, H.; Jeon, S.; Seo, J.; MacDonald, A. H.; Bernevig, B. A.; Yazdani, A. Science 2014, 346, 602. doi:10.1126/science.1259327
    Return to citation in text: [1] [2] [3]
  40. Pawlak, R.; Kisiel, M.; Klinovaja, J.; Maier, T.; Kawai, S.; Glatzel, T.; Loss, D.; Meyer, E. npj Quantum Inf. 2016, 2, 16035. doi:10.1038/npjqi.2016.35
    Return to citation in text: [1] [2] [3]
  41. Ruby, M.; Pientka, F.; Peng, Y.; von Oppen, F.; Heinrich, B. W.; Franke, K. J. Phys. Rev. Lett. 2015, 115, 197204. doi:10.1103/PhysRevLett.115.197204
    Return to citation in text: [1] [2]
  42. Jeon, S.; Xie, Y.; Li, J.; Wang, Z.; Bernevig, B. A.; Yazdani, A. Science 2017, 358, 772. doi:10.1126/science.aan3670
    Return to citation in text: [1] [2] [3] [4] [5]
  43. Maśka, M. M.; Gorczyca-Goraj, A.; Tworzydło, J.; Domański, T. Phys. Rev. B 2017, 95, 045429. doi:10.1103/PhysRevB.95.045429
    Return to citation in text: [1]
  44. Liu, X.; Li, X.; Deng, D.-L.; Liu, X.-J.; Das Sarma, S. Phys. Rev. B 2016, 94, 014511. doi:10.1103/PhysRevB.94.014511
    Return to citation in text: [1]
  45. Stanescu, T. D.; Tewari, S. J. Phys.: Condens. Matter 2013, 25, 233201. doi:10.1088/0953-8984/25/23/233201
    Return to citation in text: [1]
  46. Sticlet, D.; Bena, C.; Simon, P. Phys. Rev. Lett. 2012, 108, 096802. doi:10.1103/PhysRevLett.108.096802
    Return to citation in text: [1]
  47. Maśka, M. M.; Domański, T. Sci. Rep. 2017, 7, 16193. doi:10.1038/s41598-017-16323-3
    Return to citation in text: [1] [2]
  48. Vernek, E.; Penteado, P. H.; Seridonio, A. C.; Egues, J. C. Phys. Rev. B 2014, 89, 165314. doi:10.1103/PhysRevB.89.165314
    Return to citation in text: [1]
  49. Deng, M. T.; Vaitiekenas, S.; Hansen, E. B.; Danon, J.; Leijnse, M.; Flensberg, K.; Nygård, J.; Krogstrup, P.; Marcus, C. M. Science 2016, 354, 1557. doi:10.1126/science.aaf3961
    Return to citation in text: [1]
  50. Liu, C.-X.; Sau, J. D.; Stanescu, T. D.; Das Sarma, S. Phys. Rev. B 2017, 96, 075161. doi:10.1103/PhysRevB.96.075161
    Return to citation in text: [1] [2]
  51. Hoffman, S.; Chevallier, D.; Loss, D.; Klinovaja, J. Phys. Rev. B 2017, 96, 045440. doi:10.1103/PhysRevB.96.045440
    Return to citation in text: [1]
  52. Ptok, A.; Kobiałka, A.; Domański, T. Phys. Rev. B 2017, 96, 195430. doi:10.1103/PhysRevB.96.195430
    Return to citation in text: [1] [2]
  53. Chevallier, D.; Szumniak, P.; Hoffman, S.; Loss, D.; Klinovaja, J. Phys. Rev. B 2018, 97, 045404. doi:10.1103/PhysRevB.97.045404
    Return to citation in text: [1]
  54. Chirla, R.; Moca, C. P. Phys. Rev. B 2016, 94, 045405. doi:10.1103/PhysRevB.94.045405
    Return to citation in text: [1]
  55. Prada, E.; Aguado, R.; San-Jose, P. Phys. Rev. B 2017, 96, 085418. doi:10.1103/PhysRevB.96.085418
    Return to citation in text: [1]
  56. Barański, J.; Kobiałka, A.; Domański, T. J. Phys.: Condens. Matter 2017, 29, 075603. doi:10.1088/1361-648X/aa5214
    Return to citation in text: [1]
  57. Domański, T.; Žonda, M.; Pokorný, V.; Górski, G.; Janiš, V.; Novotný, T. Phys. Rev. B 2017, 95, 045104. doi:10.1103/PhysRevB.95.045104
    Return to citation in text: [1]
  58. Mourik, V.; Zuo, K.; Frolov, S. M.; Plissard, S. R.; Bakkers, E. P. A. M.; Kouwenhoven, L. P. Science 2012, 336, 1003. doi:10.1126/science.1222360
    Return to citation in text: [1]
  59. Gül, Ö.; Zhang, H.; Bommer, J. D. S.; de Moor, M. W. A.; Car, D.; Plissard, S. R.; Bakkers, E. P. A. M.; Geresdi, A.; Watanabe, K.; Taniguchi, T.; Kouwenhoven, L. P. Nat. Nanotechnol. 2018, 13, 192. doi:10.1038/s41565-017-0032-8
    Return to citation in text: [1]
  60. Golub, A. Phys. Rev. B 2015, 91, 205105. doi:10.1103/PhysRevB.91.205105
    Return to citation in text: [1]
  61. Bauer, J.; Oguri, A.; Hewson, A. C. J. Phys.: Condens. Matter 2007, 19, 486211. doi:10.1088/0953-8984/19/48/486211
    Return to citation in text: [1]
  62. Yamada, Y.; Tanaka, Y.; Kawakami, N. Phys. Rev. B 2011, 84, 075484. doi:10.1103/PhysRevB.84.075484
    Return to citation in text: [1] [2] [3]
  63. Martín-Rodero, A.; Levy Yeyati, A. Adv. Phys. 2011, 60, 899. doi:10.1080/00018732.2011.624266
    Return to citation in text: [1]
  64. Barański, J.; Domański, T. J. Phys.: Condens. Matter 2013, 25, 435305. doi:10.1088/0953-8984/25/43/435305
    Return to citation in text: [1] [2]
  65. Liu, D. E.; Cheng, M.; Lutchyn, R. M. Phys. Rev. B 2015, 91, 081405. doi:10.1103/PhysRevB.91.081405
    Return to citation in text: [1]
  66. Elliott, S. R.; Franz, M. Rev. Mod. Phys. 2015, 87, 137. doi:10.1103/RevModPhys.87.137
    Return to citation in text: [1]
  67. Ruiz-Tijerina, D. A.; Vernek, E.; Dias da Silva, L. G. G. V.; Egues, J. C. Phys. Rev. B 2015, 91, 115435. doi:10.1103/PhysRevB.91.115435
    Return to citation in text: [1] [2]
  68. Lee, M.; Lim, J. S.; López, R. Phys. Rev. B 2013, 87, 241402. doi:10.1103/PhysRevB.87.241402
    Return to citation in text: [1] [2]
  69. Cheng, M.; Becker, M.; Bauer, B.; Lutchyn, R. M. Phys. Rev. X 2014, 4, 031051. doi:10.1103/PhysRevX.4.031051
    Return to citation in text: [1] [2]
  70. van Beek, I. J.; Braunecker, B. Phys. Rev. B 2016, 94, 115416. doi:10.1103/PhysRevB.94.115416
    Return to citation in text: [1] [2]
  71. Weymann, I. J. Phys.: Condens. Matter 2017, 29, 095301. doi:10.1088/1361-648X/aa5526
    Return to citation in text: [1] [2]
  72. Weymann, I.; Wójcik, K. P. Phys. Rev. B 2017, 95, 155427. doi:10.1103/PhysRevB.95.155427
    Return to citation in text: [1] [2]
  73. López, R.; Lee, M.; Serra, L.; Lim, J. S. Phys. Rev. B 2014, 89, 205418. doi:10.1103/PhysRevB.89.205418
    Return to citation in text: [1]
  74. Tanaka, Y.; Kawakami, N.; Oguri, A. J. Phys. Soc. Jpn. 2007, 76, 074701. doi:10.1143/JPSJ.76.074701
    Return to citation in text: [1]
  75. Vecino, E.; Martín-Rodero, A.; Levy Yeyati, A. Phys. Rev. B 2003, 68, 035105. doi:10.1103/PhysRevB.68.035105
    Return to citation in text: [1]
  76. Béri, B.; Cooper, N. R. Phys. Rev. Lett. 2012, 109, 156803. doi:10.1103/PhysRevLett.109.156803
    Return to citation in text: [1]
  77. Galpin, M. R.; Mitchell, A. K.; Temaismithi, J.; Logan, D. E.; Béri, B.; Cooper, N. R. Phys. Rev. B 2014, 89, 045143. doi:10.1103/PhysRevB.89.045143
    Return to citation in text: [1]
  78. Plugge, S.; Zazunov, A.; Eriksson, E.; Tsvelik, A. M.; Egger, R. Phys. Rev. B 2016, 93, 104524. doi:10.1103/PhysRevB.93.104524
    Return to citation in text: [1]
  79. Béri, B. Phys. Rev. Lett. 2017, 119, 027701. doi:10.1103/PhysRevLett.119.027701
    Return to citation in text: [1]
  80. Jacob, D.; Soriano, M.; Palacios, J. J. Phys. Rev. B 2013, 88, 134417. doi:10.1103/PhysRevB.88.134417
    Return to citation in text: [1]

© 2018 Głodzik et al.; licensee Beilstein-Institut.
This is an Open Access article under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The license is subject to the Beilstein Journal of Nanotechnology terms and conditions: (

Back to Article List