No Arabic abstract
Corrections to scaling in the 3D Ising model are studied based on non-perturbative analytical arguments and Monte Carlo (MC) simulation data for different lattice sizes L. Analytical arguments show the existence of corrections with the exponent (gamma-1)/nu (approximately 0.38), the leading correction-to-scaling exponent being omega =< (gamma-1)/nu. A numerical estimation of omega from the susceptibility data within 40 =< L =< 2048 yields omega=0.25(33). It is consistent with the statement omega =< (gamma-1)/nu, as well as with the value omega = 1/8 of the GFD theory. We reconsider the MC estimation of omega from smaller lattice sizes to show that it does not lead to conclusive results, since the obtained values of omega depend on the particular method chosen. In particular, estimates ranging from omega =1.274(72) to omega=0.18(37) are obtained by four different finite-size scaling methods, using MC data for thermodynamic average quantities, as well as for partition function zeros. We discuss the influence of omega on the estimation of exponents eta and nu.
Corrections to scaling in the two-dimensional scalar phi^4 model are studied based on non-perturbative analytical arguments and Monte Carlo (MC) simulation data for different lattice sizes L (from 4 to 1536) and different values of the phi^4 coupling constant lambda, i.~e., lambda = 0.1, 1, 10. According to our analysis, amplitudes of the nontrivial correction terms with the correction-to-scaling exponents omega_l < 1 become small when approaching the Ising limit (lambda --> infinity), but such corrections generally exist in the 2D phi^4 model. Analytical arguments show the existence of corrections with the exponent 3/4. The numerical analysis suggests that there exist also corrections with the exponent 1/2 and, very likely, also corrections with the exponent about 1/4, which are detectable at lambda = 0.1. The numerical tests clearly show that the structure of corrections to scaling in the 2D phi^4 model differs from the usually expected one in the 2D Ising model.
By considering the appropriate finite-size effect, we explain the connection between Monte Carlo simulations of two-dimensional anisotropic Heisenberg antiferromagnet in a field and the early renormalization group calculation for the bicritical point in $2+epsilon$ dimensions. We found that the long length scale physics of the Monte Carlo simulations is indeed captured by the anisotropic nonlinear $sigma$ model. Our Monte Carlo data and analysis confirm that the bicritical point in two dimensions is Heisenberg-like and occurs at T=0, therefore the uncertainty in the phase diagram of this model is removed.
Energy eigenvalues and order parameters are calculated by exact diagonalization for the transverse Ising model on square lattices of up to 6x6 sites. Finite-size scaling is used to estimate the critical parameters of the model, confirming universality with the three-dimensional classical Ising model. Critical amplitudes are also estimated for both the energy gap and the ground-state energy.
Critical phenomena and Goldstone mode effects in spin models with O(n) rotational symmetry are considered. Starting with the Goldstone mode singularities in the XY and O(4) models, we briefly review different theoretical concepts as well as state-of-the art Monte Carlo simulation results. They support recent results of the GFD (grouping of Feynman diagrams) theory, stating that these singularities are described by certain nontrivial exponents, which differ from those predicted earlier by perturbative treatments. Furthermore, we present the recent Monte Carlo simulation results of the three-dimensional Ising model for very large lattices with linear sizes up to L=1536. These results are obtained, using a parallel OpenMP implementation of the Wolff single cluster algorithm. The finite-size scaling analysis of the critical exponent eta, assuming the usually accepted correction-to-scaling exponent omega=0.8, shows that eta is likely to be somewhat larger than the value 0.0335 +/- 0.0025 of the perturbative renormalization group (RG) theory. Moreover, we have found that the actual data can be well described by different critical exponents: eta=omega=1/8 and nu=2/3, found within the GFD theory.
Using Monte Carlo simulations, finite-size effects of interfacial properties in the rough phase of the Ising on a cubic lattice with $Ltimes Ltimes R$ sites are studied. In particular, magnetization profiles perpendicular to the flat interface of size L$times$R are studied, with $L$ being considerably larger than $R$, in the (pre)critical temperature range. The resulting $R$-dependences are compared with predictions of the standard capillary-wave theory, in the Gaussian approximation, and with a field theory based on effective string actions, for $L$=$infty$.