Next Article in Journal
Long-Term Optical Monitoring of Broad-Line AGNs (LoTerm AGN): Case Study of NGC 3516
Next Article in Special Issue
The Normal Casimir Force for Lateral Moving Planes with Isotropic Conductivities
Previous Article in Journal
A Thermodynamic Comparison of Nanotip and Nanoblade Geometries for Ultrafast Laser Field Emission via the Finite Element Method
Previous Article in Special Issue
Vacuum Interaction of Topological Strings at Short Distances
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Casimir–Lifshitz Frictional Heating in a System of Parallel Metallic Plates

Institute of Informatics, Electronics and Robotics, Kabardino-Balkarian State University, Chernysheskogo 174, Nalchik 360004, Russia
Physics 2024, 6(1), 13-30; https://doi.org/10.3390/physics6010002
Submission received: 2 September 2023 / Revised: 21 November 2023 / Accepted: 28 November 2023 / Published: 27 December 2023
(This article belongs to the Special Issue 75 Years of the Casimir Effect: Advances and Prospects)

Abstract

:
The Casimir–Lifshitz force of friction between neutral bodies in relative motion, along with the drag effect, causes their heating. This paper considers this frictional heating in a system of two metal plates within the framework of fluctuation electromagnetic theory. Analytical expressions for the friction force in the limiting cases of low (zero) temperature and low and high speeds, as well as general expressions describing the kinetics of heating, have been obtained. It is shown that the combination of low temperatures (T < 10 K) and velocities of 10–103 m/s provides the most favorable conditions when measuring the Casimir–Lifshitz friction force from heat measurements. In particular, the friction force of two coaxial disks of gold 10 cm in diameter and 500 nm in thickness, one of which rotates at a frequency of 10–103 rps (revolutions per second), can be measured using the heating effect of 1–2 K in less than 1 min. A possible experimental layout is discussed.

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 103–104 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, V , 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, T . 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, η = F x / V , 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,   a n d   c are the reduced Planck constant and the speed of light in vacuum, respectively, T 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, ε 1 , ε 2 , and permeabilities, μ 1 ,   μ 2 , depending on the frequency, ω , and local temperatures, T 1 and T 2 (Figure 1).
In line with Refs. [31,32], the power, F x V , of the friction force x-component, F x , per unit surface area applied to plate 2 in the laboratory coordinate system associated with plate 1 is calculated using
F x V = P 1 + P 2 / γ .
Here, P 1 and P 2 are the heat fluxes of the plates from a unit surface area per unit time, and γ = 1 V 2 / c 2 1 / 2 . For all quantities, indices 1 and 2 here and in what follows correspond to the numbering in Figure 1. Moreover, P 1 and P 2 are calculated in the rest frames of the plates. General relativistic expressions for P 1 and P 2 were obtained in Ref. [31]. In the nonrelativistic case, V / c 1 , but taking retardation into account, a more compact form of P 1 and P 2 reads [32]:
P 1 = 4 π 3 0 d ω ω d 2 k q 2 Q μ 2 I m q 1 μ 1 I m q ~ 2 μ 2 c o t h ω 2 T 1 c o t h ω 2 T 2 + μ 1,2 ε 1,2 ,
P 2 = 4 π 3 0 d ω d 2 k ω q 2 Q μ 2 I m q 1 μ 1 I m q ~ 2 μ 2 c o t h ω 2 T 1 c o t h ω 2 T 2 + μ 1,2 ε 1,2 ,
Q μ = q + q 1 / μ 1 ) ( q + q ~ 2 / μ ~ 2 ) e x p ( q a q q 1 / μ 1 ) ( q q ~ 2 / μ ~ 2 ) ( q a ,   Q μ Q ε .
Here, ω = ω k x V , q = k 2 ω 2 / c 2 , q 1,2 = k 2 ε 1,2 μ 1,2 ω 2 / c 2 , and a is the gap width in Figure 1. Variables with a tilde, such as q ~ 2 , should be used replacing ω ω . The terms μ 1,2 ε 1,2 are defined by the same expressions with appropriate replacements. In the general case, the expressions depending on ε 1,2 and μ 1,2 correspond to the contributions of electromagnetic modes with P (transverse magnetic, TM) and S (transverse electric, TE) polarizations. The quantities P 1 and P 2 are directly related to the heating (cooling) rates of the plates: d Q 1 / d t = P 1 and d Q 2 / d t = P 2 , where t denotes the time.
Using Equations (1)–(4), the power of the friction force F x V = P 1 + P 2 takes the form:
F x V = 4 π 3 0 d ω d 2 k k x V q 2 Q μ 2 I m q 1 μ 1 I m q ~ 2 μ 2 c o t h ω 2 T 1 c o t h ω 2 T 2 + μ 1,2 ε 1,2 .
Formula (5) can be also recast into a more familiar form in terms of the Fresnel reflection coefficients [9,32].
At T 1 = T 2 = T , due to the symmetry of the system, the heating rates of identical plates are equal. One then has F x V = 2 P 1,2 , and the friction force can be determined using the heating rate of any plate. For T 1 T 2 , it follows that P 1 P 2 , but P 1 T 1 , T 2 = P 2 T 2 , T 1 and, correspondingly, P 1 T 1 , T 2 + P 2 T 1 , T 2 = P 1 T 1 , T 2 + P 1 T 2 , T 1 = P 2 T 1 , T 2 + P 2 T 2 , T 1 . 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, ω p , and damping parameter, ν T = ω p 2 ρ ( T ) / 4 π , with ρ ( T ) being the resistivity:
ε ω = 1 ω p 2 ω ω + i · ν T .
Figure 2 plots the dependences ρ ( T ) 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 ρ 0 = 2.3 × 10 10 Ω · m (see Figure 2).
Hereinafter, for simplicity, the plates are assumed been made of similar nonmagnetic metal ( μ 1 = μ 2 = 1 ) with the same plasma frequency ω p , but different ν T dependence.
Since ε ω 1 for good conductors and the inequality becomes stronger as T 0 , the terms with ε 1,2 in Equations (2), (3) and (5), corresponding to modes with P-polarization, are negligible compared to the terms with μ 1,2 , 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 ω = ν m T 1 , T 2 t , with ν m T 1 , T 2 = m a x ( ν 1 T 1 , ν 2 T 2 ) and ν i T i ,   i = 1 ,   2 , being the damping parameters of plates 1 and 2 depending on their temperatures T 1 and T 2 , respectively. The absolute value k of the two-dimensional wave vector (using the polar coordinates k , ϕ in the plane ( k x , k y ) ) is expressed as k = ( ω p / c ) y 2 + β m 2 t 2 in the evanescent sector k > ω / c ( 0 y < ) and k = ω p / c β m 2 t 2 y 2 in the radiation sector k < ω / c ( 0 y β m t ) . Here, the parameters β m = ν m / ω p , α i = ν i / T i , γ i = ν i / ν m , λ = ω p a / c , ζ = ( V / c ) β m 1 , and K = h ν m 2 ω p / c 2 2 π 3 are introduced. With these definitions, for k > ω / c , Equations (2), (3) and (5) take the form:
P 1 = K 0 d t 0 d y y 3 f 1 t , y ,
P 2 = K 0 d t 0 d y y 3 f 2 t , y ,
F x V = K 0 d t 0 d y y 3 y 2 + β m 2 t 2 f 3 t , y ,
f 1 t , y = t 0 π d ϕ   I m w 1 I m w 2 D 2 Z t , y , ϕ ,
f 2 t , y = 0 π d ϕ t I m w 1 I m w 2 D 2 Z t , y , ϕ ,
f 3 t , y = ζ 0 π d ϕ c o s ϕ I m w 1 I m w 2 D 2 Z t , y , ϕ ,
Z t , y , ϕ = c o t h α 1 t 2 c o t h α 2 t 2 ,
w 1 = y 2 + t t + i · γ 1   ,   w 2 = y 2 + t t + i · γ 2   ,   t = t ζ c o s ϕ y 2 + β m 2 t 2 ,
D = y + w 1 y + w 2 exp λ y y w 1 y w 2 e x p ( λ y ) .
In the sector k < ω / c , Formulas (14) and (15) should be modified by replacing y i y and substituting β m t for in Equations (7)–(9) in the integrals over y . The expressions for I m w 1,2 can be additionally simplified. For example, it follows that
I m w 1 = γ 1 2 y 2 + 1 + y 2 t 2 2 + γ 1 2 t 2 γ 1 2 y 2 1 + y 2 t 2 1 / 2 s g n ( t ) 2 γ 1 2 + t 2 .
The I m w 2 is defined by the same expression (16), substituting γ 2 for γ 1 and t for t . For two identical plates at quasithermal equilibrium, it follows that γ 1 = γ 2 = 1 , 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:
I m w 1 t · s g n ( t ) 2 1 + t 2 y 2 + 1 + y 2 t 2 .
In this case, an approximate analytical calculations can be done.

2.3. Quantum Friction

In the case T 1 = T 2 = 0 , corresponding to the conditions of quantum friction, the main role is played by the evanescent modes k > ω / c . At finite temperatures, the evanescent modes make the dominant contribution at a < 1 µm. This range of distances is highly promising experimentally. For this reason, hereinafter, let us consider only evanescent modes, omitting the small term β m 2 t 2 in Equations (9), (14) and other formulas. Therefore, at zero temperature, substituting the identity Z t , y , ϕ = s g n t s g n t ζ y c o s ϕ into Equation 9 yields:
F x V = 2 K ζ 0 d y y 4 0 π 2 d ϕ c o s ϕ 0 ζ y c o s ϕ d t I m w 1 I m w 2 D 2 .
The most straightforward asymptotics of Equation (18) can be worked out for two identical plates in the limit of low velocities, ζ 1 . Using Equations (15) and (17), one obtains:
I m w 1 I m w 2 t ( ζ y c o s ϕ t ) 4 y 2 ,   D 2 1 16 y 4 exp 2 λ y .
Inserting Equation (19) into Equation (18) yields:
F x V = π K 2 12 ζ 4 λ 2 = ω p 2 2 13 π 2 ω p c 2 V c 2 V a ν 0 2 = 1 2 9 ρ 0 2 a 2 V c 4 ,
where ρ 0 is the residual resistivity corresponding to the zero-temperature damping factor ν 0 = ν 0 . The limit of high velocities, ζ 1 is more laborious. A reasonable representation of the double integral in Equation (9) can be worked out using an approximate expression for I m w 1 I m w 2 , based on Equation (17):
I m w 1 I m w 2 = t b t s g n t s g n ( b t ) 4 1 + t 2 1 + ( t b ) 2 y 2 + t 2 ( 1 + y 2 ) y 2 + ( t b ) 2 ( 1 + y 2 ) 1 / 2 .
where b = ζ y c o s ϕ . The product I m w 1 I m w 2 , as a function of t in the range 0 t b , reaches its maximum close to the point t = b / 2 , with zeroing at the end points t = 0 ,   t = b of the integration domain of the inner integral in Equation (9). At the same time, the dependence on t in D 2 is much weaker. By virtue of this, t = b / 2 was inserted into the denominator of Equation (21) and into D 2 (in the latter case, c o s ϕ 1   i s a l s o u s e d ) . Expression (21) then takes the form
I m w 1 I m w 2 t b t 4 1 + b 2 / 4 y 2 + b 2 ( 1 + y 2 / 4 ) .
With these transformations, it follows that (see Appendix A):
F x V = K ζ 2 0 d y y 4 D 2 0 π / 2 d ϕ c o s ϕ ψ y , ϕ 0 ζ y c o s ϕ d t t ζ y c o s ϕ t ,
where D 2 and ψ y , ϕ are calculated using Equations (A2) and (A3). The integrals over t and ϕ are calculated explicitly, and finally, we obtain (see Equations (A4) and (A6)):
F x V = ν 0 2 3 π 2 ω p c 2 0 d y y 5 exp 2 λ y y + 1 + y 2 4 1 y 2 1 + y 2 2 y 2 y 2 ζ 2 + 4 + 2 ( 1 + y 2 ) ( 1 + y 2 ) ζ 2 + 4 .
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 ζ 1 implies ρ 0 ω p 4 π V / c , and along with V / c 1 , it can only be satisfied if ρ 0 is more than three orders of magnitude smaller than the MBG value shown in Figure 2. Interestingly, for ζ 1 , 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, T 1 = T 2 = T , for two identical metal plates in the linear in velocity approximation, Equations (5) and (9) can be recast into the form of [19,20]:
F x V = V 2 8 π 2 ω p c 4 α 1 Y 1 λ , α ,
Y 1 ( λ , α ) = α 2 0 d t s i n h 2 ( α t / 2 ) 0 d y y 5 I m w 1 2 D 2 .
In this limit, the friction parameter, η = F x / V , does not depend on V . It is the dependence F x α 1 in Equation (25) that leads to a large enhancement of friction at low temperatures, when α = ν T / T 0 , because the function Y 1 λ , α weakly depends on α . The main contribution to Y 1 λ , α in this case makes the values t < 1 , y ~ 1 / 2 λ ~ 1 , and one can again use Equation (17) for I m w 1 . Meantime, α 2 s i n h 2 ( α t / 2 ) 4 / t 2 (this is a suitable approximation at α < 0.3 )   a n d   D 2 16 y 4 exp 2 λ y . Making use of these simplifications in Equation (26), one arrives at (see Appendix B)
Y 1 λ , α χ λ = π 32 π 4 λ H 1 2 λ N 1 ( 2 λ ) 1 4 λ 2 ,
where H 1 ( x ) and N 1 x are the Struve and Neumann functions [35], respectively. Using the series representations of these functions yields:
χ λ = π 64 1 λ 1 2 λ 2 + 1 4 λ 3 + , λ 1 , π 32 1 4 + 2 λ 3 2 λ l n λ + 0.577 , λ 1 .
A more straightforward and physically transparent low-temperature representation of Equation (25) is obtained by using the relation ν T = ω p 2 ρ ( T ) / 4 π between the damping factor and resistivity, yielding
F x V = 1 2 π V c 2 ω p c 2 T ρ T χ λ .
Combining the relation α 1 , which implies ν T T , and ζ 1 , which implies the limit of low velocities V / c ν ( T ) / ω p , one concludes that the Formula (29) holds at
ω p V / c ν T T .
As a result, the conditions of a low-temperature increase in friction and the applicability of the low-speed approximation are met at V / c T / ω p . For gold, at T = 1   K , this implies V / c 1.5 × 10 5 .
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 ζ 1 at finite but low temperatures ( α 1 ) 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 t ~ 1 b = ζ y c o s ϕ . Due to this, let us take an advantage of the substitution t b b in the denominator of Equation (21). For the second term in Equation (13), a new variable t = t b is introduced and the substitution t + b b is made in the denominator of Equation (21), while the integral (9) is then determined using the large exponential factor exp α t 1 1 at t ~ 1 b . 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)
I ζ 2 2 0 d y y 4 D 2 0 π 2 d ϕ ψ 1 ( y , ϕ ) 0 d t t e α t 1 ψ 2 ( y , ϕ ) ,
where
D 2 y + 1 + y 2 4 exp 2 λ y ,
ψ 1 ( y , ϕ ) = c o s 2 ϕ 1 + ζ 2 y 2 c o s 2 ϕ 1 / 2 1 + ζ 2 ( 1 + y 2 ) c o s 2 ϕ 1 / 2 ,
and
ψ 2 y , t = 1 1 + t 2 y 2 + t 2 ( 1 + y 2 1 / 2 .
To proceed further, we replace the function t / e α t 1 with 1 / α in the inner integral (31), which is again a good approximation for α < 0.3 . The remaining integral yields:
0 d t t e α t 1 ψ 2 y , t 1 α K q 1 + y 2 ,   q = 1 + y 2 1 / 2 ,
where K q 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):
c o s 2 ϕ 1 + ζ 2 1 + y 2 c o s 2 ϕ 1 < ψ 1 y , ϕ < c o s 2 ϕ 1 + ζ 2 y 2 c o s 2 ϕ 1 .
Substituting Equation (A20) into Equations (31) and (9) finally yields:
F x V = T 8 π 2 V ω p c 3 Y 2 λ , ζ ,
where Y 2 ( λ , ζ ) is calculated using:
Y 2 λ , ζ = 0 d y y 4 e 2 λ y y + y 2 + 1 4 1 + y 2 1 2 K 1 1 + y 2 1 + ζ 2 y 2 1 y 2 1 + ζ 2 y 2 + 1 + ζ 2 1 + y 2 1 1 + y 2 1 + ζ 2 1 + y 2 .
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 ( T = 0 ) . Equation (20) has the same low-speed dependence, F x ~ V 3 , but the opposite dependence on resistivity (Equation (16) in Ref. [22]): F x ~ ρ 0 2 (yet the additional condition ρ 0 ω p 4 π V / c should be met). At high but nonrelativistic velocities ( ρ 0 ω p 4 π V / c ) , Equation (24) yields F x ~ ρ 0 2 in contrast to F x ~ ρ 0 (Equation (19) in Ref. [22]). The dependence on velocity in Equation (24) is more moderate, F x ~ V 1 , which qualitatively agrees with that in Ref. [22]. The case T > 0 , 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 F x ~ 1 / ρ 0 (assuming the condition (30)), which is different from that in Ref. [23]: F x ~   ρ ( T ) 2 .
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
P 1 T 1 , T 2 t = h 1 ρ 1 c 1 T 1 T 1   a n d   P 2 T 1 , T 2 t = h 2 ρ 2 c 2 T 2 T 2 .
With c i T i being the specific heat capacities, h i and ρ i are the thicknesses and densities of materials, P 1 T 1 , T 2 and P 2 T 1 , T 2 are defined using Equations (2) and (3), and the temperature gains T i correspond to the interval of time t . The dependences T 2 ( T 1 ) and T 1 T 2 can be determined using the equation
d T 2 d T 1 = P 1 ( T 1 , T 2 ) c 1 ( T 1 ) h 1 ρ 1 P 2 ( T 1 , T 2 ) c 2 ( T 2 ) h 2 ρ 2 .
For identical plates, in Equation (40), one can use the replacements P 2 ( T 1 , T 2 ) P 1 ( T 2 , T 1 ) , P 1 ( T 1 , T 2 ) P 2 ( T 2 , T 1 ) . 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, T / t = a 2 2 T / z 2 , the characteristic time of the heat diffusion necessary to reach thermal quasiequilibrium, is τ = h 2 / a 2 (where a 2 = κ / c ρ , and κ is the thermal conductivity). Then it follows that τ = h 2 c ρ / κ and in the case of gold at T = 10   K and h 1,2 = h = 500   µ m   c = 2.2 J/(kg·K), κ = 3200 W/(m·K), ρ = 19.8 · 10 3 kg/m3 [36]) one obtains τ 3   μ s . 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 P i T , T = 0.5 η ( T , V ) V 2 , from Equation (39) one obtains:
t = 2 h ρ V 2 T 0 T c ( T ) η ( T , V ) d T ,
the heating time from the initial temperature T 0 to the final temperature T . In the simplest case of η = c o n s t and c T = a 1 T + a 2 T 3 (this is a typical low-temperature dependence for metals), it follows from Equation (41) that
T t = β + β 2 + T 0 4 + 2 β T 0 2 + 2 η V 2 t / h ρ a 2 1 / 2 ,
where β = a 1 / a 2 . At T 1 T 2 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., P 1 T 1 , T 2 = P 2 T 1 , T 2 . 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:
T = β + β 2 + β T 1 2 + T 2 2 + T 1 4 + T 2 4 1 / 2 ,
where T 1 and T 2 are their initial temperatures. After establishing quasithermal equilibrium, the temperature of the plates will increase according to Equations (41) and (42).

3. Numerical Results

For an ideal metal without impurities and defects, within the BG model, the damping frequency ν T in Equation (6) is defined using the formula [33]:
ν T = 0.0212 Θ / T 5 0 Θ / T d x x 5 s i n h 2 ( x / 2 )   e V .
Numerical calculations were performed using Equation (44) and the MBG approximation shown in Figure 2 (BG scaled). The used plasma frequency of gold is ω p = 9.03 eV. All calculations were performed with a gap width of a = 10 nm (Figure 1) unless another value is indicated. It should be noted that at distances of a > 10 nm, the processes of radiative heat transfer and friction due to tunneling of electrons and phonons [37,38] seem do not occur or become insignificant [39,40,41].

3.1. Quantum Friction

Figure 3 shows the velocity-dependent quantum friction force between the plates of gold, calculated using Formulas (20) (green line), (9) (red line) and (24) (blue line).
The curves in Figure 3a,b were calculated at residual resistances of 2.13 × 10 13 Ω·m and 2.3 × 10 10 Ω·m, which correspond to the BG model (44) at T = 5 K and the MBG model at T = 0 K. Note that in the latter case, the residual resistance coincides with that defined using Formula (44) at T = 20.9 K.

3.2. Temperature-Dependent Friction at Thermal Quasiequlibrium

Figure 4 shows the plots of the friction parameter, η = F x / V , depending on the temperature T of the gold plates, corresponding to the BG and MBG models. The curves with symbols were calculated using Equation (9) for V = 1 m/s. Solid curves were plotted using approximation (25) along with Equation (26) (green lines) or Equation (28) at λ 1 (blue lines). In Figure 4a, both solid lines merge. The presence of maxima and their positions on the curves agree with Equations (29) and (30), respectively. These results show that the linear-in-velocity approximation is valid only to the right of the maxima of the η ( T ) dependences.
Figure 5a,b demonstrates the velocity dependences of η in the BG and MBG models. The red, blue, and green lines correspond to quasiequilibrium temperatures of 5 K, 10 K, and 77 K, respectively. The different (temperature) order of lines in Figure 5a compared to Figure 5b is explained by the high residual resistance of gold in the MBG model: the condition ν T < T , which is necessary for the law-temperature increase in friction, is violated at T = 5 and 10 K.
Table 1 shows the calculated values of the friction parameter η of the gold plates at V = 1 m/s, depending on the temperature, T , and separation distance, a. Similar to that in Figure 4, one can note the effect of increasing friction (up to a maximum) with decreasing temperature at T < θ D , which is more better expressed in the BG model. The height of this maximum depends on the velocity-to-resistivity ratio. When the temperature becomes sufficiently low, condition (30) is violated, and the coefficient of friction decreases.
The dependence of η on the separation distance a in all the cases is close to inverse proportionality ( η a 1 ). This is clearly seen from the data in Table 1 and agrees with our previous results [19,20,32].

3.3. Friction and Heating under Different Conditions

Figure 6 and Figure 7 show the calculated heating rates of plate 1 (Figure 6a and Figure 7a) and friction parameters (Figure 6b and Figure 7b), depending on the velocity V of plate 2 for various thermal configurations.
One can see that at V < 10 m/s (Figure 6a) and V < 10 2 m/s (Figure 7a), the heating rates of plates 1 and 2 are equal in absolute value, differing in sign. According to their temperatures, T 1 = 4 K and T 2 = 6 K, plate 1 heats up and plate 2 cools down, realizing the “normal” heat exchange regime. At the same time, the friction parameters weakly depend on the temperature (Figure 6b and Figure 7b). When the speed of plate 2 increases, both plates heat up faster. Then, one can see the effect of the “anomalous” heating of plate 2 for some time, when it continues to heat up despite the higher temperature. This is similar to the case of heating a hotter metal particle moving above a cold surface [21]. However, due to different absolute values of the heating rates (cf. the upper and lower lines shown with open squares (□) in Figure 6a and Figure 7a), the temperature of plate 1 “catches up” with the temperature of plate 2, and further on, both plates heat up at the same rate.
The drop in friction parameters for high velocities of plate 2 (Figure 6b and Figure 7b) is explained by the change in sign of the Doppler-shifted frequency ω = ω k x V = ω k V c o s ϕ in Equation (5). This occurs at V > ν T a because the characteristic absorption frequency is ω ~ ν T and the characteristic wave vector is k ~ 1 / a . The positions of the “kinks” on the curves η ( V ) in Figure 6 and Figure 7 correlate with resistivity because ν T ~ ρ ( T ) . Indeed, it follows from Figure 2 that ρ M B G / ρ B G = 10 2 10 3 at T = 4–6 K. At the same time, the ratio η M B G / η B G in this case is inversely proportional to resistivities (see Equation (29) and Table 1).
In general, as follows from the calculations for all considered temperatures and velocities (Figure 5, Figure 6 and Figure 7, Table 1), the maximum friction parameter in the BG and MBG models (at a = 10 nm) is 10 6 10 3 kg/(m2·s).
Figure 8 shows the heating time of the plates versus the velocity of plate 2, calculated using numerical integration of Equation (40) from 4 K to 5 K and from 4 K to 8 K. In these calculations, the fitting parameters a 1 = 0.0035   J / ( k g · K 2 ) and a 2 = 0.0023   J / ( k g · K 4 ) of the dependence c T = a 1 T + a 2 T 3 were determined using the data [36] for gold at T < 20   K .
As follows from Figure 8, quite comfortable (from the experimental point of view) values of the plate heating times (1–100 s) can be obtained in the velocity range 1 10 3 m/s. On the contrary, heating by 1 K at T 0 = 300 K, a = 10 nm, and V = 10 3 m/s will take about 2 h. Thus, low-temperature thermal measurements have exceptional advantages over measurements under normal conditions due to a significant reduction in measurement time and the elimination of noise and other undesirable effects.

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 7 × 10 3 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 π D w . 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 T 0 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 V < σ 0 / ρ 0.5 , with σ 0 = 1.5 × 10 8 N/m2 and ρ = 1.93 × 10 4 kg/m3 being the tensile strength and density of gold [35]. By plugging the numerical numbers into the above condition, one obtains V ~ 10 2 m/s, or n   ~   300 rps, which seems to be a well acceptable value.

5. Concluding Remarks

The Casimir–Lifshitz friction force mediated by the fluctuating electromagnetic field between metal plates moving with constant velocity relative to each other causes their heating. In a state out of thermal equilibrium, “anomalous” heating of the moving plate can be observed when it is heated for some time despite the higher temperature. However, the system rapidly reaches a state of thermal quasiequilibrium. At low temperatures T   θ D , the Casimir–Lifshitz friction and heating of metal plates increase significantly (see Equations (29) and (30)) while the heat capacity decreases. In combination with a fairly high speed of movement, this provides a fairly short heating time, which is convenient for experiments (see Equation (41)).

Funding

This research received no external funding.

Data Availability Statement

Data can be obtained upon reasonable request.

Acknowledgments

I am grateful to Carsten Henkel for fruitful remarks at the preliminary stage of the work and for providing data on the low-temperature resistivity of gold.

Conflicts of Interest

The author declares no conflicts of interest.

Appendix A. Evaluation of the Integral (24)

Substituting t = b / 2 = ζ y c o s ϕ / 2 into Equation (22) takes into account that typically, ζ y c o s ϕ / 2 1 (since ζ 1 , y ~ 1 / 2 λ ~ 1 ) and w 1,2 y 2 + 1 . Then D 2 takes the form
D 2 y + y 2 + 1 4 e x p 2 λ y .
With these simplifications, Formula (18) reads:
F x V = K ζ 2 0 d y y 2 y + y 2 + 1 4 e 2 λ y 0 π / 2 d ϕ c o s ϕ ψ y , ϕ 0 ζ y c o s ϕ d t t ζ y c o s ϕ t ,
ψ y , ϕ = 1 + ζ 2 y 2 c o s 2 ϕ / 4 1 + ζ 2 ( 1 + y 2 ) c o s 2 ϕ / 4 .
The t -integral in Equation (A2) is just ζ 3 y 3 c o s 3 ϕ / 6 , while the integral over ϕ is
I ϕ = ζ 3 y 3 6 0 π / 2 d ϕ c o s 4 ϕ ψ y , ϕ = 8 y 3 3 ζ 0 π / 2 d ϕ c o s 4 ϕ u 2 + y 2 c o s 2 ϕ u 2 + ( 1 + y 2 ) c o s 2 ϕ ,
where u = 2 / ζ . The integral in Equation (A4) is calculated explicitly using the table integral [35]
0 π / 2 d ϕ a 2 + b 2 c o s 2 ϕ = π 2 1 a a 2 + b 2 .
Using Equation (A5) yields:
I ϕ = 4 π y 3 3 ζ 1 y 2 ( 1 + y 2 ) 2 y 2 y 2 ζ 2 + 4 + 2 ( 1 + y 2 ) ( 1 + y 2 ) ζ 2 + 4 .
Substituting Equation (A6) into Equation (A2) yields Equation (24).

Appendix B. Evaluation of the Integral (26)

In the case α 1 , the main contribution to Equation (26) makes the values t < 1 , y ~ 1 / 2 λ ~ 1 . Then, from Equation (14), it follows that w 1,2 = y 2 + t / t + i 1 / 2 y . Using this, one finds:
D 2 16 y 4 e x p 2 λ y .
Meantime, from Equation (21) it follows that
I m w 1,2 2 = t 2 4 1 + t 2 y 2 + t 2 1 + y 2 .
Substituting Equations (A7) and (A8) into Equation (26) yields:
Y λ , α = α 2 64 0 d y y e 2 λ y 0 d t t 2 s i n h ( α t / 2 ) 2 1 ( 1 + t 2 ) y 2 + t 2 ( 1 + y 2 ) .
Using the approximation t 2 s i n h ( α t / 2 ) 2 4 / α 2 and the table integral [35]
0 d x 1 a 2 + x 2 1 b 2 + x 2 = π 2 a b ( a + b ) ,
one obtains:
Y 1 λ , α χ λ = π 32 0 d y e 2 λ y y + 1 + y 2 = π 32 π 4 λ H 1 2 λ N 1 2 λ 1 4 λ 2 ,
where H 1 x and N 1 x are the Struve and Neumann functions [35].

Appendix C. Evaluation of the Integral (31)

Let us rewrite Equation (13) in the form
Z t , y , ϕ = 2 exp α t 1 2 exp α t 1   ,   t = t ζ y c o s ϕ .
The integral in Equation (9) includes two exponential factors, defined by Equation (A12). By changing the order of integration in the first term, one obtains:
I = 2 ζ 0 d y y 4 0 π d ϕ c o s ϕ 0 d t I m w 1 I m w 2 D 2 1 exp α t 1 .
Similar to Appendix B, one can again take advantage of the behavior of the t - integral in Equation (A13) for α 1 , and ζ 1 , substituting 1 exp α t 1 1 / α t and using Equation (A1) for D 2 . For I m w 1 I m w 2 , let us use Equation (21) with the replacement t b ζ y c o s ϕ . Then Equation (21) takes the form
I m w 1 I m w 2 = t · ζ c o s ϕ · s i g n ( t ζ y c o s ϕ ) 4 1 + ζ 2 y 2 c o s 2 ϕ 1 + ζ 2 1 + y 2 c o s 2 ϕ 1 + t 2 y 2 + t 2 1 + y 2 1 / 2 .
Inserting Equation (A14) into Equation (A13) yields:
I = ζ 2 α 0 d y y 4 e 2 λ y y + y 2 + 1 4 1 + y 2 1 / 2 0 π / 2 d ϕ ψ 1 ( y , ϕ ) 0 d t ψ 2 t , y ,
where
ψ 1 ϕ , y = c o s 2 ϕ 1 + ζ 2 y 2 c o s 2 ϕ 1 + ζ 2 1 + y 2 c o s 2 ϕ 1 / 2 ,
ψ 2 t , y = 1 1 + t 2 y 2 1 + y 2 + t 2 1 / 2 .
Substituting Equations (A16) and (A17) into Equation (A15) and taking into account Equations (A10) and (35), the inner integrals are calculated yielding
I ϕ ( y ) = 0 d ϕ c o s 2 ϕ ψ 1 ϕ , y π 2 ζ 2 1 + ζ 2 y 2 1 y 2 1 + ζ 2 y 2 + 1 + ζ 2 1 + y 2 1 1 + y 2 1 + ζ 2 1 + y 2 ,
I t ( y ) = 0 d t ψ 2 t , y = 1 1 + y 2 K 1 1 + y 2 ,
where K x is the elliptic integral. Finally, substituting Equations (A18) and (A19) into Equation (A15) yields:
I = π 4 α 0 d y y 4 e 2 λ y y + y 2 + 1 4 1 + y 2 1 / 2 K 1 1 + y 2 1 + ζ 2 y 2 1 y 2 1 + ζ 2 y 2 + 1 + ζ 2 1 + y 2 1 1 + y 2 1 + ζ 2 1 + y 2 .
The second integral in Equation (9) by including t in Equation (A12) gets to the same result (but ultimately having the opposite sign) by introducing a new variable t = t b , and using the substitution t + b ζ y c o s ϕ in Equation (21).

References

  1. Casimir, H.B.G. On the attraction between two perfectly conducting plates. Proc. Kon. Ned. Akad. Wetensch. B 1948, 51, 793–795. Available online: https://dwc.knaw.nl/DL/publications/PU00018547.pdf (accessed on 25 November 2023).
  2. Lifshitz, E.M. The theory of molecular attractive forces between solids. Sov. Phys. JETP 1956, 2, 73–83. Available online: http://jetp.ras.ru/cgi-bin/e/index/e/2/1/p73?a=list (accessed on 25 November 2023).
  3. Yablonovitch, E. Accelerating reference frame for electromagnetic waves in a rapidly growing plasma: Unruh-Davies-Fulling-DeWitt radiation and the nonadiabatic Casimir effect. Phys. Rev. Lett. 1956, 62, 1742–1746. [Google Scholar] [CrossRef] [PubMed]
  4. Dodonov, V.V.; Klimov, A.B.; Man’ko, V.I. Nonstationary Casimir effect and oscillator energy level shift. Phys. Rev. Lett. A 1989, 142, 511–513. [Google Scholar] [CrossRef]
  5. Schwinger, J. Casimir energy for dielectrics. Proc. Nat. Acad. Sci. USA 1989, 89, 4091–4093. [Google Scholar] [CrossRef] [PubMed]
  6. Dodonov, V. Fifty years of the dynamic Casimir effect. Physics 2020, 2, 67–104. [Google Scholar] [CrossRef]
  7. Mostepanenko, V.M. Casimir puzzle and Casimir conundrum: Discovery and search for resolution. Universe 2021, 7, 84. [Google Scholar] [CrossRef]
  8. Reiche, D.; Intravaia, F.; Busch, K. Wading through the void: Exploring quantum friction and vacuum fluctuations. APL Photonics 2022, 7, 030902. [Google Scholar] [CrossRef]
  9. Volokitin, A.I.; Persson, B.N.J. Near-field radiation heat transfer and noncontact friction. Rev. Mod. Phys. 2007, 79, 1291–1329. [Google Scholar] [CrossRef]
  10. Milton, K.A.; Høye, J.S.; Brevik, I. The reality of Casimir friction. Symmetry 2016, 8, 29–83. [Google Scholar] [CrossRef]
  11. Pendry, J.B. Shearing the vacuum–quantum friction. J. Phys. C Condens. Matter 1997, 9, 10301–10320. [Google Scholar] [CrossRef]
  12. Stipe, B.C.; Stowe, T.D.; Kenny, T.W.; Rugar, D. Noncontact friction and force fluctuations between closely spaced bodies. Phys. Rev. Lett. 2001, 87, 096801. [Google Scholar] [CrossRef] [PubMed]
  13. Volokitin, A.I. Casimir friction force between a SiO2 probe and a graphene-coated SiO2 substrate. JETP Lett. 2016, 104, 534–539. [Google Scholar] [CrossRef]
  14. Volokitin, A.I. Effect of an electric field in the heat transfer between metals in the extreme near field. JETP Lett. 2019, 110, 749–754. [Google Scholar] [CrossRef]
  15. Viotti, L.; Farias, M.B.; Villar, P.I.; Lombardo, F.C. Thermal corrections to quantum friction and decoherence: A closed time-path approach to atom-surface interaction. Phys. Rev. D 2019, 99, 105005. [Google Scholar] [CrossRef]
  16. Farias, M.B.; Lombardo, F.C.; Soba, A.A.; Villar, P.I.; Decca, R.S. Towards detecting traces of non-contact quantum friction in the corrections of the accumulated geometric phase. npj Quantum Inf. 2020, 6, 25. [Google Scholar] [CrossRef]
  17. Lombardo, F.C.; Decca, R.S.; Viotti, L.; Villar, P.I. Detectable signature of quantum friction on a sliding particle in vacuum. Adv. Quant. Technol. 2021, 4, 2000155. [Google Scholar] [CrossRef]
  18. Gurudev Dutt, M.V.; Childress, L.; Jiang, L.; Togan, E.; Maze, J.; Jelezko, F.; Zibrov, A.S.; Hemmer, P.R.; Lukin, D. Quantum register based on individual electronic and nuclear spin qubits in diamond. Science 2007, 316, 1312. [Google Scholar] [CrossRef]
  19. Dedkov, G.V. Low-temperature increase in the van der Waals friction force with the relative motion of metal plates. JETP Lett. 2021, 114, 779–784. [Google Scholar] [CrossRef]
  20. Dedkov, G.V. Puzzling low-temperature behavior of the van der Waals friction force between metallic plates in relative motion. Universe 2021, 7, 427. [Google Scholar] [CrossRef]
  21. Dedkov, G.V. Nonequilibrium Casimir-Lifshitz friction force and anomalous radiation heating of a small particle. Appl. Phys. Lett. 2022, 121, 231603. [Google Scholar] [CrossRef]
  22. Reiche, D.; Intravaia, F.; Hsiang, J.-T.; Busch, K.; Hu, B.-L. Nonequilibrium thermodynamics of quantum friction. Phys. Rev. 2020, 102, 050203. [Google Scholar] [CrossRef]
  23. Oelschläger, M.; Reiche, D.; Egerland, C.H.; Busch, K.; Intravaia, F. Electromagnetic viscosity in complex structured environments: From black-body to quantum friction. Phys. Rev. A 2022, 106, 052205. [Google Scholar] [CrossRef]
  24. Brevik, I.; Shapiro, B.; Silveirinha, M. Fluctuational electrodynamics in and out of equilibrium. Int. J. Mod. Phys. A 2022, 37, 2241012. [Google Scholar] [CrossRef]
  25. Intravaia, F. How modes shape Casimir physics. Int. J. Mod. Phys. A 2022, 37, 2241014. [Google Scholar] [CrossRef]
  26. Jentschura, U.D.; Lach, G.; De Kiviet, M.; Pachucki, K. One-loop dominance in the imaginary part of the polarizability: Application to blackbody and noncontact van der Waals friction. Phys. Rev. Lett. 2015, 114, 043001. [Google Scholar] [CrossRef]
  27. Klatt, J.; Bennett, R.; Buhmann, S.Y. Spectroscopic signatures of quantum friction. Phys. Rev. A 2016, 94, 063803. [Google Scholar] [CrossRef]
  28. Guo, X.; Milton, K.A.; Kennedy, G.; McNulty, W.P.; Pourtolami, N.; Li, Y. Energetics of quantum friction: Field fluctuations. Phys. Rev. D 2021, 104, 116006. [Google Scholar] [CrossRef]
  29. Guo, X.; Milton, K.A.; Kennedy, G.; McNulty, W.P.; Pourtolami, N.; Li, Y. Energetics of quantum friction. II: Dipole fluctuations and field fluctuations. Phys. Rev. D 2022, 106, 016008. [Google Scholar] [CrossRef]
  30. Milton, K.A.; Guo, X.; Kennedy, G.; Pourtolami, N.; DelCol, D.M. Vacuum torque, propulsive forces, and anomalous tangential forces: Effects of nonreciprocal media out of thermal equilibrium. arXiv 2023, arXiv:2306.02197. [Google Scholar] [CrossRef]
  31. Polevoi, V.G. Tangential molecular forces caused between moving bodies by a fluctuating electromagnetic field. Sov. Phys. JETP 1990, 71, 1119–1124. Available online: http://jetp.ras.ru/cgi-bin/e/index/e/71/6/p1119?a=list (accessed on 25 November 2023).
  32. Dedkov, G.V.; Kyasov, A.A. Friction and radiative heat exchange in a system of two parallel plate moving sideways: Levin-Polevoi-Rytov theory revisited. Chin. Phys. 2018, 56, 3002. [Google Scholar] [CrossRef]
  33. Condon, E.U.; Odishaw, H. (Eds.) Handbook of Physics; McGraw-Hill Book Company, Inc.: New York, NY, USA, 1967. [Google Scholar]
  34. Baptiste, J. Resistivity of gold. In The Physics Factbook; Elert, G., Ed.; 2004; Available online: https://hypertextbook.com/facts/2004/JennelleBaptiste.shtml (accessed on 25 November 2023).
  35. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series and Products; Academic Press/Elsevier Inc.: Waltham, MA, USA, 2014. [Google Scholar] [CrossRef]
  36. Grigoriev, I.S.; Meilikhov, E.Z. (Eds.) Handbook of Physical Quantities; CRC Press/Taylor & Francis Group: Boca Raton, FL, USA, 1996; Available online: https://archive.org/details/handbookofphysic0000unse_b1e4 (accessed on 25 November 2023).
  37. Biehs, S.-A.; Kittel, A.; Ben-Abdallah, P. Fundamental limitations of the mode temperature concept in strongly coupled system. Z. Naturforsch. A 2020, 75, 803–807. [Google Scholar] [CrossRef]
  38. Viloria, M.G.; Guo, Y.; Merabia, S.; Ben-Abdallah, P.; Messina, R. Role of Nottingham effect in the heat transfer in extreme near field regime. Phys. Rev. 2023, 107, 125414. [Google Scholar] [CrossRef]
  39. Pendry, J.B.; Sasihithlu, K.; Kraster, R.V. Phonon-assisted heat transfer between vacuum-separated surfaces. Phys. Rev. B 2016, 94, 075414. [Google Scholar] [CrossRef]
  40. Sasihithlu, K.; Pendry, J.B. Van der Waals force assisted heat transfer. Z. Naturforsch. A 2017, 72, 181–188. [Google Scholar] [CrossRef]
  41. Kuehn, S.; Loring, R.F.; Marohn, J.A. Dielectric fluctuations and the origins of noncontact friction. Phys. Rev. 2006, 96, 156103. [Google Scholar] [CrossRef]
  42. Dedkov, G.V. Casimir-Lifshitz friction force and kinetics of radiative heat transfer between metal plates in relative motion. JETP Lett. 2023, 117, 952–957. [Google Scholar] [CrossRef]
Figure 1. Configuration of parallel plates in relative motion. See text for details.
Figure 1. Configuration of parallel plates in relative motion. See text for details.
Physics 06 00002 g001
Figure 2. Resistivity of gold [34]. To obtain resistivity in the Gaussian units, one should use the relation Ω · m = 1 / 9 10 9 s . The ‘zero RR’ stands for zero residual resistance, T D is the Debye temperature, and ρ ~ T denotes a linear fit. See text for more details.
Figure 2. Resistivity of gold [34]. To obtain resistivity in the Gaussian units, one should use the relation Ω · m = 1 / 9 10 9 s . The ‘zero RR’ stands for zero residual resistance, T D is the Debye temperature, and ρ ~ T denotes a linear fit. See text for more details.
Physics 06 00002 g002
Figure 3. Quantum friction force of the plates of gold as a function of the velocity of a moving plate 2: (a) residual resistance of gold corresponds to Bloch–Grüneisen (BG) model at T = 5   K ; (b) residual resistance corresponds to modified BG (MBG) model at T = 0 (see Figure 2). The red lines represent complete numerical integration in Equation (9), the green and blue lines are calculations using Formulas (20) and (24). The positions of the maxima and the corresponding velocities are shown, respectively, by vertical lines and the numbers indicated.
Figure 3. Quantum friction force of the plates of gold as a function of the velocity of a moving plate 2: (a) residual resistance of gold corresponds to Bloch–Grüneisen (BG) model at T = 5   K ; (b) residual resistance corresponds to modified BG (MBG) model at T = 0 (see Figure 2). The red lines represent complete numerical integration in Equation (9), the green and blue lines are calculations using Formulas (20) and (24). The positions of the maxima and the corresponding velocities are shown, respectively, by vertical lines and the numbers indicated.
Physics 06 00002 g003
Figure 4. Friction parameter of gold plates as a function of their quasiequilibrium temperature for (a) BG and (b) MBG models. The curves with symbols were calculated using complete numerical integration in Equation (9). Solid lines correspond to the calculations using Equation (25) with Equation (26) (green lines) and (28) (blue lines); in (a), the blue and green lines merge. The vertical numbered lines show the temperatures corresponding to the maxima of the curves.
Figure 4. Friction parameter of gold plates as a function of their quasiequilibrium temperature for (a) BG and (b) MBG models. The curves with symbols were calculated using complete numerical integration in Equation (9). Solid lines correspond to the calculations using Equation (25) with Equation (26) (green lines) and (28) (blue lines); in (a), the blue and green lines merge. The vertical numbered lines show the temperatures corresponding to the maxima of the curves.
Physics 06 00002 g004
Figure 5. Friction parameter of gold plates as a function of the velocity of plate 2 for the (a) BG and (b) MBG models. The solid lines represent the calculations using Equation (9), dashed lines—Equations (37) and (38). The red, blue, and green lines correspond to quasistationary temperatures of 5 K, 10 K, and 77 K, respectively, for both plates. The different temperature order of the curves in (b) is explained by a different sequence of parameters α = ν / T : α 10 < α 5 < α 77 . The plateau in the curves corresponds to the linear velocity dependence of the friction force.
Figure 5. Friction parameter of gold plates as a function of the velocity of plate 2 for the (a) BG and (b) MBG models. The solid lines represent the calculations using Equation (9), dashed lines—Equations (37) and (38). The red, blue, and green lines correspond to quasistationary temperatures of 5 K, 10 K, and 77 K, respectively, for both plates. The different temperature order of the curves in (b) is explained by a different sequence of parameters α = ν / T : α 10 < α 5 < α 77 . The plateau in the curves corresponds to the linear velocity dependence of the friction force.
Physics 06 00002 g005
Figure 6. Heating rate of plate 1 (a) and friction parameter η = F x / V (b) as a function of velocity V of plate 2 in the BG model. The temperature configuration for the plates are indicated on the curves as follows: e.g., 6/4 denotes T 1 = 6   K ,   T 2 = 4   K . Thermal configurations T 1 = 6   K ,     T 2 = 4   K and T 1 = 4   K ,     T 2 = 6   K have the same friction parameters, configurations T = 4 and 6   K correspond to a quasiequilibrium thermal mode. The data shown by open triangles (∆ (a)) and by open diamonds (◊ (b)) are multiplied by 3 (cf. [42], where all numerical data to be reduced by π times).
Figure 6. Heating rate of plate 1 (a) and friction parameter η = F x / V (b) as a function of velocity V of plate 2 in the BG model. The temperature configuration for the plates are indicated on the curves as follows: e.g., 6/4 denotes T 1 = 6   K ,   T 2 = 4   K . Thermal configurations T 1 = 6   K ,     T 2 = 4   K and T 1 = 4   K ,     T 2 = 6   K have the same friction parameters, configurations T = 4 and 6   K correspond to a quasiequilibrium thermal mode. The data shown by open triangles (∆ (a)) and by open diamonds (◊ (b)) are multiplied by 3 (cf. [42], where all numerical data to be reduced by π times).
Physics 06 00002 g006
Figure 7. Same as in Figure 6 but in the MBG model. No additional numerical factors for the data are used.
Figure 7. Same as in Figure 6 but in the MBG model. No additional numerical factors for the data are used.
Physics 06 00002 g007
Figure 8. Heating time of gold plates as a function of the velocity of plate 2 at h = 500 μm according to BG (a) and MBG (b) models. Two upper lines correspond to heating from 4 K to 8 K at a = 20 nm (crimpson) and a = 10 nm (blue), and two lower lines correspond to heating from 4 K to 5 K at a = 20 nm (green) and a = 10 nm (red).
Figure 8. Heating time of gold plates as a function of the velocity of plate 2 at h = 500 μm according to BG (a) and MBG (b) models. Two upper lines correspond to heating from 4 K to 8 K at a = 20 nm (crimpson) and a = 10 nm (blue), and two lower lines correspond to heating from 4 K to 5 K at a = 20 nm (green) and a = 10 nm (red).
Physics 06 00002 g008
Figure 9. A possible setup for measuring Casimir–Lifshitz friction force (side view). The thermal protection layer is shown in blue, the metal coating is shown in brown. When the upper disk rotates, the circular sections of disks locating at a distance a move at a linear velocity of 0.5 Ω D relative to each other. At rotation frequencies n = 1–104 rps (revolutions per second) and disk diameter D = 0.1 m, the velocity range to be 0.3–3000 m/s.
Figure 9. A possible setup for measuring Casimir–Lifshitz friction force (side view). The thermal protection layer is shown in blue, the metal coating is shown in brown. When the upper disk rotates, the circular sections of disks locating at a distance a move at a linear velocity of 0.5 Ω D relative to each other. At rotation frequencies n = 1–104 rps (revolutions per second) and disk diameter D = 0.1 m, the velocity range to be 0.3–3000 m/s.
Physics 06 00002 g009
Table 1. Friction parameter η (in kg/(m2 · s)) of gold plates for velocity V = 1 m/s at thermal quasiequilibrium, Equation (9).
Table 1. Friction parameter η (in kg/(m2 · s)) of gold plates for velocity V = 1 m/s at thermal quasiequilibrium, Equation (9).
Temperature of Plates, Ka = 10 nma = 20 nma = 10 nma = 20 nm
Model BGModel MBG
14.81 × 10−62.77 × 10−65.60 × 10−82.80 × 10−8
22.63 × 10−41.47 × 10−41.22 × 10−75.98 × 10−8
31.10 × 10−35.73 × 10−41.87 × 10−79.13 × 10−8
53.44 × 10−41.67 × 10−43.08 × 10−71.50 × 10−7
102.15 × 10−51.04 × 10−55.63 × 10−72.73 × 10−7
154.30 × 10−52.09 × 10−67.77 × 10−73.76 × 10−7
201.52 × 10−67.35 × 10−77.19 × 10−73.33 × 10−7
502.04 × 10−79.90 × 10−82.56 × 10−71.25 × 10−7
1001.30 × 10−76.30 × 10−81.77 × 10−78.63 × 10−8
2001.14 × 10−75.54 × 10−81.42 × 10−76.81 × 10−8
3001.11 × 10−75.41 × 10−81.39 × 10−76.81 × 10−8
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dedkov, G.V. Casimir–Lifshitz Frictional Heating in a System of Parallel Metallic Plates. Physics 2024, 6, 13-30. https://doi.org/10.3390/physics6010002

AMA Style

Dedkov GV. Casimir–Lifshitz Frictional Heating in a System of Parallel Metallic Plates. Physics. 2024; 6(1):13-30. https://doi.org/10.3390/physics6010002

Chicago/Turabian Style

Dedkov, George V. 2024. "Casimir–Lifshitz Frictional Heating in a System of Parallel Metallic Plates" Physics 6, no. 1: 13-30. https://doi.org/10.3390/physics6010002

APA Style

Dedkov, G. V. (2024). Casimir–Lifshitz Frictional Heating in a System of Parallel Metallic Plates. Physics, 6(1), 13-30. https://doi.org/10.3390/physics6010002

Article Metrics

Back to TopTop