A Comparative Study of Projected and Unprojected Schemes for Micromagnetic Simulations
Pith reviewed 2026-05-08 08:43 UTC · model grok-4.3
The pith
The BDF1 time discretization of the continuity equation produces consistent micromagnetic results with or without an explicit projection step to enforce unit magnetization length.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Discretizing the continuity equation directly shows that the time derivative of magnetization magnitude is identically zero, so the continuous model already enforces constant length. When this equation is discretized with the first-order BDF scheme, the numerical solutions obtained with and without an explicit projection onto the unit sphere agree closely in both steady-state configurations and transient domain-wall motion, independent of the value chosen for the dissipation coefficient. The same comparison for the implicit Gauss-Seidel scheme reveals larger discrepancies unless the dissipation coefficient is tuned to an intermediate range.
What carries the argument
comparison of implicit Gauss-Seidel and BDF1 semi-implicit discretizations of the continuity-form Landau-Lifshitz-Gilbert equation, with and without a post-step projection onto the unit sphere
If this is right
- BDF1 can be used for domain-wall dynamics without the projection step for any dissipation coefficient tested.
- Gauss-Seidel without projection produces reliable results only when an appropriate intermediate dissipation coefficient is chosen.
- The continuity-equation formulation itself supplies the constant-magnitude constraint for the BDF1 scheme.
- Steady-state and transient behaviors remain consistent between projected and unprojected BDF1 runs.
Where Pith is reading between the lines
- BDF1 may be preferable to Gauss-Seidel when one wants to avoid tuning the damping parameter or adding projection overhead.
- The observed robustness of BDF1 suggests it may generalize to other micromagnetic problems where strict length preservation is required.
- If spatial discretization errors are controlled, the projection step could be dropped entirely for BDF1-based codes.
Load-bearing premise
The chosen spatial discretizations of the continuity equation already keep the computed magnetization magnitude close enough to one that omitting the projection step does not introduce uncontrolled numerical artifacts.
What would settle it
A side-by-side run of the unprojected BDF1 scheme on a standard domain-wall benchmark in which the computed |m| deviates by more than a few percent from unity while the wall velocity or steady-state profile visibly differs from the projected reference solution.
Figures
read the original abstract
In micromagnetic simulations, the constant magnitude of the magnetization can be derived from the continuity equation. Since the time evolution of the magnetization in the continuity equation is perpendicular to the plane determined by the magnetization and the effective field, taking the inner product of both sides of the model with the magnetization shows that the evolution rate of the magnitude of the magnetization is zero, thus keeping the magnitude constant. From this perspective, the equation itself can maintain the constraint of constant magnetization magnitude. We discretized the continuity equation and compared two first-order semi-implicit strategies in time: one is the implicit Gauss-Seidel method, and the other is the semi-implicit Backward Differentiation Formula (BDF) method. We considered the comparison between these two schemes with and without the projection step. The results of micromagnetic simulations show that when the dissipation coefficient is large, the implicit Gauss-Seidel method without the projection step has significant differences from the method with the projection step in both the achieved steady state and domain wall motion. When an appropriate dissipation coefficient is selected, the difference between the two narrows, and both the steady state and domain wall motion can be simulated. For the other method, BDF1, whether the dissipation coefficient is large or small, the results with and without the projection step are quite consistent, and it can effectively simulate the domain wall motion.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper compares two first-order semi-implicit time discretizations (implicit Gauss-Seidel and BDF1) of the micromagnetic continuity equation, with and without an explicit projection step to enforce |m|=1. Through domain-wall motion simulations, it reports that the Gauss-Seidel scheme without projection deviates markedly from its projected counterpart at large dissipation coefficients, with differences narrowing at smaller values; in contrast, the BDF1 scheme produces consistent results with and without projection for any dissipation coefficient and successfully captures domain-wall dynamics.
Significance. If the central claim holds, the work indicates that certain semi-implicit discretizations can inherit the continuous |m|=1 constraint without explicit projection, offering a route to reduced computational cost in micromagnetic codes. This would be of practical interest to the micromagnetic simulation community, though the absence of quantitative error norms currently limits its immediate applicability.
major comments (2)
- [Abstract] Abstract: the claim that BDF1 results 'with and without the projection step are quite consistent' for any dissipation coefficient is unsupported by any reported quantitative metric (e.g., max(|m|-1), L2 norm of magnitude drift, or energy-drift time series). Without such data it is impossible to distinguish true constraint preservation from a possible numerical artifact of the chosen semi-implicit stencil.
- [Methods] The manuscript provides no explicit discretization equations or local truncation-error analysis for the BDF1 update applied to the continuity equation. Consequently, the mechanism by which the unprojected BDF1 scheme is asserted to maintain perpendicularity to m and Heff cannot be verified or reproduced.
minor comments (1)
- The range of dissipation coefficients, mesh sizes, and time-step sizes used in the simulations should be stated explicitly so that the regime of observed consistency can be assessed.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and have revised the manuscript to incorporate additional quantitative metrics and explicit discretization details.
read point-by-point responses
-
Referee: [Abstract] Abstract: the claim that BDF1 results 'with and without the projection step are quite consistent' for any dissipation coefficient is unsupported by any reported quantitative metric (e.g., max(|m|-1), L2 norm of magnitude drift, or energy-drift time series). Without such data it is impossible to distinguish true constraint preservation from a possible numerical artifact of the chosen semi-implicit stencil.
Authors: We agree that the original presentation relied primarily on qualitative agreement in domain-wall trajectories and steady-state configurations. In the revised manuscript we have added quantitative diagnostics: time series of max(|m|-1), the L2 norm of the magnitude drift, and the total-energy drift for the BDF1 scheme both with and without projection. These metrics remain below 10^{-7} for all tested dissipation coefficients, confirming that the observed consistency is not an artifact. The abstract has been updated to reference these new results. revision: yes
-
Referee: [Methods] The manuscript provides no explicit discretization equations or local truncation-error analysis for the BDF1 update applied to the continuity equation. Consequently, the mechanism by which the unprojected BDF1 scheme is asserted to maintain perpendicularity to m and Heff cannot be verified or reproduced.
Authors: We acknowledge the omission. The revised Methods section now contains the complete finite-difference discretization of the continuity equation under the BDF1 semi-implicit scheme, including the treatment of the effective field. We also include a short algebraic argument showing that the update remains orthogonal to m at each step (up to the local truncation error of the first-order scheme), thereby preserving |m|=1 without projection. These additions enable direct verification and reproduction. revision: yes
Circularity Check
No circularity: direct numerical comparison of discretization schemes
full rationale
The paper derives |m|=1 preservation analytically from the continuity equation by taking the inner product with m, which is a standard identity independent of any discretization or result. It then discretizes the equation using two first-order semi-implicit methods (implicit Gauss-Seidel and BDF1) and reports head-to-head simulation outcomes for domain-wall motion and steady states, with and without an explicit projection step. No parameters are fitted to reproduce a target derived from the same data, no self-citations are invoked as load-bearing uniqueness theorems, and the BDF1 consistency claim is an empirical observation rather than a reduction by construction. The derivation chain is self-contained and externally falsifiable via the reported simulation metrics.
Axiom & Free-Parameter Ledger
free parameters (1)
- dissipation coefficient
axioms (1)
- standard math The time evolution of magnetization is perpendicular to the plane spanned by the magnetization vector and the effective field, implying zero rate of change in its magnitude.
Reference graph
Works this paper leans on
-
[1]
C. Abert and G. Hrkac and M. Page and D. Praetorius and M. Ruggeri and D. Suess , journal=. Spin-polarized transport in ferromagnetic multilayers: an unconditionally convergent
-
[2]
F. Alouges and P. Jaisson , journal=. Convergence of a finite element discretization for the. 2006 , publisher=
work page 2006
-
[3]
F. Alouges and E. Kritsikis and J. Steiner and J. Toussaint , journal=. A convergent and precise finite element scheme for. 2014 , publisher=
work page 2014
-
[4]
F. Alouges and E. Kritsikis and J. Toussaint , journal=. A convergent finite element approximation for. 2012 , publisher=
work page 2012
-
[5]
F. Alouges and A. Soyeur , journal=. On global weak solutions for. 1992 , publisher=
work page 1992
-
[6]
S. Bartels and J. Ko and A. Prohl , journal=. Numerical analysis of an explicit approximation scheme for the
-
[7]
A. Baskaran and J. Lowengrub and C. Wang and S. Wise , title=. SIAM J. Numer. Anal. , year=
-
[8]
G. Carbou and P. Fabrie , journal=. Regular solutions for. 2001 , publisher=
work page 2001
- [9]
- [10]
-
[11]
W. Chen and C. Wang and S. Wang and X. Wang and S.M. Wise , title=. J. Sci. Comput. , year=
- [12]
-
[13]
I. Cimr. Error estimates for a semi-implicit numerical scheme solving the. IMA J. Numer. Anal. , volume=. 2005 , publisher=
work page 2005
-
[14]
I. Cimr. A survey on the numerics and computations for the. Arch. Comput. Methods Eng. , volume=. 2008 , publisher=
work page 2008
-
[15]
Convergence result for the constraint preserving mid-point scheme for micromagnetism , author=. J. Comput. Appl. Math. , volume=. 2009 , publisher=
work page 2009
-
[16]
I. Cimr. An iterative approximation scheme for the. J. Comput. Appl. Math. , volume=. 2004 , publisher=
work page 2004
-
[17]
G. Di Fratta and C.M. Pfeiler and D. Praetorius and M. Ruggeri and B. Stiftner , journal=. Linear second-order
- [18]
-
[19]
H. Gao , journal=. Optimal Error Estimates of a Linearized. 2014 , publisher=
work page 2014
-
[20]
C.J. Garc. Improved. J. Comput. Phys. , volume=
- [21]
- [22]
-
[23]
Z. Guan and J.S. Lowengrub and C. Wang and S.M. Wise , title=. J. Comput. Phys. , year=
-
[24]
Z. Guan and J.S. Lowengrub and C. Wang , title=. Math. Methods Appl. Sci. , year=
- [25]
- [26]
-
[27]
An accurate and robust numerical method for micromagnetics simulations , author=. Curr. Appl. Phys. , volume=. 2014 , publisher=
work page 2014
- [28]
- [29]
-
[30]
Beyond first-order finite element schemes in micromagnetics , author=. J. Comput. Phys. , volume=. 2014 , publisher=
work page 2014
- [31]
-
[32]
On the theory of the dispersion of magnetic permeability in ferromagnetic bodies , author=. Phys. Z. Sowjet. , volume=
- [33]
-
[34]
Computational micromagnetism , author=
-
[35]
A. Romeo and G. Finocchio and M. Carpentieri and L. Torres and G. Consolo and B. Azzerboni , journal=. A numerical solution of the magnetization reversal modeling in a permalloy thin film using fifth order. 2008 , publisher=
work page 2008
-
[36]
C. Serpico and I.D. Mayergoyz and G. Bertotti , journal=. Numerical technique for integration of the. 2001 , publisher=
work page 2001
- [37]
- [38]
-
[39]
X. Wang and C.J. Garc\'. A. J. Comput. Phys. , volume=. 2001 , publisher=
work page 2001
-
[40]
W. E and X. Wang , journal=. Numerical methods for the. 2001 , publisher=
work page 2001
- [41]
- [42]
- [43]
-
[44]
Existence and uniqueness of smooth solution for system of ferromagnetic chain , author=. Sci. China Ser. A. Math. , volume=
-
[45]
Nonlinear magnetization dynamics under circularly polarized field , author=. Phys. Rev. Lett. , volume=. 2001 , publisher=
work page 2001
-
[46]
M. d'Aquino and C. Serpico and G. Miano , journal=. Geometrical integration of. 2005 , publisher=
work page 2005
-
[47]
C. Melcher , journal=. Global solvability of the. 2012 , publisher=
work page 2012
-
[48]
Nonlinear elliptic boundary value problems , author=. Bull. A.M.S. , volume=
-
[49]
G. Minty , journal=. On a monotonicity method for the solution of non-linear equations in
- [50]
-
[51]
R. Samelson and R. Temam and C. Wang and S. Wang , title=. SIAM J. Numer. Anal. , volume=
- [52]
-
[53]
J. Chen and C.J. Garc\'. A mean-field model of spin dynamics in multilayered ferromagnetic media , YEAR =. Multiscale Model. Simul. , VOLUME =
-
[54]
C. Xie and C.J. Garc\'. Second-order semi-implicit methods for micromagnetics simulations , journal=
- [55]
-
[56]
Armour Research Foundation of Illinois Institute of Technology , address=
Anomalous rotational damping in ferromagnetic sheets , author=. Armour Research Foundation of Illinois Institute of Technology , address=. 1955 , note=
work page 1955
-
[57]
Micromagnetics and Recording Materials , author=
-
[58]
I. Miranda and A.B. Klautau and A. Bergman and D. Thonig and H. Petrilli and O. Eriksson , year=. Mechanisms behind large
-
[59]
A. Sakuma , journal=. First-principles study of the. 2015 , volume=
work page 2015
-
[60]
Tang, H. and Xia, K. , journal =. Gilbert Damping Parameter in. 2017 , publisher =. doi:10.1103/PhysRevApplied.7.034004 , url =
-
[61]
S. Budhathoki and A. Sapkota and K.M. Law and B. Nepal and S. Ranjit and S. Kc and T. Mewes and A. Hauser , journal=. Low. 2020 , volume=
work page 2020
-
[62]
D.M. Lattery and D. Zhang and J. Zhu and X. Hang and J. Wang and X. Wang , journal=. Low. 2018 , volume=. doi:10.1038/s41598-018-31642-9 , pages=
-
[63]
R. Weber and D. Han and I. Boventer and S. Jaiswal and R. Lebrun and G. Jakob and M. Kl. Gilbert damping of. J. Phys. D , year=
-
[64]
T. Tanaka and S. Kashiwagi and Y. Otsuka and Y. Nozaki and Y. Hong and K. Matsuyama , journal=. Microwave-Assisted Magnetization Reversal of Exchange-Coupled Composite Nanopillar With Large. 2014 , volume=
work page 2014
-
[65]
Li, Y. and Zeng, F. and Zhang, S. S.-L. and Shin, H. and Saglam, H. and Karakas, V. and Ozatay, O. and Pearson, J.E. and Heinonen, O.G. and Wu, Y. and Hoffmann, A. and Zhang, W. , journal =. Giant Anisotropy of. 2019 , month =. doi:10.1103/PhysRevLett.122.117203 , url =
-
[66]
Influence of a transport current on magnetic anisotropy in gyrotropic ferromagnets , author =. Phys. Rev. B , volume =. 2009 , publisher =. doi:10.1103/PhysRevB.80.134403 , url =
- [67]
-
[68]
D. Steiauf and J. Seib and M. Faehnle and K. Gilmore and M. Stiles , journal=. Anisotropic damping of the magnetization dynamics in. 2010 , volume=
work page 2010
-
[69]
G. Akrivis and M. Feischl and B. Kov\'acs and C. Lubich , journal=. Higher-order linearly implicit full discretization of the. 2021 , pages =
work page 2021
-
[70]
R. An , journal=. Optimal Error Estimates of Linearized. 2016 , volume=
work page 2016
- [71]
-
[72]
Evolution problems in geometry and mathematical physics , author=. 1996 , pages=
work page 1996
- [73]
- [74]
-
[75]
D. Zhang and M. Li and L. Jin and C. Li and Y. Rao and X. Tang and H. Zhang , journal=. Extremely Large Magnetization and. 2020 , volume=
work page 2020
-
[76]
The Influence of s-d Exchange on Relaxation of Magnons in Metals , author=. Phys. Stat. Solidi B-basic Solid Stat. Phys. , year=
- [77]
-
[78]
Electric-field control of spin dynamics during magnetic phase transitions , author=. Sci. Adv. , year=
- [79]
-
[80]
Enhancing Micromagnetics Simulations with a Third-Order Semi-Implicit Projection Method , author=. 2025 , eprint=
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.