Nonequilibrium Casimir-Polder Interaction Between Nanoparticles and Substrates Coated with Gapped Graphene

The out-of-thermal-equilibrium Casimir-Polder force between nanoparticles and dielectric substrates coated with gapped graphene is considered in the framework of the Dirac model using the formalism of the polarization tensor. This is an example of physical phenomena violating the time-reversal symmetry. After presenting the main points of the used formalism, we calculate two contributions to the Casimir-Polder force acting on a nanoparticle on the source side of a fused silica glass substrate coated with gapped graphene, which is either cooler or hotter than the environment. The total nonequilibrium force magnitudes are computed as a function of separation for different values of the energy gap and compared with those from an uncoated plate and with the equilibrium force in the presence of graphene coating. According to our results, the presence of a substrate increases the magnitude of the nonequlibrium force. The force magnitude becomes larger with higher and smaller with lower temperature of the graphene-coated substrate as compared to the equilibrium force at the environmental temperature. It is shown that with increasing energy gap the magnitude of the nonequilibrium force becomes smaller, and the graphene coating makes a lesser impact on the force acting on a nanoparticle from the uncoated substrate. Possible applications of the obtained results are discussed.


Introduction
Much use is made of the physical systems which parts are kept at different temperatures, i.e., in out-of-thermal-equilibrium conditions.The processes occurring in thermally nonequilibrium systems are usually irreversable, i.e., they lack of symmetry inherent in systems characterized by some unique temperature equal to that of the environment.The theoretical description of thermally nonequilibrium systems is a rather complicated subject.However, the demands of both fundamental physics and its applications in nanotechnology lent impetus to a successful search for new theoretical approaches.
The Casimir-Polder interaction is an example on how the standard theory of a thermally equilibrium phenomenon can be extended to the out-of-thermal-equilibrium conditions.Casimir and Polder [1] derived an expression for the force acting between a small polarizable particle and an ideal metal plane kept at zero temperature.In the framework of the Lifshitz theory [2], this result was generalized to the case when an ideal metal plane is replaced with thick material plate kept at the same temperature T E as the environment [3,4].The Casimir-Polder force is a thermally equilibrium quantum phenomenon determined by the zero-point and thermal fluctuations of the electromagnetic field.It finds extensive applications in atomic physics and condensed matter physics (see the monographs [5][6][7][8] for a large body of examples).Specifically, much attention is given to the Casimir-Polder interaction between nanoparticles and material surfaces of different nature including biomembranes [9][10][11][12][13][14][15][16][17][18].
The extension of the Lifshitz theory to physical systems out of thermal equilibrium, which covers both cases of two macroscopic bodies and a small particle near the surface of a macroscopic body, was developed in [19][20][21][22][23][24] and further elaborated in [25][26][27][28][29][30][31][32].This was done under an assumption that each body is in the state of local thermal equilibrium [22].The developed theory was confirmed experimentally by measuring the Casimir-Polder force between 87 Rb atoms and a fused silica glass plate heated up to the temperatures much higher than T E [33].The case when the dielectric properties of a plate depend heavily on the temperature offers a fertile field for the study of the nonequilibrium Casimir-Polder force.This is true the for metallic plates [28,30] and for the plates made of a material, which undergoes the phase transition with increasing temperature [29].
The material, which demonstrates a profound effect of temperature on its dielectric properties, is graphene, i.e., a one-atom-thick layer of carbon atoms packed in the hexagonal lattice [34].At low energies, graphene is well-described by the Dirac model [35,36].Due to its relatively simple structure, the spatially nonlocal dielectric properties of graphene can be expressed via the polarization tensor in (2+1) dimensions [37,38] and found on the solid basis of quantum electrodynamics at any temperature [39][40][41][42].This presents the immediate possibility to calculate the Casimir-Polder force between nanoparticles and either heated or cooled graphene sheets using the extension of the Lifshits theory to the out-of-thermal-equilibrium situations.Calculations of this kind are of prime interest not only for fundamental physics, but for numerous applications harnessing interaction of nanoparticles with graphene and graphene-coated substrates as well (see, e.g., [43][44][45][46][47][48]).
The investigation of the nonequilibrium Casimir-Polder interaction between nanoparticles and graphene originated in [49], where the freestanding in vacuum pristine graphene sheet was considered.The pristine character of graphene means that there is no energy gap in the spectrum of quasiparticles, i.e., they are massless, and that it possesses the perfect crystal lattice with no foreign atoms.Both these assumptions are in the basis of the original Dirac model of graphene [34][35][36].According to [49], the nanoparticles have the same temperature as the environment, whereas the graphene sheet may be either cooler or hotter than the environment, which takes the nanoparticle-graphene system out of the state of thermal equilibrium.It was shown [49] that the nonequilibrium conditions have strong effect on the nanoparticle-graphene force and can even change the Casimir-Polder attraction with repulsion if the temperature of a graphene sheet is lower than T E .
Taking into account that real graphene sheets are characterized by some energy gap, i.e., small but nonzero mass of quasiparticles [34], an impact of the energy gap on the nonequilibrium Casimir-Polder force between a nanoparticle and a freestanding gapped graphene sheet was investigated in [50].It was shown that by varying the energy gap it is possible to control the force value.Unlike the case of a pristine graphene, the force acting on a nanoparticle from the sheet of gapped graphene preserves its attractive character.
When employing graphene sheets in physical experiments or in micro-and nanodevices, they are usually not freestanding in vacuum, but deposited on some dielectric substrate.Because of this, it is important to determine an impact of substrate underlying the gapped graphene sheet on the nonequilibrium Casimir-Polder force acting on a nanoparticle.In this article, we investigate the force on a spherical nanoparticle from the source side of dielectric substrate coated with gapped graphene, which is either cooled or heated as compared to the environmental temperature T E .As to the temperature of nanoparticles, it is assumed to be equal to T E .
After presenting the main points of the theoretical formalism, we continue with a comparison between the equilibrium Casimir-Polder forces on a nanoparticle from the uncoated and coated with gapped graphene dielectric substrate.Then, the nonequilibrium Casimir-Polder force acting on a nanoparticle from an uncoated substrate is considered.Next, we turn our attention to the calculation of two contributions of different nature to the nonequilibrium force on a nanoparticle from the graphene-coated substrates with various values of the energy gap of graphene coating kept at different temperatures.Finally, the total nonequilibrium force on a nanoparticle is studied as the function of separation at different temperatures of a graphene-coated substrate and typical values of the energy gap.
It is shown that the presence of a substrate increases the magnitude of the nonequilibrium force, which is larger and smaller than the equilibrium one for the graphene-substrate temperature higher and lower than T E , respectively.In all cases, an increase of the energy gap leads to smaller force magnitudes and to a lesser impact of the graphene coating on the nonequilibrium force on a nanoparticle from the uncoated dielectric substrate.
The article is organized as follows.Section 2 presents the necessary information from the theory of nonequilibrium Casimir-Polder force acting on nanoparticles from the substrates coated with gapped graphene.In Section 3, both the equilibrium and nonequilibrium forces on nanoparticles from the fused silica glass plate are considered and compared with the equilibrium one in the presence of graphene coating.In Section 4, the nonequilibrium Casimir-Polder force from a fused silica plate coated with gapped graphene is investigated.Section 5 contains the discussion.In Section 6, the reader will find our conclusions.

Theoretical Description of Nonequilibrium Casimir-Polder Interaction from Graphene-Coated Substrates
Here, we briefly present the formalism required for calculation of the Casimir-Polder force between nanoparticles and dielectric substrates coated with gapped graphene in out-of-thermal-equilibrium conditions.We assume that nanoparticles and the environment have the temperature T E (in computations we use T E = 300 K), whereas the graphenecoated substrate has the temperature T g , which may be either lower or higher than T E .The separation distance between nanoparticles and the plane surface of area S of a graphenecoated substrate is a.The energy gap of a graphene coating is denoted ∆.Below we consider the simplest case of spherical nanoparticles whose radius R satisfies the conditions where k B is the Boltzmann constant.It is suggested also that a ≪ √ S. Nanoparticles of arbitrary shape spaced at any separation from the surface, which should not be necessarily plane, can be considered using the more complicated scattering matrix approach [24,27].
It is convenient to represent the nonequilibrium Casimir-Polder force as a sum of two contributions [22,24] Here, F M somewhat resembles the equilibrium Casimir-Polder force expressed as a sum over the discrete Matsubara frequencies, whereas F r is the proper nonequilibrium contribution.
Actually, both terms in (2) describe some part of the effects of nonequilibrium and their explicit forms depend on which of the items responsible for these effects are included in F M and which in F r .Below we employ F M and F r in the forms derived in [49] and used in [50].Thus, F M has the form where α(0) is the static polarizability of a nanoparticle, k is the component of the wave vector parallel to the plane of graphene, k = |k|, the prime on the sum divides the term with l = 0 by 2, and the Matsubara frequencies and q l are defined as The quantities R sub TM,TE are the reflection coefficients of the electromagnetic waves on the graphene-coated substrate for the transverse magnetic and transverse electric polarizations.In fact (3) has the same form as the Lifshitz formula for the equilibrium Casimir-Polder force between a nanoparticle and a plate, but with one important difference.Here, the Matsubara frequencies are defined at the environmental temperature T E , but the reflection coefficients -at the temperature of graphene-coated substrate T g (in the Lifshitz formula T E = T g ).Note also that (3) contains the static polarizability α(0) in front of the summation sign in l, whereas in the Lifshitz formula for an atom-wall interaction α(iξ l ) appears under the sum in l.This is because the dynamic polarizability of nanoparticles at several first Matsubara frequencies contributing to the force under the condition (1) reduces to α(iξ l ) ≈ α(0) [26].
The proper nonequilibrium contribution F r in ( 2) is given by [49,50] where It is presented as an integral over the real frequency axis.In so doing, only k > ω/c, i.e., only the evanescent waves, contribute to (5).
In the literature, the other forms of F M and F r have been considered, where F r contains contributions from both the evanescent and propagating waves [22].The advantage of ( 5) is a presence of the exponential factor with real q and negative power which secures the quick convergence of the integral.
The reflection coefficients on the graphene-coated substrate were expressed via the frequency-dependent dielectric permittivity of substrate ε(ω) and the polarization tensor of graphene where For an uncoated substrate Π 00 = Π = 0 and (7) transforms into the familiar Fresnel reflection coefficients.To obtain the reflection coefficients at the pure imaginary frequencies entering (3), one should put in ( 7) ω = iξ E,l .Here, we also assume that the dielectric permittivity of a substrate material does not depend on temperature, contrary to the polarization tensor of graphene.In the case of materials (metals, for instance) with the temperature-dependent dielectric permittivity, one should change ε(ω) in (7) with ε(ω, T g ).
Computations of the nonequilibrium Casimir-Polder force between nanoparticles and graphene-coated substrate can be performed by Equations ( 2), ( 3), (5), and (7) if one knows the explicit expressions for the components of the polarization tensor of graphene.For the gapped graphene, these components defined along both real and imaginary frequency axes were found in [41] and presented more specifically in [49,50] in terms of the variables ω and k.It should be noted, however, that the numerical computations of (5) along the real frequency axis for substrates coated with gapped graphene are much more involved than of (3) and should be performed using the convenient choice of dimensionless variables.That is why, below we introduce the necessary dimensionless quantities and present the results ( 5)-( 8) and the expressions for the dimensionless polarization tensor in these terms.
Let us define the following dimensionless quantities Using ( 9) and ( 10), the reflection coefficients ( 7) can be rewritten as Below in this section, we rewrite the results of [41,50] for the polarization tensor in the region of evanescent waves in terms of dimensionless variables ( 9) and (10).It is common to present the polarization tensor as the sum of two contributions The first terms on the right-hand side of (13) refer to the case of zero temperature, whereas the second have the meaning of the thermal corrections to them.The polarization tensor in the area of evanescent waves k > ω/c has different forms in the so-called plasmonic region [52] and in the region where v F ≈ c/300 is the Fermi velocity for graphene and ṽF = v F /c.We begin with the plasmonic region (14).Here, the polarization tensor is also defined differently depending on the fulfilment of some condition.Thus, if the following inequality is satisfied: the first contributions to the polarization tensor ( 13) are given by where α = e 2 /(hc) is the fine structure constant and Under the condition ( 16), the second contributions to (13) take the form where From ( 17) and (19), it is seen that under the inequality ( 16) the polarization tensor is real.
Next, we continue the consideration of the plasmonic region ( 14) but under the inequality opposite to ( 16) In this case, the first contributions to the polarization tensor (13) take the form where The second contributions to (13) under the inequality ( 21) take a more complicated form where Note that the polarization tensor in the plasmonic region under the inequality ( 21) is the complex quantity.The first contributions to (13) defined in (22) are complex due to (23).As to the second contributions to (13) defined in (24), they are complex due to the square root in the denominators of (24) which have an imaginary part in the z interval near z (−) .Much care must be taken to this interval in numerical computations.
We are coming now to the region (15) and rewrite the results of [41,50] for the polarization tensor in terms of the dimensionless quantities ( 9) and (10).In this region, the first contributions to (13) are where The second contributions to (13) in the region (15) are expressed as From ( 26) and ( 28) it is seen that in the region (15) the first contributions to the polarization tensor, Π 00 and Π (0) , are real whereas the second, Π 00 and Π (1) , are complex.In the end of this section, we rewrite the proper nonequilibrium contribution (5) to the Casimir-Polder force it terms of the dimensionless variables ( 9), (10), and ( 12) where and the reflection coefficients R sub TM,TE expressed via the dimensionless variables are contained in (11).

Equilibrium and Nonequilibrium Casimir-Polder Forces from Fused Silica Plate Compared to Equilibrium Force in the Presence of Graphene Coating
In this and in the next sections, we use the formalism presented in Section 2 for numerical computations of the Casimir-Polder force acting on a nanoparticle on the source side of the graphene-coated substrate.As the substrate material, we choose fused silica glass, SiO 2 , which is in often use both in applications of graphene in nanoelectronics [53][54][55][56] and in physical experiments on measuring the Casimir force [57,58].
For better understanding of the relative roles of substrate, of nonequilibrium effects, and of graphene coating, in this section we start from calculation of the equilibrium Casimir-Polder forces acting on a nanoparticles from the uncoated and coated with gapped graphene SiO 2 substrate.Next, the nonequilibrium force from an uncoated SiO 2 substrate will be computed.
The equilibrium Casimir-Polder force acing on a nanoparticle from both uncoated and graphene-coated SiO 2 substrates is given by (3), where we put T E = T g .If the substrate is uncoated, the reflection coefficients R sub TM,TE are given by (7) where ω = iξ E,l and Π 00 = Π = 0.If the substrate is coated with graphene, one should use (7) with T g = T E .In this case, the polarization tensor of graphene is given by expressions ( 26) and (28) found in the frequency region (15), where, after returning to the dimensional variables, one can immediately put ω = iξ E,l (the explicit expressions for the polarization tensor of graphene at the pure imaginary frequencies iξ E,l are contained in [59]).
Computations of the equilibrium force F eq and the contribution F M to the nonequilibrium force require the dielectric permittivity of SiO 2 along the imaginary frequency axis at our disposal, whereas computations of the contribution F r to the nonequilibrium force require ε SiO 2 (ω) along the axis of real frequencies.In Figure 1, we present the real (a) and imaginary (b) parts of the dielectric permittivity of SiO 2 along the real frequency axis plotted by the tabulated optical data for the complex index of refraction of SiO 2 [60].The data collected in [60] extend from 0.0025 to 2000 eV, whereas in Figure 1 the most important region for our computations is presented.The values of the dielectric permittivity along the imaginary frequency axis are obtained from Im ε SiO 2 by means of the Kramers-Kronig relations and can be found in many literature sources (see, e.g., [7]).Now we compute the equilibrium Casimir-Polder force (31) acting on a nanoparticle from either an uncoated or graphene-coated SiO 2 plate as a function of separation at the plate temperature T p = T E = 300 K. Keeping in mind that the force values strongly depend on separation, we normalize them to the Casimir-Polder force on a nanoparticle from the ideal metal plane at zero temperature [7] Taking into account that the Dirac model is applicable at energies below 3 eV [61], we consider the nanoparticle-plate separations exceeding 200 nm where the characteristic energy hω c is much less than 3 eV.From the above, we restrict the separation region by 2 µm, where the condition ( 1) is yet applicable.
In Figure 2, the computational results for F eq /F (0) CP at T p = T E = 300 K are shown as the function of separation by the three lines counted from bottom to top for an uncoated SiO 2 plate and for coated by a graphene sheet with the energy gap ∆ = 0.2 eV and 0.1 eV, respectively (in the two latter cases the temperature of graphene is T g = T p = T E ).As is seen in Figure 2, the presence of graphene coating increases the magnitude of the Casimir-Polder force acting on a nanoparticle.This increase, however, is larger for a graphene sheet with smaller energy gap.
The ratio of the equilibrium Casimir-Polder force acting on a nanoparticle from a plate at T p = T = 300 K to that from an ideal metal plane at T p = T E = 0 is shown as the function of separation by the three lines plotted for a uncoated SiO 2 plate and coated by a graphene sheet with the energy gap ∆ = 0.2 eV and 0.1 eV.
To obtain the absolute values of force acting on a nanoparticle one can use the data of Figure 2 in combination with the following expressions for the static polarizabilities: which are valid for dielectric and metallic nanoparticles, respectively [26].Thus, for SiO 2 one obtains ε(0) = 3.81 from Figure 1(a).Now we calculate the nonequilibrium Casimir-Polder force F SiO 2 neq acting on a nanoparticle from the uncoated SiO 2 plate.In so doing we assume that the plate temperature is either lower, T p = 77 K, or higher, T p = 500 K than the environmental temperature T E = 300 K.The computations are performed by Equations ( 2), (3), and (29) with the reflection coefficients (7) where Π 00 = Π = 0 in the absence of graphene coating.
Note that the computation of the proper nonequilibrium contribution F r to the nonequlibrium Casimir-Polder force, which is performed along the real frequency axis, is much more labor and time consuming, especially in the presence of graphene coating (see in the next section).These computations were realized by using a program written in the C++ program-ming language [50] at the Supercomputer Center of the Peter the Great Saint Petersburg Polytechnic University.
In Figure 3(a), we present the computational results for the ratio F SiO 2 neq /F (0) CP as the function of separation by the bottom and top lines for the plate temperatures T p = 77 K and 500 K, respectively.For comparison purposes, the equilibrium Casimir-Polder force ratio CP for a plate kept at the environmental temperature T p = T E = 300 K is shown by the middle line which reproduces the bottom line from Figure 2. In the inset, the region of small separations from 0.2 to 0.6 µm is shown on an enlarged scale for better visualization.a (µm) Figure 3.The ratio of the nonequilibrium Casimir-Polder force acting on a nanoparticle from a SiO 2 plate to the equilibrium one from an ideal metal plane (a) at T p = T E = 0 and (b) at T p = T E = 300 K (the classical limit) is shown as the function of separation.The bottom and top lines are for the plate temperatures T p = 77 K and 500 K, respectively.The middle lines demonstrate similar ratio when T p = T E = 300 K.In the inset, the region of short separations is shown on an enlarged scale.
From Figure 3(a), it is seen that the effects of nonequilibrium have a strong impact on the equilibrium Casimir-Polder force from a SiO 2 plate making its magnitude larger if T p > T E and smaller if T p < T E .From Figure 3(a) one can also see that at separations exceeding 1 µm the Casimir-Polder force (32) from an ideal metal plane poorly reproduces the dependence of the nonequilibrium force on the separation distance.Because of this, in the region from 1 to 2 µm, we consider the same computational results for F SiO 2  neq , but normalize them to the thermal Casimir-Polder force from an ideal metal plane at large separations (the so-called classical limit) [7] In Figure 3(b), the computational results for the ratio F SiO 2 neq /F cl are shown as the function of separation by the bottom and top lines for the plate temperatures T p = 77 K and 500 K, respectively.The middle line demonstrates the ratio F SiO 2 eq /F cl in the case of thermal equilibrium T p = T E = 300 K. Figure 3(b) confirms the conclusions already made from Figure 3(a).It is seen also that the effects of nonequilibrium due to higher temperature than the environmental one make a greater impact on the equilibrium force than those due to a decrease of temperature to below the environmental one.Further decrease of temperature to below 77 K leads to only an insignificant decrease of the magnitude of nonequilibrium Casimir-Polder force.

Nonequilibrium Casimir-Polder Force from Fused Silica Plate Coated with Gapped Graphene
In the previous section, we already considered the nonequilibrium Casimir-Polder force from a SiO 2 plate.The main difference of the case of a graphene-coated plate considered now is that the response of graphene to the electromagnetic field described by the polarization tensor strongly depends on temperature.This is not the case for a SiO 2 plate whose dielectric permittivity is temperature-independent.As a result, for an uncoated SiO 2 plate the effects of nonequilibrium influence the force only through the factor Θ in ( 5) defined in (6).
According to (2), the total nonequilibrium force is the sum of two contributions F M and F r .The first of them has the same form as the equilibrium Casimir-Polder force, whereas the second is the proper nonequilibrium contribution.In order to understand physics of the process, we consider each of them separately starting from F M (for an uncoated SiO 2 plate F M is equal to the equilibrium force at T = 300 K shown by the bottom line in Figure 2 for any temperature of the plate).
Computations of F M were performed by ( 3), (7), and the polarization tensor defined at the pure imaginary frequencies iξ E,i (see Section 3 and [59]).The computational results for the ratio F M /F (0) CP are presented in Figure 4(a,b) as the function of separation for (a) graphene coating with the energy gap ∆ = 0.1 eV and (b) graphene coating with ∆ = 0.2 eV by the bottom and top lines plotted for the graphene (and plate) temperatures T p = 77 K and 500 K.The middle line shows the ratio F eq /F (0) CP when the graphene temperature is equal to that of the environment, T g = T E = 300 K.
Unlike the case of an uncoated SiO 2 plate, in the presence of graphene coating F M is not equal to the equilibrium force at T g = T E = 300 K.As is seen in Figure 4(a), for the relatively small energy gap of the graphene coating ∆ = 0.1 eV, the quantity F M at T g = 500 K differs little from the equilibrium force at T g = T E = 300 K over the entire separation region from 0.2 to 2 µm.Thus, the relative deviation varies in this case from 1.21% at a = 0.2 µm to 1.14% at a = 2 µm, i.e., is almost independent on separation.
From Figure 4(a) it is also seen that both the middle and top lines grow with separation.This means that for a substrate coated by graphene with ∆ = 0.1 eV considerable contribution to both F M at T g = 500 K and F eq at T g = T E = 300 K is given by the term of the Lifshitz formula with l = 0.At the same time, for a graphene coating with ∆ = 0.1 eV at T g = 77 K the deviation of F M from F eq is much larger and it increases with increasing separation [see the bottom and middle lines in Figure 4  The comparison of Figure 4(a) plotted for a graphene coating with ∆ = 0.1 eV with Figure 4(b) plotted for ∆ = 0.1 eV shows that the values of F M computed at T g = 77 K almost coincide (the difference of 0.6% at a = 0.2 µm and 1% at a = 2 µm).This means that at the relatively low T g = 77 K the thermal corrections are rather small and depend only slightly on the value of ∆ in the range of separations from 0.2 to 2 µm.The comparison of the values of F M at T g = 500 K for graphene coatings with ∆ = 0.1 and 0.2 eV [top lines in Figures 4(a) and 4(b)] also demonstrates rather small deviations (1% at a = 0.2 µm and 1.9% at a = 2 µm).This is because at so high temperature the thermal corrections are rather large and have only a weak dependence on ∆.As to the equilibrium Casimir-Polder force at T g = T E = 300 K, it depends on ∆ more distinctly (compare with Figure 2).Note also that the magnitudes of F M at both 77 K and 500 K are larger than the magnitudes of F neq from an uncoated SiO 2 plate at the same respective temperatures.This is seen from the comparison of Figure 4(a,b) and Figure 3.
Next, we consider the second contribution, F r , to the nonequilibrium Casimir-Polder force (2) acting on a nanoparticle from a graphene-coated substrate.The numerical computations of F r were performed in the dimensionless variables by ( 11), (29) using the polarization tensor in ( 17), ( 19), ( 22), ( 24), (26), and (28).The computational results for the ratio F r /F (0) CP are presented in Figure 5 as the function of separation for (a) graphene-plate temperature T g = 77 K and (b) graphene-plate temperature T g = 500 K.The solid line in Figure 5(a) is plotted for a graphene coating with ∆ = 0.1 eV, whereas another (dashed) line presents the coinciding results for the graphene coating with ∆ = 0.2 eV and for the uncoated SiO 2 plate.The top and bottom solid lines in Figure 5(b) are for a graphene coating with ∆ = 0.1 and 0.2 eV, respectively, and the dashed line is for the uncoated SiO 2 plate.The form of presentation in Figure 5(a,b) provides a way to determine the comparative contributions of the regions ( 14) and ( 15) to F r .As is seen in Figure 5(a,b), the quantity F r , unlike F M , substantially depends on the value of ∆ at both T g = 77 and 500 K.The point is that for T g < T E the contribution to F r from the regions ( 14) and ( 15) are positive, i.e., decrease the force magnitude.The opposite situation occurs for T g > T E , i.e., the contribution of regions ( 14) and (15) to F r are negative leading to the increase of force magnitude.
For the nonequilibrium Casimir-Poder force acting on a nanoparticle from a freestanding in vacuum graphene sheet, different contributions to F r were investigated in [50].It was shown that for a graphene with ∆ = 0.2 eV at T g = 77 K it holds F r ≈ 0, whereas for ∆ = 0.1 eV the main contribution to F r is given by the region (15) which leads to the increase of the force magnitude.If the graphene with ∆ = 0.1 eV is heated up to T g = 500 K, at short separations (a ≲ 0.4 µm) the contribution from the region ( 15) is dominant.With increasing separation, however, the region (14) takes the main role leading to a minor increase of force magnitude.For a freestanding graphene sheet with ∆ = 0.2 eV, the main contribution to F r at all separations is given by the region (15) [50].
For a graphene-coated substrate one has a more complicated situation.The point is that for an uncoated SiO 2 plate it is just F r which determines a distinction of the nonequilibrium force from the equilibrium one.As s result, for an uncoated plate F r > 0 for T g < T E and F r < 0 for T g > T E .As is seen in Figure 5(a), for T g = 77 K < T E the graphene coating with ∆ = 0.2 eV does not influence the value of F r , which is fully determined by the properties of a substrate.However, the graphene coating with ∆ = 0.1 eV makes a sizable effect on F r .Now we are in a position to present the computational results for the total nonequilibrium Casimir-Polder force F neq between a nanoparticle and a graphene-coated substrate which is the sum of contributions F M given in Figure 4 and F r given in Figure 5.In a (µm) In each figure, the bottom, middle, and top lines are plotted for the graphene-plate temperature T g = 77 K, 300 K, and 500 K, respectively.In Figures 6(a  CP from an ideal metal plane (32), whereas in Figures 6(b) and 7(b) -to the classical limit of the Casimir-Polder force F cl from an ideal metal plane (34).The bottom and top lines are for the graphene-plate temperatures T g = 77 K and 500 K, respectively.The middle lines demonstrate similar ratio where T g = T E = 300 K.In the inset, the region of short separations is shown on an enlarged scale.
Thus, the middle lines in Figures 6(a) and 7(a) reproduce the top and middle lines in Figure 2, respectively, which are also plotted for the equilibrium Casimir-Polder force from the graphene-coated substrate with ∆ = 0.1 and 0.2 eV.Figures 6(b) and 7(b) can be compared with Figures 1 and 2 in [50], plotted there for a freestanding in vacuum graphene sheet, in order to determine an impact of the SiO 2 substrate on the nonequilibrium force.By and large, Figures 6 and 7 are in analogy to Figure 3 related to the case of an uncoated substrate.
From Figures 6 and 7 it is seen that the nonequilibrium force F neq from the graphenecoated SiO 2 plate at both T g = 77 K and 500 K is larger in magnitude than the force F SiO 2 neq from an uncoated SiO 2 plate and from a freestanding graphene sheet.Thus, for a SiO 2 plate coated with a graphene sheet with ∆ = 0.1 eV at T g = 77 K the ratio F neq /F SiO 2  neq is equal to 1.05 and 0.94 at separations a = 0.2 and 2 µm, respectively.The same ratio for the same graphene coating but at T g = 500 K is equal to 1.12 and 1.33 at the same respective separations.For a graphene coating with ∆ = 0.2 eV at T g = 77 K, the ratio F neq /F SiO 2 neq is equal to 1.04 and 1.01 at a = 0.2 and 2 µm, and, if the graphene-coated plate is at the temperature T g = 500 K, -to 1.04 and 1.32, respectively.
If we compare the Casimir-Polder force from the graphene-coated SiO 2 plate with that from a freestanding graphene sheet, F free , it is seen that the presence of a substrate significantly increase the magnitudes of both equilibrium and nonequilibrium forces.This increase is the most pronounced at short separations.As an example, for a graphene coating on a SiO 2 plate and a freestanding graphene sheet with ∆ = 0.1 eV, we obtain F eq /F free eq = 4.3 and 1.4 at a = 0.2 and 2 µm, respectively.For the nonequilibrium forces with the same value of ∆, the ratio F neq /F free neq is equal to 8.2 and 2.2 at T g = 77 K and to 4.7 and 1.6 at T g = 500 K at the same respective separations.
The above results allow to control the value of the nonequilibrium Casimir-Polder force acting on nanoparticles from the graphene-coated substrates.

Discussion
The out-of-thermal-equilibrium Casimir and Casimir-Polder forces is a rather novel subject which is investigated only during the last 25 years.Despite of this, considerable advances have already been made in understanding the underlying physics and mathematical description of the effects of nonequilibrium.Specifically, as noted in Section 1, the fundamental Lifshitz theory of the Casimir and Casimir-Polder forces was generalized for the out-of-thermal-equilibrium conditions [19][20][21][22][23][24][25].At first this was made for the case when the material properties are temperature-independent, but in succeeding years generalized for materials whose response functions to the electromagnetic field explicitly depend on temperature as a parameter [29,30].
Graphene is a unique novel material whose response functions described by the polarization tensor essentially depend on the temperature.Because of this, the effects of nonequilibrium in the Casimir-Polder force acting on nanoparticles from a graphene sheet are best suited to both the theoretical study and practical utility in nanoelectronics.Previously these effects were considered only for a freestanding graphene sheets [49,50], which is a configuration not easily accessible in a laboratory.In this article, the nonequilibrium Casimir-Polder force on a nanoparticle from a substrate coated with gapped graphene sheet was investigated as the function of separation, temperature, and energy gap in the case of the most often used fused silica glass substrate.This opens a way to a practical implementation of this physical phenomenon in the field effect transistors and other devices of bioelectronics [62][63][64].The obtained results can be also generalized to other two-dimensional materials and van der Waals heterostructures employing different 2D crystals [65][66][67].

Conclusions
To conclude, in this article the formalism of the Lifshitz theory generalized to outof-thermal-equilibrium conditions with temperature-dependent material properties was used to investigate the nonequilibrium Casimir-Polder force between nanoparticles and either cooled or heated fused silica plate substrate coated with gapped graphene sheet.The response of graphene to the electromagnetic field was described by the polarization tensor in the framework of the Dirac model.
We investigated two different contributions to the nonequilibrium Casimir-Polder force and determined their physical meaning and relative role depending on the frequency and separation regions.The total nonequilibrium force from a graphene-coated fused silica glass substrate was compared with the equilibrium one from the same source, as well as with both equilibrium and nonequilibrium forces from an uncoated silica glass plate.A comparison with the nonequilibrium Casimir-Polder force from a freestanding in vacuum graphene sheet has also been made.
It was shown that the nonequilibrium force from the graphene-coated silica glass substrate kept at both lower and higher temperature than in the environment is larger in magnitude than the nonequilibrium force from an uncoated silica glass plate and from a freestanding graphene sheet.According to the results obtained, an increase of the energy gap of graphene coating leads to smaller force magnitudes and to a lesser impact of the graphene coating on the nonequilibrium force acting on a nanoparticle on the source side of an uncoated silica glass plate.By and large, we determined the impact of temperature of a graphene sheet, the role of a substrate and of the nonzero energy gap of graphene coating on the nonequilibrium Casimir-Polder force between a nanoparticle and a graphene-coated substrate spaced at different separations.
The above results may find application in the rapidly progressing areas of nanotechnology dealing with integrated nanoparticle-biomolecular systems.

Figure 1 .
Figure 1.Optical data for the (a) real and (b) imaginary parts of the dielectric permittivity of fused silica glass are shown as the function of frequency along the real frequency axis.

Figure 4 .
Figure 4.The ratio of the first contribution to the nonequilibrium Casimir-Polder force acting on a nanoparticle from a SiO 2 plate coated by a graphene sheet with the energy gap (a) ∆ = 0.1 eV and (b) ∆ = 0.2 eV to the equilibrium one from an ideal metal plane at T p = T E = 0 is shown as the function of separation.The bottom and top lines are for the graphene-plate temperatures T g = 77 K and 500 K, respectively.The middle lines demonstrate similar ratio where T g = T E = 300 K.

Figure 5 .
Figure 5.The ratio of the second contribution to the nonequilibrium Casimir-Polder force acting on a nanoparticle from a SiO 2 plate coated by a graphene sheet (a) at T g = 77 K and (b) T g = 500 K to the equilibrium one from an ideal metal plane at T p = T E = 0 is shown as the function of separation (a) by the two lines for ∆ = 0.1 and 0.2 eV, where the latter coincides with that for an uncoated plate shown by the dashed line, and (b) by the two solid lines for ∆ = 0.1 and 0.2 eV, where the latter coincides with that for an uncoated plate shown by the dashed line Figures 6(a,b) and 7(a,b) the results for F neq are shown for the graphene coating with ∆ = 0.1 and 0.2 eV, respectively.

Figure 6 .
Figure 6.The ratio of the nonequilibrium Casimir-Polder force acting on a nanoparticle from a SiO 2 plate coated by a graphene sheet with ∆ = 0.1 eV to the equilibrium one from an ideal metal plane (a) at T p = T E = 0 and (b) at T p = T E = 300 K (the classical limit) is shown as the function of separation.The bottom and top lines are for the graphene-plate temperatures T g = 77 K and 500 K, respectively.The middle lines demonstrate similar ratio where T g = T E = 300 K.
) and 7(a), the values of

F
neq are normalized to the zero-temperature Casimir-Polder force F (0)

Figure 7 .
Figure 7.The ratio of the nonequilibrium Casimir-Polder force acting on a nanoparticle from a SiO 2 plate coated by a graphene sheet with ∆ = 0.2 eV to the equilibrium one from an ideal metal plane (a) at T p = T E = 0 and (b) at T p = T E = 300 K (the classical limit) is shown as the function of separation.The bottom and top lines are for the graphene-plate temperatures T g = 77 K and 500 K, respectively.The middle lines demonstrate similar ratio where T g = T E = 300 K.In the inset, the region of short separations is shown on an enlarged scale.

Funding:
The work of O.Yu.T. was supported by the Russian Science Foundation under Grant No. 21-72-20029.G.L.K. was 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.