Comment on "Critical properties of highly frustrated pyrochlore antiferromagnets"

We argue that the analysis of Reimers {\it et al.} [ Phys. Rev. B {\bf 45}, 7295 (1992)] of their Monte Carlo data on the Heisenberg pyrochlore antiferromagnet, which suggests a new universality class, is not conclusive. By re-analysis of their data, we demonstrate asymptotic volume dependence in some thermodynamic quantities, which suggests the possibility that the transition may be first order.Comment: 5 pages (RevTex 3.0), 3 figures available upon request, CRPS-93-0

Critical exponents of a three dimensional O(4) spin model

By Monte Carlo simulation we study the critical exponents governing the transition of the three-dimensional classical O(4) Heisenberg model, which is considered to be in the same universality class as the finite-temperature QCD with massless two flavors. We use the single cluster algorithm and the histogram reweighting technique to obtain observables at the critical temperature. After estimating an accurate value of the inverse critical temperature \Kc=0.9360(1), we make non-perturbative estimates for various critical exponents by finite-size scaling analysis. They are in excellent agreement with those obtained with the $4-\epsilon$ expansion method with errors reduced to about halves of them.Comment: 25 pages with 8 PS figures, LaTeX, UTHEP-28

Magnetic Phase Diagram of the Ferromagnetically Stacked Triangular XY Antiferromagnet: A Finite-Size Scaling Study

Histogram Monte-Carlo simulation results are presented for the magnetic-field -- temperature phase diagram of the XY model on a stacked triangular lattice with antiferromagnetic intraplane and ferromagnetic interplane interactions. Finite-size scaling results at the various transition boundaries are consistent with expectations based on symmetry arguments. Although a molecular-field treatment of the Hamiltonian fails to reproduce the correct structure for the phase diagram, it is demonstrated that a phenomenological Landau-type free-energy model contains all the esstential features. These results serve to complement and extend our earlier work [Phys. Rev. B {\bf 48}, 3840 (1993)].Comment: 5 pages (RevTex 3.0), 6 figures available upon request, CRPS 93-

Tricritical behavior of the frustrated XY antiferromagnet

Extensive histogram Monte-Carlo simulations of the XY antiferromagnet on a stacked triangular lattice reveal exponent estimates which strongly favor a scenario of mean-field tricritical behavior for the spin-order transition. The corresponding chiral-order transition occurs at the same temperature but appears to be decoupled from the spin-order. These results are relevant to a wide class of frustrated systems with planar-type order and serve to resolve a long-standing controversy regarding their criticality.Comment: J1K 2R1 4 pages (RevTex 3.0), 4 figures available upon request, Report# CRPS-94-0

Scaling Properties of Antiferromagnetic Transition in Coupled Spin Ladder Systems Doped with Nonmagnetic Impurities

We study effects of interladder coupling on critical magnetic properties of spin ladder systems doped with small concentrations of nonmagnetic impurities, using the scaling theory together with quantum Monte Carlo (QMC) calculations. Scaling properties in a wide region in the parameter space of the impurity concentration x and the interladder coupling are governed by the quantum critical point (QCP) of the undoped system for the transition between antiferromagnetically ordered and spin-gapped phases. This multi-dimensional and strong-coupling region has characteristic power-law dependences on x for magnetic properties such as the N\'eel temperature. The relevance of this criticality for understanding experimental results of ladder compounds is stressed.Comment: 4 pages LaTeX including 3 PS figure

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

Order-Disorder Transition in a Two-Layer Quantum Antiferromagnet

We have studied the antiferromagnetic order -- disorder transition occurring at $T=0$ in a 2-layer quantum Heisenberg antiferromagnet as the inter-plane coupling is increased. Quantum Monte Carlo results for the staggered structure factor in combination with finite-size scaling theory give the critical ratio $J_c = 2.51 \pm 0.02$ between the inter-plane and in-plane coupling constants. The critical behavior is consistent with the 3D classical Heisenberg universality class. Results for the uniform magnetic susceptibility and the correlation length at finite temperature are compared with recent predictions for the 2+1-dimensional nonlinear $\sigma$-model. The susceptibility is found to exhibit quantum critical behavior at temperatures significantly higher than the correlation length.Comment: 11 pages (5 postscript figures available upon request), Revtex 3.

Monte Carlo Study of Cluster-Diameter Distribution: A New Observable to Estimate Correlation Lengths

We report numerical simulations of two-dimensional $q$-state Potts models with emphasis on a new quantity for the computation of spatial correlation lengths. This quantity is the cluster-diameter distribution function $G_{diam}(x)$, which measures the distribution of the diameter of stochastically defined cluster. Theoretically it is predicted to fall off exponentially for large diameter $x$, $G_{diam} \propto \exp(-x/\xi)$, where $\xi$ is the correlation length as usually defined through the large-distance behavior of two-point correlation functions. The results of our extensive Monte Carlo study in the disordered phase of the models with $q=10$, 15, and $20$ on large square lattices of size $300 \times 300$, $120 \times 120$, and $80 \times 80$, respectively, clearly confirm the theoretically predicted behavior. Moreover, using this observable we are able to verify an exact formula for the correlation length $\xi_d(\beta_t)$ in the disordered phase at the first-order transition point $\beta_t$ with an accuracy of about $1$ for all considered values of $q$. This is a considerable improvement over estimates derived from the large-distance behavior of standard (projected) two-point correlation functions, which are also discussed for comparison.Comment: 20 pages, LaTeX + 13 postscript figures. See also http://www.cond-mat.physik.uni-mainz.de/~janke/doc/home_janke.htm

Critical Exponents of the Classical 3D Heisenberg Model: A Single-Cluster Monte Carlo Study

We have simulated the three-dimensional Heisenberg model on simple cubic lattices, using the single-cluster Monte Carlo update algorithm. The expected pronounced reduction of critical slowing down at the phase transition is verified. This allows simulations on significantly larger lattices than in previous studies and consequently a better control over systematic errors. In one set of simulations we employ the usual finite-size scaling methods to compute the critical exponents $\nu,\alpha,\beta,\gamma, \eta$ from a few measurements in the vicinity of the critical point, making extensive use of histogram reweighting and optimization techniques. In another set of simulations we report measurements of improved estimators for the spatial correlation length and the susceptibility in the high-temperature phase, obtained on lattices with up to $100^3$ spins. This enables us to compute independent estimates of $\nu$ and $\gamma$ from power-law fits of their critical divergencies.Comment: 33 pages, 12 figures (not included, available on request). Preprint FUB-HEP 19/92, HLRZ 77/92, September 199