Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0607250 · Physical Review B

Hidden zero-temperature bicritical point in the two-dimensional anisotropic Heisenberg model: Monte Carlo simulations and proper finite-size scaling

Chenggang Zhou, David Paul Landau, T. C. Schulthess

Published Aug 11, 2006• 16 citations
arXivHTMLPDF

Abstract

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+ϵ$ dimensions. We found that the long-length-scale physics of the Monte Carlo simulations is indeed captured by the anisotropic nonlinear $\ensuremath{\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.

DOI
10.1103/physrevb.74.064407

Related papers

Bicritical universality of the anisotropic Heisenberg model in a crystal field

Rodrigo Teixeira Santos Freire, J. A. Plascak

The bicritical properties of the three-dimensional classical anisotropic Heisenberg model in a crystal field are investigated through extensive Monte Carlo simulations on a simple cubic lattice, using Metropolis and Wolff algorithms. Field-mixing and multidimensional histogram techniques were employed in order to compute the probability distribution function of the extensive conjugate variables of interest and, using finite-size scaling analysis, the first-order transition line of the model was precisely located. The fourth-order cumulant of the order parameter was then calculated along this line and the bicritical point located with good precision from the cumulant crossings. The bicritical properties of this point were further investigated through the measurement of the universal probability distribution function of the order parameter. The results lead us to conclude that the studied bicritical point belongs in fact to the three-dimensional Heisenberg universality class.

Physical Review E•Mar 27, 2015•4 citations

Puzzle of bicriticality in the XXZ antiferromagnet

Amnon Aharony A, Ora Entin-wohlman A

Renormalization-group theory predicts that the XXZ antiferromagnet in a magnetic field along the easy $Z$ axis has asymptotically either a tetracritical phase diagram or a triple point in the field-temperature plane. Neither experiments nor Monte Carlo simulations procure such phase diagrams. Instead, they find a bicritical phase diagram. Here, this discrepancy is resolved: After generalizing a ubiquitous condition identifying the tetracritical point, we employ different renormalization-group recursion relations near the isotropic fixed point, exploiting group-theoretical considerations and using accurate exponents at three dimensions. These show that the results from experiments and simulations can only be understood if their trajectories flow towards the fluctuation-driven first-order transition (and the associated triple point), but reach this limit only for prohibitively large system sizes or correlation lengths. In the crossover region one expects a bicritical phase diagram, as indeed is observed. A similar scenario may explain puzzling discrepancies between simulations and renormalization-group predictions for a variety of other phase diagrams with competing order parameters.

arXiv:2203.01747•Physical review. B./Physical review. B•Sep 21, 2022•5 citations

High-resolution Monte Carlo study of the multicritical point in the three-dimensional XXZ Heisenberg antiferromagnet

Siyan Hu, Shan-Ho Tsai, D. P. Landau

We use Monte Carlo simulations to study the XXZ Heisenberg antiferromagnet in a field in order to clearly determine the nature of the multicritical point. We use a hybrid sampling method with Metropolis and Wolff-cluster algorithms, along with histogram reweighting techniques. Staggered magnetization susceptibilities, Binder cumulants, and finite-size scaling are considered in an effort to detect a possible biconical phase. An analysis of the probability distribution of the magnetization allowed us to conclude that the multicritical point is bicritical and it is in the three-dimensional Heisenberg universality class.

Physical Review E•Mar 17, 2014•6 citations

Bicritical or tetracritical: the 3D anisotropic Heisenberg antiferromagnet

Shan-Ho Tsai, Siyan Hu, David Paul Landau

The classical uniaxially anisotropic Heisenberg antiferromagnet on the simple cubic lattice, in the presence of an external magnetic field, is believed to have a multicritical point; however, there has been controversy whether it is a bicritical or a tetracritical point. We perform Monte Carlo simulations of this model and analyze the components of the staggered magnetization, the susceptibilities and the probability distribution of the magnetization to conclude that the multicritical point is bicritical and it is in the three-dimensional Heisenberg universality class.

Journal of Physics Conference Series•Mar 11, 2014•3 citations

Asymptotic freedom and finite-size scaling of the two-dimensional classical Heisenberg model

Dingyun Yao, Chao Zhang, Z. Y. Xie +2 more

The classical Heisenberg model is one of the most fundamental models in statistical and condensed matter physics. Extensive theoretical and numerical studies suggest that, in two dimensions, this model does not exhibit a finite-temperature phase transition but instead manifests asymptotic freedom. However, some research has also proposed the possibility of a Berezinskii-Kosterlitz-Thouless (BKT) phase transition over the years. In this study, we revisit the classical two-dimensional (2D) Heisenberg model through large-scale simulations with linear system sizes up to $L=16384$. Our Monte Carlo data, without any extrapolation, clearly reveal an exponential divergence of the correlation length $\ensuremath{\xi}$ as a function of the inverse temperature $\ensuremath{\beta}$, a hallmark of asymptotic freedom. Moreover, extrapolating $\ensuremath{\xi}$ to the thermodynamic limit in the low-temperature regime achieves close agreement with the three-loop perturbative calculations. We further propose a finite-size scaling (FSS) ansatz for $\ensuremath{\xi}$, demonstrating that the pseudocritical point ${\ensuremath{\beta}}_{L}$ diverges logarithmically with $L$. The thermodynamic and finite-size scaling behaviors of the magnetic susceptibility $\ensuremath{\chi}$ are also investigated and corroborate the predictions of asymptotic freedom. Our paper provides solid evidence for asymptotic freedom in the 2D Heisenberg model and advances understanding of finite-size scaling in such systems.

arXiv:2503.24226•Physical review. B./Physical review. B•Jun 2, 2025•4 citations

Monte Carlo renormalization-group analysis of the classical Heisenberg model in two dimensions

Stephen H. Shenker, Jan Tobochnik

Monte Carlo renormalization-group methods were applied to the classical three-component Heisenberg model on a two-dimensional lattice. Expectation values of local correlations of spins and various sized block spins were computed using traditional Monte Carlo methods. By matching quantities at different length scales generated by different Monte Carlo Hamiltonians we directly determined the renormalization of the nearest-neighbor coupling. Using these data and the results of high-temperature expansions and low-temperature renormalization-group calculations we have determined that this model does not have a phase transition. We have also obtained the amplitude for the low-temperature divergence of the susceptibility and the correlation length.

Physical review. B, Condensed matter•Nov 1, 1980•259 citations
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+ϵ2+ϵ2+ϵ dimensions. We found that the long-length-scale physics of the Monte Carlo simulations is indeed captured by the anisotropic nonlinear \ensuremathσ\ensuremath{\sigma}\ensuremathσ model. Our Monte Carlo data and analysis confirm that the bicritical point in two dimensions is Heisenberg-like and occurs at T=0T=0T=0; therefore the uncertainty in the phase diagram of this model is removed.
The bicritical properties of the three-dimensional classical anisotropic Heisenberg model in a crystal field are investigated through extensive Monte Carlo simulations on a simple cubic lattice, using Metropolis and Wolff algorithms. Field-mixing and multidimensional histogram techniques were employed in order to compute the probability distribution function of the extensive conjugate variables of interest and, using finite-size scaling analysis, the first-order transition line of the model was precisely located. The fourth-order cumulant of the order parameter was then calculated along this line and the bicritical point located with good precision from the cumulant crossings. The bicritical properties of this point were further investigated through the measurement of the universal probability distribution function of the order parameter. The results lead us to conclude that the studied bicritical point belongs in fact to the three-dimensional Heisenberg universality class.
Renormalization-group theory predicts that the XXZ antiferromagnet in a magnetic field along the easy ZZZ axis has asymptotically either a tetracritical phase diagram or a triple point in the field-temperature plane. Neither experiments nor Monte Carlo simulations procure such phase diagrams. Instead, they find a bicritical phase diagram. Here, this discrepancy is resolved: After generalizing a ubiquitous condition identifying the tetracritical point, we employ different renormalization-group recursion relations near the isotropic fixed point, exploiting group-theoretical considerations and using accurate exponents at three dimensions. These show that the results from experiments and simulations can only be understood if their trajectories flow towards the fluctuation-driven first-order transition (and the associated triple point), but reach this limit only for prohibitively large system sizes or correlation lengths. In the crossover region one expects a bicritical phase diagram, as indeed is observed. A similar scenario may explain puzzling discrepancies between simulations and renormalization-group predictions for a variety of other phase diagrams with competing order parameters.
We use Monte Carlo simulations to study the XXZ Heisenberg antiferromagnet in a field in order to clearly determine the nature of the multicritical point. We use a hybrid sampling method with Metropolis and Wolff-cluster algorithms, along with histogram reweighting techniques. Staggered magnetization susceptibilities, Binder cumulants, and finite-size scaling are considered in an effort to detect a possible biconical phase. An analysis of the probability distribution of the magnetization allowed us to conclude that the multicritical point is bicritical and it is in the three-dimensional Heisenberg universality class.
The classical uniaxially anisotropic Heisenberg antiferromagnet on the simple cubic lattice, in the presence of an external magnetic field, is believed to have a multicritical point; however, there has been controversy whether it is a bicritical or a tetracritical point. We perform Monte Carlo simulations of this model and analyze the components of the staggered magnetization, the susceptibilities and the probability distribution of the magnetization to conclude that the multicritical point is bicritical and it is in the three-dimensional Heisenberg universality class.
The classical Heisenberg model is one of the most fundamental models in statistical and condensed matter physics. Extensive theoretical and numerical studies suggest that, in two dimensions, this model does not exhibit a finite-temperature phase transition but instead manifests asymptotic freedom. However, some research has also proposed the possibility of a Berezinskii-Kosterlitz-Thouless (BKT) phase transition over the years. In this study, we revisit the classical two-dimensional (2D) Heisenberg model through large-scale simulations with linear system sizes up to L=16384L=16384L=16384. Our Monte Carlo data, without any extrapolation, clearly reveal an exponential divergence of the correlation length $\ensuremath{\xi}$ as a function of the inverse temperature $\ensuremath{\beta}$, a hallmark of asymptotic freedom. Moreover, extrapolating $\ensuremath{\xi}$ to the thermodynamic limit in the low-temperature regime achieves close agreement with the three-loop perturbative calculations. We further propose a finite-size scaling (FSS) ansatz for $\ensuremath{\xi}$, demonstrating that the pseudocritical point ${\ensuremath{\beta}}_{L}$ diverges logarithmically with LLL. The thermodynamic and finite-size scaling behaviors of the magnetic susceptibility $\ensuremath{\chi}$ are also investigated and corroborate the predictions of asymptotic freedom. Our paper provides solid evidence for asymptotic freedom in the 2D Heisenberg model and advances understanding of finite-size scaling in such systems.
Monte Carlo renormalization-group methods were applied to the classical three-component Heisenberg model on a two-dimensional lattice. Expectation values of local correlations of spins and various sized block spins were computed using traditional Monte Carlo methods. By matching quantities at different length scales generated by different Monte Carlo Hamiltonians we directly determined the renormalization of the nearest-neighbor coupling. Using these data and the results of high-temperature expansions and low-temperature renormalization-group calculations we have determined that this model does not have a phase transition. We have also obtained the amplitude for the low-temperature divergence of the susceptibility and the correlation length.