Pith. sign in

REVIEW 3 major objections 5 minor 31 references

Wave dynamics in a macroscopic square artificial spin ice

T0 review · 3 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash

Pith's one-line read A Dirac string in a macroscopic artificial spin ice traps waves in a distinct 4.9 Hz mode below the propagation band.

desk verdict Plausible but unproven claim of a 4.9 Hz Dirac-string-localized mode; the missing eigenmode analysis is load-bearing and a referee should demand it. read the letter →

arxiv 2501.09820 v1 pith:HVVGTRN7 submitted 2025-01-16 cond-mat.mes-hall

classification cond-mat.mes-hall
keywords macroscopicartificialspiniceDiracstringmechano-magneticdynamicswavelocalizationmagneticmonopolemodelscatteringbandstructurenonlinear
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper numerically simulates a macroscopic artificial spin ice—an array of hinged bar magnets in a square lattice—to ask what a line defect, a Dirac string, does to mechanical waves moving through the lattice. The authors find that the defect does not act as a simple barrier: waves scatter among modes throughout the lattice and lose coherence, but this scattering feeds a distinct resonance localized on the Dirac string itself. That mode sits near 4.9 Hz, below the 5.7–14.8 Hz band of propagating waves, as expected for a spatially confined oscillation. The result matters because it shows that mechano-magnetic macro-ices reproduce the defect-localized dynamics seen in nanoscopic spin ices, and it suggests a route to channeling mechanical waves along defects rather than through the bulk.

What carries the argument

The central object is the Dirac string: a chain of type-II vertices where both horizontal and vertical magnets share the same orientation, carrying higher energy than the surrounding type-I vortex configuration and terminating in emergent monopoles at the lattice edges. The argument is carried by the magnetic-monopole model of Eq. (1), in which each bar magnet is reduced to two magnetic charges interacting via a Coulomb-like potential, with rotary friction and an external drive; parameters (moment of inertia $I$, friction $\eta$, effective charge $q = 2.08\,\mathrm{A\,m}$) are inherited from prior macro-ASI work. The localized mode is identified in forced-dynamics simulations by Fourier-transforming each magnet's angle time trace and convolving the amplitudes with a Gaussian to visualize the mode volume, and by comparing output spectra with and without the Dirac string. The mechanism is inter-mode scattering throughout the lattice that resonantly excites the Dirac-string resonance below the propagation band.

What would settle it

Time-resolve the orientation of every magnet in the physical macro-ASI while driving at 11.2 Hz with a ~1 mT field: if no spectral peak near 4.9 Hz appears with amplitude concentrated along the Dirac string, the central claim is falsified.

Watch

Extended reading notes

Core claim

The central claim is that a Dirac string in a macroscopic square artificial spin ice supports a localized vibrational mode near 4.9 Hz, below the lowest propagating band at about 5.7 Hz. The mode is not excited by direct scattering of an incoming wave off the string; instead, the wave scatters throughout the lattice, and that distributed scattering resonantly drives the Dirac string's own oscillation. This is shown by the mode profile, which concentrates amplitude along the string, and by the nonlinear threshold of about 0.8 mT at 11.2 Hz drive. The paper argues that this is the macroscopic analogue of low-frequency modes tied to emergent monopoles in nanoscopic spin ices, with the caveat that permanent magnets here prevent edge modes, so localization is due to the differing coupling at type-II vertices.

Load-bearing premise

The load-bearing premise is that the magnetic-monopole model of Eq. (1), with effective charge $q = 2.08\,\mathrm{A\,m}$ and parameters taken from prior macro-ASI studies, faithfully represents the physical bar-magnet lattice, so that the predicted 4.9 Hz Dirac-string mode is a real property rather than a numerical artifact.

Editorial extensions

If this is right

  • A Dirac string cannot serve as a simple gate for mechano-magnetic waves, because the incoming wave scatters throughout the lattice rather than being blocked or transmitted by the string.
  • The localized mode below the band could channel excitations along the Dirac string, analogous to spin-wave propagation along domain walls, provided coherence losses from inter-mode scattering are controlled.
  • Driving the system directly at the Dirac-string resonance near 4.9 Hz may produce more coherent channeling than driving at in-band frequencies.
  • The onset of the localized mode near 0.8 mT implies a nonlinear activation threshold for defect-localized dynamics in macro-ASIs.
  • The two-band dispersion with degenerate flat regions along $k_x=0$ and $k_y=0$ and preferential diagonal propagation is a handle for steering waves in macro-ASI devices.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the monopole model's coarse-graining is the only simplification, the exact 4.9 Hz value is likely to shift in a physical realization; the robust prediction is the existence of a sub-band mode localized on the string.
  • The same scattering-driven mechanism may apply to other defect geometries, such as grain boundaries or isolated monopole pairs, making macro-ASIs a general testbed for defect-engineered mechanical wave control.
  • A direct experiment driving at 4.9 Hz and mapping the resulting amplitude should show energy concentrating on the Dirac string, which would test the channeling hypothesis independently of the in-band driving protocol used here.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper numerically studies a macroscopic square artificial spin ice composed of hinged bar magnets, using the magnetic-monopole model of Eq. (1) with parameters taken from the authors' prior work. In an extended 25x25 ground-state array it computes a dispersion relation with two bands between about 5.7 and 14.8 Hz, degenerate along kx=0 and ky=0. In a confined 60-magnet array it then stabilizes a Dirac string and drives two corner magnets with an oscillating field. The main reported results are (i) strong inter-mode scattering and loss of coherence, and (ii) a mode near 4.9 Hz, below the propagating band, whose Fourier amplitude (averaged over 4-6 Hz) is localized along the Dirac string. The authors conclude that wave scattering throughout the lattice, rather than direct scattering off the string, excites this localized mode, and they suggest that direct excitation at the resonant frequency could channel waves with better coherence.

Significance. If the localized Dirac-string mode is real, the paper provides a valuable mechanical analogue of defect-localized modes in nanoscopic artificial spin ices and demonstrates a concrete mechanism for spatial confinement in mechano-magnetic metamaterials. A particular strength is that the central result is an emergent numerical observation, not a fitted input, and the simulation recipe is fully specified by Eq. (1) and the stated parameters. However, the evidence for the existence of the localized mode is currently indirect: it is identified from a nonlinear drive and a broad spectral average, and the dispersion extraction relies on an interpolation procedure that is not independently validated. Because the paper's novelty claim depends on this mode being a genuine eigenmode of the defective lattice, the current evidence is not yet sufficient to establish the central conclusion.

major comments (3)
  1. [Section IV, Fig. 5] The existence of a distinct Dirac-string-localized mode at approximately 4.9 Hz is inferred from the response to a nonlinear 11.2 Hz drive, using Fourier amplitudes averaged over 4-6 Hz. This does not rule out the alternatives that the 4.9 Hz feature is a combination tone between in-band modes, spectral leakage from the drive, or a static deformation of type-II vertices captured by the broad spectral average. The authors should demonstrate that the linearized dynamics of Eq. (1) about the relaxed defective state has an eigenmode near 4.9 Hz whose participation is concentrated on the Dirac string, or equivalently that a weak drive at 4.9 Hz excites the same localized profile without nonlinear mixing. This check is load-bearing for the central claim.
  2. [Section III, Fig. 3] The dispersion relation is obtained by interpolating the checkerboard-sampled angle field with MATLAB's natural interpolation and averaging over 90-degree rotations. Because the central conclusion that the 4.9 Hz mode lies 'under the band' is defined relative to the lower band edge of about 5.7 Hz, a spurious low-frequency feature created by the interpolation mask could affect this comparison. The authors should validate the extraction procedure, for example by computing the Fourier transform directly on the lattice sites without interpolation, or by comparing the numerical dispersion with the linearized equations of motion on the periodic lattice.
  3. [Section IV, Fig. 6] The claim that the onset of the localized mode occurs at approximately B = 0.8 mT is not supported by a quantitative criterion or by error estimates. The spectra in Fig. 6 are single realizations with no stated convergence checks for the ode15s solver tolerances, finite-size effects, or number of drive cycles; a 20 s simulation corresponds to only about 98 cycles at 4.9 Hz, and spectral leakage from the 11.2 Hz drive could contribute to the 4-6 Hz band. Reporting the amplitude of the 4.9 Hz peak as a function of drive amplitude with a defined threshold would make the onset claim falsifiable.
minor comments (5)
  1. [Section IV, first paragraph] The text says the minimum wave propagation frequency was 'determined in Section II,' but the dispersion is presented in Section III; the cross-reference should be corrected.
  2. [Section II, Eq. (1)] The friction coefficient eta is written as eta = 10^-7 without units; specifying the SI units (e.g., N m s or kg m^2 s^-1) would improve reproducibility.
  3. [Introduction, paragraph 3] The statement that 'permanent magnets preclude localized dynamical modes' appears to contradict the paper's main observation of a localized Dirac-string mode; the sentence should be rephrased, for example as 'preclude edge modes,' for internal consistency.
  4. [Fig. 5 caption and text] Please clarify what 'average of modes amplitudes' means quantitatively, including whether the field is root-mean-square amplitude and how the width of the Gaussian convolution was chosen.
  5. [Section IV, Fig. 4] The axis labels and captions should state units explicitly for the magnetic field amplitude (mT) and for the frequency axes (Hz), and the location of the 'output magnet' should be defined in the schematic or text.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the Dirac-string mode at ~4.9 Hz is an emergent numerical observation from a fixed physical model, not a fitted input or a renamed prior result.

full rationale

The paper's derivation chain is: adopt a magnetic-monopole model Eq. (1) with parameters I, η, and q fixed by the physical bar magnets (Ms = 1050 kA/m) and prior work; relax the ground and defective states; compute the dispersion in a 25×25 array; then simulate forced dynamics and observe a spectral peak at ≈4.9 Hz whose spatial profile is localized on the Dirac string. The central claim—that a localized mode exists below the 5.7–14.8 Hz band—is not equivalent to any input. The frequency is obtained from the simulation output, not imposed. The below-band position is checked after the fact as a consistency condition with localization theory, not used to define the mode. The model parameters are physical inputs from a previous experimental/numerical study (Ref. [28]), and that prior work is an externally falsifiable source, so the self-citation is not load-bearing in a circular sense. The paper does not fit any parameter to produce the 4.9 Hz feature, nor does it invoke a uniqueness theorem or rename a known result. The skeptical concern that the 4.9 Hz feature is not shown to be a linear eigenmode (e.g., it could be a combination tone or spectral leakage) is a correctness/evidence question, not a circularity. Under the stated criteria, there is no step where a 'prediction' reduces by construction to its inputs.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The central claim rests on a coarse-grained mechano-magnetic model from prior work by the same group (Refs [26-28]). No new physical entities are introduced; the Dirac string is an established configurational defect. The free parameters are material and damping constants inherited from the experimental macro-ASI, not fitted to the predicted mode.

free parameters (3)
  • effective monopole charge q = 2.08 A m
    Represents the nominal saturation magnetization Ms = 1050 kA/m of the bar magnets; used in all dipolar interactions in Eq. (1). Inherited from Refs [26-28], not fitted to the predicted mode.
  • moment of inertia I = 2.03e-8 kg m^2
    Physical inertia of the 2.54 cm long bar magnets; sets the frequency scale of the dynamics.
  • rotary friction coefficient eta = 1e-7
    Damping coefficient in Eq. (1); inherited from the physical setup of Ref [28]. It affects mode quality and onset of scattering but is not tuned to produce the localized mode.
assumptions (5)
  • domain assumption Magnetic-monopole representation of each bar magnet as two Coulomb-interacting charges (Eq. 1).
    Section II: the torque law replaces each magnet by charges +/-q, ignoring shape anisotropy and internal magnetization dynamics. The central claim is obtained within this approximation.
  • domain assumption Planar rotation constraint (mu_i x k) imposed by the rotary mount.
    Section II: the torque is projected onto this axis, reducing each magnet to a single angular degree of freedom theta_i.
  • domain assumption Stabilization protocol (noise, relaxation, field pulse) reaches the intended ground state or Dirac string state.
    Section II: stable initial states are asserted after this three-step procedure; no energy minimization or uniqueness proof is given.
  • ad hoc to paper Natural interpolation of the checkerboard-sampled theta(x,y,t) yields an accurate dispersion relation.
    Section III: the discrete magnet positions create a mask; the authors use MATLAB natural interpolation and average 90-degree rotated datasets to mitigate artifacts, but this is a numerical assumption.
  • ad hoc to paper The 25x25 extended array and 20 s simulations provide adequate wavevector and frequency resolution.
    Sections III and IV: finite size and finite time set the resolution; no convergence study is reported.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Wave dynamics in a macroscopic square artificial spin ice." pith.science (2026). https://pith.science/paper/HVVGTRN7

@misc{pith2026250109820,
  author       = {Pith},
  title        = {Pith review of: Wave dynamics in a macroscopic square artificial spin ice},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HVVGTRN7}},
  note         = {Machine review of arXiv:2501.09820}
}
read the original abstract

A macroscopic square artificial spin ice, or macro-ASI, is a collection of bar magnets placed in a square lattice arrangement. Each magnet is supported by hinges that allow their mechanical rotation. Previous investigations in these structures have shown ground-state configurations and driven dynamics similar to those of their nanosized counterparts. Here, we numerically investigate the impact of a defect, a Dirac string, on the driven dynamics. We find that waves quickly lose coherence by scattering between modes in this system. In addition, we observe a distinct mode associated with the isolated oscillation of the Dirac string. As expected from spatial localization, this resonant mode is under the band of propagating waves in the macro-ASI. The results are analogous to the development of low-frequency edge modes in nanoscopic spin ices in the presence of defects. Our results provide valuable insight into the physics of mechano-magnetic systems, demonstrating the existence of wave scattering and spatial confinement phenomena in macro-ASIs.

Figures

Figures reproduced from arXiv: 2501.09820 by the authors.

Figure 1
Figure 1. FIG. 1. Schematic of the macro-ASI, with arrows pointing [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Stabilized ASI states for the (a) ground state and [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Dispersion relation for the macro-ASI. The contour [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Frequency spectra of one output magnet for varying [PITH_FULL_IMAGE:figures/full_fig_p003_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. The Dirac-string mode profile as an average of modes [PITH_FULL_IMAGE:figures/full_fig_p004_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Frequency spectra of one output magnet as a function [PITH_FULL_IMAGE:figures/full_fig_p004_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

31 extracted references · 28 canonical work pages

  1. [1]

    S. H. Skjærvø, C. H. Marrows, R. L. Stamps, and L. J. Heyderman, Advances in artificial spin ice, Nat. Rev. Phys. 2, 13 (2020)

  2. [2]

    L. J. Heyderman and R. L. Stamps, Artificial ferroic sys- tems: novelfunctionalityfromstructure, interactionsand dynamics, Journal of Physics: Condensed Matter 25, 363201 (2013)

  3. [3]

    R. Wang, C. Nisoli, R. Freitas, J. Li, W. McConville, B. Cooley, M. Lund, N. Samarth, C. Leighton, V. Crespi, and P. Schiffer, Artificial spin ice in a geometrical frus- trated lattice of nanoscale ferromagnetic islands, Nature 439, 303 (2006)

  4. [4]

    Gilbert, G.-W

    I. Gilbert, G.-W. Chern, S. Zhange, L. O’Brien, B. Foe, C. Nisoli, and P. Schiffer, Emergent ice rule and magnetic charge screening from vertex frustration in artificial spin ice, Nature Physics10, 670 (2014)

  5. [5]

    Gliga, E

    S. Gliga, E. Iacocca, and O. G. Heinonen, Dynamics of reconfigurable artificial spin ice: Toward magnonic functional materials, APL Materials 8, 040911 (2020), https://doi.org/10.1063/1.5142705

  6. [6]

    Lendinez and M

    S. Lendinez and M. B. Jungfleisch, Magnetization dy- namics in artificial spin ice, J. Phys.: Condens. Matter 32, 013001 (2019)

  7. [7]

    Iacocca, S

    E. Iacocca, S. Gliga, R. L. Stamps, and O. Heinonen, Reconfigurablewavebandstructureofanartificialsquare ice, Phys. Rev. B93, 134420 (2016)

  8. [8]

    Iacocca and O

    E. Iacocca and O. Heinonen, Topologically nontriv- ial magnon bands in artificial square spin ices with dzyaloshinskii-moriya interaction, Phys. Rev. Applied8, 034015 (2017)

Show all 31 references
  1. [9]

    Iacocca, S

    E. Iacocca, S. Gliga, and O. G. Heinonen, Tailoring spin- wave channels in a reconfigurable artificial spin ice, Phys. Rev. Applied 13, 044047 (2020)

  2. [10]

    M. B. Jungfleisch, W. Zhang, E. Iacocca, J. Sklenar, J. Ding, W. Jiang, S. Zhang, J. E. Pearson, V. Novosad, J. B. Ketterson, O. Heinonen, and A. Hoffmann, Dy- namic response of an artificial square spin ice, Phys. Rev. B 93, 100401 (2016)

  3. [11]

    V. S. Bhat, F. Heimbach, I. Stasinopoulos, and D. Grundler, Magnetization dynamics of topological de- fects and the spin solid in a kagome artificial spin ice, Phys. Rev. B93, 140401 (2016)

  4. [12]

    Negrello, F

    R. Negrello, F. Montoncello, M. T. Kaffash, M. B. Jungfleisch, and G. Gubbiotti, Dynamic coupling and spin-wave dispersions in a magnetic hybrid system made of an artificial spin-ice structure and an extended nife underlayer, APL Materials10, 091115 (2022)

  5. [13]

    D. M. Arroo, J. C. Gartside, and W. R. Branford, Sculpt- ing the spin-wave response of artificial spin ice via mi- crostate selection, Phys. Rev. B100, 214425 (2019)

  6. [14]

    T. Dion, D. M. Arroo, K. Yamanoi, T. Kimura, J. C. Gartside, L. F. Cohen, H. Kurebayashi, and W. R. Bran- ford, Tunable magnetization dynamics in artificial spin ice via shape anisotropy modification, Phys. Rev. B100, 054433 (2019)

  7. [15]

    J. C. Gartside, A. Vanstone, T. Dion, K. D. Stenning, D. M. Arroo, H. Kurebayashi, and W. R. Branford, Re- configurable magnonic mode-hybridisation and spectral control in a bicomponnt artificial spin ice, Nature Com- munications 12, 2488 (2021)

  8. [16]

    Goryca, X

    M. Goryca, X. Zhang, J. Li, A. L. Balk, J. D. Watts, C. Leighton, C. Nisoli, P. Schiffer, and S. A. Crooker, Field-induced magnetic monopole plasma in artificial spin ice, Phys. Rev. X11, 011042 (2021)

  9. [17]

    Vanstone, J

    A. Vanstone, J. C. Gartside, K. D. Stenning, T. Dion, D. M. Arroo, and W. R. Branford, Spectral fingerprint- 6 ing: microstate readout via remanence ferromagnetic res- onance in artificial spin ice, New Journal of Physics24, 043017 (2022)

  10. [18]

    J. C. Gartside, K. D. Stenning, A. Vanstone, H. H. Holder, D. M. Arroo, T. Dion, F. Caravelli, H. Kure- bayashi, and W. R. Branford, Reconfigurable training and reservoir computing in an artificial spin-vortex ice via spin-wave fingerprinting, Nature Nanotechnology17, 460 (2022)

  11. [19]

    K. D. Stenning, J. C. Gartside, L. Manneschi, C. T. S. Cheung, T. Chen, A. Vanstone, J. Love, H. Holder, F.Caravelli, H.Kurebayashi, K.Everschor-Sitte, E.Vasi- laki, and W. R. Branford, Neuromorphic overparameteri- sation and few-shot learning in multilayer physical neural net...

  12. [20]

    Lendinez, M

    S. Lendinez, M. T. Kaffash, O. G. Heinonen, S. Gliga, E. Iacocca, and M. B. Jungfleisch, Nonlinear multi- magnon scattering in artificial spin ice, Nature Commu- nications 14, 3419 (2023)

  13. [21]

    A. May, M. Saccone, A. van den Berg, J. Askey, M. Hunt, and S. Ladak, Magnetic charge propagation upon a 3d ar- tificial spin ice, Nature Communications12, 3217 (2021)

  14. [22]

    Sahoo, A

    S. Sahoo, A. May, A. van Den Berg, A. K. Mondal, S. Ladak, and A. Barman, Observation of coherent spin waves in a three-dimensional artificial spin ice structure, Nano Letters 21, 4629 (2021)

  15. [23]

    T. Dion, K. D. Stenning, A. Vanstone, H. H. Holder, R. Sultana, G. Alatteili, V. Martinez, M. T. Kaffash, K.T., H.Kurebayashi, W.R.Branford, E.Iacocca, B.M. Jungfleisch, and J. C. Gartside, Ultrastrong magnon- magnon coupling and chiral spin texture control in a dipolar 3d mul...

  16. [24]

    Berchialla, G

    L. Berchialla, G. M. Macauley, and L. J. Hey- derman, Focus on three-dimensional artificial spin ice, Applied Physics Letters 125, 220501 (2024), https://pubs.aip.org/aip/apl/article- pdf/doi/10.1063/5.0229120/20269806/220501_1_5.0229120.pdf

  17. [25]

    Olive and P

    E. Olive and P. Molho, Thermodynamic study of a lattice of compass needles in dipolar interaction, Phys. Rev. B 58, 9238 (1998)

  18. [26]

    Mellado, A

    P. Mellado, A. Concha, and L. Mahadevan, Macro- scopicmagneticfrustration,Phys.Rev.Lett. 109,257203 (2012)

  19. [27]

    Teixeira, M

    H. Teixeira, M. Bernardo, F. Nascimento, M. Saccone, F. Caravelli, C. Nisoli, and C. de Araujo, Macroscopic magnetic monopoles in a 3d-printed mechano-magnet, Journal of Magnetism and Magnetic Materials 596, 171929 (2024)

  20. [28]

    Peroor, L

    R. Peroor, L. Scafuri, D. A. Bozhko, and E. Iacocca, Frequency comb in a macroscopic mechano-magnetic ar- tificial spin ice, , arXiv:2409.13658

  21. [29]

    Gliga, A

    S. Gliga, A. Kákay, R. Hertel, and O. G. Heinonen, Spec- tral analysis of topological defects in an artificial spin-ice lattice, Phys. Rev. Lett.110, 117205 (2013)

  22. [30]

    Venkat, D

    G. Venkat, D. Kumar, M. Franchin, O. Dmytriiev, M. Mruczkiewicz, H. Fangohr, A. Barman, M. Krawczyk, and A. Prabhakar, Proposal for a standard micromag- neticproblem: spinwavedispersioninamagnonicwaveg- uide, IEEE Trans. Magn.49, 524 (2013)

  23. [31]

    Schultheiss, Magnetic domain walls as re- configurable spin-wave nanochannels, Nat

    K.Wagner, A.Kákay, K.Schultheiss, A.Henschke, T.Se- bastian, and H. Schultheiss, Magnetic domain walls as re- configurable spin-wave nanochannels, Nat. Nanotechnol. 11, 432 (2016)

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.