The critical Casimir force and its fluctuations in lattice spin models: exact and Monte Carlo results

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 $=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.Comment: 21 pages, 7 figures, to appear in PR

Critical Casimir amplitudes for nn-component ϕ4\phi^4 models with O(n)-symmetry breaking quadratic boundary terms

Euclidean $n$-component $\phi^4$ theories whose Hamiltonians are O(n) symmetric except for quadratic symmetry breaking boundary terms are studied in films of thickness $L$. The boundary terms imply the Robin boundary conditions $\partial_n\phi_\alpha =\mathring{c}^{(j)}_\alpha \phi_\alpha$ at the boundary planes $\mathfrak{B}_{j=1,2}$ at $z=0$ and $z=L$. Particular attention is paid to the cases in which $m_j$ of the $n$ variables $\mathring{c}^{(j)}_\alpha$ take the special value $\mathring{c}_{m_j\text{-sp}}$ corresponding to critical enhancement while the remaining ones are subcritically enhanced. Under these conditions, the semi-infinite system bounded by $\mathfrak{B}_j$ has a multicritical point, called $m_j$-special, at which an $O(m_j)$ symmetric critical surface phase coexists with the O(n) symmetric bulk phase, provided $d$ is sufficiently large. The $L$-dependent part of the reduced free energy per area behaves as $\Delta_C/L^{d-1}$ as $L\to\infty$ at the bulk critical point. The Casimir amplitudes $\Delta_C$ are determined for small $\epsilon=4-d$ in the general case where $m_{c,c}$ components $\phi_\alpha$ are critically enhanced at both boundary planes, $m_{c,D} + m_{D,c}$ components are enhanced at one plane but satisfy asymptotic Dirichlet boundary conditions at the respective other, and the remaining $m_{D,D}$ components satisfy asymptotic Dirichlet boundary conditions at both $\mathfrak{B}_j$. Whenever $m_{c,c}>0$, these expansions involve integer and fractional powers $\epsilon^{k/2}$ with $k\ge 3$ (mod logarithms). Results to $O(\epsilon^{3/2})$ for general values of $m_{c,c}$, $m_{c,D}+m_{D,c}$, and $m_{D,D}$ are used to estimate the $\Delta_C$ of 3D Heisenberg systems with surface spin anisotropies when $(m_{c,c}, m_{c,D}+ m_{D,c}) = (1,0)$, $(0,1)$, and $(1,1)$.Comment: Latex source file with 5 eps files; version with minor amendments and corrected typo

Universality of the thermodynamic Casimir effect

Recently a nonuniversal character of the leading spatial behavior of the thermodynamic Casimir force has been reported [X. S. Chen and V. Dohm, Phys. Rev. E {\bf 66}, 016102 (2002)]. We reconsider the arguments leading to this observation and show that there is no such leading nonuniversal term in systems with short-ranged interactions if one treats properly the effects generated by a sharp momentum cutoff in the Fourier transform of the interaction potential. We also conclude that lattice and continuum models then produce results in mutual agreement independent of the cutoff scheme, contrary to the aforementioned report. All results are consistent with the {\em universal} character of the Casimir force in systems with short-ranged interactions. The effects due to dispersion forces are discussed for systems with periodic or realistic boundary conditions. In contrast to systems with short-ranged interactions, for $L/\xi \gg 1$ one observes leading finite-size contributions governed by power laws in $L$ due to the subleading long-ranged character of the interaction, where $L$ is the finite system size and $\xi$ is the correlation length.Comment: 11 pages, revtex, to appear in Phys. Rev. E 68 (2003

Casimir force in the rotor model with twisted boundary conditions

We investigate the three dimensional lattice XY model with nearest neighbor interaction. The vector order parameter of this system lies on the vertices of a cubic lattice, which is embedded in a system with a film geometry. The orientations of the vectors are fixed at the two opposite sides of the film. The angle between the vectors at the two boundaries is $\alpha$ where $0 \le \alpha \le \pi$. We make use of the mean field approximation to study the mean length and orientation of the vector order parameter throughout the film---and the Casimir force it generates---as a function of the temperature $T$, the angle $\alpha$, and the thickness $L$ of the system. Among the results of that calculation are a Casimir force that depends in a continuous way on both the parameter $\alpha$ and the temperature and that can be attractive or repulsive. In particular, by varying $\alpha$ and/or $T$ one controls \underline{both} the sign \underline{and} the magnitude of the Casimir force in a reversible way. Furthermore, for the case $\alpha=\pi$, we discover an additional phase transition occurring only in the finite system associated with the variation of the orientations of the vectors.Comment: 14 pages, 9 figure

Thermodynamic Casimir effects involving interacting field theories with zero modes

Systems with an O(n) symmetrical Hamiltonian are considered in a $d$-dimensional slab geometry of macroscopic lateral extension and finite thickness $L$ that undergo a continuous bulk phase transition in the limit $L\to\infty$. The effective forces induced by thermal fluctuations at and above the bulk critical temperature $T_{c,\infty}$ (thermodynamic Casimir effect) are investigated below the upper critical dimension $d^*=4$ by means of field-theoretic renormalization group methods for the case of periodic and special-special boundary conditions, where the latter correspond to the critical enhancement of the surface interactions on both boundary planes. As shown previously [\textit{Europhys. Lett.} \textbf{75}, 241 (2006)], the zero modes that are present in Landau theory at $T_{c,\infty}$ make conventional RG-improved perturbation theory in $4-\epsilon$ dimensions ill-defined. The revised expansion introduced there is utilized to compute the scaling functions of the excess free energy and the Casimir force for temperatures T\geqT_{c,\infty} as functions of $\mathsf{L}\equiv L/\xi_\infty$, where $\xi_\infty$ is the bulk correlation length. Scaling functions of the $L$-dependent residual free energy per area are obtained whose $\mathsf{L}\to0$ limits are in conformity with previous results for the Casimir amplitudes $\Delta_C$ to $O(\epsilon^{3/2})$ and display a more reasonable small-$\mathsf{L}$ behavior inasmuch as they approach the critical value $\Delta_C$ monotonically as $\mathsf{L}\to 0$.Comment: 23 pages, 10 figure

Exact Three Dimensional Casimir Force Amplitude, CC-function and Binder's Cumulant Ratio: Spherical Model Results

The three dimensional mean spherical model on a hypercubic lattice with a film geometry $L\times \infty ^2$ under periodic boundary conditions is considered in the presence of an external magnetic field $H$. The universal Casimir amplitude $\Delta$ and the Binder's cumulant ratio $B$ are calculated exactly and found to be $\Delta =-2\zeta (3)/(5\pi)\approx -0.153051$ and $B=2\pi /(\sqrt{5}\ln ^3[(1+\sqrt{5})/2]).$ A discussion on the relations between the finite temperature $C$-function, usually defined for quantum systems, and the excess free energy (due to the finite-size contributions to the free energy of the system) scaling function is presented. It is demonstrated that the $C$-function of the model equals 4/5 at the bulk critical temperature $T_c$. It is analytically shown that the excess free energy is a monotonically increasing function of the temperature $T$ and of the magnetic field $|H|$ in the vicinity of $T_c.$ This property is supposed to hold for any classical $d$-dimensional $O(n),n>2,$ model with a film geometry under periodic boundary conditions when $d\leq 3$. An analytical evidence is also presented to confirm that the Casimir force in the system is negative both below and in the vicinity of the bulk critical temperature $T_c.$Comment: 12 pages revtex, one eps figure, submitted to Phys. Rev E A set of references added with the text needed to incorporate them. Small changes in the title and in the abstrac

Excess free energy and Casimir forces in systems with long-range interactions of van-der-Waals type: General considerations and exact spherical-model results

We consider systems confined to a $d$-dimensional slab of macroscopic lateral extension and finite thickness $L$ that undergo a continuous bulk phase transition in the limit $L\to\infty$ and are describable by an O(n) symmetrical Hamiltonian. Periodic boundary conditions are applied across the slab. We study the effects of long-range pair interactions whose potential decays as $b x^{-(d+\sigma)}$ as $x\to\infty$, with $2<\sigma<4$ and $2<d+\sigma\leq 6$, on the Casimir effect at and near the bulk critical temperature $T_{c,\infty}$, for $2<d<4$. For the scaled reduced Casimir force per unit cross-sectional area, we obtain the form L^{d} {\mathcal F}_C/k_BT \approx \Xi_0(L/\xi_\infty) + g_\omega L^{-\omega}\Xi\omega(L/\xi_\infty) + g_\sigma L^{-\omega_\sigm a} \Xi_\sigma(L \xi_\infty). The contribution $\propto g_\sigma$ decays for $T\neq T_{c,\infty}$ algebraically in $L$ rather than exponentially, and hence becomes dominant in an appropriate regime of temperatures and $L$. We derive exact results for spherical and Gaussian models which confirm these findings. In the case $d+\sigma =6$, which includes that of nonretarded van-der-Waals interactions in $d=3$ dimensions, the power laws of the corrections to scaling $\propto b$ of the spherical model are found to get modified by logarithms. Using general RG ideas, we show that these logarithmic singularities originate from the degeneracy $\omega=\omega_\sigma=4-d$ that occurs for the spherical model when $d+\sigma=6$, in conjunction with the $b$ dependence of $g_\omega$.Comment: 28 RevTeX pages, 12 eps figures, submitted to PR

Casimir force in O(n) lattice models with a diffuse interface

On the example of the spherical model we study, as a function of the temperature $T$, the behavior of the Casimir force in O(n) systems with a diffuse interface and slab geometry $\infty^{d-1}\times L$, where $2<d<4$ is the dimensionality of the system. We consider a system with nearest-neighbor anisotropic interaction constants $J_\parallel$ parallel to the film and $J_\perp$ across it. The model represents the $n\to\infty$ limit of O(n) models with antiperiodic boundary conditions applied across the finite dimension $L$ of the film. We observe that the Casimir amplitude $\Delta_{\rm Casimir}(d|J_\perp,J_\parallel)$ of the anisotropic $d$-dimensional system is related to that one of the isotropic system $\Delta_{\rm Casimir}(d)$ via $\Delta_{\rm Casimir}(d|J_\perp,J_\parallel)=(J_\perp/J_\parallel)^{(d-1)/2} \Delta_{\rm Casimir}(d)$. For $d=3$ we find the exact Casimir amplitude $\Delta_{\rm Casimir}= [ {\rm Cl}_2 (\pi/3)/3-\zeta (3)/(6 \pi)](J_\perp/J_\parallel)$, as well as the exact scaling functions of the Casimir force and of the helicity modulus $\Upsilon(T,L)$. We obtain that $\beta_c\Upsilon(T_c,L)=(2/\pi^{2}) [{\rm Cl}_2(\pi/3)/3+7\zeta(3)/(30\pi)] (J_\perp/J_\parallel)L^{-1}$, where $T_c$ is the critical temperature of the bulk system. We find that the effect of the helicity is thus strong that the Casimir force is repulsive in the whole temperature region.Comment: 15 pages, 3 figure

Interplay of critical Casimir and dispersion forces

Using general scaling arguments combined with mean-field theory we investigate the critical ($T \simeq T_c$) and off-critical ($T\ne 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) (\beta\Delta \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|)^{\nu/\beta}$, with $\nu$ 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.Comment: 26 pages, 14 figure