No Arabic abstract
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 $ntoinfty$ 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_cUpsilon(T_c,L)=(2/pi^{2}) [{rm Cl}_2(pi/3)/3+7zeta(3)/(30pi)] (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.
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 proceeding comment raises a few points concerning our paper Dantchev textit{et al.}, Phys. Rev. E. {bf 89}, 042116 (2014). In this reply we stress that while Refs. Diehl textit{et al.} EPL {bf 100}, 10004 (2012) and Phys. Rev. E. {bf 89}, 062123 (2014) use three different models to study the the Casimir force for the $O(n rightarrow infty)$ model with free boundary conditions we study a single model over the entire range of temperatures, from above the bulk critical temperature, $T_c$, to absolute temperatures down to $T=0$. The use of a single model renders more transparent the crossover from effects dominated by critical fluctuations in the vicinity of the bulk transition temperature to effects controlled by Goldstone modes at low temperatures. Contrary to the assertion in the comment, we make no claim for the superiority of our model over any of those considered by Diehl textit{et al}. We also present additional evidence supporting our conclusion in Dantchev textit{et al.}, Phys. Rev. E. {bf 89}, 042116 (2014) that the temperature range in which our low-temperature analytical expansion for the Casimir force increases as $L$ grows and remains accurate for values of the ratio $T/T_c$ that become closer and closer to unity, while $T$ remains well outside of the critical region.
We consider two disordered lattice models on the square lattice: on the medial lattice the random field Ising model at T=0 and on the direct lattice the random bond Potts model in the large-q limit at its transition point. The interface properties of the two models are known to be related by a mapping which is valid in the continuum approximation. Here we consider finite random samples with the same form of disorder for both models and calculate the respective equilibrium states exactly by combinatorial optimization algorithms. We study the evolution of the interfaces with the strength of disorder and analyse and compare the interfaces of the two models in finite lattices.
Monte Carlo (MC) analysis of the Goldstone mode singularities for the transverse and the longitudinal correlation functions, behaving as G_{perp}(k) simeq ak^{-lambda_{perp}} and G_{parallel}(k) simeq bk^{-lambda_{parallel}} in the ordered phase at k -> 0, is performed in the three-dimensional O(n) models with n=2, 4, 10. Our aim is to test some challenging theoretical predictions, according to which the exponents lambda_{perp} and lambda_{parallel} are non-trivial (3/2<lambda_{perp}<2 and 0<lambda_{parallel}<1 in three dimensions) and the ratio bM^2/a^2 (where M is a spontaneous magnetization) is universal. The trivial standard-theoretical values are lambda_{perp}=2 and lambda_{parallel}=1. Our earlier MC analysis gives lambda_{perp}=1.955 pm 0.020 and lambda_{parallel} about 0.9 for the O(4) model. A recent MC estimation of lambda_{parallel}, assuming corrections to scaling of the standard theory, yields lambda_{parallel} = 0.69 pm 0.10 for the O(2) model. Currently, we have performed a similar MC estimation for the O(10) model, yielding lambda_{perp} = 1.9723(90). We have observed that the plot of the effective transverse exponent for the O(4) model is systematically shifted down with respect to the same plot for the O(10) model by Delta lambda_{perp} = 0.0121(52). It is consistent with the idea that 2-lambda_{perp} decreases for large $n$ and tends to zero at n -> infty. We have also verified and confirmed the expected universality of bM^2/a^2 for the O(4) model, where simulations at two different temperatures (couplings) have been performed.
We present exact results on the behavior of the thermodynamic Casimir force and the excess free energy in the framework of the $d$-dimensional spherical model with a power law long-range interaction decaying at large distances $r$ as $r^{-d-sigma}$, where $sigma<d<2sigma$ and $0<sigmaleq2$. For a film geometry and under periodic boundary conditions we consider the behavior of these quantities near the bulk critical temperature $T_c$, as well as for $T>T_c$ and $T<T_c$. The universal finite-size scaling function governing the behavior of the force in the critical region is derived and its asymptotics are investigated. While in the critical and under critical region the force is of the order of $L^{-d}$, for $T>T_c$ it decays as $L^{-d-sigma}$, where $L$ is the thickness of the film. We consider both the case of a finite system that has no phase transition of its own, when $d-1<sigma$, as well as the case with $d-1>sigma$, when one observes a dimensional crossover from $d$ to a $d-1$ dimensional critical behavior. The behavior of the force along the phase coexistence line for a magnetic field H=0 and $T<T_c$ is also derived. We have proven analytically that the excess free energy is always negative and monotonically increasing function of $T$ and $H$. For the Casimir force we have demonstrated that for any $sigma ge 1$ it is everywhere negative, i.e. an attraction between the surfaces bounding the system is to be observed. At $T=T_c$ the force is an increasing function of $T$ for $sigma>1$ and a decreasing one for $sigma<1$. For any $d$ and $sigma$ the minimum of the force at $T=T_c$ is always achieved at some $H e 0$.