1. Introduction
Over the past two decades, much effort has been spent investigating the static [
1,
2] and dynamic [
3,
4] Casimir effect in various geometric configurations, including a system of two parallel metal (dielectric) plates separated by a narrow vacuum gap. The main objectives of these studies are the properties of a fluctuating electromagnetic field and its interaction with matter on the nanoscale. Measurement of these effects paves the way to the core of nonequilibrium quantum field theory [
5,
6,
7,
8].
In addition to the attractive (in most cases) Casimir forces between electrically neutral bodies at rest, a dissipative tangential force arises when one or both bodies move relative to each other. The sources of these forces are dissipative effects within the plates (Joule losses). In this case, the corresponding fluctuation electromagnetic forces are called “van der Waals” [
9], “Casimir” [
10], or “quantum” [
11] forces of friction. As it looks, it is convenient to use the general name “Casimir-Lifshitz” (CL) friction force, which incorporates all the features of these dissipative forces regarding their distance, temperature, and material properties.
It is worth noting that, despite of many intense efforts, no convincing experimental measurements of CL friction forces have been carried out to date. This is due not only to the small magnitude of these forces relative to the “ordinary” Casimir forces (forces of attraction) but also to the imperfections of the measurement layout. In particular, the effective interaction area and relative velocity are significantly limited in the “pendulum” measurement scheme used in Ref. [
12]. Other experimental scenarios [
9,
13,
14,
15,
16,
17] seem to be more exotic. Recently, in Refs. [
15,
16,
17], to measure traces of quantum friction, the authors suggested a scenario in which the nitrogen vacancy center in diamond acquires the geometric phase during rotation at a frequency of 10
3–10
4 rps (revolutions per second) near the Si- or Au-coated surface. Nitrogen vacancy centers have been proposed for use as the main components of quantum computer processors [
18].
Nearly all experiments to measure Casimir–Lifshitz forces (both conservative and dissipative) have been performed with well-conducting materials (metals like gold) under near-normal temperature conditions. Regarding Casimir–Lifshitz friction forces, it has usually been assumed that they decrease with decreasing temperature as the resistivity of metals and ohmic losses decrease. Therefore, at first glance, the friction force also does. The conclusion that for metals, the temperature behavior of CL friction is not that simple was first made in Ref. [
19] and later discussed in [
20,
21]. It has been shown that at temperatures
T <<
θD (
θD is the Debye temperature), the force of friction can increase by several orders of magnitude compared with normal conditions. However, several issues have not been elucidated, in particular, the relation between friction and heating effects at thermal nonequilibrium, the relation between quantum friction and friction at close to zero temperature, and the kinetics of radiation heating, etc. In particular, the interplay between nonequilibrium dynamics, the quantum and thermal properties of the radiation, and the confinement of light at the vacuum-surface interface may lead to several intriguing features caused by nonequilibrium thermodynamics of quantum friction [
22,
23,
24,
25]. Some other effects were considered in Refs. [
26,
27,
28,
29,
30].
The main objective of this paper, in addition to studying CL friction and heating in a system of parallel metallic plates of nonmagnetic metals like gold, is to substantiate the possibility of determining the friction force from thermal measurements. In the calculations, the general results of fluctuation electrodynamics [
31,
32] are used, without a linear expansion in velocity in the basic expressions. It is shown that identical metal plates with different initial temperatures, moving with a constant nonrelativistic velocity,
, relative to each other, rapidly reach a state of quasithermal equilibrium and continue to heat up further. The heating rate is then determined by the power of the friction force.
The outline of this paper is as follows. In
Section 2, the general relations between radiative heating and friction force for parallel plates in relative nonrelativistic motion are given. In
Section 2.2,
Section 2.3,
Section 2.4,
Section 2.5 I consider the simplest case of identical plates of Drude metals having the same material parameters and temperature,
. Analytical expressions are obtained for the friction force of metal plates in the limiting cases of low (zero) temperature and low and high speeds, as well as general expressions describing the kinetics of heating. In
Section 3, the results of the numerical calculations (heating rates of plate 1 and friction parameter,
with
Fx the
x-component of the friction force) are given for different thermal configurations and velocities. The analytical results of
Section 2 are compared with the results of numerical integration according to the general formulas.
Section 4 is devoted to a brief discussion of a possible layout of an experiment for determining the CL friction force by measuring the heating rates of gold plates. Concluding remarks are given in
Section 5.
Appendix A,
Appendix B and
Appendix C contain the details of the analytical calculations. All formulas are written in the Gaussian units,
are the reduced Planck constant and the speed of light in vacuum, respectively,
denotes the absolute temperature and is given in units of energy.
2. General Results
2.1. Radiative Heating and Friction Force for Parallel Plates in Relative Motion
Here, we use the standard formulation of the problem, in which the plates are assumed to be made of homogeneous and isotropic materials with permittivities,
,
, and permeabilities,
, depending on the frequency,
and local temperatures,
and
(
Figure 1).
In line with Refs. [
31,
32], the power,
, of the friction force
x-component,
per unit surface area applied to plate 2 in the laboratory coordinate system associated with plate 1 is calculated using
Here,
and
are the heat fluxes of the plates from a unit surface area per unit time, and
. For all quantities, indices 1 and 2 here and in what follows correspond to the numbering in
Figure 1. Moreover,
and
are calculated in the rest frames of the plates. General relativistic expressions for
and
were obtained in Ref. [
31]. In the nonrelativistic case,
, but taking retardation into account, a more compact form of
and
reads [
32]:
Here,
,
, and
is the gap width in
Figure 1. Variables with a tilde, such as
should be used replacing
. The terms
are defined by the same expressions with appropriate replacements. In the general case, the expressions depending on
and
correspond to the contributions of electromagnetic modes with P (transverse magnetic, TM) and S (transverse electric, TE) polarizations. The quantities
and
are directly related to the heating (cooling) rates of the plates:
and
, where
t denotes the time.
Using Equations (1)–(4), the power of the friction force
takes the form:
Formula (5) can be also recast into a more familiar form in terms of the Fresnel reflection coefficients [
9,
32].
At , due to the symmetry of the system, the heating rates of identical plates are equal. One then has , and the friction force can be determined using the heating rate of any plate. For , it follows that , but and, correspondingly, . This means that when measuring the CL friction force, it is sufficient to control the temperature of only one plate.
2.2. Metal Plates in the Drude Model
In order to treat the problem of temperature-dependent CL friction force between ordinary metals, they are described using the Drude model in terms of plasma frequency,
and damping parameter,
, with
being the resistivity:
Figure 2 plots the dependences
corresponding to the Bloch–Grüneisen (BG) model [
33] and the modified Bloch–Grüneisen (MBG) model (BG scaled in
Figure 2 to the data from Ref. [
34]). In the former case, the residual resistance is zero or can be specified by indicating the effective temperature, below which it is constant. In the MBG model, the residual resistivity is
(see
Figure 2).
Hereinafter, for simplicity, the plates are assumed been made of similar nonmagnetic metal () with the same plasma frequency , but different dependence.
Since for good conductors and the inequality becomes stronger as the terms with in Equations (2), (3) and (5), corresponding to modes with P-polarization, are negligible compared to the terms with , corresponding to modes with S-polarization. Therefore, in what follows, the contributions of P modes are omitted.
When calculating the integrals in (2), (3), and (5), it is convenient to introduce a new frequency variable
, with
and
being the damping parameters of plates 1 and 2 depending on their temperatures
and
, respectively. The absolute value
of the two-dimensional wave vector (using the polar coordinates
in the plane
is expressed as
in the evanescent sector
(
and
in the radiation sector
(
. Here, the parameters
,
,
,
,
, and
are introduced. With these definitions, for
, Equations (2), (3) and (5) take the form:
In the sector
, Formulas (14) and (15) should be modified by replacing
and substituting
for
in Equations (7)–(9) in the integrals over
. The expressions for
can be additionally simplified. For example, it follows that
The
is defined by the same expression (16), substituting
for
and
for
. For two identical plates at quasithermal equilibrium, it follows that
and a simpler useful expression is obtained by expanding the square root in Equation (16) and leaving the expansion terms up to the second order:
In this case, an approximate analytical calculations can be done.
2.3. Quantum Friction
In the case
, corresponding to the conditions of quantum friction, the main role is played by the evanescent modes
. At finite temperatures, the evanescent modes make the dominant contribution at
µm. This range of distances is highly promising experimentally. For this reason, hereinafter, let us consider only evanescent modes, omitting the small term
in Equations (9), (14) and other formulas. Therefore, at zero temperature, substituting the identity
into Equation
yields:
The most straightforward asymptotics of Equation (18) can be worked out for two identical plates in the limit of low velocities,
. Using Equations (15) and (17), one obtains:
Inserting Equation (19) into Equation (18) yields:
where
is the residual resistivity corresponding to the zero-temperature damping factor
The limit of high velocities,
is more laborious. A reasonable representation of the double integral in Equation (9) can be worked out using an approximate expression for
, based on Equation (17):
where
. The product
, as a function of
in the range
reaches its maximum close to the point
, with zeroing at the end points
of the integration domain of the inner integral in Equation (9). At the same time, the dependence on
in
is much weaker. By virtue of this,
was inserted into the denominator of Equation (21) and into
(in the latter case,
. Expression (21) then takes the form
With these transformations, it follows that (see
Appendix A):
where
and
are calculated using Equations (A2) and (A3). The integrals over
and
are calculated explicitly, and finally, we obtain (see Equations (A4) and (A6)):
As follows from Equation (24), in this approximation, the power of the quantum friction force does not depend on the velocity. However, it is worth noting that the condition
implies
, and along with
, it can only be satisfied if
is more than three orders of magnitude smaller than the MBG value shown in
Figure 2. Interestingly, for
, Equation (24) also agrees quite well with the numerical calculations and approximation (20) (see
Section 3.2 below).
2.4. Low Temperatures, Linear in Velocity Approximation
In the quasiequilibrium thermal regime,
, for two identical metal plates in the linear in velocity approximation, Equations (5) and (9) can be recast into the form of [
19,
20]:
In this limit, the friction parameter,
does not depend on
. It is the dependence
in Equation (25) that leads to a large enhancement of friction at low temperatures, when
because the function
weakly depends on
The main contribution to
in this case makes the values
,
, and one can again use Equation (17) for
. Meantime,
(this is a suitable approximation at
Making use of these simplifications in Equation (26), one arrives at (see
Appendix B)
where
and
are the Struve and Neumann functions [
35], respectively. Using the series representations of these functions yields:
A more straightforward and physically transparent low-temperature representation of Equation (25) is obtained by using the relation
between the damping factor and resistivity, yielding
Combining the relation
, which implies
, and
, which implies the limit of low velocities
one concludes that the Formula (29) holds at
As a result, the conditions of a low-temperature increase in friction and the applicability of the low-speed approximation are met at . For gold, at , this implies .
According to Refs. [
19,
20], the dependence (29) is associated with a growing penetration depth of S-polarized electromagnetic modes and an increase in their density at low temperatures. A significant low-temperature increase in the friction parameter was also noted in the case of the movement of a metal particle above the metal surface [
21].
2.5. Low Temperatures, High-Velocity Limit
The limit
at finite but low temperatures (
can be analyzed similarly to the case of zero temperatures using the properties of the function (21). When substituting Equation (21) into Equation (9) with allowance for Equation (13), the first exponential term in Equation (13) makes the dominant contribution at
. Due to this, let us take an advantage of the substitution
in the denominator of Equation (21). For the second term in Equation (13), a new variable
is introduced and the substitution
is made in the denominator of Equation (21), while the integral (9) is then determined using the large exponential factor
at
. Then, taking into account these transformations in Equation (9), and summing both contributions, the double integral in Equation (9) finally takes the form (see
Appendix C)
where
and
To proceed further, we replace the function
with
in the inner integral (31), which is again a good approximation for
. The remaining integral yields:
where
is a complete elliptic integral [
35]. Taking this into account, the
-integral in Equation (31) can be evaluated as the arithmetic mean between the integrals calculated with the limit functions on the left and right sides of the inequality (see
Appendix C):
Substituting Equation (A20) into Equations (31) and (9) finally yields:
where
is calculated using:
Similar to the case of quantum friction (24), the power of the friction force (37) does not depend on velocity.
To date, there are no other relevant calculations for the friction forces between metal plates, corresponding to low-temperature conditions. However, it is interesting to compare the results obtained here with those in the case of an atom moving above the metal surface [
22,
23]. Let us compare the dependences on the velocity and resistivity of metal for quantum friction (
Equation (20) has the same low-speed dependence,
, but the opposite dependence on resistivity (Equation (16) in Ref. [
22]):
(yet the additional condition
should be met). At high but nonrelativistic velocities (
, Equation (24) yields
in contrast to
(Equation (19) in Ref. [
22]). The dependence on velocity in Equation (24) is more moderate,
, which qualitatively agrees with that in Ref. [
22]. The case
when the friction force is linear in velocity, is less informative, because the results of Refs. [
22,
23] correspond to room conditions. Yet, Formula (29) yields
(assuming the condition (30)), which is different from that in Ref. [
23]:
.
In general, one should not expect close qualitative similarity between the plate–plate and atom–plate configurations because in the latter case, the radiative energy exchange processes, according to Refs. [
22,
23], are determined by the specific thermal nonequilibrium in the system. In the case of macroscopic bodies, such as two plates, the system must reach a state of thermal quasiequilibrium; see
Section 2.6 just below.
2.6. Kinetics of Heating of Plates
The heat transfer of plates is described using the equations
With
being the specific heat capacities,
and
are the thicknesses and densities of materials,
and
are defined using Equations (2) and (3), and the temperature gains
correspond to the interval of time
. The dependences
and
can be determined using the equation
For identical plates, in Equation (40), one can use the replacements , . In what follows, only this case is considered.
When writing Equations (39) and (40), it is also assumed that the heat exchange due to radiative heat transfer occurs much slower than under thermal diffusion, and the plates acquire equal temperature at all points because of high thermal conductivity. Using the thermal diffusion equation along the normal to the plates,
, the characteristic time of the heat diffusion necessary to reach thermal quasiequilibrium, is
(where
, and
is the thermal conductivity). Then it follows that
and in the case of gold at
and
J/(kg·K),
W/(m·K),
kg/m
3 [
36]) one obtains
. In turn, the kinetics of heating induced by friction takes dozens of seconds or minutes (see
Section 3.3 below), depending on the velocity and other parameters. Assuming that
, from Equation (39) one obtains:
the heating time from the initial temperature
to the final temperature
. In the simplest case of
and
(this is a typical low-temperature dependence for metals), it follows from Equation (41) that
where
. At
and relatively low velocities of plate 2, as follows from numerical calculations (see
Section 3.1 below), the heating/cooling rates of metal plates differ only in sign, i.e.,
. This is the normal mode of heat transfer, when a hotter body cools down, and a colder one heats up. Then the left sides of Equation (39) can be equated, and the corresponding quasistationary temperature of the plates is calculated using:
where
and
are their initial temperatures. After establishing quasithermal equilibrium, the temperature of the plates will increase according to Equations (41) and (42).
4. Experimental Proposal
Initiated by the advantage of the experimental design [
15,
16,
17] to measure the quantum friction force, I suggested [
42] using another experimental layout, as shown schematically in
Figure 9. Unlike in Ref. [
17], where the setup includes a disk 10 cm in diameter rotating with an angular frequency of up to
rps, it is proposed to use two identical disks placed in one thermostat, one of which rotates at a controlled speed. In the peripheral region, the disks have an annular metal coating with an effective area
. The non-inertiality of the reference system of disk 2 does not appear in this case because the rotation frequency is small compared to the characteristic frequencies of the fluctuation electromagnetic field. Accordingly, the original expressions (2) and (3), for heating rates remain valid.
A possible measurement scenario in this case is the quasiequilibrium thermal mode, in which the temperatures of plates increase from the initial temperature
at the same rate. It should be noted that the experimental design must take into account possible limitations on angular velocity imposed by the tensile strength of the material used. Assuming that the main body of the plate is made of gold, a quite moderate assessment of the linear velocity of the far-distant annular parts of the plate yields
, with
N/m
2 and
kg/m
3 being the tensile strength and density of gold [
35]. By plugging the numerical numbers into the above condition, one obtains
m/s, or
rps, which seems to be a well acceptable value.