Next Article in Journal
Nuclear Equation of State in the Relativistic Point-Coupling Model Constrained by Excitations in Finite Nuclei
Next Article in Special Issue
Casimir Effect for Fermion Condensate in Conical Rings
Previous Article in Journal / Special Issue
Role of Electronic Relaxation Rates in the Casimir Force between High-Tc Superconductors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Low-Temperature Expansion of the Casimir-Polder Free Energy of an Atom with Graphene

by
Nail Khusnutdinov
* and
Natalia Emelianova
Centro de Matemática, Computação e Cognição, Universidade Federal do ABC, Santo André 09210-170, SP, Brazil
*
Author to whom correspondence should be addressed.
Universe 2021, 7(3), 70; https://doi.org/10.3390/universe7030070
Submission received: 4 February 2021 / Revised: 8 March 2021 / Accepted: 12 March 2021 / Published: 16 March 2021
(This article belongs to the Special Issue The Casimir Effect: From a Laboratory Table to the Universe)

Abstract

:
We consider the low-temperature expansion of the Casimir-Polder free energy for an atom and graphene by using the Poisson representation of the free energy. We extend our previous analysis on the different relations between chemical potential μ and mass gap parameter m. The key role plays the dependence of graphene conductivities on the μ and m. For simplicity, we made the manifest calculations for zero values of the Fermi velocity. For μ > m , the thermal correction T 2 , and for μ < m , we confirm the recent result of Klimchitskaya and Mostepanenko, that the thermal correction T 5 . In the case of exact equality μ = m , the correction T . This point is unstable, and the system falls to the regime with μ > m or μ < m . The analytical calculations are illustrated by numerical evaluations for the Hydrogen atom/graphene system.

1. Introduction

The Casimir [1] and Casimir-Polder [2] dispersion forces play an important role in different phenomena [3,4]. The Casimir-Polder force is usually referred to as the van der Waals force on large distances between micro-particles and macro-objects when the retardation of interaction is taken into account. The Casimir-Polder force essentially depends on the material of the macro-objects, its dimension, shapes, conductivity, and temperature [3,4], and it is important for the interaction of graphene with micro-particles [5,6,7,8,9]. The Casimir-Polder force and torque for anisotropic molecules have been the subject of investigations in the recent years from the theoretical point of view [10,11,12,13,14], as well as experimentally [15,16]. The non-ground state of an atom gives a specific resonant term contribution to the Casimir-Polder force [17,18], but, here, we restrict ourself by an atom in the ground state, only.
The thermal corrections to the Casimir-Polder interaction for micro-particle/graphene were considered in Refs. [19,20,21,22,23,24,25]. In the last few years, much attention has been given to the low-temperature expansion of the Casimir-Polder free energy for the atom/graphene system [7,23,24,25]. In the case of an atom/ideal plane, the low-temperature correction to the Casimir-Polder free energy is proportional to the fourth degree of temperature T 4 . As opposed to the ideal case, the conductivity of graphene depends on the chemical potential and temperature, and it has temporal and spatial dispersion [26,27,28,29]. The low-temperature expansion depends on the relations between these macro-parameters.
It was shown in Ref. [23] that the low-temperature expansion reveals the unusual quadratic T 2 behavior. Next, detailed considerations [24,25] showed a more rich picture of low-temperature expansion depending on the relation between chemical potential μ and mass gap parameter m of the Dirac electron. For μ > m , the same quadratic behavior was confirmed, but, in the case μ < m the T 5 , dependence was obtained. In the case of the very specific exact relation μ = m , the linear T dependence was observed. To obtain these results, the authors of Refs. [24,25] made the sophisticated treatment of the Matsubara series.
In the present paper, we extend our analysis made in Ref. [23] in the framework of the Poisson representation of the Matsubara series to all relations between μ and m and confirm results by numerical analysis. The Poisson representation is more suitable for the low-temperature expansions—the zeroth term of the series coincides with a zero-temperature contribution (with possible temperature and chemical potential dependence via conductivity) and the rest of the series gives temperature correction. In Ref. [23], we considered conductivity of graphene with cutting scattering rate, as in Refs. [30,31], where the Kubo approach was used. In this case, the conductivity has a constant value at zero frequencies which depends on the scattering rate parameter γ . It leads to the low-temperature dependence T 2 for any relation between chemical potential μ and mass gap m. In the framework of the polarization tensor approach [27,29], there is no scattering parameter, and the behavior of the conductivity at zero frequency strongly depends on the relation between μ and m.
In the framework of our approach, we show that, for μ > m , the temperature corrections to the free energy T 2 , and, for μ < m , we obtain correction T 5 , and, in the point μ = m , the linear dependence T appears. Therefore, we confirm expansions obtained in Ref. [24]. The main contribution for μ m to the low-temperature expansion of the free energy comes from the expression for energy at zero temperatures (zeroth term of Poisson expansion) via temperature dependence of conductivity. For μ < m , the main contribution comes from the rest part of the series. The point μ = m is a very specific point where the free energy is linear over temperature; therefore, the entropy is constant, and the Nernst theorem is no longer valid. But, for an infinitely small deviation from equality, μ = m shows different regimes where the Nernst theorem holds. The numerical evaluations reveal the singular “beak-shaped” form of the temperature correction of the free energy, which confirms this conclusion. The point μ = m is the point of changing a regime of the temperature dependence of the free energy. The numerical evaluations show that, for any small value of the deviations | μ m | , there is a domain of temperatures T close to the origin where the derivative of the free energy with respect to the temperature, the entropy, is zero for T = 0 , and the Nernst theorem is valid. The method developed here may be used to obtain the low-temperature expansion in many other situations.
Throughout the paper, the units = c = k B = 1 are used.

2. The Casimir-Polder Free Energy

Taking into account the Poisson summation formula (see details in Ref. [23]), the free energy may be represented in the following form
F tm te E CP = 8 3 l = 0 0 z 3 e 2 z d z 0 1 d x cos z x l a T α ( λ ) α ( 0 ) { x 2 2 | x 2 } r tm te ,
normalized to the E CP = 3 α ( 0 ) / 8 π a 4 —the Casimir-Polder (CP) energy for an ideal plane/atom at large distance a. Here, the prime means factor 1 / 2 for l = 0 , and we have to use λ = z x a , k = z a 1 x 2 for imaginary frequency ω = i λ and wave-vector k in conductivities η tm = 2 π σ tm and η te = 2 π σ te . The α is the polarizability of atom or molecule at the imaginary frequency, and the refraction coefficients of TE and TM modes are
r te = 1 1 + 1 x η te , r tm = 1 1 + x η tm .
The form of the zero terms, l = 0 , in (1) coincides exactly with that obtained for zero temperatures but, in general, with temperature and chemical potential dependence through the conductivities η tm | te = η tm | te ( λ , k , μ , m , T ) . We extract the zero, l = 0 , term
F = F tm + F te = F 0 + Δ F ,
and we consider the low-temperature expansion for Δ F and F 0 separately. We extract the temperature contribution from F 0 . Then,
F = F T = 0 + Δ F 0 + Δ F ,
where Δ F 0 = F 0 F T = 0 . Therefore, the total temperature correction, Δ T F , consists of two parts,
Δ T F = Δ F 0 + Δ F .
In general (for the graphene case, for example), the F T = 0 depends on the chemical potential.
The expansion crucially depends on the behavior of the conductivities at zero frequencies. In Ref. [23], we considered in detail the different models of conductivities with a constant value of conductivity at zero frequencies. In particular, we have taken into account the graphen’s conductivity with finite scattering factor γ , which means that the η tm | te | ω 0 = η tm | te | ω = γ . As a result, the free energy has the main low-temperature term T 2 for any relation between μ , m and T. With zero scattering factor, we have to consider this expansion more carefully.
In the framework of the polarization tensor approach [26], the conductivities of the TM and TE modes read [27,29]
η i = η i 0 + Δ η i , Δ η i = η gr m d y f i ( y ) Ξ ( y , μ , T ) ,
where
η te 0 η gr = 4 m π λ 1 + k F 2 4 m 2 2 m k F arctan k F 2 m , η tm 0 η te 0 = λ 2 k F 2 , f te ( y ) = 8 π λ Re ( 4 m 2 + q 2 ) ( q 2 k F 2 + 4 m 2 k 2 v F 2 ) q 2 k F 2 λ 2 r ( q 2 k F 2 + 4 m 2 k 2 v F 2 + q λ r ) , q = λ 2 i y , f tm ( y ) = 8 π Re q ( q 2 + k 2 v F 2 + 4 m 2 ) λ r r ( r + q λ ) , Ξ = 1 e y + μ T + 1 + 1 e y μ T + 1 , k F = λ 2 + v F 2 k 2 , r = k F 2 ( q 2 + k 2 v F 2 ) + 4 m 2 k 2 v F 2 ,
and η gr = 2 π σ gr c = π e 2 2 c = 0.0114 with σ gr = e 2 4 being the graphene universal conductivity.
Let us consider, for simplicity, the conductivity in the zero approximation over the Fermi velocity v F = 1 / 300 1 . In this approximation, we obtain more simple expressions ( x = tm , te )
η x 0 η gr = 4 m π λ 1 + λ 2 4 m 2 2 m λ arctan λ 2 m , f x ( y ) = 16 m 2 + y 2 π λ λ 2 + 4 y 2 ,
and we observe that the conductivities have no dependence on k, which should be the case because the Fermi velocity and wave-vector come in the single combination k v F .

2.1. Expansion of the Δ F

The sum with l 1 , the Δ F , may be represented in the following form [23]:
Δ F E CP = 8 3 Re l = 1 0 d z e i Λ l z Y tm + Y te ,
where
Y tm ( z ) = α ( λ ) α ( 0 ) z e 2 s s ( 2 s 2 z 2 ) s + z / η tm d s , Y te ( z ) = α ( λ ) α ( 0 ) z e 2 s z 3 z + s / η te d s ,
and Λ l = l / a T , λ = z / a , and k = s 2 z 2 / a . This representation is suitable for T 0 analysis, which means Λ l .
In the case v F 0 , the conductivities do not depend on k; therefore, we can calculate integral over s:
Y tm ( z ) = α z a e 2 z α ( 0 ) 2 η tm 3 η tm η tm 2 + η tm 2 2 η tm + 2 z 2 + η tm ( 2 η tm 1 ) z 2 η tm 2 2 z 3 e 2 1 + η tm 1 z Ei 2 z 1 + η tm 1 , Y te ( z ) = α z a α ( 0 ) η te z 3 e 2 η te z Ei 2 z ( 1 + η te ) ,
where Ei [ x ] is the exponential logarithm function. The function Ei [ x ] has the following representation as a series
Ei ( x ) = γ E + ln ( x ) + n 1 x n n · n ! ,
that contains polynomials, as well as logarithmic contributions.
Then, we make expansion over z,
Y x = n 0 A n x z n + n 1 B n x z n ln z ,
and use the Lemmas Erdélyi (These lemmas are sometimes called etalon integrals in the asymptotic methods of the stationary phase.) (see Ref. [32], Equations (1.13) and (1.35)), and Ref. [23]) to calculate asymptotic Λ l for integral over z in Equation (8) for each term of series. The manifest form of the coefficients depends on the specific model of conductivity. We obtain the series (12) in which we have to make replacements
z n n ! e i π 2 ( n + 1 ) Λ l n 1 , z n ln z n ! e i π 2 ( n + 1 ) i π 2 + ψ ( n + 1 ) ln Λ l Λ l n 1 ,
where ψ ( x ) is the digamma function, and then we take the real part (see Equation (8)) and obtain the following replacements
z 2 n + 1 ( 1 ) n + 1 ( 2 n + 1 ) ! Λ l 2 n 2 , z 2 n 0 , z 2 n + 1 ln z ( 1 ) n + 1 ( 2 n + 1 ) ! ( ψ ( n + 1 ) ln Λ l ) Λ l 2 n 2 , z 2 n ln z ( 1 ) n + 1 ( 2 n ) ! π 2 Λ l 2 n 1 .
Then, we make a summation over l 1 and arrive with relation
Δ F E CP = 8 3 X tm + X te ,
where
X x = n 0 A 2 n + 1 x ( 1 ) n + 1 ( 2 n + 1 ) ! ζ R ( 2 n + 2 ) ( a T ) 2 n + 2 + n 1 B 2 n x ( 1 ) n + 1 ( 2 n ) ! π 2 ζ R ( 2 n + 1 ) ( a T ) 2 n + 1 + n 0 B 2 n + 1 x ( 1 ) n + 1 ( 2 n + 1 ) ! ( [ ψ ( n + 1 ) + ln ( a T ) ] ζ R ( 2 n + 2 ) ζ R ( 2 n + 2 ) ) ( a T ) 2 n + 2 ,
and ζ R is the Riemann zeta function. We observe that the polynomial contributions come from odd A 2 n + 1 x and B n x . The logarithmic contributions come from odd B 2 n + 1 i . The main contribution reads
Δ F x E CP = 4 π 2 9 A 1 x B 1 x γ E + 6 ζ R ( 2 ) π 2 ( a T ) 2 + 8 π ζ R ( 3 ) 3 B 2 x ( a T ) 3 + 8 π 4 45 A 3 x + B 3 x 1 γ E 90 ζ R ( 4 ) π 4 ( a T ) 4 32 π ζ R ( 5 ) B 4 x ( a T ) 5 32 π 6 189 2 A 5 x + B 5 x 3 2 γ E 1890 ζ R ( 6 ) π 6 ( a T ) 6 + + ln ( a T ) 4 π 2 9 B 1 x ( a T ) 2 + 8 π 4 45 B 3 x ( a T ) 4 64 π 6 189 B 5 x ( a T ) 6 + ,
where γ E is the Euler constant. Note that, in general, the coefficients A n x and B n x depend on m , μ , and T.
For the constant conductivity case B 1 x = B 2 x = 0 and A 1 te = 0 , A 1 tm = 1 / 2 η tm . Therefore, the main contribution comes from TM mode and reads [23]
Δ F E CP = ( 2 π a T ) 2 9 A 1 tm = 2 π 2 ( a T ) 2 9 η tm .
For graphene case, the conductivities are expanded in the following series over z:
η x η gr = 1 z ( b 0 + z 2 b 2 + z 4 b 4 + ) ,
where the coefficients b n are functions of m , μ , and T:
b 0 = 4 a π m d y y 2 ( m 2 + y 2 ) Ξ , b 2 = 4 3 π m a 1 π a m d y y 4 ( m 2 + y 2 ) Ξ , b 4 = 2 15 π m 3 a 3 + 1 4 π a 3 m d y y 6 ( m 2 + y 2 ) Ξ .
The zero term b 0 crucially depends on m and μ for low temperatures:
I. μ > m , T μ m
b 0 = m μ g ( y ) d y + π 2 6 T 2 g ( μ ) + 7 π 4 360 T 4 g ( μ ) + O ( e μ m T ) .
II. μ = m , T m
b 0 = T ln 2 g ( m ) + π 2 12 T 2 g ( m ) + 3 4 ζ R ( 3 ) T 3 g ( m ) + 7 π 4 720 T 4 g ( m ) + O ( e m T ) .
III. μ < m , T m μ
b 0 = O ( e m μ T ) .
Here, g ( y ) = 4 a π m 2 + y 2 y 2 . These asymptotic are illustrated by numerical evaluation in Figure 1.
In the first case (19a) with μ > m , the expansion b 0 over T starts from the constant term (see Figure 1, right panel), and we obtain the following non-zero coefficients
A 3 tm = 1 3 , A 3 te = 1 , A 5 tm = α ( 0 ) 6 a 2 α ( 0 ) + 2 3 η gr b 0 3 η gr 2 b 0 2 2 15 , A 5 te = α ( 0 ) 2 a 2 α ( 0 ) 2 3 η gr b 0 1 3 η gr 2 b 0 2 2 3 , B 4 tm = 1 η gr b 0 , B 6 tm = α ( 0 ) 2 a 2 α ( 0 ) η gr b 0 + b 2 2 η gr 2 b 0 2 + 2 η gr 3 b 0 3 .
The main contribution to the b 0 comes from the first term of expansion in Equation (19a):
b 0 = 4 a π μ 2 m 2 μ .
In the last case (19c) with μ < m , the b 0 is exponentially small (see Figure 1, left panel), and we set it at zero, and the conductivities are expanded as the following:
η x η gr = z ( b 2 + z 2 b 4 + ) ,
and the firsts non-zero coefficients read
A 5 tm = 1 15 η gr b 2 ( η gr b 2 + 2 ) , A 5 te = 2 η gr b 2 , B 4 te = η gr b 2 .
For the specific case (19b), when μ = m , we consider the main contribution for T 0 and can set b 0 = 0 , as in the case (19c). Therefore, for all cases, we have the following expansion up to T 5 :
Δ F E CP = 16 π 4 135 ( a T ) 4 + 8 π 2 ζ R ( 3 ) μ a η gr ( μ 2 m 2 ) ( a T ) 5 + , ( μ > m , T μ m ) , Δ F E CP = 128 η gr ζ R ( 5 ) 3 a m ( a T ) 5 + , ( μ = m , T m ) , Δ F E CP = 128 η gr ζ R ( 5 ) 3 a m ( a T ) 5 + , ( μ < m , T m μ ) .

2.2. Expansion of the F 0

The zero term reads
F tm | te 0 E CP = 4 3 0 z 3 e 2 z d z 0 1 d x α ( z x a ) α ( 0 ) 2 x 2 1 + x η tm x 2 1 + 1 x η te .
The conductivities (6) depend on the temperature, and Δ η x have the following low-temperature expansions [23]
I. μ > m , T μ m
Δ η x η gr = m μ f x ( y ) d y + π 2 6 T 2 f x ( μ ) + 7 π 4 360 T 4 f x ( μ ) + O ( e μ m T ) .
II. μ = m , T m
Δ η x η gr = T ln 2 f x ( m ) + π 2 12 T 2 f x ( m ) + 3 4 ζ R ( 3 ) T 3 f x ( m ) + 7 π 4 720 T 4 f x ( m ) + O ( e m T ) .
III. μ < m , T m μ
Δ η x η gr = O ( e m μ T ) .
Taking these expansions into account, we obtain the following low-temperature corrections to zero term
Δ F x 0 E CP = ( a T ) 2 G x , ( μ > m , T μ m ) , Δ F x 0 E CP = ( a T ) H x , ( μ = m , T m ) , Δ F x 0 E CP = 0 , ( μ < m , T μ m ) ,
where
G tm | te = 2 π 2 9 a 2 0 z 3 e 2 z d z 0 1 d x α ( z x a ) α ( 0 ) ( 2 x 2 ) x f tm ( μ ) ( x + η tm 0 + Δ 0 η tm ) 2 x 3 f te ( μ ) ( 1 + x η te 0 + Δ 0 η te ) 2 , H tm | te = 4 ln 2 3 a 0 z 3 e 2 z d z 0 1 d x α ( z x a ) α ( 0 ) ( 2 x 2 ) x f tm ( m ) ( x + η tm 0 ) 2 x 3 f te ( m ) ( 1 + x η te 0 ) 2 .
The functions f x are given by Equation (6) and
Δ 0 η x = η gr m μ f x ( y ) d y .
Therefore, taking into account (24), we obtain the low-temperature expansion of the free energy
Δ T F E CP = ( a T ) 2 ( G tm + G te ) , ( μ > m , T μ m ) ,
Δ T F E CP = ( a T ) ( H tm + H te ) , ( μ = m , T m ) ,
Δ T F E CP = 128 η gr ζ R ( 5 ) 3 a m ( a T ) 5 + ( μ < m , T m μ ) .

3. Numerical Evaluations

We evaluated numerically the total temperature correction (5) by using the expression for the free energy in the form of Matsubara sum. We considered the Hydrogen atom at distance a = 100 nm from the graphene sheet. The polarizability of the Hydrogen atom in the single oscillator approximation may be found in Ref. [8], for example. The graphene conductivities are given by Equations (6). We used the Fermi velocity v F = 1 / 300 and the mass gap m = 0.1 eV. To visualize the dynamic of μ dependence, we made calculations for the value of μ close to m = 0.1 eV. The free energy is normalized to the E CP = 3 α ( 0 ) / 8 π a 4 —the CPenergy for an ideal plane/atom at large distance a.
The zero-temperature free energy, F T = 0 , depends on the chemical potential μ , and this dependence is shown in Figure 2. For μ m , it has the constant value, and it grows up starting with mass gap μ > m .
We proceed now to the consideration of the temperature correction, Δ T F to the free energy. First of all, let us consider the functions G x , which define low temperature expansion for μ > m case (30a). They are plotted in Figure 3. We observe that they are negative, and contribution from the TE mode is 100 times smaller.
One comment is in order. Because G x < 0 , then the temperature correction Δ T F > 0 close to the T = 0 . At the same time, in Refs. [7,24], the positive value of correction was observed. The disagreement is connected with that in which we use distance a = 100 nm for numerical evaluations, which is out of distances considered in Refs. [7,24].
The numerical evaluation of H x ( a ) for a = 100 nm gives the following values H tm = 3 × 10 6 and H te = 0.18 . Again, the main contribution comes from TM mode. The value of H tm strongly depends on the value of the Fermi velocity. For v F 0 , the H tm .
The temperature contribution for μ m is shown in Figure 4. We observe that, for low temperatures, the free energy has the form of parabola, ∼ G tm ( a T ) 2 , in agreement with (30) with negative parameter G tm (see Figure 3).
The closer μ > m to m, the smaller domain of temperature T where this approximation is valid, and the greater value of the parameter of parabola G tm , in agreement with Figure 3. If μ = m , this domain becomes zero, and the free energy drastically changes its form. If μ 0 , the part of the curve which is out of this domain (the vertical part of the green curve, for example) goes to free energy for this very special position with μ = m (black curve). Therefore, for any infinitely small difference μ m 0 , the derivative of free energy with respect temperature T is zero for T = 0 , and the Nernst theorem is valid. The experimental realization of the exact equality μ = m cannot be realized, and we conclude that the Nernst theorem is valid for this system.
The temperature contribution for μ m is shown in Figure 5. We observe the completely different dependence of the energy on the μ .
For zero chemical potential (brown curve), the temperature correction is, in fact, zero for the large domain of temperatures. The closer μ to m, the smaller domain in which temperature correction is zero. According to (30b),
Δ T F E CP = 128 η gr ζ R ( 5 ) 3 a m ( a T ) 5 = 10 24 T 5 ( K ) ,
in this domain.
Figure 6 shows the temperature contribution to the free energy as a function of chemical potential. The function has a very sharp form with the maximum for μ = m with a different slope at the left and the right of this point.
The point μ = m looks like a phase transition point between different regimes from Δ T F T 4 to Δ T F T 2 . In fact, it is an unstable point—infinitely small deviation μ from m changes regime.
From Figure 4 and Figure 5 and relations (30), we observe the different signs of the entropy, S, for μ < m and μ > m . The entropy is the negative derivative of the free energy with respect of temperature. Therefore, S μ < m < 0 (see Figure 5), and S μ > m > 0 (see Figure 4). In both cases, the Nernst theorem is valid, S T 0 0 . The negative entropy of the dispersion forces has already been observed in Refs. [19,33] for plain and spherical configurations in the framework of plasma model and was also discussed recently in Ref. [34].
The domains of temperatures where the low-temperature expansions (30) over dimensionless parameter ( a T ) 1 are valid depend on the chemical potential μ , mass gap m, and the distance a between an atom and plain. Let us make some numerical estimations. The mass gap m = 0.1 eV corresponds to the temperature T m = 1160 K. Therefore, the restriction T | μ m | in Equations (30a) and (30c) plays the role for | μ m | < 0.1 eV, in agreement with Figure 4 and Figure 5. The restriction in the second case (30b) plays the role if T m < 0.1 eV. The more strong limitation appears due to the huge value of parameter H tm 10 6 (for a = 100 nm): T ( H tm a ) 1 = 0.14 K for a = 100 nm. It means that, in this case, there is no domain of temperatures where this case is valid, according to above discussion. The restriction ( a T ) 1 is important starting from a = 10 μm. In this case, the effective temperature T a = 1 / a = 1438 K.

4. Conclusions

We considered the low-temperature correction to the Casimir-Polder free energy for atom/graphene system by using the Poisson representation of the free energy, which is more suitable for low-temperature analysis. The analysis is naturally broken into three different regions, (i) μ > m , (ii) μ = m , and (iii) μ < m , for chemical potential. This division is the consequence of the same regions for the conductivity of graphene (see Equation (26). The conductivities have completely different expansion in these regions. It starts from the constant in the first region, linear on the temperature in the second one, and exponentially small in the third region.
The free energy may be divided into the two parts (3). The first one, F 0 , has the form of the free energy at zero temperature but with μ , m , and T dependence via the conductivity dependence on these parameters. The main contribution in the low-temperature expansion in the first (i) and second (ii) regimes comes from this first term, and it is quadratic and linear on the temperature, correspondingly (see Equation (30). In the third (iii) regime, the main contribution T 5 comes from the rest part Δ F .

Author Contributions

Investigation, both authors; writing, both authors; numerical calculations, both authors. All authors have read and agreed to the published version of the manuscript.

Funding

The NK was supported in part by the grants 2019/10719-9, 2016/03319-6 of São Paulo Research Foundation (FAPESP) and by the Russian Foundation for Basic Research Grant No. 19-02-00496-a.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All necessary data are contained in this paper.

Acknowledgments

We are grateful to Dmitri Vassilevich, Galina Klimchitskaya, and Vladimir Mostepanenko for fruitful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Casimir, H.B.G. On the attraction between two perfectly conducting plates. Kon. Ned. Akad. Wetensch. Proc. 1948, 51, 793–795. [Google Scholar]
  2. Casimir, H.B.G.; Polder, D. The influence of retardation on the London-van der Waals forces. Phys. Rev. 1948, 73, 360–372. [Google Scholar] [CrossRef]
  3. Parsegian, A.V. Van der Waals Forces. A Handbook for Biologists, Chemists, Engineers, and Physicists; Cambridge University Press: Cambridge, UK, 2006; p. 380. [Google Scholar]
  4. Bordag, M.; Klimchitskaya, G.L.; Mohideen, U.; Mostepanenko, V.M. Advances in the Casimir Effect; Cambridge University Press: Cambridge, UK, 2009; pp. 1–768. [Google Scholar] [CrossRef]
  5. Bondarev, I.V.; Lambin, P. Van der Waals coupling in atomically doped carbon nanotubes. Phys. Rev. B 2005, 72, 35451. [Google Scholar] [CrossRef] [Green Version]
  6. Blagov, E.V.; Klimchitskaya, G.L.; Mostepanenko, V.M. Van der Waals interaction between a microparticle and a single-wall carbon nanotube. Phys. Rev. B 2007, 75, 235413. [Google Scholar] [CrossRef] [Green Version]
  7. Klimchitskaya, G.L.; Mostepanenko, V.M. Casimir and Casimir-Polder Forces in Graphene Systems: Quantum Field Theoretical Description and Thermodynamics. Universe 2020, 6, 150. [Google Scholar] [CrossRef]
  8. Khusnutdinov, N.; Kashapov, R.; Woods, L.M. Casimir-Polder effect for a stack of conductive planes. Phys. Rev. A 2016, 94, 12513. [Google Scholar] [CrossRef] [Green Version]
  9. Khusnutdinov, N.; Woods, L.M. Casimir Effects in 2D Dirac Materials (Mini-review). JETP Lett. 2019, 110, 1–10. [Google Scholar] [CrossRef]
  10. Babb, J.F. Long-range atom-surface interactions for cold atoms. J. Phys. Conf. Ser. 2005, 19, 1–9. [Google Scholar] [CrossRef] [Green Version]
  11. Marachevsky, V.N.; Pis’mak, Y.M. Casimir-Polder effect for a plane with Chern-Simons interaction. Phys. Rev. D 2010, 81, 65005. [Google Scholar] [CrossRef] [Green Version]
  12. Shajesh, K.V.; Schaden, M. Repulsive long-range forces between anisotropic atoms and dielectrics. Phys. Rev. A 2012, 85, 012523. [Google Scholar] [CrossRef] [Green Version]
  13. Thiyam, P.; Parashar, P.; Shajesh, K.V.; Persson, C.; Schaden, M.; Brevik, I.; Parsons, D.F.; Milton, K.A.; Malyi, O.I.; Boström, M. Anisotropic contribution to the van der Waals and the Casimir-Polder energies for CO2 and CH4 molecules near surfaces and thin films. Phys. Rev. A 2015, 92, 052704. [Google Scholar] [CrossRef] [Green Version]
  14. Antezza, M.; Fialkovsky, I.; Khusnutdinov, N. Casimir-Polder force and torque for anisotropic molecules close to conducting planes and their effect on CO2. Phys. Rev. B 2020, 102, 195422. [Google Scholar] [CrossRef]
  15. Obrecht, J.M.; Wild, R.J.; Antezza, M.; Pitaevskii, L.P.; Stringari, S.; Cornell, E.A. Measurement of the Temperature Dependence of the Casimir-Polder Force. Phys. Rev. Lett. 2007, 98, 063201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Laliotis, A.; de Silans, T.P.; Maurin, I.; Ducloy, M.; Bloch, D. Casimir-Polder interactions in the presence of thermally excited surface modes. Nat. Commun. 2014, 5. [Google Scholar] [CrossRef] [Green Version]
  17. Wylie, J.M.; Sipe, J.E. Quantum electrodynamics near an interface. II. Phys. Rev. A 1985, 32, 2030–2043. [Google Scholar] [CrossRef] [PubMed]
  18. Buhmann, S.; Welsch, D. Dispersion forces in macroscopic quantum electrodynamics. Prog. Quantum Electron. 2007, 31, 51–130. [Google Scholar] [CrossRef] [Green Version]
  19. Bezerra, V.B.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Romero, C. Lifshitz theory of atom-wall interaction with applications to quantum reflection. Phys. Rev. A 2008, 78, 042901. [Google Scholar] [CrossRef] [Green Version]
  20. Chaichian, M.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Tureanu, A. Thermal Casimir-Polder interaction of different atoms with graphene. Phys. Rev. A 2012, 86, 12515. [Google Scholar] [CrossRef] [Green Version]
  21. Bordag, M. Low Temperature Expansion in the Lifshitz Formula. Adv. Math. Phys. 2014, 2014, 1–34. [Google Scholar] [CrossRef]
  22. Khusnutdinov, N.; Kashapov, R.; Woods, L.M. Thermal Casimir and Casimir–Polder interactions in N parallel 2D Dirac materials. 2D Mater. 2018, 5, 35032. [Google Scholar] [CrossRef] [Green Version]
  23. Khusnutdinov, N.; Emelianova, N. Low-temperature expansion of the Casimir-Polder free energy for an atom interacting with a conductive plane. Int. J. Mod. Phys. A 2019, 34, 1950008. [Google Scholar] [CrossRef] [Green Version]
  24. Klimchitskaya, G.L.; Mostepanenko, V.M. Nernst heat theorem for an atom interacting with graphene: Dirac model with nonzero energy gap and chemical potential. Phys. Rev. D 2020, 101, 116003. [Google Scholar] [CrossRef]
  25. Klimchitskaya, G.L.; Mostepanenko, V.M. Quantum field theoretical description of the Casimir effect between two real graphene sheets and thermodynamics. Phys. Rev. D 2020, 102, 016006. [Google Scholar] [CrossRef]
  26. Bordag, M.; Fialkovsky, I.V.; Gitman, D.M.; Vassilevich, D.V. Casimir interaction between a perfect conductor and graphene described by the Dirac model. Phys. Rev. B 2009, 80, 245406. [Google Scholar] [CrossRef] [Green Version]
  27. Fialkovsky, I.V.; Marachevsky, V.N.; Vassilevich, D.V. Finite-temperature Casimir effect for graphene. Phys. Rev. B 2011, 84, 35446. [Google Scholar] [CrossRef] [Green Version]
  28. Bordag, M.; Klimchitskaya, G.L.; Mostepanenko, V.M.; Petrov, V.M. Quantum field theoretical description for the reflectivity of graphene. Phys. Rev. D 2015, 91, 045037. [Google Scholar] [CrossRef] [Green Version]
  29. Bordag, M.; Fialkovskiy, I.; Vassilevich, D. Enhanced Casimir effect for doped graphene. Phys. Rev. B 2016, 93, 075414, Erratum in Phys. Rev. B 2017, 95, 119905. doi:10.1103/physrevb.95.119905. [Google Scholar] [CrossRef] [Green Version]
  30. Falkovsky, L.A.; Varlamov, A.A. Space-time dispersion of graphene conductivity. Eur. Phys. J. B 2007, 56, 281–284. [Google Scholar] [CrossRef] [Green Version]
  31. Gusynin, V.P.; Sharapov, S.G.; Carbotte, J.P. Magneto-optical conductivity in Graphene. J. Phys. Condens. Matter 2007, 19, 26222. [Google Scholar] [CrossRef] [Green Version]
  32. Fedoryuk, M.V. The Saddle-Point Method; Nauka: Moscow, Russia, 1977. (In Russian) [Google Scholar]
  33. Khusnutdinov, N.R. The thermal Casimir–Polder interaction of an atom with a spherical plasma shell. J. Phys. A Math. Theor. 2012, 45, 265301. [Google Scholar] [CrossRef] [Green Version]
  34. Li, Y.; Milton, K.; Parashar, P.; Hong, L. Negativity of the Casimir Self-Entropy in Spherical Geometries. Entropy 2021, 23, 214. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The functions b 0 for a = 100 nm, m = 0.1 eV, and different values of chemical potential μ . For μ < m , the function b 0 0 for low temperatures; for μ = m , it is linear, and, for μ > m , it starts from constant values in agreement with Equation (19).
Figure 1. The functions b 0 for a = 100 nm, m = 0.1 eV, and different values of chemical potential μ . For μ < m , the function b 0 0 for low temperatures; for μ = m , it is linear, and, for μ > m , it starts from constant values in agreement with Equation (19).
Universe 07 00070 g001
Figure 2. The zero temperature free energy, F T = 0 , as a function of chemical potential μ .
Figure 2. The zero temperature free energy, F T = 0 , as a function of chemical potential μ .
Universe 07 00070 g002
Figure 3. The functions G x (30a) for different values of the chemical potential μ > m .
Figure 3. The functions G x (30a) for different values of the chemical potential μ > m .
Universe 07 00070 g003
Figure 4. The temperature contribution to the free energy in the different intervals of temperatures and μ m .
Figure 4. The temperature contribution to the free energy in the different intervals of temperatures and μ m .
Universe 07 00070 g004
Figure 5. The temperature contribution to the free energy in the different intervals of temperatures and μ m .
Figure 5. The temperature contribution to the free energy in the different intervals of temperatures and μ m .
Universe 07 00070 g005
Figure 6. The temperature contribution to the free energy as a function of the chemical potential for distance between Hydrogen atom and graphene a = 100 nm.
Figure 6. The temperature contribution to the free energy as a function of the chemical potential for distance between Hydrogen atom and graphene a = 100 nm.
Universe 07 00070 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khusnutdinov, N.; Emelianova, N. The Low-Temperature Expansion of the Casimir-Polder Free Energy of an Atom with Graphene. Universe 2021, 7, 70. https://doi.org/10.3390/universe7030070

AMA Style

Khusnutdinov N, Emelianova N. The Low-Temperature Expansion of the Casimir-Polder Free Energy of an Atom with Graphene. Universe. 2021; 7(3):70. https://doi.org/10.3390/universe7030070

Chicago/Turabian Style

Khusnutdinov, Nail, and Natalia Emelianova. 2021. "The Low-Temperature Expansion of the Casimir-Polder Free Energy of an Atom with Graphene" Universe 7, no. 3: 70. https://doi.org/10.3390/universe7030070

APA Style

Khusnutdinov, N., & Emelianova, N. (2021). The Low-Temperature Expansion of the Casimir-Polder Free Energy of an Atom with Graphene. Universe, 7(3), 70. https://doi.org/10.3390/universe7030070

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