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:2112.03783 · Physical review. B./Physical review. B

Three-dimensional O(N) -invariant ϕ4 models at criticality for N≥4

Martin Hasenbusch B

Published Feb 24, 2022• 23 citations
arXivHTMLPDF

Abstract

We study the $O(N)$-invariant ${\ensuremath{\phi}}^{4}$ model on the simple cubic lattice by using Monte Carlo simulations. By using a finite-size scaling analysis, we obtain accurate estimates for the critical exponents $\ensuremath{\nu}$ and $\ensuremath{\eta}$ for $N=4$, 5, 6, 8, 10, and 12. We study the model for each $N$ for at least three different values of the parameter $\ensuremath{\lambda}$ to control leading corrections to scaling. We compare our results with those obtained by other theoretical methods.

DOI
10.1103/physrevb.105.054428

Related papers

Cubic fixed point in three dimensions: Monte Carlo simulations of the ϕ4 model on the simple cubic lattice

Martin Hasenbusch B

We study the cubic fixed point for $N=3$ and 4 by using finite-size scaling applied to data obtained from Monte Carlo simulations of the $N$-component ${\ensuremath{\phi}}^{4}$ model on the simple cubic lattice. We generalize the idea of improved models to a two-parameter family of models. The two-parameter space is scanned for the point, where the amplitudes of the two leading corrections to scaling vanish. To this end, a dimensionless quantity is introduced that monitors the breaking of the $O(N)$ invariance. For $N=4$, we determine the correction exponents ${\ensuremath{\omega}}_{1}=0.763(24)$ and ${\ensuremath{\omega}}_{2}=0.082(5)$. In the case of $N=3$, we obtain ${Y}_{4}=0.0142(6)$ for the renormalization group exponent of the cubic perturbation at the $O(3)$-invariant fixed point, while the correction exponent ${\ensuremath{\omega}}_{2}=0.0133(8)$ at the cubic fixed point. Simulations close to the improved point result in the estimates $\ensuremath{\nu}=0.7202(7)$ and $\ensuremath{\eta}=0.0371(2)$ of the critical exponents of the cubic fixed point for $N=4$. For $N=3$, at the cubic fixed point, the $O(3)$ symmetry is only mildly broken and the critical exponents differ only by little from those of the $O(3)$-invariant fixed point. We find $\ensuremath{-}0.000\phantom{\rule{0.16em}{0ex}}01⪅{\ensuremath{\eta}}_{\text{cubic}}\ensuremath{-}{\ensuremath{\eta}}_{O(3)}⪅0.000\phantom{\rule{0.16em}{0ex}}07$ and ${\ensuremath{\nu}}_{\text{cubic}}\ensuremath{-}{\ensuremath{\nu}}_{O(3)}=\ensuremath{-}0.000\phantom{\rule{0.16em}{0ex}}61(10)$.

arXiv:2211.16170•Physical review. B./Physical review. B•Jan 10, 2023•23 citations

Finite size scaling study of lattice models in the three-dimensional Ising universality class

Martin Hasenbusch B

We study the spin-1/2 Ising model and the Blume-Capel model at various values of the parameter $D$ on the simple cubic lattice. To this end we perform Monte Carlo simulations using a hybrid of the local Metropolis, the single cluster and the wall cluster algorithm. Using finite size scaling we determine the value ${D}^{\ensuremath{\ast}}=0.656(20)$ of the parameter $D$, where leading corrections to scaling vanish. We find $\ensuremath{\omega}=0.832(6)$ for the exponent of leading corrections to scaling. In order to compute accurate estimates of critical exponents, we construct improved observables that have a small amplitude of the leading correction for any model. Analyzing data obtained for $D=0.641$ and 0.655 on lattices of a linear size up to $L=360$ we obtain $\ensuremath{\nu}=0.63002(10)$ and $\ensuremath{\eta}=0.03627(10)$. We compare our results with those obtained from previous Monte Carlo simulations and high-temperature series expansions of lattice models, by using field-theoretic methods and experiments.

arXiv:1004.4486•Physical Review B•Nov 23, 2010•333 citations

Monte Carlo study of an improved clock model in three dimensions

Martin Hasenbusch B

We study a generalized clock model on the simple cubic lattice. The parameter of the model can be tuned such that the amplitude of the leading correction to scaling vanishes. In the main part of the study, we simulate the model with ${Z}_{8}$ symmetry. At the transition, with increasing length scale, $O(2)$ symmetry emerges. We perform Monte Carlo simulations using a hybrid of local Metropolis and cluster algorithms of lattices with a linear size up to $L=512$. The field variable requires less memory and the updates are faster than for a model with $O(2)$ symmetry at the microscopic level. Our finite-size scaling analysis yields accurate estimates for the critical exponents of the three-dimensional $XY$-universality class. In particular, we get $\ensuremath{\eta}=0.03810(8),\ensuremath{\nu}=0.67169(7)$, and $\ensuremath{\omega}=0.789(4)$. Furthermore, we obtain estimates for fixed point values of phenomenological couplings and critical temperatures.

arXiv:1910.05916•Physical review. B./Physical review. B•Dec 26, 2019•107 citations

Monte Carlo study of a generalized icosahedral model on the simple cubic lattice

Martin Hasenbusch B

We study the critical behavior of a generalized icosahedral model on the simple cubic lattice. The field variable of the icosahedral model might take 1 of 12 vectors of unit length which are given by the normalized vertices of the icosahedron as value. Similar to the Blume-Capel model, where in addition to $\ensuremath{-}1$ and 1, as in the Ising model, the spin might take the value 0, we add in the generalized model (0,0,0) as allowed value. There is a parameter $D$ that controls the density of these voids. For a certain range of $D$, the model undergoes a second-order phase transition. On the critical line, $\text{O}(3)$ symmetry emerges. Furthermore, we demonstrate that within this range, similar to the Blume-Capel model on the simple cubic lattice, there is a value of $D$, where leading corrections to scaling vanish. We perform Monte Carlo simulations for lattices of a linear size up to $L=400$ by using a hybrid of local Metropolis and cluster updates. The motivation to study this particular model is mainly of technical nature. Less memory and CPU time are needed than for a model with $\text{O}(3)$ symmetry at the microscopic level. As the result of a finite-size scaling analysis we obtain $\ensuremath{\nu}=0.711\phantom{\rule{0.16em}{0ex}}64(10)$, $\ensuremath{\eta}=0.037\phantom{\rule{0.16em}{0ex}}84(5)$, and $\ensuremath{\omega}=0.759(2)$ for the critical exponents of the three-dimensional Heisenberg universality class. The estimate of the irrelevant renormalization group eigenvalue that is related with the breaking the $\text{O}(3)$ symmetry is ${y}_{\mathrm{ico}}=\ensuremath{-}2.19(2)$.

arXiv:2005.04448•Physical review. B./Physical review. B•Jul 6, 2020•52 citations

High-precision Monte Carlo study of the 3DXY-universality class

M Hasenbusch, T Török

We present a Monte Carlo study of the two-component 4 model on the simple cubic lattice in three dimensions. By suitable tuning of the coupling constant, , we eliminate leading-order corrections to scaling. High-statistics simulations using finite-size scaling techniques yield = 0.6723(3)[8] and = 0.0381(2)[2], where the statistical and systematical errors are given in the first and second bracket, respectively. These results are more precise than any previous theoretical estimate of the critical exponents for the 3D XY universality class.

arXiv:cond-mat/9904408•Journal of Physics A Mathematical and General•Aug 26, 1999•113 citations

Monte Carlo test of critical exponents and amplitudes in 3D Ising and phi^4 lattice models

J. Kaupužs

We have tested the leading correction-to-scaling exponent omega in O(n)-symmetric models on a three-dimensional lattice by analysing the recent Monte Carlo (MC) data. We have found that the effective critical exponent, estimated at finite sizes of the system L and L/2, decreases remarkably within the range of the simulated L values. This shows the incorrectness of some claims that omega has a very accurate value 0.845(10) at n=1. A selfconsistent infinite volume extrapolation yields row estimates omega=0.547, omega=0.573, and omega=0.625 at n=1, 2, and 3, respectively, in approximate agreement with the corresponding exact values 1/2, 5/9, and 3/5 predicted by our recently developed GFD (grouping of Feynman diagrams) theory. We have fitted the MC data for the susceptibility of 3D Ising model at criticality showing that the effective critical exponent eta tends to increase well above the usually accepted values around 0.036. We have fitted the data within [L;8L], including several terms in the asymptotic expansion with fixed exponents, to obtain the effective amplitudes depending on L. This method clearly demonstrates that the critical exponents of GFD theory are correct (the amplitudes converge to certain asymptotic values at L tending to infinity), whereas those of the perturbative renormalization group (RG) theory are incorrect (the amplitudes diverge). A modification of the standard Ising model by introducing suitable "improved" action (Hamiltonian) does not solve the problem in favour of the perturbative RG theory.

arXiv:cond-mat/0106271•Jun 14, 2001•0 citations
We study the O(N)O(N)O(N)-invariant \ensuremathϕ4{\ensuremath{\phi}}^{4}\ensuremathϕ4 model on the simple cubic lattice by using Monte Carlo simulations. By using a finite-size scaling analysis, we obtain accurate estimates for the critical exponents \ensuremathν\ensuremath{\nu}\ensuremathν and \ensuremathη\ensuremath{\eta}\ensuremathη for N=4N=4N=4, 5, 6, 8, 10, and 12. We study the model for each NNN for at least three different values of the parameter \ensuremathλ\ensuremath{\lambda}\ensuremathλ to control leading corrections to scaling. We compare our results with those obtained by other theoretical methods.
We study the cubic fixed point for N=3N=3N=3 and 4 by using finite-size scaling applied to data obtained from Monte Carlo simulations of the NNN-component ${\ensuremath{\phi}}^{4}$ model on the simple cubic lattice. We generalize the idea of improved models to a two-parameter family of models. The two-parameter space is scanned for the point, where the amplitudes of the two leading corrections to scaling vanish. To this end, a dimensionless quantity is introduced that monitors the breaking of the O(N)O(N)O(N) invariance. For N=4N=4N=4, we determine the correction exponents ${\ensuremath{\omega}}_{1}=0.763(24)$ and ${\ensuremath{\omega}}_{2}=0.082(5)$. In the case of N=3N=3N=3, we obtain Y4=0.0142(6){Y}_{4}=0.0142(6)Y4​=0.0142(6) for the renormalization group exponent of the cubic perturbation at the O(3)O(3)O(3)-invariant fixed point, while the correction exponent ${\ensuremath{\omega}}_{2}=0.0133(8)$ at the cubic fixed point. Simulations close to the improved point result in the estimates $\ensuremath{\nu}=0.7202(7)$ and $\ensuremath{\eta}=0.0371(2)$ of the critical exponents of the cubic fixed point for N=4N=4N=4. For N=3N=3N=3, at the cubic fixed point, the O(3)O(3)O(3) symmetry is only mildly broken and the critical exponents differ only by little from those of the O(3)O(3)O(3)-invariant fixed point. We find $\ensuremath{-}0.000\phantom{\rule{0.16em}{0ex}}01⪅{\ensuremath{\eta}}_{\text{cubic}}\ensuremath{-}{\ensuremath{\eta}}_{O(3)}⪅0.000\phantom{\rule{0.16em}{0ex}}07$ and ${\ensuremath{\nu}}_{\text{cubic}}\ensuremath{-}{\ensuremath{\nu}}_{O(3)}=\ensuremath{-}0.000\phantom{\rule{0.16em}{0ex}}61(10)$.
We study the spin-1/2 Ising model and the Blume-Capel model at various values of the parameter DDD on the simple cubic lattice. To this end we perform Monte Carlo simulations using a hybrid of the local Metropolis, the single cluster and the wall cluster algorithm. Using finite size scaling we determine the value ${D}^{\ensuremath{\ast}}=0.656(20)$ of the parameter DDD, where leading corrections to scaling vanish. We find $\ensuremath{\omega}=0.832(6)$ for the exponent of leading corrections to scaling. In order to compute accurate estimates of critical exponents, we construct improved observables that have a small amplitude of the leading correction for any model. Analyzing data obtained for D=0.641D=0.641D=0.641 and 0.655 on lattices of a linear size up to L=360L=360L=360 we obtain $\ensuremath{\nu}=0.63002(10)$ and $\ensuremath{\eta}=0.03627(10)$. We compare our results with those obtained from previous Monte Carlo simulations and high-temperature series expansions of lattice models, by using field-theoretic methods and experiments.
We study a generalized clock model on the simple cubic lattice. The parameter of the model can be tuned such that the amplitude of the leading correction to scaling vanishes. In the main part of the study, we simulate the model with Z8{Z}_{8}Z8​ symmetry. At the transition, with increasing length scale, O(2)O(2)O(2) symmetry emerges. We perform Monte Carlo simulations using a hybrid of local Metropolis and cluster algorithms of lattices with a linear size up to L=512L=512L=512. The field variable requires less memory and the updates are faster than for a model with O(2)O(2)O(2) symmetry at the microscopic level. Our finite-size scaling analysis yields accurate estimates for the critical exponents of the three-dimensional XYXYXY-universality class. In particular, we get $\ensuremath{\eta}=0.03810(8),\ensuremath{\nu}=0.67169(7)$, and $\ensuremath{\omega}=0.789(4)$. Furthermore, we obtain estimates for fixed point values of phenomenological couplings and critical temperatures.
We study the critical behavior of a generalized icosahedral model on the simple cubic lattice. The field variable of the icosahedral model might take 1 of 12 vectors of unit length which are given by the normalized vertices of the icosahedron as value. Similar to the Blume-Capel model, where in addition to $\ensuremath{-}1$ and 1, as in the Ising model, the spin might take the value 0, we add in the generalized model (0,0,0) as allowed value. There is a parameter DDD that controls the density of these voids. For a certain range of DDD, the model undergoes a second-order phase transition. On the critical line, O(3)\text{O}(3)O(3) symmetry emerges. Furthermore, we demonstrate that within this range, similar to the Blume-Capel model on the simple cubic lattice, there is a value of DDD, where leading corrections to scaling vanish. We perform Monte Carlo simulations for lattices of a linear size up to L=400L=400L=400 by using a hybrid of local Metropolis and cluster updates. The motivation to study this particular model is mainly of technical nature. Less memory and CPU time are needed than for a model with O(3)\text{O}(3)O(3) symmetry at the microscopic level. As the result of a finite-size scaling analysis we obtain $\ensuremath{\nu}=0.711\phantom{\rule{0.16em}{0ex}}64(10)$, $\ensuremath{\eta}=0.037\phantom{\rule{0.16em}{0ex}}84(5)$, and $\ensuremath{\omega}=0.759(2)$ for the critical exponents of the three-dimensional Heisenberg universality class. The estimate of the irrelevant renormalization group eigenvalue that is related with the breaking the O(3)\text{O}(3)O(3) symmetry is ${y}_{\mathrm{ico}}=\ensuremath{-}2.19(2)$.
We present a Monte Carlo study of the two-component 4 model on the simple cubic lattice in three dimensions. By suitable tuning of the coupling constant, , we eliminate leading-order corrections to scaling. High-statistics simulations using finite-size scaling techniques yield = 0.6723(3)[8] and = 0.0381(2)[2], where the statistical and systematical errors are given in the first and second bracket, respectively. These results are more precise than any previous theoretical estimate of the critical exponents for the 3D XY universality class.
We have tested the leading correction-to-scaling exponent omega in O(n)-symmetric models on a three-dimensional lattice by analysing the recent Monte Carlo (MC) data. We have found that the effective critical exponent, estimated at finite sizes of the system L and L/2, decreases remarkably within the range of the simulated L values. This shows the incorrectness of some claims that omega has a very accurate value 0.845(10) at n=1. A selfconsistent infinite volume extrapolation yields row estimates omega=0.547, omega=0.573, and omega=0.625 at n=1, 2, and 3, respectively, in approximate agreement with the corresponding exact values 1/2, 5/9, and 3/5 predicted by our recently developed GFD (grouping of Feynman diagrams) theory. We have fitted the MC data for the susceptibility of 3D Ising model at criticality showing that the effective critical exponent eta tends to increase well above the usually accepted values around 0.036. We have fitted the data within [L;8L], including several terms in the asymptotic expansion with fixed exponents, to obtain the effective amplitudes depending on L. This method clearly demonstrates that the critical exponents of GFD theory are correct (the amplitudes converge to certain asymptotic values at L tending to infinity), whereas those of the perturbative renormalization group (RG) theory are incorrect (the amplitudes diverge). A modification of the standard Ising model by introducing suitable "improved" action (Hamiltonian) does not solve the problem in favour of the perturbative RG theory.