Hyperuniversality and the renormalization group for finite systems

Statistical field theory with constraints: Application to critical Casimir forces in the canonical ensemble

The effect of imposing a constraint on a fluctuating scalar order parameter field in a system of finite volume is studied within statistical field theory. The canonical ensemble, corresponding to a fixed total integrated order parameter (e.g., the total number of particles), is obtained as a special case of the theory. A perturbative expansion is developed which allows one to systematically determine the constraint-induced finite-volume corrections to the free energy and to correlation functions. In particular, we focus on the Landau-Ginzburg model in a film geometry (i.e., in a rectangular parallelepiped with a small aspect ratio) with periodic, Dirichlet, or Neumann boundary conditions in the transverse direction and periodic boundary conditions in the remaining, lateral directions. Within the expansion in terms of ε=4-d, where d is the spatial dimension of the bulk, the finite-size contribution to the free energy of the confined system and the associated critical Casimir force are calculated to leading order in ε and are compared to the corresponding expressions for an unconstrained (grand canonical) system. The constraint restricts the fluctuations within the system and it accordingly modifies the residual finite-size free energy. The resulting critical Casimir force is shown to depend on whether it is defined by assuming a fixed transverse area or a fixed total volume. In the former case, the constraint is typically found to significantly enhance the attractive character of the force as compared to the grand canonical case. In contrast to the grand canonical Casimir force, which, for supercritical temperatures, vanishes in the limit of thick films, in the canonical case with fixed transverse area the critical Casimir force attains for thick films a negative value for all boundary conditions studied here. Typically, the dependence of the critical Casimir force both on the temperaturelike and on the fieldlike scaling variables is different in the two ensembles.

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.

Asymmetric fluid criticality. II. Finite-size scaling for simulations

The vapor-liquid critical behavior of intrinsically asymmetric fluids is studied in finite systems of linear dimensions L focusing on periodic boundary conditions, as appropriate for simulations. The recently propounded "complete" thermodynamic (L--> infinity) scaling theory incorporating pressure mixing in the scaling fields as well as corrections to scaling [Phys. Rev. E 67, 061506 (2003)] is extended to finite L, initially in a grand canonical representation. The theory allows for a Yang-Yang anomaly in which, when L--> infinity, the second temperature derivative (d2musigma/dT2) of the chemical potential along the phase boundary musigmaT diverges when T-->Tc-. The finite-size behavior of various special critical loci in the temperature-density or (T,rho) plane, in particular, the k-inflection susceptibility loci and the Q-maximal loci--derived from QL(T,<rho>L) is identical with <m2>2L/<m4>L where m is identical with rho-<rho>L--is carefully elucidated and shown to be of value in estimating Tc and rhoc. Concrete illustrations are presented for the hard-core square-well fluid and for the restricted primitive model electrolyte including an estimate of the correlation exponent nu that confirms Ising-type character. The treatment is extended to the canonical representation where further complications appear.

Exact amplitude ratio and finite-size corrections for the MxN square lattice Ising model

Let f, U, and C represent, respectively, the free energy, the internal energy, and the specific heat of the critical Ising model on the MxN square lattice with periodic boundary conditions, and f(infinity) represents f for fixed M/N and N-->infinity. We find that f, U, and C can be written as N(f-f(infinity))= summation operator(infinity)(i=1)f(2i-1)/N(2i-1), U=-square root of [2]+ summation operator(infinity)(i=1)u(2i-1)/N(2i-1), and C=8 ln N/pi+ summation operator(infinity)(i=0)c(i)/N(i), i.e., Nf and U are odd functions of N(-1). We also find that u(2i-1)/c(2i-1)=1/square root of [2] and u(2i)/c(2i)=0 for 1 < or = i <infinity and obtain closed form expressions for f, U, and C up to orders 1/N(5), 1/N(5), and 1/N(3), respectively, which implies an analytic equation for c(5).