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

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.

Exact expressions for the partition function of the one-dimensional Ising model in the fixed-M ensemble

We obtain exact closed-form expressions for the partition function of the one-dimensional Ising model in the fixed-M ensemble, for three commonly used boundary conditions: periodic, antiperiodic, and Dirichlet. These expressions allow for the determination of fluctuation-induced forces in the canonical ensemble, which we term Helmholtz forces. The thermodynamic expressions and the calculations flowing from them should provide insights into the nature and behavior of fluctuation-induced forces in interesting and as-yet unexplored regimes.

Exact Critical Casimir Amplitude of Anisotropic Systems from Conformal Field Theory and Self-Similarity of Finite-Size Scaling Functions in d≥2 Dimensions

The exact critical Casimir amplitude is derived for anisotropic systems within the d=2 Ising universality class by combining conformal field theory with anisotropic φ^{4} theory. Explicit results are presented for the general anisotropic scalar φ^{4} model and for the fully anisotropic triangular-lattice Ising model in finite rectangular and infinite strip geometries with periodic boundary conditions. These results demonstrate the validity of multiparameter universality for confined anisotropic systems and the nonuniversality of the critical Casimir amplitude. We find an unexpected complex form of self-similarity of the anisotropy effects near the instability where weak anisotropy breaks down. This can be traced back to the property of modular invariance of isotropic conformal field theory for d=2. More generally, for d>2 we predict the existence of self-similar structures of the finite-size scaling functions of O(n)-symmetric systems with planar anisotropies and periodic boundary conditions both in the critical region for n≥1 as well as in the Goldstone-dominated low-temperature region for n≥2.

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.

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.

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.

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.

Universal anisotropic finite-size critical behavior of the two-dimensional Ising model on a strip and of d-dimensional models on films

Anisotropy effects on the finite-size critical behavior of a two-dimensional Ising model on a general triangular lattice in an infinite-strip geometry with periodic, antiperiodic, and free boundary conditions (bc) in the finite direction are investigated. Exact results are obtained for the scaling functions of the finite-size contributions to the free energy density. With ξ(>) the largest and ξ(<) the smallest bulk correlation length at a given temperature near criticality, we find that the dependence of these functions on the ratio ξ(<)/ξ(>) and on the angle parametrizing the orientation of the correlation volume is of geometric nature. Since the scaling functions are independent of the particular microscopic realization of the anisotropy within the two-dimensional Ising model, our results provide a limited verification of universality. We explain our observations by considering finite-size scaling of free energy densities of general weakly anisotropic models on a d-dimensional film (i.e., in an L×∞(d-1) geometry) with bc in the finite direction that are invariant under a shear transformation relating the anisotropic and isotropic cases. This allows us to relate free energy scaling functions in the presence of an anisotropy to those of the corresponding isotropic system. We interpret our results as a simple and transparent case of anisotropic universality, where, compared to the isotropic case, scaling functions depend additionally on the shape and orientation of the correlation volume. We conjecture that this universality extends to cases where the geometry and/or the bc are not invariant under the shear transformation and argue in favor of validity of two-scale factor universality for weakly anisotropic systems.

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

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.