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.
We study critical Casimir forces (CCF) $f_{mathrm C}$ 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 (RSF) on both surfaces. We consider the case that, in the absence of RSF, 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 CCF still exhibit scaling, acquiring a random field scaling variable $w$ which is zero for pure systems. We confirm these analytic predictions by MC simulations. Moreover, our MC data show that $f_{mathrm C}$ varies as $f_{mathrm C}(wto 0)-f_{mathrm C}(w=0)sim w^2$. Asymptotically, for large $L$, $w$ scales as $w sim L^{-0.26} to 0$ indicating that this type of disorder is an irrelevant perturbation of the ordinary surface universality class. However, for thin films such that $w simeq 1$, we find that the presence of RSF with vanishing mean value increases significantly the strength of CCF, as compared to systems without them, and shifts the extremum of the scaling function of $f_{mathrm C}$ towards lower temperatures. But $f_{mathrm C}$ remains attractive.
We present a new Monte Carlo method to calculate Casimir forces acting on objects in a near-critical fluid, considering the two basic cases of a wall and a sphere embedded in a two-dimensional Ising medium. During the simulation, the objects are moved through the system with appropriate statistical weights, and consequently are attracted or repelled from the system boundaries depending on the boundary conditions. The distribution function of the object position is utilized to obtain the residual free energy, or Casimir potential, of the configuration as well as the corresponding Casimir force. The results are in perfect agreement with known exact results. The method can easily be generalized to more complicated geometries, to higher dimensions, and also to colloidal suspensions with many particles.
The confinement of critical fluctuations in soft media induces critical Casimir forces acting on the confining surfaces. The temperature and geometry dependences of such forces are characterized by universal scaling functions. A novel approach is presented to determine them for films via Monte Carlo simulations of lattice models. The method is based on an integration scheme of free energy differences. Our results for the Ising and the XY universality class compare favourably with corresponding experimental results for wetting layers of classical binary liquid mixtures and of 4He, respectively.
Using general scaling arguments combined with mean-field theory we investigate the critical ($T simeq T_c$) and off-critical ($T e T_c$) 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 $T_c$. If, however, the direct substrate-substrate contribution is subtracted, the force is repulsive everywhere except near the bulk critical point $(T_c,mu_c)$, where critical density fluctuations arise, or except at low temperatures and $(L/a) (betaDelta 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}$. Essential deviations from the standard finite-size scaling behavior are observed within the finite-size critical region $L/xi=O(1)$ for films with thicknesses $L lesssim L_{rm crit}$, where $L_{rm crit}=xi_0^pm (16 |s|)^{ u/beta}$, with $ u$ 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.
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 f_C. We calculate the scaling functions vartheta of f_C within models representing the corresponding universality classes. For the mixtures our analysis provides an understanding of the rich behavior of vartheta 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 f_C below the tricritical point as observed experimentally. Near the tricritical point we find logarithmic corrections ~L^(-3)(ln L)^(1/2) for the leading behavior of vartheta dominating the contributions from the background dispersion forces.