A Monte Carlo study of leading order scaling corrections of phi4theory on a three-dimensional lattice

Finite-size scaling at fixed renormalization-group invariant

Finite-size scaling at fixed renormalization-group invariant is a powerful and flexible technique to analyze Monte Carlo data at a critical point. It consists in fixing a given renormalization-group invariant quantity to a given value, thereby trading its statistical fluctuations with those of a parameter driving the transition. One remarkable feature is the observed significant improvement of statistical accuracy of various quantities, as compared to a standard analysis. We review the method, discussing in detail its implementation, the error analysis, and a previously introduced covariance-based optimization. Comprehensive benchmarks on the Ising model in two and three dimensions show large gains in the statistical accuracy, which are due to cross-correlations between observables. As an application, we compute an accurate estimate of the inverse critical temperature of the improved O(2) ϕ^{4} model on a three-dimensional cubic lattice.

Boundary Critical Behavior of the Three-Dimensional Heisenberg Universality Class

We study the boundary critical behavior of the three-dimensional Heisenberg universality class, in the presence of a bidimensional surface. By means of high-precision Monte Carlo simulations of an improved lattice model, where leading bulk scaling corrections are suppressed, we prove the existence of a special phase transition, with unusual exponents, and of an extraordinary phase with logarithmically decaying correlations. These findings contrast with naïve arguments on the bulk-surface phase diagram, and allow us to explain some recent puzzling results on the boundary critical behavior of quantum spin models.

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.

Improvement of Monte Carlo estimates with covariance-optimized finite-size scaling at fixed phenomenological coupling

In the finite-size scaling analysis of Monte Carlo data, instead of computing the observables at fixed Hamiltonian parameters, one may choose to keep a renormalization-group invariant quantity, also called phenomenological coupling, fixed at a given value. Within this scheme of finite-size scaling, we exploit the statistical covariance between the observables in a Monte Carlo simulation in order to reduce the statistical errors of the quantities involved in the computation of the critical exponents. This method is general and does not require additional computational time. This approach is demonstrated in the Ising model in two and three dimensions, where large gain factors in CPU time are obtained.

Error estimation and reduction with cross correlations

Besides the well-known effect of autocorrelations in time series of Monte Carlo simulation data resulting from the underlying Markov process, using the same data pool for computing various estimates entails additional cross correlations. This effect, if not properly taken into account, leads to systematically wrong error estimates for combined quantities. Using a straightforward recipe of data analysis employing the jackknife or similar resampling techniques, such problems can be avoided. In addition, a covariance analysis allows for the formulation of optimal estimators with often significantly reduced variance as compared to more conventional averages.

Equilibriumlike invaded cluster algorithm: critical exponents and dynamical properties

We present a detailed study of the equilibrium-like invaded cluster algorithm, recently proposed as an extension of the invaded cluster algorithm, designed to drive the system to criticality while still preserving the equilibrium ensemble. We perform extensive simulations on two special cases of the Potts model and examine the precision of critical exponents by including the leading corrections. We show that both thermal and magnetic critical exponents can be obtained with high accuracy compared to the best available results. The choice of the auxiliary parameters of the algorithm is discussed in context of dynamical properties. We also discuss the relation to the Li-Sokal bound for the dynamical exponent z.

Cross correlations in scaling analyses of phase transitions

Thermal or finite-size scaling analyses of importance sampling Monte Carlo computer simulations in the vicinity of phase transition points often combine different estimates for the same quantity, such as a critical exponent, with the intent to reduce statistical fluctuations. We point out that the origin of such estimates in the same time series often results in pronounced cross correlations which are usually ignored even in high-precision studies, generically leading to significant underestimation of statistical fluctuations. We suggest to use a simple extension of the conventional analysis taking correlation effects into account, which leads to improved estimators with often substantially reduced statistical fluctuations at almost no extra cost in terms of computation time.

Universal dependence on disorder of two-dimensional randomly diluted and random-bond +/-J Ising models

We consider the two-dimensional randomly site diluted Ising model and the random-bond +/-J Ising model (also called the Edwards-Anderson model), and study their critical behavior at the paramagnetic-ferromagnetic transition. The critical behavior of thermodynamic quantities can be derived from a set of renormalization-group equations, in which disorder is a marginally irrelevant perturbation at the two-dimensional Ising fixed point. We discuss their solutions, focusing in particular on the universality of the logarithmic corrections arising from the presence of disorder. Then, we present a finite-size scaling analysis of high-statistics Monte Carlo simulations. The numerical results confirm the renormalization-group predictions, and in particular the universality of the logarithmic corrections to the Ising behavior due to quenched dilution.

Three-dimensional randomly dilute Ising model: Monte Carlo results

We perform a high-statistics simulation of the three-dimensional randomly dilute Ising model on cubic lattices L3 with L< or =256. We choose a particular value of the density, x=0.8, for which the leading scaling corrections are suppressed. We determine the critical exponents, obtaining nu=0.683(3), eta=0.035(2), beta=0.3535(17), and alpha=-0.049(9), in agreement with previous numerical simulations. We also estimate numerically the fixed-point values of the four-point zero-momentum couplings that are used in field-theoretical fixed-dimension studies. Although these results somewhat differ from those obtained using perturbative field theory, the field-theoretical estimates of the critical exponents do not change significantly if the Monte Carlo result for the fixed point is used. Finally, we determine the six-point zero-momentum couplings, relevant for the small-magnetization expansion of the equation of state, and the invariant amplitude ratio R(+)(xi) that expresses the universality of the free-energy density per correlation volume. We find R(+)(xi)=0.2885(15).

25th-order high-temperature expansion results for three-dimensional Ising-like systems on the simple-cubic lattice

25th-order high-temperature series are computed for a general nearest-neighbor three-dimensional Ising model with arbitrary potential on the simple cubic lattice. In particular, we consider three improved potentials characterized by suppressed leading scaling corrections. Critical exponents are extracted from high-temperature series specialized to improved potentials, obtaining gamma=1.2373(2), nu=0.63012(16), alpha=0.1096(5), eta=0.036 39(15), beta=0.326 53(10), and delta=4.78 93(8). Moreover, biased analyses of the 25th-order series of the standard Ising model provide the estimate Delta=0.52(3) for the exponent associated with the leading scaling corrections. By the same technique, we study the small-magnetization expansion of the Helmholtz free energy. The results are then applied to the construction of parametric representations of the critical equation of state, using a systematic approach based on a global stationarity condition. Accurate estimates of several universal amplitude ratios are also presented.

Critical structure factors of bilinear fields in O(N) vector models

We compute the two-point correlation functions of general quadratic operators in the high-temperature phase of the three-dimensional O(N) vector model by using field-theoretical methods. In particular, we study the small- and large-momentum behavior of the corresponding scaling functions, and give general interpolation formulas based on a dispersive approach. Moreover, we determine the crossover exponent phi(T) associated with the traceless tensorial quadratic field, by computing and analyzing its six-loop perturbative expansion in fixed dimension. We find phi(T)=1.184(12), phi(T)=1.271(21), and phi(T)=1.40(4) for N=2,3,5, respectively.

Monte Carlo renormalization-group analysis of the lattice phi 4 model in D = 3,4

We present a simple, sophisticated method to capture renormalization-group flow in Monte Carlo simulation, which provides important information of critical phenomena. We applied the method to the D = 3,4 lattice phi 4 model and obtained a renormalization flow diagram that well reproduces theoretically predicted behavior of the continuum phi 4 model. We also show that the method can be easily applied to much more complicated models, such as frustrated spin models.

Improved high-temperature expansion and critical equation of state of three-dimensional Ising-like systems

High-temperature series are computed for a generalized three-dimensional Ising model with arbitrary potential. Three specific "improved" potentials (suppressing leading scaling corrections) are selected by Monte Carlo computation. Critical exponents are extracted from high-temperature series specialized to improved potentials, achieving high accuracy; our best estimates are gamma=1.2371(4), nu=0.630 02(23), alpha=0.1099(7), eta=0.0364(4), beta=0.326 48(18). By the same technique, the coefficients of the small-field expansion for the effective potential (Helmholtz free energy) are computed. These results are applied to the construction of parametric representations of the critical equation of state. A systematic approximation scheme, based on a global stationarity condition, is introduced (the lowest-order approximation reproduces the linear parametric model). This scheme is used for an accurate determination of universal ratios of amplitudes. A comparison with other theoretical and experimental determinations of universal quantities is presented.