Critical Casimir forces between homogeneous and chemically striped surfaces

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

Boundary Criticality of the 3D O(N) Model: From Normal to Extraordinary

It was recently realized that the three-dimensional O(N) model possesses an extraordinary boundary universality class for a finite range of N≥2. For a given N, the existence and universal properties of this class are predicted to be controlled by certain amplitudes of the normal universality class, where one applies an explicit symmetry breaking field to the boundary. In this Letter, we study the normal universality class for N=2, 3 using Monte Carlo simulations on an improved lattice model and extract these universal amplitudes. Our results are in good agreement with direct Monte Carlo studies of the extraordinary universality class serving as a nontrivial quantitative check of the connection between the normal and extraordinary classes.

Boundary Critical Behavior of the Three-Dimensional Heisenberg Universality Class

We study the boundary critical behavior of the three-dimensional Heisenberg universality class, in the presence of a bidimensional surface. By means of high-precision Monte Carlo simulations of an improved lattice model, where leading bulk scaling corrections are suppressed, we prove the existence of a special phase transition, with unusual exponents, and of an extraordinary phase with logarithmically decaying correlations. These findings contrast with naïve arguments on the bulk-surface phase diagram, and allow us to explain some recent puzzling results on the boundary critical behavior of quantum spin models.

Two-dimensional critical systems with mixed boundary conditions: Exact Ising results from conformal invariance and boundary-operator expansions

With conformal-invariance methods, Burkhardt, Guim, and Xue studied the critical Ising model, defined on the upper half plane y>0 with different boundary conditions a and b on the negative and positive x axes. For ab=-+ and f+, they determined the one- and two-point averages of the spin σ and energy ε. Here +,-, and f stand for spin-up, spin-down, and free-spin boundaries, respectively. The case +-+-+⋯, where the boundary condition switches between + and - at arbitrary points, ζ_{1},ζ_{2},⋯ on the x axis was also analyzed. In the first half of this paper a similar study is carried out for the alternating boundary condition +f+f+⋯ and the case -f+ of three different boundary conditions. Exact results for the one- and two-point averages of σ,ε, and the stress tensor T are derived with conformal-invariance methods. From the results for 〈T〉, the critical Casimir interaction with the boundary of a wedge-shaped inclusion is derived for mixed boundary conditions. In the second half of the paper, arbitrary two-dimensional critical systems with mixed boundary conditions are analyzed with boundary-operator expansions. Two distinct types of expansions-away from switching points of the boundary condition and at switching points-are considered. Using the expansions, we express the asymptotic behavior of two-point averages near boundaries in terms of one-point averages. We also consider the strip geometry with mixed boundary conditions and derive the distant-wall corrections to one-point averages near one edge due to the other edge. Finally we confirm the consistency of the predictions obtained with conformal-invariance methods and with boundary-operator expansions, in the the first and second halves of the paper.

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.

On the critical Casimir interaction between anisotropic inclusions on a membrane

Using a lattice model and a versatile thermodynamic integration scheme, we study the critical Casimir interactions between inclusions embedded in a two-dimensional critical binary mixtures. For single-domain inclusions we demonstrate that the interactions are very long range, and their magnitudes strongly depend on the affinity of the inclusions with the species in the binary mixtures, ranging from repulsive when two inclusions have opposing affinities to attractive when they have the same affinities. When one of the inclusions has no preference for either of the species, we find negligible critical Casimir interactions. For multiple-domain inclusions, mimicking the observations that membrane proteins often have several domains with varying affinities to the surrounding lipid species, the presence of domains with opposing affinities does not cancel the interactions altogether. Instead we can observe both attractive and repulsive interactions depending on their relative orientations. With increasing number of domains per inclusion, the range and magnitude of the effective interactions decrease in a similar fashion to those of electrostatic multipoles. Finally, clusters formed by multiple-domain inclusions can result in an effective affinity patterning due to the anisotropic character of the Casimir interactions between the building blocks.

Manipulation and amplification of the Casimir force through surface fields using helicity

We present both exact and numerical results for the behavior of the Casimir force in O(n) systems with a finite extension L in one direction when the system is subjected to surface fields that induce helicity in the order parameter. We show that for such systems, the Casimir force in certain temperature ranges is of the order of L^{-2}, both above and below the critical temperature, T_{c}, of the bulk system. An example of such a system would be one with chemically modulated bounding surfaces, in which the modulation couples directly to the system's order parameter. We demonstrate that, depending on the parameters of the system, the Casimir force can be either attractive or repulsive. The exact calculations presented are for the one-dimensional XY and Heisenberg models under twisted boundary conditions resulting from finite surface fields that differ in direction by a specified angle, and the three-dimensional Gaussian model with surface fields in the form of plane waves that are shifted in phase with respect to each other. Additionally, we present exact and numerical results for the mean-field version of the three-dimensional O(2) model with finite surface fields on the bounding surfaces. We find that all significant results are consistent with the expectations of finite-size scaling.

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.

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 force in the presence of random local adsorption preference

We study the critical Casimir force for a film geometry in the Ising universality class. We employ a homogeneous adsorption preference on one of the confining surfaces, while the opposing surface exhibits quenched random disorder, leading to a random local adsorption preference. Disorder is characterized by a parameter p, which measures, on average, the portion of the surface that prefers one component, so that p=0,1 correspond to homogeneous adsorption preference. By means of Monte Carlo simulations of an improved Hamiltonian and finite-size scaling analysis, we determine the critical Casimir force. We show that by tuning the disorder parameter p, the system exhibits a crossover between an attractive and a repulsive force. At p=1/2, disorder allows to effectively realize Dirichlet boundary conditions, which are generically not accessible in classical fluids. Our results are relevant for the experimental realizations of the critical Casimir force in binary liquid mixtures.

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.

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.