Quantum Phase Transition Induced by Magnetic Impurity

The Yu–Shiba–Rusinov bound states can arise around magnetic impurities in conventional superconductors. Competition between screened and unscreened states can lead to the quantum phase transition. We discuss basic properties of this quantum phase transition in the case of the triangular lattice with on-site pairing. We show that the quantum phase transistion results from the change of the ground state of the system.


Introduction
Interplay between impurities and superconductors is the focus of many studies [1]. Bound states near magnetic impurities in the conventional s-wave superconductor are an example of such interplay. This behavior was proposed by Yu [2], Shiba [3] and Rusinov [4] in the 1960s. In those pioneering papers, authors describe the possibility of pair breaking by a single magnetic impurity, that leads to emergence of in-gap quasi-particles states.
Recent experiments [5][6][7] show the possibility of realization of this bound state, called Yu-Shiba-Rusinov (YSR) Andrzej [8,9]. A change in the coupling strength between the magnetic impurity and superconductor, can lead to the quantum phase transition (QPT) [10]. The main idea of this QPT is shown schematically in Fig. 1. When the coupling is weak, there are numerous bound Cooper pairs (left part). However, after QPT, when the coupling is strong enough, the local magnetic field, generated by the impurity, leads to pair breaking and formation of bound broken pairs (right panel).
In this paper, we discuss basic properties of the QPT in the case of triangular lattice by studying the evolution of YSR states using the Bogoliubov-de Gennes (BdG) equations. Our microscopic model and the BdG technique are shown in Section 2. Numerical results and discussion are presented in Section 3, and we summarize our work in Section 4.

Model and Methods
We study a 2D triangular lattice, described by attractive Hubbard model (U < 0). The Hamiltonian is given as: Here is the free-electron kinetic term. c † iσ (c iσ ) denotes creation (annihilation) of electron with spin σ = {↑, ↓} at i-th lattice site, t is the hopping integral between nearest neighbors and μ is the chemical potential. For simplicity, we neglect the local modification of chemical potential at impurity site. For a large spin number S, the magnetic impurity can be treated classically In this case, the potential J plays the role of local magnetic field at the impurity site.
On-site pairing is given by the Coulomb term This part can be decoupled by means of the mean-field approximation, which gives: where we define i = c i↓ c i↑ as the superconducting order parameter (SOP). For the inhomogeneous system, the Hamitonian (1) can be diagonalized by Bogoliubov-Valatin transformation [11] where γ n and γ † n are the fermionic quasi-particle operators, whereas u and v are the BdG eigenvectors. Then, the selfconsistent BdG equations in real space take the form where where N is the total number of states.
The BdG (7) can be solved self-consistently, with respect to the SOP, which can be found as where f (ω) = 1/ (1 + exp(ω/k B T )) is the Fermi-Dirac distribution. Similar, the local density of states (LDOS) for particles with given spin σ in i-th site [12] can be calculated from whereas (total) density of states (DOS) ρ(ω) = iσ ρ iσ (ω).

Numerical Results and Discussion
In this section, we turn our attention to numerical results for a 2D triangular lattice with N x × N y = 41 × 41 sites in real space. In our calculations we assume U/t = −3, μ/t = 0 and k B T /t = 0.
In the first step, we calculate the total DOS as a function of J (Fig. 2a) (|ω|/t < 0.4). For any value of |J |, eigenvalues corresponding to those bound states are located symmetrically around ω/t = 0. Different spectral weight and opposite spin polarisation of positive and negative eigenvalues stems from the fact that YSR states are induced by magnetic impurity. In this sense, we can tell about spin polarized bound states [30], corresponding to nonzero value of ρ i↑ (ω) − ρ i↓ (ω).
As we can see, in some critical value of J (show by yellow dashed line), the polarization of bound state are changed. This critical values J c correspond to QPT induced by magnetic impurity. One can observe that when J → J c eigenvalues corresponding to YSR states are approaching ω/t = 0. Moreover, during the QPT energy of bound states changes discontinuously.
Other signature of the QPT in this system is the sign change of the SOP at the impurity site imp (Fig. 2b). On other hand, like in the case of a scalar impurity, increasing J leads to a decrease of | imp |, i.e., when J /t → ±∞ then | imp | → 0.
When describing a QPT induced by magnetic impurity, it is important to look into the total energy of the system, by calculating the grand canonical potential (8). For fixed parameter, the energy is only a function of the spatial decomposition of SOP in real space, i.e., ≡ ( ). Here, = ( 1 , 2 , · · · ) is a N x ×N y -component vector of the SOP at every site of real space. To compare the of a system before and after QPT, we calculate the grand canonical potential in two cases: when the SOP at the impurity site is (i) positive and (ii) negative. We denote the grand canonical potential in these two cases¯ + and¯ − , respectively. The comparison of those quantities is shown in Fig. 3, where the QPTs are marked by black arrows. As we can see the minimum energy of the system is given by¯ + for |J | < |J c |, while in other cases by¯ − .
The QPTs induced by magnetic impurities are similar to the classical phase transitions e.g. from superconducting to normal state. On one hand, energy of the system is given by a global minimum of the energy of different states (bound Cooper pairs or broken pair bound states). When a change of ground state is more favorable energetically, we observe a phase transition. On the other hand, QPTs are discontinuous. This can be deduced from the fact that the eigenvalues of YSR states change discontinuously and from the sign change of the SOP at the impurity site.

Summary
The magnetic impurity in a conventional superconductor can lead to a quantum phase transition. In this paper, we described this possibility in the case of triangular lattice with on-site pairing. We described basic properties of the quantum phase transition induced by the the magnetic impurity, i.e., the change of polarization of the bound states, discontinuous change of energy of bound states induced by impurity and the sign change of the superconducting order parameter at the impurity site. We have shown that the quantum phase transition is connected to a change of ground state of the system, which chooses the more energetically favorable state.