Self-stabilizing temperature driven crossover between topological and non-topological ordered phases in one-dimensional conductors

Self-stabilizing temperature driven crossover between topological and non-topological ordered phases in one-dimensional conductors

Bernd Braunecker SUPA, School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews KY16 9SS, UK    Pascal Simon Laboratoire de Physique des Solides, CNRS UMR-8502, Université de Paris Sud, 91405 Orsay Cedex, France
September 26, 2019

We present a self-consistent analysis of the topological superconductivity arising from the interaction between self-ordered localized magnetic moments and electrons in one-dimensional conductors in contact with a superconductor. We show that due to a gain in entropy there exists a magnetically ordered yet non-topological phase at finite temperatures that is relevant for systems of magnetic adatom chains on a superconductor. Spin-orbit interaction is taken into account, and we show that it causes a modification of the magnetic order yet without affecting the topological properties.

74.20.Mn, 71.10.Pm, 75.30.Hx, 75.75.-c

Introduction. Topological superconductors have received much attention recently, partly because they host exotic low energy excitations such as Majorana bound states kane_RMP ; alicea ; leijnse , whose non-Abelian statistics are attractive for topological quantum computation nayak ; pachos:2012 . As a remarkable feature, topological superconductivity can be created artificially by contacting specific materials with a conventional -wave superconductor. For instance, it arises at the interface between the surface states of a three-dimensional topological insulator and a -wave superconductor Fu-Kane ; in one-dimensional (1D) semiconducting wires with a strong spin-orbit interaction (SOI) and a Zeeman magnetic field with proximitized superconductivity lutchyn:2010 ; oreg:2010 ; mourik:2012 ; das:2012 ; deng:2012 ; or in arrays of magnetic nanoparticles or magnetic adatoms on top of a superconducting surface choy:2011 ; nadj-perge:2013 ; pientka:2013 ; pientka:2014 ; poyhonen:2014 ; rontynen:2014 ; kim:2014 ; heimes:2014 ; brydon:2015 ; ojanen:2015 ; schecter:2015b such as iron adatoms on lead nadj-perge:2014 ; pawlak:2015 ; ruby:2015 .

Figure 1: Zoom on a 1D conductor with embedded magnetic moments on top of a superconductor. The magnetic moments self-order in the form a spiral order with spatial period .

The systems we consider in this letter exhibit a topological phase emerging from self-organization of magnetic moments embedded in 1D conductors with proximity induced superconductivity. This situation may apply to semiconducting wires with extrinsic magnetic impurities or intrinsic moments such as nuclear spins, or a conducting wire made of magnetic adatoms on a superconducting surface. Due to the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction mediated through the electrons, the magnetic moments can undergo an ordering transition below a temperature and form a spiral with a spatial period characterized by the wave number (see Fig. 1) such that , for the Fermi momentum. This ordering mechanism was first demonstrated for normal conductors braunecker:2009a ; braunecker:2009b , then conjectured nadj-perge:2013 and self-consistently demonstrated braunecker:2013 ; klinovaja:2013 ; vazifeh:2013 for the superconducting case. These results were corroborated recently by showing that the spiral order persists beyond the RKKY limit, and stays close to , as long as is away from commensurate band fillings and the coupling strength between magnetic moments and electrons remains smaller than the electron bandwidth schecter:2015 ; hu:2015 .

The locking of to has important consequences. The magnetic spiral forms a periodic superstructure that causes a part of the electrons to undergo a spin-selective Peierls transition braunecker:2010 to a non-conducting spiral electron spin density wave, whereas the remaining conducting electron states become helical (spin-filtered). The induced superconductivity then becomes of the topological -wave type, and Majorana bound states appear at the two ends of the 1D wire. A system with such a spiral magnetic order is indeed equivalent braunecker:2010 ; choy:2011 ; kjaergaard:2012 ; martin:2012 ; jelena:2012 ; nadj-perge:2013 ; braunecker:2013 ; klinovaja:2013 ; vazifeh:2013 ; pientka:2013 ; pientka:2014 ; poyhonen:2014 ; rontynen:2014 ; kim:2014 ; ojanen:2015 ; jelena:2015 to the original proposals for topological superconductivity in nanowires lutchyn:2010 ; oreg:2010 . Remarkably, by this mechanism, the topological superconducting phase emerges naturally as the ground state without any fine tuning. Although both the RKKY based and the nonperturbative approaches consistently predict the locking condition , it must be stressed that the former results are based on the further condition of large magnetic and superconducting gap energies, whereas the latter apply only at zero temperature.

In this letter, we provide a general analysis which incorporates entropy and thermal fluctuations in the non-zero temperature regime and we show that there exists a previously unknown crossover to a magnetically ordered yet non-topological phase. Furthermore, we show that the spin-orbit interaction (SOI), which is genuinely present in such systems either intrinsically or through interface effects, causes a modification of the magnetic spiral but has no effect on the topological properties.

Heuristic considerations. The physical origin of the non-topological phase is illustrated in Fig. 2.

Figure 2: Example of dispersion relations in the transformed basis obtained by the spin dependent shift of momenta for which the spiral effective magnetic field becomes ferromagnetic and the Zeeman like gap opens at braunecker:2010 ; note_offset . In the plots, the chemical potential remains constant, but varies. Top panels for the normal state (), with red and blue colors corresponding to opposite spin projections perpendicular to the spiral field. Bottom panels for the induced superconducting state (, with ). For , the chemical potential (dashed line) lies in the middle of the gap (b) and the superconducting system is fully gapped (e). For smaller or larger , the gap lies at lower or higher energies and eventually touches the upper (a) or lower gap edge (c). At both touching points the superconducting gap closes (d,f) and the state becomes non-topological.

At low temperatures , the thermodynamic ground state is determined by the gain in electronic energy obtained by maintaining large magnetic () and superconducting () gaps, and the system adjusts to . As is raised, however, the ground state is dictated by the free energy , and the entropy can play a decisive role. Indeed, if is lowered or raised to the values as indicated in the left and right panels of Fig. 2, such that the chemical potential touches the bottom or top of a band, the induced superconducting gap closes (for ) because the touched bands are fully spin polarized. The effective dispersion arising from the superconducting case becomes gapless, with a larger entropy than in the gapped case. As a result, if is large enough, typically (with the Boltzmann constant), the minimization of can be dominated by the enhancement of , and the thermodynamic ground state corresponds to situations (d) or (f) in Fig. 2, a topologically trivial yet magnetically ordered phase.

Not yet taken into account in this argument is the stability of this phase upon thermal fluctuations of the magnetic moments. As shown in braunecker:2009a ; braunecker:2009b ; braunecker:2013 , for both the ungapped and gapped cases a mean-field description of the spin-wave fluctuations captures the correct value of the ordering temperature and , with the transverse spin susceptibility at momentum . Since increases for a gapless dispersion, closing the superconducting gap by moving away from causes furthermore an enhanced stability against thermal fluctuations, provided that the effect occurs at . Obviously, the latter condition depends on the considered material as specified below. For practical implementations, the condition is a priori required for the topological self-tuning phase to be accessible at high enough temperatures, and for close to the non-topological ordered phase may indeed be favored. Yet we find that for semiconductor bands with an effective mass as in Fig. 2, the value of remains generally still too low, such that would dominate only at temperatures where no order can persist anyway.

This situation changes drastically for tight binding systems such as shown in Fig. 3, which are the natural description for adatom chains. Due to the cosine nature of the dispersion, magnetic gaps appear at two points in the Brillouin zone, and can self-adjust such that the superconducting-magnetic gaps at both points become equal and fulfill the condition [see Fig. 3 (f)]. At the effective doubling of thermally accessible states provides a doubling of the value of . This is sufficient to push the transition to the non-topological to in precisely the systems that are most attractive for realizing a self-sustained topological phase. Furthermore, the equality of the two gaps leads to (with the lattice spacing), which corresponds to an antiferromagnetic arrangement of the magnetic moments if the latter are on the same lattice sites.

Figure 3: Dispersion relations as in Fig. 2 for the tight binding model associated with Eq. (1) and corresponding to Fig. 4. Parts (a,d) represent the low-temperature topological phase; parts (b,e) the crossover to the non-topological phase where the gap closes; parts (c,f) the final high-temperature phase where and the system is magnetically ordered yet non-topological.

Quantitative analysis. For a quantitative investigation we consider a quantum wire with induced superconductivity and embedded magnetic moments, described by the Hamiltonian


Here are the operators for electron with spin and dispersion relation . is the chemical potential, the induced superconducting gap, the vector of Pauli matrices, and the vector the effective SOI in the system, arising from the sum of SOI contributions due to the internal structure of the wire or to interface effects with the substrate. The vectors are the Fourier transforms of the chain of magnetic scatterers coupling to the electron spin, where are the magnetic moments and is the coupling strength. The are placed at positions that can be irregular but are sufficiently dense with respect to such that they can be considered as a continuum.

For and , the ground state energy is minimized if the vectors are confined to an arbitrary two-dimensional plane, spanned by orthogonal unit vectors , in which they rotate as a spiral as a function of , choy:2011 ; kjaergaard:2012 ; martin:2012 ; nadj-perge:2013 ; braunecker:2013 ; klinovaja:2013 ; vazifeh:2013 ; kim:2014 ; schecter:2015 ; hu:2015 . Choosing such that , the corresponding term in the Hamiltonian becomes . Letting , braunecker:2010 produces a unitary transformation diagonalizing the Hamiltonian, in which forms a uniform ferromagnetic coupling along the spin- direction. For a parabolic dispersion with the band mass we obtain Fig. 2.

Since the SOI term is linear in , a produces a similar spin-dependent momentum shift. For a quadratic dispersion with the band mass, and spin axes such that , the SOI can be absorbed by letting , with braunecker:2010 . If is not perpendicular to the plane, the two shifts by and are not compatible, and the diagonalization of would generally mix all momenta. Such modification of the long-ranged wave functions by the spiral order would cause an extensive energy cost which is not favored energetically. This can be avoided, however, by an alignment of the spiral to the plane perpendicular to . We thus define the spin directions such that , and , for which and are parallel and can be directly added. Remarkably, to maintain the optimal by minimizing the free energy, the spiral undergoes an adjustment of to a such that . While the “” sign arises from the choice of spin axes, can have either sign, and so two spirals with opposite helicities and different periods, , are possible. Therefore even a large SOI has no further influence than the pinning of the plane of the magnetic spiral together with the adjustment of , provided that does not become smaller than the electron lattice spacing or the average spacing between the . As long as (up to dependent corrections that can be included), a measurement of the period and plane of the magnetic spiral could therefore give a direct measurement of .

Due to the extensive energy cost, there are furthermore no conical deformations out of the spiral plane kim:2014 ; reis:2014 . In a spin-Nambu matrix representation spanned by the vectors the Hamiltonian takes then the form


for , , , and the restriction of the summation to to avoid state overcounting. The energy offset is and due to its dependence must be kept for comparison of different . The diagonalization of the matrix in Eq. (2), for , leads to the energies , for , with for . This leads to the ground state energy


and the entropy


for the Fermi function. Notice that the sums are restricted to to avoid overcounting. To proceed, we exclude the well-known scattering effects occurring at commensurate filling factors that can cause orderings different from spirals schecter:2015 ; hu:2015 .

Minimizing the free energy determines the ordering vector . However, analyzing only is incomplete to assure the stability of the ordered phase since the long-wavelength spin-wave fluctuations smooth any magnetization at finite for a system of finite size braunecker:2009a ; braunecker:2009b ; braunecker:2013 ; klinovaja:2013 ; vazifeh:2013 . Taking this condition into account, it was demonstrated braunecker:2009b ; braunecker:2013 that for any realistic system size the mean field result,


provides the ordering temperature for both the gapped and gapless cases, where is the static transverse spin susceptibility, expressed in terms of the Fermi energy , the Fermi velocity , and the short distance cutoff (limited by lattice spacing or distance between the ). If is the gap as indicated in Figs. 2 and 3 the energy is roughly set by and its value reflects the transition between the gapped () and gapless regimes (). As a consequence, if departs from and shrinks, the value of initially grows but then saturates at a self-consistent value where . Electron interactions further modify the logarithm in Eq. (5) to a temperature dependent power-law in the gapless case and can lead to a considerable further increase of braunecker:2009a ; braunecker:2009b . We note that Eq. (5) results from analyzing fluctuations arising from the RKKY interaction. While the RKKY limit is insufficient to characterize the ground state and needs to be used, it correctly encodes the fluctuations away from the ground state configuration for .

Application to various systems. Since the transition to the non-topological phase occurs at , materials with large shall be considered. With typical meV, most systems with nuclear magnetic moments have , hence a guaranteed topological phase, albeit a low transition temperature . For most 1D wire situations (including InAs braunecker:2013 ), we indeed find that remains too low to allow for a strong impact of the entropy, and the self-sustained topological phases remain stable.

Dense chains of magnetic adatoms on a superconducting substrate have a larger coupling constant . If neighboring adatom orbitals hybridize, the chains become conductors of the type of Eq. (1). SOI effects are generally strong in such systems nadj-perge:2014 ; pawlak:2015 ; ruby:2015 . However, this causes here just a mere rearrangement of the magnetic helix . While systems such as in nadj-perge:2014 ; pawlak:2015 ; ruby:2015 likely depend much on the direct exchange interaction between neighboring moments, we focus here instead on the case where RKKY dominates over the latter (therefore our results do not a priori apply to nadj-perge:2014 ; pawlak:2015 ; ruby:2015 ).

As demonstrated in braunecker:2013 , with meV menzel:2012 , we obtain meV for the gapped phase. Requirement for such large is, as seen in Eq. (5) a large prefactor , which means a rather small bandwidth or a renormalization of the Fermi velocity Pientka:2015 . A tight binding model is therefore suitable, in which the factor is replaced by the hopping integral , and the dispersion relation is . Consequently, a minimum of at a shift has a particle-hole reversed minimum corresponding to the shift . Away from half filling, , both minima are inequivalent, yet for not too far from they can lie energetically close enough together such that further entropy can be gained by tuning through the topological boundary and pin it to , an antiferromagnetic order at which the system is non-topological and has small gaps but with both minima contributing to . In combination with a small bandwidth, a large , and , the entropy can become large enough to dominate . An example is given in Fig. 4, showing as a function of and . The values minimizing are indicated by the red line. At low temperatures lies near, but not on since and cause a significant band deformation due to the small band width. At we observe a crossing into the non-topological region (indicated by the hatching), well below , and the pinning to .

Figure 4: (a) Free energy for an adatom chain conductor of length 1 m as a function of total spiral momentum (including SOI contributions) and temperature note_offset . Contour lines complement the color coding. The minima are marked by the red line and correspond to the ground state configuration. The values , , and are marked by horizontal dashed lines. Parameters for the tight binding model are meV, Å, meV, meV, m, (7/20 filling). Stability of the order is ensured up to (at the right plot limits). The gap closes at the phase boundary between the topological (clear) and non-topological (hatched) regions. The system becomes non-topological at and then stabilizes at , corresponding to an antiferromagnetic order. (b) Values of minimizing the free energy [same as in (a)], the entropy , and the energy . While alone favors a topological phase, the entropy favors a gapless phase, but is at higher temperatures further enhanced by tuning , and eventually dominates the minimum of at . For the two identical gaps indicated in Fig. 3 (f) are .

Conclusions. We have analyzed a 1D conductor with spin-orbit interaction coupled to a 1D chain of magnetic moments. Through a self-consistent analysis taking into account the full electronic free energy and the fluctuations about the ordered magnetic ground states, we have determined the stability of the topological superconducting phase at finite temperature. We showed that spin-orbit interaction causes only a pinning of the plane of the magnetic spiral and an adjustment of its spatial period. Furthermore, in some situations especially met in systems of magnetic adatoms, we demonstrated that there is a significant temperature range, in which a magnetic order persist but the electronic state is non-topological.

Acknowledgments. We thank C. Carroll, J. Klinovaja, D. Loss and D. Morr for helpful discussions. PS acknowledges support by the French Agence Nationale de la Recherche through the contract ANR Mistral.


  • (1) M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010).
  • (2) J. Alicea, Rep. Prog. Phys. 75, 076501 (2012).
  • (3) M. Leijnse and K. Flensberg, Semicond. Sci. Technol. 27, 124003 (2012).
  • (4) C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das Sarma, Rev. Mod. Phys. 80, 1083 (2008).
  • (5) J. K. Pachos, Introduction to Topological Quantum Computation (Cambridge University Press, Cambridge, 2012).
  • (6) L. Fu and C. L. Kane, Phys. Rev. Lett. 100, 096407 (2008).
  • (7) R. M. Lutchyn, J. D. Sau, and S. Das Sarma, Phys. Rev. Lett. 105, 077001 (2010).
  • (8) Y. Oreg, G. Refael, and F. Oppen Phys. Rev. Lett. 105, 177002 (2010).
  • (9) V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. A. M. Bakkers, and L. P. Kouwenhoven, Science 336, 1003 (2012).
  • (10) A. Das, Y. Ronen, Y. Most, Y. Oreg, M. Heiblum, and H. Shtrikman, Nano Lett. 8, 887 (2012).
  • (11) M. T. Deng, C. L. Yu, G. Y. Huang, M. Larsson, P. Caroff, and H. Q. Xu, Nano Lett. 12, 6414 (2012).
  • (12) T.-P. Choy, J. M. Edge, A. R. Akhmerov, and C. W. J. Beenakker, Phys. Rev. B 84, 195442 (2011).
  • (13) S. Nadj-Perge, I. K. Drozdov, B. A. Bernevig, and A. Yazdani, Phys. Rev. B 88, 020407(R) (2013).
  • (14) F. Pientka, L. I. Glazman, and F. von Oppen, Phys. Rev. B 88, 155420 (2013).
  • (15) F. Pientka, L. I. Glazman, and F. von Oppen, Phys. Rev. B 89, 180505(R) (2014).
  • (16) K. Pöyhönen, A. Westström, J. Röntynen, and T. Ojanen, Phys. Rev. B 89, 115109 (2014).
  • (17) J. Röntynen and T. Ojanen, Phys. Rev. B 90, 180503 (2014).
  • (18) Y. Kim, M. Cheng, B. Bauer, R. M. Lutchyn, and S. Das Sarma, Phys. Rev. B 90, 060401(R) (2014).
  • (19) A Heimes, P Kotetes, G Schön, Phys. Rev. B 90, 060507(R) (2014).
  • (20) P. M. R. Brydon, S. Das Sarma, H.-Y. Hui, and J. D. Sau, Phys. Rev. B 91, 064505 (2015).
  • (21) K. Pöyhönen, A. Westström, T. Ojanen, arXiv:1509.05223.
  • (22) M. Schecter, K. Flensberg, M. H. Christensen, B. M. Andersen, J. Paaske, arXiv:1509.07399.
  • (23) S. Nadj-Perge, I. K. Drozdov, J. Li, H. Chen, S. Jeon, J. Seo, A. H. MacDonald, B. A. Bernevig, and A. Yazdani, Science 346, 602 (2014).
  • (24) R. Pawlak, M. Kisiel, J. Klinovaja, T. Meier, S. Kawai, T. Glatzel, D. Loss, and E. Meyer, arXiv:1505.06078.
  • (25) M. Ruby, F. Pientka, Y. Peng, F. von Oppen, B. W. Heinrich, and K. J. Franke, arXiv:1507.03104.
  • (26) B. Braunecker, P. Simon, and D. Loss, Phys. Rev. Lett. 102, 116403 (2009).
  • (27) B. Braunecker, P. Simon, and D. Loss, Phys. Rev. B 80, 165119 (2009).
  • (28) B. Braunecker and P. Simon, Phys. Rev. Lett. 111, 147202 (2013).
  • (29) J. Klinovaja, P. Stano, A. Yazdani, and D. Loss, Phys. Rev. Lett. 111, 186805 (2013).
  • (30) M. M. Vazifeh and M. Franz, Phys. Rev. Lett. 111, 206802 (2013).
  • (31) M. Schecter, M. S. Rudner, and K. Flensberg, Phys. Rev. Lett 114, 247205 (2015).
  • (32) W. Hu, R. T. Scalettar, and R. R. P. Singh, arXiv:1506.04809.
  • (33) B. Braunecker, G. I. Japaridze, J. Klinovaja, and D. Loss, Phys. Rev. B 82, 045127 (2010).
  • (34) M. Kjaergaard, K. Wölms, and K. Flensberg, Phys. Rev. B 85, 020503(R) (2012).
  • (35) I. Martin and A. Morpurgo, Phys. Revf. B 85, 144505 (2012).
  • (36) J. Klinovaja, P. Stano, and D. Loss, Phys. Rev. Lett. 109, 236801 (2012).
  • (37) C.-H. Hsu, P. Stano, J. Klinovaja, and D. Loss, arXiv:1509.01685.
  • (38) For a magnetic spiral, the shift diagonalizes the Hamiltonian without changing the energy, and hence the results for different are directly comparable. The same type of transformation also diagonalizes the SOI but adds an energy offset . However, since the SOI is fixed for a given system, this offset is the same for all further magnetic spirals and has no effect on the discussion.
  • (39) I. Reis, D. J. J. Marchand, and M. Franz, Phys. Rev. B 90, 085124 (2014).
  • (40) M. Menzel, Y. Mokrousov, R. Wieser, J. E. Bickel, E. Vedmedenko, S. Blügel, S. Heinze, K. von Bergmann, A. Kubetzka, and R. Wiesendanger, Phys. Rev. Lett. 108, 197204 (2012).
  • (41) Y. Peng, F. Pientka, L. I. Glazman, F. von Oppen, Phys. Rev. Lett. 115, 087001 (2015).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description