Universality of the thermodynamic Casimir effect

Fluctuations of the critical Casimir force

The critical Casimir force (CCF) arises from confining fluctuations in a critical fluid and thus it is a fluctuating quantity itself. While the mean CCF is universal, its (static) variance has previously been found to depend on the microscopic details of the system which effectively set a large-momentum cutoff in the underlying field theory, rendering it potentially large. This raises the question how the properties of the force variance are reflected in experimentally observable quantities, such as the thickness of a wetting film or the position of a suspended colloidal particle. Here, based on a rigorous definition of the instantaneous force, we analyze static and dynamic correlations of the CCF for a conserved fluid in film geometry for various boundary conditions within the Gaussian approximation. We find that the dynamic correlation function of the CCF is independent of the momentum cutoff and decays algebraically in time. Within the Gaussian approximation, the associated exponent depends only on the dynamic universality class but not on the boundary conditions. We furthermore consider a fluid film, the thickness of which can fluctuate under the influence of the time-dependent CCF. The latter gives rise to an effective non-Markovian noise in the equation of motion of the film boundary and induces a distinct contribution to the position variance. Within the approximations used here, at short times, this contribution grows algebraically in time whereas, at long times, it saturates and contributes to the steady-state variance of the film thickness.

Ensemble dependence of critical Casimir forces in films with Dirichlet boundary conditions

In a recent study [Phys. Rev. E 94, 022103 (2016)2470-004510.1103/PhysRevE.94.022103] it has been shown that, for a fluid film subject to critical adsorption, the resulting critical Casimir force (CCF) may significantly depend on the thermodynamic ensemble. Here we extend that study by considering fluid films within the so-called ordinary surface universality class. We focus on mean-field theory, within which the order parameter (OP) profile satisfies Dirichlet boundary conditions and produces a nontrivial CCF in the presence of external bulk fields or, respectively, a nonzero total order parameter within the film. Additionally, we study the influence of fluctuations by means of Monte Carlo simulations of the three-dimensional Ising model. We show that, in the canonical ensemble, i.e., when fixing the so-called total mass within the film, the CCF is repulsive for large absolute values of the total OP, instead of attractive as in the grand canonical ensemble. Based on the Landau-Ginzburg free energy, we furthermore obtain analytic expressions for the order parameter profiles and analyze the relation between the total mass in the film and the external bulk field.

Critical and near-critical phase behavior and interplay between the thermodynamic Casimir and van der Waals forces in a confined nonpolar fluid medium with competing surface and substrate potentials

We study, using general scaling arguments and mean-field type calculations, the behavior of the critical Casimir force and its interplay with the van der Waals force acting between two parallel slabs separated at a distance L from each other, confining some fluctuating fluid medium, say a nonpolar one-component fluid or a binary liquid mixture. The surfaces of the slabs are coated by thin layers exerting strong preference to the liquid phase of the fluid, or one of the components of the mixture, modeled by strong adsorbing local surface potentials ensuring the so-called (+,+) boundary conditions. The slabs, on the other hand, influence the fluid by long-range competing dispersion potentials, which represent irrelevant interactions in renormalization-group sense. Under such conditions, one usually expects attractive Casimir force governed by universal scaling function, pertinent to the extraordinary surface universality class of Ising type systems, to which the dispersion potentials provide only corrections to scaling. We demonstrate, however, that below a given threshold thickness of the system L(crit) for a suitable set of slabs-fluid and fluid-fluid coupling parameters the competition between the effects due to the coatings and the slabs can result in sign change of the Casimir force acting between the surfaces confining the fluid when one changes the temperature T, the chemical potential of the fluid μ, or L. The last implies that by choosing specific materials for the slabs, coatings, and the fluid for L≲L(crit) one can realize repulsive Casimir force with nonuniversal behavior which, upon increasing L, gradually turns into an attractive one described by a universal scaling function, depending only on the relevant scaling fields related to the temperature and the excess chemical potential, for L≫L(crit). We present arguments and relevant data for specific substances in support of the experimental feasibility of the predicted behavior of the force. It can be of interest, e.g., for designing nanodevices and for governing behavior of objects, say colloidal particles, at small distances. We formulate the corresponding criterion for determination of L(crit). The universality is regained for L≫L(crit). We also show that for systems with L≲L(crit), the capillary condensation phase diagram suffers modifications which one does not observe in systems with purely short-ranged interactions.

Critical Casimir interactions around the consolute point of a binary solvent

Spatial confinement of a near-critical medium changes its fluctuation spectrum and modifies the corresponding order parameter distribution, resulting in effective, so-called critical Casimir forces (CCFs) acting on the confining surfaces. These forces are attractive for like boundary conditions of the order parameter at the opposing surfaces of the confinement. For colloidal particles dissolved in a binary liquid mixture acting as a solvent close to its critical point of demixing, one thus expects the emergence of phase segregation into equilibrium colloidal liquid and gas phases. We analyze how such phenomena occur asymmetrically in the whole thermodynamic neighborhood of the consolute point of the binary solvent. By applying field-theoretical methods within mean-field approximation and the semi-empirical de Gennes-Fisher functional, we study the CCFs acting between planar parallel walls as well as between two spherical colloids and their dependence on temperature and on the composition of the near-critical binary mixture. We find that for compositions slightly poor in the molecules preferentially adsorbed at the surfaces, the CCFs are significantly stronger than at the critical composition, thus leading to pronounced colloidal segregation. The segregation phase diagram of the colloid solution following from the calculated effective pair potential between the colloids agrees surprisingly well with experiments and simulations.

Casimir force in the rotor model with twisted boundary conditions

We investigate the three-dimensional lattice XY model with nearest neighbor interaction. The vector order parameter of this system lies on the vertices of a cubic lattice, which is embedded in a system with a film geometry. The orientations of the vectors are fixed at the two opposite sides of the film. The angle between the vectors at the two boundaries is α where 0≤α≤π. We make use of the mean field approximation to study the mean length and orientation of the vector order parameter throughout the film--and the Casimir force it generates--as a function of the temperature T, the angle α, and the thickness L of the system. Among the results of that calculation are a Casimir force that depends in a continuous way on both the parameter α and the temperature and that can be attractive or repulsive. In particular, by varying α and/or T one controls both the sign and the magnitude of the Casimir force in a reversible way. Furthermore, for the case α=π, we discover an additional phase transition occurring only in the finite system associated with the variation of the orientations of the vectors.

Aspect-ratio dependence of thermodynamic Casimir forces

We consider the three-dimensional Ising model in a L(⊥)×L(∥)×L(∥) cuboid geometry with a finite aspect ratio ρ=L(⊥)/L(∥) and periodic boundary conditions along all directions. For this model the finite-size scaling functions of the excess free energy and thermodynamic Casimir force are evaluated numerically by means of Monte Carlo simulations. The Monte Carlo results compare well with recent field theoretical results for the Ising universality class at temperatures above and slightly below the bulk critical temperature T(c). Furthermore, the excess free energy and Casimir force scaling functions of the two-dimensional Ising model are calculated exactly for arbitrary ρ and compared to the three-dimensional case. We give a general argument that the Casimir force vanishes at the critical point for ρ=1 and becomes repulsive in periodic systems for ρ>1.

Finite-size effects in film geometry with nonperiodic boundary conditions: Gaussian model and renormalization-group theory at fixed dimension

Finite-size effects are investigated in the Gaussian model with isotropic and anisotropic short-range interactions in film geometry with nonperiodic boundary conditions (bc) above, at, and below the bulk critical temperature Tc. We have obtained exact results for the free energy and the Casimir force for antiperiodic, Neumann, Dirichlet, and Neumann-Dirichlet mixed bc in 1<d<4 dimensions. For the Casimir force, finite-size scaling is found to be valid for all bc. For the free energy, finite-size scaling is valid in 1<d<3 and 3<d<4 dimensions for antiperiodic, Neumann, and Dirichlet bc, but logarithmic deviations from finite-size scaling exist in d=3 dimensions for Neumann and Dirichlet bc. This is explained in terms of the borderline dimension d*=3 , where the critical exponent 1-α-ν=(d-3)∕2 of the Gaussian surface energy density vanishes. For Neumann-Dirichlet bc, finite-size scaling is strongly violated above Tc for 1<d<4 because of a cancelation of the leading scaling terms. For antiperiodic, Dirichlet, and Neumann-Dirichlet bc, a finite film critical temperature Tc,film(L)<Tc exists at finite film thickness L . Our results include an exact description of the dimensional crossover between the d -dimensional finite-size critical behavior near bulk Tc and the (d-1) -dimensional critical behavior near Tc,film(L). This dimensional crossover is illustrated for the critical behavior of the specific heat. Particular attention is paid to an appropriate representation of the free energy in the region Tc,film(L)≤T≤Tc. For 2<d<4 , the Gaussian results are renormalized and reformulated as one-loop contributions of the φ4 field theory at fixed dimension d and are then compared with the ε=4-d expansion results at ε=1 as well as with d=3 Monte Carlo data. For d=2 , the Gaussian results for the Casimir force scaling function are compared with those for the Ising model with periodic, antiperiodic, and free bc; unexpected exact relations are found between the Gaussian and Ising scaling functions. For both the d -dimensional Gaussian model and the two-dimensional Ising model it is shown that anisotropic couplings imply nonuniversal scaling functions of the Casimir force that depend explicitly on microscopic couplings. Our Gaussian results provide the basis for the investigation of finite-size effects of the mean spherical model in film geometry with nonperiodic bc above, at, and below the bulk critical temperature.

Pair correlation function decay in models of simple fluids that contain dispersion interactions

We investigate the intermediate-and longest-range decay of the total pair correlation function h(r) in model fluids where the inter-particle potential decays as -r(-6), as is appropriate to real fluids in which dispersion forces govern the attraction between particles. It is well-known that such interactions give rise to a term in q(3) in the expansion of [Formula: see text], the Fourier transform of the direct correlation function. Here we show that the presence of the r(-6) tail changes significantly the analytic structure of [Formula: see text] from that found in models where the inter-particle potential is short ranged. In particular the pure imaginary pole at q = iα(0), which generates monotonic-exponential decay of rh(r) in the short-ranged case, is replaced by a complex (pseudo-exponential) pole at q = iα(0)+α(1) whose real part α(1) is negative and generally very small in magnitude. Near the critical point α(1)∼-α(0)(2) and we show how classical Ornstein-Zernike behaviour of the pair correlation function is recovered on approaching the mean-field critical point. Explicit calculations, based on the random phase approximation, enable us to demonstrate the accuracy of asymptotic formulae for h(r) in all regions of the phase diagram and to determine a pseudo-Fisher-Widom (pFW) line. On the high density side of this line, intermediate-range decay of rh(r) is exponentially damped-oscillatory and the ultimate long-range decay is power-law, proportional to r(-6), whereas on the low density side this damped-oscillatory decay is sub-dominant to both monotonic-exponential and power-law decay. Earlier analyses did not identify the pseudo-exponential pole and therefore the existence of the pFW line. Our results enable us to write down the generic wetting potential for a 'real' fluid exhibiting both short-ranged and dispersion interactions. The monotonic-exponential decay of correlations associated with the pseudo-exponential pole introduces additional terms into the wetting potential that are important in determining the existence and order of wetting transitions.

Casimir force in O(n) systems with a diffuse interface

We study the behavior of the Casimir force in O(n) systems with a diffuse interface and slab geometry infinity;{d-1}xL , where 2<d<4 is the dimensionality of the system. We consider a system with nearest-neighbor anisotropic interaction constants J_{ parallel} parallel to the film and J_{ perpendicular} across it. We argue that in such an anisotropic system the Casimir force, the free energy, and the helicity modulus will differ from those of the corresponding isotropic system, even at the bulk critical temperature, despite that these systems both belong to the same universality class. We suggest a relation between the scaling functions pertinent to the both systems. Explicit exact analytical results for the scaling functions, as a function of the temperature T , of the free energy density, Casimir force, and the helicity modulus are derived for the n-->infinity limit of O(n) models with antiperiodic boundary conditions applied along the finite dimension L of the film. We observe that the Casimir amplitude Delta_{Casimir}(dmid R:J_{ perpendicular},J_{ parallel}) of the anisotropic d -dimensional system is related to that of the isotropic system Delta_{Casimir}(d) via Delta_{Casimir}(dmid R:J_{ perpendicular},J_{ parallel})=(J_{ perpendicular}J_{ parallel});{(d-1)2}Delta_{Casimir}(d) . For d=3 we derive the exact Casimir amplitude Delta_{Casimir}(3,mid R:J_{ perpendicular},J_{ parallel})=[Cl_{2}(pi3)3-zeta(3)(6pi)](J_{ perpendicular}J_{ parallel}) , as well as the exact scaling functions of the Casimir force and of the helicity modulus Upsilon(T,L) . We obtain that beta_{c}Upsilon(T_{c},L)=(2pi;{2})[Cl_{2}(pi3)3+7zeta(3)(30pi)](J_{ perpendicular}J_{ parallel})L;{-1} , where T_{c} is the critical temperature of the bulk system. We find that the contributions in the excess free energy due to the existence of a diffuse interface result in a repulsive Casimir force in the whole temperature region.

Diversity of critical behavior within a universality class

We study spatial anisotropy effects on the bulk and finite-size critical behavior of the O(n) symmetric anisotropic phi;{4} lattice model with periodic boundary conditions in a d -dimensional hypercubic geometry above, at, and below Tc. The absence of two-scale factor universality is discussed for the bulk order-parameter correlation function, the bulk scattering intensity, and for several universal bulk amplitude relations. The anisotropy parameters are observable by scattering experiments at Tc. For the confined system, renormalization-group theory within the minimal subtraction scheme at fixed dimension d for 2<d<4 is employed. In contrast to the epsilon=4-d expansion, the fixed- d finite-size approach keeps the exponential form of the order-parameter distribution function unexpanded. For the case of cubic symmetry and for n=1 , our perturbation approach yields excellent agreement with the Monte Carlo (MC) data for the finite-size amplitude of the free energy of the three-dimensional Ising model at Tc by Mon [Phys. Rev. Lett. 54, 2671 (1985)]. The epsilon expansion result is in less good agreement. Below Tc, a minimum of the scaling function of the excess free energy is found. We predict a measurable dependence of this minimum on the anisotropy parameters. The relative anisotropy effect on the free energy is predicted to be significantly larger than that on the Binder cumulant. Our theory agrees quantitatively with the nonmonotonic dependence of the Binder cumulant on the ferromagnetic next-nearest-neighbor (NNN) coupling of the two-dimensional Ising model found by MC simulations of Selke and Shchur [J. Phys. A 38, L739 (2005)]. Our theory also predicts a nonmonotonic dependence for small values of the antiferromagnetic NNN coupling and the existence of a Lifshitz point at a larger value of this coupling. The nonuniversal anisotropy effects in the finite-size scaling regime are predicted to satisfy a kind of restricted universality. The tails of the large- L behavior at T++Tc violate both finite-size scaling and universality even for isotropic systems as they depend on the bare four-point coupling of the phi4 theory, on the cutoff procedure, and on subleading long-range interactions.

Interplay of critical Casimir and dispersion forces

Using general scaling arguments combined with mean-field theory we investigate the critical (T approximately Tc) and off-critical (T not equal Tc) behavior of the Casimir forces in fluid films of thickness L governed by dispersion forces and exposed to long-ranged substrate potentials which are taken to be equal on both sides of the film. We study the resulting effective force acting on the confining substrates as a function of T and of the chemical potential mu. We find that the total force is attractive both below and above Tc. If, however, the direct substrate-substrate contribution is subtracted, the force is repulsive everywhere except near the bulk critical point (Tc, mu(c)), where critical density fluctuations arise, or except at low temperatures and (L/a)(beta(Delta)(mu))=O(1), with Delta(mu)=mu-mu(c)<0 and a the characteristic distance between the molecules of the fluid, i.e., in the capillary condensation regime. While near the critical point the maximal amplitude of the attractive force if of order of L(-d) in the capillary condensation regime the force is much stronger with maximal amplitude decaying as L(-1). In the latter regime we observe that the long-ranged tails of the fluid-fluid and the substrate-fluid interactions further increase that amplitude in comparison with systems with short-range interactions only. Although in the critical region the system under consideration asymptotically belongs to the Ising universality class with short-ranged forces, we find deviations from the standard finite-size scaling for xi(ln)(xi/xi0(+/-)) >>L even for xi, L>>xi0(+/-), where xi[t=(T-Tc)/Tc-->+/-0,Delta(mu)=0]=xi0(+/-)/t/-nu, is the bulk correlation length. In this regime the dominant finite-size contributions to the free energy and to the force stem from the long-ranged algebraically decaying tails of the interactions; they are not exponentially small in L, as it is the case there in systems governed by purely short-ranged interactions, but exhibit a power law decay in L. Essential deviations from the standard finite-size scaling behavior are observed also within the finite-size critical region L/xi=O(1) for films with thicknesses L less than or approximately equal Lcrit, where Lcrit=xi0(+/-)(16/s/)nu/beta, with nu and beta as the standard bulk critical exponents and with s=O(1) as the dimensionless parameter that characterizes the relative strength of the long-ranged tail of the substrate fluid over the fluid-fluid interaction. We present the modified finite-size scaling pertinent for such a case and analyze in detail the finite-size behavior in this region. The standard finite-size scaling behavior is recovered only for L>>Lcrit.

Finite-size effects on the behavior of the susceptibility in van der Waals films bounded by strongly absorbing substrates

We study critical point finite-size effects in the case of the susceptibility of a film in which interactions are characterized by a van der Waals-type power law tail. The geometry is appropriate to a slablike system with two bounding surfaces. Boundary conditions are consistent with surfaces that both prefer the same phase in the low temperature, or broken symmetry, state. We take into account both interactions within the system and interactions between the constituents of the system and the material surrounding it. Specific predictions are made with respect to the behavior of 3He and 4He films in the vicinity of their respective liquid-vapor critical points.

Excess free energy and Casimir forces in systems with long-range interactions of van der Waals type: general considerations and exact spherical-model results

We consider systems confined to a d-dimensional slab of macroscopic lateral extension and finite thickness L that undergo a continuous bulk phase transition in the limit L --> infinity and are describable by an O(n) symmetrical Hamiltonian. Periodic boundary conditions are applied across the slab. We study the effects of long-range pair interactions whose potential decays as bx-(d+sigma) as x --> infinity, with 2<sigma<4 and 2<d+sigma< or =6, on the Casimir effect at and near the bulk critical temperature Tc,infinity, for 2<d<4. These interactions decay sufficiently fast to leave bulk critical exponents and other universal bulk quantities unchanged--i.e., they are irrelevant in the renormalization group (RG) sense. Yet they entail important modifications of the standard scaling behavior of the excess free energy and the Casimir force Fc. We generalize the phenomenological scaling Ansätze for these quantities by incorporating these long-range interactions. For the scaled reduced Casimir force per unit cross-sectional area, we obtain the form LdFc/kBt approximately Xi0(L/xi infinity) + g omegaL -omega Xi omega (L/Xi infinity) + g sigma L -omega sigma Xi sigma (L/Xi infinity). Here Xi0, Xi omega, and Xi sigma are universal scaling functions; g omega and g sigma are scaling fields associated with the leading corrections to scaling and those of the long-range interaction, respectively; omega and omega sigma = sigma + eta - 2 are the associated correction-to-scaling exponents, where eta denotes the standard bulk correlation exponent of the system without long-range interactions; xi infinity is the (second-moment) bulk correlation length (which itself involves corrections to scaling). The contribution proportional variant g sigma decays for T not = Tc,infinity algebraically in L rather than exponentially, and hence becomes dominant in an appropriate regime of temperatures and L. We derive exact results for spherical and Gaussian models which confirm these findings. In the case d + sigma = 6, which includes that of nonretarded van der Waals interactions in d = 3 dimensions, the power laws of the corrections to scaling proportional to b of the spherical model are found to get modified by logarithms. Using general RG ideas, we show that these logarithmic singularities originate from the degeneracy omega = omega sigma = 4 - d that occurs for the spherical model when d + sigma = 6, in conjunction with the b dependence of g omega.

Critical Casimir forces for O(n) systems with long-range interaction in the spherical limit

We present exact results on the behavior of the thermodynamic Casimir force and the excess free energy in the framework of the d -dimensional spherical model with a power law long-ranged interaction decaying at large distances r as r(-d-sigma) , where sigma<d<2sigma and 0<sigma< or =2 . For a film geometry and under periodic boundary conditions we consider the behavior of these quantities near the bulk critical temperature T(c) , as well as for T> T(c) and T< T(c) . The universal finite-size scaling function governing the behavior of the force in the critical region is derived and its asymptotics are investigated. While in the critical and subcritical region the force is of the order of L(-d) , for T> T(c) it decays as L(-d-sigma) , where L is the thickness of the film. We consider both the case of a finite system that has no phase transition of its own, when d-1<sigma , as well as the case with d-1>sigma , when one observes a dimensional crossover from d to a d-1 dimensional critical behavior. The behavior of the force along the phase coexistence line for a magnetic field H=0 and T< T(c) is also derived. We have proven analytically that the excess free energy is always negative and monotonically increasing function of T and H . For the Casimir force we have demonstrated that for any sigma > or =1 it is everywhere negative, i.e., an attraction between the surfaces bounding the system is to be observed. At T= T(c) the force is an increasing function of T for sigma>1 and a decreasing one for sigma<1 . For any d and sigma the minimum of the force at T= T(c) is always achieved at some H unequal to 0 .

Solvation force for long-ranged wall–fluid potentials

The solvation force of a simple fluid confined between identical planar walls is studied in two model systems with short ranged fluid-fluid interactions and long-ranged wall-fluid potentials decaying as -Az(-p),z--> infinity, for various values of p. Results for the Ising spins system are obtained in two dimensions at vanishing bulk magnetic field h=0 by means of the density-matrix renormalization-group method; results for the truncated Lennard-Jones (LJ) fluid are obtained within the nonlocal density functional theory. At low temperatures the solvation force f(solv) for the Ising film is repulsive and decays for large wall separations L in the same fashion as the boundary field f(solv) approximately L(-p), whereas for temperatures larger than the bulk critical temperature f(solv) is attractive and the asymptotic decay is f(solv) approximately L(-(p+1)). For the LJ fluid system f(solv) is always repulsive away from the critical region and decays for large L with the the same power law as the wall-fluid potential. We discuss the influence of the critical Casimir effect and of capillary condensation on the behavior of the solvation force.