pith. sign in

arxiv: 2509.05195 · v2 · submitted 2025-09-05 · ❄️ cond-mat.stat-mech

Constraint effective action and critical correlation functions at fixed magnetization

Pith reviewed 2026-05-18 18:36 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech
keywords functional renormalization groupconstraint effective actionIsing universality classcritical correlation functionsderivative expansionfixed magnetizationrate function
0
0 comments X

The pith

The functional renormalization group at second order in the derivative expansion computes the constraint effective action and momentum-dependent correlation functions at fixed magnetization for the Ising universality class.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper extends the functional renormalization group framework from computing critical probability distributions of the order parameter to also obtaining momentum-dependent observables. It focuses on the constraint effective action at fixed magnetization in the Ising universality class, derives the exact flow equations, and solves them numerically at the second order of the derivative expansion. The solutions yield universal rate functions and correlation functions that are benchmarked against Monte Carlo simulations. In three dimensions the results recover the rate function and reproduce the first few Fourier modes accurately, while in two dimensions the second-order approximation succeeds where the local potential approximation fails.

Core claim

The exact flow equations for the constraint effective action at fixed magnetization are derived and solved numerically within the second-order derivative expansion, producing universal rate functions and momentum-dependent correlation functions that agree with Monte Carlo simulations in three dimensions and show qualitative agreement in two dimensions.

What carries the argument

The constraint effective action at fixed magnetization, with its exact renormalization group flow equations solved at second order of the derivative expansion to extract rate functions and Fourier modes of the correlation function.

If this is right

  • In three dimensions the rate function is recovered and the first few Fourier modes of the constrained correlation function are accurately reproduced.
  • The numerical solutions demonstrate convergence of the method.
  • In two dimensions at least the second-order derivative expansion is required to describe the critical point.
  • The results confirm the robustness of the approach for both zero- and finite-momentum critical observables at fixed magnetization.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same flow-equation setup could be used to extract higher-order correlation functions at fixed magnetization.
  • The framework may extend directly to other O(N) universality classes by changing the symmetry of the effective action.
  • Finite-size scaling corrections to the rate function could be incorporated by adjusting the boundary conditions in the numerical integration.

Load-bearing premise

The second-order derivative expansion truncation remains accurate enough to describe the critical point, including in two dimensions where the local potential approximation fails.

What would settle it

High-precision Monte Carlo data for the first few Fourier modes of the constrained correlation function in three dimensions that disagree with the numerically extracted modes beyond numerical error bars.

Figures

Figures reproduced from arXiv: 2509.05195 by Adam Ran\c{c}on, F\'elix Rose, Ivan Balog.

Figure 1
Figure 1. Figure 1: FIG. 1. Sketches of the small-momentum behavior of the vertices [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Main panel: critical rate function [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Main panel: critical rate function [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Critical rate function [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Critical rate function [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Function [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Function [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
read the original abstract

We present an extension of the functional renormalization group (FRG) framework developed to compute critical probability distributions of the order parameter to momentum-dependent observables. Focusing on the constraint effective action at fixed magnetization for the Ising universality class, we derive its exact flow equations and solve them at the second order of the derivative expansion (DE2). We solve these flow equations numerically for two- and three-dimensional systems, extract universal rate functions and momentum-dependent correlation functions, and benchmark them against Monte Carlo simulations. In three dimensions, we recover the rate function and accurately reproduce the first few Fourier modes of the constrained correlation function and demonstrate the convergence of the method. In two dimensions, the lowest order approximations such as local potential approximation (LPA) fail, and it is required to consider at least the DE2 to describe the critical point. Our results are in qualitative agreement with the numerics. We confirm the robustness of the FRG approach for calculating both zero- and finite-momentum critical observables at fixed magnetization.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript extends the functional renormalization group (FRG) framework to momentum-dependent observables by deriving exact flow equations for the constraint effective action at fixed magnetization in the Ising universality class. These are solved numerically at second order in the derivative expansion (DE2) for two- and three-dimensional systems. Universal rate functions and momentum-dependent correlation functions are extracted and benchmarked against Monte Carlo simulations, with accurate reproduction of the rate function and first few Fourier modes in 3D and qualitative agreement in 2D where LPA fails.

Significance. If the central numerical claims hold, the work provides a systematic FRG route to both zero- and finite-momentum critical observables under fixed-magnetization constraints, extending beyond local-potential approximations. The derivation of exact flow equations within the FRG and direct benchmarking against independent Monte Carlo data are clear strengths that could aid studies of constrained critical phenomena.

major comments (2)
  1. [Two-dimensional results section] Two-dimensional results: The claim that DE2 is required and sufficient because LPA fails rests on qualitative agreement with Monte Carlo for the rate function and constrained correlation functions. No quantitative assessment of truncation error is provided via comparison to DE3 or to exact 2D Ising correlation functions, leaving open whether the regulator-dependent integrals have converged for finite-momentum observables.
  2. [Numerical implementation and results] Numerical solution and regulator dependence: The flow equations for the momentum-dependent parts involve regulator-dependent integrals whose specific form and cutoff sensitivity are not detailed with convergence checks or error estimates. This affects the robustness of the reported Fourier modes in 3D and the overall demonstration of method convergence.
minor comments (1)
  1. [Abstract] The abstract states that the method demonstrates 'convergence' in 3D; this should be defined explicitly (e.g., with respect to regulator variation or truncation order) in the main text.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments. We address the two major comments point by point below.

read point-by-point responses
  1. Referee: [Two-dimensional results section] Two-dimensional results: The claim that DE2 is required and sufficient because LPA fails rests on qualitative agreement with Monte Carlo for the rate function and constrained correlation functions. No quantitative assessment of truncation error is provided via comparison to DE3 or to exact 2D Ising correlation functions, leaving open whether the regulator-dependent integrals have converged for finite-momentum observables.

    Authors: We agree that a direct comparison to a DE3 truncation would allow a more quantitative estimate of truncation errors. However, the flow equations at DE3 become substantially more involved, with additional momentum-dependent vertices and a significantly higher computational cost for the numerical solution. Our present work is focused on establishing that the DE2 truncation is already sufficient to capture the qualitative features where LPA fails, as demonstrated by the improvement over LPA and the agreement with Monte Carlo data. Exact analytic expressions for the momentum-dependent constrained correlation functions of the 2D Ising model at fixed magnetization are not available in closed form for direct comparison. Monte Carlo simulations therefore remain the most practical benchmark. In a revised manuscript we will add an explicit paragraph discussing the truncation level, the practical limitations of extending to DE3, and the role of the Monte Carlo comparison. revision: partial

  2. Referee: [Numerical implementation and results] Numerical solution and regulator dependence: The flow equations for the momentum-dependent parts involve regulator-dependent integrals whose specific form and cutoff sensitivity are not detailed with convergence checks or error estimates. This affects the robustness of the reported Fourier modes in 3D and the overall demonstration of method convergence.

    Authors: We thank the referee for this observation. The manuscript states the regulator choice and outlines the numerical procedure, yet we acknowledge that the explicit form of the momentum integrals and systematic convergence tests are not presented in sufficient detail. In the revised version we will add a short appendix (or extended methods section) that (i) writes out the regulator-dependent integrals appearing in the flow equations for the momentum-dependent vertices, (ii) describes the discretization and quadrature scheme used to evaluate them, and (iii) reports convergence checks obtained by varying the ultraviolet cutoff scale, the infrared regulator parameter, and the momentum-grid resolution. These checks will be shown to leave the first few Fourier modes in three dimensions stable within the reported precision, thereby strengthening the demonstration of numerical robustness. revision: yes

Circularity Check

0 steps flagged

FRG derivation and DE2 solution are self-contained with external Monte Carlo benchmarks

full rationale

The paper states it derives exact flow equations for the constraint effective action within the standard FRG framework, then truncates to DE2 and solves numerically for rate functions and momentum-dependent correlations in 2D and 3D Ising class. These outputs are benchmarked against independent Monte Carlo data rather than fitted to the same observables or derived by construction from inputs. No self-definitional reductions, fitted parameters renamed as predictions, or load-bearing self-citations appear in the provided derivation chain; the central results rest on the FRG equations plus numerical solution plus external validation.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the standard FRG flow equation and the assumption that the DE2 truncation captures the essential physics at criticality. No new particles or forces are introduced.

axioms (2)
  • standard math The Wetterich equation provides the exact flow for the effective average action
    Invoked to derive the flow equations for the constraint effective action as stated in the abstract.
  • domain assumption The derivative expansion truncated at second order is sufficient to describe the critical fixed point
    Explicitly required for the numerical solution in both 2D and 3D; the abstract notes that LPA fails in 2D.

pith-pipeline@v0.9.0 · 5705 in / 1498 out tokens · 54894 ms · 2026-05-18T18:36:45.042597+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Cost/FunctionalEquation washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    We solve these flow equations numerically for two- and three-dimensional systems, extract universal rate functions and momentum-dependent correlation functions... In two dimensions, the lowest order approximations such as local potential approximation (LPA) fail, and it is required to consider at least the DE2

  • IndisputableMonolith/Foundation/AlexanderDuality alexander_duality_circle_linking unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    the two-point vertex in constant field is expected to take the form Γ̂(2)ₖ(p;s)=I''ₖ(s)+(1−δ_{p,0})Δ̂₂,ₖ(s)+Ẑₖ(s)p²

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

58 extracted references · 58 canonical work pages

  1. [1]

    All FRG computations in three dimensions use the exponential regulator, with optimal prefactorαopt≈4.65at LPA andαopt ≈1.30at DE2

    We therefore define the most conservative error bar at DE2 order to be one-quarter of the difference between the LPA and DE2 rate function curves. All FRG computations in three dimensions use the exponential regulator, with optimal prefactorαopt≈4.65at LPA andαopt ≈1.30at DE2. At these parameters, the DE2 approximation yieldsη≈0.0455andν≈0.6275, compared ...

  2. [2]

    PHC COGITO

    IB wishes to acknowledge the support of the INFaR and FrustKor projects financed by the EU through the National Recovery and Resilience Plan (NRRP) 2021- 2026. AR has benefited from the financial support of the Grant No. ANR-24-CE30-6695 FUSIoN. This work was partially supported by an IEA CNRS project and by the “PHC COGITO” program (project number: 49149...

  3. [3]

    Bouchaud and A

    J.-P. Bouchaud and A. Georges, Anomalous diffusion in disordered media: Statistical mechanisms, models and physical applications, Phys. Rep.195, 127 (1990)

  4. [4]

    A. D. Bruce, T. Schneider, and E. Stoll, Cluster Phenomena and the Displacive to Order-Disorder Crossover, Phys. Rev. Lett.43, 1284 (1979)

  5. [5]

    Binder, Critical Properties from Monte Carlo Coarse Graining and Renormalization, Phys

    K. Binder, Critical Properties from Monte Carlo Coarse Graining and Renormalization, Phys. Rev. Lett.47, 693 (1981)

  6. [6]

    Binder, Finite Size Scaling Analysis of Ising Model Block Distribution Functions, Z

    K. Binder, Finite Size Scaling Analysis of Ising Model Block Distribution Functions, Z. Phys. B43, 119 (1981)

  7. [7]

    Rudnick, H

    J. Rudnick, H. Guo, and D. Jasnow, Finite-size scaling and the renormalization group, J. Stat. Phys.41, 353 (1985)

  8. [8]

    Eisenriegler and R

    E. Eisenriegler and R. Tomaschitz, Helmholtz free energy of finite spin systems near criticality, Phys. Rev. B35, 4876 (1987)

  9. [9]

    Kuti and Y

    J. Kuti and Y. Shen, Supercomputing the Effective Action, Phys. Rev. Lett.60, 85 (1988)

  10. [10]

    Nicolaides and A

    D. Nicolaides and A. D. Bruce, Universal configurational structure in two-dimensional scalar models, J. Phys. A21, 233 (1988)

  11. [11]

    Brankov and D

    J. Brankov and D. Danchev, A probabilistic view on finite-size scaling in infinitely coordinated spherical models, Physica A158, 842 (1989)

  12. [12]

    A. D. Bruce and N. B. Wilding, Scaling fields and universality of the liquid-gas critical point, Phys. Rev. Lett.68, 193 (1992)

  13. [13]

    M. M. Tsypin, Universal Effective Potential for Scalar Field Theory in Three Dimensions by Monte Carlo Computation, Phys. Rev. Lett.73, 2015 (1994)

  14. [14]

    Hilfer and N

    R. Hilfer and N. B. Wilding, Are critical finite-size scaling functions calculable from knowledge of an appropriate critical exponent?, J. Phys. A28, L281 (1995)

  15. [15]

    Esser, V

    A. Esser, V. Dohm, and X. Chen, Field theory of finite-size effects for systems with a one-component order parameter, Physica A222, 355 (1995)

  16. [16]

    X. S. Chen, V. Dohm, and N. Schultka, Order-Parameter Distribution Function of FiniteO(n)Symmetric Systems, Phys. Rev. Lett.77, 3641 (1996)

  17. [17]

    A. D. Bruce, Finite-size critical behavior in the Gibbs ensemble, Phys. Rev. E55, 2315 (1997)

  18. [18]

    Rudnick, W

    J. Rudnick, W. Lay, and D. Jasnow, Effective potential, critical point scaling, and the renormalization group, Phys. Rev. E58, 2902 (1998)

  19. [19]

    M. M. Tsypin and H. W. J. Blöte, Probability distribution of the order parameter for the three-dimensional Ising-model universality class: A high-precision Monte Carlo study, Phys. Rev. E62, 73 (2000)

  20. [20]

    Portelli and P

    B. Portelli and P. C. W. Holdsworth, Relevance of soft modes for order parameter fluctuations in the two-dimensional XY model, J. Phys. A35, 1231 (2002)

  21. [21]

    Hilfer, B

    R. Hilfer, B. Biswal, H. G. Mattutis, and W. Janke, Multicanonical Monte Carlo study and analysis of tails for the order-parameter distribution of the two-dimensional Ising model, Phys. Rev. E68, 46123 (2003)

  22. [22]

    Hilfer, B

    R. Hilfer, B. Biswal, H.-G. Mattutis, and W. Janke, Multicanonical simulations of the tails of the order-parameter distri- bution of the two-dimensional Ising model, Comput. Phys. Commun.169, 230 (2005)

  23. [23]

    Malakis, N

    A. Malakis, N. G. Fytas, and G. Gülpinar, Critical Binder cumulant and universality: Fortuin-Kasteleyn clusters and order-parameter fluctuations, Phys. Rev. E89, 042103 (2014)

  24. [24]

    J. Xu, A. M. Ferrenberg, and D. P. Landau, High-resolution Monte Carlo study of the order-parameter distribution of the three-dimensional Ising model, Phys. Rev. E101, 023315 (2020)

  25. [25]

    Balog, A

    I. Balog, A. Rançon, and B. Delamotte, Critical Probability Distributions of the Order Parameter from the Functional Renormalization Group, Phys. Rev. Lett.129, 210602 (2022)

  26. [26]

    Rançon, B

    A. Rançon, B. Delamotte, L. Šaravanja, and I. Balog, Probability distributions of the order parameter of the O(N) model, Phys. Rev. E111, 034131 (2025)

  27. [27]

    S. Sahu, B. Delamotte, and A. Rançon, Generalization of the central limit theorem to critical systems: Revisiting pertur- bation theory, Phys. Rev. E111, 034128 (2025)

  28. [28]

    Balog, B

    I. Balog, B. Delamotte, and A. Rancon, Universal and non-universal large deviations in critical systems, SciPost Physics 18, 119 (2025)

  29. [29]

    Sahu, Probability distributions of the order parameter of the Ising Model at two loops, arXiv:2501.18615 (2025)

    S. Sahu, Probability distributions of the order parameter of the Ising Model at two loops, arXiv:2501.18615 (2025). 17

  30. [30]

    Sahu, Critical probability distributions of the order parameter at two loops II:O(n)universality class, arXiv:2503.21799 (2025)

    S. Sahu, Critical probability distributions of the order parameter at two loops II:O(n)universality class, arXiv:2503.21799 (2025)

  31. [31]

    S. T. Bramwell, P. C. W. Holdsworth, and J. F. Pinton, Universality of rare fluctuations in turbulence and critical phenomena, Nature (London)396, 552 (1998)

  32. [32]

    Joubaud, A

    S. Joubaud, A. Petrosyan, S. Ciliberto, and N. B. Garnier, Experimental Evidence of Non-Gaussian Fluctuations near a Critical Point, Phys. Rev. Lett.100, 180601 (2008)

  33. [33]

    K. A. Takeuchi and M. Sano, Universal Fluctuations of Growing Interfaces: Evidence in Turbulent Liquid Crystals, Phys. Rev. Lett.104, 230601 (2010)

  34. [34]

    Allemand, G

    M. Allemand, G. Dupuy, P. Paquiez, N. Dupuis, A. Rançon, T. Roscilde, T. Chalopin, and D. Clément, Observation of universal non-Gaussian statistics of the order parameter across a continuous phase transition, arXiv:2508.21623 (2025)

  35. [35]

    Jona-Lasinio, The renormalization group: A probabilistic view, Il Nuovo Cimento B (1971-1996)26, 99 (1975)

    G. Jona-Lasinio, The renormalization group: A probabilistic view, Il Nuovo Cimento B (1971-1996)26, 99 (1975)

  36. [36]

    Gallavotti and H

    G. Gallavotti and H. J. F. Knops, Block-spins interactions in the Ising model, Comm. Math. Phys.36, 171 (1974)

  37. [37]

    Gallavotti and A

    G. Gallavotti and A. Martin-Löf, Block-spin distributions for short-range attractive Ising models, Il Nuovo Cimento B (1971-1996)25, 425 (1975)

  38. [38]

    Dupuis, L

    N. Dupuis, L. Canet, A. Eichhorn, W. Metzner, J. Pawlowski, M. Tissier, and N. Wschebor, The nonperturbative functional renormalization group and its applications, Phys. Rep.910, 1 (2021)

  39. [39]

    Fukuda and E

    R. Fukuda and E. Kyriakopoulos, Derivation of the effective potential, Nucl. Phys. B85, 354 (1975)

  40. [40]

    O’Raifeartaigh, A

    L. O’Raifeartaigh, A. Wipf, and H. Yoneyama, The constraint effective potential, Nucl. Phys. B271, 653 (1986)

  41. [41]

    Hasenfratz, K

    A. Hasenfratz, K. Jansen, J. Jersák, C. B. Lang, T. Neuhaus, and H. Yoneyama, Study of the 4-componentϕ4 model, Nucl. Phys. B317, 81 (1989)

  42. [42]

    Göckeler and H

    M. Göckeler and H. Leutwyler, Shape of the constraint effective potential, Nucl. Phys. B350, 228 (1991)

  43. [43]

    Göckeler and H

    M. Göckeler and H. Leutwyler, Constraint correlation functions in the O(N) model, Nucl. Phys. B361, 392 (1991)

  44. [44]

    Mukaida and Y

    H. Mukaida and Y. Shimada, Convex effective potential ofO(N)-symmetricϕ4 theory for largeN, Nucl. Phys. B479, 663 (1996)

  45. [45]

    Camia, C

    F. Camia, C. Garban, and C. M. Newman, Planar Ising magnetization field II. Properties of the critical and near-critical scaling limits, Ann. inst. Henri Poincare (B) Probab. Stat.52, 146 (2016)

  46. [46]

    Blaizot, R

    J.-P. Blaizot, R. Méndez-Galain, and N. Wschebor, A new method to solve the non-perturbative renormalization group equations, Phys. Lett. B632, 571 (2006)

  47. [47]

    Dupuis,Field Theory of Condensed Matter and Ultracold Gases(World Scientific (Europe), 2023)

    N. Dupuis,Field Theory of Condensed Matter and Ultracold Gases(World Scientific (Europe), 2023)

  48. [48]

    Rançon, Nonperturbative renormalization group approach to quantum XY spin models, Phys

    A. Rançon, Nonperturbative renormalization group approach to quantum XY spin models, Phys. Rev. B89, 214418 (2014)

  49. [49]

    Benitez, J.-P

    F. Benitez, J.-P. Blaizot, H. Chaté, B. Delamotte, R. Méndez-Galain, and N. Wschebor, Nonperturbative renormalization group preserving full-momentum dependence: Implementation and quantitative evaluation, Phys. Rev. E85, 026707 (2012)

  50. [50]

    Bervillier, Universal relations among critical amplitude

    C. Bervillier, Universal relations among critical amplitude. Calculations up to orderϵ2 for systems with continuous sym- metry, Phys. Rev. B14, 4964 (1976)

  51. [51]

    Privman and M

    V. Privman and M. E. Fisher, Universal critical amplitudes in finite-size scaling, Phys. Rev. B30, 322 (1984)

  52. [52]

    Balog, H

    I. Balog, H. Chaté, B. Delamotte, M. Marohnić, and N. Wschebor, Convergence of Nonperturbative Approximations to the Renormalization Group, Phys. Rev. Lett.123, 240604 (2019)

  53. [53]

    De Polsi, I

    G. De Polsi, I. Balog, M. Tissier, and N. Wschebor, Precision calculation of critical exponents in the O(N) universality classes with the nonperturbative renormalization group, Phys. Rev. E101, 42113 (2020)

  54. [54]

    F. Kos, D. Poland, D. Simmons-Duffin, and A. Vichi, Precision islands in the ising andO(N)models, JHEP2016(8), 36

  55. [55]

    R. H. Swendsen and J.-S. Wang, Nonuniversal critical dynamics in monte carlo simulations, Phys. Rev. Lett.58, 86 (1987)

  56. [56]

    Newman and G

    M. Newman and G. Barkema,Monte Carlo methods in statistical physics(Clarendon Press, Oxford : New York, 1999)

  57. [57]

    Balog and A

    I. Balog and A. Rançon, to be published

  58. [58]

    Canet, B

    L. Canet, B. Delamotte, D. Mouhanna, and J. Vidal, Optimization of the derivative expansion in the nonperturbative renormalization group, Phys. Rev. D67, 065004 (2003)