Specific heat of critical films, the Casimir force, and wetting films near critical end points

Critical Casimir forces in soft matter

We review recent advances in the theoretical, numerical, and experimental studies of critical Casimir forces in soft matter, with particular emphasis on their relevance for the structures of colloidal suspensions and on their dynamics. Distinct from other interactions which act in soft matter, such as electrostatic and van der Waals forces, critical Casimir forces are effective interactions characterised by the possibility to control reversibly their strength via minute temperature changes, while their attractive or repulsive character is conveniently determined via surface treatments or by structuring the involved surfaces. These features make critical Casimir forces excellent candidates for controlling the equilibrium and dynamical properties of individual colloids or colloidal dispersions as well as for possible applications in micro-mechanical systems. In the past 25 years a number of theoretical and experimental studies have been devoted to investigating these forces primarily under thermal equilibrium conditions, while their dynamical and non-equilibrium behaviour is a largely unexplored subject open for future investigations.

Multiparameter universality and intrinsic diversity of critical phenomena in weakly anisotropic systems

Recently a unified hypothesis of multiparameter universality for the critical behavior of bulk and confined anisotropic systems has been formulated [V. Dohm, Phys. Rev. E 97, 062128 (2018)2470-004510.1103/PhysRevE.97.062128]. We prove the validity of this hypothesis in d≥2 dimensions on the basis of the principle of two-scale-factor universality for isotropic systems at vanishing external field. We introduce an angular-dependent correlation vector and a generalized shear transformation that transforms weakly anisotropic systems to isotropic systems. As examples we consider the O(n)-symmetric φ^{4} model, Gaussian model, and n-vector model. By means of the inverse of the shear transformation we determine the general structure of the bulk order-parameter correlation function, of the singular bulk part of the critical free energy, and of critical bulk amplitude relations of anisotropic systems at and away from T_{c}. It is shown that weakly anisotropic systems exhibit a high degree of intrinsic diversity due to d(d+1)/2-1 independent parameters that cannot be determined by thermodynamic measurements. Exact results are derived for the d=2 Ising universality class and for the spherical and Gaussian universality classes in d≥2 dimensions. For the d=3 Ising universality class we identify the universal scaling function of the isotropic bulk correlation function from the nonuniversal result of the functional renormalization group. A proof is presented for the validity of multiparameter universality of the exact critical free energy and critical Casimir amplitude in a finite rectangular geometry of weakly anisotropic systems with periodic boundary conditions in the Ising universality class. This confirms the validity of recent predictions of self-similar structures of finite-size effects in the (d=2,n=1) universality class at T=T_{c} derived from conformal field theory [V. Dohm and S. Wessel, Phys. Rev. Lett. 126, 060601 (2021)PRLTAO0031-900710.1103/PhysRevLett.126.060601]. This also substantiates the previous notion of an effective shear transformation for anisotropic two-dimensional Ising models. Our theory paves the way for a quantitative theory of nonuniversal critical Casimir forces in anisotropic superconductors for which experiments have been proposed by G. A. Williams [Phys. Rev. Lett. 92, 197003 (2004)PRLTAO0031-900710.1103/PhysRevLett.92.197003].

Effective pair interaction of patchy particles in critical fluids

We study the critical Casimir interaction between two spherical colloids immersed in a binary liquid mixture close to its critical demixing point. The surface of each colloid prefers one species of the mixture with the exception of a circular patch of arbitrary size, where the other species is preferred. For such objects, we calculate, within the Derjaguin approximation, the scaling function describing the critical Casimir potential, and we use it to derive the scaling functions for all components of the forces and torques acting on both colloids. The results are compared with available experimental data. Moreover, the general relation between the scaling function for the potential and the scaling functions for the force and the torque is derived.

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.

Crossover from low-temperature to high-temperature fluctuations: Universal and nonuniversal Casimir forces of isotropic and anisotropic systems

We study the crossover from low-temperature to high-temperature fluctuations including Goldstone-dominated and critical fluctuations in confined isotropic and weakly anisotropic O(n)-symmetric systems on the basis of a finite-size renormalization-group approach at fixed dimension d introduced previously [V. Dohm, Phys. Rev. Lett. 110, 107207 (2013)PRLTAO0031-900710.1103/PhysRevLett.110.107207]. Our theory is formulated within the φ^{4} lattice model in a d-dimensional block geometry with periodic boundary conditions. We calculate the finite-size scaling functions F^{ex} and X of the excess free-energy density and the thermodynamic Casimir force, respectively, for 1≤n≤∞, 2<d<4. Exact results are derived for n→∞. Applications are given for L_{∥}^{d-1}×L slab geometry with an aspect ratio ρ=L/L_{∥}>0 and for film geometry (ρ=0). Good overall agreement is found with Monte Carlo (MC) data for isotropic spin models with n=1,2,3. For ρ=0, the low-temperature limits of F^{ex} and X vanish for n=1, whereas they are finite for n≥2. For ρ>0 and n=1, we find a finite low-temperature limit of F^{ex}, which deviates from that of the Ising model. We attribute this deviation to the nonuniversal difference between the φ^{4} model with continuous variables and the Ising model with discrete variables. For n≥2 and ρ>0, a logarithmic divergence of F^{ex} in the low-temperature limit is predicted, in excellent agreement with MC data. For 2≤n≤∞ and ρ<ρ_{0}=0.8567 the Goldstone modes generate a negative low-temperature Casimir force that vanishes for ρ=ρ_{0} and becomes positive for ρ>ρ_{0}. For anisotropic systems a unified hypothesis of multiparameter universality is introduced for both bulk and confined systems. The dependence of their scaling functions on d(d+1)/2-1 microscopic anisotropy parameters implies a substantial reduction of the predictive power of the theory for anisotropic systems as compared to isotropic systems. An exact representation is derived for the nonuniversal large-distance behavior of the bulk correlation function of anisotropic systems and quantitative predictions are made. The validity of multiparameter universality is proven analytically for the d=2,n=1 universality class. A nonuniversal anisotropy-dependent minimum of the Casimir force scaling function X is found. Both the sign and magnitude of X and the shift of the film critical temperature are affected by the lattice anisotropy.

Tricritical Casimir forces and order parameter profiles in wetting films of ^{3}He-^{4}He mixtures

Tricritical Casimir forces in ^{3}He-^{4}He wetting films are studied, within mean field theory, in terms of a suitable lattice gas model for binary liquid mixtures with short-ranged surface fields. The proposed model takes into account the continuous rotational symmetry O(2) of the superfluid degrees of freedom associated with ^{4}He and it allows, inter alia, for the occurrence of a vapor phase. As a result, the model facilitates the formation of wetting films, which provides a strengthened theoretical framework to describe available experimental data for tricritical Casimir forces acting in ^{3}He-^{4}He wetting films.

Critical Casimir interactions between Janus particles

Recently there has been strong experimental and theoretical interest in studying the self-assembly and the phase behavior of patchy and Janus particles, which form colloidal suspensions. Although in this quest a variety of effective interactions have been proposed and used in order to achieve a directed assembly, the critical Casimir effect stands out as being particularly suitable in this respect because it provides both attractive and repulsive interactions as well as the potential of a sensitive temperature control of their strength. Specifically, we have calculated the critical Casimir force between a single Janus particle and a laterally homogeneous substrate as well as a substrate with a chemical step. We have used the Derjaguin approximation and compared it with results from full mean field theory. A modification of the Derjaguin approximation turns out to be generally reliable. Based on this approach we have derived the effective force and the effective potential between two Janus cylinders as well as between two Janus spheres.

Hydrodynamics in bridging and aggregation of two colloidal particles in a near-critical binary mixture

We investigate bridging and aggregation of two colloidal particles in a near-critical binary mixture when the fluid far from the particles is outside the coexistence (CX) curve and is rich in the component disfavored by the colloid surfaces. In such situations, the adsorption-induced interaction is enhanced, leading to bridging and aggregation of the particles. We realize bridging firstly by changing the temperature with a fixed interparticle separation and secondly by letting the two particles aggregate. The interparticle attractive force dramatically increases upon bridging. The dynamics is governed by hydrodynamic flow around the colloid surfaces. In aggregation, the adsorption layers move with the particles and squeezing occurs at narrow separation. These results suggest relevance of bridging in the reversible colloid aggregation observed so far. We use the local functional theory [J. Chem. Phys., 2012, 136, 114704] to take into account the renormalization effect and the simulation method [Phys. Rev. Lett., 2000, 85, 1338] to calculate the hydrodynamic flow around the colloidal particles.

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.

Line contribution to the critical Casimir force between a homogeneous and a chemically stepped surface

Recent experimental realizations of the critical Casimir effect have been implemented by monitoring colloidal particles immersed in a binary liquid mixture near demixing and exposed to a chemically structured substrate. In particular, critical Casimir forces have been measured for surfaces consisting of stripes with periodically alternating adsorption preferences, forming chemical steps between them. Motivated by these experiments, we analyze the contribution of such chemical steps to the critical Casimir force for the film geometry and within the Ising universality class. By means of Monte Carlo simulations, mean-field theory and finite-size scaling analysis we determine the universal scaling function associated with the contribution to the critical Casimir force due to individual, isolated chemical steps facing a surface with homogeneous adsorption preference or with Dirichlet boundary condition. In line with previous findings, these results allow one to compute the critical Casimir force for the film geometry and in the presence of arbitrarily shaped, but wide stripes. In this latter limit the force decomposes into a sum of the contributions due to the two homogeneous parts of the surface and due to the chemical steps between the stripes. We assess this decomposition by comparing the resulting sum with actual simulation data for the critical Casimir force in the presence of a chemically striped substrate.

Critical Casimir forces between planar and crenellated surfaces

We study critical Casimir forces between planar walls and geometrically structured substrates within mean-field theory. As substrate structures, crenellated surfaces consisting of periodic arrays of rectangular crenels and merlons are considered. Within the widely used proximity force approximation, both the top surfaces of the merlons and the bottom surfaces of the crenels contribute to the critical Casimir force. However, for such systems the full, numerically determined critical Casimir forces deviate significantly from the pairwise addition formalism underlying the proximity force approximation. A first-order correction to the proximity force approximation is presented in terms of a step contribution arising from the critical Casimir interaction between a planar substrate and the right-angled steps of the merlons consisting of their upper and lower edges as well as their sidewalls.

Critical Casimir forces in the presence of random surface fields

We study critical Casimir forces (CCFs) fC for films of thickness L which in the three-dimensional bulk belong to the Ising universality class and which are exposed to random surface fields (RSFs) on both surfaces. We consider the case in which, in the absence of RSFs, the surfaces of the film belong to the surface universality class of the so-called ordinary transition. We carry out a finite-size scaling analysis and show that for weak disorder, CCFs still exhibit scaling, acquiring a random field scaling variable w that is zero for pure systems. We confirm these analytic predictions by Monte Carlo (MC) simulations. Moreover, our MC data show that fC varies as fC(w→0)-fC(w=0)∼w2. Asymptotically, for large L, w scales as w∼L-0.26→0, indicating that this type of disorder is an irrelevant perturbation of the ordinary surface universality class. However, for thin films such that w≃1, we find that the presence of RSFs with vanishing mean value increases significantly the strength of CCFs, as compared to systems without them, and it shifts the extremum of the scaling function of fC toward lower temperatures. But fC remains attractive.

Thermodynamic Casimir effect in films: the exchange cluster algorithm

We study the thermodynamic Casimir force for films with various types of boundary conditions and the bulk universality class of the three-dimensional Ising model. To this end, we perform Monte Carlo simulations of the improved Blume-Capel model on the simple cubic lattice. In particular, we employ the exchange or geometric cluster cluster algorithm [Heringa and Blöte, Phys. Rev. E 57, 4976 (1998)]. In a previous work, we demonstrated that this algorithm allows us to compute the thermodynamic Casimir force for the plate-sphere geometry efficiently. It turns out that also for the film geometry a substantial reduction of the statistical error can achieved. Concerning physics, we focus on (O,O) boundary conditions, where O denotes the ordinary surface transition. These are implemented by free boundary conditions on both sides of the film. Films with such boundary conditions undergo a phase transition in the universality class of the two-dimensional Ising model. We determine the inverse transition temperature for a large range of thicknesses L(0) of the film and study the scaling of this temperature with L(0). In the neighborhood of the transition, the thermodynamic Casimir force is affected by finite size effects, where finite size refers to a finite transversal extension L of the film. We demonstrate that these finite size effects can be computed by using the universal finite size scaling function of the free energy of the two-dimensional Ising model.

Phase diagram of fluid phases in (3)He-(4)He mixtures

Fluid parts of the phase diagram of (3)He-(4)He mixtures are obtained from a mean-field analysis of a suitable lattice gas model for binary liquid mixtures. The proposed model takes into account the continuous rotational symmetry O(2) of the superfluid degrees of freedom associated with (4)He and includes the occurrence of vacancies. This latter degree of freedom allows the model to exhibit a vapor phase and hence can provide the theoretical framework to describe the experimental conditions for measurements of tricritical Casimir forces in (3)He-(4)He wetting films.

Pronounced minimum of the thermodynamic Casimir forces of O(n) symmetric film systems: analytic theory

Thermodynamic Casimir forces of film systems in the O(n) universality classes with Dirichlet boundary conditions are studied below bulk criticality. Substantial progress is achieved in resolving the long-standing problem of describing analytically the pronounced minimum of the scaling function observed experimentally in ^{4}He films (n=2) by Garcia and Chan [Phys. Rev. Lett. 83, 1187 (1999)] and in Monte Carlo simulations for the three-dimensional Ising model (n=1) by O. Vasilyev et al. [Europhys. Lett. 80, 60009 (2007)]. Our finite-size renormalization-group approach describes the film systems as the limit of finite-slab systems with vanishing aspect ratio. This yields excellent agreement with the depth and the position of the minimum for n=1 and semiquantitative agreement with the minimum for n=2. Our theory also predicts a pronounced minimum for the n=3 Heisenberg universality class.

Alignment of cylindrical colloids near chemically patterned substrates induced by critical Casimir torques

Recent experiments have demonstrated a fluctuation-induced lateral trapping of spherical colloidal particles immersed in a binary liquid mixture near its critical demixing point and exposed to chemically patterned substrates. Inspired by these experiments, we study this kind of effective interaction, known as the critical Casimir effect, for elongated colloids of cylindrical shape. This adds orientational degrees of freedom. When the colloidal particles are close to a chemically structured substrate, a critical Casimir torque acting on the colloids emerges. We calculate this torque on the basis of the Derjaguin approximation. The range of validity of the latter is assessed via mean-field theory. This assessment shows that the Derjaguin approximation is reliable in experimentally relevant regimes, so that we extend it to Janus particles endowed with opposing adsorption preferences. Our analysis indicates that critical Casimir interactions are capable of achieving well-defined, reversible alignments both of chemically homogeneous and of Janus cylinders.

Casimir force in the O(n→∞) model with free boundary conditions

We present results for the temperature behavior of the Casimir force for a system with a film geometry with thickness L subject to free boundary conditions and described by the n→∞ limit of the O(n) model. These results extend over all temperatures, including the critical regime near the bulk critical temperature Tc, where the critical fluctuations determine the behavior of the force, and temperatures well below it, where its behavior is dictated by the Goldstone mode contributions. The temperature behavior when the absolute temperature, T, is a finite distance below Tc, up to a logarithmic-in-L proximity of the bulk critical temperature, is obtained both analytically and numerically; the critical behavior follows from numerics. The results resemble-but do not duplicate-the experimental curve behavior for the force obtained for He4 films.

Critical Casimir forces between homogeneous and chemically striped surfaces

Recent experiments have measured the critical Casimir force acting on a colloid immersed in a binary liquid mixture near its continuous demixing phase transition and exposed to a chemically structured substrate. Motivated by these experiments, we study the critical behavior of a system, which belongs to the Ising universality class, for the film geometry with one planar wall chemically striped, such that there is a laterally alternating adsorption preference for the two species of the binary liquid mixture, which is implemented by surface fields. For the opposite wall we employ alternatively a homogeneous adsorption preference or homogeneous Dirichlet boundary conditions, which within a lattice model are realized by open boundary conditions. By means of mean-field theory, Monte Carlo simulations, and finite-size scaling analysis we determine the critical Casimir force acting on the two parallel walls and its corresponding universal scaling function. We show that in the limit of stripe widths small compared with the film thickness, on the striped surface the system effectively realizes Dirichlet boundary conditions, which generically do not hold for actual fluids. Moreover, the critical Casimir force is found to be attractive or repulsive, depending on the width of the stripes of the chemically patterned surface and on the boundary condition applied to the opposing surface.

Critical Casimir effect for colloids close to chemically patterned substrates

Colloids immersed in a critical or near-critical binary liquid mixture and close to a chemically patterned substrate are subject to normal and lateral critical Casimir forces of dominating strength. For a single colloid, we calculate these attractive or repulsive forces and the corresponding critical Casimir potentials within mean-field theory. Within this approach we also discuss the quality of the Derjaguin approximation and apply it to Monte Carlo simulation data available for the system under study. We find that the range of validity of the Derjaguin approximation is rather large and that it fails only for surface structures which are very small compared to the geometric mean of the size of the colloid and its distance from the substrate. For certain chemical structures of the substrate, the critical Casimir force acting on the colloid can change sign as a function of the distance between the particle and the substrate; this provides a mechanism for stable levitation at a certain distance which can be strongly tuned by temperature, i.e., with a sensitivity of more than 200 nm/K.

Critical Casimir forces for Ising films with variable boundary fields

Monte Carlo simulations based on an integration scheme for free energy differences is used to compute critical Casimir forces for three-dimensional Ising films with various boundary fields. We study the scaling behavior of the critical Casimir force, including the scaling variable related to the boundary fields. Finite size corrections to scaling are taken into account. We pay special attention to that range of surface field strengths within which the force changes from repulsive to attractive upon increasing the temperature. Our data are compared with other results available in the literature.

Critical free energy and Casimir forces in rectangular geometries

We study the critical behavior of the free energy and the thermodynamic Casimir force in a L(∥)(d-1) × L block geometry in 2<d<4 dimensions with aspect ratio ρ=L/L(∥) on the basis of the O(n) symmetric ϕ4 lattice model with periodic boundary conditions and with isotropic short-range interactions. Exact results are derived in the large-n limit describing the geometric crossover from film (ρ=0) over cubic (ρ=1) to cylindrical (ρ=∞) geometries. For n=1, three perturbation approaches in the minimal renormalization scheme at fixed d are presented that cover both the central finite-size regime near T(c) for 1/4≲ρ≲3 and the region well above and below T(c). At bulk T(c), we predict the critical Casimir force in the vertical (L) direction to be negative (attractive) for a slab (ρ<1), positive (repulsive) for a rod (ρ>1), and zero for a cube (ρ=1). Our results for finite-size scaling functions agree well with Monte Carlo data for the three-dimensional Ising model by Hasenbusch for ρ=1 and by Vasilyev et al. for ρ=1/6 above, at, and below T(c).

Boundary conditions and the critical Casimir force on an Ising model film: exact results in one and two dimensions

Finite-size effects in certain critical systems can be understood as universal Casimir forces. Here, we compare the Casimir force for free, fixed, periodic, and antiperiodic boundary conditions in the exactly calculable case of the ferromagnetic Ising model in one and two dimensions. We employ a procedure which allows us to calculate the Casimir force with the aforementioned boundary conditions analytically in a transparent manner. Among other results, we find an attractive Casimir force for the case of periodic boundary conditions and a repulsive Casimir force in the antiperiodic case.

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.

Crossover of critical Casimir forces between different surface universality classes

In confined systems near a continuous phase transition the long-ranged fluctuations of the corresponding order parameter are subject to boundary conditions. These constraints result in so-called critical Casimir forces acting as effective forces on the confining surfaces. For systems belonging to the Ising bulk universality class corresponding to a scalar order parameter the critical Casimir force is studied for the film geometry in the crossover regime characterized by different surface fields at the two surfaces. The scaling function of the critical Casimir force is calculated within mean-field theory. Within our approach, the scaling functions of the critical Casimir force and of the order parameter profile for finite surface fields can be mapped by rescaling, except for a narrow crossover regime, onto the corresponding scaling function of the so-called normal fixed point of strong surface fields. In the crossover regime, the critical Casimir force as function of temperature exhibits more than one extremum and for certain ranges of surface field strengths it changes sign twice upon varying temperature. Monte Carlo simulation data obtained for a three-dimensional Ising film show similar trends. The sign of the critical Casimir force can be inferred from the comparison of the order parameter profiles in the film and in the semi-infinite geometry.

Critical Casimir effect in classical binary liquid mixtures

If a fluctuating medium is confined, the ensuing perturbation of its fluctuation spectrum generates Casimir-like effective forces acting on its confining surfaces. Near a continuous phase transition of such a medium the corresponding order parameter fluctuations occur on all length scales and therefore close to the critical point this effect acquires a universal character, i.e., to a large extent it is independent of the microscopic details of the actual system. Accordingly it can be calculated theoretically by studying suitable representative model systems. We report on the direct measurement of critical Casimir forces by total internal reflection microscopy with femtonewton resolution. The corresponding potentials are determined for individual colloidal particles floating above a substrate under the action of the critical thermal noise in the solvent medium, constituted by a binary liquid mixture of water and 2,6-lutidine near its lower consolute point. Depending on the relative adsorption preferences of the colloid and substrate surfaces with respect to the two components of the binary liquid mixture, we observe that, upon approaching the critical point of the solvent, attractive or repulsive forces emerge and supersede those prevailing away from it. Based on the knowledge of the critical Casimir forces acting in film geometries within the Ising universality class and with equal or opposing boundary conditions, we provide the corresponding theoretical predictions for the sphere-planar wall geometry of the experiment. The experimental data for the effective potential can be interpreted consistently in terms of these predictions and a remarkable quantitative agreement is observed.

Another method to compute the thermodynamic Casimir force in lattice models

We discuss a method that allows us to compute the thermodynamic Casimir force at a given temperature in lattice models by performing a single Monte Carlo simulation. It is analogous to the one used by de Forcrand and co-workers in the study of 't Hooft loops and the interface tension in SU(N) lattice gauge models in four dimensions. We test the method at the example of thin films in the XY universality class. In particular we simulate the improved two-component phi4 model on the simple cubic lattice. This allows us to compare with our previous study, where we have computed the Casimir force by numerically integrating energy densities over the inverse temperature.

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.

Universal scaling functions of critical Casimir forces obtained by Monte Carlo simulations

Effective Casimir forces induced by thermal fluctuations in the vicinity of bulk critical points are studied by means of Monte Carlo simulations in three-dimensional systems for film geometries and within the experimentally relevant Ising and XY universality classes. Several surface universality classes of the confining surfaces are considered, some of which are relevant for recent experiments. An approach introduced previously [O. Vasilyev, EPL 80, 60009 (2007)], based inter alia on an integration scheme of free-energy differences, is utilized to compute the universal scaling functions of the critical Casimir forces in the critical range of temperatures above and below the bulk critical temperature. The resulting predictions are compared with corresponding experimental data for wetting films of fluids and with available theoretical results.

Off-critical Casimir effect in Ising slabs with symmetric boundary conditions in d=3

Extended de Gennes-Fisher (EdGF) local-functional method has been applied to the thermodynamic Casimir effect away from the critical point for systems in the Ising universality class confined between parallel plane plates with symmetric boundary conditions [denoted (ab)=(++)]. Results on the universal scaling functions of the Casimir force W++(y) (y is a temperature-dependent scaling variable) and Gibbs adsorption G[over ](y) are presented in spatial dimension d=3. Also, the mean-field form of the universal scaling function of the Gibbs adsorption G[over ](y) is derived within the local functional theory. Asymptotic behavior of W++(y) for large values of the scaling variable y is analyzed in general dimension d.

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.

Thermodynamic Casimir effect in 4He films near Tlambda: Monte Carlo results

The universal finite-size scaling function of the critical Casimir force for the three dimensional XY universality class with Dirichlet boundary conditions is determined using Monte Carlo simulations. The results are in excellent agreement with recent experiments on 4He Films at the superfluid transition and with available theoretical predictions.

Critical Casimir effect in superfluid wetting films

Recent experimental data for the complete wetting behavior of pure 4He and of 3He-4He mixtures exposed to solid substrates show that there is a change of the corresponding film thicknesses L upon approaching thermodynamically the lambda transition and the tricritical end point, respectively, which can be attributed to critical Casimir forces fC. We calculate the scaling functions theta of fC within models representing the corresponding universality classes. For the mixtures our analysis provides an understanding of the rich behavior of theta deduced from the experimental data and predicts the crossover behavior between the tricritical point and the lambda transition of pure 4He which are connected by a line of critical points. The formation of a "soft-mode" phase within the wetting films gives rise to a pronounced maximum of fC below the tricritical point as observed experimentally. Near the tricritical point we find logarithmic corrections approximately L(-3)(ln L)1/2 for the leading behavior of theta dominating the contributions from the background dispersion forces.

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.

Forces between chemically structured substrates mediated by critical fluids

We consider binary liquid mixtures close to their critical points confined by two parallel, geometrically flat, and chemically structured substrates. Universal order parameter profiles are calculated within mean field theory for periodic patterns of stripes with alternating preferences for the two species of the mixture and with different relative positions of the two substrates. From the order parameter profiles the effective forces between the two plates are derived. The tuning of Casimir amplitudes is discussed.

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.

Universality of the thermodynamic Casimir effect

Recently a nonuniversal character of the leading spatial behavior of the thermodynamic Casimir force has been reported [X. S. Chen and V. Dohm, Phys. Rev. E 66, 016102 (2002)]. We reconsider the arguments leading to this observation and show that there is no such leading nonuniversal term in the systems with short-ranged interactions if one treats properly the effects generated by a sharp momentum cutoff in the Fourier transform of the interaction potential. We also conclude that lattice and continuum models then produce results in mutual agreement independent of the cutoff scheme, contrary to the aforementioned report. All results are consistent with the universal character of the Casimir force in the systems with short-ranged interactions. The effects due to dispersion forces are discussed for the systems with periodic or realistic boundary conditions. In contrast to the systems with short-ranged interactions, for L/xi>>1, one observes leading finite-size contributions governed by power laws in L due to the subleading long-ranged character of the interaction, where L is the finite system size and xi is the correlation length.

Colloidal aggregation in polymer blends

We consider here a low-density assembly of colloidal particles immersed in a critical polymer mixture of two chemically incompatible polymers. We assume that, close to the critical point of the free mixture, the colloids prefer to be surrounded by one polymer (critical adsorption). As result, one is assisted to a reversible colloidal aggregation in the nonpreferred phase, due the existence of a long-range attractive Casimir force between particles. This aggregation is a phase transition driving the colloidal system from dilute to dense phases, as the usual gas-liquid transition. We are interested in a quantitative investigation of the phase diagram of the immersed colloids. We suppose that the positions of particles are disordered, and the disorder is quenched and follows a Gaussian distribution. To apprehend the problem, use is made of the standard phi(4) theory, where the field phi represents the composition fluctuation (order parameter), combined with the standard cumulant method. First, we derive the expression of the effective free energy of colloids and show that this is of Flory-Huggins type. Second, we find that the interaction parameter u between colloids is simply a linear combination of the isotherm compressibility and specific heat of the free mixture. Third, with the help of the derived effective free energy, we determine the complete shape of the phase diagram (binodal and spinodal) in the (Psi,u) plane, with Psi as the volume fraction of immersed colloids. The continuous "gas-liquid" transition occurs at some critical point K of coordinates (Psi(c) = 0.5,u(c) = 2). Finally, we emphasize that the present work is a natural extension of that, relative to simple liquid mixtures incorporating colloids.

Critical Casimir effect in three-dimensional Ising systems: measurements on binary wetting films

The critical Casimir force (CF) is observed in thin wetting films of a binary liquid mixture close to the liquid/vapor coexistence. X-ray reflectivity shows thickness (L) enhancement near the bulk consolute point. The extracted Casimir amplitude Delta(+-)=3+/-1 agrees with the theoretical universal value for the antisymmetric 3D Ising films. The onset of CF in the one-phase region occurs at L/xi approximately 5 regardless of whether the bulk correlation length xi is varied with temperature or composition. The shape of the Casimir scaling function depends monotonically on the dimensionality.

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 .

Nonuniversal finite-size scaling in anisotropic systems

We study the bulk and finite-size critical behavior of the O(n) symmetric phi4 theory with spatially anisotropic interactions of noncubic symmetry in d<4 dimensions. In such systems of a given (d,n) universality class, two-scale factor universality is absent in bulk correlation functions, and finite-size scaling functions including the Privman-Fisher scaling form of the free energy, the Binder cumulant ratio, and the Casimir amplitude are shown to be nonuniversal. In particular it is shown that, for anisotropic confined systems, isotropy cannot be restored by an anisotropic scale transformation.

Critical adsorption and Casimir torque in wedges and at ridges

Geometrical structures of confining surfaces profoundly influence the adsorption of fluids upon approaching a critical point T(c) in their bulk phase diagram, i.e., for t= (T- T(c) ) / T(c) -->+/-0 . Guided by general scaling considerations, we calculate, within the mean-field theory, the temperature dependence of the order parameter profile in a wedge with opening angle gamma<pi and close to a ridge (gamma>pi) for T>< T(c) and in the presence of surface fields. For a suitably defined reduced excess adsorption Gamma (+/-) (gamma,t-->+/-0) approximately Gamma (+/-) (gamma) mid R:tmid R:(beta-2nu) we compute the universal amplitudes Gamma (+/-) (gamma) , which diverge as Gamma (+/-) (gamma-->0) approximately 1/gamma for small opening angles, vary linearly close to gamma=pi for gamma<pi , and increase exponentially for gamma-->2pi . There is evidence that, within the mean-field theory, the ratio Gamma (+) (gamma) / Gamma (-) (gamma) is independent of gamma . We also discuss the critical Casimir torque acting on the sides of the wedge as a function of the opening angle and temperature.