RKKY interaction in graphene at finite temperature

TL;DR

This paper extends previous zero-temperature calculations of RKKY interaction in graphene to finite temperatures using Matsubara formalism, providing a more comprehensive understanding of magnetic impurity interactions in graphene.

Contribution

It generalizes the earlier zero-temperature approach to finite temperatures, enabling more accurate modeling of magnetic interactions in graphene under realistic conditions.

Findings

RKKY interaction can be calculated at finite temperature using the generalized approach.

Finite temperature effects modify the strength and range of magnetic interactions.

The method simplifies previous calculations by using Matsubara formalism.

Abstract

In our publication from 8 years ago (Phys. Rev. B {\bf 84}, 115119 (2011)) we calculated 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.

Equations

$$H_T=H+H_{int}=H-J{\mathbf{S}}_i{\cdot \mathbf{s}}_i-J{\mathbf{S}}_j{\cdot \mathbf{s}}_j,$$

$$\Delta \Omega =-T\ln ⟨S⟩\equiv -T\ln \mathrm{tr}\left\{S\cdot e^{-H/T}/Z\right\},$$

$$S=T_{\tau }\exp \left\{-{\int }_0^{1/T}{H}_{int}(\tau )d\tau \right\}.$$

$${\mathbf{s}}_i=\frac{1}{2}{c}_{i\alpha }^{†}{\sigma }_{\alpha \beta }{c}_{i\beta },$$

$$\Delta \Omega =\frac{J^{2}T}{4}\sum _{\alpha \beta \gamma \delta }{\mathbf{S}}_i{\cdot \sigma }_{\alpha \beta }{\mathbf{S}}_j{\cdot \sigma }_{\gamma \delta }$$

$${\int }_0^{1/T}{\int }_0^{1/T}d{\tau }_{1}d{\tau }_2⟨T_{\tau }\left\{c_{i\alpha }^{†}({\tau }_1)c_{i\beta }({\tau }_1)c_{j\gamma }^{†}({\tau }_2)c_{j\delta }({\tau }_2)\right\}⟩.$$

$$⟨T_{\tau }\left\{c_{i\alpha }^{†}({\tau }_1)c_{i\beta }({\tau }_1)c_{j\gamma }^{†}({\tau }_2)c_{j\delta }({\tau }_2)\right\}⟩$$

$$=-{\mathcal{G}}_{\beta \gamma }(i,j;{\tau }_1-{\tau }_2){\mathcal{G}}_{\delta \alpha }(j,i;{\tau }_2-{\tau }_1),$$

$${\mathcal{G}}_{\beta \gamma }(i,j,{\tau }_1-{\tau }_2)=-⟨T_{\tau }\left\{c_{i\beta }({\tau }_1)c_{j\gamma }^{†}({\tau }_2)\right\}⟩$$

$${\mathcal{G}}_{\beta \gamma }(i,j,{\tau }_1-{\tau }_2)=-{\delta }_{\beta \gamma }⟨T_{\tau }\left\{c_i({\tau }_1)c_j^{†}({\tau }_2)\right\}⟩.$$

$$\sum _{\alpha \beta }{\mathbf{S}}_i{\cdot \sigma }_{\alpha \beta }{\mathbf{S}}_j{\cdot \sigma }_{\beta \alpha }={\mathbf{S}}_i{\cdot \mathbf{S}}_j,$$

$$\Delta \Omega =-J^2{\chi }_{ij}{\mathbf{S}}_i{\cdot \mathbf{S}}_j,$$

$${\chi }_{ij}=-\frac{1}{4}{\int }_0^{1/T}\mathcal{G}(i,j;\tau )\mathcal{G}(j,i;-\tau )d\tau$$

$$H_{RKKY}=-J^2{\chi }_{ij}{\mathbf{S}}_i{\cdot \mathbf{S}}_j,$$

$$\left(H-E_n\right)u_n=0.$$

$$\mathcal{G}(i,j;\tau )=\sum _n{u}_n^{\ast }(i)u_n(j)e^{-{\xi }_n\tau }$$

$$\displaystyle\times\left\{\begin{array}[]{ll}-\left(1-n_{F}(\xi_{n})\right),&\tau>0\\ n_{F}(\xi_{n}),&\tau<0\end{array}\right.,$$

$$u_n(i)=e^{i{\mathbf{p}\cdot \mathbf{R}}_i}{\psi }_{\mathbf{p}},$$

$${\psi }_{\nu ,\mathbf{K}}(\mathbf{k})$$

$${\psi }_{\nu ,{\mathbf{K}}^{\prime }}(\mathbf{k})$$

$${\mathcal{G}}^{AA}(i,j;\tau >0)=-\frac{1}{2}\left[e^{i{\mathbf{K}\cdot \mathbf{R}}_{ij}}+e^{i{{\mathbf{K}}^{\prime }\cdot \mathbf{R}}_{ij}}\right]$$

$$\frac{a^2}{(2\pi {)}^2}\int d^2\mathbf{k}{e}^{i{\mathbf{k}\cdot \mathbf{R}}_{ij}-E_{+}(\mathbf{k})\tau },$$

$${\mathcal{G}}^{AB}(i,j;\tau >0)=\frac{1}{2}\frac{a^2}{(2\pi {)}^2}\int d^2\mathbf{k}{e}^{-E_{+}(\mathbf{k})\tau }$$

$$\times \left[e^{i(\mathbf{K}+\mathbf{k}){\cdot \mathbf{R}}_{ij}-i{\theta }_k}-e^{i({\mathbf{K}}^{\prime }+\mathbf{k}){\cdot \mathbf{R}}_{ij}+i{\theta }_k}\right].$$

$$\int d^2\mathbf{k}{e}^{i{\mathbf{k}\cdot \mathbf{R}}_{ij}-E(k)\tau }={\int }_0^{\infty }dkk{J}_0(kR)e^{-E(k)\tau }$$

$$\int d^2\mathbf{k}{e}^{i{\mathbf{k}\cdot \mathbf{R}}_{ij}\pm i{\theta }_k-E(k)\tau }=e^{\pm i{\theta }_{\mathbf{R}}}{\int }_0^{\infty }dkk{J}_1(kR)e^{-E(k)\tau }$$

$$E_{\pm }(k)=\pm v_{F}k,$$

$${\int }_0^{\infty }{x}^{n-1}{e}^{-px}{J}_{\nu }(cx)dx$$

$$=(-1{)}^{n-1}{c}^{-\nu }\frac{{\partial }^{n-1}}{\partial p^{n-1}}\frac{{\left(\sqrt{p^2+c^2}-p\right)}^{\nu }}{\sqrt{p^2+c^2}},$$

$${\mathcal{G}}^{AA}(i,j;\tau >0)=-\frac{a^2}{4\pi }\frac{v\tau }{{\left(v^2{\tau }^2+R^2\right)}^{3/2}}$$

$$\left[e^{i{\mathbf{K}\cdot \mathbf{R}}_{ij}}+e^{i{{\mathbf{K}}^{\prime }\cdot \mathbf{R}}_{ij}}\right]$$

$${\mathcal{G}}^{AB}(i,j;\tau >0)=\frac{a^2}{4\pi }\frac{R}{{\left(v^2{\tau }^2+R^2\right)}^{3/2}}$$

$$\left[e^{i(\mathbf{K}+\mathbf{k}){\cdot \mathbf{R}}_{ij}-i{\theta }_k}-e^{i({\mathbf{K}}^{\prime }+\mathbf{k}){\cdot \mathbf{R}}_{ij}+i{\theta }_k}\right].$$

$${\chi }_T^{AA}\left({\mathbf{R}}_{ij}\right)={\chi }_{T=0}^{AA}\left({\mathbf{R}}_{ij}\right)\frac{16}{\pi }{\int }_0^{v/RT}\frac{x^{2}dx}{(x^2+1{)}^3}$$

$${\chi }_T^{AB}\left({\mathbf{R}}_{ij}\right)={\chi }_{T=0}^{AB}\left({\mathbf{R}}_{ij}\right)\frac{16}{3\pi }{\int }_0^{v/RT}\frac{dx}{(x^2+1{)}^3},$$

$${\chi }_{T=0}^{AA}\left({\mathbf{R}}_{ij}\right)=\frac{a^4}{256v_F{R}^3}\left[1+\cos ((\mathbf{K}-{\mathbf{K}}^{\prime }){\cdot \mathbf{R}}_{ij})\right]$$

Full text

I Introduction

More than 60 years ago it was understood that localized spins in metals can interact by means of Ruderman-Kittel-Kasuya-Yosida (RKKY) mechanism ruderman ; kasuya ; yosida . 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 zyuzin , superconductors abr ; aristov ; gali , topological insulators liu ; biswas ; garate ; zhu ; abanin ; check ; zyuz ; pank , carbon nanotubes braun ; klin , semiconducting wires simon , in Weyl and Dirac semimetals sun1 ; chang ; hosseini ; sun2 ; yudson ; zyuzin0 , but most thoroughly in graphene vozmediano ; dugaev ; brey ; saremi ; hwang ; black ; sherafati ; uchoa ; power ; patrone ; kogan2 ; gorman ; stano ; salo ; gumbs ; xiao ; shirakawa ; zare ; frans ; hanini ; agarwal ; sousa . RKKY interaction in graphene is also the subject of the present publication.

Initial motivation for our study of the RKKY interaction in 2011 kogan was 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 to eliminate this problem. It turned out that using Matsubara formalism abrikosov 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 kogan .

Though Matsubara formalism is quite appropriate for finite temperature calculations (it was invented for that purpose) in our previous publications kogan ; kogan2 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 uchoa ; ruckenstein . 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 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

[TABLE]

where $H$ is the Hamiltonian of the electron system, ${\bf S}_{i}$ is the spins of the impurity and ${\bf s}_{i}$ is the spin of itinerant electrons at site $i$.

II Theoretical methods

The consideration is based on the perturbation theory for the thermodynamic potential abrikosov . The correction to the thermodynamic potential due to interaction is

[TABLE]

where the $S$–matrix is given by the equation

[TABLE]

Writing down ${\bf s}_{i}$ in the second quantization representation

[TABLE]

the second order term of the expansion with respect to the interaction is

[TABLE]

Notice that we have ignored the terms proportional to ${\bf S}_{i}^{2}$ and ${\bf 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 kogan3 ; kogan4 , and its competition with the RKKY interaction pruser ; kroha1 ; kroha2 , 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 Eq. (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 Eq. (5) in the form

[TABLE]

where

[TABLE]

is the Matsubara Green’s function abrikosov . We can connect ${\cal G}_{\beta\gamma}$ with the Green’s function of spinless electron

[TABLE]

Presence of delta-symbols allows to perform summation with respect to spin indices in Eq. (5)

[TABLE]

which gives

[TABLE]

where

[TABLE]

is the free electrons static real space spin susceptibility cheianov ; kogan .

Thus we obtain

[TABLE]

The Green’s function can be easily written down using representation of eigenvectors and eigenvalues of the operator $H$

[TABLE]

It is

[TABLE]

where $\xi_{n}=E_{n}-\mu$, and $n_{F}(\xi)=\left(e^{\beta\xi}+1\right)^{-1}$ is the Fermi distribution function.

In calculations of the RKKY interaction in graphene the $\sum_{n}$ in Eq. (II) turns into $\frac{a^{2}}{(2\pi)^{2}}\int d^{2}{\bf 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

[TABLE]

where $\psi_{\bf p}$ is the appropriate component of spinor electron wave-function (depending upon which sublattice the magnetic adatom belongs to) in momentum representation.

Further on the integration with respect to $d^{2}{\bf p}$ we’ll treat as the integration in the vicinity of two Dirac points $K,K^{\prime}$ and present ${\bf p}={\bf K}({\bf K}^{\prime})+{\bf k}$. The wave function for the momentum around Dirac points $K$ and $K^{\prime}$ has respectively the form

[TABLE]

where $\nu=\pm 1$ corresponds to electron and hole band castro ; 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 kogan2 we consider the case of doped graphene. But here, like in our first publication on the subject kogan , we consider only the case of undoped graphene, with the chemical potential at the Dirac points. the quantities $E_{+}({\bf k})$ and $E_{-}({\bf k})$ in this case would be electron and hole energy. Then Eq. (II) takes the form: for $i$ and $j$ belonging to the same sublattice

[TABLE]

and for $i$ and $j$ belonging to different sublattices

[TABLE]

For $\tau<0$ we should change the sign of the Green’s functions and substitute $E_{-}$ for $E_{+}$.

For isotropic dispersion law $E({\bf k})=E(k)$ we can perform the angle integration in Eqs. (II) and (II) to get

[TABLE]

($J_{0}$ and $J_{1}$ are the Bessel function of zero and first order respectively, and $\theta_{\bf R}$ is the angle between the vectors ${\bf K}-{\bf K}^{\prime}$ and ${\bf R}_{ij}$; $R=|{\bf R}_{ij}|$).

For the linear dispersion law

[TABLE]

using mathematical identity prudnikov

[TABLE]

we can explicitly calculate the Green’s functions to get

[TABLE]

III Results

Now we have substitute the results obtained into Eq. (11). In our previous publication kogan only the case $T=0$ was considered. But 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$

[TABLE]

where $\chi^{AA}_{T=0}$ are the zero temperature susceptibilities calculated in our previous publication kogan

[TABLE]

Actually, while rederiving Eqs. (33), (III) we have found additional multiplier $1/2\pi$, but since we were already quite sloppy with the numerical multiplier in going from summation to integration in Eq. (II), we decided to leave the equations in this paper as they were in the previously published one kogan . In any case, Eqs. (III) and (III) are valid independently of the presence of additional multiplier for $\chi_{T=0}$.

Integrals in Eqs. (III), (III) can be easily calculated, and we obtain for the intrasublattice interaction

[TABLE]

and for the intersublattice interaction

[TABLE]

where $z=v/R\tau$. The limiting cases of Eqs. (III) and (III) are particularly simple. For $T\ll v/R$ the first term in the braces both in Eq. (III) and in Eq. (III) is equal to $\pi/2$ and the other two terms can be neglected. Thus we obtain the previous ($T=0$) results. The opposite limiting case $T\gg v/R$ is easier to get directly from Eqs. (III) and (III). Expanding the integrand we get for $T\gg v/R$

[TABLE]

We must mention that comparing our results with those obtained earlier for the case of doped graphene klier , one should be aware of the fact that the exponential decrease of the RKKY interaction with the distance at high temperatures obtained in Ref. klier, , was obtained for $k_{F}R\gg 1$ (in our case $k_{F}=0$).

IV 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.

Acknowledgements.

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.