# Quantum phase diagram of Shiba impurities from bosonization

###### Abstract

A characteristic feature of Shiba impurities is the existence of a spin- and parity-changing quantum phase transition (known as “0-” transition) which has been observed in scanning tunneling microscopy (STM) and transport experiments. Using the Abelian bosonization technique, here we analyze the ground-state properties and the quantum phase diagram of a classical (i.e., Ising-like) Shiba impurity. In particular, we analyze the cases of an impurity in a three- and a one-dimensional superconductor. Within the bosonization framework, the ground-state properties are determined by simple soliton-like solutions of the classical equations of motion of the bosonic fields, whose topological charge is related to the spin and parity quantum numbers. Our results indicate that the quantum phase diagram of the superconductor can be strongly affected by geometrical and dimensional effects. Exploiting this fact, in the one-dimensional case we propose an experimental superconducting nanodevice in which a novel parity-preserving, spin-changing “0-0” transition is predicted.

###### pacs:

PACS ACAIntroduction.— Yu-Shiba-Rusinov (or simply Shiba) states are localized subgap states arising in a superconductor due to the local pair-breaking processes induced by a magnetic impurity Yu (1965); Shiba (1968); Rusinov (1968); Balatsky et al. (2006), or by a quantum dot in superconductor-hybrid nanodevices, where they are more commonly known as “Andreev bound states” De Franceschi et al. (2010); Deacon et al. (2010); Martín-Rodero and Yeyati (2011); Maurand et al. (2012). Recently, these systems have become the focus of intensive research due to the prediction that Majorana zero-modes (i.e., non-Abelian quasiparticles with potential uses in fault-tolerant topological quantum computation Nayak et al. (2008)) could be realized in a linear chain of magnetic impurities deposited on top of a superconductor. In these systems, the Shiba states can overlap and form a topological “Shiba band” Nadj-Perge et al. (2013); Klinovaja et al. (2013); Braunecker and Simon (2013); Pientka et al. (2013) that mimics the physics of the Kitaev one-dimensional model Kitaev (2001). Subsequent STM studies of chains of Fe atoms deposited on top of clean Pb surfaces have shown compelling evidence for the Majorana scenario Nadj-Perge et al. (2014); Pawlak et al. (2016); Ruby et al. (2015), generating a lot of excitement.

A salient feature of Shiba systems is the existence of an experimentally accessible quantum phase transition, known as the “0- transition”, determined by the position of the level inside the superconducting gap: if the Shiba state falls below the Fermi energy, it becomes occupied and the nature of the collective ground state changes from an even-parity BCS-like singlet to an odd-parity spin-1/2 doublet Sakurai (1970). This parity- and spin-changing transition has been experimentally observed both in adatom/superconductor systems with STM techniques Franke et al. (2011); Bauer et al. (2013), and in superconductor-quantum dot devices via quantum transport experiments Deacon et al. (2010); Lee et al. (2012); Maurand et al. (2012); Schindele et al. (2014); Island et al. (2017).

Recent progress in nanofabrication has allowed the realization of ultrathin superconducting epitaxial nanostructures [i.e., one-dimensional nanowires (1DNWs)] using the proximity effect Krogstrup et al. (2015); Taupin et al. (2016). These advances pave the way for novel technological applications, and can potentially enable the study of Shiba states in systems with reduced dimensionality. A question which naturally arises in this context is whether the different dimensionality or geometrical properties of the superconductor can qualitatively modify the properties of Shiba systems and the above-mentioned phase diagram. For instance, scattering processes which are crucial in the one-dimensional (1D) geometry can be profoundly affected by interactions Fisher and Kane (1992); Kane and Fisher (1992). In a 1D system (such as a proximitized 1DNW), the enhanced effect of correlations might indeed change the properties of a superconductor, thus modifying the features of the induced quantum phases.

In this work we implement the Abelian bosonization formalism to study the zero-temperature phase diagram of a Shiba impurity, both in a 1D and in a three dimensional (3D) superconductor. Interestingly, within the bosonic representation the ground-state properties in both geometries can be studied in an unified way. In particular for the 1D case, we propose a novel experimental device based on proximitized semiconducting 1DNWs with a ferromagnetic (FM) nanocontact grown ontop to induce a controllable Shiba state. We predict that this device could host a novel parity-conserving,spin-changing “0-0” transition. Our results might have impact in recent theoretical and experimental developments where Shiba states have been observed, and in recent works where superconductivity is induced in semiconducting nanowires by the means of proximity effect.

Theoretical model.— We begin by analyzing the 1D geometry. We describe a proximitized single-channel superconducting 1DNW of length with the Hamiltonian , where

(1) |

describes the (interacting) semiconductor wire. Here, anhilates a fermion with spin , and is a repulsive Hubbard-type interaction parameter (here we have used ). Linearization of the 1DNW spectrum in a region of width around the Fermi energy allows to write the Fermi fields as , where is the Fermi wavevector, and () are left/right moving chiral fields slowly-varying on the scale . We now introduce the bosonization formalism Giamarchi (2004); Gogolin et al. (1988), and represent the chiral fermion fields as where are chiral bosonic fields obeying the Kac-Moody commutation relations , is the short-distance cutoff of the continuum theory, and are standard Klein factors which ensure the proper anticommutation relations of the Fermi fields.

It is customary to introduce the new fields , where () refers to charge- (spin-) type operators, and where the convention on the r.h.s has been used for the branch, and similarly, for . The new fields satisfy canonical commutation relations , and physically, the field is related to the spin density operator , while is related to the Josephson phase field of the superconductor and to the current density . In terms of these fields, the Hamiltonian takes a Luttinger liquid form with decoupled charge and spin bosonic sectors, i.e., , where

(2) |

for . The parameter encodes the interactions in each sector, and physically controls the decay of the correlation function . In our particular case, , and due to the SU(2) symmetry of the model, the value of is constrained to . On the other hand, is the velocity of the 1D charge plasmons, and is the velocity of the 1D spinons.

Next, the proximity-induced superconducting pairing in the 1DNW is described by , where we have neglected the rapidly oscillating terms proportional to since they average to zero. Physically, the induced pairing potential emerges from the integration of the bulk superconductor, and strongly depends on the transparency and disorder of the superconductor/nanowire contact Takei et al. (2013). In terms of the bosonic fields, the pairing term writes . A simple scaling analysis where we change the cutoff indicates that is a relevant perturbation in the RG sense, i.e., , and flows to strong coupling as when .

Finally, Shiba states emerge due to the presence of a localized magnetic field or impurity in the superconductor. In an ultrathin mesoscopic 1DNW, atomic-sized impurities are hard to control and manipulate experimentally. Therefore, in order to induce a controllable Shiba state, we assume the presence of a FM insulating contact of width grown ontop of the 1DNW. This Hamiltonian writes , where is the exchange field induced by the FM contact, and is its dimensionless magnetization profile, assumed to be oriented along the -axis Sau et al. (2010). For concreteness, we assume a Gaussian profile , where the magnetization in the bulk can take values depending on the temperature, size, material details, etc. Under the reasonable assumption that , with the coherence length in the 1DNW, the FM contact becomes a “point-like” perturbation from the perspective of the chiral fields , i.e., . Therefore, the Hamiltonian writes , where we have neglected oscillatory terms ). In fermion language, this effective delta-like potential can be eliminated with a gauge transformation , where we have defined the phase shift , and where is the unit step function. The bosonized Hamiltonian of the FM contact then writes . Experimental values nm Jespersen et al. (2009) and nm Albrecht et al. (2016) are encouraging for the realization of this device. Note that in the absence of pairing the Hamiltonian is akin to the X-ray edge problem Gogolin et al. (1988); Giamarchi (2004).

An interesting feature of the bosonic representation is that it allows to obtain the total spin and parity of the ground state from simple expressions

(3) | ||||

(4) |

The physical significance of can be clarified imposing periodic boundary conditions in the problem, i.e., . In this case, can only take integer values which are related to the total number of particles through the relation . Then, since , we conclude that determines the fermion-parity of the ground state Haldane (1981).

Note that in the terms and , only the commuting fields and appear. This indicates that the Hamiltonian has a well-defined “classical” limit when flows to strong coupling, i.e., when it becomes of the order of the ultraviolet cutoff . This classical limit is more conventiently captured defining the fields , and their canonical conjugates , in terms of which we have

(5) | ||||

(6) | ||||

(7) |

where we have defined , , and the helicity parameters and . Even though the classical limit of can be studied in the generic case, here we concentrate on the simpler but more insightful point (i.e., ), where the classical static configurations of decouple. The generic case will be discussed elsewhere Bortolin et al. (). At the point , we obtain the decoupled equations

(8) |

Here , and is the derivative of the Dirac delta-function, which has to be interpreted in terms of the approximant as .

3D Case.—It is illuminating to briefly turn to the well-known case of a magnetic impurity in a 3D superconductor Yu (1965); Shiba (1968); Rusinov (1968). For a spherically-symmetric problem, only the -wave component of the 3D Fermi field couples to a point-like impurity placed at the origin. Then, linearizing the spectrum around the Fermi energy, we can write , where is the radial coordinate. Equivalently, the half-line can be unfolded to the whole line by defining

(9) |

and keeping a single chirality for each fermion. This definition includes the boundary condition at , which guarantees that is finite at the origin. This procedure is standard and we refer the reader to Refs. Giamarchi, 2004; Gogolin et al., 1988 for details. In what follows, when treating the 3D case, we will keep only the branches and . Imposing a hard-wall boundary condition in an sphere of radius (i.e., ), induces periodic boundary conditions on the effective 1D chiral fermions, . Therefore, we can write

(10) |

which upon bosonization results , where only the field appears. Interestingly, Eq. (10) is the same Hamiltonian that describes the edge states of a spin-Hall insulator Wu et al. (2006). In addition, since the (repulsive) interactions are efficiently screened in a 3D superconductor, we have and , and the system is automatically at the point . On the other hand, it is easy to see that due to the elimination of the chiral fields , the pairing Hamiltonian in the -wave channel writes , where here is an intrinsic interaction. Finally, the presence of a classical magnetic impurity is usually modelled using the - Hamiltonian Yu (1965); Shiba (1968); Rusinov (1968): , where is the exchange coupling and is the -component of the magnetic impurity , assumed as an Ising spin. In terms of the chiral fields, this term writes , where we have used Eq. (9) and the result . Again, can be absorbed via the gauge transformation , where is the -wave phase shift (note that we have used the definition of the density of states at the Fermi level )Balatsky et al. (2006). Then, in bosonic language we obtain the expression . In this way, the connection between and can be clearly seen: the 3D Hamiltonian is identical to the 1D Hamiltonian at the point with half the degrees of freedom (i.e., the field is absent, and the system becomes a helical liquid Wu et al. (2006). Consequently, the classical configuration of the field is also given by Eq. (8).

Results—We now turn to the resolution of Eq. (8). To that end, we write the fields as

(11) |

which allows to eliminate the term . Here is a continuous and smooth function at . We have found that the soliton/antisoliton-like functions ()

(12) |

are exact solutions of Eq. (8). Here, is an integer which determine the different minima of the potential, and is a parameter which must be chosen so that the continuity condition at the origin is verified. From this condition we obtain the equation

(13) |

from where is obtained. Once the value of and the nature of the solution (i.e., soliton or anti-soliton ) are fixed, the integer becomes restricted to only two possibilities since the left hand side of Eq. (13) can only vary between and .

Of special importance here is the physical meaning associated to the topological charge of the solitonic solutions (11) and (12). Replacing them into Eqs. (3) and (4), we obtain and . In addition, replacing them into the classical part of the Hamiltonian (5)-(7) allows to obtain their associated energy, i.e., , where , where we have substracted a constant background term so that the ground-state energy for (i.e., decoupled impurity) is zero. The numbers and determine a particular “energy branch” associated to that solution (see Fig. 1). In the limit (or ) we obtain the analytical result

(14) |

where we have defined .

It is instructive to analyze first the results for a Shiba impurity in a 3D host. For a decoupled impurity (i.e., when ), the ground-state configuration of the superconductor is described by a constant , which trivially minimizes the pairing term . This situation physically corresponds to the “classical” BCS ground state. Our solutions (11) and (12) precisely recover this behavior for =, for which Eq. (13) yields , and Eq. (14) yields the energy (see continous blue line in Fig. 1(a)). This ground state has topological numbers and (i.e., even-parity singlet), as expected physically for a BCS superconductor. On the other hand, the solution with () and corresponds to a topological kink (anti-kink) connecting the minima at and () at , and corresponds to an excited state with () and () (odd-parity doublet). We interpret this state as having an extra (one less) fermion with respect to the BCS ground state and has energy (dashed purple line in Fig. 1(a)). Other excited branches are shown in that figure. As the value of is increased, the branches and approach to each other and eventually cross at the value , signalling a quantum phase transition (i.e., the transition) at which the ground state abruptly changes spin and parity. Interestingly, our formalism allows to exactly recover the critical value of the transition obtained from the crossing of Shiba levels (where ) Shiba (1968); Balatsky et al. (2006) (note that precisely when , we obtain ).

In the case of the 1DNW, one can immediately determine the quantum phase diagram exploiting the degeneracy of the solutions at . The ground-state energy is given by the sum of the lowest degenerate branches at a particular value of . For , the ground-state energy is given by and corresponds to an even-parity singlet phase with topological numbers and (see blue line in Fig. 1(b)). On the other hand, for the ground state is magnetic (topological number ) and has even parity (). Therefore, in this case the value indicates a “0-0 transition” that preserves parity but not the spin, in stark contrast to the 3D case. This result can be understood due to the presence of two independent bosonic modes in a 1D geometry. Physically, each one of these modes are related to the fermionic bilinears and , each one capable of binding one Shiba state with at the critical point. Therefore, in the 1D geometry the transition occurs when these two independent Shiba states simultaneously cross the Fermi energy and become occupied. Thus, the spin quantum number of the ground state jumps from to , whereas the parity remains unchanged. Experimentally, this effect would imply the absence of sign-reversal of a supercurrent flowing through the 1DNW.

Note that for this effect to occur it is crucial to avoid the mixing of the fields . To that end, backwards scattering terms in the Hamiltonian (7) should be absent. Although in principle backscattering terms are present in the expression of the spin density , due the assumption of a mesoscopic FM contact of width this mechanism can be strongly suppressed. Note that this would not be the case in the more usual superconductor-quantum dot-superconductor geometry. Another requirement for the observation of the 0-0 transition is that the system obeys . Away from this point the fields become mixed, enabling the study of potentially interesting phenomena within the theoretical framework of bosonization.

###### Acknowledgements.

This work is partially supported by ANPCyT, Argentina. AML acknowledges the use of Toko Cluster facilities from FCEN-UNCuyo, which is part of the SNCAD-MinCyT, Argentina.## References

- Yu (1965) L. Yu, Acta Phys. Sin. 21, 75 (1965).
- Shiba (1968) H. Shiba, Prog. Theor. Phys. 40, 435 (1968).
- Rusinov (1968) A. I. Rusinov, Zh. Eksp. Teor. Fiz. Pisma. Red. 9, 146 (1968), [JETP Lett. 9, 85 (1969)].
- Balatsky et al. (2006) A. V. Balatsky, I. Vekhter, and J.-X. Zhu, Rev. Mod. Phys. 78, 373 (2006).
- De Franceschi et al. (2010) S. De Franceschi, L. Kouwenhoven, C. Schönenberger, and W. Wernsdorfer, Nature Nanotechnology 5, 703 (2010).
- Deacon et al. (2010) R. S. Deacon, Y. Tanaka, A. Oiwa, R. Sakano, K. Yoshida, K. Shibata, K. Hirakawa, and S. Tarucha, Phys. Rev. Lett. 104, 076805 (2010).
- Martín-Rodero and Yeyati (2011) A. Martín-Rodero and A. L. Yeyati, Advances in Physics 60, 899 (2011).
- Maurand et al. (2012) R. Maurand, T. Meng, E. Bonet, S. Florens, L. Marty, and W. Wernsdorfer, Phys. Rev. X 2, 011009 (2012).
- Nayak et al. (2008) C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das Sarma, Rev. Mod. Phys. 80, 1083 (2008).
- Nadj-Perge et al. (2013) S. Nadj-Perge, I. K. Drozdov, B. A. Bernevig, and A. Yazdani, Phys. Rev. B 88, 020407 (2013).
- Klinovaja et al. (2013) J. Klinovaja, P. Stano, A. Yazdani, and D. Loss, Phys. Rev. Lett. 111, 186805 (2013).
- Braunecker and Simon (2013) B. Braunecker and P. Simon, Phys. Rev. Lett. 111, 147202 (2013).
- Pientka et al. (2013) F. Pientka, L. I. Glazman, and F. von Oppen, Phys. Rev. B 88, 155420 (2013).
- Kitaev (2001) A. Y. Kitaev, Phys. Usp. 44, 131 (2001).
- Nadj-Perge et al. (2014) 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).
- Pawlak et al. (2016) R. Pawlak, M. Kisiel, J. Klinovaja, T. Meier, S. Kawai, T. Glatzel, D. Loss, and E. Meyer, npj Quantum Information 2, 16035 (2016).
- Ruby et al. (2015) M. Ruby, F. Pientka, Y. Peng, F. von Oppen, B. W. Heinrich, and K. J. Franke, Phys. Rev. Lett. 115, 197204 (2015).
- Sakurai (1970) A. Sakurai, Prog. Theor. Phys. 44, 1472 (1970).
- Franke et al. (2011) K. J. Franke, G. Schulze, and J. I. Pascual, Science 332, 940 (2011).
- Bauer et al. (2013) J. Bauer, J. I. Pascual, and K. J. Franke, Phys. Rev. B 87, 075125 (2013).
- Lee et al. (2012) E. J. H. Lee, X. Jiang, R. Aguado, G. Katsaros, C. M. Lieber, and S. De Franceschi, Phys. Rev. Lett. 109, 186802 (2012).
- Schindele et al. (2014) J. Schindele, A. Baumgartner, R. Maurand, M. Weiss, and C. Schönenberger, Phys. Rev. B 89, 045422 (2014).
- Island et al. (2017) J. O. Island, R. Gaudenzi, J. de Bruijckere, E. Burzurí, C. Franco, M. Mas-Torrent, C. Rovira, J. Veciana, T. M. Klapwijk, R. Aguado, and H. S. J. van der Zant, Phys. Rev. Lett. 118, 117001 (2017).
- Krogstrup et al. (2015) P. Krogstrup, N. L. B. Ziino, W. Chang, S. M. Albrecht, M. H. Madsen, E. Johnson, J. Nygård, C. . M. Marcus, and T. S. Jespersen, Nature Materials 14, 400 (2015), article.
- Taupin et al. (2016) M. Taupin, E. Mannila, P. Krogstrup, V. F. Maisi, H. Nguyen, S. M. Albrecht, J. Nygård, C. M. Marcus, and J. P. Pekola, Phys. Rev. Applied 6, 054017 (2016).
- Fisher and Kane (1992) M. P. Fisher and C. Kane, Phys. Rev. Lett. 68, 1824 (1992).
- Kane and Fisher (1992) C. L. Kane and M. P. A. Fisher, Phys. Rev. B 46, 15233 (1992).
- Giamarchi (2004) T. Giamarchi, Quantum Physics In One Dimension (Oxford University Press, 2004).
- Gogolin et al. (1988) A. O. Gogolin, A. A. Nersesyan, and A. M. Tsvelik, Bosonization and strongly correlated systems (Cambridge, 1988).
- Takei et al. (2013) S. Takei, B. M. Fregoso, H.-Y. Hui, A. M. Lobos, and S. Das Sarma, Phys. Rev. Lett. 110, 186803 (2013).
- Sau et al. (2010) J. D. Sau, S. Tewari, R. M. Lutchyn, T. D. Stanescu, and S. D. Sarma, Phys. Rev. B 82, 214509 (2010).
- Jespersen et al. (2009) T. S. Jespersen, M. L. Polianski, C. B. Sørensen, K. Flensberg, and J. Nygård, New Journal of Physics 11, 113025 (2009).
- Albrecht et al. (2016) S. M. Albrecht, A. P. Higginbotham, M. Madsen, F. Kuemmeth, T. S. Jespersen, J. Nygård, P. Krogstrup, and C. M. Marcus, Nature 531, 206 EP (2016).
- Haldane (1981) F. D. M. Haldane, Phys. Rev. Lett. 47, 1840 (1981).
- (35) T. Bortolin, A. Lobos, and A. Iucci, In preparation.
- Wu et al. (2006) C. Wu, B. A. Bernevig, and S.-C. Zhang, Phys. Rev. Lett. 96, 106401 (2006).