No Arabic abstract
We present a systematic method to calculate the universal scaling functions for the critical Casimir force and the according potential of the two-dimensional Ising model with various boundary conditions. Therefore we start with the dimer representation of the corresponding partition function $Z$ on an $Ltimes M$ square lattice, wrapped around a torus with aspect ratio $rho=L/M$. By assuming periodic boundary conditions and translational invariance in at least one direction, we systematically reduce the problem to a $2times2$ transfer matrix representation. For the torus we first reproduce the results by Kaufman and then give a detailed calculation of the scaling functions. Afterwards we present the calculation for the cylinder with open boundary conditions. All scaling functions are given in form of combinations of infinite products and integrals. Our results reproduce the known scaling functions in the limit of thin films $rhoto 0$. Additionally, for the cylinder at criticality our results confirm the predictions from conformal field theory.
Based on the results published recently [SciPost Phys. 7, 026 (2019)], the influence of surfaces and boundary fields are calculated for the ferromagnetic anisotropic square lattice Ising model on finite lattices as well as in the finite-size scaling limit. Starting with the open cylinder, we independently apply boundary fields on both sides which can be either homogeneous or staggered, representing different combinations of boundary conditions. We confirm several predictions from scaling theory, conformal field theory and renormalisation group theory: we explicitly show that anisotropic couplings enter the scaling functions through a generalised aspect ratio, and demonstrate that open and staggered boundary conditions are asymptotically equal in the scaling regime. Furthermore, we examine the emergence of the surface tension due to one antiperiodic boundary in the system in the presence of symmetry breaking boundary fields, again for finite systems as well as in the scaling limit. Finally, we extend our results to the antiferromagnetic Ising model.
We have dramatically extended the zero field susceptibility series at both high and low temperature of the Ising model on the triangular and honeycomb lattices, and used these data and newly available further terms for the square lattice to calculate a number of terms in the scaling function expansion around both the ferromagnetic and, for the square and honeycomb lattices, the antiferromagnetic critical point.
We present general arguments and construct a stress tensor operator for finite lattice spin models. The average value of this operator gives the Casimir force of the system close to the bulk critical temperature $T_c$. We verify our arguments via exact results for the force in the two-dimensional Ising model, $d$-dimensional Gaussian and mean spherical model with $2<d<4$. On the basis of these exact results and by Monte Carlo simulations for three-dimensional Ising, XY and Heisenberg models we demonstrate that the standard deviation of the Casimir force $F_C$ in a slab geometry confining a critical substance in-between is $k_b T D(T)(A/a^{d-1})^{1/2}$, where $A$ is the surface area of the plates, $a$ is the lattice spacing and $D(T)$ is a slowly varying nonuniversal function of the temperature $T$. The numerical calculations demonstrate that at the critical temperature $T_c$ the force possesses a Gaussian distribution centered at the mean value of the force $<F_C>=k_b T_c (d-1)Delta/(L/a)^{d}$, where $L$ is the distance between the plates and $Delta$ is the (universal) Casimir amplitude.
The critical Casimir force (CCF) arises from confining fluctuations in a critical fluid and thus it is a fluctuating quantity itself. While the mean CCF is universal, its (static) variance has previously been found to depend on the microscopic details of the system which effectively set a large-momentum cutoff in the underlying field theory, rendering it potentially large. This raises the question how the properties of the force variance are reflected in experimentally observable quantities, such as the thickness of a wetting film or the position of a suspended colloidal particle. Here, based on a rigorous definition of the instantaneous force, we analyze static and dynamic correlations of the CCF for a conserved fluid in film geometry for various boundary conditions within the Gaussian approximation. We find that the dynamic correlation function of the CCF is independent of the momentum cutoff and decays algebraically in time. Within the Gaussian approximation, the associated exponent depends only on the dynamic universality class but not on the boundary conditions. We furthermore consider a fluid film, the thickness of which can fluctuate under the influence of the time-dependent CCF. The latter gives rise to an effective non-Markovian noise in the equation of motion of the film boundary and induces a distinct contribution to the position variance. Within the approximations used here, at short times, this contribution grows algebraically in time whereas, at long times, it saturates and contributes to the steady-state variance of the film thickness.
We present results of a Monte Carlo study for the ferromagnetic Ising model with long range interactions in two dimensions. This model has been simulated for a large range of interaction parameter $sigma$ and for large sizes. We observe that the results close to the change of regime from intermediate to short range do not agree with the renormalization group predictions.