arXiv:2405.15233 · AIP Advances

In-depth investigation of phase transition phenomena in network models derived from lattice models

Yonglong Ding

Published Aug 1, 2024• 1 citation

Abstract

Lattice models exhibit significant potential in investigating phase transitions, yet they encounter numerous computational challenges. To address these issues, this study introduces a Monte Carlo-based approach that transforms lattice models into a network model with intricate inter-node correlations. This framework enables a profound analysis of Ising, JQ, and XY models. By decomposing the network into a maximum entropy component and a conservative component, under the constraint of detailed balance, this work derives an estimation formula for the temperature-dependent magnetic induction in Ising models. Notably, the critical exponent β in the Ising model aligns well with the established results, and the predicted phase transition point in the three-dimensional Ising model exhibits a mere 0.7% deviation from numerical simulations.

DOI: 10.1063/5.0219207

Related papers

Pushing the limits of Monte Carlo simulations for the three-dimensional Ising model

Alan M. Ferrenberg, Jiahao Xu, David Paul Landau

While the three-dimensional Ising model has defied analytic solution, various numerical methods like Monte Carlo, Monte Carlo renormalization group, and series expansion have provided precise information about the phase transition. Using Monte Carlo simulation that employs the Wolff cluster flipping algorithm with both 32-bit and 53-bit random number generators and data analysis with histogram reweighting and quadruple precision arithmetic, we have investigated the critical behavior of the simple cubic Ising Model, with lattice sizes ranging from 16^{3} to 1024^{3}. By analyzing data with cross correlations between various thermodynamic quantities obtained from the same data pool, e.g., logarithmic derivatives of magnetization and derivatives of magnetization cumulants, we have obtained the critical inverse temperature K_{c}=0.221654626(5) and the critical exponent of the correlation length ν=0.629912(86) with precision that exceeds all previous Monte Carlo estimates.

arXiv:1806.03558Physical review. EApr 4, 2018191 citations

On the Ising model for the simple cubic lattice

Roland Häggkvist, Anders Rosengren, P. H. Lundow +3 more

The Ising model was introduced in 1920 to describe a uniaxial system of magnetic moments, localized on a lattice, interacting via nearest-neighbour exchange interaction. It is the generic model for a continuous phase transition and arguably the most studied model in theoretical physics. Since it was solved for a two-dimensional lattice by Onsager in 1944, thereby representing one of the very few exactly solvable models in dimensions higher than one, it has served as a testing ground for new developments in analytic treatment and numerical algorithms. Only series expansions and numerical approaches, such as Monte Carlo simulations, are available in three dimensions. This review focuses on Monte Carlo simulation. We build upon a data set of unprecedented size. A great number of quantities of the model are estimated near the critical coupling. We present both a conventional analysis and an analysis in terms of a Puiseux series for the critical exponents. The former gives distinct values of the high- and low-temperature exponents; by means of the latter we can get these exponents to be equal at the cost of having true asymptotic behaviour being found only extremely close to the critical point. The consequences of this for simulations of lattice systems are discussed at length.

Advances In PhysicsSep 1, 200740 citations

Cluster Monte Carlo and network flow algorithms for critical phenomena

Ilija Dukovski

Phase transitions and critical phenomena in several spin models are studied. These models include a two and three dimensional XY model and Ising models in structured and random field. The methods of study are cluster Monte Carlo simulations and network flow optimization algorithms. The XY model and a related soft spin model belonging to the O(2) universality class are studied using a modified version of the Invaded Cluster algorithm. The static critical properties of the models and the dynamic characteristics of the algorithm are studied. The critical exponent η for the O(2) universality class in two and three dimensions is measured. The Invaded Cluster algorithm for the XY model does not show any critical slowing for global variables. Benzene adsorption in a nanoporous solid is modeled with an Ising model with both internal and external magnetic fields. The first order phase transition and the critical properties of the system are studied. The system is simulated using the two replica cluster algorithm. The critical temperature of the system is found to decrease with increasing difference of the values of the internal fields. It is found that the system undergoes a phase transition for a relatively wide, experimentally obtainable, region of values of its parameters. The random field Ising model in three dimensions is studied at zero temperature. The nature of randomness induced phase transition is studied using a numerical network flow optimization method. It is found that the latent heat at the phase transition vanishes in the limit of infinite system size which can be taken as a proof of a second order phase transition.

Jan 1, 20010 citations

Monte-Carlo simulation of phase transition in 2D and 3D ising model

Rajendra Neupane, Hari Ram Krishna Gauli, Krishna Bahadur Rai +1 more

In this work, Markov Chain-Monte Carlo technique was used to study the phase transition in two and three dimensional Ising Model (IM) in a square and cubic lattice. The study of temperature dependence of average magnetization and specific heat in different magnetic fields has been carried out in the 3x3 and 3x3x3 lattice with periodic boundary. Critical temperature point kBTc / J for 2D and 3D Ising Model has been observed at around 2.2 and 4.3 respectively at zero field. Our work satisfies Onsager’s critical value in 2D IM. The simulation suggests bifurcation in average magnetization below critical temperature Tc. Temperature plays the role of increasing randomness of spins. We found that Ising Model in small lattice size still retains interesting features like spontaneous magnetization and symmetry breaking below Tc at B = 0. At a non-zero field, the likelihood of spins to prefer certain alignment depends on the direction of the external field and magnitude of magnetization depends on magnitude of field ±B. Specific heat Cv, which gives us fluctuation of energy at particular temperature, has been found to attain maximum value at critical regions.

Scientific WorldAug 4, 20237 citations

Ising model on three-dimensional random lattices: A Monte Carlo study

Wolfhard Janke E, Ramon Villanova

We report single-cluster Monte Carlo simulations of the Ising model on three-dimensional Poissonian random lattices with up to $128000\ensuremath{\approx}{50}^{3}$ sites which are linked together according to the Voronoi-Delaunay prescription. For each lattice size quenched averages are performed over 96 realizations. By using reweighting techniques and finite-size scaling analyses we investigate the critical properties of the model in the close vicinity of the phase transition point. Our random lattice data provide strong evidence that, for the available system sizes, the resulting effective critical exponents are indistinguishable from recent high-precision estimates obtained in Monte Carlo studies of the Ising model and ${\ensuremath{\varphi}}^{4}$ field theory on three-dimensional regular cubic lattices.

arXiv:cond-mat/0211446Physical review. B, Condensed matterOct 31, 200231 citations

Accurate Monte Carlo critical exponents for Ising lattices

J. G. García, Julio A. Gonzalo, Manuel I. Marqués

A careful Monte Carlo investigation of the phase transition very close to the critical point (T -> Tc, H -> 0) in relatively large d = 3, s = 1/2 Ising lattices did produce critical exponents beta = 0.3126(4) =~ 5/16, delta^{-1} = 0.1997(4) =~ 1/5 and gamma_{3D} = 1.253(4) =~ 5/4. Our results indicate that, within experimental error, they are given by simple fractions corresponding to the linear interpolations between the respective two-dimensional (Onsager) and four-dimensional (mean field) critical exponents. An analysis of our inverse susceptibility data chi^{-1}(T) vs. /T - Tc/ shows that these data lead to a value of gamma_{3D} compatible with gamma' = gamma and Tc = 4.51152(12), while gamma values obtained recently by high and low temperature series expansions and renormalization group methods are not.

arXiv:cond-mat/0211270Nov 14, 20020 citations