Tracing Monopoles and Anti-monopoles in a Magnetic Hedgehog Lattice

Satoru Hayami; Shun Okumura; Yasuyuki Kato; Yukitoshi Motome

REVIEW 2 major objections 5 minor 1 cited by

Tracing Monopoles and Anti-monopoles in a Magnetic Hedgehog Lattice

T0 review · 2 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read On a discrete lattice, magnetic monopoles and anti-monopoles repel before they annihilate.

desk verdict A solid numerical study with a new repulsion-before-annihilation trajectory, undermined only by the unqualified abstract wording and the lack of a stable-branch check. read the letter →

arxiv 1909.01316 v3 pith:7MGVJ4LF submitted 2019-09-03 cond-mat.str-el

Shun Okumura , Satoru Hayami , Yasuyuki Kato , Yukitoshi Motome This is my paper · ORCID

classification cond-mat.str-el

keywords magnetichedgehoglatticemonopoleanti-monopolescalarspinchiralitysimulatedannealing3Qstatetopologicaltexture

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 asks what happens to the magnetic monopoles and anti-monopoles living in a 3Q hedgehog lattice when a magnetic field is applied, and whether lattice discreteness changes the motion predicted by continuum theory. By simulated annealing with a field sweep on an effective spin model, the authors track each monopole and anti-monopole through the unit cell. They establish that, on the discrete lattice, a monopole and an anti-monopole approaching from different layers do not collide; they repel and return to their original planes before finally annihilating at a higher field. The result matters because real short-period hedgehog lattices are far from the continuum limit, so the annihilation process relevant to experiments could be governed by lattice-scale repulsion rather than direct collision.

What carries the argument

The argument is carried by three tools: an effective spin Hamiltonian with RKKY, biquadratic, and Dzyaloshinskii-Moriya interactions; a local scalar spin chirality defined in vector form that marks the hedgehog cores; and a monopole charge computed from the solid angles of eight spins around each unit cube. The monopole charge, which takes values +1 and -1 when a monopole or anti-monopole occupies the cube, locates the defects unambiguously at interstitial positions. Simulated annealing with a field sweep lets the authors follow the same defects from zero field through the metastable 3Q branch until the total monopole number vanishes, and the trajectories of the Qm = +1 and Qm = -1 positions form the central evidence for the repulsion claim.

What would settle it

Repeat the calculation along the stable ground-state branch obtained by direct energy minimization at each field rather than annealing from the previous field step, and track the positions where Qm equals +1 and -1: if the monopole and anti-monopole that approach from different layers collide instead of repelling before annihilation, the claimed lattice repulsion is an artifact of the metastable sweep. Alternatively, place a single monopole and anti-monopole in an otherwise ordered 3Q background and compute the pair force as a function of separation to see whether the interaction is repulsive at short range.

Watch

Extended reading notes

Core claim

The paper's central discovery is that, in a lattice model of the 3Q hedgehog lattice, monopoles and anti-monopoles move under an increasing [001] field, and those approaching from different layers repel before pair annihilation. This is stated in Section 4 as a difference from the continuum approximation: the defects do not collide with each other but are repelled by other monopoles and anti-monopoles before pair annihilations. The motion is traced by computing the monopole charge in every unit cube, and the trajectories show monopoles shifting to upper layers and anti-monopoles to lower layers while also moving in the xy plane, with a repulsion event around h approximately 0.66 followed by annihilation near h approximately 0.78. The same trajectories connect the motion to the field dependence of the uniform scalar spin chirality: the fictitious fluxes incline toward the field direction, increasing the absolute value of the chirality, and the final annihilation rapidly reduces it. These results are obtained on a metastable 3Q branch accessed by the field sweep, as the paper explicitly notes.

Load-bearing premise

The field sweep follows a metastable 3Q state rather than the stable ground-state branch, and if the repulsion seen near h approximately 0.66 is an artifact of that metastable path, the central claim would not hold for the true equilibrium evolution.

Editorial extensions

If this is right

In short-period hedgehog lattices, pair annihilation is preceded by a repulsion step that changes the defect trajectories, so annihilation is not a simple head-on collision.
The uniform scalar spin chirality grows as fictitious fluxes tilt toward the field direction and drops sharply at annihilation, giving an observable transport signature tied to the defect motion.
The same field-sweep method can track defect motion for other field directions and for other types of hedgehog lattices, so the repulsion mechanism is a candidate general feature of discrete topological spin textures.
At fields just below the transition, the system carries a metastable 3Q state whose monopole arrangement has shifted from the initial layers, which may affect the topological Hall response in field sweeps.

Reading between the lines

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

If the repulsion is a robust lattice effect, continuum descriptions should be corrected with a short-range repulsive interaction between monopoles and anti-monopoles belonging to different layers; a lattice Landau theory with defect-defect interactions could test this.
The metastable branch caveat suggests a direct comparison with the stable ground-state path is needed; if the repulsion disappears there, the phenomenon is a property of the relaxation protocol rather than of the equilibrium 3Q hedgehog lattice.
The relation between defect trajectories and scalar chirality implies that measurements of the topological Hall effect during field sweeps may carry a fingerprint of the repulsion event, such as a two-step change in the chirality before the transition.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

Colleague,

Quick take: this is a legitimate numerical study that finds a genuinely new trajectory for monopoles in a lattice hedgehog lattice—repulsion before pair annihilation—but the result is explicitly obtained on the metastable field-sweep branch, and the abstract presents it without that caveat. That overstatement is the main thing to fix, not a reason to sink the paper.

The new physics is in Fig. 5: as the field increases along [001], monopoles and anti-monopoles move between layers, and at h~0.66 a monopole is repelled by an anti-monopole approaching from a different layer, then the pair annihilates at h~0.78. This contrasts with the continuum calculation of Zhang et al. [6], which predicted collision. The connection to uniform scalar spin chirality is a nice extra, and the monopole detection via solid-angle charge is standard and clearly specified.

What the paper does well: the model and parameters (K=0.7, D=0.3) come from the authors' prior work, so nothing is fitted to the target claim. The methods—simulated annealing, field sweep, solid-angle monopole counting—are standard and described well enough to reproduce. They also deserve credit for explicitly stating in Sec. 2.2 that the field sweep may follow a metastable state beyond first-order transitions, and the concluding remarks say the same. That disclosure is real.

The soft spots, in order of weight:

First, the central claim is only demonstrated on the metastable continuation. The stable ground-state branch is never checked. If a sweep on the stable branch gives colliding/annihilating trajectories like the continuum prediction, then the claimed lattice-discretization effect is actually a property of the metastable continuation, not of the 3Q-HL in equilibrium. The authors acknowledge this possibility indirectly, but they never test it. This makes the abstract's unqualified 'we show' broader than the evidence.

Second, there are no error bars or independent runs; the trajectories are from a single annealing protocol. This is minor for a deterministic-looking result, but it would strengthen the paper to show robustness to annealing schedule or initial conditions.

Third, only one parameter set and one field direction are studied. Fine for a short report, but it limits the generality of the 'repulsive interaction' claim.

Who is this for: people working on topological spin textures, especially lattice effects on monopole dynamics. It is a useful counterpoint to the continuum picture and relevant to interpreting the short-period HL in MnSi1-xGex, though the authors are appropriately cautious about quantitative comparison with experiment.

My recommendation: send it to peer review. It deserves referee time. A good referee should push for a clear statement in the abstract that the observed repulsion occurs on the metastable branch, and ideally a check of the stable branch or at least a discussion of why that branch is not the relevant one for the annihilation process. With that revision, the paper is a solid contribution.

Referee Report

2 major / 5 minor

Summary. This manuscript studies the motion of magnetic monopoles and anti-monopoles in a lattice model of the 3Q magnetic hedgehog lattice, using simulated annealing with a magnetic field sweep along [001]. The authors detect monopole charges via solid-angle computations, trace their positions as a function of field, and report that, contrary to the continuum approximation, the defects move and repel each other before pair annihilation. They also connect the trajectories to the field dependence of the uniform scalar spin chirality. The paper explicitly notes that the field-sweep protocol may follow a metastable branch rather than the equilibrium ground-state sequence.

Significance. If the central claim holds, the paper provides a lattice-discretization correction to the continuum prediction of monopole collision and annihilation, which is directly relevant to short-period hedgehog lattices such as in MnSi1-xGex. The methods are transparent and standard: classical Monte Carlo with simulated annealing, solid-angle monopole detection, and trajectory visualization. The model parameters are taken from prior work rather than fitted to the target claim, and the field dependence is an emergent output. These features make the study a useful and falsifiable numerical contribution, provided the metastable-branch dependence of the central result is properly bounded.

major comments (2)

[Sec. 2.2 / Sec. 3.2 / Fig. 5] The central claim that monopoles and anti-monopoles move and repel before pair annihilation is supported only on the metastable field-sweep branch. Section 2.2 states that the sweep 'may follow a metastable state beyond the first-order phase transitions,' and Section 3.2 notes that the phase boundary differs from the stable ground state obtained in Ref. [7]. The observed repulsion at h≈0.66, including the event where monopole 4 is repelled by anti-monopole 5, could therefore be a property of the metastable continuation rather than a general lattice-discretization effect. The abstract, however, presents the repulsion result without this caveat. To make the conclusion load-bearing, the authors should either compute a stable-branch field sweep (for example, by re-annealing from equilibrium states at each field or by sweeping the field down from the polarized state) or explicitly restrict the abstract and conclusions to the metastable 3Q-HL branch.
[Sec. 3.2 / Fig. 5] The procedure for constructing monopole trajectories is underspecified. The text says only that trajectories are drawn by tracing the positions where Qm=±1; it does not state how monopole charges are matched between successive field steps, how ambiguities are resolved when two charges approach each other near h≈0.66, or whether the identity labels in Fig. 5 remain well-defined under that matching. Because the 'repelled by anti-monopole 5' event is the central evidence, the tracking criterion should be described explicitly and, ideally, validated against alternative matching rules to show that the repulsion is not an artifact of the labeling procedure.

minor comments (5)

[Sec. 2.2] The number of Monte Carlo sweeps per field step is not specified; only the total for the zero-field annealing (10^5–10^6 sweeps) is given. Please state the schedule used at each field value so that the protocol is reproducible.
[Figs. 4 and 5] The results appear to come from a single simulated-annealing run with no statistical uncertainties. Please state explicitly that these are single-run results and comment on the expected run-to-run variation, particularly for the trajectory details near the repulsion event.
[Sec. 2.1] Typo: 'repersents' should be 'represents'.
[Sec. 2.3] Typo: 'amd' should be 'and'.
[Abstract] The abstract should carry the metastable-state caveat that appears in Sec. 2.2 and Sec. 4; currently it states the result as a general property of the lattice system.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the monopole/anti-monopole repulsion is an emergent simulation result, not an input or a self-citation chain.

full rationale

The paper's central claim, that monopoles and anti-monopoles move and repel before pair annihilation, is obtained from simulated annealing of an explicit spin Hamiltonian followed by a field sweep, with the defect positions detected through the monopole charge defined in Eqs. (2)-(5). No parameter is fitted to the claimed repulsion behavior, and no quantity used in the detection is defined in terms of the trajectories or the repulsion outcome. The model parameters K=0.7 and D=0.3 are taken from the authors' prior work (Ref. [7]) to stabilize the 3Q hedgehog lattice, but this self-citation is not load-bearing for the new claim: the zero-field 3Q state used as the starting point is independently reproduced in Sec. 3.1 from the same Hamiltonian, rather than merely imported. The field dependence, the kink in the magnetization, the vanishing of the monopole number, and the trajectories shown in Figs. 4 and 5 are all emergent outputs of the simulation. The acknowledged metastability of the field-sweep protocol (Sec. 2.2: 'in this method we may follow a metastable state beyond the first-order phase transitions') is a scope or correctness caveat about how general the repulsion result is, not a circularity: the prediction is not equivalent by construction to the protocol, and the paper explicitly discloses that the phase boundary differs from the stable ground-state result of Ref. [7]. The comparison with the continuum approximation result [6] is a comparison of two independent calculations, not a renaming of a known result. No circular step can be identified, so the appropriate finding is a non-finding with score 0.

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

The central claim rests on model parameters K, D, and Q chosen from prior work, the metastable field-sweep protocol, and standard monopole detection; no new entities or fitted outputs are introduced.

free parameters (3)

biquadratic coupling K = 0.7
Set by hand to stabilize the 3Q hedgehog lattice at zero field, following Ref. [7]. Not derived from first principles.
Dzyaloshinskii-Moriya coupling D = 0.3
Same as above; chosen with K to reproduce the hedgehog lattice at zero field.
modulation wave vector Q = pi/12
Chooses the 24-site magnetic unit cell to mimic a short-period hedgehog lattice; not derived from material parameters.

assumptions (4)

domain assumption The effective spin Hamiltonian (Eq. 1) with classical unit-length spins faithfully models the 3Q hedgehog lattice in B20 compounds.
The Hamiltonian is taken from perturbation expansions and prior work (Ref. [7]); the paper does not justify it against experimental data.
ad hoc to paper The field sweep follows a metastable state that is the physically relevant branch for observing pair annihilation.
Section 2.2 explicitly acknowledges the phase boundary differs from the stable ground state; the central claim relies on this protocol.
domain assumption The solid-angle monopole charge Qm in Eq. (5) correctly locates monopoles and anti-monopoles on the discrete lattice.
The standard solid-angle definition is used without validation against other detection methods; the trajectories depend on this choice.
domain assumption Periodic boundary conditions with N=24^3 are free of finite-size effects.
The unit cell matches the ordering wave vectors, but no finite-size scaling is checked in the paper.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Tracing Monopoles and Anti-monopoles in a Magnetic Hedgehog Lattice." pith.science (2026). https://pith.science/paper/7MGVJ4LF

@misc{pith2026190901316,
  author       = {Pith},
  title        = {Pith review of: Tracing Monopoles and Anti-monopoles in a Magnetic Hedgehog Lattice},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7MGVJ4LF}},
  note         = {Machine review of arXiv:1909.01316}
}

abstract

The magnetic hedgehog lattice (HL), which was recently discovered in the $B$20-type chiral magnet MnSi$_{1-x}$Ge$_x$, is a topological spin texture with a periodic array of magnetic monopoles and anti-monopoles. Within the continuum approximation, the monopoles and anti-monopoles are predicted to move, collide, and pair annihilate in an applied magnetic field, but it remains unclear how the lattice discretization affects their motions. Here, we study the trajectories of monopoles and anti-monopoles in a lattice system by simulated annealing with field sweep. We show that the monopoles and anti-monopoles move and repel before pair annihilations. We also clarify that their motions are closely related with the field dependence of the scalar spin chirality.

Figures

Figures reproduced from arXiv: 1909.01316 by the authors.

**Figure 2.** Schematic pictures of (a) the local scalar spin chirality χ γ sc(rl) in Eq. (2) and (b) the solid angle Ωα (rc + δα 2 xˆα) in Eq. (3). The gray spheres represent the lattice sites on the cubic lattice, and the blue arrows are the spins at each site. The red and green triangles represent the local scalar spin chirality and the solid angle, respectively. In (b), the gray arrows denote the order of the outer products i… view at source ↗

**Figure 3.** (a) displays a MC snapshot on an xy plane after the simulated annealing. The spin configuration is a noncoplanar one similar to [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

Figures from the paper (2 more)

**Figure 4.** Figure 4: Magnetic field sweep along the [001] direction, starting from the 3Q-HL at zero field: (a) the magnetization m, the uniform scalar spin chirality χsc, the number of monopoles and anti-monopoles Nm, and (b) the magnetic moments with wave vector Qη, mQη . The results ar…

**Figure 5.** Figure 5: Trajectories of monopoles and anti-monopoles in the [001] field: (a) the z positions as functions of the magnetic field h and (b) the projection onto the xy plane. The red (blue) lines represent the trajectories of monopoles (anti-monopoles). The solid squares, diamond…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

Magnetic hedgehog lattices in noncentrosymmetric metals
cond-mat.str-el 2019-08 conditional novelty 7.0 of 10

An effective spin model for noncentrosymmetric metals stabilizes both 3Q and 4Q magnetic hedgehog lattices at zero field when Dzyaloshinskii-Moriya and biquadratic interactions act together.

Reference graph

Works this paper leans on

19 extracted references · 16 canonical work pages · cited by 1 Pith paper

[1]

Kanazawa, J.-H

N. Kanazawa, J.-H. Kim, D. S. Inosov, J. S. White, N. Egete nmeyer, J. L. Gavilano, S. Ishiwata, Y . Onose, T. Arima, B. Keimer, and Y . Tokura, Phys. Rev. B86, 134425 (2012)

work page 2012
[6]

Zhang, A

X.-X. Zhang, A. S. Mishchenko, G. De Filippis, and N. Naga osa, Phys. Rev. B 94, 174428 (2016)

work page 2016
[7]

Okumura, S

S. Okumura, S. Hayami, Y . Kato, and Y . Motome, arXiv:1908 .05044

work page 1908
[2]

Tanigaki, K

T. Tanigaki, K. Shibata, N. Kanazawa, X. Y u, Y . Onose, H. S . Park, D. Shindo, and Y . Tokura, Nano Lett. 15, 5438 (2015)

work page 2015
[3]

Fujishiro, N

Y . Fujishiro, N. Kanazawa, T. Nakajima, X. Z. Y u, K. Ohishi, Y . Kawamura, K. Kakurai, T. Arima, H. Mi- tamura, A. Miyake, K. Akiba, M. Tokunaga, A. Matsuo, K. Kindo , T. Koretsune, R. Arita, and Y . Tokura, Nat. Commun. 10, 1059 (2019)

work page 2019
[4]

Kanazawa, Y

N. Kanazawa, Y . Nii, X. X. Zhang, A. S. Mishchenko, G. De Fi lippis, F. Kagawa, Y . Iwasa, N. Nagaosa, and Y . Tokura, Nat. Commun.7, 11622 (2016)

work page 2016
[5]

Kanazawa, Y

N. Kanazawa, Y . Onose, T. Arima, D. Okuyama, K. Ohoyama, S . Wakimoto, K. Kakurai, S. Ishiwata, and Y . Tokura, Phys. Rev. Lett.106, 156603 (2011)

work page 2011
[8]

M. A. Ruderman and C. Kittel, Phys. Rev. 96, 99 (1954)

1954

Show all 19 references

[9]

Kasuya, Prog

T. Kasuya, Prog. Theor. Phys. 16, 45 (1956)

1956
[10]

Y osida, Phys

K. Y osida, Phys. Rev. 106, 893 (1957)

1957
[11]

Akagi, M

Y . Akagi, M. Udagawa, and Y . Motome, Phys. Rev. Lett. 108, 096401 (2012)

2012
[12]

Hayami and Y

S. Hayami and Y . Motome, Phys. Rev. B 90, 060402 (2014)

2014
[13]

Hayami, R

S. Hayami, R. Ozawa, and Y . Motome, Phys. Rev. B 95, 224424 (2017)

2017
[14]

Dzyaloshinsky, J

I. Dzyaloshinsky, J. Phys. Chem. Solids 4, 241 (1958)

1958
[15]

Moriya, Phys

T. Moriya, Phys. Rev. 120, 91 (1960)

1960
[16]

Hayami and Y

S. Hayami and Y . Motome, Phys. Rev. Lett. 121, 137202 (2018)

2018
[17]

Y ang, Y .-H

S.-G. Y ang, Y .-H. Liu, and J. H. Han, Phys. Rev. B 94, 054420 (2016)

2016
[18]

Shiomi, N

Y . Shiomi, N. Kanazawa, K. Shibata, Y . Onose, and Y . Tokura, Phys. Rev. B 88, 064409 (2013)

2013
[19]

Fujishiro, N

Y . Fujishiro, N. Kanazawa, T. Shimojima, A. Nakamura, K . Ishizaka, T. Koretsune, R. Arita, A. Miyake, H. Mitamura, K. Akiba, M. Tokunaga, J. Shiogai, S. Kimura, S. A waji, A. Tsukazaki, A. Kikkawa, Y . Taguchi, and Y . Tokura, Nat. Commun.9, 408 (2018). 6

2018

Pith tools

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

[1] [1]

Kanazawa, J.-H

N. Kanazawa, J.-H. Kim, D. S. Inosov, J. S. White, N. Egete nmeyer, J. L. Gavilano, S. Ishiwata, Y . Onose, T. Arima, B. Keimer, and Y . Tokura, Phys. Rev. B86, 134425 (2012)

work page 2012

[2] [6]

Zhang, A

X.-X. Zhang, A. S. Mishchenko, G. De Filippis, and N. Naga osa, Phys. Rev. B 94, 174428 (2016)

work page 2016

[3] [7]

Okumura, S

S. Okumura, S. Hayami, Y . Kato, and Y . Motome, arXiv:1908 .05044

work page 1908

[4] [2]

Tanigaki, K

T. Tanigaki, K. Shibata, N. Kanazawa, X. Y u, Y . Onose, H. S . Park, D. Shindo, and Y . Tokura, Nano Lett. 15, 5438 (2015)

work page 2015

[5] [3]

Fujishiro, N

Y . Fujishiro, N. Kanazawa, T. Nakajima, X. Z. Y u, K. Ohishi, Y . Kawamura, K. Kakurai, T. Arima, H. Mi- tamura, A. Miyake, K. Akiba, M. Tokunaga, A. Matsuo, K. Kindo , T. Koretsune, R. Arita, and Y . Tokura, Nat. Commun. 10, 1059 (2019)

work page 2019

[6] [4]

Kanazawa, Y

N. Kanazawa, Y . Nii, X. X. Zhang, A. S. Mishchenko, G. De Fi lippis, F. Kagawa, Y . Iwasa, N. Nagaosa, and Y . Tokura, Nat. Commun.7, 11622 (2016)

work page 2016

[7] [5]

Kanazawa, Y

N. Kanazawa, Y . Onose, T. Arima, D. Okuyama, K. Ohoyama, S . Wakimoto, K. Kakurai, S. Ishiwata, and Y . Tokura, Phys. Rev. Lett.106, 156603 (2011)

work page 2011

[8] [8]

M. A. Ruderman and C. Kittel, Phys. Rev. 96, 99 (1954)

1954