FluctuationInduced Interactions Between Ellipsoidal Particle and Planar Substrate Immersed in Critical Medium
Abstract
In our study we investigate the behaviour of the net force (NF) emerging between an ellipsoidal particle and a thick plate (slab), when the interaction takes place in a near critical fluid medium with account for the omnipresent van der Waals forces (vdWF). Here we consider the case of complete wetting of the objects surfaces by the medium, due to strong adsorbing local surface potentials, exerted by thin solid coating films. The influence of the bulk inner regions of the particle and the slab on the constituents of the fluid results in longranged competing dispersion potentials. As a consequence from the critical fluctuations of the medium, the system experiences an additional effective interaction, traditionally termed critical Casimir force (CCF). The forces of interest are evaluated numerically from integral expressions obtained utilizing general scaling arguments and meanfield type calculations in combination with the socalled ”surface integration approach” (SIA). Within the scenario considered here, this technique is applicable if one has knowledge of the forces between two parallel semiinfinite plates, confining in between some fluctuating fluid medium characterized by its temperature and chemical potential . It is demonstrated that for a suitable set of particlefluid, slabfluid, and fluidfluid coupling parameters the competition between the effects due to the coatings and the core regions of the objects result, when one changes or , in sign change of the NF acting between the ellipsoid and the slab.
aff1]G. S. Valchev aff1]D. M. Dantchev
Introduction
When a confined fluid medium is at the vicinity of its bulk critical point, as first suggested by Fisher and de Gennes [1], an additional component adds up to the already acting forces if any, resulting from the correlated critical fluctuation of the medium. This new contribution is of a longranged character and strongly depends on the boundary conditions which the confining objects, immersed in the medium, impose on it at their surfaces. Since this fluctuationinduced force depends additionally only on some gross features of the fluid medium [2, 3] it can be treated as universal in nature. One can make a parallel with the force between neutral bodies due to the quantum and temperature fluctuations of the charge distributions in them, i.e., of the electromagnetic field, which force is known today under the general name of quantum electrodynamic (QED) Casimir force [4, 5]. Because of this analogy it became customary to term the fluctuation part of the net force critical Casimir force [2, 3].
In a medium where both quantum and thermodynamical fluctuations are present, can be thought as a sum of a regular background and a singular contributions. The first depends in an analytic way on the parameters characterizing the medium, whereas the latter emerges due to the critical fluctuations of the medium, respectively. Therefore
(0) 
It is normal to equate
(0) 
since the dispersion van der Waals interactions, ubiquitous for any system, are not influenced by the thermodynamic fluctuations.
Following Refs. [6, 7, 8], here we are going to be interested in the simultaneous manifestation of both the critical Casimir force (CCF) and the van der Waals force (vdWF), which gives rise to the net force between an ellipsoidal in shape particle and a thick smooth plate, say the handling arm of a microgripper – Figure 1(b). It will be demonstrated that by proper choice of the materials (cores) of the colloid particle and the plate it is indeed possible to achieve control over the net force (NF) by simply changing the temperature and chemical potential of the fluid medium. So far there is only one article [9], which addresses the geometry of the particle studied here. Moreover, the interactions accounted there are of pure shortrange type.
The content of the article is arranged as follows. In the following subsection we recall and comment on the finitesize behaviour of the CCF, vdWF and NF when they act between pair of parallel plates. Here we also identify the main coupling parameters characterizing the dispersion interactions in the systems, and give the general expressions used to calculate the CCF and NF. Next, the ”surface integration approach” (SIA), within which we study the commented spectrum of forces, is introduced. The exact equations used to evaluate numerically the investigated forces are derived in the last section, which ends with a discussion on the observed behaviour of the NF.
Influence of the dispersion forces on the thermodynamic Casimir one in nonpolar fluid systems with film geometry
If an object of certain shape is introduced in a medium, it alters its thermodynamic behaviour in such a way that any quantitative occurrence shows dependance on the penetration depth of the symmetry breaking effect into the volume. The range to which this effect is felt within the system depends both on the scope of the interactions and on that of the correlations between the fluctuations in the fluid, which mediate the interactions between the bodies.
The range of the correlations is set by the correlation length of the socalled order parameter of the medium. This quantity becomes large, and theoretically diverges, at the vicinity of the bulk critical point : , where , and , where . Here and are the usual critical exponents which, for classical fluids, are those of the threedimensional Ising model, and and are the corresponding nonuniversal amplitudes of the correlation length along the and axes. When becomes comparable to the characteristic dimension of the system, say the separation between the objects, the size dependence of the thermodynamic functions enters into the thermodynamic potentials through the ratio , and takes a scaling form given by the finitesize scaling theory [3].
Let us consider the system depicted on Figure 1(a). Upon approaching the vicinity of the bulk critical point (critical region) of , following Ref. [6] for the occurring force , per cross section area and , is customary to write the following expression
(0) 
Here is the dimensionality of the system, – an exponent describing the decay of the dispersion interactions, is a microscopic length introduced for dimensional reasons (see the text below Eq. (9) in Ref. [7]), is the Hamaker term, whose dependence from the temperature and chemical potential is given by the socalled Hamaker constant [10, 11] (for details see the Appendix in Ref. [6]), is dimensionless, universal scaling function, and are the temperature and field relevant (in renormalization group sense) scaling variables, whereas , and are irrelevant scaling variables (see the text below Eq. (2.6) in Ref. [6]) associated with the dispersion interactions in the system. The latter are introduced in the theory through the dimensionless coupling constants (see Section III in Refs. [6, 7])
(0) 
In the expressions shown in Eq. (Influence of the dispersion forces on the thermodynamic Casimir one in nonpolar fluid systems with film geometry), and designate the long and, respectively, the shortranged components of the Londonvan der Waals potentials, between the constituents of the core regions and the fluid and within the medium itself ; and are the number densities of and the critical one of , respectively, in units .
According to the scaling hypothesis of the CCF one expects that near the bulk critical point
(0) 
where is a scaling function, that for large enough with fixed and approaches the scaling function of the shortranged system (for details see Eqs. (2.12) and (4.10) in Ref. [6]). From Eqs. (INTRODUCTION) and (INTRODUCTION), together with Eq. (Influence of the dispersion forces on the thermodynamic Casimir one in nonpolar fluid systems with film geometry) follows that is proportional to the sum of and the singular part of the Hamaker term . The last implies that in order to determine the CCF in systems with dispersion interaction one has to decompose the contribution captured through the Hamaker term in a singular and a regular parts, i.e.
(0) 
Thus, for one has
(0) 
while equals to the second term of Eq. (Influence of the dispersion forces on the thermodynamic Casimir one in nonpolar fluid systems with film geometry), with .
In what follows, we are going to present results for the CCF, vdWF and NF in the cases of ellipsoidplate system, utilizing the knowledge gained from studies of the corresponding interactions between parallel plates.
The Surface Integration Approach and Its Application for the EllipsoidPlate System
In 1934 the soviet scientist B. Derjaguin was the first to propose an approach [12] for calculating geometry dependent interactions in systems where at least one of the objects has a nonplanar geometry. In particular, this approach focuses on relating the interaction force/potential between two gently curved colloidal particles with the knowledge for that between a pair of parallel plates . An important feature of this technique is that it is only applicable if the separation distance between the interacting objects is much smaller than their geometrical characteristics. In order to overcome this inconvenient condition the cocalled ”surface integration approach” (SIA) reported in Ref. [13] was developed. Here we must note that both DA and SIA are strictly valid if the interactions involved can be described by pair potentials.
Within this new technique the interaction force between an object (say a colloid particle) of arbitrary shape and a flat surface bounded by the plane of a Cartesian coordinate system, is determined by subtracting from the contributions stemming from the surface regions of the particle that ”face towards” the plane those from the regions that ”face away” from it, i.e.,
(0) 
Here and are the projections of the corresponding parts of the surface of the body on the plane, with . The expression Eq. (THE SURFACE INTEGRATION APPROACH AND ITS APPLICATION FOR THE ELLIPSOIDPLATE SYSTEM) also takes into account that the force on a given point of is along the normal to the surface at that point (for details see Section 2 in Ref. [13]). Therefore one can make use of the SIA to calculate the CCF, since under the assumption for mechanical equilibrium of the particle in the fluid, this force acts at any point of along the normal to the surface. Thus, one can get a reasonably good approximation to the effect of that force by keeping just the integration over the part of the surface of the body that faces the plane. This leads to
(0) 
In what follows we are going to make use of the SIA to estimate the vdWF and the CCF between a planar substrate and an ellipsoidal particle with halfradii and center positioned at a point with respect to the substrate’s surface, where and . The shape of the particle allows two cases to be considered:
 (i)

when the distance is fixed;
 (ii)

when the distance of closest approach between the surfaces of the interacting objects is fixed.
Here is the separation between the particle and the substrate when points toward the plane.
Hence, using Eq. (THE SURFACE INTEGRATION APPROACH AND ITS APPLICATION FOR THE ELLIPSOIDPLATE SYSTEM) one can write the following general expression in case (i)
(0) 
whereas in case (ii)
(0) 
In Eqs. (THE SURFACE INTEGRATION APPROACH AND ITS APPLICATION FOR THE ELLIPSOIDPLATE SYSTEM) and (THE SURFACE INTEGRATION APPROACH AND ITS APPLICATION FOR THE ELLIPSOIDPLATE SYSTEM) the following notations are introduced
(0) 
In the definition for , and denote two of the Euler angles to which, in general, the particle can be orientated with respect to a righthanded Cartesian coordinate system, attached to the surface of the planar substrate (see Figure 4 in Ref. [13]). Substituting in the above expressions, for the regular (van der Waals) contribution to the NF between an ellipsoid and a plate one obtains in cases (i) and (ii) the following
(0) 
and
(0) 
respectively.
The corresponding expressions in for the CCF arising between the studied objects, within the two considered cases, follow directly from the substitution of Eq. (Influence of the dispersion forces on the thermodynamic Casimir one in nonpolar fluid systems with film geometry) in Eq. (THE SURFACE INTEGRATION APPROACH AND ITS APPLICATION FOR THE ELLIPSOIDPLATE SYSTEM), and have the form
(0) 
and
(0) 
These two equations are strictly valid only when: a) , i.e., points toward the plane; b) , i.e., points toward the plane and c) any value of if , i.e., when points toward the plane. This is because the CCF originates from the confinement of the correlated fluctuations of the fluid at its critical point. When the particle is arbitrarily orientated the projected infinitesimal area lies either between two nonconcentric ellipses, or when the ellipses intersect is their disjoint area. The latter has a part which faces the plate, and hence contributes to the CCF, and such that is orientated away from it (the direction of the component of the unit normal vector coincides with that of the coordinate – see Figure 4 in Ref. [13] and the text therein). Hence Eqs. (THE SURFACE INTEGRATION APPROACH AND ITS APPLICATION FOR THE ELLIPSOIDPLATE SYSTEM) and (THE SURFACE INTEGRATION APPROACH AND ITS APPLICATION FOR THE ELLIPSOIDPLATE SYSTEM) can not be used outside cases a)  c), within which the spectrum of occurred fluctuations is restricted only between the lower half of the ellipsoid’s surface and that of the plate.
Results
Within the meanfield theory the and dependance of the CCF and hence the NF between any pair of objects immersed in near critical fluid medium is given by the temperature and field scaling variables and , respectively [see the text beneath Eq. (Influence of the dispersion forces on the thermodynamic Casimir one in nonpolar fluid systems with film geometry)], with and , i.e., for . In our numerical treatment we take these variables to range in the intervals: and . To study the scaling function of the CCF, the separation between the set of parallel plates is varied from 20 to 120 with step 2, from 23 to 45 with step 1, from 100 to 120 with step 1, from 120 to 220 with step 5 and from 180 to 220 with step 2. In order to demonstrate the effect of the sign change of the considered forces and in view of potential experimental realization of the predicted effects we choose the following values of the coupling parameters: , and . The physical argumentation of the so chosen values is explained in Section IV C of Ref. [7].
The van der Waals and critical Casimir forces between an ellipsoid and a plate in within the SIA
In , considering genuine van der Waals interactions between all constituents comprising the system ellipsoidfluidplate, Eqs. (THE SURFACE INTEGRATION APPROACH AND ITS APPLICATION FOR THE ELLIPSOIDPLATE SYSTEM) and (THE SURFACE INTEGRATION APPROACH AND ITS APPLICATION FOR THE ELLIPSOIDPLATE SYSTEM) become
(0) 
respectively. Here the following notations are introduced: , , with interpreted either as or depending on the considered case, . The quantities normalized by contain an additional prime symbol . Accordingly for the CCF from Eqs. (THE SURFACE INTEGRATION APPROACH AND ITS APPLICATION FOR THE ELLIPSOIDPLATE SYSTEM) and (THE SURFACE INTEGRATION APPROACH AND ITS APPLICATION FOR THE ELLIPSOIDPLATE SYSTEM) one can write
(0) 
where
(0) 
with ; is dimensionless variable; is the largest system for which numerical data are available, the arguments of the scaling function , and is the Heaviside step function with the convention . For the present study . In Eq. (The van der Waals and critical Casimir forces between an ellipsoid and a plate in within the SIA) with respect to case (i) we have
(0) 
while when one considers case (ii)
(0) 
In Eq. (The van der Waals and critical Casimir forces between an ellipsoid and a plate in within the SIA) the critical component of the scaling function calculated within the meanfield theory is normalized by the procedure described inandout in Section V of Ref. [6], so that it contributes properly to the CCF and hence to the net interaction force when . The need of such normalization is explained in details in Ref. [14] (see there Sections IV.A.1 and IV.A.3).
After presenting the mathematical means to calculate the CCF and vdWF, we now pass to the detailed discussion of the results and argumentation of the experimental feasibility of the parameters used in the model calculations utilizing the SIA.
Discussion and concluding remarks
In the present section for physical reasons related to future experimental studies we will assume that the ellipsoidal particle has either ruthenium (Ru) or platinum (Pt) core, while the planar substrate is made either of carbon (C) or silica (SiO) aerogels. The contact surface of any of the interacting objects is considered coated by monolayer of lead (Pb) or thallium nitride (TlN) to ensure the (+,+) boundary conditions. Finally, we complete the so assembled system by adopting xenon (Xe) as the critical fluctuating medium in which the interacting objects are immersed. For detail comment and argumentation on the choice of the presented substances the reader can refer to Ref. [7].
Some results about the behavior of the net force in an ellipsoidplate system are presented on Figure 2. The effect of the particle size is included in the study through the consideration of particles with dimensions: and . For the sake of simplicity we choose . The analysis of the numerical data showed that the critical component of the Casimir force is negative for any value of and , i.e., it corresponds to attraction between the ellipsoid and the plate. On the other hand at , has positive values for , i.e. in the liquid phase of , and is zero for , since it is proportional to the bulk order parameter of the medium. Therefore, the behavior of Casimir interaction coincides with that of its critical component above , but one is able observes a single sign change in the ”liquid” phase of the fluid. The superposition between the and , leads to interaction which is positive for any temperature and shows explicit decrease with the increase of the temperature like (see additionally Eqs. (33) and (34) in Ref. [8]). As a result becomes repulsive outside the critical region, changes sign twice, after which turns attractive only near the critical point. The nature of this observation is dictated by the increased correlations between the fluctuations in the fluid close to . The so described behaviour is well illustrated on subfigured (a) and (c).
Now, at as a function of , is attractive in the ”gas” phase of the fluid, i.e., for and repulsive otherwise with an infinite slope at the critical point, where a sign change occurs. Hence, the CCF changes sign once in the ”liquid” phase () and its minimum increases in magnitude slightly in comparison to . Since does not depend on the chemical potential, when the corresponding vdWF is added to , the resultant interaction ”shift up” and change sign for . Therefore, as presented on subfigured (b) and (d), the net force exhibits behavior similar to that of , but now the sign change occurs closer to the critical point.
With respect to the mutual orientation between the interacting objects, when is fixed, as expected is maximal [subfigured (a’) and (b’)] at any and when , and minimal for any value of and . Additionally, one observes that in the critical region the net interactions for and are very close both quantitatively and qualitatively, when the halfradii or are orientated towards the plate. This observation can be attributed to the increased minimal separation [entering explicitly in Eq. (The van der Waals and critical Casimir forces between an ellipsoid and a plate in within the SIA)] in these two cases in comparison to the one when is towards the plate. Hence irrespective of the initial orientation, if one is able to maintain fixed, the fluctuations in the system are going to drive the particle in a position where plane.
In contrast to the above, for fixed , the behavior of with respect to its angular dependance is in reverse – the force is maximal for any value of and and minimal when [subfigured (c’) and (d’)]. One also observes that in each of the three orientations the minimum of the NF is deeper for than for . This is easily understood having that: the confining surface of is bigger than that of and when plane maximal part of it is close to the plates surface, thus minimazing the separation effect. Therefore, this is the equilibrium orientation which the system will strive to achieve.
Although the magnitude of the NF may seems rather negligible, a comparison with say the weight of a single particle proofs otherwise. For instance, the weight of a platinum ellipsoid in shape particle with dimensions is approximately . In the liquid phase of xenon at and one finds that , and . Therefore the repulsive part of the net force in one such concrete system is indeed capable of levitating a single particle or a system of such. Thus, one can indeed make use of the interplay between the quantum and thermodynamical fluctuations for governing the behavior of objects, say colloidal particles, at small distances. It can also provide a strategy for solving problems with handling, feeding, trapping and fixing of micro parts in nanotechnology resolving the issues related to sticking of the particles on the surface of the mechanical manipulator utilizing, e.g., the reversible dependence on the forces under minute changes of the temperature of the critical medium. One can perform grabbing of particles for small values of , where the force is attractive and release them at a given spacial position after slightly increasing or decreasing of temperature achieving in that way a value of with a repulsive NF.
Acknowledgements
The authors gratefully acknowledges the financial support via Contract No. DN 02/8 of Bulgarian NSF.
References
 Fisher and de Gennes [1978] M. E. Fisher and P. G. de Gennes, C. R. Seances Acad. Sci. Paris Ser. B 287, p. 207 (1978).
 Krech [1994] M. Krech, The Casimir Effect in Critical Systems (World Scientific, Singapore, 1994).
 Brankov, Dantchev, and Tonchev [2000] J. G. Brankov, D. M. Dantchev, and N. S. Tonchev, The Theory of Critical Phenomena in FiniteSize Systems  Scaling and Quantum Effects (World Scientific, Singapore, 2000).
 Casimir [1948] H. Casimir, Proc. K. Ned. Akad. Wet. 51, p. 793 (1948).
 Casimir and Polder [1948] H. Casimir and D. Polder, Phys. Rev. 73, 360–372 (1948).
 Valchev and Dantchev [2015] G. Valchev and D. Dantchev, Phys. Rev. E 92, p. 012119 (2015).
 Valchev and Dantchev [2017] G. Valchev and D. Dantchev, Phys. Rev. E 96, p. 022107 (2017).
 Valchev [2018] G. Valchev, MATEC Web of Conferences 145, p. 01008 (2018).
 Kondrat, Harnau, and Dietrich [2009] S. Kondrat, L. Harnau, and S. Dietrich, J. Chem. Phys. 131, p. 204902 (2009).
 Parsegian [2006] V. A. Parsegian, Van der Waals Forces (Cambridge University Press, New York, 2006).
 Israelachvili [2011] J. N. Israelachvili, Intermolecular and surface forces (Academic, London, 2011).
 Derjaguin [1934] B. Derjaguin, Kolloid Z. 69, p. 155 (1934).
 Dantchev and Valchev [2012] D. Dantchev and G. Valchev, 372, 148 – 163 (2012).
 Dantchev, Schlesener, and Dietrich [2007] D. Dantchev, F. Schlesener, and S. Dietrich, Phys. Rev. E 76, p. 011121 (2007).