Order-Parameter Distribution Function of Finite O(n) Symmetric Systems

Exact Critical Casimir Amplitude of Anisotropic Systems from Conformal Field Theory and Self-Similarity of Finite-Size Scaling Functions in d≥2 Dimensions

The exact critical Casimir amplitude is derived for anisotropic systems within the d=2 Ising universality class by combining conformal field theory with anisotropic φ^{4} theory. Explicit results are presented for the general anisotropic scalar φ^{4} model and for the fully anisotropic triangular-lattice Ising model in finite rectangular and infinite strip geometries with periodic boundary conditions. These results demonstrate the validity of multiparameter universality for confined anisotropic systems and the nonuniversality of the critical Casimir amplitude. We find an unexpected complex form of self-similarity of the anisotropy effects near the instability where weak anisotropy breaks down. This can be traced back to the property of modular invariance of isotropic conformal field theory for d=2. More generally, for d>2 we predict the existence of self-similar structures of the finite-size scaling functions of O(n)-symmetric systems with planar anisotropies and periodic boundary conditions both in the critical region for n≥1 as well as in the Goldstone-dominated low-temperature region for n≥2.

Crossover from low-temperature to high-temperature fluctuations: Universal and nonuniversal Casimir forces of isotropic and anisotropic systems

We study the crossover from low-temperature to high-temperature fluctuations including Goldstone-dominated and critical fluctuations in confined isotropic and weakly anisotropic O(n)-symmetric systems on the basis of a finite-size renormalization-group approach at fixed dimension d introduced previously [V. Dohm, Phys. Rev. Lett. 110, 107207 (2013)PRLTAO0031-900710.1103/PhysRevLett.110.107207]. Our theory is formulated within the φ^{4} lattice model in a d-dimensional block geometry with periodic boundary conditions. We calculate the finite-size scaling functions F^{ex} and X of the excess free-energy density and the thermodynamic Casimir force, respectively, for 1≤n≤∞, 2<d<4. Exact results are derived for n→∞. Applications are given for L_{∥}^{d-1}×L slab geometry with an aspect ratio ρ=L/L_{∥}>0 and for film geometry (ρ=0). Good overall agreement is found with Monte Carlo (MC) data for isotropic spin models with n=1,2,3. For ρ=0, the low-temperature limits of F^{ex} and X vanish for n=1, whereas they are finite for n≥2. For ρ>0 and n=1, we find a finite low-temperature limit of F^{ex}, which deviates from that of the Ising model. We attribute this deviation to the nonuniversal difference between the φ^{4} model with continuous variables and the Ising model with discrete variables. For n≥2 and ρ>0, a logarithmic divergence of F^{ex} in the low-temperature limit is predicted, in excellent agreement with MC data. For 2≤n≤∞ and ρ<ρ_{0}=0.8567 the Goldstone modes generate a negative low-temperature Casimir force that vanishes for ρ=ρ_{0} and becomes positive for ρ>ρ_{0}. For anisotropic systems a unified hypothesis of multiparameter universality is introduced for both bulk and confined systems. The dependence of their scaling functions on d(d+1)/2-1 microscopic anisotropy parameters implies a substantial reduction of the predictive power of the theory for anisotropic systems as compared to isotropic systems. An exact representation is derived for the nonuniversal large-distance behavior of the bulk correlation function of anisotropic systems and quantitative predictions are made. The validity of multiparameter universality is proven analytically for the d=2,n=1 universality class. A nonuniversal anisotropy-dependent minimum of the Casimir force scaling function X is found. Both the sign and magnitude of X and the shift of the film critical temperature are affected by the lattice anisotropy.

Pronounced minimum of the thermodynamic Casimir forces of O(n) symmetric film systems: analytic theory

Thermodynamic Casimir forces of film systems in the O(n) universality classes with Dirichlet boundary conditions are studied below bulk criticality. Substantial progress is achieved in resolving the long-standing problem of describing analytically the pronounced minimum of the scaling function observed experimentally in ^{4}He films (n=2) by Garcia and Chan [Phys. Rev. Lett. 83, 1187 (1999)] and in Monte Carlo simulations for the three-dimensional Ising model (n=1) by O. Vasilyev et al. [Europhys. Lett. 80, 60009 (2007)]. Our finite-size renormalization-group approach describes the film systems as the limit of finite-slab systems with vanishing aspect ratio. This yields excellent agreement with the depth and the position of the minimum for n=1 and semiquantitative agreement with the minimum for n=2. Our theory also predicts a pronounced minimum for the n=3 Heisenberg universality class.

Crossover from Goldstone to critical fluctuations: Casimir forces in confined O(n)-symmetric systems

We study the crossover between thermodynamic Casimir forces arising from long-range fluctuations due to Goldstone modes and those arising from critical fluctuations. Both types of forces exist in the low-temperature phase of O(n)-symmetric systems for n>1 in a d-dimensional L(||)(d-1) × L slab geometry with a finite aspect ratio ρ = L/L(||). Our finite-size renormalization-group treatment for periodic boundary conditions describes the entire crossover from the Goldstone regime with a nonvanishing constant tail of the finite-size scaling function far below T(c) up to the region far above T(c) including the critical regime with a minimum of the scaling function slightly below T(c). Our analytic result for ρ << 1 agrees well with Monte Carlo data for the three-dimensional XY model. A quantitative prediction is given for the crossover of systems in the Heisenberg universality class.

Diversity of critical behavior within a universality class

We study spatial anisotropy effects on the bulk and finite-size critical behavior of the O(n) symmetric anisotropic phi;{4} lattice model with periodic boundary conditions in a d -dimensional hypercubic geometry above, at, and below Tc. The absence of two-scale factor universality is discussed for the bulk order-parameter correlation function, the bulk scattering intensity, and for several universal bulk amplitude relations. The anisotropy parameters are observable by scattering experiments at Tc. For the confined system, renormalization-group theory within the minimal subtraction scheme at fixed dimension d for 2<d<4 is employed. In contrast to the epsilon=4-d expansion, the fixed- d finite-size approach keeps the exponential form of the order-parameter distribution function unexpanded. For the case of cubic symmetry and for n=1 , our perturbation approach yields excellent agreement with the Monte Carlo (MC) data for the finite-size amplitude of the free energy of the three-dimensional Ising model at Tc by Mon [Phys. Rev. Lett. 54, 2671 (1985)]. The epsilon expansion result is in less good agreement. Below Tc, a minimum of the scaling function of the excess free energy is found. We predict a measurable dependence of this minimum on the anisotropy parameters. The relative anisotropy effect on the free energy is predicted to be significantly larger than that on the Binder cumulant. Our theory agrees quantitatively with the nonmonotonic dependence of the Binder cumulant on the ferromagnetic next-nearest-neighbor (NNN) coupling of the two-dimensional Ising model found by MC simulations of Selke and Shchur [J. Phys. A 38, L739 (2005)]. Our theory also predicts a nonmonotonic dependence for small values of the antiferromagnetic NNN coupling and the existence of a Lifshitz point at a larger value of this coupling. The nonuniversal anisotropy effects in the finite-size scaling regime are predicted to satisfy a kind of restricted universality. The tails of the large- L behavior at T++Tc violate both finite-size scaling and universality even for isotropic systems as they depend on the bare four-point coupling of the phi4 theory, on the cutoff procedure, and on subleading long-range interactions.

Dynamics of fluids near the consolute critical point: A molecular-dynamics study of the Widom–Rowlinson mixture

Molecular-dynamics simulations are presented for the dynamic behavior of the Widom-Rowlinson mixture [B. Widom, and J. S. Rowlinson, J. Chem. Phys. 52, 1670 (1970)] at its critical point. This model consists of two components where like species do not interact and unlike species interact via a hard-core potential. Critical exponents are obtained from a finite-size scaling analysis. The self-diffusion coefficient shows no anomalous behavior near the critical point. The shear viscosity and thermal conductivity show no divergent behavior for the system sizes considered, although there is a significant critical enhancement. The mutual diffusion coefficient, D(AB), vanishes as D(AB) approximately xi(-1.26 +/- 0.08), where xi is the correlation length. This is different from the renormalization-group (D(AB) approximately xi(-1.065)) mode coupling theory (D(AB) approximately xi(-1)) predictions. The theories and simulations can be reconciled if we assume that logarithmic corrections to scaling are important.

Scaling function for the critical diffusion coefficient of a critical fluid in a finite geometry

The long-wavelength diffusion coefficient of a critical fluid confined between two parallel plates, separated by a distance L, is strongly affected by the finite size. Finite size scaling leads us to expect that the vanishing of the diffusion coefficient as xi(-1) for xi<<L, xi being the correlation length, would crossover to L-1 for xi>>L. We show that this is not strictly true. There is a logarithmic scaling violation. We construct a Kawasaki-like scaling function that connects the thermodynamic regime to the extreme critical (xi>>L) regime.

Anomalous fluctuations in phases with a broken continuous symmetry

It is shown that the Goldstone modes associated with a broken continuous symmetry lead to anomalously large fluctuations of the zero field order parameter at any temperature below T(c). In dimensions 2<d<4, the variance of the extensive spontaneous magnetization scales as L4 with the system size L, independent of the order parameter dynamics. The anomalous scaling is a consequence of the 1/q(4-d) divergence of the longitudinal susceptibility. For ground states in two dimensions with Goldstone modes vanishing linearly with momentum, the dynamical susceptibility contains a singular contribution (q(2)-omega(2)/c(2))(-1/2). The dynamic structure factor thus exhibits a critical continuum above the undamped spin wave pole, which may be detected by neutron scattering in the Néel phase of 2D quantum antiferromagnets.

Probability distribution of the order parameter for the three-dimensional ising-model universality class: A high-precision monte carlo study

We study the probability distribution P(M) of the order parameter (average magnetization) M, for the finite-size systems at the critical point. The systems under consideration are the 3-dimensional Ising model on a simple cubic lattice, and its 3-state generalization known to have remarkably small corrections to scaling. Both models are studied in a cubic box with periodic boundary conditions. The model with reduced corrections to scaling makes it possible to determine P(M) with unprecedented precision. We also obtain a simple, but remarkably accurate, approximate formula describing the universal shape of P(M).