Emergent QCD Kondo effect in two-flavor color superconducting phase

We show that effective coupling strengths between ungapped and gapped quarks in the two-flavor color superconducting (2SC) phase are renormalized by logarithmic quantum corrections. We obtain a set of coupled renormalization-group (RG) equations for two distinct effective coupling strengths arising from gluon exchanges carrying different color charges. The diagram of RG flow suggests that both of the coupling strengths evolve into a strong-coupling regime as we decrease the energy scale toward the Fermi surface. This is a characteristic behavior observed in the Kondo effect, which has been known to occur in the presence of impurity scatterings via non-Abelian interactions. We propose a novel Kondo effect emerging without doped impurities, but with the gapped quasiexcitations and the residual SU(2) color group intrinsic in the 2SC phase, which we call the 2SC Kondo effect.


I. INTRODUCTION
The Kondo effect gives rise to rich physics from an emergent strong-coupling regime in the low-energy dynamics. In 1964, Jun Kondo pointed out that a long-standing problem about an anomalous behavior in the resistivity of alloy originates from the quantum corrections to the impurity-scattering amplitudes [1]. Namely, the coupling strength between conduction electrons and an impurity is enhanced due to contributions of the next-to-leading order scattering processes. The modification of the coupling strength is best captured by the concept of the renormalization group (RG), which inspired subsequent developments in this important concept [2][3][4]. The scaling argument clearly indicates that the existence of a Fermi surface is crucial for the Kondo effect to occur [5,6] as we will explain in a brief review part in the next section.
In recent years, the Kondo effect was applied to nuclear physics [7,8]. Especially, possible realization of the QCD Kondo effect in dense quark matter was discussed when heavy-quark impurities are embedded in light-quark matter [8]. The results of the perturbative RG analyses indicate that the effective interaction strength between light and heavy quarks evolves into a Landau pole despite the small value of the QCD coupling constant at high density. Subsequent studies investigated further consequences of the QCD Kondo effect, including interplay/competition with color superconductivity [9], formation of "Kondo condensates" and modification of the QCD phase diagram [10][11][12][13], nonperturbative aspects near the IR fixed point by conformal field theory [14,15], estimates of transport coefficients [16], and QCD equation of state for an application to neutron/quark star physics [17]. It is also remarkable that light quarks have the same scaling dimensions in a dense system and in a strong magnetic field as a consequence of analogous effective dimensional reductions [6], and that a strong magnetic field alone induces the QCD Kondo effect even at zero density [18].
In this paper, we propose a realization of the Kondo effect without (doped) impurities in QCD. We consider the RG evolution of effective coupling strengths between gapped and ungapped quarks appearing in the two-flavor color superconducting (2SC) phase [19,20]. In the 2SC phase with three colors, two of three color states of quarks are involved in the Cooper pairing and acquire a gap above the Fermi surface (cf. Fig. 3), while the other color state remains gapless with a finite density of states at the Fermi surface (see, e.g., Refs. [21,22] for review articles). Therefore, the gapped quarks may play a role of impurities in the low-energy dynamics below the gap size, and the 2SC phase intrinsically has the necessary setup of the Kondo effect.
The formation of the diquark condensates breaks the SU(3) color symmetry group down to SU (2). Gluons belonging to the unbroken SU (2) subgroup are coupled only to the gapped quarks, so that they are decoupled from all quarks in the aforementioned low-energy regime, realizing a pure gluodynamics [23]. Here, we can focus on the remaining five gluons which mediate the interactions between ungapped and gapped quarks. Those gluons have different properties depending on if they are associated with the diagonal or off-diagonal Gell-Mann matrices, so that we introduce two distinct effective coupling strengths.
We will derive coupled RG equations for these two effective coupling strengths from the next-to-leading order scattering amplitudes (cf. Fig. 4) and obtain an RG-flow diagram in Fig. 5. The fate of the RG flow depends on the initial conditions for the RG equations. We take the treelevel coupling strengths as initial conditions which depend on the magnitudes of the Debye and Meissner masses in the 2SC phase [24,25]. Plugging the initial conditions evaluated with those quantities, we find that both of the coupling strengths evolve into strong-coupling regimes, but have different signs for attraction and repulsion.
This paper is organized as follows. In the next section, we provide a brief review of the QCD Kondo effect which will be useful to identify the basic mechanism of the Kondo effect. Then, we investigate a novel Kondo effect occurring in the 2SC phase in Sec. III. We conclude the paper in Sec. IV with discussions. Some useful properties of the high-density and heavy-quark effective field theories are briefly summarized in the Appendix.

II. BRIEF REVIEW OF QCD KONDO EFFECT
We first provide a brief review of the QCD Kondo effect [8], highlighting the essential points of the discussions given in a review article [6]. These preliminary discussions will be useful to identify the necessary ingredients of the Kondo effect.

A. High-density effective field theory
The Kondo effect is induced by low-energy excitations near the Fermi surface scattering off dilute impurities. Therefore, we use the high-density effective field theory (HD-EFT) to extract such low-energy degrees of freedom (d.o.f.) from the full theory [26][27][28][29][30].
We decompose the fermion momentum into a sum of the large Fermi momentum and the small residual momentum l μ ¼ ðl 0 ; ℓ Þ as where μ and v F are the chemical potential and the Fermi velocity, respectively, and l 0 , jlj ≪ μ. The energy p 0 is measured from the Fermi surface. Then, at the leading order (LO) of 1=μ expansion, the low-energy d.o.f. are extracted from the Dirac Lagrangian as (see the Appendix) where v μ FAE ≡ ð1; AEv F Þ. We introduced the low-energy field ψ þ ðl; v F Þ ≡ P þ ψðl; v F Þ for particle and hole excitations around the Fermi momentum μv F , by the use of projection operators P AE ≡ 1 2 ð1 AE γ 0 v F · γÞ. From this expansion, the dispersion relation of the lowenergy excitations near the Fermi surface is read off as We define l k ≡ v F · ℓ for later use. This is a linear dispersion relation in the (1 þ 1)-dimensional phase space normal to the Fermi surface and does not depend on the residual two-dimensional momentum tangential to the Fermi surface (cf. Fig. 1). This means that an effective dimensional reduction occurs in the low-energy excitation near the large Fermi sphere, and the phase space is degenerated in the residual two dimensions. From Eq. (2), the free propagator is found to be where we have = v FAE ¼ 2P AE γ 0 . At the LO, only the temporal and the parallel components of the gauge field, A 0 and v F · A, are coupled to the low-energy fermion excitations. The gamma matrix is not involved in these couplings, because the spin direction is frozen along the Fermi velocity.

B. Scaling dimensions in QCD Kondo effect
One can determine the scaling dimension of the lowenergy excitation field ψ þ assuming that the kinetic term (2) is invariant under the scaling transformation, l 0 → sl 0 (t → s −1 t) with s < 1. According to the (1 þ 1)-dimensional dispersion relation, only the l k scales as l k → sl k , and the tangential momentum ℓ ⊥ is intact. Therefore, the ψ þ scales with a factor of s −1=2 when the energy scale is reduced toward the Fermi energy (l 0 ¼ 0). 1 FIG. 1. Effective dimensional reduction near a large Fermi surface. Excitations cost energy only in the normal direction, so that there is a two-dimensional degeneracy in the tangential directions. 1 Since there is a degeneracy in the momentum space, it is useful to count the scaling dimensions in a mixed representation ψ þ ðt; ℓ Þ.
In addition, we introduce a heavy-quark impurity embedded in the dense light-quark matter. One may use the heavy-quark effective field theory (HQ-EFT) which is organized with an expansion with respect to the inverse heavy-quark mass 1=m H . This expansion is analogous to that in the HD-EFTas we summarize in the Appendix. At the LO, the kinetic term of the particle state Ψ þ is given by We consider a static heavy quark with a vanishing spatial velocity. In this case, the spatial derivative is not contained in the kinetic term, and we find that the heavy-quark field Ψ þ and its spatial momentum k do not scale when t → s −1 t. This is reasonable since the static impurity acts as a scattering center in the same way at any energy scale of light particles. Now that we have determined the scaling dimensions of both the light and heavy quark fields, we look into the four-Fermi operator composed of the light and heavy quarks: Plugging the scaling dimensions of the fields and of the momenta discussed above, we find that the light-heavy four-Fermi operator has a marginal scaling dimension . This result suggests that the four-Fermi interaction acquires logarithmic quantum corrections from the scattering of the light quark off the heavy-quark impurity with loop diagrams. We will see how the logarithmic enhancement arises from the second-order heavy-light scattering amplitude and determine the sign of the logarithmic correction. The following computation by the HD-EFT and HQ-EFT confirms the result in Ref. [8].
Note again that the lower scaling dimension of ψ þ due to the effective dimensional reduction is crucial for the Kondo effect. It is worth mentioning that the BCS instability, leading to superconductivity, can be understood as a consequence of the same dimensional reduction (see, e.g., Refs. [5,26,27]). The RG analyses were performed for color superconductivity [31][32][33][34][35]. Also, the scaling argument can be applied to the low-energy dynamics in a strong magnetic field where an analogous effective dimensional reduction occurs in the lowest Landau level [6]. Consequences of this dimensional reduction are known as the magnetic catalysis of chiral symmetry breaking [36][37][38] and the magnetically induced QCD Kondo effect [18].
C. Effective interaction and the leading-order scattering amplitude The gluon propagator at high density is split into two transverse structures: where ξ is a gauge parameter and the diagonal color structure is suppressed for notational simplicity. The longitudinal and transverse projections are given by In this section, we consider a normal phase (without Cooper pairing) and use the gluon self-energy from the hard dense loop approximation [30,39,40]. In the heavyquark limit, we only need the electric component which is screened by the Debye screening mass We consider a light quark scattering off a static heavy quark and identify the effective coupling G in Eq. (6) with the S-wave projection of the color electric interaction [6,18]: We have put the initial and final momenta of the light quark on the Fermi surface since the Kondo effect occurs in such a low-energy regime. This leads to a spacelike momentum transfer k 2 ¼ ðp ð3Þ − p ð1Þ Þ 2 ¼ −2μ 2 ð1 − cosθÞ, and we have integrated out the scattering angle θ ≡ arccosðv F =μ 2 Þ. The gauge term proportional to ξ is suppressed in this kinematics. A similar S-wave projection was performed for the Cooper pairing [32,33,35].
Using the effective coupling (9), the leading-order scattering amplitude is given by where N c is the number of colors. For a notational simplicity, we have suppressed the spinor structure, ½ū þ γ 0 u þ ½Ū þ U þ with the projected spinors u þ ¼ P þ u and U þ ¼ Q þ U for the light and heavy quarks, respectively.

D. Kondo scale emerging from the RG evolution
Since the four-Fermi operator, composed of light and heavy quarks, has a marginal scaling dimension, we anticipate emergence of a dynamical infrared scale. Below, we shall see how the logarithmic correction from the loop integral drives the RG evolution of the effective coupling to an infrared Landau pole. At the one-loop level, there are two relevant diagrams for the Kondo effect (cf. Fig. 2). In terms of the effective coupling G, the propagators (4) and (A13), those amplitudes are written down as where the overall minus signs come from the fermionic statistics and S H is the heavy-quark propagator given in Eq. (A13). Again, we put the initial and final momenta of the light quark on the Fermi surface, i.e., l ðiÞ μ ¼ l ðfÞ μ ¼ 0, and suppressed the spinor structures which are found to be the same as that of the tree-level amplitude (10). We assume a static heavy quark with a vanishing velocity. As we will see shortly, it is important to have noncommutative color matrices Performing the integrals in Eq. (11), we have Note that, since the integral regions in Eq. (13) are restricted to the above and below of the Fermi surface, diagrams ð2aÞ and ð2bÞ of Fig. 2 provide a particle and hole contribution, respectively, as specified by the pole positions in Eq. (11). The density of states on the Fermi surface ρ F has been obtained from its area: We now examine an increment when the energy scale Λ is reduced to Λ − dΛ. The sum of the two one-loop amplitudes, integrated over a thin shell of a thickness dΛ, is obtained as The relative minus sign in the curly brackets originates from the fact that the particle and hole contributions in Eq. (13) have opposite signs. The logarithms from the two distinct diagrams would cancel each other if the interaction were an Abelian type. Therefore, the non-Abelian nature of the interaction plays an essential role for the logarithmic correction to survive in the total amplitude. By the use of an identify T ðaÞ ij;lm − T ðbÞ ij;lm ¼ −N c =2 P r ðt r Þ ij ðt r Þ lm , we find a logarithmic correction to the total one-loop amplitude: It is this logarithm that renormalizes the effective coupling G. Now, combining the results in Eqs. (10) and (16), we obtain the RG equation The solution to this RG equation is found to be FIG. 2. One-loop diagrams that give rise to logarithmic corrections. The single and double lines denote the light and heavy quarks. Both initial and final light-quark momenta are put on the Fermi surface, i.e., l μ ðiÞ ¼ l μ ðfÞ ¼ 0.
Here, Λ 0 is the initial energy scale, and the initial condition of G is given by the tree-level result in Eq. (9). Namely, It is clear that the effective coupling (18) is enhanced according to a negative beta function, when GðΛ 0 Þ > 0. We can read off the location of the Landau pole that is called the Kondo scale [6,8]: where we took the initial energy scale at the hard scale Λ 0 ¼ μ. When the temperature is reduced below the Kondo scale, the system becomes nonperturbative no matter how small the initial coupling GðΛ 0 Þ or α s is. Intuitively speaking, the light particles (or carriers of transport phenomena) are trapped around the impurity due to the strong-coupling nature of the low-energy dynamics. As a consequence, the (electrical) resistance is enhanced below the Kondo temperature T K of which the scale is given by T K ∼ Λ K , and there emerges a minima at T K .

III. 2SC KONDO EFFECT
We here highlight four important ingredients for the Kondo effect discussed in the last section. First of all, the Kondo effect needs (i) impurities. Next, as implied by the essential d.o.f. in HD-EFT and the scaling argument, the (1 þ 1)-dimensional dispersion relation plays a crucial role. Therefore, the Kondo effect needs (ii) the existence of the Fermi surface so that the effective dimensional reduction occurs in the low-energy dynamics. For this low dimensionality, the four-Fermi operator, composed of the heavy and light particles, acquires a marginal scaling dimension. Consequently, the effective coupling strength is renormalized due to (iii) the logarithmic quantum correction from the loop integrals. Finally, the logarithms from the two distinct one-loop diagrams do not cancel out (iv) only when the interaction is a non-Abelian type.
Once these ingredients are identified, one may consider extensions of the Kondo effect. As already discussed, the QCD Kondo effect in dense quark matter is a straightforward extension since the fourth ingredient (iv) is provided by the color exchange interaction [8]. One may also replace the second ingredient (ii), the Fermi surface, by a strong external magnetic field which also causes an effective dimensional reduction in the low-lying state, i.e., the lowest Landau level [18].
In this section, we propose a Kondo effect which occurs without impurities. Namely, we do not introduce the most essential ingredient (i) externally, but investigate a situation in which the gapped excitations (i.e., the "impurities") emerge dynamically through a spontaneous symmetry breaking. The physical system we will consider is the 2SC phase of dense quark matter. The gapped quarks and the broken generators of the color symmetry will play a role of a heavy impurity and a non-Abelian interaction, respectively. We anticipate that a novel Kondo effect emerges with all the necessary ingredients inherent in the 2SC phase.

A. Gapped quarks
We have examined the quark propagator in the normal phase in the previous section. Here, we prepare a quark propagator for a gapped quark. Without losing generality, we hereafter choose blue quarks (i ¼ 3) to be ungapped ones, so that red and green quarks are gapped above the Fermi surface (cf. Fig. 3). Then, the color structure of the quark propagator reads where S and S Δ are the quark propagators with and without an energy gap, respectively. There would be off-diagonal components in the flavor space (f, g) if one considers interactions between the quasiexcitations and the hudi condensate in the 2SC phase. However, they are suppressed with a small value of the condensate. As shown in Refs. [24,41], the propagator of the gapped quasiexcitations is given by The dispersion relations read where the energy gaps Δ andΔ are generated as a consequence of the diquark and diantiquark condensate formation, respectively. As in Eq. (1) for the HD-EFT, we decompose the momentum into a large Fermi momentum and a small fluctuation near the Fermi surface. Then, we have jpj 2 ∼ μ 2 þ 2μl k . Therefore, we find jpj ∼ μ þ l k and At the leading order in the 1=μ expansion, the propagator reads The imaginary displacements are explicitly shown for the quasiparticle and quasihole excitations. This propagator has the same projection operator P þ as that of the ungapped quark (4), meaning that the highly suppressed antiparticle excitations are neglected and that the coupling to a gluon field is again simplified as in the leading-order Lagrangian of the HD-EFT (2). Namely, the gamma matrix is replaced by the Fermi velocity v μ F .

B. Two effective coupling strengths from screened gluon exchanges
The diquark condensation modifies not only the quark dispersion relation but also properties of the gluons. In the 2SC phase, the SUð3Þ color symmetry is broken down to the SUð2Þ subgroup, and the three gluons remain massless and do not interact with the ungapped blue quarks. They are also decoupled from the other two gapped quarks in the low-energy region below the size of the gap (cf. the yellow band in Fig. 3). Therefore, a pure gluodynamics of SUð2Þ theory was investigated in such a low-energy region [23]. Here, we consider the other five gluons which are screened by the Meissner mass ½Π T ð0Þ ¼ m 2 M as well as the Debye mass ½Π L ð0Þ ¼ m 2 D [24,25]. 2 Those gluons mediate the interactions between the gapped and ungapped quarks.
As we have seen in the last section, the energy gap of an impurity should be a dominant scale so that an impurity field is invariant under the scaling and plays the role of a heavy scattering center. In the current case, the gap size Δ should be the dominant scale for a gapped-quark field to be invariant under the scaling of an ungapped-quark field. This is realized in a scattering between a pair of gapped and ungapped quarks moving in almost the same direction within a relative angle Δ=μ, where the momentum scale of the ungapped quark is smaller than Δ in the comoving frame of the gapped quark. Therefore, we focus on those pairs whose incoming momenta are labeled by the same Fermi velocity v F . Then, dynamics is dominated by small fluctuations near the Fermi surface.
Since the gamma matrix on the interaction vertex is replaced by the four Fermi velocity v μ F , we identify the effective coupling strength from the S-wave scattering [cf. Eq. (9)]: As in the previous section, we have put the initial and final momenta of the ungapped quark on the Fermi surface, and integrated out the scattering angle θ. The magnitudes of m D and m M depend on the color index of gluons. Since the red and green quarks play the same role, we consider the scattering between the red (1) and blue (3) quarks. Then, there are only three relevant gluons A 4;5;8 which mediate the interactions between the red and blue quarks. 3 Quoting the results in Table I of Ref. [24], we have m M =m D ¼ 1= ffiffi ffi 3 p for A 4;5 and m M =m D ¼ 1=3 for A 8 in the 2SC phase at T ¼ 0. Therefore, we define two effective coupling strengths: It is useful to compare the setup in the 2SC phase to the anisotropic Kondo effect discussed in Refs. [2,3]. Our coupling strengths G 2SC andḠ 2SC play similar roles of J AE and J z there. A negative J z corresponds to the ferromagnetic Kondo problem of which the fate at low energy depends on the relative strength between J AE and jJ z j. In order to determine the low-energy physics, we will investigate how those effective couplings G 2SC andḠ 2SC are renormalized with the next-to-leading order scattering processes between gapped and ungapped excitations.

C. Color flows in the next-to-leading order scattering diagrams
We write the scattering diagrams in terms of the effective coupling constants defined above. At the leading order, we simply have where we defined a shorthand spinor notation U ab ≡ū a þ γ 0 u b þ . Next, we shall identify next-to-leading order diagrams which renormalize the effective coupling strengths. In Fig. 4, we show flows of color charges where the color matrices t 4;5 mix the red and blue quarks, while the diagonal color matrix t 8 does not. Similar to the diagrams in Fig. 2, we should consider a pair of diagrams ð2aÞ and ð2bÞ of Fig. 4. However, we do not need to consider the diagrams with the diagonal matrix t 8 on all the vertices, since possible logarithmic corrections will cancel out. On the other hand, we should include diagrams ð2 0 Þ and ð2 00 Þ of Fig. 4 which have each of t 4;5 and t 8 . Neither of these two diagrams has a relevant cross channel, since its cross channel is a disconnected diagram, indicating an annihilation between an ungapped particle and gapped hole or an ungapped hole and gapped particle. Diagrams ð2 0 Þ and ð2 00 Þ of Fig. 4 contain a factor of G 2SCḠ2SC , and give rise to a mixing between the two coupling strengths.
Similar to the previous section, the scattering amplitudes of diagrams ð2aÞ and ð2bÞ of Fig. 4 are written down as On the external lines, we put all the spatial momenta and energies of ungapped quarks on the Fermi surface. Only one finite component of the external momentum is the energy of a gapped quark, i.e., l ð2Þ μ¼0 ¼ l ð4Þ μ¼0 ¼ Δ. The density of states at the Fermi surface ρ F is obtained as in the previous section. The longitudinal components of the integrals are given by Here, we clearly see that Performing the contour integral with respect to l 0 and keeping only the singular terms when l k → 0, 4 we have Integrating over a shin shell Λ − dΛ ≤ l k ≤ Λ, we find a logarithmic contribution FIG. 4. Next-to-leading order scattering diagrams which renormalize the effective coupling constants. Single and double lines show ungapped and gapped quarks, respectively. 4 One can enclose the contour either in the upper or lower half planes. The results are the same in both cases.
The subsequent terms are a polynomial of Λ which provides only irrelevant corrections. Plugging this result into M 2 2SC , we obtain The remaining two diagrams can be computed in the same manner. Actually, those diagrams have the same kinematics as diagram ð2aÞ of Fig. 4, so that we only need to take care of the color structures. Writing down those two contributions, we have where the products of the color matrices read

D. Coupled RG equations and RG-flow diagram
We are now in position to derive RG equations for G 2SC andḠ 2SC . Plugging the leading-order amplitude (26) and the next-to-leading order amplitudes (33) and (36), we immediately obtain the coupled RG equations Correspondingly, the right-hand sides of the RG equations provide two distinct beta functions. In Fig. 5, we draw the RG flow driven by the "velocity field" identified with those beta functions.
To understand the RG-flow profile, we write the RG equations (37) as This means that the RG flow evolves along parabolic curves where the constant C is determined by the initial conditions at Λ ¼ Λ 0 . We take the tree-level coupling strengths (25) as the initial conditions and have a positive constant C > 0. In Fig. 5, we start out at a point G 2SC ðΛ 0 Þ > 0 and G 2SC ðΛ 0 Þ < 0, and we find that the RG flow goes into the lower right corner. This means that the G 2SC ðΛÞ and G 2SC ðΛÞ evolve toward positive and negative infinity, respectively, away from the weak-coupling regime near the origin. Thus, the RG evolution driven by the interaction between the gapped and ungapped quarks gives rise to a strongly coupled system in the low-energy regime. This is a characteristic behavior of the Kondo effect. From the RG equations (37), we get and its solution where c ¼ ffiffiffiffiffiffiffiffi ffi C=2 p . This solution hits a Landau pole when the argument of the tangent approaches π=2, giving rise to an emergent scale (see Ref. [38] for the same form of solution). According to the relation (39), the other coupling strength G 2SC ðΛÞ also hits the Landau pole at the same energy scale (cf. Fig. 5).
Noting that cρ F ∼ Oðg 2 Þ andḠ 2SC ðΛ 0 Þ=c ∼ Oð1Þ in Eq. (41), the parametric dependence of the Kondo temperature is found to be The explicit form of the order-one constant c 1 can be obtained easily, but is suppressed for a simple parametric estimate. It is natural to take the initial scale at the gap size, Λ 0 ∼ Δ ∼ μ expð−c 2 =gÞ, where c 2 ∼ Oð1Þ. Importantly, the gap size has a weaker exponential suppression because of the enhancement arising from the unscreened color-magnetic interaction [33,35,43,44]. However, the gluons involved in the 2SC Kondo effect are screened by both the Debye and Meissner masses. Therefore, in the weakcoupling limit g ≪ 1, we get a parametric estimate This hierarchy confirms our basic picture of the 2SC Kondo effect (cf. Fig. 3). Note that the two dynamical scales T 2SC K and Δ emerge in the different color sectors. The emergence of the strong-coupling regime implies formation of a bound state or condensation between the gapped and ungapped quarks. In condensed matter systems, such a bound state between the conduction electron and impurity has been known as the Kondo singlet [45][46][47]. While the impurity magnetic moment is localized in such systems, the gapped quarks are thermally excited in the bulk. Therefore, in the present system, we may think of it as a condensate rather than a localized bound state. Then, the formation of condensation breaks the residual SU(2) color symmetry, and the associated gluons become massive via the Higgs mechanism. 5 Nevertheless, whether such a "Higgs phase" emerges depends on the gapped-quark distribution which reduces as we decrease temperature.
At this moment, we conclude that the residual color symmetry is broken as long as the 2SC Kondo phase manifests itself in the QCD phase diagram. More generally, one may ask how a possible phase structure depends on the impurity distribution as an axis of an extended phase diagram (see Refs. [10,13] for the chiral symmetry breaking in the presence of a homogeneous distribution of the heavy-quark impurities in quark matter). We leave those issues to future works.

IV. CONCLUSIONS AND DISCUSSIONS
In this paper, we investigated the RG evolutions of the coupling strengths between gapped and ungapped quarks in the two-flavor superconducting phase. The next-to-leading order diagrams generate logarithmic quantum corrections and have the effective coupling strengths renormalized. We obtained coupled RG equations for the two coupling strengths associated with distinct color channels. The RG-flow diagram indicates that both of the coupling strengths evolve into a strong-coupling regime as the energy scale is reduced toward the Fermi energy. This is a characteristic behavior of the Kondo effect, so that we call it the "2SC Kondo effect." Once the system approaches the strong-coupling regime, the fate of the RG evolution needs to be investigated with nonperturbative methods. For example, a mapping from the Kondo problem in the vicinity of an infrared fixed point to conformal field theory has been known as a useful method (see recent works [14,15] and references therein). In the present system, one could ask if there is an infrared fixed point and what the ground state of the system is.
The magnitudes of the Kondo effect on bulk quantities, such as transport coefficients, depend on the concentration of impurities, i.e., the Fermi-Dirac distribution of the gapped quarks in the present case. At strict zero temperature, there are no gapped quark excitations. Therefore, effects of the 2SC Kondo effect will be most prominent in between the transition temperature of the 2SC phase and zero temperature. This contrasts with the conventional Kondo effect which remains important at zero temperature.
If the 2SC phase exists in the neutron-star cores (see Ref. [22]), we expect that the occurrence of the 2SC Kondo effect may have implications for the neutron-star physics. In particular, the transport properties of the neutron-star matter could be altered in a way similar to how the conventional Kondo effect affects the low-temperature resistivity of magnetic alloys. One typical example is the neutrino emissivity which is important as a cooling mechanism of neutron stars. It is conventionally known that, in the 2SC phase, the existence of ungapped quarks opens a large phase space for neutrino emission, leading to a too-fast cooling as compared to the astrophysical data [49]. The 2SC Kondo effect may serve as a mechanism to suppress the emissivity and give a cooling rate closer to the 5 The formation of a bound state/condensate may give rise to a pole in the Bethe-Salpeter amplitude (cf. an analogous structure in the Cooper pairing [48]), suggesting that the scattering ungapped quarks have off-shell momenta on the external lines. The natural scales of the off-shellness and binding energy are given by the Kondo scale (cf. Ref. [47]) which is the emergent and smallest scale in our analysis as shown in Eq. (43). The offshellness or deviation from the Fermi surface below the emergent scale do not change the form of the logarithm that drives the renormalization flow from the ultraviolet to the Landau pole. This observation, a posteriori, justifies our setup of the kinematics with the ungapped-quark momenta on the Fermi surface. data. This calls for a more detailed investigation in the future.
Besides, we would like to mention a new possibility that ultracold atoms could serve as a designed platform for studying quantum many-body physics. Realization of "color superconductivity" has been discussed with fermionic atoms carrying colorlike internal d.o.f. [50][51][52][53]. Those models, however, do not have "gluons." Although the 2SC Kondo effect does not necessarily need dynamical gluon exchanges, which can be substituted by contact interactions, a non-Abelian matrix on each interaction vertex is an indispensable ingredient. While Cooper pairing with a "color-flipping" effect was recently discussed [54], a noncommutative property is yet more demanding. Nevertheless, the direction for ultracold atoms deserves further study. It may also be worth mentioning that realization of the Kondo effect [55][56][57][58] and Shiba state [59][60][61], a mixture of impurity and superconducting states, was discussed in terms of ultracold atoms.

APPENDIX A: EFFECTIVE FIELD THEORIES
In this appendix, we briefly summarize the basic properties of the high-density and heavy-quark EFTs at the leading order of the expansions with respect to a large chemical potential 1=μ and a heavy-quark mass 1=m H , respectively.

High-density EFT
For a given Fermi velocity, the corresponding plane wave is factorized as where ψ v F ðlÞ is the field for low-energy excitations (in the momentum space). Plugging the above decomposition into the Lagrangian, one gets where v μ FAE ≡ ð1; AEv F Þ. The kinetic term in Eq. (A2) yields not only the one in Eq. (2) but also the other three termsψ where ψ AE ðx; v F Þ ≡ P AE ψðx; v F Þ, γ μ k ≡ ðγ 0 ; ðv F · γÞv F Þ, and γ μ ⊥ ≡ γ μ − γ μ k . From Eq. (A3a), the antiparticle states are gapped by 2μ, so that those excitations are highly suppressed in the dense system. When l ⊥ is smaller than the gap 2μ in Eq. (A3b), the mixing between the particle and antiparticle states is also suppressed.
Here are some basic properties of the projection operators: By using the identities, we find

Heavy-quark EFT
We also briefly summarize the heavy-quark effective field theory at the leading order [62]. We shall decompose the heavy-quark momentum as where the velocity is normalized as v 2 H ¼ 1. Since excitations of the antiparticle states are highly suppressed by 1=m H , it is natural to introduce operators which, in the rest frame (v H ¼ 0), project out the particle and antiparticle states. Here are some basic properties of the projection operators To get these identities, we used Assuming that the on-shell momentum m H v μ H does not change in the low-energy dynamics, we factorize the corresponding plane wave as Ψ AE ðxÞ ≡ e im H v μ H x μ Q AE ΨðxÞ: Then, we have Therefore, by using the identities (A8), the kinetic term of the heavy quark is decomposed as At the leading order, the interaction vertex does not contain the gamma matrix, because the magnetic moment is suppressed by 1=m H . The free HQ propagator is read off as