Testing fixed points in the 2D O(3) non-linear sigma model

Using high statistic numerical results we investigate the properties of the O(3) non-linear 2D sigma-model. Our main concern is the detection of an hypothetical Kosterlitz-Thouless-like (KT) phase transition which would contradict the asymptotic freedom scenario. Our results do not support such a KT-like phase transition.Comment: Latex, 7 pgs, 4 eps-figures. Added more analysis on the KT-transition. 4-loop beta function contains corrections from D.-S.Shin (hep-lat/9810025). In a note-added we comment on the consequences of these corrections on our previous reference [16

Critical behavior of the planar magnet model in three dimensions

We use a hybrid Monte Carlo algorithm in which a single-cluster update is combined with the over-relaxation and Metropolis spin re-orientation algorithm. Periodic boundary conditions were applied in all directions. We have calculated the fourth-order cumulant in finite size lattices using the single-histogram re-weighting method. Using finite-size scaling theory, we obtained the critical temperature which is very different from that of the usual XY model. At the critical temperature, we calculated the susceptibility and the magnetization on lattices of size up to $42^3$. Using finite-size scaling theory we accurately determine the critical exponents of the model and find that $\nu$=0.670(7), $\gamma/\nu$=1.9696(37), and $\beta/\nu$=0.515(2). Thus, we conclude that the model belongs to the same universality class with the XY model, as expected.Comment: 11 pages, 5 figure

O(N) and RP^{N-1} Models in Two Dimensions

I provide evidence that the 2D $RP^{N-1}$ model for $N \ge 3$ is equivalent to the $O(N)$-invariant non-linear $\sigma$-model in the continuum limit. To this end, I mainly study particular versions of the models, to be called constraint models. I prove that the constraint $RP^{N-1}$ and $O(N)$ models are equivalent for sufficiently weak coupling. Numerical results for their step-scaling function of the running coupling $\bar{g}^2= m(L) L$ are presented. The data confirm that the constraint $O(N)$ model is in the samei universality class as the $O(N)$ model with standard action. I show that the differences in the finite size scaling curves of $RP^{N-1}$i and $O(N)$ models observed by Caracciolo et al. can be explained as a boundary effect. It is concluded, in contrast to Caracciolo et al., that $RP^{N-1}$ and $O(N)$ models share a unique universality class.Comment: 14 pages (latex) + 1 figure (Postscript) ,uuencode

Physical tests for Random Numbers in Simulations

We propose three physical tests to measure correlations in random numbers used in Monte Carlo simulations. The first test uses autocorrelation times of certain physical quantities when the Ising model is simulated with the Wolff algorithm. The second test is based on random walks, and the third on blocks of n successive numbers. We apply the tests to show that recent errors in high precision simulations using generalized feedback shift register algorithms are due to short range correlations in random number sequences. We also determine the length of these correlations.Comment: 16 pages, Post Script file, HU-TFT-94-

Self-Averaging, Distribution of Pseudo-Critical Temperatures and Finite Size Scaling in Critical Disordered Systems

The distributions $P(X)$ of singular thermodynamic quantities in an ensemble of quenched random samples of linear size $l$ at the critical point $T_c$ are studied by Monte Carlo in two models. Our results confirm predictions of Aharony and Harris based on Renormalization group considerations. For an Ashkin-Teller model with strong but irrelevant bond randomness we find that the relative squared width, $R_X$, of $P(X)$ is weakly self averaging. $R_X\sim l^{\alpha/\nu}$, where $\alpha$ is the specific heat exponent and $\nu$ is the correlation length exponent of the pure model fixed point governing the transition. For the site dilute Ising model on a cubic lattice, known to be governed by a random fixed point, we find that $R_X$ tends to a universal constant independent of the amount of dilution (no self averaging). However this constant is different for canonical and grand canonical disorder. We study the distribution of the pseudo-critical temperatures $T_c(i,l)$ of the ensemble defined as the temperatures of the maximum susceptibility of each sample. We find that its variance scales as $(\delta T_c(l))^2 \sim l^{-2/\nu}$ and NOT as $\sim l^{-d}. We find that$R_\chi$is reduced by a factor of$\sim 70$with respect to$R_\chi (T_c)$by measuring$\chi$of each sample at$T_c(i,l)$. We analyze correlations between the magnetization at criticality$m_i(T_c,l)$and the pseudo-critical temperature$T_c(i,l)$in terms of a sample independent finite size scaling function of a sample dependent reduced temperature$(T-T_c(i,l))/T_c$. This function is found to be universal and to behave similarly to pure systems.Comment: 31 pages, 17 figures, submitted to Phys. Rev.

Low Temperature Static and Dynamic Behavior of the Two-Dimensional Easy-Axis Heisenberg Model

We apply the self-consistent harmonic approximation (SCHA) to study static and dynamic properties of the two-dimensional classical Heisenberg model with easy-axis anisotropy. The static properties obtained are magnetization and spin wave energy as functions of temperature, and the critical temperature as a function of the easy-axis anisotropy. We also calculate the dynamic correlation functions using the SCHA renormalized spin wave energy. Our analytical results, for both static properties and dynamic correlation functions, are compared to numerical simulation data combining cluster-Monte Carlo algorithms and Spin Dynamics. The comparison allows us to conclude that far below the transition temperature, where the SCHA is valid, spin waves are responsible for all relevant features observed in the numerical simulation data; topological excitations do not seem to contribute appreciably. For temperatures closer to the transition temperature, there are differences between the dynamic correlation functions from SCHA theory and Spin Dynamics; these may be due to the presence of domain walls and solitons.Comment: 12 pages, 14 figure

Ising spins coupled to a four-dimensional discrete Regge skeleton

Regge calculus is a powerful method to approximate a continuous manifold by a simplicial lattice, keeping the connectivities of the underlying lattice fixed and taking the edge lengths as degrees of freedom. The discrete Regge model employed in this work limits the choice of the link lengths to a finite number. To get more precise insight into the behavior of the four-dimensional discrete Regge model, we coupled spins to the fluctuating manifolds. We examined the phase transition of the spin system and the associated critical exponents. The results are obtained from finite-size scaling analyses of Monte Carlo simulations. We find consistency with the mean-field theory of the Ising model on a static four-dimensional lattice.Comment: 19 pages, 7 figure

Cluster Monte Carlo Simulations of the Nematic--Isotropic Transition

We report the results of simulations of the Lebwohl-Lasher model of the nematic-isotropic transition using a new cluster Monte Carlo algorithm. The algorithm is a modification of the Wolff algorithm for spin systems, and greatly reduces critical slowing down. We calculate the free energy in the neighborhood of the transition for systems up to linear size 70. We find a double well structure with a barrier that grows with increasing system size, obeying finite size scaling for systems of size greater than 35. We thus obtain an estimate of the value of the transition temperature in the thermodynamic limit.Comment: 4 figure