pith. machine review for the scientific record. sign in

arxiv: 2604.08777 · v1 · submitted 2026-04-09 · ❄️ cond-mat.mtrl-sci · physics.comp-ph

Recognition: unknown

Including sample shape in micromagnetics with 3D periodic boundary conditions

Frederik Laust Durhuus , Andrea Roberto Insinga , Rasmus Bj{\o}rk
Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:42 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci physics.comp-ph
keywords micromagneticsperiodic boundary conditionsshape effectsdemagnetizing fieldaverage magnetizationcomputational methodsmagnetic simulationsquasi-periodic macrogeometry
0
0 comments X

The pith

For large magnetic samples, only average magnetization produces non-negligible shape effects in periodic simulations.

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

The paper establishes a formal proof that shape-dependent magnetic fields in periodic micromagnetic simulations depend, for sufficiently large samples, solely on the average magnetization rather than on its spatial variations. This result addresses the practical problem that standard quasi-periodic macrogeometry methods require many replicated copies of the simulated domain, which becomes inefficient or impossible when the domain is incommensurate with the target sample shape. By isolating the average-magnetization contribution, the authors derive a minimal correction that existing PBC implementations can adopt with negligible added cost. The approach therefore lets researchers include realistic finite-sample demagnetizing effects while retaining the computational advantages of periodic boundary conditions.

Core claim

The central claim is that, for sufficiently large magnetic samples, the shape contribution to the effective field is exhausted by the term arising from the spatially averaged magnetization; all higher-order spatial fluctuations of the magnetization produce vanishingly small corrections to the demagnetizing field. This mathematical reduction justifies a direct, low-cost modification to existing quasi-periodic PBC codes that adds the missing shape term without enlarging the computational cell.

What carries the argument

The formal proof that shape effects reduce to the average magnetization for large samples, which supplies the correction term inserted into standard PBC field calculations.

If this is right

  • Standard PBC micromagnetic codes can be updated with a single average-magnetization correction to recover finite-sample shape effects.
  • Domains whose periodicity does not match the desired overall shape can now be simulated without prohibitive replication overhead.
  • Magnetization patterns with strong spatial variation can be treated periodically while still capturing the dominant demagnetizing field from the sample boundary.
  • The method preserves the speed advantage of PBC while removing the systematic error that arises when shape is ignored.

Where Pith is reading between the lines

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

  • The same reduction argument may apply to other periodic physical fields, such as electrostatic potentials in repeating dielectric structures.
  • Empirical calibration of the minimal sample size for a given material and discretization would turn the 'sufficiently large' criterion into a practical rule of thumb.
  • Hybrid schemes could combine the corrected PBC interior with a boundary-element treatment of the outermost shape.

Load-bearing premise

The sample must be large enough that higher-order spatial variations of the magnetization make negligible contributions to the overall shape effects.

What would settle it

A numerical test comparing the corrected PBC method against a full macrogeometry simulation (many replicated cells) for the same magnetization pattern while systematically increasing sample size until the two results converge within a chosen error tolerance.

Figures

Figures reproduced from arXiv: 2604.08777 by Andrea Roberto Insinga, Frederik Laust Durhuus, Rasmus Bj{\o}rk.

Figure 1
Figure 1. Figure 1: FIG. 1: Illustration of micromagnetic problem in 3D [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Illustration of the simulated system, i.e. a [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Convergence rate vs. number of macrogeometry [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Normalised magnetisation vs. applied field [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Coercive field vs. aspect ratio and intergrain [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
read the original abstract

Periodic boundary conditions (PBCs) for computing magnetic fields in repeating magnetic structures, e.g. in micromagnetic simulations, are typically imposed using the quasi periodic macrogeometry approach, where many copies of the simulated domain are introduced. This can be computationally problematic, especially if the simulated domain is incommensurate with the desired sample shape. In this work, we present a formal proof that for sufficiently large magnetic samples, only the average magnetisation gives non-negligible shape effects. Using this insight, we develop a simple, computationally efficient modification of existing implementations which incorporates shape effects in PBC methods.

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 paper claims that for sufficiently large magnetic samples, shape effects in the magnetostatic field reduce to a contribution determined solely by the spatially averaged magnetization, with higher-order spatial variations of M becoming negligible. It uses this to derive a simple, computationally efficient correction to standard 3D periodic-boundary-condition (PBC) implementations that incorporates finite-sample shape effects without needing an incommensurate quasi-periodic macrogeometry.

Significance. If the central approximation holds with quantifiable error bounds, the method would allow accurate micromagnetic modeling of finite-shaped periodic structures at far lower cost than full quasi-periodic copies, addressing a common practical limitation in simulations of real devices. The formal proof, if rigorous and accompanied by validation, would provide a stronger foundation than existing heuristic corrections and could be adopted in standard micromagnetic codes.

major comments (2)
  1. [Abstract] Abstract: the central claim of a 'formal proof' that only the average magnetisation produces non-negligible shape effects for 'sufficiently large' samples is load-bearing for the subsequent PBC modification, yet the abstract (and by extension the manuscript) supplies neither an explicit length-scale criterion (e.g., sample size versus exchange length or magnetostatic length) nor an a-priori error bound when the condition is violated.
  2. [Method / Results] The proposed modification discards all higher-order spatial variations of the magnetization when computing the shape correction; without numerical benchmarks comparing the corrected PBC result against a full finite-shape calculation for samples near the 'sufficiently large' threshold, it is impossible to assess the practical range of validity.
minor comments (2)
  1. [Introduction] Notation for the average magnetization and the shape-correction term should be introduced with explicit definitions and units in the first appearance.
  2. [Results] The manuscript would benefit from a short table or figure quantifying the computational speedup relative to standard quasi-periodic PBC for representative sample sizes.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed review. The comments highlight important points regarding the presentation of our central approximation and the need for practical validation. We address each major comment below and will incorporate revisions to strengthen the manuscript.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim of a 'formal proof' that only the average magnetisation produces non-negligible shape effects for 'sufficiently large' samples is load-bearing for the subsequent PBC modification, yet the abstract (and by extension the manuscript) supplies neither an explicit length-scale criterion (e.g., sample size versus exchange length or magnetostatic length) nor an a-priori error bound when the condition is violated.

    Authors: We agree that the abstract would benefit from greater specificity on the length-scale criterion and error bound to make the claim self-contained. The formal proof (detailed in Section 3 of the manuscript) derives that the shape demagnetizing field contribution from spatial variations of M decays as O((a/L)^2) where L is the sample linear size and a the discretization cell size; this becomes negligible for L ≫ a with the leading error term independent of the exchange length for the magnetostatic interaction. We will revise the abstract to include an explicit statement of this criterion (e.g., sample dimensions exceeding ~10 cell sizes) together with a reference to the derived error scaling in the main text. revision: yes

  2. Referee: [Method / Results] The proposed modification discards all higher-order spatial variations of the magnetization when computing the shape correction; without numerical benchmarks comparing the corrected PBC result against a full finite-shape calculation for samples near the 'sufficiently large' threshold, it is impossible to assess the practical range of validity.

    Authors: The analytical proof establishes the asymptotic validity, and the results section already illustrates the corrected PBC method on large samples. However, we acknowledge that direct numerical comparisons near the threshold are needed to quantify practical errors. We will add a new subsection and accompanying figure that benchmarks the shape-corrected PBC results against reference calculations on finite (non-periodic) samples with dimensions 5–20 times the characteristic magnetostatic length, reporting the relative error in the demagnetizing field as a function of sample size. revision: yes

Circularity Check

0 steps flagged

No significant circularity; formal proof and modification are self-contained

full rationale

The paper's central derivation is a formal proof that magnetostatic shape effects reduce to the average magnetization for sufficiently large samples, followed by a direct modification to existing PBC implementations. No load-bearing step reduces by construction to a fitted parameter, self-defined quantity, or chain of self-citations whose validity depends on the present work. The proof is presented as an independent asymptotic argument on the magnetostatic kernel, and the modification applies the resulting insight without re-deriving or fitting the input quantities. The lack of an explicit size threshold is an applicability limitation rather than a circularity in the derivation chain.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The proof and correction rest on the domain assumption that the sample is large enough for higher multipole contributions to vanish; no free parameters or new entities are introduced.

axioms (1)
  • domain assumption For sufficiently large samples only the average magnetization produces non-negligible shape effects
    Invoked in the abstract as the basis for both the proof and the practical modification; no size threshold or error bound is supplied.

pith-pipeline@v0.9.0 · 5407 in / 1247 out tokens · 31699 ms · 2026-05-10T16:42:49.927280+00:00 · methodology

discussion (0)

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

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Exact demagnetisation field for periodic one-dimensional array of rectangular prisms

    cond-mat.mtrl-sci 2026-04 accept novelty 7.0

    An exact closed-form expression is derived for the on-axis demagnetization field of a periodically repeating infinite array of rectangular prisms, becoming exact in the thin-prism limit.

Reference graph

Works this paper leans on

38 extracted references · 1 canonical work pages · cited by 1 Pith paper

  1. [1]

    BecauseH i is linear in magnetisation, we can separate the average magnetisation from the rest by writing Hi =H avg i + ∆Hi (A3) whereH avg i and ∆H i are also given by eq

    Decomposing magnetic field The field from domain Γ i at positionris given by Hi(r) = 1 4π Z Γ 3[M(r′)·∆ˆ r]∆ˆ r−M(r′) ∆r3 dr′,(A2) wherer ′ is position relative to the center of Γ i denoted ri, while ∆r=r−r ′ −r i and ∆ˆ r=∆r ∆r with ∆r=|∆r|. BecauseH i is linear in magnetisation, we can separate the average magnetisation from the rest by writing Hi =H av...

  2. [2]

    Rewriting the field from non-uniform magnetisation The magnetic field resulting from non-uniform mag- netisation is given by ∆H i from eq. (A3), which we may 7 write as ∆Hi =− 1 4π Z Γ Ndip(∆r)∆M(r′)dr′,(A4) where Ndip αβ (∆r) = 1 ∆r3 δαβ − 3∆rβ∆rα ∆r5 (A5) is the demagnetisation tensor for a dipole,α, β=x, y, z while ∆r x = ∆xand similarly for ∆r y,∆r z....

  3. [3]

    To this end, let Λ be the union of Γ and itsNclosest copies and let¬Λ =R 3 −Λ be the rest of 3D space

    Proof of convergence for higher order contributions In the remainder of this proof, we show that the com- bined ∆Hfield from all domain copies beyond a certain distance is negligible. To this end, let Λ be the union of Γ and itsNclosest copies and let¬Λ =R 3 −Λ be the rest of 3D space. We are interested in an upper bound for the field magnitude ∆H¬Λ = 1 4...

  4. [4]

    Exchanging sum and limit Technically to take theR Λ − → ∞limit term-by-term in eq. (A11), one first has to prove the infinite sum has a finite value, so we note that whenR Λ ≥2R Γ ∞X n=1 ∥Gn∆M∥max n+ 1 n2 + 3n+ 2 R2 ΛRn Γ (RΛ −R Γ)n+2 = ∞X n=1 An∥Gn∆M∥max < ∞X n=1 ∥Gn∆M∥max, sinceA n ≤A 1 ≤1 withA 1 = 1 forR Λ = 2RΓ. Because RΓ is fixed whileR Λ can be ma...

  5. [5]

    Novel magnets through interdisiplinarity and nanocomposites

    Necessity of macroscopically uniform magnetisation With PBCs the sample magnetisation is modelled as uniform on scales much larger than the characteristic length of the simulated domain,R Γ. Here we show the necessity of this macroscopic uniformity for the preceding proof to hold. Consider the magnetisation M=MΘ(r≥R 1)Θ(r≤R 2)ˆ z ( 1θ < π 4 ∨θ > 3π 4 −1θ∈...

  6. [6]

    Oommf user’s guide: Version 1.0, 1999

    Michael Joseph Donahue and Donald G Porter. Oommf user’s guide: Version 1.0, 1999

    1999

  7. [7]

    The design and verification of mu- max3.AIP advances, 4(10), 2014

    Arne Vansteenkiste, Jonathan Leliaert, Mykola Dvornik, Mathias Helsen, Felipe Garcia-Sanchez, and Bartel Van Waeyenberge. The design and verification of mu- max3.AIP advances, 4(10), 2014

    2014

  8. [8]

    Bjørk, E.B

    R. Bjørk, E.B. Poulsen, K.K. Nielsen, and A.R. Insinga. MagTense: A micromagnetic framework using the an- alytical demagnetization tensor.Journal of Magnetism and Magnetic Materials, 535:168057, October 2021

    2021

  9. [9]

    OOMMF user’s guide, version 1.0

    M J Donahue and D G Porter. OOMMF user’s guide, version 1.0. Technical Report NIST IR 6376, National In- stitute of Standards and Technology, Gaithersburg, MD, 1999

    1999

  10. [10]

    Pe- riodic boundary conditions for demagnetization interac- tions in micromagnetic simulations.Journal of Physics D: Applied Physics, 41(17):175005, September 2008

    K M Lebecki, M J Donahue, and M W Gutowski. Pe- riodic boundary conditions for demagnetization interac- tions in micromagnetic simulations.Journal of Physics D: Applied Physics, 41(17):175005, September 2008

    2008

  11. [11]

    Two-dimensional pe- riodic boundary conditions for demagnetization interac- tions in micromagnetics.Computational Materials Sci- ence, 49(1):84–87, June 2010

    Weiwei Wang, Congpu Mu, Bin Zhang, Qingfang Liu, Jianbo Wang, and Desheng Xue. Two-dimensional pe- riodic boundary conditions for demagnetization interac- tions in micromagnetics.Computational Materials Sci- ence, 49(1):84–87, June 2010

    2010

  12. [12]

    Wysocki and Vladimir P

    Aleksander L. Wysocki and Vladimir P. Antropov. Mi- cromagnetic simulations with periodic boundary condi- tions: Hard-soft nanocomposites.Journal of Magnetism and Magnetic Materials, 428:274–286, April 2017

    2017

  13. [13]

    Strayfield calculation for micro- magnetic simulations using true periodic boundary con- ditions.Scientific Reports, 11(1):9202, April 2021

    Florian Bruckner, Amil Ducevic, Paul Heistracher, Claas Abert, and Dieter Suess. Strayfield calculation for micro- magnetic simulations using true periodic boundary con- ditions.Scientific Reports, 11(1):9202, April 2021

    2021

  14. [14]

    Peri- odic micromagnetic finite element method.Journal of Magnetism and Magnetic Materials, 615:172776, March 2025

    Fangzhou Ai, Jiawei Duan, and Vitaliy Lomakin. Peri- odic micromagnetic finite element method.Journal of Magnetism and Magnetic Materials, 615:172776, March 2025

    2025

  15. [15]

    Micromagnetic simulations with periodic strayfield calculation of soft magnetic composite-materials.Scientific Reports, 15(1):17782, May 2025

    Amil Ducevic, Florian Bruckner, Claas Abert, Christoph Vogler, and Dieter Suess. Micromagnetic simulations with periodic strayfield calculation of soft magnetic composite-materials.Scientific Reports, 15(1):17782, May 2025

    2025

  16. [16]

    Almind, Matti Knaapila, Ziwei Ouyang, Marco Beleggia, and Cathrine Frandsen

    Mathias Zambach, Miriam Var´ on, Mads R. Almind, Matti Knaapila, Ziwei Ouyang, Marco Beleggia, and Cathrine Frandsen. The effect of demagnetization on the susceptibility of single-domain particles and assemblies. (arXiv:2502.03394), February 2026

    work page arXiv 2026

  17. [17]

    Hans Fangohr, Giuliano Bordignon, Matteo Franchin, Andreas Knittel, Peter A. J. De Groot, and Thomas Fis- chbacher. A new approach to (quasi) periodic bound- ary conditions in micromagnetics: The macrogeometry. Journal of Applied Physics, 105(7):07D529, April 2009

    2009

  18. [18]

    The design and verification of Mu- Max3.AIP Advances, 4(10), October 2014

    Arne Vansteenkiste, Jonathan Leliaert, Mykola Dvornik, Mathias Helsen, Felipe Garcia-Sanchez, and Bartel Van Waeyenberge. The design and verification of Mu- Max3.AIP Advances, 4(10), October 2014

    2014

  19. [19]

    Demagnetiz- ing Field in Micromagnetic Simulation Under Periodic Boundary Conditions.IEEE Transactions on Magnetics, 47(5):902–905, May 2011

    Tetsuji Matsuo and Yuya Yamazaki. Demagnetiz- ing Field in Micromagnetic Simulation Under Periodic Boundary Conditions.IEEE Transactions on Magnetics, 47(5):902–905, May 2011

    2011

  20. [20]

    Mag- netic susceptibility of soft magnetic composite mate- rials.Journal of Magnetism and Magnetic Materials, 551:169103, June 2022

    Graham Weir, J´ erˆ ome Leveneur, and Nick Long. Mag- netic susceptibility of soft magnetic composite mate- rials.Journal of Magnetism and Magnetic Materials, 551:169103, June 2022

    2022

  21. [21]

    Arzbacher, P

    S. Arzbacher, P. Amann, B. Weidenfeller, T. Loert- ing, A. Ostermann, and J. Petrasch. Tomography based numerical simulation of the demagnetizing field in soft magnetic composites.Journal of Applied Physics, 117(16):163905, April 2015

    2015

  22. [22]

    Baars, Albert P

    Alessio Caciagli, Roel J. Baars, Albert P. Philipse, and Bonny W.M. Kuipers. Exact expression for the magnetic field of a finite cylinder with arbitrary uniform magneti- zation.Journal of Magnetism and Magnetic Materials, 456:423–432, June 2018

    2018

  23. [23]

    Durhuus, Ellen Fogh, Thomas Jauho, and Marco Beleggia

    Frederik L. Durhuus, Ellen Fogh, Thomas Jauho, and Marco Beleggia. Demagnetization factors and field of hemispherical objects.Journal of Magnetism and Mag- netic Materials, 626:173047, August 2025

    2025

  24. [24]

    Smith, K

    A. Smith, K. K. Nielsen, D. V. Christensen, C. R. H. Bahl, R. Bjørk, and J. Hattel. The demagnetizing field of a nonuniform rectangular prism.Journal of Applied Physics, 107(10):103910, May 2010

    2010

  25. [25]

    The Stray and Demagnetizing Field of a Homogeneously Magnetized Tetrahedron.IEEE Mag- netics Letters, 10:1–5, 2019

    Kaspar Kirstein Nielsen, Andrea Roberto Insinga, and Rasmus Bjork. The Stray and Demagnetizing Field of a Homogeneously Magnetized Tetrahedron.IEEE Mag- netics Letters, 10:1–5, 2019

    2019

  26. [26]

    J. A. Osborn. Demagnetizing Factors of the General El- lipsoid.Physical Review, 67(11-12):351–357, June 1945

    1945

  27. [27]

    Beleggia, M

    M. Beleggia, M. De Graef, and Y. Millev. Demagnetiza- tion factors of the general ellipsoid: An alternative to the Maxwell approach.Philosophical Magazine, 86(16):2451– 2466, June 2006

    2006

  28. [28]

    Beleggia and M

    M. Beleggia and M. De Graef. On the computation of the demagnetization tensor field for an arbitrary par- ticle shape using a Fourier space approach.Journal of Magnetism and Magnetic Materials, 263(1-2):L1–L9, July 2003

    2003

  29. [29]

    R. I. Joseph and E. Schl¨ omann. Demagnetizing Field in Nonellipsoidal Bodies.Journal of Applied Physics, 36(5):1579–1593, May 1965

    1965

  30. [30]

    Moskowitz and E

    R. Moskowitz and E. Della Torre. Theoretical aspects of demagnetization tensors.IEEE Transactions on Mag- netics, 2(4):739–744, December 1966

    1966

  31. [31]

    Clarendon Press ; Oxford University Press, third edition, 1891

    James Clerk Maxwell.A Treatise on Electricity and Mag- netism, volume 2. Clarendon Press ; Oxford University Press, third edition, 1891. Chapter V, part 437

  32. [32]

    K. L. Pisane, Sobhit Singh, and M. S. Seehra. Unusual enhancement of effective magnetic anisotropy with de- creasing particle size in maghemite nanoparticles.Ap- plied Physics Letters, 110(22):222409, May 2017

    2017

  33. [33]

    Ox- ford University Press, first edition edition, 2001

    Stephen Blundell.Magnetism in Condensed Matter. Ox- ford University Press, first edition edition, 2001. ap- pendix D

    2001

  34. [34]

    Reinert, A

    J. Reinert, A. Brockmeyer, and R.W.A.A. De Doncker. Calculation of losses in ferro- and ferrimagnetic mate- rials based on the modified Steinmetz equation.IEEE 10 Transactions on Industry Applications, 37(4):1055–1061, July-Aug./2001

    2001

  35. [35]

    Dover, New York, 1990

    Konrad Knopp.Theory and application of infinite series. Dover, New York, 1990. Chapter III,§13, pages 116-120

    1990

  36. [36]

    K. Halbach. Design of permanent multipole magnets with oriented rare earth cobalt material.Nuclear Instruments and Methods, 169(1):1–10, February 1980

    1980

  37. [37]

    Bjørk, A

    R. Bjørk, A. Smith, and C.R.H. Bahl. Analysis of the magnetic field, force, and torque for two-dimensional Hal- bach cylinders.Journal of Magnetism and Magnetic Ma- terials, 322(1):133–141, January 2010

    2010

  38. [38]

    Bjørk, A

    R. Bjørk, A. Smith, and C.R.H. Bahl. The efficiency and the demagnetization field of a general Halbach cylinder. Journal of Magnetism and Magnetic Materials, 384:128– 132, June 2015

    2015