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

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics

Monte Carlo simulation of the ferromagnetic order-disorder transition in a Heisenberg fluid

M. J. P. Nijmeijer, J. J. Weis

Published Jan 1, 1996• 43 citations
Publisher

Abstract

We have simulated a fluid of particles which carry a Heisenberg spin and interact via a Heisenberg-like interaction. The fluid displays a magnetic order-disorder transition, which we studied by varying the temperature at a fixed particle density. We located the transition and estimated critical exponents by means of a finite-size scaling analysis. The exponents differ by a small but significant amount from the ones for the lattice Heisenberg model. \textcopyright{} 1996 The American Physical Society.

DOI
10.1103/physreve.53.591

Related papers

Monte Carlo Simulation of the Ferromagnetic Order-Disorder Transition in a Heisenberg Fluid

M. J. P. Nijmeijer, J. J. Weis

We have performed Monte Carlo simulations of the ferromagnetic order-disorder transition in an off-lattice version of the Heisenberg model. Varying the temperature at a fixed density, we located the transition for three different densities by means of a finite-size scaling analysis. The obtained critical exponents differ from those for the lattice Heisenberg model. Although the limited size of our systems warrants caution, the results question whether the lattice and the off-lattice model are in the same universality class.

Physical Review Letters•Oct 9, 1995•39 citations

Ferromagnetic Phase Transition in a Heisenberg Fluid: Monte Carlo Simulations and Fisher Corrections to Scaling

I. M. Mryglod, I. P. Omelyan, R. Folk

The magnetic phase transition in a Heisenberg fluid is studied by means of the finite size scaling technique. We find that even for larger systems, considered in an ensemble with fixed density, the critical exponents show deviations from the expected lattice values similar to those obtained previously. This puzzle is clarified by proving the importance of the leading correction to the scaling that appears due to Fisher renormalization with the critical exponent equal to the absolute value of the specific heat exponent alpha. The appearance of such new corretions to scaling is a general feature of systems with constraints.

arXiv:cond-mat/0012110•Physical Review Letters•Apr 2, 2001•28 citations

Ferromagnetic order-disorder transition in an Ising fluid

M. J. P. Nijmeijer, Alberto Parola, L. Reatto

The Ising fluid we have investigated is a system of hard spheres which carry an Ising spin. A short range ferromagnetic potential acts between particles. A renormalization group study of the model, based on the hierarchical reference theory of fluids, gives quantitative predictions on the phase diagram and the phase boundaries. The model can show two distinct phase transitions: a magnetic order-disorder transition and a liquid-vapor phase separation which meet at a tricritical point. The critical properties of the magnetic transition have been studied by Monte Carlo simulations at two densities in systems with up to 4000 particles. The locations of the critical points are in good agreement with the theoretical estimates. However, a finite-size scaling analysis shows effective critical exponents which deviate from the expected Ising values. These discrepancies resemble the analogous result of a previous numerical study of Heisenberg fluids. Here they are interpreted as due to crossover phenomena, partly related to Fisher exponent renormalization induced by the presence of noncritical density fluctuations.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jan 1, 1998•24 citations

Phase transition in the Ising ferromagnetic model with fixed spins

A. Labarta, Joaquín Marro, Benjamín Martínez +1 more

The ferromagnetic Ising spin-(1/2 model in a finite simple-cubic lattice \ensuremath{\Lambda} is studied by Monte Carlo methods when two subsets of the lattice sites in \ensuremath{\Lambda}, say ${\ensuremath{\Omega}}^{+}$ and ${\ensuremath{\Omega}}^{\mathrm{\ensuremath{-}}}$, contain (the same number of) spins fixed at \ifmmode\pm\else\textpm\fi{}1, respectively, the global defect concentration being x\ensuremath{\le}0.25. We study the thermodynamic properties of the model for different choices of \ensuremath{\Omega}=${\ensuremath{\Omega}}^{+}$\ensuremath{\cup}${\ensuremath{\Omega}}^{\mathrm{\ensuremath{-}}}$. A finite-size-scaling analysis reveals that the transition remains second order with pure critical exponents for regularly spaced defects, the critical temperature varying with the symmetry of \ensuremath{\Omega}. Any small randomness in \ensuremath{\Omega}, however, makes the transition weakly first order; the transition becomes more abrupt for defects located fully at random, and the long-range order is suppressed when the numbers of defects in ${\ensuremath{\Omega}}^{+}$ and ${\ensuremath{\Omega}}^{\mathrm{\ensuremath{-}}}$ differ from each other. We also discuss our findings in relation to the random-field and frustration problems.

Physical review. B, Condensed matter•Jul 1, 1988•3 citations

Possibility of Fisher renormalization of the critical exponents in an Ising fluid

W. Fenz, R. Folk, I. M. Mryglod +1 more

Using Monte Carlo simulation techniques, we study the ferromagnetic order-disorder phase transition in Ising spin fluids with hard-core Yukawa interaction truncated at various cutoff radii r{c}. We focus our interest on the dependence of critical quantities such as the Binder cumulant and various exponent ratios on the value of r{c}, and on the question whether the Fisher-renormalized exponents expected for such systems can be observed in the simulations. It turns out that the corrections to scaling decaying with a rather small exponent prevent reaching the asymptotic region with the computational power available. Thus, we observe only effective exponents, with different (nonuniversal) values depending on the cutoff radius. The same behavior is also found for the critical Binder cumulant. Nevertheless, an exact investigation of the effective susceptibility exponent gamma{eff} as a function of temperature seems to point towards a Fisher-renormalized value. For two selected cutoff radii, the critical temperature is determined more accurately using, in addition to the cumulant crossing technique, the scanning technique and the shifting technique, taking into account corrections to scaling. Simulations of Ising fluids with constant cutoff radius and varying Yukawa-tail screening lengths lambda also show a nonuniversal dependence of U{c} on lambda. Finally, we have performed simulations of the Ising lattice model with increasing number of couplings which show the expected asymptotic behavior, independent of the range of interactions.

Physical Review E•Jun 21, 2007•10 citations

Static critical behavior of three-dimensional classical Heisenberg models: A high-resolution Monte Carlo study

Kun Chen, Alan M. Ferrenberg, David Paul Landau

Using both recently developed cluster-algorithm and histogram methods, we have carried out a high-resolution Monte Carlo study of static critical properties of classical ferromagnetic Heisenberg models. Extensive Monte Carlo simulations were performed at several temperatures in the critical region, using an improved cluster-updating scheme, on L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L simple-cubic and body-centered-cubic systems with L\ensuremath{\le}40. Thermodynamic quantities as a function of temperature in the vicinity of the critical point were obtained by an optimized multiple-histogram method, and the critical temperature and static critical exponents were extracted using finite-size scaling. Our best estimates for the inverse critical temperatures are 0.693 035(37) for the simple-cubic system and 0.486 798(12) for the body-centered-cubic system. Estimated static critical exponents for both systems agree with each other within their respective error bars, and the mean estimates \ensuremath{\nu}=0.7048(30) and \ensuremath{\gamma}=1.3873(85) are also in excellent agreement with field-theoretic predictions 0.705(3) and 1.386(4).

Physical review. B, Condensed matter•Aug 1, 1993•331 citations
We have simulated a fluid of particles which carry a Heisenberg spin and interact via a Heisenberg-like interaction. The fluid displays a magnetic order-disorder transition, which we studied by varying the temperature at a fixed particle density. We located the transition and estimated critical exponents by means of a finite-size scaling analysis. The exponents differ by a small but significant amount from the ones for the lattice Heisenberg model. \textcopyright{} 1996 The American Physical Society.
We have performed Monte Carlo simulations of the ferromagnetic order-disorder transition in an off-lattice version of the Heisenberg model. Varying the temperature at a fixed density, we located the transition for three different densities by means of a finite-size scaling analysis. The obtained critical exponents differ from those for the lattice Heisenberg model. Although the limited size of our systems warrants caution, the results question whether the lattice and the off-lattice model are in the same universality class.
The magnetic phase transition in a Heisenberg fluid is studied by means of the finite size scaling technique. We find that even for larger systems, considered in an ensemble with fixed density, the critical exponents show deviations from the expected lattice values similar to those obtained previously. This puzzle is clarified by proving the importance of the leading correction to the scaling that appears due to Fisher renormalization with the critical exponent equal to the absolute value of the specific heat exponent alpha. The appearance of such new corretions to scaling is a general feature of systems with constraints.
The Ising fluid we have investigated is a system of hard spheres which carry an Ising spin. A short range ferromagnetic potential acts between particles. A renormalization group study of the model, based on the hierarchical reference theory of fluids, gives quantitative predictions on the phase diagram and the phase boundaries. The model can show two distinct phase transitions: a magnetic order-disorder transition and a liquid-vapor phase separation which meet at a tricritical point. The critical properties of the magnetic transition have been studied by Monte Carlo simulations at two densities in systems with up to 4000 particles. The locations of the critical points are in good agreement with the theoretical estimates. However, a finite-size scaling analysis shows effective critical exponents which deviate from the expected Ising values. These discrepancies resemble the analogous result of a previous numerical study of Heisenberg fluids. Here they are interpreted as due to crossover phenomena, partly related to Fisher exponent renormalization induced by the presence of noncritical density fluctuations.
The ferromagnetic Ising spin-(1/2 model in a finite simple-cubic lattice \ensuremath{\Lambda} is studied by Monte Carlo methods when two subsets of the lattice sites in \ensuremath{\Lambda}, say ${\ensuremath{\Omega}}^{+}$ and ${\ensuremath{\Omega}}^{\mathrm{\ensuremath{-}}}$, contain (the same number of) spins fixed at \ifmmode\pm\else\textpm\fi{}1, respectively, the global defect concentration being x\ensuremath{\le}0.25. We study the thermodynamic properties of the model for different choices of \ensuremath{\Omega}=${\ensuremath{\Omega}}^{+}$\ensuremath{\cup}${\ensuremath{\Omega}}^{\mathrm{\ensuremath{-}}}$. A finite-size-scaling analysis reveals that the transition remains second order with pure critical exponents for regularly spaced defects, the critical temperature varying with the symmetry of \ensuremath{\Omega}. Any small randomness in \ensuremath{\Omega}, however, makes the transition weakly first order; the transition becomes more abrupt for defects located fully at random, and the long-range order is suppressed when the numbers of defects in ${\ensuremath{\Omega}}^{+}$ and ${\ensuremath{\Omega}}^{\mathrm{\ensuremath{-}}}$ differ from each other. We also discuss our findings in relation to the random-field and frustration problems.
Using Monte Carlo simulation techniques, we study the ferromagnetic order-disorder phase transition in Ising spin fluids with hard-core Yukawa interaction truncated at various cutoff radii r{c}. We focus our interest on the dependence of critical quantities such as the Binder cumulant and various exponent ratios on the value of r{c}, and on the question whether the Fisher-renormalized exponents expected for such systems can be observed in the simulations. It turns out that the corrections to scaling decaying with a rather small exponent prevent reaching the asymptotic region with the computational power available. Thus, we observe only effective exponents, with different (nonuniversal) values depending on the cutoff radius. The same behavior is also found for the critical Binder cumulant. Nevertheless, an exact investigation of the effective susceptibility exponent gamma{eff} as a function of temperature seems to point towards a Fisher-renormalized value. For two selected cutoff radii, the critical temperature is determined more accurately using, in addition to the cumulant crossing technique, the scanning technique and the shifting technique, taking into account corrections to scaling. Simulations of Ising fluids with constant cutoff radius and varying Yukawa-tail screening lengths lambda also show a nonuniversal dependence of U{c} on lambda. Finally, we have performed simulations of the Ising lattice model with increasing number of couplings which show the expected asymptotic behavior, independent of the range of interactions.
Using both recently developed cluster-algorithm and histogram methods, we have carried out a high-resolution Monte Carlo study of static critical properties of classical ferromagnetic Heisenberg models. Extensive Monte Carlo simulations were performed at several temperatures in the critical region, using an improved cluster-updating scheme, on L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L simple-cubic and body-centered-cubic systems with L\ensuremath{\le}40. Thermodynamic quantities as a function of temperature in the vicinity of the critical point were obtained by an optimized multiple-histogram method, and the critical temperature and static critical exponents were extracted using finite-size scaling. Our best estimates for the inverse critical temperatures are 0.693 035(37) for the simple-cubic system and 0.486 798(12) for the body-centered-cubic system. Estimated static critical exponents for both systems agree with each other within their respective error bars, and the mean estimates \ensuremath{\nu}=0.7048(30) and \ensuremath{\gamma}=1.3873(85) are also in excellent agreement with field-theoretic predictions 0.705(3) and 1.386(4).