# Magnetic Quantum Phase Diagram of Magnetic Impurities in 2 Dimensional Disordered Electron Systems

## Abstract

The quantum phase diagram of disordered electron systems as function of the concentration of magnetic impurities and the local exchange coupling is studied in the dilute limit. We take into account the Anderson localisation of the electrons by a nonperturbative numerical treatment of the disorder potential. The competition between RKKY interaction and the Kondo effect, as governed by the temperature scale , is known to gives rise to a rich magnetic quantum phase diagram, the Doniach diagram. Our numerical calculations show that in a disordered system both the Kondo temperature and are widely distributed. Accordingly, also their ratio, is widely distributed as shown in Fig. 1 (a). However, we find a sharp cutoff of that distribution, which allows us to define a critical density of magnetic impurities below which Kondo screening wins at all sites of the system above a critical coupling , forming the Kondo phase [see Fig. 1 (b)]. As disorder is increased, increases and a spin coupled phase is found to grow at the expense of the Kondo phase. From these distribution functions we derive the magnetic susceptibility which show anomalous power law behavior. In the Kondo phase that power is determined by the wide distribution of the Kondo temperature, while in the spin coupled phase it is governed by the distribution of . At low densities and small we identify a paramagnetic phase. We also report results on a honeycomb lattice, graphene, where we find that the spin coupled phase is more stable against Kondo screening, but is more easily destroyed by disorder into a PM phase.

## I Introduction

Phenomena which emerge from the interplay of strong correlations and disorder remain a challenge for condensed matter theory. Spin correlations and disorder effects are however relevant for a wide range of materials, including doped semiconductors like Si:P close to metal-insulator transitions,v. Löhneysen (2000) and heavy Fermion systems, materials with 4f or 5f atoms, notably Ce, Yb, or U.Löhneysen et al. (2007) Many of these materials show a remarkable magnetic quantum phase transition which can be understood by the competition between indirect exchange interaction, the Ruderman-Kittel-Kasuya-Yoshida (RKKY) interaction between localised magnetic momentsRuderman and Kittel (1954); Kasuya (1956); Yosida (1957) and their Kondo screening.

Thereby, one finds a suppression of long range magnetic order when exchange coupling is increased and Kondo screening succeeds. This results in a typical quantum phase diagram with a quantum critical point where the of the magnetic phase is vanishing, the Doniach diagram.Doniach (1977)

Recently, controlled studies of magnetic adatoms on the surface of metals,Zhou et al. (2010) on graphene,Chen et al. (2007) and on the conducting surface of topological insulatorsHsieh et al. (2009); Chen et al. (2010); Scholz et al. (2012); Seibel et al. (2012) with surface sensitive experimental methods like spin resolved STM and ARPES became possible. This demands a theoretical study of the Doniach diagram for magnetically doped disordered electron systems (DES), in particular 2D systems.

In any material there is some degree of disorder. In doped semiconductors it arises from the random positioning of the dopants themselves, in heavy Fermion metals and in 2D metals it may arise from structural defects or impurities. Disorder is known to cause Anderson localisation, which therefore has to be taken into account when deriving the Doniach diagram of disordered electrons systems. Moreover, as noted already early,Anderson (1977) the physics of random systems is fully described by probability distributions, not just averages. This must be particularly true for systems with random local magnetic impurities (MIs),Mott, N. F. (1976) since the magnetic impurities are exposed to the local density of states of the conduction electrons, which is widely distributed itself. In fact, it has been noticed that a wide distribution of the Kondo temperature of MIs in disordered host metals gives rise to non-Fermi liquid behavior, such as the low temperature power-law divergence of the magnetic susceptibility.Mott, N. F. (1976); Miranda and Dobrosavljevic (2005); Miranda et al. (1996); Bhatt and Fisher (1992); Langenfeld and WÃ¶lfle (1995); Castro Neto and Jones (2000); Cornaglia et al. (2006); Kettemann and Mucciolo (2007); Tran and Kim (2010) Nonmagnetic disorder quenches the Kondo screening of MIs due to Anderson-localisation and the formation of local pseudogaps at the Fermi energy,Kettemann and Mucciolo (2007); Zhuravlev et al. (2007); Kettemann and Raikh (2003); Kettemann et al. (2012) resulting in bimodal distributions of and a finite concentration of free, paramagnetic moments (PMs). However, in these studies the RKKY interaction between different MIs has not yet been taken into account. is mediated by the conduction electrons, and aligns the spins of the MIs ferromagnetically or antiferromagnetically, depending on their distance . This is a long-ranged interaction, with a power law decay , where d is the dimension, and its typical value is not changed by weak disorder.Lerner (1993); Bulaevskii and Panyukov (1986); Bergmann (1987); Lee et al. (2012a) However, its amplitude has a wide log-normal distribution in disordered metals.Lerner (1993); Lee et al. (2012b) In this article we therefore intend to study the competition between RKKY interaction and the Kondo effect in disordered electron systems.

In the next section we introduce the model, and provide the equations for the Kondo temperature and the RKKY coupling. In section III, we derive numerically the distribution function of , and compare it with an analytical result, based on a perturbative expansion of the nonlinear sigma model. We derive next numerically the distribution function of finding excellent agreement with approximate analytical results which were obtained, taking into account the multifracatlity and power law correlations of wave functions. In section IV we present the main results, the distribution function of the Ratio of TK and RKKY Interaction, for various distances between magnetic impurities R. From that we show how to derive the zero temperature magnetic quantum phase diagram as function of magnetic impurity density and exchange coupling, for 2D disordered electronic systems. At low densities and small we identify a paramagnetic phase. For graphene we find that the spin coupled phase is more stable against Kondo screening, but is more easily destroyed by disorder into a PM phase. In section V we derive from the distribution functions the magnetic susceptibility as function of temperature, which show anomalous power law behavior. In the Kondo phase that power is found to be determined by the wide distribution of the Kondo temperature, while at small exchange coupling there we identify spin coupled phase where the magnetic susceptibility is governed by the distribution of . In the final section we conclude and discuss the relevance and limitations of our results.

## Ii Model

In order to obtain the Doniach diagram of random electron systems we extend the approach of DoniachDoniach (1977) by calculating the distribution functions of and and their ratio. Thus, in our approach we try to draw conclusions on the quantum phase diagram of an electron system with a finite density of magnetic impurities, by considering the Kondo temperature of single impurities and the RKKY coupling of pairs of magnetic moments.

We start from a microscopic description of the MIs, the Anderson impurity model coupled to a non-interacting disordered electronic Hamiltonian with on-site disorder. Then, we map it with the Schrieffer-Wolff transformation on a model of Kondo impurity spins coupled to the disordered host electron spins by the local coupling .Kettemann et al. (2012) We consider the single-impurity and the coupling between a pair of spins. For the numerical calculations we employ the single-band Anderson tight-binding model on a square lattice of size and lattice spacing ,

(1) |

where is the hopping energy between nearest neighbours , is the on-site disorder potential distributed in the interval . , where is the Fermi energy measured from the band edge, in 2D . We use periodic boundary conditions.

In the dilute limit, one can calculate the of a single magnetic impurity at position from the Nagaoka-Suhl one-loop equation,Nagaoka (1965); Suhl (1965)

(2) |

with band width . is the local density of states (LDOS). The RKKY coupling between two MIs located at positions , is in the zero temperature limit () given byRoche and Mayou (1999); Lee et al. (2012a)

(3) |

where , and is the magnitude of the MI spin.

## Iii Distribution Functions

Using the Kernel Polynomial method (KPM),Roche and Mayou (1999); Weiße et al. (2006) one can evaluate the matrix elements of the density matrix Lee et al. (2012a); Weiße et al. (2006); Mucciolo (2010) with a polynomial expansion of order . Here, we increase the cutoff degree linearly with the linear system size based on our analysis for the convergence of RKKY interaction with respect in Ref. Lee et al., 2012b. It has been also carefully discussed in Ref. Wenk and Bouzerar, 2013 that the choice of , not , gives proper DOS and LDOS results avoiding finite size effect.

Eq. (3) yields in a clean 2D system

in the asymptotic limit with effective electron mass , and Fermi wave vector . Ruderman and Kittel (1954) Its geometrical average is close to the clean limit for distances smaller than localisation length , and decays exponentially at larger distances,Sobota et al. (2007); Lee et al. (2012a) As shown in Figs. 2 a, b, the distribution of the absolute value of is well fitted by a log-normal,

where and the fitting gives for and , and width increasing with the disorder strength . This is qualitatively consistent with analytical results obtained at weak disorder,Lerner (1993) while analytical calculations at strong disorder have not been performed yet. This distribution width hardly depends on the distance . We used disorder configurations.

The distribution of is shown in Fig. 2 c, as obtained from the numerical solution of Eq. (2) for , . Since for every sample only one single site is taken to avoid a distortion of the distribution due to intersite correlations, we had to use a huge number of different random disorder configurations to get sufficient statistics. It has a strongly bimodal shape where the low - peak becomes more distinctive with larger disorder amplitude .Cornaglia et al. (2006); Kettemann and Mucciolo (2007, 2006) In Fig. 2 d we show these results for fixed disorder strength for various exchange couplings . Recently, an analytical derivation of the low -tail of was done, using the multifractal distribution and correlations of intensities.Kettemann et al. (2012) These correlations are in 2D logarithmic with an amplitude of order , where . For weak disorder, , it corresponds to a power law correlation with power The correlation energy is of the order of the elastic scattering rate . Thus, for ,Kettemann et al. (2012)

(4) |

where , the ratio of free PMs. Eq. (4) has a power law tail with power in good agreement with the numerical results, Fig. 2 d, for . For one findsKettemann et al. (2012)

(5) |

where . This expression is in agreement with the numerical results, see Fig. 2 d, using , and , fitting only and the prefactor. Thus, we confirm that the power law tail is governed by the multifractal correlation with power .

The quantum phase transition between the free paramagnetic moment phase (PM) and a Kondo screened phase can be studied by calculating the critical exchange coupling above which there is no more than one free magnetic moment in the sample volume .Zhuravlev et al. (2007) From the multifractality of the eigenfunction intensities it is found to be related to the power of the power law correlations in the 2D DES as and thus to increase in 2D linearly with disorder strength as,Kettemann et al. (2012)

(6) |

In Fig. 3 a, Eq. (6) is plotted together with numerical results as function of disorder strength . We find good agreement. There are only deviations at large disorder, , where the expansion breaks down. We plot as function of in Fig. 3 b, together with the density of states (DOS). We find that is increasing towards the band edge as in agreement with Eq. (6). Far outside of of the clean system it increases as due to the gap in the DOS.

## Iv Magnetic Phase Diagram at

In clean systems the critical density above which the MIs are coupled with each other can be obtained from the condition that .Doniach (1977) Thus, in 2D with and , , one finds .

In disordered systems, of an MI at a given site competes with the RKKY coupling to another MI at distance . Thus, the distribution function of the ratio of these two energy scales for a given disordered sample with density of MIs , where is the average distance between the MIs, is crucial to determine its magnetic state. The distribution of for and is shown for several distances in Fig. 1 (a) (, , and ). Likewise and , the distribution of has an exponentially wide width characterized by a small- tails and a sharp upper cutoff in as shown in Fig. 1 (a). As increasing the distance between the magnetic impurities the distribution is shifted to the left (smaller ), since the RKKY interaction decreases with . The sharp upper cutoff in allows us to define a critical density below which the Kondo effect dominates in the competition with RKKY interaction at all sites. is plotted in Fig. 1 (b) for various values of disorder strength . When the MI density exceeds , magnetic clusters start to form at some sites and the MIs may be coupled by . We see that this coupled moment phase (CM) expands at the expense of the Kondo phase with increasing . When is larger than localisation length the coupling is exponentially small and there is a paramagnetic phase (PM) below where MIs remain free up to exponentially small temperatures.

In graphene the pseudogap at the Dirac point quenches the Kondo effect below , independently on disorder amplitude . Thus, in graphene there is a larger parameter space where the MIs are coupled (CM) than in a normal 2DES, see Fig. 4. However, short range disorder localises the electrons, cutting off the RKKY-interaction and for there is a PM phase. Thus, the magnetic phase in graphene is more stable against Kondo screening but is more easily destroyed by disorder.

## V Doniach Phase Diagram of disordered 2DES and Graphene.

We find, that the Kondo phase splits at finite temperature into a Kondo Fermi-liquid (FL) phase, where all MIs are screened, and a Kondo Non-Fermi-liquid (NFL) phase, at , where some MIs remain unscreened and contribute to the magnetic susceptibility with an anomalous temperature dependence, given by,Kettemann et al. (2012)

(7) |

The temperature , plotted schematically in Fig. 5 (blue line), is given by the position of the low peak in the distribution , see Fig. 2 c. We note that may be distributed itself and may add a nonuniversal, material dependent contribution to the distribution of Langenfeld and WÃ¶lfle (1995) and .

For there is a succession of phases, starting with the RKKY phase where clusters are formed locally due to the widely distributed RKKY coupling. Anomalous power laws are observed when clusters are broken up successively as temperature is raised. From the log-normal distribution one obtains for the magnetic susceptibility,

(8) | |||||

where width increases with disorder strength . Accordingly, the excess specific heat is

(9) |

The detailed analysis of the quantum phase diagram at higher concentrations requires to go beyond our present analysis. One expects that at a spin-glass phase appears, where the magnetic susceptibility shows a peak at spin glass temperature as studied in Refs. Binder and Young, 1986; Coqblin et al., 2003. Above a critical density a phase with long range order may form below a critical temperature .Coqblin et al. (2003); Varma (1976); Magalhaes et al. (2006); Bouzerar et al. (2005)

## Vi Conclusions and Discussion

We conclude that it is the full distribution function of the ratio of the RKKY coupling and the Kondo temperature which determines the magnetic phase diagram of magnetic moments in disordered electron systems, especially at low concentrations. We identified a critical density of magnetic impurities below which Kondo wins at all positions in a disordered sample above a critical coupling , which increases with the disorder amplitude. As a result, the Kondo phase is diminished as the disorder is increased, favoring a phase where the MI spins are coupled. The magnetic susceptibility obeys an anomalous power law behavior, which crosses over as function of from the Kondo regime where that power is determined by the wide distribution of the Kondo temperature , to a spin coupled phase where it is governed by the log-normal distribution of . At low densities and small , we identify a paramagnetic phase. The distribution function of is expected to determine also the magnetic phase diagram of magnetically doped graphene and the surface of topological insulators with magnetic adatoms, see Fig. 4. This distribution function may also be crucial to explain the anomalous magnetic properties of doped semiconductors in the vicinity of metal-insulator transition,v. Löhneysen (2000) where we expect that is replaced by the universal value , with the universal multifractality parameter .

In this work we considered the distribution function of the Kondo temperature of single impurities and the RKKY coupling of pairs of magnetic moments and extracted information on the quantum phase diagram of systems with finite concentrations of MIs. While this approach has its limitations, for example at finite concentration the RKKY coupling can reduce as has been already found by Tsay and Klein in the 70s.Tsay and Klein (1973, 1975) However, they concluded that this reduction is minor. More importantly, later work revealed that the Kondo lattice of a finite density of magnetic moments, which is coupled to the conduction electrons, has a coherent low temperature heavy fermion phase, and a Kondo insulator phase at half filling of the magnetic moment levels. More recently, the Kondo lattice in 1 dimension was studied more rigorously (see the review by Tsunetsugu et. alTsunetsugu et al. (1997)), and it was shown that, at least in 1D, the groundstate of this system can not be understood by the mere extension of the single and two-magnetic impurity problem, where the physics is governed by the competition between these two energy scales, the Kondo temperature and the RKKY coupling. However, the higher temperature behavior was found to be still governed by the competition between these two energy scales. Therefore, we expect that the consideration of the reduced problem of two impurity spins, will give important information on the physics of disordered electron systems at finite concentration of magnetic moments, which becomes more meaningful the lower the density and the higher the temperature is. Going beyond the limitations of this approach, one will have to study the disordered Kondo lattice where a finite density of magnetic moments is coupled to the conduction electrons. For a clean Kondo lattice it is known that a coherent low temperature heavy fermion phase, and a Kondo insulator phase at half filling of the magnetic moment levels appears.Coleman (1984); Newns and Read (1987) It remains to see how these low temperature phases are modified by the presence of nonmagnetic disorder.

###### Acknowledgements.

We gratefully acknowledge useful discussions with Georges Bouzerar, Ki-Seok Kim, Eduardo Mucciolo and Keith Slevin, as well as the support by the BK21 Plus funded by the Ministry of Education, Korea (10Z20130000023).### References

- H. v. Löhneysen, Adv. in Solid State Phys 40, 143 (2000).
- H. v. Löhneysen, A. Rosch, M. Vojta, and P. Wölfle, Rev. Mod. Phys. 79, 1015 (2007).
- M. A. Ruderman and C. Kittel, Phys. Rev. 96, 99 (1954).
- T. Kasuya, Progress of Theoretical Physics 16, 45 (1956).
- K. Yosida, Phys. Rev. 106, 893 (1957).
- S. Doniach, Physica B+C 91, 231 (1977), ISSN 0378-4363.
- L. Zhou, J. Wiebe, S. Lounis, E. Vedmedenko, F. Meier, S. Blugel, P. H. Dederichs, and R. Wiesendanger, Nat. Phys. 6, 187 (2010).
- J.-H. Chen, L. Li, W. G. Cullen, E. D. Williams, and M. S. Fuhrer, Nat. Phys. 7, 1745 (2007).
- D. Hsieh, Y. Xia, L. Wray, D. Qian, A. Pal, J. H. Dil, J. Osterwalder, F. Meier, G. Bihlmayer, C. L. Kane, et al., Science 323, 919 (2009).
- Y. L. Chen, J.-H. Chu, J. G. Analytis, Z. K. Liu, K. Igarashi, H.-H. Kuo, X. L. Qi, S. K. Mo, R. G. Moore, D. H. Lu, et al., Science 329, 659 (2010).
- M. R. Scholz, J. Sánchez-Barriga, D. Marchenko, A. Varykhalov, A. Volykhov, L. V. Yashina, and O. Rader, Phys. Rev. Lett. 108, 256810 (2012).
- C. Seibel, H. Maaß, M. Ohtaka, S. Fiedler, C. Jünger, C.-H. Min, H. Bentmann, K. Sakamoto, and F. Reinert, Phys. Rev. B 86, 161105 (2012).
- P. W. Anderson, Nobel Lectures in Physics 1980, 376 (1977).
- Mott, N. F., J. Phys. Colloques 37, C4 (1976).
- E. Miranda and V. Dobrosavljevic, Reports on Progress in Physics 68, 2337 (2005).
- E. Miranda, V. Dobrosavljevic, and G. Kotliar, Journal of Physics: Condensed Matter 8, 9871 (1996).
- R. N. Bhatt and D. S. Fisher, Phys. Rev. Lett. 68, 3072 (1992).
- A. Langenfeld and P. WÃ¶lfle, Annalen der Physik 507, 43 (1995), ISSN 1521-3889.
- A. Castro Neto and B. Jones, Phys. Rev. B 62, 14975 (2000).
- P. S. Cornaglia, D. R. Grempel, and C. A. Balseiro, Phys. Rev. Lett. 96, 117209 (2006).
- S. Kettemann and E. R. Mucciolo, Phys. Rev. B 75, 184407 (2007).
- M.-T. Tran and K.-S. Kim, Phys. Rev. Lett. 105, 116403 (2010).
- A. Zhuravlev, I. Zharekeshev, E. Gorelov, A. I. Lichtenstein, E. R. Mucciolo, and S. Kettemann, Phys. Rev. Lett. 99, 247202 (2007).
- S. Kettemann and M. E. Raikh, Phys. Rev. Lett. 90, 146601 (2003).
- S. Kettemann, E. R. Mucciolo, I. Varga, and K. Slevin, Phys. Rev. B 85, 115112 (2012).
- I. V. Lerner, Phys. Rev. B 48, 9462 (1993).
- L. N. Bulaevskii and S. V. Panyukov, JETP Lett. 43, 240 (1986).
- G. Bergmann, Phys. Rev. B 36, 2469 (1987).
- H. Lee, J. Kim, E. R. Mucciolo, G. Bouzerar, and S. Kettemann, Phys. Rev. B 85, 075420 (2012a).
- H. Lee, E. R. Mucciolo, G. Bouzerar, and S. Kettemann, Phys. Rev. B 86, 205427 (2012b).
- Y. Nagaoka, Phys. Rev. 138, A1112 (1965).
- H. Suhl, Phys. Rev. 138, A515 (1965).
- S. Roche and D. Mayou, Phys. Rev. B 60, 322 (1999).
- A. Weiße, G. Wellein, A. Alvermann, and H. Fehske, Rev. Mod. Phys. 78, 275 (2006).
- E. R. Mucciolo, unpublished (2010).
- P. Wenk and G. Bouzerar, unpublished (2013).
- J. A. Sobota, D. Tanasković, and V. Dobrosavljević, Phys. Rev. B 76, 245106 (2007).
- S. Kettemann and E. Mucciolo, JETP Letters 83, 240 (2006), ISSN 0021-3640.
- K. Binder and A. P. Young, Rev. Mod. Phys. 58, 801 (1986).
- B. Coqblin, C. Lacroix, M. A. Gusmão, and J. R. Iglesias, Phys. Rev. B 67, 064417 (2003).
- C. M. Varma, Rev. Mod. Phys. 48, 219 (1976).
- S. G. Magalhaes, F. M. Zimmer, P. R. Krebs, and B. Coqblin, Phys. Rev. B 74, 014427 (2006).
- G. Bouzerar, T. Ziman, and J. Kudrnovsky, Eur. Phys. Lett. 69, 812 (2005).
- Y. C. Tsay and M. W. Klein, Phys. Rev. B 7, 352 (1973).
- Y. C. Tsay and M. W. Klein, Phys. Rev. B 11, 318 (1975).
- H. Tsunetsugu, M. Sigrist, and K. Ueda, Rev. Mod. Phys. 69, 809 (1997).
- P. Coleman, Phys. Rev. B 29, 3035 (1984).
- D. Newns and N. Read, Advances in Physics 36, 799 (1987).