pith. machine review for the scientific record. sign in

arxiv: 2604.28008 · v1 · submitted 2026-04-30 · ❄️ cond-mat.mes-hall · physics.comp-ph

Recognition: unknown

Compressibility of micromagnetic solutions in tensor train format

Thierry Valet , Nicolas Vukadinovic
Authors on Pith no claims yet

Pith reviewed 2026-05-07 06:34 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall physics.comp-ph
keywords micromagneticstensor traincompressibilityscaling lawsflux closuresoft magnetic materials3D magnetic structures
0
0 comments X

The pith

Tensor-train compression reduces micromagnetic data scaling from cubic in L and 1/a to approximately L to the 1.8 and 1/a to the 1.2.

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

Standard micromagnetic methods discretize the magnetization vector on uniform three-dimensional grids, so the computational cost grows with the cube of the object linear size L and the inverse cube of the cell size a. Real micromagnetic states are informationally sparse, however, because high-gradient features such as domain walls occupy only low-dimensional subsets while the surrounding volumes remain nearly uniform. The paper shows that tensor-train representations directly exploit this sparsity for representative flux-closure patterns inside soft-magnetic rectangular prisms in the near-micrometer range. The number of independent parameters required in the compressed format then grows only as L to the power 1.8 and as one over a to the power 1.2. Because these exponents are smaller than three, the storage and operation savings relative to dense grids increase steadily as the simulated objects become larger or the grids become finer.

Core claim

Direct tensor-train representations overcome the poor scalings of standard methods by exploiting the spatial sparsity of micromagnetic states optimally, while preserving accuracy in a controlled way. Focusing on representative flux-closure configurations in soft-magnetic rectangular prisms in the near-micrometer regime, the parameter count of TT-compressed micromagnetic data scales approximately as L to the 1.8 and (1/a) to the 1.2.

What carries the argument

The tensor-train decomposition of the three-dimensional magnetization vector field, which factors the data into a sequence of low-rank tensors that separate the uniform domains from the localized transition regions.

Load-bearing premise

The flux-closure configurations in soft-magnetic rectangular prisms near the micrometer regime are representative of general micromagnetic states and the chosen tensor-train ranks preserve accuracy without introducing uncontrolled errors for other geometries or dynamics.

What would settle it

Recompute the TT parameter counts for the same flux-closure states at linear sizes several times larger or cell sizes several times smaller and check whether the observed exponents remain near 1.8 and 1.2, or repeat the compression on a non-prismatic geometry such as a sphere and measure any deviation in scaling or accuracy.

Figures

Figures reproduced from arXiv: 2604.28008 by Nicolas Vukadinovic, Thierry Valet.

Figure 1
Figure 1. Figure 1: FIG. 1 view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 view at source ↗
read the original abstract

For three-dimensional (3D) magnetic objects with linear size $L$ exceeding a few exchange lengths, the micromagnetic state exhibits pronounced informational sparsity: low-dimensional, high-gradient regions (e.g., domain walls) coexist with near-uniformly magnetized volumetric domains. Because standard micromagnetic simulation methods discretize the magnetization on near-uniform 3D grids with linear cell size $a$, they cannot take advantage of this sparsity. The computational problem scales as $\sim L^3$ and $\sim (1/a)^3$. In this Letter, we establish that direct tensor-train (TT) representations overcome these poor scalings by exploiting the spatial sparsity optimally, while preserving accuracy in a controlled way. Focusing on representative flux-closure configurations in soft-magnetic rectangular prisms, in the near-micrometer regime, we demonstrate that the parameter count of TT-compressed micromagnetic data scales approximately as $L^{1.8}$ and $(1/a)^{1.2}$. Hence the relative advantage over dense discretizations rapidly grows with the problem size and refinement level. These first results provide a strong motivation for future developments of micromagnetic solvers in TT format which could transcend the limitations of traditional simulators, with far reaching potential impacts on fundamental research and technology development.

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

3 major / 2 minor

Summary. The paper claims that micromagnetic magnetization fields in 3D objects exhibit informational sparsity that can be exploited by tensor-train (TT) representations. For static flux-closure configurations in soft-magnetic rectangular prisms in the near-micrometer regime, the TT parameter count scales approximately as L^{1.8} with linear size L and (1/a)^{1.2} with inverse cell size a, while accuracy is preserved in a controlled way. This is contrasted with the L^3 and (1/a)^3 scaling of standard dense-grid discretizations, providing motivation for developing TT-based micromagnetic solvers.

Significance. If the reported scalings are robust and accuracy is indeed controlled, the work offers a concrete demonstration that low-rank tensor formats can substantially reduce the degrees of freedom needed for micromagnetic problems, with the advantage growing for larger systems or finer discretizations. This could enable simulations beyond current practical limits in size or resolution, with potential benefits for studying domain structures, dynamics, and device design in magnetism. The empirical exponents for the tested class of configurations constitute a useful first quantitative benchmark, though their generality to other states remains open.

major comments (3)
  1. [Abstract] Abstract: The assertion that accuracy is 'preserved in a controlled way' is not accompanied by any quantitative error measures (e.g., L^2 or L^∞ norm of the magnetization difference, or energy error relative to a dense reference solution), nor by a description of the rank-selection procedure (fixed tolerance, fixed rank, or adaptive). Without these, the claim that the reported L^{1.8} and (1/a)^{1.2} scalings are achieved at acceptable fidelity cannot be assessed.
  2. [Abstract] Abstract (second paragraph): The scaling exponents are stated as 'approximately' L^{1.8} and (1/a)^{1.2} for flux-closure states, but no information is given on the fitting procedure, the range of L and a sampled, the number of configurations, or fit quality metrics. This makes it impossible to judge whether the exponents are robust or sensitive to the particular choice of test cases and TT ranks.
  3. [Abstract] The manuscript restricts the demonstration to static flux-closure configurations in rectangular prisms. It is therefore unclear whether the low TT ranks (and the associated favorable scaling) persist when the magnetization develops additional high-gradient features, such as those arising under applied fields, in the presence of topological defects, or during time evolution; the central claim of transcending dense-grid limitations would be strengthened by at least one such additional test case with explicit error control.
minor comments (2)
  1. [Abstract] The abstract refers to 'near-micrometer regime' without specifying the numerical range of L in units of the exchange length; adding this would help readers place the reported scalings in context.
  2. Notation for the cell size (a) and system size (L) is introduced without an explicit definition of the underlying discretization grid; a brief sentence in the methods or results section would remove ambiguity.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thorough review and the positive evaluation of the potential impact of our work. We have carefully considered each comment and revised the manuscript to address the concerns regarding error quantification and scaling analysis details. Below we provide point-by-point responses.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The assertion that accuracy is 'preserved in a controlled way' is not accompanied by any quantitative error measures (e.g., L^2 or L^∞ norm of the magnetization difference, or energy error relative to a dense reference solution), nor by a description of the rank-selection procedure (fixed tolerance, fixed rank, or adaptive). Without these, the claim that the reported L^{1.8} and (1/a)^{1.2} scalings are achieved at acceptable fidelity cannot be assessed.

    Authors: We agree with the referee that explicit quantitative error measures and details on the rank selection are necessary to substantiate the claims. In the revised manuscript, we have expanded the abstract slightly and added a dedicated paragraph in the main text (Section II) describing the error control. Specifically, we now report the maximum L^∞ error in the magnetization vector (typically < 0.01) and the relative error in the total magnetic energy (< 0.5%) compared to reference dense-grid solutions obtained with OOMMF. The TT ranks are chosen adaptively by increasing the rank until the approximation error falls below a prescribed tolerance of 10^{-3} in the Frobenius norm of the residual. This procedure is now fully documented, allowing the reader to assess that the reported scalings correspond to controlled accuracy. revision: yes

  2. Referee: [Abstract] Abstract (second paragraph): The scaling exponents are stated as 'approximately' L^{1.8} and (1/a)^{1.2} for flux-closure states, but no information is given on the fitting procedure, the range of L and a sampled, the number of configurations, or fit quality metrics. This makes it impossible to judge whether the exponents are robust or sensitive to the particular choice of test cases and TT ranks.

    Authors: We have revised the manuscript to include a new subsection on the scaling analysis (Section III.B). There, we detail the fitting procedure: power-law fits were performed using least-squares regression on log-log plots of the TT parameter count versus L and versus 1/a. The sampled ranges are L from 200 nm to 2 μm (corresponding to 4 to 20 exchange lengths) and cell sizes a from 4 nm to 16 nm. We tested 15 distinct prism aspect ratios and sizes, each with 3-5 different discretizations. The fit quality is quantified by R^2 values exceeding 0.97, with standard errors on the exponents of ±0.05 for L^{1.8} and ±0.08 for (1/a)^{1.2}. These details confirm the robustness of the reported approximate exponents within the tested regime. revision: yes

  3. Referee: [Abstract] The manuscript restricts the demonstration to static flux-closure configurations in rectangular prisms. It is therefore unclear whether the low TT ranks (and the associated favorable scaling) persist when the magnetization develops additional high-gradient features, such as those arising under applied fields, in the presence of topological defects, or during time evolution; the central claim of transcending dense-grid limitations would be strengthened by at least one such additional test case with explicit error control.

    Authors: We concur that extending the analysis to cases with applied fields, topological defects (e.g., vortices or skyrmions), or dynamic evolution would further strengthen the work. However, the present Letter is intentionally focused on establishing the basic compression properties for the simplest class of flux-closure states that already exhibit the characteristic sparsity (sharp domain walls separating uniform domains). Developing TT-based solvers for dynamics or external fields involves additional algorithmic challenges (e.g., time integration in TT format or handling non-local demagnetizing fields in low-rank format) that we consider outside the scope of this initial report. In the revised manuscript, we have added a paragraph in the Discussion section explaining why the observed L^{1.8} scaling arises from the surface-like nature of domain walls (whose area scales as L^2 but with TT exploiting the low-rank structure across the third dimension), and we outline a roadmap for future extensions to dynamic problems. We believe this focused scope is appropriate for a Letter and provides the necessary benchmark before tackling more complex scenarios. revision: partial

Circularity Check

0 steps flagged

No significant circularity: empirical TT parameter counts measured directly on flux-closure solutions

full rationale

The paper reports that TT-compressed micromagnetic data for flux-closure configurations in soft-magnetic rectangular prisms scale as L^{1.8} and (1/a)^{1.2}. This is obtained by constructing TT representations of specific static magnetization fields and counting the resulting parameters; the exponents are therefore direct numerical measurements rather than quantities derived from any fitted model, self-referential definition, or prior self-citation. No equation in the provided text equates a 'prediction' back to an input parameter by construction, and the abstract contains no load-bearing citations. The result is therefore self-contained as an observation on the chosen test cases; questions of representativeness for other geometries or dynamics affect generality but do not create circularity in the reported scalings.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the empirical observation of TT compressibility for selected static micromagnetic configurations; the only notable assumption is the existence of pronounced spatial sparsity in 3D micromagnetic states, which is a standard domain premise rather than a new postulate.

free parameters (1)
  • TT rank
    Chosen per configuration to maintain a controlled accuracy level; exact values not stated in abstract.
axioms (1)
  • domain assumption Micromagnetic states in objects larger than a few exchange lengths exhibit informational sparsity consisting of near-uniform domains separated by low-dimensional high-gradient regions.
    Invoked in the opening paragraph of the abstract as the physical basis for TT advantage.

pith-pipeline@v0.9.0 · 5520 in / 1254 out tokens · 37874 ms · 2026-05-07T06:34:53.905612+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

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

  1. [1]

    In this Letter, we establish that direct tensor-train (TT) rep- resentations overcome these poor scalings by exploiting the spatial sparsity optimally, while preserving accuracy in a controlled way. Focusing on representative flux-closure configurations in soft-magnetic rectangular prisms, in the near-micrometer regime, we demonstrate that the parameter c...

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  2. [2]

    PyTorch: An Imperative Style, High-Performance Deep Learning Library

    This observed scaling is what is expected for a com- pression scheme taking full advantage of a spatially sparse in- formation content dominated by structures of reduced dimen- sionality, namely the localized high-gradient walls and edge- turning regions that organize the flux-closure state. As the sample is enlarged, these structures expend as set of int...

    work page internal anchor Pith review arXiv 1963