pith. machine review for the scientific record. sign in

arxiv: 2605.11573 · v1 · submitted 2026-05-12 · ✦ hep-th · cond-mat.stat-mech· hep-lat

Recognition: 2 theorem links

· Lean Theorem

A Monte Carlo Study of the Dipolar Universality Class in Three Dimensions

Akira Matsumoto , Yu Nakayama , Toshiki Onagi , Slava Rychkov
Authors on Pith no claims yet

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

classification ✦ hep-th cond-mat.stat-mechhep-lat
keywords dipolar universality classMonte Carlo simulationcritical exponentsphase transitionrotation invarianceBinder ratiothree dimensionsferromagnets
0
0 comments X

The pith

Monte Carlo simulations on a transversely constrained lattice model provide estimates of critical exponents for the three-dimensional dipolar universality class and demonstrate the emergence of rotation invariance at criticality.

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

The paper introduces a lattice model and Markov Chain Monte Carlo algorithm to simulate the phase transition in three-dimensional systems with strong dipolar interactions. This universality class was previously studied mainly through renormalization group methods and experiments, with little numerical lattice data available. The authors enforce the transverse constraint on the order parameter and use local Metropolis updates plus global zero-mode updates to sample configurations on lattices up to 48 cubed. They observe a continuous transition and compute universal quantities including critical exponents and the Binder ratio, comparing them to theoretical results. Additionally, they examine how full rotation invariance appears at the critical point despite the underlying anisotropy from dipolar forces.

Core claim

By implementing a lattice model that faithfully captures the transverse constraint on the magnetization and employing a hybrid Monte Carlo algorithm with local constraint-preserving moves and global updates, simulations reveal a continuous phase transition in the dipolar universality class. Estimates for the main critical exponents and Binder ratio are obtained and found to be consistent with renormalization group calculations, while the critical point exhibits restored rotational symmetry.

What carries the argument

A lattice model enforcing the transverse constraint on the order parameter combined with a Markov Chain Monte Carlo algorithm using local Metropolis updates and global zero-mode updates.

If this is right

  • The phase transition between disordered and ordered phases is continuous.
  • Critical exponents and the Binder ratio agree with renormalization group results.
  • Rotation invariance emerges at the critical point.
  • The algorithm enables access to universal quantities on cubic lattices up to volume 48 cubed.
  • The method bridges the gap between theory and direct numerical simulation for this class.

Where Pith is reading between the lines

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

  • The constrained update technique could apply to simulations of other systems with similar order-parameter restrictions.
  • These exponent estimates offer a benchmark for refining experimental probes of dipolar ferromagnets.
  • Larger-scale runs or alternative update schemes might tighten the numerical values further.
  • The observed symmetry restoration illustrates how fluctuations can isotropize anisotropic interactions near criticality in related models.

Load-bearing premise

The lattice model and update rules correctly implement the transverse constraint and properly sample the dipolar interaction ensemble without introducing artifacts.

What would settle it

A simulation result showing a first-order transition or critical exponents differing substantially from renormalization group predictions would falsify the model's validity for the dipolar class.

read the original abstract

The dipolar universality class describes the phase transition in 3D ferromagnets with strong dipolar interactions, as first discussed by Aharony and Fisher in the 1970s. While this universality class has been studied theoretically using renormalization group methods, as well as experimentally, little is known about it from Monte Carlo simulations. In this paper we aim to bridge this gap. We introduce a lattice model that faithfully implements the transverse constraint on the order parameter. We introduce a Markov Chain Monte Carlo algorithm which involves a combination of local Metropolis updates preserving the constraint, and a global update of the zero mode. We perform simulations on cubic lattices up to volume $48\times 48 \times 48$. We observe a continuous phase transition between the disordered and ordered phases. We obtain estimates of universal quantities such as the main critical exponents and the Binder ratio, and compare them with results from other techniques. We also investigate the emergence of rotation invariance at the critical point.

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 / 2 minor

Summary. The manuscript introduces a lattice model implementing the transverse constraint for the dipolar universality class in 3D ferromagnets. It presents a hybrid MCMC algorithm combining local Metropolis updates (which preserve the constraint) with global zero-mode updates, performs simulations on cubic lattices up to volume 48^3, reports a continuous phase transition, extracts estimates of critical exponents and the Binder ratio, compares these to renormalization-group and other results, and examines the emergence of rotational invariance at criticality.

Significance. If the sampling procedure is validated, the work supplies the first Monte Carlo estimates of universal quantities in this long-studied but numerically unexplored universality class, enabling quantitative confrontation with Aharony-Fisher renormalization-group predictions and with experimental data on dipolar magnets. The 48^3 volumes and the direct test of rotational symmetry restoration constitute concrete strengths that would strengthen the literature if the central numerical claims hold.

major comments (2)
  1. [Algorithm section (around the definition of the global update)] The description of the global zero-mode update states that it preserves the transverse constraint, yet no explicit demonstration of detailed balance for this move, no test of ergodicity on the constrained manifold, and no numerical diagnostics (e.g., comparison of cumulants on small volumes against exact enumeration or autocorrelation times) are supplied. Because every reported exponent, Binder ratio, and scaling observation rests on the assumption that the hybrid MCMC correctly samples the constrained ensemble, this verification is load-bearing for the central claims.
  2. [Results and comparison paragraphs] No error bars are quoted on the extracted critical exponents, no details of the finite-size scaling procedure (which observables, which fitting windows, how corrections-to-scaling are controlled) are given, and no quantitative comparison table with RG values or prior estimates appears. These omissions prevent assessment of the precision and consistency of the reported agreement with other techniques.
minor comments (2)
  1. [Model definition] The abstract states that the model 'faithfully implements the transverse constraint' but the precise lattice definition (how the constraint is enforced site-by-site) should be written explicitly with an equation.
  2. [Figures] Figure captions and axis labels for Binder-ratio crossings or scaling plots should include the precise definition of the Binder cumulant used and the range of volumes shown.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments on our manuscript. The points raised are important for strengthening the presentation and validation of our numerical results. We respond to each major comment below and indicate the revisions we will make.

read point-by-point responses

  1. Referee: [Algorithm section (around the definition of the global update)] The description of the global zero-mode update states that it preserves the transverse constraint, yet no explicit demonstration of detailed balance for this move, no test of ergodicity on the constrained manifold, and no numerical diagnostics (e.g., comparison of cumulants on small volumes against exact enumeration or autocorrelation times) are supplied. Because every reported exponent, Binder ratio, and scaling observation rests on the assumption that the hybrid MCMC correctly samples the constrained ensemble, this verification is load-bearing for the central claims.

    Authors: We agree that explicit verification of the sampling procedure is essential. In the revised manuscript we will expand the algorithm section with a detailed demonstration that the global zero-mode update satisfies detailed balance while preserving the transverse constraint. We will also address ergodicity of the hybrid algorithm on the constrained manifold and supply numerical diagnostics, including autocorrelation times for the primary observables together with comparisons of cumulants on small volumes against exact enumeration where feasible. These additions will directly support the reliability of the reported exponents and scaling results. revision: yes

  2. Referee: [Results and comparison paragraphs] No error bars are quoted on the extracted critical exponents, no details of the finite-size scaling procedure (which observables, which fitting windows, how corrections-to-scaling are controlled) are given, and no quantitative comparison table with RG values or prior estimates appears. These omissions prevent assessment of the precision and consistency of the reported agreement with other techniques.

    Authors: We acknowledge that the current version lacks error bars, a full description of the finite-size scaling analysis, and a quantitative comparison table. In the revision we will report statistical error bars on all critical exponents and the Binder ratio. We will describe the finite-size scaling procedure in detail, specifying the observables, fitting windows, and the treatment of corrections to scaling. We will also add a table providing direct numerical comparisons between our Monte Carlo estimates and renormalization-group predictions as well as other available results. These changes will enable a clearer evaluation of precision and consistency. revision: yes

Circularity Check

0 steps flagged

No circularity: results are direct numerical outputs from lattice Monte Carlo sampling

full rationale

The paper defines a lattice model implementing the transverse constraint, introduces a hybrid MCMC (local Metropolis + global zero-mode updates), runs it on volumes up to 48^3, and reports measured critical exponents, Binder ratio, and rotation invariance. These quantities are generated by the dynamics of the simulation rather than obtained by algebraic reduction, parameter fitting to the target observables, or load-bearing self-citation. No derivation chain exists that equates outputs to inputs by construction; comparisons to RG or experiment are external benchmarks, not internal tautologies. The sampling correctness is an assumption whose validity is independent of the reported numbers.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the correctness of the newly introduced lattice model and MCMC algorithm; no free parameters are fitted in the abstract, and no new entities are postulated.

axioms (1)
  • domain assumption Standard assumptions of statistical mechanics that a continuous phase transition belongs to a universality class whose exponents are independent of microscopic details once the symmetry and range of interactions are fixed.
    Invoked when the authors compare their numerical exponents to renormalization-group predictions for the dipolar class.

pith-pipeline@v0.9.0 · 5475 in / 1304 out tokens · 43699 ms · 2026-05-13T01:18:14.101655+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/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    We introduce a lattice model that faithfully implements the transverse constraint... combination of local Metropolis updates preserving the constraint, and a global update of the zero mode... cubic lattices up to volume 48×48×48... estimates of ... critical exponents ... Binder ratio ... emergence of rotation invariance

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

44 extracted references · 44 canonical work pages · 2 internal anchors

  1. [1]

    The Renormalization group and the epsilon expansion,

    K. Wilson and J. B. Kogut, “The Renormalization group and the epsilon expansion,” Phys.Rept.12(1974) 75–200

    work page 1974

  2. [2]

    Critical phenomena and renormalization group theory,

    A. Pelissetto and E. Vicari, “Critical phenomena and renormalization group theory,”Phys. Rept.368(2002) 549–727,arXiv:cond-mat/0012164 [cond-mat]

    work page arXiv 2002

  3. [3]

    Critical exponents and equation of state of the three-dimensional Heisenberg universality class,

    M. Campostrini, M. Hasenbusch, A. Pelissetto, P. Rossi, and E. Vicari, “Critical exponents and equation of state of the three-dimensional Heisenberg universality class,”Phys. Rev. B 65(2002) 144520,arXiv:cond-mat/0110336

    work page arXiv 2002

  4. [4]

    Precision Islands in the Ising and O(N) Models,

    F. Kos, D. Poland, D. Simmons-Duffin, and A. Vichi, “Precision Islands in the Ising and O(N) Models,”JHEP08(2016) 036,arXiv:1603.04436 [hep-th]

    work page arXiv 2016

  5. [5]

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

    M. Hasenbusch, “Monte Carlo study of a generalized icosahedral model on the simple cubic lattice,”Phys. Rev. B102no. 2, (2020) 024406,arXiv:2005.04448 [cond-mat.stat-mech]

    work page arXiv 2020

  6. [6]

    Bootstrapping Heisenberg magnets and their cubic instability,

    S. M. Chester, W. Landry, J. Liu, D. Poland, D. Simmons-Duffin, N. Su, and A. Vichi, “Bootstrapping Heisenberg magnets and their cubic instability,”Phys. Rev. D104no. 10, (2021) 105013,arXiv:2011.14647 [hep-th]

    work page arXiv 2021

  7. [7]

    ϕ 4 lattice model with cubic symmetry in three dimensions: Renormalization group flow and first-order phase transitions,

    M. Hasenbusch, “ϕ 4 lattice model with cubic symmetry in three dimensions: Renormalization group flow and first-order phase transitions,”Phys. Rev. B109no. 5, (2024) 054420,arXiv:2307.05165 [hep-lat]

    work page arXiv 2024

  8. [8]

    M. F. Collins,Magnetic Critical Scattering. Oxford University Press, Oxford, 1989

    work page 1989

  9. [9]

    Dependence of universal critical behaviour on symmetry and range of interaction,

    A. Aharony, “Dependence of universal critical behaviour on symmetry and range of interaction,” inPhase Transitions and Critical Phenomena, C. Domb and M. S. Green, eds., vol. 6, pp. 357–424. Academic Press, London and New York, 1976

    work page 1976

  10. [10]

    Dipolar Interactions at Ferromagnetic Critical Points,

    M. E. Fisher and A. Aharony, “Dipolar Interactions at Ferromagnetic Critical Points,”Phys. Rev. Lett.30(1973) 559–562

    work page 1973

  11. [11]

    Critical Behavior of Magnets with Dipolar Interactions. I. Renormalization Group near Four Dimensions,

    A. Aharony and M. E. Fisher, “Critical Behavior of Magnets with Dipolar Interactions. I. Renormalization Group near Four Dimensions,”Phys. Rev. B8(1973) 3323–3341

    work page 1973

  12. [12]

    Critical behavior of magnets with dipolar interactions. ii. feynman-graph expansion for ferromagnets near four dimensions,

    A. Aharony, “Critical behavior of magnets with dipolar interactions. ii. feynman-graph expansion for ferromagnets near four dimensions,”Phys. Rev. B8(1973) 3342–3348

    work page 1973

  13. [13]

    Scale without conformal invariance in dipolar ferromagnets,

    A. Gimenez-Grau, Y. Nakayama, and S. Rychkov, “Scale without conformal invariance in dipolar ferromagnets,”Phys. Rev. B110no. 2, (2024) 024421,arXiv:2309.02514 [hep-th]

    work page arXiv 2024

  14. [14]

    The Conformal Bootstrap: Theory, Numerical Techniques, and Applications

    D. Poland, S. Rychkov, and A. Vichi, “The Conformal Bootstrap: Theory, Numerical Techniques, and Applications,”Rev. Mod. Phys.91(2019) 015002,arXiv:1805.04405 [hep-th]. – 24 –

    work page Pith review arXiv 2019

  15. [15]

    Critical exponents of ferromagnets with dipolar interactions: Second-orderϵexpansion,

    A. D. Bruce and A. Aharony, “Critical exponents of ferromagnets with dipolar interactions: Second-orderϵexpansion,”Phys. Rev. B10(1974) 2078–2087

    work page 1974

  16. [16]

    Critical behavior of isotropic systems with strong dipole-dipole interaction: Three-loop study,

    A. Kudlis and A. Pikelner, “Critical behavior of isotropic systems with strong dipole-dipole interaction: Three-loop study,”Nuclear Physics B985(2022) 115990,arXiv:2204.02838 [cond-mat.stat-mech]

    work page arXiv 2022

  17. [17]

    Functional renormalization group approach to dipolar fixed point which is scale invariant but nonconformal,

    Y. Nakayama, “Functional renormalization group approach to dipolar fixed point which is scale invariant but nonconformal,”Phys. Rev. D110no. 2, (2024) 025020, arXiv:2309.15307 [hep-th]

    work page arXiv 2024

  18. [18]

    Critical behavior of isotropic systems with strong dipole-dipole interaction from the functional renormalization group

    G. Kalagov and N. Lebedev, “Critical behavior of isotropic systems with strong dipole-dipole interaction from the functional renormalization group,”arXiv:2602.04313 [cond-mat.stat-mech]

    work page internal anchor Pith review Pith/arXiv arXiv

  19. [19]

    Ferromagnetic Europium Compounds,

    T. R. McGuire and M. W. Shafer, “Ferromagnetic Europium Compounds,”J. Applied Phys. 35no. 3, (1964) 984–988

    work page 1964

  20. [20]

    Faraday rotation near the Curie point of EuO,

    C.-C. Huang and J. T. Ho, “Faraday rotation near the Curie point of EuO,”Phys. Rev. B12 (1975) 5255–5260

    work page 1975

  21. [21]

    Critical order-parameter relaxation in EuS, a nearly isotropic ferromagnet,

    J. Kotzler, G. Kamleiter, and G. Weber, “Critical order-parameter relaxation in EuS, a nearly isotropic ferromagnet,”J. Phys. C9no. 13, (1976) L361

    work page 1976

  22. [22]

    Depression of Longitudinal Fluctuations aboveT c of Heisenberg Ferromagnets Observed by Polarized Neutrons,

    J. K¨ otzler, D. G¨ orlitz, F. Mezei, and B. Farago, “Depression of Longitudinal Fluctuations aboveT c of Heisenberg Ferromagnets Observed by Polarized Neutrons,”Europhys. Lett.1 no. 12, (1986) 675

    work page 1986

  23. [23]

    Monte Carlo study of the Heisenberg model with local dipolar interaction,

    E. Itou, A. Matsumoto, Y. Nakayama, and T. Onagi, “Monte Carlo study of the Heisenberg model with local dipolar interaction,”Phys. Rev. B112no. 13, (2025) 134409, arXiv:2503.15874 [hep-th]

    work page arXiv 2025

  24. [24]

    The “Coulomb Phase

    C. L. Henley, “The “Coulomb Phase” in Frustrated Systems,”Annual Review of Condensed Matter Physics1(Apr., 2010) 179–210,arXiv:0912.4531 [cond-mat.stat-mech]

    work page arXiv 2010

  25. [25]

    Universal conductivity of dirty bosons at the superconductor-insulator transition,

    E. S. Sørensen, M. Wallin, S. M. Girvin, and A. P. Young, “Universal conductivity of dirty bosons at the superconductor-insulator transition,”Phys. Rev. Lett.69(Aug, 1992) 828–831

    work page 1992

  26. [26]

    Exchange Monte Carlo Method and Application to Spin Glass Simulations

    K. Hukushima and K. Nemoto, “Exchange Monte Carlo Method and Application to Spin Glass Simulations,”J. Phys. Soc. Jap.65no. 6, (1996) 1604–1608, arXiv:cond-mat/9512035

    work page Pith review arXiv 1996

  27. [27]

    Parallel tempering and decorrelation of topological charge in full QCD,

    E.-M. Ilgenfritz, W. Kerler, M. Muller-Preussker, and H. Stuben, “Parallel tempering and decorrelation of topological charge in full QCD,”arXiv:hep-lat/0007039

    work page internal anchor Pith review arXiv

  28. [28]

    Cluster Monte Carlo algorithm for the quantum rotor model,

    F. Alet and E. S. Sørensen, “Cluster Monte Carlo algorithm for the quantum rotor model,” Phys.Rev.E67no. 1, (Jan., 2003) 015701,arXiv:cond-mat/0211262 [cond-mat.str-el]

    work page arXiv 2003

  29. [29]

    Directed geometrical worm algorithm applied to the quantum rotor model,

    F. Alet and E. S. Sørensen, “Directed geometrical worm algorithm applied to the quantum rotor model,”Phys.Rev.E68no. 2, (Aug., 2003) 026702,arXiv:cond-mat/0303080 [cond-mat.str-el]

    work page arXiv 2003

  30. [30]

    Finite size scaling analysis of Ising model block distribution functions,

    K. Binder, “Finite size scaling analysis of Ising model block distribution functions,”Z. Phys. B43(1981) 119–140

    work page 1981

  31. [31]

    A measure of data collapse for scaling,

    S. M. Bhattacharjee and F. Seno, “A measure of data collapse for scaling,”Journal of Physics A: Mathematical and General34no. 33, (Aug., 2001) 6375–6380. – 25 –

    work page 2001

  32. [32]

    Critical behavior of the three-dimensional±j model in a magnetic field,

    N. Kawashima and N. Ito, “Critical behavior of the three-dimensional±j model in a magnetic field,”Journal of the Physical Society of Japan62no. 2, (1993) 435–438, https://doi.org/10.1143/JPSJ.62.435.https://doi.org/10.1143/JPSJ.62.435

    work page doi:10.1143/jpsj.62.435.https://doi.org/10.1143/jpsj.62.435 1993

  33. [33]

    Low-temperature behavior of two-dimensional gaussian ising spin glasses,

    J. Houdayer and A. Hartmann, “Low-temperature behavior of two-dimensional gaussian ising spin glasses,”Physical Review B70no. 1, (July, 2004) . http://dx.doi.org/10.1103/PhysRevB.70.014418

    work page doi:10.1103/physrevb.70.014418 2004

  34. [34]

    Critical Behavior of Anisotropic Cubic Systems,

    A. Aharony, “Critical Behavior of Anisotropic Cubic Systems,”Phys. Rev. B8(1973) 4270–4273

    work page 1973

  35. [35]

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

    M. Hasenbusch, “Cubic fixed point in three dimensions: Monte Carlo simulations of theϕ4 model on the simple cubic lattice,”Phys. Rev. B107no. 2, (2023) 024409, arXiv:2211.16170 [cond-mat.stat-mech]

    work page arXiv 2023

  36. [36]

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

    M. Hasenbusch, “Finite size scaling study of lattice models in the three-dimensional Ising universality class,”Phys. Rev. B82(2010) 174433,arXiv:1004.4486 [cond-mat.stat-mech]

    work page arXiv 2010

  37. [37]

    Two- and three-point functions at criticality: Monte Carlo simulations of the improved three-dimensional Blume-Capel model,

    M. Hasenbusch, “Two- and three-point functions at criticality: Monte Carlo simulations of the improved three-dimensional Blume-Capel model,”Phys. Rev. E97no. 1, (2018) 012119, arXiv:1711.10946 [hep-lat]

    work page arXiv 2018

  38. [38]

    Dynamic critical exponent z of the three-dimensional Ising universality class: Monte Carlo simulations of the improved Blume-Capel model,

    M. Hasenbusch, “Dynamic critical exponent z of the three-dimensional Ising universality class: Monte Carlo simulations of the improved Blume-Capel model,”Phys. Rev. E101 no. 2, (2020) 022126,arXiv:1908.01702 [cond-mat.stat-mech]

    work page arXiv 2020

  39. [39]

    Monte Carlo study of an improved clock model in three dimensions,

    M. Hasenbusch, “Monte Carlo study of an improved clock model in three dimensions,”Phys. Rev. B100no. 22, (2019) 224517,arXiv:1910.05916 [cond-mat.stat-mech]

    work page arXiv 2019

  40. [40]

    Eliminating leading and subleading corrections to scaling in the three-dimensional XY universality class,

    M. Hasenbusch, “Eliminating leading and subleading corrections to scaling in the three-dimensional XY universality class,”Phys. Rev. B112no. 18, (2025) 184512, arXiv:2507.19265 [cond-mat.stat-mech]

    work page arXiv 2025

  41. [41]

    Neutron scattering from the heisenberg ferromagnets euo and eus. ii. static critical properties,

    J. Als-Nielsen, O. W. Dietrich, and L. Passell, “Neutron scattering from the heisenberg ferromagnets euo and eus. ii. static critical properties,”Phys. Rev. B14(1976) 4908–4922

    work page 1976

  42. [42]

    Critical magnetic properties and exchange interactions in euo,

    N. Menyuk, K. Dwight, and T. B. Reed, “Critical magnetic properties and exchange interactions in euo,”Phys. Rev. B3(1971) 1689–1698

    work page 1971

  43. [43]

    Heat capacity of euo near the curie temperature,

    A. Kornblit and G. Ahlers, “Heat capacity of euo near the curie temperature,”Phys. Rev. B 11(1975) 2678–2688

    work page 1975

  44. [44]

    Heat capacity of eus near the ferromagnetic transition,

    A. Kornblit, G. Ahlers, and E. Buehler, “Heat capacity of eus near the ferromagnetic transition,”Phys. Rev. B17(1978) 282–292. – 26 –

    work page 1978