1. Introduction
The Casimir–Polder force is a typical example of dispersion interactions caused by a fluctuating electromagnetic field which is zero in the mean but possesses a nonzero dispersion. Discovered by Casimir and Polder [
1] as a relativistic generalization of the van der Waals interaction between two atoms (molecules) or an atom and a material surface, it has received wide recognition in both fundamental physics dealing with precision experiments and in nanotechnology (for more, see the monographs [
2,
3,
4,
5,
6,
7]).
The modern theoretical description of the Casimir–Polder force is based on quantum electrodynamics. In the original article [
1], Casimir and Polder considered the interaction of two neutral atoms and an atom with an ideal metal plate. The case of an atom spaced at large separation from a plate made of real material characterized by the frequency-dependent dielectric permittivity was considered in [
8]. This was done on the basis of the Lifshitz theory [
9] describing the Casimir force acting between two thick plates by rarefying the material of one of them.
A more general formula for the potential energy of a molecule of solute in a dilute solution of low concentration as a function of the interface distance was obtained in [
10] using the same rarefying procedure of one of the plates (see also [
11]). At a later time, the Lifshitz formula for the Casimir–Polder interaction of microparticles with material plates was used in several hundreds papers (see [
2,
3,
4,
5,
6,
7] and references therein). This formula allows calculation of the Casimir–Polder interaction given the dynamic atomic (molecular) polarizability and the dielectric permittivity of plate material. Formulas of this kind are applicable in the case of nanoparticles interacting with planar or gently curved surfaces under the condition that the nanoparticle radius
R is much less than the distance to the surface [
12,
13,
14,
15].
Among the new materials which have appeared, the one-atom-thick layer of carbon atoms called graphene seems to hold the lead [
16]. Thanks to its two-dimensional structure, graphene admits a theoretical description in the framework of the Dirac formalism of quantum electrodynamics [
17,
18] without resorting to the phenomenological assumptions and far-reaching extrapolations which are often used in condensed matter physics when dealing with ordinary three-dimensional materials. This enables the response of graphene to the electromagnetic field to be found as a function of frequency, wave vector, and temperature with the use of the polarization tensor [
19,
20,
21,
22]. In actual fact, the two independent components of the polarization tensor are equivalent to the spatially nonlocal dielectric permittivities of graphene [
23,
24]. Because of this, by using this formalism one can calculate the Casimir–Polder force acting from the source side of graphene on atoms, molecules, and nanoparticles.
Calculations of the Casimir–Polder force from graphene acting on atoms (see, e.g., [
25,
26,
27,
28,
29,
30,
31,
32,
33,
34]) and nanoparticles (see, e.g., [
35,
36,
37,
38,
39,
40,
41,
42]) find applications in both fundamental physics and nanotechnology. An expression for the force takes an especially simple form in the limiting case of large separations (high temperatures). At room temperature
T = 300 K, this is reached at a few micrometers. In [
43], the large separation behavior of the Casimir–Polder force was investigated for a graphene sheet with no foreign atoms possessing zero chemical potential and sufficiently small energy gap in the spectrum of quasiparticles.
Real graphene sheets are characterized by nonzero values of both the energy gap
and chemical potential
. In a recent article [
44], the large-separation behavior of the Casimir–Polder force was explored for a real graphene sheet freestanding in vacuum as a function of
and
. It was shown that it is possible to control the force value by varying the values of
and
and deal with the analytic asymptotic expressions for the force rather than with more complicated computations by means of the Lifshitz formula.
In practical implementations, graphene is usually deposited on a dielectric substrate. Here, we investigate the large-separation behavior of the Casimir–Polder force acting on atoms or nanoparticles from real graphene sheets characterized by an energy gap and chemical potential which are deposited on thick silica glass substrate. We find that the range of large separations for a graphene-coated substrate starts at approximately 5.6 m from the surface and is almost independent of the values of and . The ratio of forces from a graphene-coated substrate and an uncoated silica glass plate decreases with increasing and increases with increasing . It is shown that the classical limit of the Casimir–Polder force from an ideal metal plane reached at the thermal length, for a graphene-coated substrate may be reached at much larger separations from the surface which essentially depend on the values of and .
The analytic asymptotic expression for the large-separation behavior of the Casimir–Polder force from a graphene-coated substrate is found as a function of separation from the surface, energy gap, chemical potential, and temperature. The asymptotic results are compared with numerical computations of the Casimir–Polder force at large separations for different values of and . According to our results, at some fixed separation an agreement between the asymptotic and numerical results becomes better with increasing chemical potential and worse with increasing energy gap. By and large, for a graphene-coated substrate the asymptotic results are found to be in slightly better agreement with the results of numerical computations than for a freestanding graphene sheet.
The structure of the article is as follows.
Section 2 is devoted to the general formalism of the Lifshitz theory for graphene deposited on a substrate. In
Section 3, we consider the impact of a substrate on the large-separation Casimir–Polder force from real graphene sheet.
Section 4 is devoted to a confrontation between the analytic asymptotic results and numerical computations in the presence of a substrate.
Section 5 provides a discussion, while
Section 6 lists our conclusions.
2. General Formalism
Here, we briefly present the formalism of the Lifshitz theory describing the Casimir–Polder force between an atom or a nanoparticle and any planar structure, e.g., a graphene-coated substrate, which are at temperature
T and separated by a distance
a. Using the dimensionless Matsubara frequencies
, where
is the Boltzmann constant,
, and
with
,
being the projection of the wave vector on the planar structure, the Lifshitz formula for the Casimir–Polder force takes the following form [
6]:
In this equation, the prime on the summation sign divides the term with
by 2 and
is the dynamic polarizability of an atom, a molecule, or a nanoparticle taken at a pure imaginary frequency that has the dimension of
.
Particular attention should be paid to the quantities
and
in (
1); these are the reflection coefficients for the transverse magnetic (TM) and transverse electric (TE), or equivalently p and s, polarized electromagnetic waves on the planar structure calculated at the pure imaginary frequencies
In our case, the planar structure is a graphene sheet deposited on a dielectric substrate. The electromagnetic response of the graphene is described by the polarization tensor in (2 + 1) dimensions, whereas the material of the substrate reacts to the electromagnetic field through its frequency-dependent dielectric permittivity. The form of reflection coefficients in this unusual case was found in [
45]. Using the dimensionless variables introduced above, these coefficients are provided by [
29]
Here, in analogy with (
2), the dielectric permittivity of a substrate is taken at the pure imaginary frequencies
The quantity
is the 00 component and
is the linear combination of components of the dimensionless polarization tensor
which are taken at the pure imaginary frequencies
.
The explicit expressions for the quantities
and
for real graphene sheets characterized by any value of the energy gap
and chemical potential
are presented in [
29]. Here, we use them in the more convenient identical form provided in [
44]. Thus, the 00 component is
The quantity
in (
7) is the fine structure constant
It should not be confused with the static polarizability
used below. The function
is defined as
where
and
is the dimensionless Fermi velocity,
.
As to the function
entering (
7), it is provided by
where
Using the same notations, we express the linear combination of the components of the polarization tensor
from Equation (
4) defining the TE reflection coefficient [
29,
44]:
As is seen from (
1), (
4)–(
13), for computation of the Casimir–Polder force acting on some particle on the source side of a graphene-coated substrate, it is desirable to know the dynamic polarizability of this particle (
2) and the dielectric permittivity of a substrate (
5) within a sufficiently wide frequency region, along with the energy gap, the chemical potential of the graphene sheet, and the temperature.
It is pertinent to note that the above formalism was developed in the framework of the Dirac model of graphene [
16,
17,
18], which applies when the energy is not too high, i.e.,
eV [
46]. Taking into consideration that the characteristic frequency of the Casimir–Polder force is provided by
, which leads to energies of less than 1 eV at all separations exceeding 100 nm [
6], it can be concluded that this formalism is very well adapted for theoretical description of the Casimir–Polder force at large separations.
3. Impact of Substrate on the Casimir–Polder Force from a Real Graphene Sheet
at Large Separations
For a particle interacting with a lonely graphene sheet which is freestanding in vacuum, it was shown in [
44] that, for separations exceeding some value
from 2.3 to
m, 99% of the total Casimir–Polder force is provided by the term of (
1) with
. This result is valid for various values of the energy gap and chemical potential of graphene. We start by finding how it is modified by the presence of a substrate.
For this purpose, we consider the term of (
1) with
where
is the static polarizability of an atom or a nanoparticle. It is seen that
does not depend on the TE reflection coefficient.
The reflection coefficient
is obtained from (
4), again by setting
,
When obtaining (
15), it was assumed that
which is valid for substrates made of both dielectric materials and metals described by the Drude model. The value of
in (
15) is provided by (
7) with
[
29]
Here, the quantity
is obtained from (
10) and
entering
is obtained from (
12)
Note that Equation (
14) of [
44] for
contains a typo: the factor
in the first term on the right-hand side of (14) in [
44] should be replaced with
y.
All numerical computations below are performed for a graphene sheet deposited on a fused silica glass (SiO
) substrate. This substrate is typical in graphene technologies (see, e.g., [
47,
48,
49,
50]) and in precision measurements of the Casimir force from graphene [
51,
52]. The optical data for the complex index of refraction of fused silica glass can be found in [
53] in the wide range of
extending from 0.0025 eV to 2000 eV. The obtained dielectric permittivity of SiO
along the imaginary frequency axis is well known and has been used in many publications (see, e.g., [
6]). It contains the two steps, one of which is due to ionic polarization and another one is due to electronic polarization. The static dielectric permittivity of SiO
is
.
Now, we compute the large-separation Casimir–Polder force from a graphene-coated SiO
substrate using the full Lifshitz formula (
1) and its term with
(
14) for various values of the energy gap and chemical potential. Taking into account that in the presence of the substrate the region of large separations starts at separations exceeding a few micrometers, full computations using (
1) can be safely performed by setting
without sacrificing precision [
6]. When using the Lifshitz formula (
1), the computations are performed by (
4)–(
13). When Equation (
14) is employed, we use (
15), (
17) and (
18) in the computations. In both cases, we use the value of
for the SiO
substrate, and in full computations we use values of
with
.
Comparing the obtained results, we find that the full values of the Casimir–Polder force
computed by (
1) differ from the force
computed by (
14) in less than 1% at all separations
m independently of the values of
and
used in computations (see below for the typical specific values of these parameters). Thus, for a graphene sheet deposited on a substrate, the large-separation behavior of the Casimir–Polder force provided by
is reached at larger separations than for a freestanding graphene sheet (i.e., from 2.3 to 3.2
m, depending on the values of
and
[
44]). What is more, for a graphene-coated substrate the value of
is essentially independent on
and
. Note that for an uncoated silica glass plate the large-separation behavior of the Casimir–Polder force is reached at
m.
In an effort to determine the impact of graphene coating on the large-separation Casimir–Polder force, we computed the force
from a graphene-coated SiO
substrate and
from an uncoated fused silica plate. To compute the latter quantity, one should set
in (
15). In
Figure 1a,b the computational results for the ratio
at
T = 300 K,
m are shown as a function of the energy gap by (a) the four lines labeled 1, 2, 3, and 4 for
, 0.05, 0.1, and 0.15 eV, respectively, and (b) the three lines labeled 4, 5, and 6 for
, 0.2, and 0.25 eV, respectively. For illustrative purposes, the dashed lines show the value of the ratio
when the graphene-coated substrate is replaced with an ideal metal plane.
As seen in
Figure 1a,b, the impact of graphene coating on the large-separation Casimir–Polder force increases with increasing chemical potential and decreasing energy gap of the graphene. As a result, the fused silica substrate coated by a graphene sheet with
eV and
eV produces almost the same Casimir–Polder force at
m separation as an ideal metal plane.
It is well known that at separations above 6
m the Casimir–Polder force from an ideal metal plane takes the so-called classical form, which does not depend on either
ℏ or
c [
6]:
In order to find how the large-separation Casimir–Polder force from a graphene-coated SiO
substrate approaches this classical limit depending on the values of
and
, we can calculate the relative difference
In
Figure 2, the computational results for
at
T = 300 K are shown for a graphene coating with
eV as a function of separation by the five lines labeled 1, 2, 3, 4, and 5 for
, 0.025, 0.05, 0.075, and 0.1 eV, respectively. The dashed line shows the lower boundary of the figure plane domain where the relative deviation between
and
is less than 1%.
As is seen in
Figure 2, for a graphene coating with
eV (line 5), the relative deviation (
20) remains within 1% over the entire separation region
m. As for the graphene coating with smaller
, 0.05, 0.025, and 0 eV, the Casimir–Polder force takes the classical form at larger separations equal to 6.5, 15.5, 35.5, and 53
m, respectively. These values are only slightly larger than for a freestanding graphene sheet [
44]. From
Figure 2, it can be concluded that the Casimir–Polder force from a graphene-coated substrate takes the classical form at larger separations for graphene with lower chemical potential.
In fact, for a graphene sheet deposited on a substrate, the initiation of an energy gap is almost unavoidable [
16]. Thus, an energy gap equal to 0.3 eV is rather typical [
51,
52]. To cover this case, in
Figure 3a,b we present the computational results for
at
T = 300 K for a graphene coating with
eV. These results are again shown as the function of separation (a) by the five lines labeled 1, 2, 3, 4, and 5 for
, 0.025, 0.05, 0.075, and 0.1 eV, respectively, over the separation interval from 5.6 to 60
m and (b) by the four lines labeled 1, 2, 3, and 4 for
, 0.025, 0.05, and 0.075 eV, respectively, over the separation interval from 60 to 200
m. The dashed lines in
Figure 3a,b again enclose the figure domain where
is less than 1%.
As is seen in
Figure 3a, only for a graphene coating with
eV (line 5) does the Casimir–Polder force take the classical form in the field of this figure (for
m). Thus, with increasing energy gap the classical form of the force is reached at larger separations. From
Figure 3b, it can be seen that for graphene sheets with
and 0.05 eV the Casimir–Polder force becomes classical starting from 63 and 157
m, respectively (lines 4 and 3). As for the graphene coating with
and 0 eV, the corresponding Casimir–Polder force does not take the classical form in the field of this figure up to
m. Calculations show that the Casimir–Polder force from the graphene-coated substrate with
eV,
and 0 eV becomes classical only at separations of 363 and 550
m, respectively. By and large,
Figure 2 and
Figure 3 demonstrate that with increasing
the classical form of the Casimir–Polder force from the graphene-coated substrate is reached at shorter and with increasing
at larger separations.
4. Analytic Asymptotic Results Confronted with Numerical Computations in the
Presence of Substrate
In this section, we obtain the analytic asymptotic expressions for the large-separation Casimir–Polder force from a graphene-coated substrate and find the measure of agreement between and the numerically computed for different values of the energy gap and chemical potential of a graphene sheet.
The asymptotic expression sought here is valid under the following condition:
which is satisfied with a large safety margin at
K,
m (in fact, as was shown in
Section 3, we consider separation distances exceeding
m).
We consider the large-separation Casimir–Polder force (
14) with the reflection coefficient (
15) under the condition (
21). This reflection coefficient can be rearranged to
According to (
17), the polarization tensor is of the order of parameter (
21); thus,
. What is more, the main contribution to the integral (
14) is provided by
. Because of these conditions, we can replace
y with unity in the denominator of (
22) and neglect
in comparison with
. Therefore, (
22) becomes
i.e., the TM reflection coefficient takes exactly the same approximate form as was found earlier for the freestanding graphene sheet [
44].
After substitution of (
23) in (
14) and integration, the desired asymptotic expression takes the same form as in [
44]
where
is defined in (
19).
The expression (
24) should be supplemented by the approximate expression for
found in [
44] under the condition (
21)
We are coming now to a detailed comparison between the asymptotic expression for the Casimir–Polder force
provided by (
24) and (
25) and the numerically computed large-separation force
from a graphene-coated substrate using Equations (
14), (
15) and (
17). This comparison can be made by considering the ratio
for different values of the graphene parameters. All computations are again performed at room temperature
K.
First, we consider an undoped graphene coating with
and compute the ratio
for two different values of the energy gap
eV and 0.2 eV. The computational results are shown in
Figure 4 as the function of separation by the top and bottom solid lines, respectively. For comparison purposes, the ratio
for a freestanding graphene sheet (i.e., with no substrate) with the same values of
is plotted by the dashed lines.
As seen in
Figure 4, for a graphene coating with
,
eV, the asymptotic results reproduce the results of the numerical computations with better than 1% accuracy at separations exceeding
m. With increasing energy gap, the accuracy of the asymptotic expressions becomes lower. Thus, for a graphene coating with
,
eV, the asymptotic expression is accurate within 1% starting from
m separation. From
Figure 4, it can be seen that in spite of the fact that the asymptotic expression used does not depend on the dielectric permittivity of a substrate it reproduces the results of numerical computations somewhat better than for a freestanding graphene sheet. This is illustrated by the solid lines, which lie slightly above the dashed ones.
Next, we consider a graphene coating with a reasonably large energy gap
eV and consider the relationship between the asymptotic and numerical results for a large-separation Casimir–Polder force from a graphene-coated substrate for different values of the chemical potential. The computational results for the ratio
in this case at
T = 300 K are presented in
Figure 5 by the three solid lines from bottom to top for the chemical potential equal to
, 0.05, and 0.075 eV, respectively. The dashed lines show the ratio
for the freestanding graphene sheets with the same values of
. In the inset, the two pairs of solid and dashed lines for
and 0.075 eV are reproduced on an enlarged scale within a narrower separation region for better visualization.
From
Figure 5, it can be seen that with increasing
the agreement between the asymptotic and numerically computed Casimir–Polder forces from the graphene-coated substrate becomes better. Thus, if for a graphene coating with
eV the 1% agreement is reached only at separations exceeding 41
m, for
eV this measure of agreement is observed at
m. As to the graphene coating with
eV, the better than 1% agreement between the asymptotic and numerical results holds over the entire range of large separations
m. According to
Figure 5 (see the inset as well), the solid lines lie above the dashed ones for all
, i.e., the asymptotic expression (
24) and (
25) is slightly more exact in the case of graphene-coated substrates than for freestanding graphene sheets.
At the end of this section, we consider the case of a graphene coating with a larger energy gap
eV, such as that used in the experiment measuring the Casimir force from a graphene-coated substrate in [
51,
52]. Here, we calculate the ratio
over the wider range of
up to
eV (the latter value was measured for the sample used in the experiment [
51,
52]). In
Figure 6a,b, we present the computational results for this ratio at
K as a function of separation. The obtained results are shown by (a) the five lines labeled 1, 2, 3, 4, and 5 for
, 0.025, 0.05, 0.075, and 0.1 eV, respectively, over the separation interval from 5.6
m to 100
m and (b) by the three lines labeled 6, 7, and 8 for
, 0.2, and 0.25 eV, respectively, over the interval from 5.6 to 30
m.
As seen in
Figure 6a, for lines 1, 2, and 3 (i.e., for the graphene coatings with
, 0.025, and 0.05 eV, respectively) the 1% agreement between the asymptotic and numerical results is not reached up to the separation of 100
m. As for the lines 4 and 5 (
, and 0.1 eV, respectively) the 1% agreement is reached at separation distances exceeding approximately 60 and 24
m. Thus, the agreement between the asymptotic and numerical results again becomes better with increasing value of the chemical potential.
From
Figure 6b, an inference can be drawn that sufficiently large values of
can fully compensate the negative role played by large
in an agreement between the asymptotic and numerical values of the large-separation Casimir–Polder force from graphene-coated substrates. According to
Figure 6b, for all the three graphene coatings with
, 0.2, and 0.25 eV the asymptotic results are within 1% agreement with the results of numerical computations at all separations exceeding the border of the large-separation region equal to 5.6
m. These results make possible the reliable use of the analytic asymptotic expression for the Casimir–Polder force from graphene-coated substrates with the proper combination of the values of
and
.