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:1502.04249 · Physical Review E

Fourier Monte Carlo renormalization-group approach to crystalline membranes

A. Tröster

Published Feb 23, 2015• 13 citations
arXivHTMLPDF

Abstract

The computation of the critical exponent η characterizing the universal elastic behavior of crystalline membranes in the flat phase continues to represent challenges to theorists as well as computer simulators that manifest themselves in a considerable spread of numerical results for η published in the literature. We present additional insight into this problem that results from combining Wilson's momentum shell renormalization-group method with the power of modern computer simulations based on the Fourier Monte Carlo algorithm. After discussing the ideas and difficulties underlying this combined scheme, we present a calculation of the renormalization-group flow of the effective two-dimensional Young modulus for momentum shells of different thickness. Extrapolation to infinite shell thickness allows us to produce results in reasonable agreement with those obtained by functional renormalization group or by Fourier Monte Carlo simulations in combination with finite-size scaling. Moreover, our method allows us to obtain a decent estimate for the value of the Wegner exponent ω that determines the leading correction to scaling, which in turn allows us to refine our numerical estimate for η previously obtained from precise finite-size scaling data.

DOI
10.1103/physreve.91.022132

Related papers

High-precision Fourier Monte Carlo simulation of crystalline membranes

A. Tröster

We report an essential improvement of the plain Fourier Monte Carlo algorithm that promises to be a powerful tool for investigating critical behavior in a large class of lattice models, in particular those containing microscopic or effective long-ranged interactions. On tuning the Monte Carlo acceptance rates separately for each wave vector, we are able to drastically reduce critical slowing down. We illustrate the resulting efficiency and unprecedented accuracy of our algorithm with a calculation of the universal elastic properties of crystalline membranes in the flat phase and derive a numerical estimate $\ensuremath{\eta}=0.795(10)$ for the critical exponent $\ensuremath{\eta}$ that challenges those derived from other recent simulations. The large system sizes accessible to our present algorithm also allow us to demonstrate that insufficiently taking into account corrections to scaling may severely hamper a finite size scaling analysis. This observation may also help to clarify the apparent disagreement of published numerical estimates of $\ensuremath{\eta}$ in the existing literature.

arXiv:1303.3726•Physical Review B•Mar 25, 2013•39 citations

Optimizing Wilson’s momentum shell renormalization group

A. Tröster

Previous attempts to accurately compute critical exponents from Wilson's momentum shell renormalization prescription suffered from the difficulties posed by the presence of an infinite number of irrelevant couplings. Taking the example of the one-dimensional long-ranged Ising model, we calculate the momentum shell renormalization flow in the plane spanned by the coupling constants $({u}_{0},{r}_{0})$ for different values of the momentum shell thickness parameter $b$ by simulation using our recently developed Fourier Monte Carlo algorithm. We report strong anomalies in the $b$ dependence of the fixed-point couplings and the resulting exponents ${y}_{\ensuremath{\tau}}$ and $\ensuremath{\omega}$ in the vicinity of a shell parameter ${b}^{\ensuremath{\ast}}<1$ characterizing a thin but finite momentum shell. Evaluation of the exponents for this value of $b$ yields a dramatic improvement of their numerical accuracy, indicating a strong damping of the influence of irrelevant couplings for $b={b}^{\ensuremath{\ast}}$.

Physical Review B•Mar 31, 2010•8 citations

Momentum-shell renormalization-group flow from simulation

A. Tröster

Our recently developed Fourier Monte Carlo algorithm permits a nonperturbative calculation of momentum-shell renormalization-group flows by simulation which despite its apparent simplicity is illustrative both numerically as well as conceptually interesting. We study the example of a varphi(4) model with long-range lattice interaction. For this model we show that the topology of the renormalization flow is globally accessible in a particularly convenient way. The nontrivial fixed point of Wilson-Fisher type is observed its accompanying critical exponents are numerically determined from fitting its surrounding flow pattern to a linearized renormalization-group transformation. The results are compared to those obtained from perturbation theory, -expansion and earlier Monte Carlo simulations. Application of our method is also expected to be rewarding in other models with long-range interactions.

Physical Review E•Mar 26, 2009•10 citations

Corrections to scaling in the 3D Ising model: a comparison between MC and MCRG results

J. Kaupužs, Roderick Melnik

Corrections to scaling in the 3D Ising model are studied based on Monte Carlo (MC) simulation results for very large lattices with linear lattice sizes up to L=3456. Our estimated values of the correction-to-scaling exponent omega tend to decrease below the usually accepted value about 0.83 when the smallest lattice sizes are discarded from the fits. This behavior apparently confirms some of the known estimates of the Monte Carlo renormalization group (MCRG) method, i.e., omega about 0.7 and omega = 0.75(5). We discuss the possibilities that omega is either really smaller than usually expected or these values of omega describe some transient behavior. We propose refining MCRG simulations and analysis to resolve this issue. In distinction from omega, our actual MC estimations of the critical exponents eta and nu provide stable values eta=0.03632(13) and nu=0.63017(31), which well agree with those of the conformal bootstrap method, i.e., eta=0.0362978(20) and nu=0.6299709(40).

arXiv:2206.04393•Jun 9, 2022•0 citations

Monte Carlo, blocking, and inference: How to measure the renormalization group flow

Luca Di Carlo

Renormalization group theory is a powerful and intriguing technique with a wide range of applications. One of the main successes of renormalization group theory is the description of continuous phase transitions and the development of scaling theory. Most courses on phase transitions focus on scaling and critical exponents, while less attention is paid to universality, renormalization group flow, and the existence of a unique fixed point, which are the ultimate reasons why scaling theory is so effective in describing continuous phase transitions. We use a combination of Monte Carlo simulations and real space renormalization group theory to determine the renormalization group flow and to show the existence of a universal fixed point in the context of the ferromagnetic Ising model.

arXiv:2401.04811•American Journal of Physics•Jan 23, 2025•0 citations

Transfer matrix and Monte Carlo tests of critical exponents in lattice models

J. Kaupužs

The corrections to finite-size scaling in the critical two-point correlation function G(r) of 2D Ising model on a square lattice have been studied numerically by means of exact transfer-matrix algorithms. The systems of square geometry with periodic boundaries oriented either along <10> or along <11> direction have been considered, including up to 800 spins. The calculation of G(r) at a distance r equal to the half of the system size L shows the existence of an amplitude correction proportional to 1/L^2. A nontrivial correction proportional to 1/L^0.25 of a very small magnitude also has been detected in agreement with predictions of our recently developed GFD (grouping of Feynman diagrams) theory. A refined analysis of the recent MC data for 3D Ising, phi^4, and XY lattice models has been performed. It includes an analysis of the partition function zeros of 3D Ising model, an estimation of the correction-to-scaling exponent omega from the Binder cumulant data near criticality, as well as a study of the effective critical exponent eta and the effective amplitudes in the asymptotic expansion of susceptibility at the critical point. In all cases a refined analysis is consistent with our (GFD) asymptotic values of the critical exponents (nu=2/3, omega=1/2, eta=1/8 for 3D Ising model and omega=5/9 for 3D XY model), while the actually accepted "conventional" exponents are, in fact, effective exponents which are valid for approximation of the finite--size scaling behavior of not too large systems.

arXiv:cond-mat/0201221•Jan 14, 2002•0 citations
The computation of the critical exponent η characterizing the universal elastic behavior of crystalline membranes in the flat phase continues to represent challenges to theorists as well as computer simulators that manifest themselves in a considerable spread of numerical results for η published in the literature. We present additional insight into this problem that results from combining Wilson's momentum shell renormalization-group method with the power of modern computer simulations based on the Fourier Monte Carlo algorithm. After discussing the ideas and difficulties underlying this combined scheme, we present a calculation of the renormalization-group flow of the effective two-dimensional Young modulus for momentum shells of different thickness. Extrapolation to infinite shell thickness allows us to produce results in reasonable agreement with those obtained by functional renormalization group or by Fourier Monte Carlo simulations in combination with finite-size scaling. Moreover, our method allows us to obtain a decent estimate for the value of the Wegner exponent ω that determines the leading correction to scaling, which in turn allows us to refine our numerical estimate for η previously obtained from precise finite-size scaling data.
We report an essential improvement of the plain Fourier Monte Carlo algorithm that promises to be a powerful tool for investigating critical behavior in a large class of lattice models, in particular those containing microscopic or effective long-ranged interactions. On tuning the Monte Carlo acceptance rates separately for each wave vector, we are able to drastically reduce critical slowing down. We illustrate the resulting efficiency and unprecedented accuracy of our algorithm with a calculation of the universal elastic properties of crystalline membranes in the flat phase and derive a numerical estimate $\ensuremath{\eta}=0.795(10)$ for the critical exponent $\ensuremath{\eta}$ that challenges those derived from other recent simulations. The large system sizes accessible to our present algorithm also allow us to demonstrate that insufficiently taking into account corrections to scaling may severely hamper a finite size scaling analysis. This observation may also help to clarify the apparent disagreement of published numerical estimates of $\ensuremath{\eta}$ in the existing literature.
Previous attempts to accurately compute critical exponents from Wilson's momentum shell renormalization prescription suffered from the difficulties posed by the presence of an infinite number of irrelevant couplings. Taking the example of the one-dimensional long-ranged Ising model, we calculate the momentum shell renormalization flow in the plane spanned by the coupling constants (u0,r0)({u}_{0},{r}_{0})(u0​,r0​) for different values of the momentum shell thickness parameter bbb by simulation using our recently developed Fourier Monte Carlo algorithm. We report strong anomalies in the bbb dependence of the fixed-point couplings and the resulting exponents ${y}_{\ensuremath{\tau}}$ and $\ensuremath{\omega}$ in the vicinity of a shell parameter ${b}^{\ensuremath{\ast}}<1$ characterizing a thin but finite momentum shell. Evaluation of the exponents for this value of bbb yields a dramatic improvement of their numerical accuracy, indicating a strong damping of the influence of irrelevant couplings for $b={b}^{\ensuremath{\ast}}$.
Our recently developed Fourier Monte Carlo algorithm permits a nonperturbative calculation of momentum-shell renormalization-group flows by simulation which despite its apparent simplicity is illustrative both numerically as well as conceptually interesting. We study the example of a varphi(4) model with long-range lattice interaction. For this model we show that the topology of the renormalization flow is globally accessible in a particularly convenient way. The nontrivial fixed point of Wilson-Fisher type is observed its accompanying critical exponents are numerically determined from fitting its surrounding flow pattern to a linearized renormalization-group transformation. The results are compared to those obtained from perturbation theory, -expansion and earlier Monte Carlo simulations. Application of our method is also expected to be rewarding in other models with long-range interactions.
Corrections to scaling in the 3D Ising model are studied based on Monte Carlo (MC) simulation results for very large lattices with linear lattice sizes up to L=3456. Our estimated values of the correction-to-scaling exponent omega tend to decrease below the usually accepted value about 0.83 when the smallest lattice sizes are discarded from the fits. This behavior apparently confirms some of the known estimates of the Monte Carlo renormalization group (MCRG) method, i.e., omega about 0.7 and omega = 0.75(5). We discuss the possibilities that omega is either really smaller than usually expected or these values of omega describe some transient behavior. We propose refining MCRG simulations and analysis to resolve this issue. In distinction from omega, our actual MC estimations of the critical exponents eta and nu provide stable values eta=0.03632(13) and nu=0.63017(31), which well agree with those of the conformal bootstrap method, i.e., eta=0.0362978(20) and nu=0.6299709(40).
Renormalization group theory is a powerful and intriguing technique with a wide range of applications. One of the main successes of renormalization group theory is the description of continuous phase transitions and the development of scaling theory. Most courses on phase transitions focus on scaling and critical exponents, while less attention is paid to universality, renormalization group flow, and the existence of a unique fixed point, which are the ultimate reasons why scaling theory is so effective in describing continuous phase transitions. We use a combination of Monte Carlo simulations and real space renormalization group theory to determine the renormalization group flow and to show the existence of a universal fixed point in the context of the ferromagnetic Ising model.
The corrections to finite-size scaling in the critical two-point correlation function G(r) of 2D Ising model on a square lattice have been studied numerically by means of exact transfer-matrix algorithms. The systems of square geometry with periodic boundaries oriented either along <10> or along <11> direction have been considered, including up to 800 spins. The calculation of G(r) at a distance r equal to the half of the system size L shows the existence of an amplitude correction proportional to 1/L^2. A nontrivial correction proportional to 1/L^0.25 of a very small magnitude also has been detected in agreement with predictions of our recently developed GFD (grouping of Feynman diagrams) theory. A refined analysis of the recent MC data for 3D Ising, phi^4, and XY lattice models has been performed. It includes an analysis of the partition function zeros of 3D Ising model, an estimation of the correction-to-scaling exponent omega from the Binder cumulant data near criticality, as well as a study of the effective critical exponent eta and the effective amplitudes in the asymptotic expansion of susceptibility at the critical point. In all cases a refined analysis is consistent with our (GFD) asymptotic values of the critical exponents (nu=2/3, omega=1/2, eta=1/8 for 3D Ising model and omega=5/9 for 3D XY model), while the actually accepted "conventional" exponents are, in fact, effective exponents which are valid for approximation of the finite--size scaling behavior of not too large systems.