Thermodynamic Casimir Effect in Films: the Exchange Cluster Algorithm

Martin Hasenbusch
DOI: https://doi.org/10.1103/PhysRevE.91.022110
2014-10-27
Abstract: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 [J.R. Heringa and H. W. J. Blöte, Phys. Rev. E 57, 4976 (1998)]. In a previous work we demonstrated that this algorithm allows 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.
Statistical Mechanics
What problem does this paper attempt to address?
This paper aims to study the thermodynamic Casimir force of thin films under different boundary conditions, especially when the bulk system belongs to the universality class of the three - dimensional Ising model. Specifically, the author carried out Monte Carlo simulations on a simple cubic lattice using an improved Blume - Capel model, especially adopting the exchange or geometric cluster algorithm. The focus of the paper is on studying thin films with ordinary surface transition boundaries, that is, thin films under free boundary conditions. Such thin films undergo a phase transition in the universality class of the two - dimensional Ising model. The author determined the inverse transition temperature within the range of film thickness and studied the variation law of this temperature with the film thickness. In addition, the author also explored the influence of finite - size effects on the thermodynamic Casimir force near the critical point and showed that these finite - size effects can be calculated by the universal finite - size scaling function of the free energy of the two - dimensional Ising model. The main contributions of the paper are as follows: 1. **Method innovation**: By adopting the exchange cluster algorithm, the statistical error is significantly reduced and the efficiency of calculating the thermodynamic Casimir force is improved. 2. **Theoretical verification**: The properties of thin films with ordinary surface transition boundaries are studied in detail, and the relationship between finite - size effects and the two - dimensional Ising model is verified. 3. **Data accuracy**: High - precision numerical results are provided, which support the theoretical interpretation of experimental data. In conclusion, this paper solves the problem of how to efficiently calculate the thermodynamic Casimir force of thin films under different boundary conditions and provides a new perspective for understanding the influence of finite - size effects on the thermodynamic Casimir force.