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.