Nonequilibrium Casimir-Polder Force between Nanoparticles and Graphene-Coated Silica Plate: Combined Effect of the Chemical Potential and Mass Gap

The Casimir-Polder force between spherical nanoparticles and a graphene-coated silica plate is investigated in situations out of thermal equilibrium, i.e., with broken time-reversal symmetry. The response of graphene coating to the electromagnetic field is described on the basis of first principles of quantum electrodynamics at nonzero temperature using the formalism of the polarization tensor in the framework of the Dirac model. The nonequilibrium Casimir-Polder force is calculated as a function of the mass-gap parameter, chemical potential of graphene and temperature of the graphene-coated plate, which can be both higher and lower than that of the environment. It is shown that the force value increases with increasing chemical potential, and this increase is more pronounced when the temperature of a graphene-coated plate is lower than that of the environment. The nonequilibrium force also increases with increasing temperature of the graphene-coated plate. This increase is larger when the plate is hotter than the environment. The effect is revealed that the combined impact of the chemical potential $\mu$ and mass gap $\Delta$ of graphene coating depends on the relationship between $\Delta$ and 2$\mu$. If $2\mu>\Delta$ the magnitude of the nonequilibrium force between nanoparticles and a cooled graphene-coated plate becomes much larger than for a graphene coating with $\mu=0$. The physical reasons explaining this effect are elucidated. Possible applications of the obtained results are discussed.


Introduction
The two-dimensional (2D) hexagonal lattice of carbon atoms called graphene [1] has attracted particular attention of many experts in both fundamental and applied physics because at low energies it adequately fits the requirements of the Dirac model.What this means is that graphene can be considered as consisting of light quasiparticles described not by the Schrödinger equation, as is customary in condensed matter physics, but by the relativistic Dirac equation in 2D space and one-dimensional (1D) time, where the speed of light c is replaced with the Fermi velocity v F ≈ c/300 [2].
For experts in fundamental physics, graphene is interesting because it possesses unusual electrical, optical and mechanical properties which offer outstanding possibilities to large-scale exploratory investigations [3,4].At the same time, graphene suggests a wide range of prospective applications in nanotechnology and material science [5,6].Specifically, certain of the nanotechnological devices rely upon interaction between nanoparticles and graphene-coated plates for their operation [7][8][9][10][11][12].
The Casimir-Polder force [13] acting between an electrically neutral microscopic particle and a material plate is the relativistic and quantum phenomenon.It is initiated by the zero-point and thermal fluctuations of the electromagnetic field.In the state of thermal equilibrium, when temperatures of the microparticle, plate, and of the environment are equal, the Casimir-Polder force is described by the Lifshitz theory [14][15][16].This theory expresses the force value via the dynamic polarizability of a particle and the frequency-dependent dielectric permittivity of the plate material.
There is an extensive literature devoted to calculations of the Casimir-Polder force between atoms, molecules, nanoparticles and diversified material surfaces, including graphene, by means of the Lifshitz theory (see, e.g., Refs.).In doing so, the spatially nonlocal dielectric permittivities and the 2D reflection coefficients of graphene sheets [45] were expressed via the polarization tensor of graphene in (2+1) dimensional space-time [46][47][48][49] in the framework of the Dirac model.
In some instances, the condition of thermal equilibrium is violated with the result that the Lifshitz theory becomes inapplicable.For example, this may happen in the applications of graphene to solar cells where the graphene-coated plate has different temperature from that of nanoparticles and of the environment [50,51].Although almost all physical theories, including the Lifshitz theory, are symmetric relatively to the operation of time reversal, in thermally nonequilibrium systems this symmetry is broken.
The generalization of the Lifshitz theory to the cases when the condition of thermal equilibrium is violated was developed under an assumption that each body in the physical system under consideration is in the state of local thermal equilibrium [52][53][54][55].As a result, it has become possible to calculate the Casimir-Polder force in situations when the temperature of a plate, a nanoparticle, an even of the environment are dissimilar.
The theory of nonequilibrium Casimir and Casimir-Polder forces was further generalized to the case of nonplanar configurations [56][57][58][59][60] and for the plate materials with temperature-dependent dielectric permittivities [61][62][63][64][65].It was also confirmed experimentally by measuring the nonequilibrium Casimir-Polder force between the Bose-Einstein condensate of 87 Rb atoms and a SiO 2 plate, which was heated as compared to the environmental temperature [66].
Investigation of the nonequilibrium Casimir-Polder interaction between spherical nanoparticles and graphene sheet using the formalism of the polarization tensor has been started by Ref. [67].In that work, the pristine graphene sheet freestanding in vacuum was considered.This means that the crystal lattice of graphene was assumed to be perfect with no foreign atoms and gapless quasiparticle spectrum.According to the results of Ref. [67], the impact of nonequilibrium conditions on the total Casimir-Polder force is the most pronounced at short separations between a nanoparticle and a graphene sheet.At these separations, the force may even change its sign and become repulsive.
Real graphene sheets are, however, not pristine and, specifically, the quasiparticles in graphene, although light, are not massless.This leads to some nonzero mass gap ∆ in the quasiparticle spectrum [3,68,69].An impact of the nonzero mass gap on the nonequilibrium Casimir-Polder interaction between nanoparticles and a freestanding in vacuum graphene sheet was investigated in Ref. [70] within the formalism of the polarization tensor.It was shown that in this case the nonequilibrium Casimir-Polder force remains attractive and it is possible to control its value by varying the mass-gap parameter.
The freestanding graphene sheets are difficult, if not impossible, to use in physical experiments and in nanodevices, where they are usually deposited on some dielectric plates.An impact of the plate material on the nonequilibrium Casimir-Polder force between nanoparticles and a silica plate coated with gapped graphene was considered in Ref. [71].It was shown that the presence of a silica plate increases the magnitude of the nonequilibrium Casimir-Polder force.With increasing mass-gap parameter, the magnitude of the nonequilibrium force decreases and the impact of graphene coating on the force becomes smaller.
It should be noted that computations of the nonequilibrium Casimir-Polder force acting on nanoparticles from real graphene sheet using the polarization tensor are rather cumbersome.They are connected with calculation of the multiple integrals of rapidly varying functions and were realized on a supercomputer [70,71].
In addition to the nonzero mass gap, real graphene sheets are unavoidably doped, i.e., contain some fraction of foreign atoms described by the nonzero value of the chemical potential µ [3,72,73].This should be taken into account in the theoretical computations intended for the comparison with the measurement data as was successfully done in the experiment on measuring the equilibrium Casimir force acting between an Au-coated sphere and a graphene-coated SiO 2 plate [74,75].
The nonequilibrium Casimir interaction between two plates coated with the gapless graphene sheets possessing the nonzero chemical potential was investigated in Ref. [76] using the phenomenological, spatially local, model for the response of graphene to the electromagnetic field [72,77].
In this paper, we investigate the Casimir-Polder force acting between spherical nanoparticles and fused silica plate coated with gapped and doped graphene sheet in situations out of thermal equilibrium using the formalism of the polarization tensor.It is assumed that the temperature of nanoparticles is the same as of the environment.However, the temperature of the graphene-coated plate can be either lower or higher than that of the environment.The spatially nonlocal response of graphene to the fluctuating electromagnetic field is described by the polarization tensor, which takes into account the dependence of graphene properties on temperature, on the mass-gap parameter, and on the chemical potential.In doing so, the explicit expressions for the polarization tensor along the real frequency axis are presented in the region of evanescent waves.The obtained results for the nonequilibrium Casimir-Polder force are compared with those for the equilibrium one in the same configuration.
It is shown that increasing of the chemical potential of graphene coating leads to the monotonous increase of both the equilibrium and nonequilibrium Casimir-Polder forces between a nanoparticle and a graphene-coated plate.In doing so, for a heated graphenecoated plate, the mass gap makes a lesser impact on the force if the chemical potential takes a nonzero value.According to our results, for a cooled graphene-coated plate the impact of chemical potential of graphene coating on the nonequilibrium Casimir-Polder force essentially depends on whether the inequality 2µ > ∆ or 2µ < ∆ is satisfied.By and large, the impact of chemical potential on the nonequilibrium force for a cooled graphene-coated plate is stronger than for a heated one.
The paper is organized as follows.In Section 2, both the formalism of the Lifshitz theory and the computational results for the equilibrium Casimir-Polder force between nanoparticles and a silica plate coated with gapped and doped graphene sheet are presented.Section 3 contains the generalization of the Lifshitz theory of Casimir-Polder force for the out-of-thermal-equilibrium conditions and the required expressions for the polarization tensor of graphene along the real frequency axis with due regard for different relationships between the mass gap and chemical potential of graphene coating.In Section 4, the computational results for the nonequilibrium Casimir-Polder force as a function of separation are presented under different values of temperature, mass gap and chemical potential of graphene coating.In Section 5, the reader will find a discussion of the obtained results.Section 6 contains our conclusions.

Equilibrium Casimir-Polder Force between Nanoparticles and Silica Plate Coated with Gapped and Doped Graphene
We consider the spherical nanoparticles of radius R spaced above large, graphenecoated fused silica (SiO 2 ) plate at height a ≫ R. The graphene coating is characterized by the mass gap ∆ and chemical potential µ.The fused silica plate chosen for computations is the most typical substrate used as a supporter for graphene [74,75,[78][79][80][81].All below equations are, however, applicable to any material substrate.
The temperature of the environment is T E = 300 K.In this section, the temperature of graphene-coated plate, T p , is equal to T E : T p = T E = 300 K. Nanoparticle temperature is assumed to be equal to T E throughout the paper.We also assume that at all temperatures T considered below the nanoparticle radius satisfies the condition R ≪ hc/(k B T), where h is the reduced Planck constant, c is the speed of light, and k B is the Boltzmann constant.At T = T E this condition reduces to the inequality R ≪ 7.6 µm.When it is satisfied, the dynamic polarizability of nanoparticle α(ω) at all frequencies ω contributing to the Casimir-Polder force is approximately equal to the static value α(0) [59].
In the situation of thermal equilibrium T p = T E , the Casimir-Polder force between a nanoparticle and a graphene-coated plate is given by the Lifshitz formula written in terms of the discrete Matsubara frequencies [14][15][16]82,83].Keeping in mind application of the same notations in the next section devoted to the situation out of thermal equilibrium, in this section we continue to denote the plate and environment temperatures as T p and T E but imply that T p = T E = 300 K.Then, the equilibrium Casimir-Polder force between a nanoparticle and a graphene-coated SiO 2 plate is given by [14][15][16]82,83] Here, k is the magnitude of the wave vector projection on the plane of graphene-coated plate, the Matsubara frequencies are ξ E,l = 2πk B T E l/h with l = 0, 1, 2, . . ., the prime on the summation sign means that the term with l = 0 is divided by 2, and The reflection coefficients on the graphene-coated plate for the electromagnetic fluctuations with transverse magnetic (TM) and transverse electric (TE) polarizations are expressed via the dielectric permittivity of the plate material ε(ω) and the polarization tensor of graphene Π βγ (ω, k, ∆, µ, T p ) calculated at the pure imaginary Matsubara frequencies [40,74,75,84] Here, and the following combination of the components of the polarization tensor is introduced: where a summation over the repeated index β = 0, 1, 2 is implied.
For numerical computations of the equilibrium Casimir-Polder force (1) between a nanoparticle and a graphene-coated substrate, one needs the values of the static polarizability of a nanoparticle, α(0), of the dielectric permittivity of plate material, ε(iξ E,l ), and of the components of the polarization tensor of graphene, Π 00 (iξ E,l , k, ∆, µ, T P ) and Π(iξ E,l , k, ∆, µ, T P ).
The dielectric permittivity of SiO 2 at the pure imaginary Matsubara frequencies is obtained from the optical data for the imaginary part of ε(ω) [85] with the help of the Kramers-Kronig relation.The obtained values were repeatedly used in calculations of the Casimir and Casimir-Polder forces [82].There are also analytic expressions for the dielectric permittivity of SiO 2 along the imaginary frequency axis [86,87].
Below we present the components of the polarization tensor of gapped and doped graphene at the Matsubara frequencies found in the framework of the Dirac model [46][47][48][49].As mentioned in Section 1, this model is applicable at low energies.In Ref. [88] the upper boundary of its application region is estimated as 3 eV.Thus, one should consider not too small separations a between the nanoparticle and the graphene-coated plate in order that all energies giving the major contribution to the Casimir-Polder interaction lie below this boundary.
It has been known that the characteristic energy of the Casimir-Polder interaction is equal to hc/(2a) [82,89].This energy is less than 1 eV at all separation distances a > 100 nm.That is why, at separations exceeding, e.g., 200 nm one can reliably use the polarization tensor of graphene derived in the framework of the Dirac model.This conclusion was confirmed by measurements of the Casimir force between an Au-coated sphere and a graphene-coated SiO 2 plate, which were found in a very good agreement with theoretical predictions using the polarization tensor of graphene [74,75].
The explicit expression for the component Π 00 of the polarization tensor is given by [40,74,75] where α = e 2 /(hc) is the fine structure constant and the following notations are introduced In a similar way, the explicit expression for the combination of the components of the polarization tensor Π takes the form [40,74,75] Computations of the equilibrium Casimir-Polder force between nanoparticles and graphene-coated SiO 2 plate were performed by Eqs. ( 1), ( 3), (6), and (8).Since the force values quickly decrease with increasing separation, they were normalized either to the Casimir-Polder force F 0 between a nanoparticle and an ideal metal plane at zero temperature or to the classical limit F cl reached in the same configuration, considered at temperature T E at large separations [82,89] In Figure 1, the computational results for (a) F SiO 2 eq /F 0 and (b) F SiO 2 eq /F cl are presented as the function of separation for a graphene coating with the mass gap ∆ = 0.2 eV and chemical potential µ = 0, 0.075, and 0.15 eV (lines 1, 2, and 3, respectively).As is seen in Figure 1, the force magnitude increases with increasing chemical potential as it should be.At all separations considered the magnitude of F SiO 2 eq decreases with separation slower than eq | decreases with separation faster than |F cl |.Next, let us consider the impact of the mass-gap parameter on the Casimir-Polder force.In Figure 2, the computational results for (a) F SiO 2  eq /F 0 and (b) F SiO 2 eq /F cl are presented as the function of separation by the pairs solid and dashed lines labeled 1 and 2. In pair 1, the mass-gap parameter ∆ = 0.2 eV for both lines, whereas the chemical potential µ = 0.075 eV for the solid and 0 for the dashed line.In pair 2, the mass-gap parameter ∆ = 0.1 eV for both lines with the same values of µ, as in pair 1, for the solid and dashed lines.
As is seen in Figure 2, the force magnitude increases with decreasing value of the energy gap.This means that the chemical potential and the energy gap of graphene coating act on the Casimir-Polder force in the opposite directions.It is seen also that if µ ̸ = 0 (see the solid lines) the dependence of the force value on ∆ is much weaker than for graphene coating with µ = 0 (see the dashed lines).
Note that numerical results presented in Figures 1 and 2 do not depend on α(0), i.e., on the nanoparticle radius.The absolute values of F SiO 2 eq can be obtained from Figures 1 and 2 by using Eq. ( 9) and the values of α(0) indicated above.
Figure 2. The equilibrium Casimir-Polder force between a nanoparticle and a graphene-coated SiO 2 plate at T p = T E = 300 K normalized to (a) the zero-temperature Casimir-Polder force from an ideal metal plane and (b) classical limit of the same force at T E = 300 K is shown as the function of separation by the two pairs solid and dashed lines.In the pair labeled 1, the mass-gap parameter ∆ = 0.2 eV for both lines, whereas the chemical potential µ = 0.075 eV for the solid line and 0 for the dashed line.In the pair labeled 2, ∆ = 0.1 eV for both lines with the same corresponding values of µ as in the pair 1.

Nonequilibrium Casimir-Polder Force between Nanoparticle and Plate Coated with Gapped and Doped Graphene: General Formalism
Now we consider the situation out of thermal equilibrium when the temperature of nanoparticles is equal to T E , i.e., is the same as of the environment, whereas the temperature of the graphene-coated plate T p can be either lower or higher than T E .
In the generalization of the Lifshitz theory to out-of-thermal-equilibrium conditions, the nonequilibrium force is usually presented as the sum of two terms, one of which is akin the equilibrium force by its form and another one is a truly nonequilibrium term [55,57].In the literature, different forms of these terms are contained leading to the same total result.Below, we represent the nonequilibrium Casimir-Polder force following Refs.[67,70,71] Here, F SiO 2 eq was already defined in Eq. ( 1) as an equilibrium Casimir-Polder force, but now it has another meaning.In Eq. ( 1), it was assumed that T p = T E = 300 K, whereas now T p ̸ = T E and, thus, F SiO 2  eq , although resembles the equilibrium force by its form, describes some part of the effects of nonequilibrium.
The truly nonequilibrium contribution to Eq. ( 10) is given by [67,70,71] where The reflection coefficients in Eq. ( 11) are similar to those in Eq. ( 3), but now they are defined along the real frequency axis where It is seen that Eq. ( 13) is obtained from Eq. ( 3) simply by putting iξ E,l = ω.However, the components of the polarization tensor ( 6) and ( 8) are rather complicated functions.They should be analytically continued from the imaginary to real frequency axis of the plane of complex frequency which is a nontrivial procedure.
The advantage of representation (10) is that the truly nonequilibrium contribution to it (11) contains an integration on only the evanescent waves.For these waves k > ω/c and the quantity q(ω, k) in the power of the exponent under the integral is real.
The analytic continuation of the polarization tensor to the real frequency axis under the condition k > ω/c is, nevertheless, rather involved.It takes different forms in the so-called plasmonic region [90] and in the region We begin with the plasmonic region (15).The analytic continuation of the polarization tensor to this region of the real frequency axis for a pristine and gapped graphene was performed in Refs.[67,70,71].Here, we present the obtained results in a more convenient analytic form and generalize them for the case of gapped graphene with nonzero chemical potential.
As was shown in Ref. [48], the form of the sought for analytic continuation depends on whether hcp < ∆ or hcp ⩾ ∆ where p ≡ p(ω, k) = ω 2 − v 2 F k 2 /c.If the inequality hcp(ω, k) < ∆ is satisfied, the results of Ref. [48] for Π 00 generalized for the case µ ̸ = 0 can be identically represented in the following form: where In the region hcp(ω, k) < ∆ under consideration A(ω, k, ∆) < 0 and the quantity Π 00 is real.
Using the same notations, for Π one obtains from Ref. [48] Π where The quantity Π is also the real function.
Now we consider the region of ω and k plane where the inequality hcp(ω, k) ⩾ ∆ is satisfied.In this region, A(ω, k, ∆) ⩾ 0 and, as a consequence, the quantities Π 00 and Π have both the real and imaginary parts.From the results of Ref. [48] generalized for the case µ ̸ = 0, after identical transformations, one obtains where and the roots of the common denominator of the functions B 1 and B 2 are These roots belonging to the integration domain D < u 1 < u 2 < ∞ are obtained for λ = −1.
In a similar manner, Π is given by Note that the imaginary parts of Π 00 and Π arise in the second integrals in Eqs. ( 21) and ( 24) when integrating with respect to u from u 1 to u 2 .
Let us continue and consider the region (16).In this region, the analytic continuation of the polarization tensor ( 6) and ( 8) to the real frequency axis can be obtained by putting iξ E,l = ω with the appropriate choice of the branch of the square root in denominator [48].For this purpose, it is first convenient to represent the real part in Eqs. ( 6) and ( 8) as half a sum of the complex conjugated quantities.Then, the quantity Π 00 in the region (16) of real frequency axis takes the form where, in accordance with Eq. ( 7), In a similar way, the analytic continuation of the quantity Π to the region ( 16) of real frequency axis is Let us discuss the obtained expressions for the polarization tensor of gapped and doped graphene along the real frequency axis.First of all let us note that the polarization tensor behaves differently depending on whether the inequality ∆ > 2µ or ∆ ⩽ 2µ is satisfied.
If ∆ > 2µ the powers of exponents in the definitions of w and w in Eqs. ( 18) and ( 26) are positive even at the lower integration limits u = D and u = D in Eqs. ( 17) and ( 19), ( 21) and ( 24), ( 25) and ( 27).This is true because from Eqs. ( 18) and ( 26) one obtains In this case, Eqs.( 18) and ( 26) lead to lim As a result, the polarization tensor at zero temperature defined as takes a much simplified form.It is equal to the first terms in Eqs. ( 17) and ( 19), in the plasmonic region (15) under the condition hcp < ∆, to the first terms in Eqs. ( 21) and ( 24), in the plasmonic region (15) under the condition hcp ⩾ ∆, and to the first terms in Eqs. ( 25) and ( 27), in the region (16).Thus, if the inequality ∆ > 2µ is satisfied, the polarization tensor at zero temperature (30) does not depend on the chemical potential µ.
It is easily seen that under the opposite inequality ∆ < 2µ there are the intervals in the integration domains where the powers of exponents in the definitions of w and w in Eqs. ( 18) and ( 26) with κ = −1 are negative.Thus, in the plasmonic region (15) the power of exponent in w with κ = −1 is negative within the interval In a similar way, in the region ( 16) the power of exponent in w with κ = −1 is negative in the interval In both cases it holds lim As a result, not only the first but all terms in Eqs. ( 17) and ( 19), ( 21) and ( 24), ( 25) and ( 27) contribute to the polarization tensor at zero temperature defined in Eq. ( 30).By way of example, in the region (16) the polarization tensor of graphene with the mass gap and chemical potential satisfying the condition ∆ < 2µ takes the following form at T = 0: where the upper integration limit ũ0 is defined in Eq. (35).Thus, under the condition ∆ < 2µ the polarization tensor at zero temperature depends on the chemical potential µ through the terms in the figure brackets in Eq. ( 37) and similar terms surviving at T = 0 in Eqs. ( 17), ( 19) and ( 21), (24) valid in the plasmonic region.
Note that in the plasmonic region (15) the expression for the polarization tensor at zero temperature (32), found under the opposite condition 2µ < ∆ ⩽ hcp contains the pure imaginary part for some values of ω and k satisfying the inequality hcp ⩾ ∆.This is in line with the fact that additional contributions to the zero-temperature polarization tensor arising under the condition ∆ < 2µ also contain the imaginary parts.It is only a quantitative difference.
Another situation arises in the region (16) where k > ω/v F .Here, under the condition ∆ > 2µ, the polarization tensor at zero temperature (33) is real, but under the condition ∆ < 2µ it gains some imaginary part depending on µ for some values of ω and k.This leads to important qualitative differences between the cases ∆ > 2µ and ∆ < 2µ, which are demonstrated in Section 4 by the results of numerical computations.

Computational Results for the Nonequilibrium Casimir-Polder Force
In this section, the nonequilibrium Casimir-Polder force F SiO 2 neq from Eq. ( 10) between nanoparticles and a silica plate coated with a gapped and doped graphene sheet is computed using the formalism presented in Sections 2 and 3.As in Section 2, the temperatures of nanoparticles and of the environment are assumed to be equal to T E = 300 K, whereas the temperature of the graphene-coated plate is either lower, T p = 77 K, or higher, T p = 500 K, than that of the environment.All computations were performed for typical values of the mass-gap parameter ∆ = 0.1 and 0.2 eV and the chemical potential µ = 0, 0.075 and 0.15 eV of the graphene coating within the separation region from 200 nm to 2 µm.
Numerical computations of the first contribution F SiO 2 eq to F SiO 2 neq given by Eq. ( 1) were performed along the imaginary frequency axis as discussed in Section 2 but with T p ̸ = T E .The computational results for the second contribution F SiO 2 r to F SiO 2 neq were obtained by Eq. ( 11) using the polarization tensor of graphene along the real frequency axis defined for different values of ω and k in Eqs. ( 17) and ( 19), ( 21) and ( 24), ( 25) and (27).Numerical data for the real and imaginary parts of the dielectric permittivity of SiO 2 are taken from Ref. [85] (see also Ref. [71] for the graphical representation of these data).Note that numerical computations of F SiO 2  neq are much more complicated than those for F SiO 2 eq made along the imaginary frequency axis owing to the quick variation of the integrands.They were performed using the program utilizing the Gauss-Kronrod and double-exponential quadrature methods for numerical integration from the GNU scientific Library [91] and Boost C ++ Libraries [92].The program was written in the C ++ programming language.It utilizes the OpenMP Library [93] for parallelism.The numerical computations were performed on a supercomputer of the Peter the Great Saint Petersburg Polytechnic University.Sufficiently high precision of the computational results was reached using the Boost Multiprecision Library [94].
We begin with the case when the SiO 2 plate coated with gapped and doped graphene is heated up to 500 K whereas the temperature of nanoparticles and of the environment here and below is always T E = 300 K.The mass-gap parameter of graphene coating is ∆ = 0.2 eV, whereas the chemical potential can be equal to zero or take the values µ = 0.075 eV and 0.15 eV.In Figure 3 the computational results for the nonequilibrium Casimir-Polder force normalized to (a) the zero-temperature Casimir-Polder force F 0 acting on a nanoparticle from an ideal metal plane and (b) the classical limit F cl of the same force at T p = T E = 300 K [see Eq. ( 9)] are shown as the function of separation by the lines labeled 1, 2, and 3 for graphene coating with µ = 0, 0.075, and 0.15 eV, respectively.
From Figure 3 it is seen that with increasing chemical potential the magnitude of the nonequilibrium Casimir-Polder force increases.Thus, in this respect the nonequilibrium force behaves in the same way as the equilibrium one.Note that the lines 1 and 2, for which 2µ < ∆ = 0.2 eV, are rather close to each other at all separations.The line 3, for which 2µ > ∆, is further apart from the lines 1 and 2 at short separations but approach them at larger separations.Thus, for a heated graphene-coated plate there is no qualitative difference between the cases 2µ < ∆ and 2µ > ∆.Now we consider the situation when the graphene-coated plate is cooled to T p = 77 K preserving unchanged all the other parameters of graphene coating listed in Figure 3.
The computational results for the normalized nonequilibrium Casimir-Polder force are presented as a function of separation in Figure 4 using the same notations for all lines as in Figure 3.
As is seen in Figure 4, the case of a cooled graphene-coated plate differ greatly from the case of a heated one.In the case of graphene-coated plate kept at T p = 77 K, an impact of the chemical potential on the nonequilibrium Casimir-Polder force essentially depends on the relationship between 2µ and ∆.If 2µ < ∆ (this holds for the line labeled 2) the impact of nonzero µ on the force value is rather moderate.If, however, 2µ > ∆ (line 3), this leads to significant deviation of the force value from that computed at µ = 0 (line 1) which increases with increasing separation between nanoparticles and a graphene-coated plate.The enhancement of nonequilibrium force for the graphene coating with 2µ > ∆ is explained by the increasing role of the polarization tensor at T = 0, which has a nonzero imaginary part in this case, and decrease of the thermal correction to it (see the discussion in the end of Section 3).By and large, although the nonequilibrium Casimir-Polder force increases with increasing µ for both the heated and cooled graphene-coated plate, in the latter case the impact of µ on the force value is much stronger.
We consider now the SiO 2 plate coated with the graphene sheet possessing the smaller mass-gap parameter ∆ = 0.1 eV and the moderate chemical potential µ = 0.075 eV.In this case, the condition 2µ > ∆ is satisfied.The nonequilibrium Casimir-Polder force acting on nanoparticles in such a configuration was computed at three different temperatures of the graphene-coated plate: T p = 77 K (the cooled plate), T p = 300 K (the situation of the thermal equilibrium), and T p = 500 K (the heated plate).
The computational results for the Casimir-Polder force normalized to (a) the force F 0 from an ideal metal plane at zero temperature and (b) the classical force F cl from an ideal metal plane at T p = T E = 300 K are shown in Figure 5 as the function of separation by the lines labeled 1, 2, and 3 for temperatures of the graphene-coated plate T p = 77 K, 300 K, and 500 K, respectively.As is seen in Figure 5, the Casimir-Polder force at all separations increases with increasing temperature.This increase is smaller when the temperature of the graphene-coated plate T p is below the temperature of the environment T E and larger when T P > T E .In doing so, the impact of temperature on the force value increases with increasing separation.
It is interesting to consider the role of the mass-gap parameter.For this purpose, the same configuration at the same three temperatures of the graphene-coated plate was considered with the larger mass-gap parameter of the graphene sheet ∆ = 0.2 eV.In this case the opposite inequality 2µ = 0.15 eV < ∆ is satisfied.The computational results for the nonequilibrium (lines 1 and 3) and equilibrium (line 2) Casimir-Polder force normalized to (a) F 0 and (b) F cl are shown on Figure 6 by the lines labeled 1, 2, and 3 for temperatures of the graphene-coated plate T p = 77 K, 300 K, and 500 K, respectively.Similar to Figure 5, the Casimir-Polder force increases with increasing temperature and this increase becomes larger at larger separations.By comparing the force values in Figure 5 and Figure 6, it is seen that at 77 K (the lines labeled 1 in both figures) the value of the nonequilibrium Casimir-Polder force for the graphene sheet with ∆ = 0.1 eV (Figure 5) is much larger than that for the graphene sheet with ∆ = 0.2 eV (Figure 6).This is because under the condition 2µ > ∆ satisfied for the configuration considered in Figure 5 the impact of the polarization tensor at T = 0 is significantly increased as discussed in the end of Section 3. As to the force values shown by the lines 2 and 3, in Figure 6, where 2µ < ∆, they are only slightly smaller than those in Figure 5.
Finally, we consider the same configuration, as in Figure 6, i.e., the SiO 2 plate coated by a graphene sheet with the mass-gap parameter ∆ = 0.2 eV, but with the twice as large chemical potential µ = 0.15 eV.In this case, the condition 2µ > ∆ is satisfied at the cost of sufficiently large µ, but not at the cost of relatively small ∆ as in Figure 5.The computational results normalized to (a) F 0 and (b) F cl are shown on Figure 7 as the function of separation by the lines 1 and 3 for the nonequilibrium Casimir-Polder force computed at T p = 77 K and 500 K, respectively, and by the line 2 for the equilibrium force computed at T p = 300 K. Similar to the previous figures, the force values in Figure 7 increase with increasing temperature, and this increase is more pronounced at larger separations between nanoparticles and a graphene-coated plate.There is only a minor impact of the increased chemical potential on the lines 2 and 3 relevant to the equilibrium situation and to the heated graphene-coated plate, respectively.These lines demonstrate only slightly larger force values in comparison with those in Figure 6.One can conclude that at sufficiently high temperature the polarization tensor is almost independent on whether 2µ > ∆ or 2µ < ∆ although the zero-temperature contribution and the thermal correction to it are different under these conditions.As to the line 1 in Figure 7, computed at T p = 77 K, it demonstrates much larger values of the nonequilibrium Casimir-Polder force than the line 1 in Figure 6.This effect is explained by the fact that at 2µ > ∆ the polarization tensor at zero temperature takes the relatively large value.Specifically, it has a nonzero imaginary part in the region (16) (see a discussion in the end of Section 3).For a cooled plate just this region contributes to the increase of the force value.

Discussion
As was shown in Section 1, the Casimir-Polder interaction of nanoparticles with material substrates and, especially, with graphene-coated plates attracts much attention in fundamental physics and finds applications in nanotechnology.Because of this, it is important to reliably predict the strength of this interaction and its dependence on all relevant parameters.Graphene is a novel material and, to a first approximation, its response to the electromagnetic field was described by using different phenomenological approaches, such as Kubo theory, 2D Drude model etc. (see Refs. [2][3][4]) for a review).
Development of the formalism of the polarization tensor [46][47][48][49] has made it possible to calculate the equilibrium Casimir and Casimir-Polder forces in graphene systems in the framework of the Dirac model on the basis of first principles of thermal quantum field theory at nonzero temperature.However, application of the same formalism to the nonequilibrium forces, although rather straightforward theoretically, presents serious computational difficulties.The reason is that the equilibrium Casimir and Casimir-Polder forces are calculated as the sum over the discrete pure imaginary Matsubara frequencies, whereas computation of the nonequilibrium forces is unavoidably involved with the integration of rapidly varying functions along the real frequency axis [52][53][54][55][56][57][58][59][60][61][62][63], which was always challenging in the Casimir physics.
This difficulty was overcome step by step.At first, the nonequilibrium Casimir-Polder force between nanoparticles and a freestanding sheet of a pristine graphene was considered [67].This case is interesting from the theoretical point of view, but does not reflect realistic experimental situations.Next, the impact of the mass gap, which exists in real graphene sheets, on the nonequilibrium force was investigated [70].Among other things, the mass gap necessary arises in graphene sheets deposited on some substrate.Because of this, the obtained results were generalized for the case of a graphene-coated dielectric plate [71] as usually takes place in experimental situations.All these computations demanded the use of a supercomputer.
One more important parameter, which was as yet disregarded, is the chemical potential reflecting the concentration of foreign atoms in a graphene sheet.As was shown previously, just the relationship between the mass gap and chemical potential determines the behavior of the Casimir and Casimir-Polder forces in graphene systems [40,42,44].In Ref. [76] the impact of chemical potential on the nonequilibrium Casimir force between two parallel graphene-coated plates was investigated, but with disregarded mass gap of graphene and using the spatially local model of its electromagnetic response.
The formalism and computational results presented in this paper take into account all the parameters important for an adequate description of the nonequilibrium Casimir-Polder interaction between nanoparticles and a graphene-coated plate in the experimental situation, i.e., the mass gap and chemical potential of a graphene sheet, the dielectric permittivity of the plate material and their temperature.As expected from the experience obtained when studying the graphene systems in thermal equilibrium, the nonequilibrium Casimir-Polder force between nanoparticles and a graphene-coated substrate essentially depends on the relationship between the energy gap and twice the chemical potential of the graphene coating.This opens opportunities for high-reliability predictions of the Casimir-Polder forces in graphene systems out of thermal equilibrium for a subsequent comparison with the experimental data.

Conclusions
To conclude, in the foregoing we have investigated the Casimir-Polder force between spherical nanoparticles and a graphene-coated SiO 2 plate in out-of-thermal-equilibrium conditions.The response of graphene to the electromagnetic field was described by the spatially nonlocal polarization tensor taking into account both the mass gap and chemical potential of the graphene coating, which was found in the framework of the Dirac model.The plate material was characterized by the frequency-dependent dielectric permittivity obtained from the tabulated optical data for the complex index of refraction of silica glass.
By computing the nonequilibrium Casimir-Polder force on a supercomputer, the following results were obtained.With increasing chemical potential µ of the graphene coating, the nonequilibrium Casimir-Polder force increases for both heated and cooled graphene-coated plate.An impact of the chemical potential on the force value is stronger when the graphene-coated plate is cooled, as compared to the environment temperature, and weaker when it is heated.An impact of the mass gap ∆ of the graphene coating is also more pronounced for a cooled graphene-coated plate.In this case, the character of the nonequilibrium Casimir-Polder force essentially depends on the relationship between the values of ∆ and 2µ.An impact of the chemical potential on the force value is rather moderate if 2µ < ∆, but if 2µ > ∆ the magnitude of the nonequilibrium force becomes much larger than that computed for a graphene coating with µ = 0.This effect is explained by the increasing role of the polarization tensor of graphene defined at zero temperature.
The magnitude of the nonequilibrium Casimir-Polder force increases with increasing temperature of the graphene-coated plate, all other conditions being equal.This increase is weaker and stronger if the temperature of the graphene-coated plate is lower and higher than that of the environment, respectively.In all cases the impact of temperature on the force value becomes stronger with increasing separation between nanoparticles and a graphene-coated plate.
To conclude, the revealed dependencies of the nonequilibrium Casimir-Polder force between nanoparticles and a graphene-coated plate from the mass-gap parameter, chemical potential, and temperature provide a way to control the force value in the graphene-based nanotechnological devices of next generations.

Figure 1 .
Figure1.The equilibrium Casimir-Polder force between a nanoparticle and a graphene-coated SiO 2 plate at T p = T E = 300 K normalized to (a) the zero-temperature Casimir-Polder force from an ideal metal plane and (b) classical limit of the same force at T E = 300 K is shown as the function of separation for a graphene coating with the mass gap ∆ = 0.2 eV and chemical potential µ = 0, 0.075, and 0.15 eV by the lines labeled 1, 2, and 3, respectively.

Figure 3 .
Figure3.The nonequilibrium Casimir-Polder force between a nanoparticle and a graphene-coated SiO 2 plate kept at T p = 500 K normalized to (a) the zero-temperature Casimir-Polder force from an ideal metal plane and (b) classical limit of the same force at T E = 300 K is shown as the function of separation for a graphene coating with the mass gap ∆ = 0.2 eV and chemical potential µ = 0, 0.075, and 0.15 eV by the lines labeled 1, 2, and 3, respectively.

Figure 4 .
Figure 4.The nonequilibrium Casimir-Polder force between a nanoparticle and a graphene-coated SiO 2 plate at T p = 77 K normalized to (a) the zero-temperature Casimir-Polder force from an ideal metal plane and (b) classical limit of the same force at T E = 300 K is shown as the function of separation for a graphene coating with the mass gap ∆ = 0.2 eV and chemical potential µ = 0, 0.075, and 0.15 eV by the lines labeled 1, 2, and 3, respectively.

Figure 5 .
Figure5.The Casimir-Polder force between a nanoparticle and a SiO 2 plate coated by a graphene sheet with the mass-gap parameter ∆ = 0.1 eV and chemical potential µ = 0.075 eV normalized to (a) the Casimir-Polder force from an ideal metal plane at zero temperature and (b) classical limit of the same force at T E = 300 K is shown as the function of separation for a graphene-coated plate temperatures T p = 77 K (thermal nonequilibrium), T p = 300 K (thermal equilibrium), and T p = 500 K (thermal nonequilibrium) by the lines 1, 2, and 3, respectively.

Figure 6 .
Figure6.The Casimir-Polder force between a nanoparticle and a SiO 2 plate coated by a graphene sheet with the mass-gap parameter ∆ = 0.2 eV and chemical potential µ = 0.075 eV normalized to (a) the Casimir-Polder force from an ideal metal plane at zero temperature and (b) classical limit of the same force at T E = 300 K is shown as the function of separation for a graphene-coated plate temperatures T p = 77 K (thermal nonequilibrium), T p = 300 K (thermal equilibrium), and T p = 500 K (thermal nonequilibrium) by the lines 1, 2, and 3, respectively.

Figure 7 .
Figure7.The Casimir-Polder force between a nanoparticle and a SiO 2 plate coated by a graphene sheet with the mass-gap parameter ∆ = 0.2 eV and chemical potential µ = 0.15 eV normalized to (a) the Casimir-Polder force from an ideal metal plane at zero temperature and (b) classical limit of the same force at T E = 300 K is shown as the function of separation for a graphene-coated plate temperatures T p = 77 K (thermal nonequilibrium), T p = 300 K (thermal equilibrium), and T p = 500 K (thermal nonequilibrium) by the lines 1, 2, and 3, respectively.

Funding:
G.L.K.and V.M.M. were partially funded by the Ministry of Science and Higher Education of Russian Federation ("The World-Class Research Center: Advanced Digital Technologies," contract No. 075-15-2022-311 dated April 20, 2022).The research of V.M.M. was partially carried out in accordance with the Strategic Academic Leadership Program "Priority 2030" of the Kazan Federal University.