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

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.

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.

Critical Casimir interactions between colloids around the critical point of binary solvents

Critical Casimir interactions between colloidal particles arise from the confinement of fluctuations of a near-critical solvent in the liquid gap between closely-spaced particles. So far, the comparison of theoretical predictions and experimental measurements of critical Casimir forces (CCFs) has focused on the critical solvent composition, while it has been lacking for off-critical compositions. We address this issue by investigating CCFs between spherical colloidal particles around the critical point of a binary solvent through a combination of experiments, previous Ising Monte Carlo simulation results and field-theoretical methods. By measuring the correlation length of the near-critical solvent and the pair potentials of the particles in terms of radial distribution functions and by determining the second virial coefficient, we test in detail theoretical predictions. Our results indicate that the critical Casimir theory gives quantitative correct predictions for the interaction potential between particles in a near critical binary mixture if weak preferential adsorption of the particle surface is taken into account.

Solvent-mediated forces in critical fluids

The effective interaction between two planar walls immersed in a fluid is investigated by use of density functional theory in the supercritical region of the phase diagram. A hard core Yukawa model of fluid is studied with special attention to the critical region. To achieve this goal a formulation of the weighted density approximation coupled with the hierarchical reference theory, able to deal with critical long wavelength fluctuations, is put forward and compared with other approaches. The effective interaction between the walls is seen to change character on lowering the temperature: The strong oscillations induced by layering of the molecules, typical of the depletion mechanism in hard core systems, are gradually smoothed and, close to the critical point, a long range attractive tail emerges leading to a scaling form which agrees with the expectations based on the critical Casimir effect. Strong corrections to scaling are seen to affect the results up to very small reduced temperatures. By use of the Derjaguin approximation, this investigation has natural implications for the aggregation of colloidal particles in critical solvents.

Critical adsorption and critical Casimir forces in the canonical ensemble

Critical properties of a liquid film between two planar walls are investigated in the canonical ensemble, within which the total number of fluid particles, rather than their chemical potential, is kept constant. The effect of this constraint is analyzed within mean-field theory (MFT) based on a Ginzburg-Landau free-energy functional as well as via Monte Carlo simulations of the three-dimensional Ising model with fixed total magnetization. Within MFT and for finite adsorption strengths at the walls, the thermodynamic properties of the film in the canonical ensemble can be mapped exactly onto a grand canonical ensemble in which the corresponding chemical potential plays the role of the Lagrange multiplier associated with the constraint. However, due to a nonintegrable divergence of the mean-field order parameter profile near a wall, the limit of infinitely strong adsorption turns out to be not well-defined within MFT, because it would necessarily violate the constraint. The critical Casimir force (CCF) acting on the two planar walls of the film is generally found to behave differently in the canonical and grand canonical ensembles. For instance, the canonical CCF in the presence of equal preferential adsorption at the two walls is found to have the opposite sign and a slower decay behavior as a function of the film thickness compared to its grand canonical counterpart. We derive the stress tensor in the canonical ensemble and find that it has the same expression as in the grand canonical case, but with the chemical potential playing the role of the Lagrange multiplier associated with the constraint. The different behavior of the CCF in the two ensembles is rationalized within MFT by showing that, for a prescribed value of the thermodynamic control parameter of the film, i.e., density or chemical potential, the film pressures are identical in the two ensembles, while the corresponding bulk pressures are not.

Critical Casimir forces for colloidal assembly

Critical Casimir forces attract increasing interest due to their opportunities for reversible particle assembly in soft matter and nano science. These forces provide a thermodynamic analogue of the celebrated quantum mechanical Casimir force that arises from the confinement of vacuum fluctuations of the electromagnetic field. In its thermodynamic analogue, solvent fluctuations, confined between suspended particles, give rise to an attractive or repulsive force between the particles. Due to its unique temperature dependence, this effect allows in situ control of reversible assembly. Both the force magnitude and range vary with the solvent correlation length in a universal manner, adjusting with temperature from fractions of the thermal energy, k B T, and nanometre range to several ten kT and micrometer length scale. Combined with recent breakthroughs in the synthesis of complex particles, critical Casimir forces promise the design and assembly of complex colloidal structures, for fundamental studies of equilibrium and out-of-equilibrium phase behaviour. This review highlights recent developments in this evolving field, with special emphasis on the dynamic interaction control to assemble colloidal structures, in and out of equilibrium.

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 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.

Spin models in three dimensions: Adaptive lattice spacing

Aiming at the study of critical phenomena in the presence of boundaries with a nontrivial shape we discuss how lattices with an adaptive lattice spacing can be implemented. Since the parameters of the Hamiltonian transform nontrivially under changes of the length scale, adapting the lattice spacing is much more difficult than in the case of the numerical solution of partial differential equations, where this method is common practice. Here we shall focus on the universality class of the three-dimensional Ising model. Our starting point is the improved Blume-Capel model on the simple cubic lattice. In our approach, the system is composed of sectors with lattice spacing a,2a,4a,.... We work out how parts of the lattice with lattice spacing a and 2a, respectively, can be coupled in a consistent way. Here, we restrict ourself to the case where the boundary between the sectors is perpendicular to one of the lattice axes. Based on the theory of defect planes one expects that it is sufficient to tune the coupling between these two regions. To this end we perform a finite-size scaling study. However, first numerical results show that slowly decaying corrections remain. It turns out that these corrections can be removed by adjusting the strength of the couplings within the boundary layers. As a benchmark, we simulate films with strongly symmetry breaking boundary conditions. We determine the magnetization profile and the thermodynamic Casimir force. For our largest thickness L0=64.5, we find that results obtained for the homogeneous system are nicely reproduced.

Crossover aspects in Ising strips under the influence of variable surface fields and a grain boundary

I use an exact variational formulation of Mikheev and Fisher to study the critical Ising strip with a grain boundary and confining surfaces characterized by arbitrary and different surface magnetic fields. Energy density profiles that serve as order parameters of the system within the used method show strong nonmonotonous behavior in the vicinity of confining surfaces. I consider short-distance expansion of energy density profiles. New universal amplitudes associated with distant-wall corrections exhibit nontrivial crossover behaviors from positive to negative values as the field's variables are continuously varied. Casimir amplitudes calculated in a self-contained manner are characterized by complex manifolds which may comprise two disconnected positive wings separated by an area of negative values. I define and determine the generalized de Gennes-Fisher amplitude, which strongly suggests that a stress tensor is present within one of the distant-wall corrections. This is an unanticipated result given that a similar discovery for standard extraordinary (E) and ordinary (O) surface universality classes was based on the conformal invariance symmetry, which is broken under the present boundary conditions. A grain boundary essentially influences both energy density profiles and Casimir amplitudes, besides surface fields with a number of accompanying features that are examined in detail. We present closed analytic forms of energy density profiles for standard (N) and (O) BCs: ɛNN(x,g),ɛOO(x,g),ɛNO(x,g) [(N) is the normal surface universality class, g is the grain boundary's strength, x is the relative distance x:=z/L with L as a film width], thus enabling us to deduce their important symmetry properties and to have a detailed insight into their general behaviors.

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.

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.

Attractive interaction and bridging transition between neutral colloidal particles due to preferential adsorption in a near-critical binary mixture

We examine the solvent-mediated interaction between two neutral colloidal particles due to preferential adsorption in a near-critical binary mixture. We take into account the renormalization effect due to the critical fluctuations using the recent local functional theory [J. Chem. Phys. 136, 114704 (2012)]. We calculate the free energy and the force between two colloidal particles as functions of the temperature T, the composition far from the colloidal particles c(∞), and the colloid separation ℓ. The interaction is much enhanced when the component favored by the colloid surfaces is poor in the reservoir. For such off-critical compositions, we find a surface of a first-order bridging transition ℓ=ℓ(cx)(T,c(∞)) in the T-c(∞)-ℓ space in a universal, scaled form, across which a discontinuous change occurs between separated and bridged states. This surface starts from the bulk coexistence surface (CX) and ends at a bridging critical line where ℓ is determined by T as ℓ=ℓ(c)(T). On approaching the critical line, the discontinuity vanishes and the derivatives of the force with respect to T and ℓ both diverge. Furthermore, bridged states continuously change into separated states if c(∞) (or T) is varied from a value on CX to a value far from CX with ℓ kept smaller than ℓ(c)(T).

Phase behavior of colloidal suspensions with critical solvents in terms of effective interactions

We study the phase behavior of colloidal suspensions the solvents of which are considered to be binary liquid mixtures undergoing phase segregation. We focus on the thermodynamic region close to the critical point of the accompanying miscibility gap. There, due to the colloidal particles acting as cavities in the critical medium, the spatial confinements of the critical fluctuations of the corresponding order parameter result in the effective, so-called critical Casimir forces between the colloids. Employing an approach in terms of effective, one-component colloidal systems, we explore the possibility of phase coexistence between two phases of colloidal suspensions, one being rich and the other being poor in colloidal particles. The reliability of this effective approach is discussed.

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 depletion

Depletion interactions and the critical Casimir effect are usually regarded as distinct phenomena in colloidal suspensions. By experimentally investigating how the Asakura-Oosawa picture, appropriate for a weakly correlated depletant, is modified when critical correlations develop within the depletion agent, we conversely show that the former merges continuously into the latter, leading to a distinctive scaling behavior solely dictated by the depletant correlation length. A model based on density functional theory provides a microscopic understanding of the phenomenon and properly accounts for the observed trends.

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.

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.

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.