Arxiv:2107.00455V2

Home , Kondo effect

arXiv:2107.00455v2 [cond-mat.str-el] 31 Aug 2021 uino h yrdzto ed eodtemean-field the beyond fields hybridization distri- the probabilistic of the bution simulate directly to proximation of physics. hybridization exploration the our of limits fluctua- richness sign severely hybridization the which the of lacking, to treatment is due tions proper not A fields do hybridization problem. usually the difficulties [44–56] probe analytical simulations directly to expan- Numerical due order hybridization low [37–43]. the only of include effective sions some but ther- to theories led some have account low-energy fluctuations into quantum take and to Kondo artificial mal Efforts an multi-impurity 35]. as in [34, systems such transition predictions mean-field phase false The (first-order) made has [36]. theory experiments ultrafast pump-probe in observed optical rise recently as give physics may new that important ignore to fields and hybridization [22–35] the of approximation fluctuations fo- mean-field studies the most on However, cus and [19–21]. by pseudofermions electrons described between conduction been impurity hybridization has the effective physics for an underlying the used and often spins pseud- is a representation Theoretically, ofermion systems. correlated studied well corre- strongly in applied the systems. widely but electronic been [5–18], lated not parameter has order identify bro- the concept to the and introduced of symmetry been knowledge ken recently explicit without has subsys- transitions It between phase correlations 4]. of [3, amount to tems used total be the may and probe variables random two depen- between statistical concept the dency key measures a theory, information, information the mutual microscopic in The of addition 2]. simple [1, a properties as viewed phenom- be collective cannot whose ena systems quantum many-body nti ok epooeasai uiir edap- field auxiliary static a propose we work, this In most the of one arguably are systems Kondo The in role key a plays subsystems between Correlation uulifrain unu hs rniin n hs co phase and transition, phase quantum information, Mutual utain n ofimtemta nomto saueu to useful a systems. as Kondo hi information in mutual These the properties. confirm and physical fluctuations on intersi influences of observable development with “Varm the so-called manifesting the information near mode mutual information Kondo phas mutual two-impurity (first-order) amplitude the gent artificial to applied the When overcomes spat and approximation. the of orders account all full takes at informat method mutual Our the use correlations. and tion problem sign the without systems epooeasai uiir edapoiaint td th study to approximation field auxiliary static a propose We 1 ejn ainlLbrtr o odne atrPyisa Physics Matter Condensed for Laboratory National Beijing 3 oghnLk aeil aoaoy ogun undn 5 Guangdong Dongguan, Laboratory, Materials Lake Songshan inJnDong, Jian-Jun 2 nvriyo hns cdm fSine,Biig100049, Beijing Sciences, of Academy Chinese of University hns cdm fSine,Biig109,China 100190, Beijing Sciences, of Academy Chinese ,2 1, oghnHuang, Dongchen Dtd etme ,2021) 1, September (Dated: systems edter a rdce rtodrtasto between mean- transition the first-order a While predicted has electrons. theory field conduction by screened are pnsnlt hl nteKnolmt( limit Kondo the in while singlet, spin a we ieetmn-oygon tts ntemag- the In ( fea- states. be- limit ground netic that competition many-body systems and different correlations simplest tween electronic the strong of ture one describes model sflisgto hi yrdzto physics. hybridization provide their to Our on models insight other point. to useful critical extended physi- easily the be on near can method particular influences phase in intersite important properties, the have cal of These measure pro- a time distance, coherence. first impurity the large Kondo for the for viding in value diminishes finite but a point, adopts regime phase fixed the “Varma-Jones” of that so-called while the tem- lowering estima- near with neural perature divergence logarithmic developed a exhibits recently tor, amplitude, the hybridization using the that calculated find of We information mutual [57–59]. the analysis renor- (NRG) numerical group exact tran- the malization with phase Kondo agreement first-order two-impurity good artificial in the sition, the to suppresses it applied As model, when problem. sign example, the using without an fields sampling hybridization the Carlo the Monte of capture the of to aspects correlations us important spatial allows be some full method reveal may The to information physics. tool hybridization mutual a that as show used We approach. where respectively, J s i H esatwt h w-muiyKnomodel, Kondo two-impurity the with start We = ,2 1, r h od n esnegecag interactions, exchange Heisenberg and Kondo the are H P epaechrnei h od regime, Kondo the in coherence phase te n ifn Yang Yi-feng and a utain ftehbiiainfields hybridization the of fluctuations ial = ǫ rniinpeitdi h mean-field the in predicted transition e hih h motneo hybridization of importance the ghlight αβ k -oe”fie on n ag phase large a and point fixed a-Jones” lt xlr h yrdzto physics hybridization the explore to ol X o omaueteitriehybridiza- intersite the measure to ion k stecnuto lcrndispersion, electron conduction the is c σ iα † ,i eel oaihial diver- logarithmically a reveals it l, J yrdzto hsc fKondo of physics hybridization e ǫ H S ~ σ k αβ 2 i c ≫ k † dIsiueo Physics, of Institute nd steiprt pnlctdat located spin impurity the is c σ iβ J c k K σ sta fcnuto lcrn.The electrons. conduction of that is ,tetoiprte r okdinto locked are impurities two the ), 30,China 23808, + J ,2 3, 2, 1, China K eec nKondo in herence X i =1 2 ∗ s i · S i + J J H K S 1 ≫ · S J 2 H R J , ,both ), K i and , and (1) 2 two phases [34, 35], NRG and the conformal field theory (a) 0.04 (b) 0.6

(CFT) suggest a crossover at ﬁnite temperatures and a χ |V|=0, | |=0 χ 0.03 P(V1,V2, ) “Varma-Jones” ﬁxed point at zero temperature for a con- 0.4 P(V ,V )

| 1 2

duction band with the particle-hole symmetry [57–61]. 1

T 0.02 The hybridization physics is best seen in the Abrikosov |V Mean Field ~σηγ 0.2 Si pseudofermion representation, = ηγ fiη† 2 fiγ . The 0.01 |V|=0, |χ|≠0 |V|≠0, |χ|≠0 Kondo and Heisenberg terms can beP decomposed using the standard Hubbard-Stratonovich factorization, result- 0 0 ing in a bilinear action of pseudofermions and conduction 0.1 0.2 0.3 0.4 0.5 0 0.01 0.02 0.03 0.04 J T electrons that are only coupled through ﬂuctuating back- K 0 5 10 15 0 1 2 0 1 2 3 ground auxiliary ﬁelds, (c) 1 J =0.1 J =0.26 J =0.4 2 2 2 K K K JK |Vi,n| JH |χn| S = β ( + ) − β λi + S1, (2) y

2 4 1 i,n i=1 0

X X V where Vi,n and χn are the auxiliary fields in Matsub- ara frequency (iωn) representing the Kondo hybridiza- -1 -1 0 1 -1 0 1 -1 0 1 tion and intersite magnetic correlation, respectively. x x x V1 V1 V1 S1 = nmσ Ψnσ† (Onm − i ωnδnm)Ψmσ, with Ψnσ = [ck , ··· ,ck ,f ,f ]T and N being the number 1σn P N0 σn 1σn 2σn 0 of k points in the Brillouin zone of conduction electrons. FIG. 1: (a) The mean-field phase diagram. (b) Comparison The matrix Onm is of the mean-field |V1| and that derived from the peak position 1 1 of p(V1,V2,χ) and p(V1,V2) as a function of temperature for ǫk1 ··· 0 V1,n m V2,n m . . .− .− JK = 0.34 (dashed line in (a)). (c) Evolution of the probabilis-  ......  x y . . . . . tic distribution p(V1) on the complex plane of V1 =(V1 ,V1 ) N0 N0 Onm =  0 ··· ǫk V V  , at T = 0.001 for three different values of JK in the magnetic N0 1,n m 2,n m  1 N − −   0 JH χn−m  regime, near the quantum critical point, and in the Kondo  V 1,m n ··· V 1,m n λ1 2  − − regime, respectively.  1 N0 JH χm−n   V 2,m n ··· V 2,m n 2 λ2   − −  · i Ri kj j JK e Vi,n−m where the auxiliary fields V and χ are random vari- where Vi,n m = √ . λi is the Lagrange mul- i − 2 N0 2 tiplier for the constraints ni = 1 and takes a real value ables satisfying p0(Vi) ∼ exp(−βJK|Vi| /2) and p0(χ) ∼ f 2 after a Wick rotation [62]. exp(−βJH|χ| /2), respectively. The model can then The above action is, however, generally unsolvable. be studied using Monte Carlo simulations [69–73]. Al- To proceed, we propose a static approximation assuming ternately, one may first eliminate the conduction electron part in O and obtain ln det(O − i ω ) = Vi,n m = Viδnm, χn m = χδnm such that Onm = Oδnm. n n − − This ignores the temporal fluctuations of the auxiliary n ln det(An − i ωn)+ S0 with P fields but takes full account of their spatial fluctuations P λ − ∆ (i ω ) JHχ − ∆ (i ω ) and probabilistic distribution beyond the mean-field ap- 1 11 n 2 12 n An = JHχ . (5) proximation [63–68]. To see how it works, we first inte- 2 − ∆21(i ωn) λ2 − ∆22(i ωn)

− · − grate out all fermionic degrees of freedom and obtain an 2 i k (Ri Rj) JKV iVj e Here ∆ (i ω ) = k . S is a effective action only of the auxiliary fields, ij n 4N0 i ωn+ǫk 0 constant from conduction electrons− and can be safely 2 2 2 P JK |Vi| JH |χ| dropped. The result is also real because An† = A n. S = β ( + −λ )−2 ln det(O−iω ). − eff 2 4 i n The hybridization physics can then be studied with Xi=1 Xn 1 (3) the probabilistic distribution p(Vi,χ) = Z− exp(−Seff), The summation over Matsubara frequency can be evalu- where Z is the partition function serving as the normal- βξl ization factor. Vi and χ are both complex numbers and, ated using n ln det(O − i ωn)= l ln(1 + e− ), where due to the high dimensionality of all variables, one may ξl denote theP eigenvalues of O andP are always real because O is Hermitian. This action can also be derived from an use the Monte Carlo and Metropolis algorithm for im- i portance sampling. Compared to the perturbation ex- effective Hamiltonian (with the constraints nf = 1): pansion or mean-field theory, all spatial fluctuations of JK H = ǫkc† ck + V c† f + H.c. the hybridization fields are included. Below we will fix eff kσ σ 2 i iσ iσ R R R Xkσ Xi,σ JH =0.1 and the impurity distance | | = | 2 − 1| =1 unless noted. For simplicity, we also set the half conduc- JH + χf † f + H.c. , (4) tion bandwidth to unity and use ǫk = −(cos k +cos k )/2 2 1σ 2σ x y Xσ with the particle-hole symmetry. The Lagrangian multi- 3 pliers are then approximated by their saddle-point value (a) 0.2 0.6 J =0.20 0.6 λ = 0 [74]. K T=0.001 i J =0.25 To better understand the fluctuation effect, we present K T=0.002 0.4 |) J =0.30

2 K 0.4 in Fig. 1(a) the mean-ﬁeld phase diagram with three dis- T=0.003 0.2 J K=0.35 |; |V 0.1 T=0.004 1 tinct regions. For small JK, the mean-ﬁeld hybridization J =0.40 0 K 0.001 0.004

I(|V 0.2 is zero (|V | = 0) and there is an artificial second order T phase transition from |χ| = 0 to |χ| 6= 0 at T = JH/4, below which the two impurity spins form a singlet due 0 0 0 0.01 0.02 0.03 0.1 0.2 0.3 0.4 0.5 to the Heisenberg interaction. With increasing JK, the T J mean-field theory predicts an artificial first-order transi- K tion to the Kondo phase (|V | 6= 0) [34, 35]. By contrast, (b) 0 10 20 30 0 10 20 30 40 50 0.8 NRG and CFT analyses suggest a “Varma-Jones” unsta- J K=0.1 J K=0.255 J K=0.4 ble fixed point at zero temperature [57–61] and argue that | fluctuations will suppress the artificial first-order transi- 2 0.4 tion and turn it into a crossover. |V This is indeed the case in our method. The mean-field 0 solution is equivalent to the saddle point approximation 0 0.4 0.8 0 0.4 0.8 0 0.4 0.8 for the effective action Seff. Figure 1(b) compares the |V 1| |V 1| |V 1| mean-field |V1| and that of maximal p(V1, V2,χ) at a cho- sen JK (the dashed line in Fig. 1(a)). As expected, the two agree well with each other. However, if we first in- FIG. 2: (a) The amplitude mutual information I(|V1|; |V2|) tegrate out χ and estimate |V1| from the maximum of as a function of the temperature and Kondo coupling. The the joint distribution p(V1, V2), the transition tempera- left inset shows a schematic diagram of the three-layer neural network and the right inset shows logarithmic divergence of ture will be greatly suppressed. Since p(V1, V2) includes the maximal mutual information at low temperatures. The the probabilistic distribution of χ, the difference reflects black solid line is a guide to the eye. (b) Comparison of the the effect of fluctuating magnetic correlations between distribution p(|V1|, |V2|) at three different JK at T = 0.001. impurities in reducing Kondo screening by conduction electrons. We will show later that once the statistical

fluctuations of Vi are also included, the first-order tran- true mutual information, I(X; Y ) ≥ Iθ(X; Y ), following sition does turn into a crossover in physical properties. the Donsker-Varadhan representation [78]. As long as Here just for illustration, we plot in Fig. 1(c) the dis- the parameter space Θ is large enough, the inequality is x y tribution p(V1) of the complex field V1 = (V1 , V1 ) after tight and becomes a good approximation. The neural integrating out V2 and χ from p(V1, V2,χ). With increas- estimator has been successfully applied to both thermal ing JK, the region of large p(V1) is seen to first expand and athermal systems [79]. We implement it here with from a spot around the origin and gradually develop into a three-layer neural network using Tensorflow and the a ring, showing a “Mexican hat” potential for the hy- Adam optimizer [80, 81]. bridization in the Kondo phase. We first integrate out the phase θi of the complex i θi To extract useful information on hybridization fluctu- Vi = |Vi| e and study the mutual information of the ations, we define the mutual information of two random amplitude |Vi|. The results are plotted in Fig. 2(a) as a variables X and Y with the joint probability p(x, y), function of the temperature and Kondo interaction, re- p(x, y) spectively. For small JK, the amplitude mutual informa- I(X; Y )= p(x, y) log dxdy, (6) tion I(|V |; |V |) is small but grows with lowering temper- p(x)p(y) 1 2 ZZ ature due to increasing intersite magnetic correlations; where p(x) = p(x, y)dy and p(y)= p(x, y)dx are the while for large JK, a peak appears at finite temperature, marginal probabilities.R Its calculationR has historically indicating the weakening of intersite magnetic correla- been challenging because we typically only have samples tions due to Kondo screening at low temperatures. Re- rather than the underlying distribution [75]. A straight- markably, the amplitude mutual information varies non- forward approach is to partition the samples into bins of monotonically with JK and exhibits a maximum whose finite size but the results are very sensitive to the bin sizes height increases rapidly and diverges logarithmically with [76]. Here we use a neural estimator recently proposed lowering temperature, manifesting the quantum critical for the mutual information [77]: behavior above the “Varma-Jones” fixed point. To get an intuitive picture, we compare in Fig. 2(b) the distri- E E fθ Iθ(X; Y ) = sup p(x,y)[fθ] − log p(x)p(y)[e ] , (7) bution p(|V1|, |V2|) at three different JK. We see that the θ Θ ∈ values of |Vi| are scattered along the diagonal direction where fθ is a function parametrized by neural networks near the critical JK. Hence the divergence comes from with the parameters θ ∈ Θ. Iθ gives a lower bound for the strong cooperative fluctuations of the hybridization fields 4

(a) 0.6 0.8 (a) 20 (b)0.16 T=0.001 0.1

JK=0.20 |R|=1 JK=0.3 0.6 T=0.002 16 ρ |R|=1 0.12 SAF 0.4 JK=0.25 ) T=0.003 MF 2

) SAF(i=1) θ

; J =0.30 12

1 K 0 0.4 T=0.004 JK=0.2 θ MF(i=1) 0.001 0.01 0.04 ρ ω=0 0.08 I( J =0.35 ( J =0.25 T K c K 0.2 SAF(i=3) A 8 J =0.3 JK=0.40 0.2 K MF(i=3) 0.04 4 0 0 0 0.01 0.02 0.03 0.1 0.2 0.3 0.4 0.5 T JK 0 0 4 6 8 2 4 6 8 2 4 0.1 0.2 0.3 0.4 0.5 0.001 0.01 0 0.5 1.0 #!# #!" !# !" 0 0.5 1.0 1.5 (b) JK T 2.0 (c)0.8 (d)0.16 JK=0.1 JK=0.255 JK=0.4 |R|=1

π 0.6 0.12 / |R|=1 |R|=8 2 1.0 θ

) |R|=8 2 300 x=0.094

θ 1.2 ; ρ 1 0.4 0.08 0.29 θ 0.8 200 I(

0 χ 0.50 0 1.0 2.0 0 1.0 2.0 0 1.0 2.0 0.4 cm/Ce)

µΩ 0.73 ( θ /π θ /π θ /π 0.2 0.04 ρ 100 1 1 1 0 0.90 0.1 0.2 0.3 0.4 0.5 1.0 JK 0 0.01 0.1 1 10 100 T(K) 0 0 0.1 0.2 0.3 0.4 0.5 2 3 4 5 6 2 3 4 FIG. 3: (a) The phase mutual information I(θ1; θ2) as a func- 0.001 0.01 JK T tion of the temperature and Kondo coupling. (b) Comparison of the distribution p(θ1,θ2) for different JK at T = 0.001. FIG. 4: (a) Comparison of the our results (SAF) and the mean-field results (MF) for the electron spectrum at (i = 1) on two impurities. Note that the maximal distribution of and away from (i = 3) the impurity sites at T = 0.001, show- |Vi| is not around the origin for small JK in polar coordi- ing strong suppression (pseudogap) due to hybridization fluc- nate, showing the presence of hybridization fluctuations tuations even in the magnetic regime. (b) Temperature evolution of the resistivity ρ at three different JK compared with even deep inside the magnetic phase. the MF result (inset). (c) Comparison of the phase mutual The mutual information of the phases is presented in information for different impurity distances |R| =1 and 8 at Fig. 3(a). Since the effective action of the two-impurity T = 0.001. The inset shows the mean-field order parameter i φi model is invariant under the transformation Vi → Vi e , χ. (d) Resistivity as a function of temperature for JK = 0.4, χ → χ e i(φ1 φ2), we fix the gauge such that χ is real showing a significant drop due to the intersite phase coher- − − R and nonnegative and study the phase mutual informa- ence for | | = 1 compared to that of the single impurity limit (|R| = 8). The inset reproduces the experimental resistivity tion I(θ1; θ2) of the hybridization fields. We see that it of Ce1−xLaxCu6 for various x [87]. is nearly zero at high temperatures where the intersite magnetic correlations are negligible. But unlike that of the amplitude, here it always grows with decreasing tem- cal density of states (LDOS), A (i,ω = 0) ≈− β G (i, τ = perature and, in the Kondo regime, varies only slightly c π c β/2), where G (i, τ) = −hT [c (τ)c†(0)]i is the Green’s with JK as shown in Fig. 3(a). Its large value reflects the c τ i i establishment of cooperative phase fluctuations between function of conduction electrons. We find a strong sup- two Kondo impurities. To see this more clearly, we plot in pression (pseudogap) at the impurity sites in the Kondo regime, in agreement with that observed in the scan- Fig. 3(b) the probabilistic distribution p(θ1,θ2) for three ning tunnelling experiments [82]. Interestingly, the sup- different values of JK. For small JK, the distribution is pression already starts in the magnetic regime owing to featureless; while for large JK, it is peaked along the diagonal direction, indicating that the two are “locked” with the hybridization fluctuations. There appears no abrupt each other. The mutual information therefore provides change across the critical point, in contrast to the mean- a useful tool to measure the intersite phase coherence of field expectation of a first-order phase transition. Fig- the hybridization fields. Note that while the pattern of ure 4(b) plots the magnetic resistivity as a function πT 2 of temperature calculated using ρ ≈ − , where p(θ1,θ2) may shift with the gauge of χ, I(θ1; θ2) is gauge gxx(τ=β/2) invariant. It is also related to the Shannon entropy of the gxx(τ) = −hTτ jx(τ)jx(0)i is the current-current corre- phase difference through I(θ ; θ ) = const − S(θ) where i 1 2 lation function and jx(τ) = 2 lσ[cl†+x,σcl,σ − h.c.] is S(θ)= − p(θ) log p(θ)dθ with θ = θ1 − θ2. the current operator [83, 84]. WeP find logarithmic diver- We mayR further study physical properties under this gence near the critical point, suppression in the magnetic scheme taking partly account of the effect of magnetic regime and saturation in the Kondo regime, indicating and hybridization fluctuations. Figure 4(a) plots the lo- non-Fermi liquid, metallic, and Kondo-like behaviors, re- 5 spectively. Again, the smooth evolution across the mean- tum Entanglement to Topological Phases of Many-Body field phase boundary implies a crossover rather than a Systems, (Springer, New York, 2019). first-order phase transition once fluctuations are included [3] T. M. Cover and J. A. Thomas, Elements of Information as in our method. Theory, (Wiley, New York, 1991). [4] M. M. Wolf, F. Verstraete, M. B. Hastings, and J. I. The importance of phase coherence is more evidently Cirac, Area Laws in Quantum Systems: Mutual Infor- seen if we increase the distance between two impurities. mation and Correlations, Phys. Rev. Lett. 100, 070502 Figure 4(c) compares the phase mutual information for (2008). |R| = 1 and 8. For large distance |R| = 8, the mutual [5] H. Matsuda, K. Kudo, R. Nakamura, O. Yamakawa, and information is reduced to zero even in the Kondo regime. T. Murata, Mutual Information of Ising Systems, Int. J. 35 A large finite phase mutual information only appears Theor. Phys. , 839 (1996). for small |R|, indicating a fundamental difference of the [6] R. T. Wicks, S. C. Chapman, and R. O. Dendy, Mu- tual information as a tool for identifying phase transi- multi-impurity Kondo physics from the single-impurity tions in dynamical complex systems with limited data, case. For comparison, the mean-field order parameter Phys. Rev. E 75, 051125 (2007). χ is also shown in the inset and found to be zero for [7] J. Wilms, M. Troyer, and F. Verstraete, Mutual infor- |R| = 8 and finite for |R| = 1. The phase coherence is mation in classical spin models, J. Stat. Mech. (2011) therefore closely associated with the intersite magnetic P10011. correlations [85, 86]. As shown in Fig. 4(d), it causes a [8] J. Wilms, J. Vidal, F. Verstraete, and S. Dusuel, Finite- significant drop of the resistivity for |R| = 1, in qual- temperature mutual information in a simple phase transition, J. Stat. Mech. (2012) P01023. itative agreement with the experimental observation in [9] J. Iaconis, S. Inglis, A. B. Kallin, and R. G. Melko, De- Ce1 xLaxCu6 at large x [87]. − tecting classical phase transitions with Renyi mutual in- We conclude that the static auxiliary field approxima- formation, Phys. Rev. B 87, 195134 (2013). tion can take good account of some fluctuation effects [10] H. W. Lau and P. Grassberger, Information theoretic as- beyond the mean-field approximation. The employment pects of the two-dimensional Ising model, Phys. Rev. E 87 of the mutual information can effectively reduce the di- , 022128 (2013). mensionality of the data and allow one to extract key [11] H. Casini, I. S. Landea, and G. Torroba, The g-theorem and quantum information theory, J. High Energy Phys. information hidden in the complicated probabilistic dis- 10 (2016) 140. tribution functions. It is thus a useful tool to probe the [12] I. Hagymási, J. Sólyom, and O.¨ Legeza, Competition quantum phase transition and phase coherence in Kondo between Hund’s coupling and Kondo effect in a one- systems. Our method can be extended easily to other dimensional extended periodic Anderson model, Phys. Kondo models to investigate spatial correlations of the Rev. B 92, 035108 (2015). hybridization physics. In the dense Kondo lattice, we [13] F. P. Toldin, T. Sato, and F. F. Assaad, Mutual informa- 99 expect that it will also suppress the artificial phase tran- tion in heavy-fermion systems, Phys. Rev. B , 155158 (2019). sitions predicted in an effective Ginzburg-Landau-Wilson [14] J. S. Hofmann, F. F. Assaad, and T. Grover, Fraction- theory [41–43] and reveal the importance of cooperative alized Fermi liquid in a frustrated Kondo lattice model, hybridization fluctuations and intersite phase coherence, Phys. Rev. B 100, 035118 (2019). which is important for understanding the heavy Fermi [15] C. Walsh, P. Sémon, D. Poulin, G. Sordi, and A. M. liquid [54]. The method may also be used to reveal the S. Tremblay, Local Entanglement Entropy and Mutual snapshot or spatial modulation of hybridization configu- Information across the Mott Transition in the Two- 122 rations [88–90] and provide novel insight on the Kondo Dimensional Hubbard Model, Phys. Rev. Lett. , 067203 (2019). physics in multi-impurity or lattice systems. [16] C. Walsh, M. Charlebois, P. Sémon, G. Sordi, and A. This work was supported by the National Key Re- M. S. Tremblay, Information-theoretic measures of super- search and Development Program of MOST of China conductivity in a two-dimensional doped Mott insulator, (Grant No. 2017YFA0303103), the National Natural Sci- Proc. Natl. Acad. Sci. USA 118, e2104114118 (2021). ence Foundation of China (NSFC Grant No. 11974397, [17] M. A. Valdez, D. Jaschke, D. L. Vargas, and L. D. Carr, No. 11774401), the Strategic Priority Research Pro- Quantifying Complexity in Quantum Phase Transitions gram of the Chinese Academy of Sciences (Grant No. via Mutual Information Complex Networks, Phys. Rev. Lett. 119, 225301 (2017). XDB33010100), and the Youth Innovation Promotion [18] A. A. Bagrov, M. Danilov, S. Brener, M. Harland, A. Association of CAS. I. Lichtenstein, and M. I. Katsnelson, Detecting quan- ′ tum critical points in the t-t Fermi-Hubbard model via complex network theory, Sci. Rep. 10, 20470 (2020). [19] A. C. Hewson, The Kondo Problem to Heavy Fermions, (Cambridge University Press, Cambridge, England, ∗ [email protected] 1997). [1] P. W. Anderson, More Is Different, Science 177, 393 [20] P. Coleman, Introduction to Many-body Physics, (Cam- (1972). bridge University Press, Cambridge, England, 2015). [2] B. Zeng, X. Chen, D.-L. Zhou, and X.-G. Wen, Quan- [21] Y.-F. Yang, D. Pines, and G. Lonzarich, Quantum crit- tum Information Meets Quantum Matter: From Quan- 6

ical scaling and fluctuations in Kondo lattice materials, [43] K. Ohara and K. Hanzawa, Nambu-Goldstone and Higgs Proc. Natl. Acad. Sci. USA 114, 6250 (2017). Modes in the Kondo Insulator, J. Phys. Soc. Jpn. 83, [22] S. Doniach, The Kondo lattice and weak antiferromag- 104604 (2014). netism, Physica B+C, 91, 231 (1977). [44] R. M. Fye and J. E. Hirsch, Quantum Monte Carlo study [23] D. M. Newns and N. Read, Mean-field theory of inter- of the two-impurity Kondo Hamiltonian, Phys. Rev. B mediate valence/heavy fermion systems, Adv. Phys. 36, 40, 4780 (1989). 799 (1987). [45] R. M. Fye, “Anomalous Fixed Point Behavior” of Two [24] C. Sire, C. M. Varma, and H. R. Krishnamurthy, The- Kondo Impurities: A Reexamination, Phys. Rev. Lett. ory of the non-Fermi-liquid transition point in the two- 72, 916 (1994). impurity Kondo model, Phys. Rev. B 48, 13833 (1993). [46] K. Hallberg and R. Egger, Two-impurity Kondo problem [25] G. M. Zhang and L. Yu, Kondo singlet state coexist- for correlated electrons, Phys. Rev. B 55, R8646 (1997). ing with antiferromagnetic long-range order: A possible [47] J. B. Silva, W. L. C. Lima, W. C. Oliveira, J. L. N. ground state for Kondo insulators, Phys. Rev. B 62, 76 Mello, L. N. Oliveira, and J. W. Wilkins, Particle-Hole (2000). Asymmetry in the Two-Impurity Kondo Model, Phys. [26] B. Coqblin, C. Lacroix, M. A. Gusmaão, and J. R. Rev. Lett. 76, 275 (1996). Iglesias, Band-filling effects on Kondo-lattice properties, [48] V. L. Campo and L. N. Oliveira, Thermodynamics for Phys. Rev. B 67, 064417 (2003). the two-impurity Kondo model, Phys. Rev. B 70, 153401 [27] T. Senthil, S. Sachdev, and M. Vojta, Fractionalized (2004). Fermi Liquids, Phys. Rev. Lett. 90, 216403 (2003). [49] R. Bulla, T. A. Costi, and T. Pruschke, Numerical renor- [28] T. Senthil, M. Vojta, and S. Sachdev, Weak magnetism malization group method for quantum impurity systems, and non-Fermi liquids near heavy-fermion critical points, Rev. Mod. Phys. 80, 395 (2008). Phys. Rev. B 69, 035111 (2004). [50] L. Zhu and J.-X. Zhu, Coherence scale of coupled Ander- [29] J. Rech, P. Coleman, G. Zarand, and O. Parcollet, son impurities, Phys. Rev. B 83, 195103 (2011). Schwinger Boson Approach to the Fully Screened Kondo [51] A. Bayat, S. Bose, P. Sodano, and H. Johannesson, En- Model, Phys. Rev. Lett. 96, 016601 (2006). tanglement Probe of Two-Impurity Kondo Physics in a [30] A. O’Bannon, I. Papadimitriou, and J. Probst, A holo- Spin Chain, Phys. Rev. Lett. 109, 066403 (2012). graphic two-impurity Kondo model, J. High Energy [52] A. Bayat, H. Johannesson, S. Bose, and P. Sodano, An Phys. 01 (2016) 103. order parameter for impurity systems at quantum criti- [31] G. Zhang, J. S. V. Dyke, and R. Flint, Cubic hastatic cality, Nat. Commun. 5, 3784 (2014). order in the two-channel Kondo-Heisenberg model, Phys. [53] D. Hu, J.-J. Dong, and Y.-F. Yang, Hybridization fluctu- Rev. B 98, 235143 (2018). ations in the half-filled periodic Anderson model, Phys. [32] A. Wugalter, Y. Komijani, and P. Coleman, Large-N ap- Rev. B 100, 195133 (2019). proach to the two-channel Kondo lattice, Phys. Rev. B [54] M. Raczkowski and F. F. Assaad, Emergent Coherent 101, 075133 (2020). Lattice Behavior in Kondo Nanosystems, Phys. Rev. [33] P. Coleman, Mixed valence as an almost broken symme- Lett. 122, 097203 (2019). try, Phys. Rev. B 35, 5072 (1987). [55] A. Cai, H. Hu, K. Ingersent, S. Paschen, and Q. Si, Dy- [34] A. Jones, B. G. Kotliar, and A. J. Millis, Mean-field anal- namical Kondo effect and Kondo destruction in effective ysis of two antiferromagnetically coupled Anderson im- models for quantum-critical heavy fermion metals, arXiv: purities, Phys. Rev. B 39, 3415 (1989). 1904.11471 (2019). [35] B. A. Jones, Pair correlation effects in heavy fermions, [56] H. Hu, A. Cai, and Q. Si, Quantum Criticality and Physica B 171, 53 (1991). Dynamical Kondo Effect in an SU(2) Anderson Lattice [36] Y. P. Liu, Y. J. Zheng, J. J. Dong, H. Lee, Z. X. Wei, W. Model, arXiv: 2004.04679 (2020). L. Zhang, C. Y. Chen, H. Q. Yuan, Y.-F. Yang, and J. Qi, [57] B. A. Jones and C. M. Varma, Study of Two Magnetic Hybridization Dynamics in CeCoIn5 Revealed by Ultra- Impurities in a Fermi Gas, Phys. Rev. Lett. 58, 843 fast Optical Spectroscopy, Phys. Rev. Lett. 124, 057404 (1987). (2020). [58] B. A. Jones, C. M. Varma, and J. W. Wilkins, Low- [37] A. Auerbach and K. Levin, Kondo Bosons and the Kondo Temperature Properties of the Two-Impurity Kondo Lattice: Microscopic Basis for the Heavy Fermi Liquid, Hamiltonian, Phys. Rev. Lett. 61, 125 (1988). Phys. Rev. Lett. 57, 877 (1986). [59] B. A. Jones and C. M. Varma, Critical point in the solu- [38] P. Coleman, J. B. Marston, and A. J. Schofield, Trans- tion of the two magnetic impurity problem, Phys. Rev. port anomalies in a simplified model for a heavy-electron B 40, 324 (1989). quantum critical point, Phys. Rev. B 72, 245111 (2005). [60] I. Affleck and A. W. W. Ludwig, Exact Critical Theory [39] I. Paul, C. Pépin, and M. R. Norman, Kondo Breakdown of the Two-Impurity Kondo Model, Phys. Rev. Lett. 68, and Hybridization Fluctuations in the Kondo-Heisenberg 1046 (1992). Lattice, Phys. Rev. Lett. 98, 026402 (2007). [61] I. Affleck, A. W. W. Ludwig, and B. A. Jones, Conformal- [40] I. Paul, C. Pépin, and M. R. Norman, Multiscale fluctu- field-theory approach to the two-impurity Kondo prob- ations near a Kondo breakdown quantum critical point, lem: Comparison with numerical renormalization-group Phys. Rev. B 78, 035109 (2008). results, Phys. Rev. B 52, 9528 (1995). [41] K. Hanzawa and K. Ohara, Asymptotic local moment [62] Y. Zhou, K. Kanoda, and T.-K. Ng, Quantum spin liquid formation in the Kondo lattice, Phys. Rev. B 76, 184407 states, Rev. Mod. Phys. 89, 025003 (2017). (2007). [63] A. Mukherjee, N. D. Patel, S. Dong, S. Johnston, A. [42] K. Ohara and K. Hanzawa, A Generalized Ginzburg- Moreo, and E. Dagotto, Testing the Monte Carlo mean Landau-Wilson Theory of the Kondo Lattice, J. Phys. field approximation in the one-band Hubbard model, Soc. Jpn. 82, 104713 (2013). Phys. Rev. B 90, 205133 (2014). 7

[64] S. Pradhan and G. V. Pai, Holstein-Hubbard model at Conference on Machine Learning, PMLR 80, 531 (2018). half filling: A static auxiliary field study, Phys. Rev. B [78] M. D. Donsker and S. R. S. Varadhan, Asymptotic Eval- 92, 165124 (2015). uation of Certain Markov Process Expectations for Large [65] M. Karmakar and P. Majumdar, Population-imbalanced Time, I, Commun. Pure Appl. Math. 28, 1 (1975). lattice fermions near the BCS-BEC crossover: Thermal [79] A. Nir, E. Sela, R. Beck, and Y. B. Sinai, Machine- physics of the breached pair and Fulde-Ferrell-Larkin- learning iterative calculation of entropy for physical sys- Ovchinnikov phases, Phys. Rev. A 93, 053609 (2016). tems, Proc. Natl. Acad. Sci. USA 117, 30234 (2020). [66] N. D. Patel, A. Mukherjee, N. Kaushal, A. Moreo, and [80] M. Abadi, P. Barham, J. Chen, Z. Chen, A. Davis, J. E. Dagotto, Non-Fermi Liquid Behavior and Continu- Dean, M. Devin, S. Ghemawat, G. Irving, M. Isard, M. ously Tunable Resistivity Exponents in the Anderson- Kudlur, J. Levenberg, R. Monga, S. Moore, D. G. Mur- Hubbard Model at Finite Temperature, Phys. Rev. Lett. ray, B. Steiner, P. Tucker, V. Vasudevan, P. Warden, M. 119, 086601 (2017). Wicke, Y. Yu, and X. Zheng, TensorFlow: A System for [67] G. Jana and A. Mukherjee, Frustration effects at finite Large-Scale Machine Learning, In 12th USENIX Sympo- temperature in the half filled Hubbard model, J. Phys. sium on Operating Systems Design and Implementation Condens. Matter 32, 365602 (2020). (OSDI 16), (Savannah, GA), 265 (2016). [68] M. Karmakar, Pauli limited d-wave superconductors: [81] D. P. Kingma and J. L. Ba, Adam: A Method For quantum breached pair phase and thermal transitions, Stochastic Optimization, In International Conference on J. Phys. Condens. Matter 32, 405604 (2020). Learning Representations (ICLR), (2015). [69] P. Farkaˇsovský, H. Cenˇcariková,ˇ and S. Mataˇs, Numeri- [82] A. Spinelli, M. Gerrits, R. Toskovic, B. Bryant, M. cal study of magnetization processes in rare-earth tetra- Ternes, and A. F. Otte, Exploring the phase diagram borides, Phys. Rev. B 82, 054409 (2010). of the two-impurity Kondo problem, Nat. Commun. 6, [70] Y. Kato, I. Martin, and C. D. Batista, Stability of 10046 (2015). the Spontaneous Quantum Hall State in the Triangu- [83] L.-Y. Wei and Y.-F. Yang, Doping-induced perturbation lar Kondo-Lattice Model, Phys. Rev. Lett. 105, 266405 and percolation in the two-dimensional Anderson lattice, (2010). Sci. Rep. 7, 46089 (2017). [71] H. Ishizuka and Y. Motome, Partial Disorder in an Ising- [84] E. W. Huang, R. Sheppard, B. Moritz, and T. P. Dev- Spin Kondo Lattice Model on a Triangular Lattice, Phys. ereaux, Strange metallicity in the doped Hubbard model, Rev. Lett. 108, 257205 (2012). Science 366, 987 (2019). [72] W.-W. Yang, J. Zhao, H.-G. Luo, and Y. Zhong, Exactly [85] S. Nakatsuji, S. Yeo, L. Balicas, Z. Fisk, P. Schlottmann, solvable Kondo lattice model in the anisotropic limit, P. G. Pagliuso, N. O. Moreno, J. L. Sarrao, and J. D. Phys. Rev. B 100, 045148 (2019). Thompson, Intersite Coupling Effects in a Kondo Lattice, [73] M. M. Maska and N. Trivedi, Temperature-driven BCS- Phys. Rev. Lett. 89, 106402 (2002). BEC crossover and Cooper-paired metallic phase in cou- [86] Y.-F. Yang, Z. Fisk, H.-O. Lee, J. D. Thompson, and D. pled boson-fermion systems, Phys. Rev. B 102, 144506 Pines, Scaling the Kondo lattice, Nature 454, 611 (2008). (2020). [87] Y. Onuki¯ and T. Komatsaubara, Heavy Fermion State [74] S. Saremi, P. A. Lee, and T. Senthil, Unifying Kondo In CeCu6, J. Magn. Magn. Mater. 63-64, 281 (1987). coherence and antiferromagnetic ordering in the honey- [88] J.-X. Zhu, I. Martin, and A. R. Bishop, Kondo Stripes comb lattice, Phys. Rev. B 83, 125120 (2011). in an Anderson-Heisenberg Model of Heavy Fermion Sys- [75] B. Poole, S. Ozair, A. V. D. Oord, A. A. Alemi, and G. tems, Phys. Rev. Lett. 100, 236403 (2008). Tucker, On Variational Bounds of Mutual Information, [89] C. Pépin, M. R. Norman, and A. Ferraz, Modulated Spin Proceedings of the 36th International Conference on Ma- Liquid: A New Paradigm for URu2Si2, Phys. Rev. Lett. chine Learning, PMLR 97, 5171 (2019). 106, 106601 (2011). [76] A. Kraskov, H. Stögbauer, and P. Grassberger, Estimat- [90] A. Thomson and S. Sachdev, Fractionalized Fermi liquid ing mutual information, Phys. Rev. E 69, 066138 (2004). on the surface of a topological Kondo insulator, Phys. [77] M. I. Belghazi, A. Baratin, S. Rajeswar, S. Ozair, Y. Ben- Rev. B 93, 125103 (2016). gio, A. Courville, and R. D. Hjelm, Mutual Information Neural Estimation, Proceedings of the 35th International