Manipulation of Majorana bound states in proximity to a quantum ring with Rashba coupling

The quest for Majorana zero modes in the laboratory is an active field of research in condensed matter physics. In this regard, there have been many theoretical proposals; however, their experimental detection remains elusive. In this article, we present a realistic setting by considering a quantum ring with Rashba spin-orbit coupling and threaded by a magnetic flux, in contact with a topological superconducting nanowire. We focus on spin-polarized persistent currents to assess the existence of Majorana zero modes. We find that the Rashba spin-orbit coupling allows for tuning the position of the zero energy crossings in the flux parameter space and has sizable effects on spin-polarized persistent currents. We believe that our results will contribute towards probing the existence of Majorana zero modes.

Majorana fermions remain a theoretical construct in the realm of particle physics. These particles, whose main feature is to be their own antiparticles, have been sought after in neutrinos, although without success thus far 1,2 . In condensed matter physics, Majorana quasiparticles emerge in p-wave topological superconductors 3,4 . Of particular relevance are the so-called Majorana zero modes (MZM) 2 , which are zero-energy Majorana quasiparticles protected by particle-hole symmetry. These quasiparticles, in contrast to their particle counterparts, exhibit yet another crucial feature. Such property is their being non-Abelian under braiding, which renders these quasiparticles as ideal candidates to perform topological quantum computing 5,6 . The reason for this stems in the large degeneracy found in the ground state of a system comprising a large number of Majorana quasiparticles. This degeneracy occurs due to these zero modes being pinned to zero energy 4 .
The archetypal model for studying Majorana modes in superconductors is that of  , where a onedimensional spinless p-wave superconductor is considered. Under suitable tuning of parameters, the model predicts that an ordinary superconductor undergoes a topological quantum phase transition into a topological superconductor. When considering open boundary conditions, the model predicts an odd number of MZMs at both ends of the chain with exponential decay into the bulk. Kitaev's model, although simple, poses significant challenges for its experimental realization. However, a great deal of experimental advances have taken place during the last few years, and a convenient setting has been put forward 4,10-12 . Such a setting has three key ingredients: a semiconductor nanowire with large spin-orbit coupling, such as InAs and InSb, a uniform magnetic field parallel to the wire to induce Zeeman splitting, and an s-wave superconductor such as Al. When placed in contact with the semiconductor and applying the magnetic field, the superconductor becomes topological when a given magnetic field is reached.
In this paper, we consider a superconducting nanowire in the topological regime placed in contact to a quantum ring. We assume the ring to be a continuation of the nanowire, as was done with a quantum dot in Ref. 13 , and therefore it displays strong spin-orbit coupling. The ring is in fact a collection of quantum dots assembled in a ring-like fashion, which can be achieved by means of gate voltages. The main results of this paper relate to the impact of the nanowire on the persistent spin currents of the ring, which are sharply different from those of the ring. As we shall show, when the MZMs interact via the quantum ring, they split from zero energy and oscillate, forming zero energy crossings. These crossings reshape the spin currents that would be present in the absence of Majoranas. Therefore, we believe our proposal could potentially be used as an additional tool to uncover the nature of MZMs.

System and model Hamiltonian
The system under consideration consists of a topological superconducting nanowire hosting MZMs at its edges, as depicted in Fig. 1a. The nanowire is placed in close proximity to a quantum ring, which is threaded by a magnetic flux. We introduce a minimal model that captures the influence of the MZMs on the persistent currents of the ring. The quantum ring consists of an array of quantum dots uniformly distributed on a circumference, as shown in Fig. 1b. Due to the similarities of our system with the one considered in Refs. 12,13 , we believe that our model might be amenable to be performed experimentally.
The model Hamiltonian is then given by where H r is the Hamiltonian of the ring, H M describes the interaction between the two MZMs at the edges of a nanowire in the topological superconducting regime, and H c contains the coupling between the MZMs and the quantum ring. Considering spin-orbit coupling and the fact that the quantum ring is threaded by a flux , H r can be written as [14][15][16] where c † n = c † n,↑ , c † n,↓ are the creation operators of the quantum ring and n runs from 1 to N, with N the number of sites of the ring, that is, the number of quantum dots that compose it. We shall measure energies in units of 2t with t the hopping energy, which implies that t = σ 0 /2 , with σ 0 the 2 × 2 identity matrix. The flux is accounted for in the phase factor θ = 2π �/N , where � = φ/φ 0 is the flux φ measured in units of the flux quantum φ 0 = h/e . The spin-orbit coupling is given by where we only consider Rashba spin-orbit interaction, disregarding a possible Dresselhaus interaction 15 . Here, σ x and σ y are the Pauli matrices, α = αN/2πR being R the radius of the ring, α the strength of the Rashba interaction and Regarding the mutual coupling between MZMs and the coupling between MZMs and the ring, we will consider a low energy theory as in Ref. 17 . Thus, the interaction term between the two MZMs is given by the minimal Kitaev Hamiltonian for the topological phase 7 where γ 1 and γ 2 are Majorana operators which satisfy γ i = γ † i and γ i , γ j = δ ij . Here ξ M ∝ e −L/ℓ 0 where L is the length of the nanowire and ℓ 0 is the superconducting coherence length 7 . The Majorana operators can also be written in terms of ordinary fermion operators, f and f † , as In terms of these operators, Eq. (5) takes the form (3) t so (n) = iα σ x cos ϕ n + σ y sin ϕ n , www.nature.com/scientificreports/ Finally, the coupling between the MZMs and the quantum ring is given by [17][18][19] where 1 and 2 are real parameters [see also Fig. 1]. We have chosen to couple the MZMs to a single spin. This can be achieved by fixing the spin canting angles of the MZMs accordingly, which can be done by means of sufficiently large magnetic fields 17 . Finally, we will not consider electron-electron interactions in the ring, which could be included within a mean field approximation 16 We can now find the spectrum of H by using the Bogoliubov-de Gennes (BdG) Hamiltonian. Indeed, if we introduce with d n,µ = d 1,− , d 1,+ , . . . , d N,− , d N,+ and likewise for d † n,µ , we can write H as The BdG Hamiltonian leads to a particle-hole redundancy and, as a result, the spectrum will be symmetrical around E = 0 . Diagonalization of H BdG provides us with the quasiparticle spectrum, E(�).

Quasiparticle energy spectra
In Figs. 2 and 3, we show the quasiparticle spectrum for rings with N = 3 and N = 4 . The spectra for larger values of N is provided in the Supplemental Material S1, and restrict to N = 3 and 4 in the remaining of the text. This is mainly for practical reasons, since the larger the number of dots, the harder an experimental implementation would prove to be.. In the following, we consider two scenarios, namely 1 = 0.1, 2 = 0 and 1 = 0.1, 2 = 0.15 , respectively. As explained in Ref. 13 , the quantity η = 2 / 1 provides a measure of locality and, as reported in that reference, a value of η ≪ 1 corresponds to highly nonlocal Majoranas, whereas the Majorana components of an Andreev bound state have η ≃ 1 . Since in our case 1 and 2 are interchangeable, we can see that the values chosen are in the nonlocal regime, that is, η < 1 . In Figs. 2a,b,e,f, we consider the case of α = 0 and coupling to a single Majorana mode, namely to γ 1 . In this case, we can observe that the spectrum is symmetric as E(�) = E(−�) and the Majorana modes remain at zero energy [see panels (a) and (b)] as long as there is no  we can see that the band that shifts to the left-hand side displays both gapless and gaped Dirac-like spectra close to zero energy, while the band that moves rightwards only shows a Dirac cone-like feature. However, for N = 3 we can see that the two shifted sides present a gaped Dirac-like spectra; this is because the Rashba spin-orbit coupling makes possible the interaction between MZMs and spin down. This shift is also observed in the oscillatory behavior of the MZMs as they are coupled together. However, these shifts follow the Dirac points corresponding to the spin-up states to which γ 1 is coupled, as we can observe in panels (g) and (h). In Fig. 3a,b,e,f, we consider the effect of coupling between the two MZMs while keeping the Rashba spin-orbit coupling turned off, α = 0 . As we can observe, the degeneracy is further lifted by coupling the ring to the two MZMs, as can be observed at zero flux. However, as we expressed above, it is not entirely lifted since there are still gapless Dirac-like features at low energy. It is interesting to notice that there is a clear distinction between even-and odd-numbered rings. Indeed, even when the MZMs are explicitly uncoupled in the Hamiltonian, when ξ M = 0 , there is hybridization due to their coupling through the ring. Although we do not show it here, this behavior can be seen to occur in rings with a larger number of sites. This effect may be due to the fact that even-numbered rings with α = 0 and at zero flux present particle-hole symmetry c n,σ → (−1) n c † n,−σ and c † n,σ → (−1) n c n,−σ , whereas odd-numbered rings do not, making these two systems inherently different. Notice that this particle-hole symmetry is a true symmetry, whereas the particle-hole symmetry of the BdG Hamiltonian is a redundancy. Finally, in panels (c), (d), (g) and (h), we show the effect of including Rashba spin-orbit coupling, which breaks the symmetry E(�) = E(−�) by shifting sideways the bands, as occurred in panels (c), (d), (g) and (h) of Fig. 2. However, in contrast to the latter, the coupling to γ 2 leads to a further reduction of the degeneracies in the problem. In fact, for N = 4 , the spectrum is wholly gaped out, even for ξ M = 0 . There is then a clear distinction with the corresponding spectrum without Rashba spin-orbit coupling, which will affect the persistent currents, as we shall show below.

Charge and spin persistent currents
Persistent currents have been studied extensively, both theoretically 16,20-25 and experimentally [26][27][28][29] . They constitute a peculiar quantum phenomenon resulting from the time-reversal symmetry breaking that occurs from coupling electrons to a gauge field, as in the Aharonov-Bohm effect.
In this section, we will study the appearance of persistent charge, I q , and spin currents, I s , in all the scenarios discussed above. We will calculate I q as and I s as where N σ = n c † n,σ c n,σ . Notice that the charge persistent current can also be obtained from the eigenenergies E ν by summing over all states I q = ν I q,ν where 16 and f (E ν ) the Fermi distribution function. The total persistent current is obtained by summing over all states 16 . After lengthy but straightforward calculations, we find that As we can see, the expressions for the persistent currents are identical, except for the coefficients. For the charge persistent current, we find where The expressions for the spin persistent current coefficients are identical by replacing sin ↔ cos in A ± n and B n and by taking α → −α in B n . Considering the fluctuation-dissipation theorem, we can write Here g h,± = (z + ε n,± ) −1 are the Green's functions for holes of the ring uncoupled to the MZMs. In order to assess the effect of having MZMs, we first show in Fig. 4 the charge and spin persistent currents for the system without MZMs. Additionally, we include the quasiparticle spectra where it can be seen that, as expected, the only effect is a lateral shift, as indicated by solid and dashed lines. Here, panels (b) and (c) correspond to the charge persistent currents for the N = 3 and N = 4 cases respectively and panels (d) and (e) correspond to the spin persistent currents for the N = 3 and N = 4 cases respectively. The different colours indicate the value of the Rashba spin-orbit interaction. As it can be noticed, the Rashba interaction leads to a secondary peak in the charge persistent current and to nonzero spin persistent currents.
n β n,− 2 g h,− + β n,+ 2 g h,+ + β n,− 2 g p,− + β n,+ 2 g p,+ .    spectra] and another one in between, ξ M = 0.05 , whose spectrum is not shown. Figure 5a,b,e,f corresponds to the settings of Fig. 2a,b,e,f. The first thing we observe is that the persistent currents display the periodicity observed in the spectra, as it should. In the case of the charge currents, we can see that upon decreasing the interaction between the two MZMs, the peaks in the current become smoother. On the other hand, the spin currents are nonzero due to the coupling of one of the two MZMs to the up-spin electrons in the ring. However, as the coupling between MZMs increases, reducing their topological robustness, the spin current tends towards zero. This is consistent with the results shown in the Fig. 4 as a greater interaction hybridizes the Majoranas. Figure 5c,d,g,h refers to the scenario considered in Fig. 2c,d,g,h, which is identical to the previous one except for a nonzero Rashba spin-orbit coupling. As we can see, in this case, the charge current is affected, as signaled by the appearance of additional peaks. Since the position of these peaks depends on the value of α , we can fine-tune the persistent current by manipulating this parameter. However, the sharpest distinction occurs when considering the spin currents. As we can observe, Rashba renders the spin currents asymmetric around zero current and, additionally, it leads to a pulse-like pattern, similar to the one we can see in Fig. 4.
In Fig. 6a,b,e,f we show the persistent currents for the parameters of Fig. 3a,b,e,f. In this case, the Rashba spin-orbit coupling is turned off, but the two MZMs are coupled to the ring. As we can observe, the results for the charge current are not dramatically altered, except for the curves being smoother in the N = 4 case and for a secondary peak at N = 3 as ξ M is increased to ξ M = 0.1 . These peaks feature in a clearer fashion in the spin currents, where a small plateau for ξ M = 0.1 is observed at certain values of flux. This can be understood from Fig. 3. Take for instance the spectrum around = 0 in panel (f) of that figure. Although it is very subtle, there is a switching of the Majorana states at zero energy and it is in the region between the two crossings where this plateau takes place. Additionally, we can observe that there is no suppression of the spin current upon increasing ξ M , as there was when 2 = 0 . Therefore, we can see that, by coupling the ring to the two MZMs we are able to suppress spin current for a range of fluxes. This could be a smoking gun in the detection of MZMs as this situation is completely different from the result in the absence of MZMs (cf. Fig. 4).
Finally, Fig. 6c,d,g,h corresponds to the case shown in Fig. 3c,d,g,h. In this scenario, the currents are similar to those of panels (a), (b), (e) and (f), but with the additional peak that was observed for N = 4 in Fig. 5c, although clearly smoothed out. In fact, the additional peak for the N = 3 case is completely suppressed. This has an effect on the spin current, where the N = 4 case still shows a pulse-like behaviour, but the N = 3 case shows plateaus for a larger range of fluxes than in panels (a), (b), (e) and (f).

Conclusions
MZMs stand out for their remarkable properties in the realm of topological quantum matter, such as their non-Abelian statistics. However, although significant progress has been made in the direction of their experimental detection in recent years, they cannot be confirmed as yet. More importantly, recent discoveries show that it is necessary to carefully examine different effects to determine their existence 31 unambiguously. In this regard, we propose to consider a topological superconducting nanowire in proximity to a quantum ring threaded by a magnetic flux. It is the natural extension of the system presented in Ref. 17 where a nanowire is coupled to a single quantum dot. In addition, we propose the ring to display Rashba spin-orbit coupling. As we have shown in this paper, the BdG spectra are entirely altered by coupling the quantum ring to one or two MZMs and having a nonzero Rashba spin-orbit coupling. The consequences are clearly observed in the persistent currents. Interestingly, spin-orbit coupling leads to an asymmetric spin current with a square pulse-like behavior. Moreover, the coupling to the two MZMs leads to plateaus for odd-numbered rings due to the switching of MZMs at zero energy. We believe that our work could serve as an additional signature towards probing the existence of Majorana zero modes.