Next Article in Journal
Preparation of MgO-Templated N-Doped Mesoporous Carbons from Polyvinylpyrrolidone: Effect of Heating Temperature on Pore Size Distribution
Previous Article in Journal
Novel Biobased Polyol Using Corn Oil for Highly Flame-Retardant Polyurethane Foams
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

RKKY Interaction in Graphene at Finite Temperature

1
Jack and Pearl Resnick Institute, Department of Physics, Bar-Ilan University, Ramat-Gan 52900, Israel
2
Max-Planck-Institut fur Physik komplexer Systeme, 01187 Dresden, Germany
Submission received: 27 February 2019 / Revised: 22 March 2019 / Accepted: 26 March 2019 / Published: 31 March 2019

Abstract

:
In our publication from eight years ago (Kogan, E. 2011, vol. 84, p. 115119), we calculated Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction between two magnetic impurities adsorbed on graphene at zero temperature. We show in this short paper that the approach based on Matsubara formalism and perturbation theory for the thermodynamic potential in the imaginary time and coordinate representation which was used then, can be easily generalized, and calculate RKKY interaction between the magnetic impurities at finite temperature.

Graphical Abstract

1. Introduction

More than 60 years ago, it was understood that localized spins in metals can interact by means of the Ruderman-Kittel-Kasuya-Yosida (RKKY) mechanism [1,2,3]. This indirect exchange between two magnetic impurities in a non–magnetic host coupling is mediated by the conduction electrons and is traditionally calculated as the second order perturbation with respect to exchange interaction between the magnetic impurity and the itinerant electrons of the host. Though analysis of the RKKY interaction in the lowest order non-zero of perturbation is simple in principle, calculation of the integrals defining the interaction (whether analytical or numerical) can pose some problems.
RKKY interaction has been investigated in materials of different nature such as disordered metals [4], superconductors [5,6,7], topological insulators [8,9,10,11,12,13,14,15], carbon nanotubes [16,17], semiconducting wires [18], in Weyl and Dirac semimetals [19,20,21,22,23,24], but most thoroughly in graphene [25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46]. RKKY interaction in graphene is also the subject of the present publication.
The initial motivation for our study of the RKKY interaction in 2011 [47] was the realization of the fact that many of the previous calculations had to deal with divergent integrals, and the complicated (and to some extent arbitrary) cut-off procedure was implemented to obtain from these integrals the finite results. So we started to look for the procedure which will allow eliminating this problem. It turned out that by using Matsubara formalism [48] and calculating the Matsubara Green’s functions in the coordinate-imaginary time representation, one is completely free from any diverging integrals in the whole calculation process [47].
Though Matsubara formalism is quite appropriate for finite temperature calculations (it was invented for that purpose) in our previous publications [35,47] we studied RKKY interaction only at zero temperature. The present short note is intended to generalize the result obtained previously to the case of finite temperature.
Adsorbed magnetic impurities can occupy different positions on graphene sheet. For example, they can be in the center of the elementary cell, or in between two adjacent graphene atoms [32,49,50,51]. However, in the previous publication we consider two magnetic impurities sitting on top of carbon atoms in graphene lattice, which is from our point of view the most interesting case. Let the two impurities sit on top of the sites i and j. Assuming a contact exchange interaction between the electrons and the magnetic impurities we can write the total Hamiltonian of the system as
H T = H + H i n t = H J S i · s i J S j · s j ,
where H is the Hamiltonian of the electron system, S i is the spins of the impurity and s i is the spin of itinerant electrons at site i.

2. Theoretical Methods

The consideration is based on the perturbation theory for the thermodynamic potential [48]. The correction to the thermodynamic potential due to interaction is
Δ Ω = T ln S T ln tr S · e H / T / Z ,
where the S—matrix is given by the equation
S = T τ exp 0 1 / T H i n t ( τ ) d τ .
Writing down s i in the second quantization representation
s i = 1 2 c i α σ α β c i β ,
the second order term of the expansion with respect to the interaction is
Δ Ω = J 2 T 4 α β γ δ S i · σ α β S j · σ γ δ 0 1 / T 0 1 / T d τ 1 d τ 2 T τ c i α ( τ 1 ) c i β ( τ 1 ) c j γ ( τ 2 ) c j δ ( τ 2 ) .
Please note that we have ignored the terms proportional to S i 2 and S j 2 , because they are irrelevant for our calculation of the effective interaction between the adatoms spins. Actually, such and similar terms of higher order can potentially lead to the Kondo effect in graphene [52,53], and its competition with the RKKY interaction [54,55,56], but we do not touch this issue, implicitly assuming that the RKKY interaction is strong enough to suppress the Kondo effect (at the temperature considered).
Leaving aside the question about the spin structure of the two–particle Green’s function standing in the r.h.s. of Equation (5) (for interacting electrons), further on we assume that the electrons are non–interacting. This will allow us to use Wick theorem and present the correlator from Equation (5) in the form
T τ c i α ( τ 1 ) c i β ( τ 1 ) c j γ ( τ 2 ) c j δ ( τ 2 ) = G β γ ( i , j ; τ 1 τ 2 ) G δ α ( j , i ; τ 2 τ 1 ) ,
where
G β γ ( i , j , τ 1 τ 2 ) = T τ c i β ( τ 1 ) c j γ ( τ 2 )
is the Matsubara Green’s function [48]. We can connect G β γ with the Green’s function of spinless electron
G β γ ( i , j , τ 1 τ 2 ) = δ β γ T τ c i ( τ 1 ) c j ( τ 2 ) .
Presence of delta-symbols allows to perform summation with respect to spin indices in Equation (5)
α β S i · σ α β S j · σ β α = S i · S j ,
which gives
Δ Ω = J 2 χ i j S i · S j ,
where
χ i j = 1 4 0 1 / T G ( i , j ; τ ) G ( j , i ; τ ) d τ
is the free electrons static real space spin susceptibility [47,57].
Thus we obtain
H R K K Y = J 2 χ i j S i · S j ,
Green’s function can be easily written down using representation of eigenvectors and eigenvalues of the operator H
H E n u n = 0 .
It is
G ( i , j ; τ ) = n u n ( i ) u n ( j ) e ξ n τ × 1 n F ( ξ n ) , τ > 0 n F ( ξ n ) , τ < 0 ,
where ξ n = E n μ , and n F ( ξ ) = e β ξ + 1 1 is the Fermi distribution function.
In calculations of the RKKY interaction in graphene the n in Equation (14) turns into a 2 ( 2 π ) 2 d 2 p , where a is the carbon–carbon distance. (Actually, there should appear a numerical multiplier, connecting the area of the elementary cell with a 2 , but we decided to discard it, which is equivalent to some numerical renormalization of J.) Also
u n ( i ) = e i p · R i ψ p ,
where ψ p is the appropriate component of spinor electron wave-function (depending upon which sublattice the magnetic adatom belongs to) in momentum representation.
Further in the integration with respect to d 2 p we will treat as the integration in the vicinity of two Dirac points K , K and present p = K ( K ) + k . The wave function for the momentum around Dirac points K and K has respectively the form
ψ ν , K ( k ) = 1 2 e i θ k / 2 ν e i θ k / 2 ψ ν , K ( k ) = 1 2 e i θ k / 2 ν e i θ k / 2 ,
where ν = ± 1 corresponds to electron and hole band [58]; the upper line of the spinor refers to the sublattice A and the lower line refers to the sublattice B.
In our publication from 2013 [35], we consider the case of doped graphene. However, here, like in our first publication on the subject [47], we consider only the case of undoped graphene, with the chemical potential at the Dirac points. The quantities E + ( k ) and E ( k ) in this case would be electron and hole energy. Then, Equation (14) takes the form: for i and j belonging to the same sublattice
G A A ( i , j ; τ > 0 ) = 1 2 e i K · R i j + e i K · R i j a 2 ( 2 π ) 2 d 2 k e i k · R i j E + ( k ) τ ,
and for i and j belonging to different sublattices
G A B ( i , j ; τ > 0 ) = 1 2 a 2 ( 2 π ) 2 d 2 k e E + ( k ) τ × e i ( K + k ) · R i j i θ k e i ( K + k ) · R i j + i θ k .
For τ < 0 we should change the sign of the Green’s functions and substitute E for E + .
For isotropic dispersion law E ( k ) = E ( k ) we can perform the angle integration in Equations (17) and (18) to get
d 2 k e i k · R i j E ( k ) τ = 0 d k k J 0 ( k R ) e E ( k ) τ d 2 k e i k · R i j ± i θ k E ( k ) τ = e ± i θ R 0 d k k J 1 ( k R ) e E ( k ) τ
( J 0 and J 1 are the Bessel function of zero and first order respectively, and θ R is the angle between the vectors K K and R i j ; R = | R i j | ).
For the linear dispersion law
E ± ( k ) = ± v F k ,
using mathematical identity [59]
0 x n 1 e p x J ν ( c x ) d x = ( 1 ) n 1 c ν n 1 p n 1 p 2 + c 2 p ν p 2 + c 2 ,
we can explicitly calculate the Green’s functions to get
G A A ( i , j ; τ > 0 ) = a 2 4 π v τ v 2 τ 2 + R 2 3 / 2 e i K · R i j + e i K · R i j G A B ( i , j ; τ > 0 ) = a 2 4 π R v 2 τ 2 + R 2 3 / 2
e i ( K + k ) · R i j i θ k e i ( K + k ) · R i j + i θ k .

3. Results

Now we have substitute the results obtained into Equation (11). In our previous publication [47], only the case T = 0 was considered. However, consideration of finite temperature is no more complicated in the formalism used. (We again realize the convenience of the imaginary time—coordinate representation of the Green’s functions for the problem at hand. The transition from zero to finite temperature result can be performed just by changing limits of integration in the single integral.) Thus we obtain for arbitrary T
χ T A A R i j = χ T = 0 A A R i j 16 π 0 v / R T x 2 d x ( x 2 + 1 ) 3
χ T A B R i j = χ T = 0 A B R i j 16 3 π 0 v / R T d x ( x 2 + 1 ) 3 ,
where χ T = 0 A A are the zero temperature susceptibilities calculated in our previous publication [47]
χ T = 0 A A R i j = a 4 256 v F R 3 1 + cos ( ( K K ) · R i j )
χ T = 0 A B R i j = 3 a 4 256 v F R 3 1 cos ( ( K K ) · R i j 2 θ R ) .
Actually, while rederiving Equations (26) and (27) we have found additional multiplier 1 / 2 π , but since we were already quite sloppy with the numerical multiplier in going from summation to integration in Equation (14), we decided to leave the equations in this paper as they were in the previously published one [47]. In any case, Equations (24) and (25) are valid independently of the presence of additional multiplier for χ T = 0 .
Integrals in Equations (24) and (25) can be easily calculated, and we obtain for the intrasublattice interaction
χ T A A R i j = χ T = 0 A A R i j 2 π tan 1 z + z z 2 + 1 2 z ( z 2 + 1 ) 2 ,
and for the intersublattice interaction
χ T A B R i j = χ T = 0 A B R i j 2 π tan 1 z + z z 2 + 1 + z 3 ( z 2 + 1 ) 2 ,
where z = v / R τ . The limiting cases of Equations (28) and (29) are particularly simple. For T v / R the first term in the braces both in Equation (28) and in Equation (29) is equal to π / 2 and the other two terms can be neglected. Thus we obtain the previous ( T = 0 ) results. The opposite limiting case T v / R is easier to get directly from Equations (24) and (25). Expanding the integrand we get for T v / R
χ T A A R i j = χ T = 0 A A R i j 16 3 π v R T 3
χ T A B R i j = χ T = 0 A B R i j 16 3 π v R T .
We must mention that comparing our results with those obtained earlier for the case of doped graphene [60], one should be aware of the fact that the exponential decrease of the RKKY interaction with the distance at high temperatures obtained in Ref. [60], was obtained for k F R 1 (in our case k F = 0 ).

4. Conclusions

We have calculated RKKY interaction between two magnetic impurities adsorbed on graphene at arbitrary temperature. Matsubara Green’s functions in coordinate and imaginary time representations were used. For zero temperature the results coinsides with those obtained by us previously. At high temperature the interaction decreases inversely proportional to to the temperature for the inter-sublattice interaction, and inversely proportional to the third power of temperature fore the intra-sublattice interaction.

Funding

This research received no external funding.

Acknowledgments

This paper was written during the author’s visit to Max-Planck-Institut fur Physik komplexer Systeme in 2019. The author cordially thanks the Institute for the hospitality extended to him during that and all the previous visits.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ruderman, M.A.; Kittel, C. Indirect Exchange Coupling of Nuclear Magnetic Moments by Conduction Electrons. Phys. Rev. 1954, 96, 99–102. [Google Scholar] [CrossRef]
  2. Kasuya, T. A Theory of Metallic Ferro- and Antiferromagnetism on Zener’s Model. Prog. Theor. Phys. 1956, 16, 45–57. [Google Scholar] [CrossRef]
  3. Yosida, K. Magnetic Properties of Cu-Mn Alloys. Phys. Rev. 1957, 106, 893–898. [Google Scholar] [CrossRef]
  4. Zyuzin, Y.A.; Spivak, B.Z. Friedel oscillations and Ruderman-Kittel interaction in disordered conductors. JETP Lett. 1986, 43, 234–237. [Google Scholar]
  5. Abrikosov, A.A. Fundamentals of the Theory of Metals; Elsevier Science Publishers: Amsterdam, The Netherlands, 1988. [Google Scholar]
  6. Aristov, D.N.; Maleyev, S.V.; Yashenkin, A.G. RKKY interaction in layered superconductors with anisotropic pairing. Z. Phys. B Condens. 1997, 102, 467–471. [Google Scholar] [CrossRef]
  7. Galitski, V.M.; Larkin, A.I. Spin glass versus superconductivity. Phys. Rev. B 2002, 66, 064526. [Google Scholar] [CrossRef]
  8. Liu, Q.; Liu, C.-X.; Xu, C.; Qi, X.-L.; Zhang, S.-C. Magnetic Impurities on the Surface of a Topological Insulator. Phys. Rev. Lett. 2009, 102, 156603. [Google Scholar] [CrossRef]
  9. Biswas, R.R.; Balatsky, A.V. Impurity-induced states on the surface of three-dimensional topological insulators. Phys. Rev. B 2010, 81, 233405. [Google Scholar] [CrossRef]
  10. Garate, I.; Franz, M. Magnetoelectric response of the time-reversal invariant helical metal. Phys. Rev. B 2010, 81, 172408. [Google Scholar] [CrossRef]
  11. Zhu, J.-J.; Yao, D.-X.; Zhang, S.-C.; Chang, K. Electrically Controllable Surface Magnetism on the Surface of Topological Insulators. Phys. Rev. Lett. 2011, 106, 097201. [Google Scholar] [CrossRef] [PubMed]
  12. Abanin, D.A.; Pesin, D.A. Ordering of Magnetic Impurities and Tunable Electronic Properties of Topological Insulators. Phys. Rev. Lett. 2011, 106, 136802. [Google Scholar] [CrossRef] [PubMed]
  13. Checkelsky, J.G.; Ye, J.; Onose, Y.; Iwasa, Y.; Tokura, Y. Dirac-fermion-mediated ferromagnetism in a topological insulator. Nat. Phys. 2012, 8, 729–733. [Google Scholar] [CrossRef]
  14. Zyuzin, A.A.; Loss, D. RKKY interaction on surfaces of topological insulators with superconducting proximity effect. Phys. Rev. B 2014, 90, 125443. [Google Scholar] [CrossRef]
  15. Klier, N.; Sharma, S.; Pankratov, O.; Shallcross, S. Tuning topological surface magnetism by bulk alloying. arXiv 2017, arXiv:1711.10760. [Google Scholar]
  16. Braunecker, B.; Simon, P.; Loss, D. Nuclear Magnetism and Electronic Order in 13C Nanotubes. Phys. Rev. Lett. 2009, 102, 116403. [Google Scholar] [CrossRef]
  17. Klinovaja, J.; Loss, D. RKKY interaction in carbon nanotubes and graphene nanoribbons. Phys. Rev. B 2013, 87, 045422. [Google Scholar] [CrossRef]
  18. Braunecker, B.; Simon, P.; Loss, D. Nuclear magnetism and electron order in interacting one-dimensional conductors. Phys. Rev. B 2009, 80, 165119. [Google Scholar] [CrossRef]
  19. Sun, J.-H.; Xu, D.-H.; Zhang, F.-C.; Zhou, Y. Magnetic impurity in a Weyl semimetal. Phys. Rev. B 2015, 92, 195124. [Google Scholar] [CrossRef]
  20. Chang, H.-R.; Zhou, J.; Wang, S.-X.; Shan, W.; Xiao, D. RKKY interaction of magnetic impurities in Dirac and Weyl semimetals. Phys. Rev. B 2015, 92, 241103. [Google Scholar] [CrossRef]
  21. Hosseini, M.V.; Askari, M. Ruderman-KittelKasuya-Yosida interaction in Weyl semimetals. Phys. Rev. B 2015, 92, 224435. [Google Scholar] [CrossRef]
  22. Sun, Y.; Wang, A. RKKY interaction of magnetic impurities in multi-Weyl semimetals. J. Phys. Condens. Matter 2017, 2, 435306. [Google Scholar] [CrossRef] [PubMed]
  23. Yevtushenko, O.M.; Yudson, V.I. Kondo Impurities Coupled to a Helical Luttinger Liquid: RKKY-Kondo Physics Revisited. Phys. Rev. Lett. 2018, 120, 147201. [Google Scholar] [CrossRef]
  24. Kaladzhyan, V.; Zyuzin, A.A.; Simon, P. RKKY interaction on the surface of 3D Dirac semimetals. arXiv 2018, arXiv:1811.05750. [Google Scholar]
  25. Vozmediano, M.A.H.; Lopez-Sancho, M.P.; Stauber, T.; Guinea, F. Local defects and ferromagnetism in graphene layers. Phys. Rev. B 2005, 72, 155121. [Google Scholar] [CrossRef]
  26. Dugaev, V.K.; Litvinov, V.I.; Barnas, J. Exchange interaction of magnetic impurities in graphene. Phys. Rev. B 2006, 74, 224438. [Google Scholar] [CrossRef]
  27. Brey, L.; Fertig, H.A.; Sarma, S.D. Diluted Graphene Antiferromagnet. Phys. Rev. Lett. 2007, 99, 116802. [Google Scholar] [CrossRef]
  28. Saremi, S. RKKY in half-filled bipartite lattices: Graphene as an example. Phys. Rev. B 2007, 76, 184430. [Google Scholar] [CrossRef]
  29. Hwang, E.H.; Sarma, S.D. Screening, Kohn Anomaly, Friedel Oscillation, and RKKY Interaction in Bilayer Graphene. Phys. Rev. Lett. 2008, 101, 156802. [Google Scholar] [CrossRef] [PubMed]
  30. Black–Schaffer, A.M. RKKY coupling in graphene. Phys. Rev. B 2010, 81, 205416. [Google Scholar] [CrossRef]
  31. Sherafati, M.; Satpathy, S. RKKY interaction in graphene from the lattice Green’s function. Phys. Rev. B 2011, 83, 165425. [Google Scholar] [CrossRef]
  32. Uchoa, B.; Rappoport, T.G.; Castro Neto, A.H. Kondo Quantum Criticality of Magnetic Adatoms in Graphene. Phys. Rev. Lett. 2011, 106, 016801. [Google Scholar] [CrossRef] [PubMed]
  33. Power, S.R.; Guimaraes, F.M.S.; Costa, A.T.; Muniz, R.B.; Ferreira, M.S. Dynamic RKKY interaction in graphene. Phys. Rev. B 2012, 85, 195411. [Google Scholar] [CrossRef]
  34. Patrone, P.N.; Einstein, T.L. Anisotropic surface-state-mediated RKKY interaction between adatoms. Phys. Rev. B 2012, 85, 045429. [Google Scholar] [CrossRef]
  35. Kogan, E. RKKY Interaction in Gapped or Doped Graphene. Graphene 2013, 2, 8–12. [Google Scholar] [CrossRef]
  36. Gorman, P.D.; Duffy, J.M.; Ferreira, M.S.; Power, S.R. RKKY interaction between adsorbed magnetic impurities in graphene: Symmetry and strain effects. Phys. Rev. B 2013, 88, 085405. [Google Scholar] [CrossRef]
  37. Stano, P.; Klinovaja, J.; Yacoby, A.; Loss, D. Local spin susceptibilities of low-dimensional electron systems. Phys. Rev. B 2013, 88, 045441. [Google Scholar] [CrossRef]
  38. Szalowski, K. Indirect RKKY interaction between localized magnetic moments in armchair graphene nanoribbons. J. Phys. Condens. Matter 2013, 25, 166001. [Google Scholar] [CrossRef]
  39. Roslyak, O.; Gumbs, G.; Huang, D. Gap-modulated doping effects on indirect exchange interaction between magnetic impurities in graphene. J. Appl. Phys. 2013, 1, 123702. [Google Scholar] [CrossRef]
  40. Xiao, X.; Liu, Y.; Wen, W. Ruderman-Kittel-Kasuya-Yosida interaction in silicene. J. Phys. Condens. Matter 2014, 26, 266001. [Google Scholar] [CrossRef] [PubMed]
  41. Shirakawa, T.; Yunoki, S. Block Lanczos density-matrix renormalization group method for general Anderson impurity models: Application to magnetic impurity problems in graphene. Phys. Rev. B 2014, 90, 195109. [Google Scholar] [CrossRef]
  42. Zare, M.; Parhizgar, F.; Asgari, R. Topological phase and edge states dependence of the RKKY interaction in zigzag silicene nanoribbon. Phys. Rev. B 2016, 94, 045443. [Google Scholar] [CrossRef]
  43. Fransson, J.; Black-Schaffer, A.M.; Balatsky, A.V. Magnon Dirac materials. Phys. Rev. B 2016, 94, 075401. [Google Scholar] [CrossRef]
  44. Rezania, H.; Azizi, F. RKKY interaction in spin polarized armchair graphene nanoribbon. JMMM 2016, 417, 272–278. [Google Scholar] [CrossRef]
  45. Agarwal, M.; Mishchenko, E.G. Long-range exchange interaction between magnetic impurities in graphene. Phys. Rev. B 2017, 95, 075411. [Google Scholar] [CrossRef]
  46. Sousa, J.F.; Amorim, B.; Castro, E.V. Dilute magnetism in graphene. arXiv 2019, arXiv:1901.08614. [Google Scholar]
  47. Kogan, E. RKKY interaction in graphene. Phys. Rev. B 2011, 84, 115119. [Google Scholar] [CrossRef]
  48. Abrikosov, A.A.; Gorkov, L.P.; Dzyloshinski, I.E. Methods of Quantum Field Theory in Statistical Physics; Pergamon Press: Oxford, UK, 1965. [Google Scholar]
  49. Wang, S.; Yang, B.; Chen, H.; Ruckenstein, E. Popgraphene: a new 2D planar carbon allotrope composed of 5–8–5 carbon rings for high-performance lithium-ion battery anodes from bottom-up programming. J. Mater. Chem. A 2018, 6, 6815–6821. [Google Scholar] [CrossRef]
  50. Wang, S.; Yang, B.; Chen, H.; Ruckenstein, E. Reconfiguring graphene for high-performance metal-ion battery anodes. Energy Storage Mater. 2019, 16, 619–624. [Google Scholar] [CrossRef]
  51. Wu, D.; Yang, B.; Ruckenstein, E.; Chen, H. Functionalization: An Effective Approach to Open and Close Channels for Electron Transfer in Nitrogenated Holey Graphene C2N Anodes in Sodium-Ion Batteries. J. Phys. Chem. Lett. 2019, 10, 721–726. [Google Scholar] [CrossRef]
  52. Kogan, E.; Noda, K.; Yunoki, S. Spin-anisotropic magnetic impurity in a Fermi gas: Integration of poor man’s scaling equations. Phys. Rev. B 2017, 95, 165412. [Google Scholar] [CrossRef]
  53. Kogan, E. Poor man’s scaling: anisotropic Kondo and Coqblin–Schrieffer models. J. Phys. Commun. 2018, 2, 085001. [Google Scholar] [CrossRef]
  54. Pruser, H.; Dargel, P.E.; Bouhassoune, M.; Ulbrich, R.G.; Pruschke, T.; Lounis, S.; Wenderoth, M. Interplay between the Kondo effect and the Ruderman–Kittel–Kasuya–Yosida interaction. Nat. Commun. 2014, 5, 5417. [Google Scholar] [CrossRef]
  55. Nejati, A.; Kroha, J. Oscillation and suppression of Kondo temperature by RKKY coupling in two-site Kondo systems. J. Phys. Conf. Ser. 2017, 807, 082004. [Google Scholar] [CrossRef]
  56. Kroha, J. Ser. Model. Simulat. 2017; 7, 12.1–12.27, (Verlag Forschungszentrum Julich, 2017).
  57. Cheianov, V.V.; Syljuasen, O.; Altshuler, B.L.; Fal’ko, V. Ordered states of adatoms on graphene. Phys. Rev. B 2009, 80, 233409. [Google Scholar] [CrossRef]
  58. Castro Neto, A.H.; Guinea, F.; Peres, N.M.R.; Novoselov, K.S.; Geim, A.K. The electronic properties of graphene. Rev. Mod. Phys. 2009, 81, 109–162. [Google Scholar] [CrossRef]
  59. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series; Gordon and Breach Science Publishers: London, UK, 1986; Volume 2. [Google Scholar]
  60. Klier, N.; Shallcross, S.; Sharma, S.; Pankratov, O. Ruderman-Kittel-Kasuya-Yosida interaction at finite temperature: Graphene and bilayer graphene. Phys. Rev. B 2015, 92, 205414. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Kogan, E. RKKY Interaction in Graphene at Finite Temperature. C 2019, 5, 14. https://0-doi-org.brum.beds.ac.uk/10.3390/c5020014

AMA Style

Kogan E. RKKY Interaction in Graphene at Finite Temperature. C. 2019; 5(2):14. https://0-doi-org.brum.beds.ac.uk/10.3390/c5020014

Chicago/Turabian Style

Kogan, Eugene. 2019. "RKKY Interaction in Graphene at Finite Temperature" C 5, no. 2: 14. https://0-doi-org.brum.beds.ac.uk/10.3390/c5020014

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop