REVIEW 3 major objections 6 minor 2 cited by
Excitation of vortex core gyration in nanopillars through driven Floquet magnons
T0 review · 3 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read The paper argues that a single high-frequency RF drive can sustain not one but several distinct steady-state gyration orbits of a magnetic vortex core, because the core's nonlinear interactions with driven Floquet magnons create multiple…
desk verdict A real extension of the self-induced Floquet magnon work to nanopillars, with a solid fixed-point analysis and a well-supported multistability claim; the main weakness is an untested linear-scaling assumption in the model. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the extended Thiele equation with a nonlinear interaction contribution f(R) added to the radial motion. The Thiele model treats the vortex core position as the only collective coordinate; the authors augment its Lagrangian with Berry-phase coupling terms describing three-particle processes |μ±1⟩ ↔ |μ⟩ ± |g⟩, which convert a velocity operator of the gyration into forces on the magnon amplitudes. When substituted into the polar-coordinate Thiele equation, these terms generate the self-consistent radial equation dR/dt = -ΓgR + f(R), whose fixed-point structure—where the NLC balances Gilbert damping—carries the entire multi-stability and bifurcation phenomenology.
What would settle it
Simulate or measure f(R) at several intermediate drive amplitudes across the predicted bistable windows (e.g., between 2.5 and 2.9 mT at 11 GHz), or perform up-and-down sweeps of brf at fixed frf and check whether the gyration radius and frequency comb jump at the predicted thresholds; if f(R) does not scale linearly in brf or the expected jumps are absent, the central multi-stability claim is refuted.
Extended reading notes
Core claim
In a 300 nm diameter, 20 nm thick ferromagnetic disk with a vortex ground state, driving one of the azimuthal (n = 0, m = ±1) spin-wave modes with an in-plane RF field produces a frequency comb spaced by the gyration frequency once the field amplitude exceeds a threshold. The paper shows that the radial dynamics of the core can be reduced to dR/dt = -ΓgR + f(R), where f(R) is a nonlinear interaction contribution from scattering between the core and Floquet magnon modes. Depending on the drive frequency, f(R) intersects the relaxation line ΓgR at multiple points, yielding up to three simultaneously stable gyration radii at a single drive (for example at frf = 12.29 GHz and brf = 4.7 mT). Each stable radius produces a distinctly different Floquet frequency comb, so the system is hysteretic under field or frequency sweeps.
Load-bearing premise
The predictions rely on the assumption that the nonlinear interaction contribution f(R) grows linearly with the radio-frequency drive amplitude, so the bifurcation curves drawn from measured f(R) at two amplitudes apply at all other amplitudes.
Editorial extensions
If this is right
- Upward and downward sweeps of the RF field amplitude (or frequency) should reveal hysteresis in the gyration radius, since the stable fixed point reached depends on the initial core position.
- The frequency-comb spectrum becomes a history-dependent readout: different steady orbits produce measurably distinct comb shapes, so the comb can be used to detect transitions between metastable gyration radii.
- The self-induced mechanism is generic to vortex-based disks: it appears in both the 300 nm nanopillars studied here and the larger 5 µm disks of earlier work, with the required drive frequency set by the azimuthal-mode splitting.
- The bifurcation scenarios include supercritical Hopf, saddle-node, and multi-fixed-point types, so vortex nanopillars constitute a small laboratory for Floquet-driven auto-oscillator physics.
Reading between the lines
- The predicted thresholds and coexistence intervals depend on the assumption that f(R) scales linearly with brf; a direct measurement of f(R) at intermediate field amplitudes (e.g., brf = 0.9, 1.5, 2.0 mT at frf = 11 GHz) would confirm or refute the entire bifurcation diagram.
- The three-particle scattering picture implies an angular-momentum bookkeeping that could be tested by measuring the m-resolved comb: the checkerboard pattern seen in S_m(f) is a direct map of azimuthal-number flow between the core and magnon modes.
- Because the full Floquet comb carries history-dependent structure, it could serve as a richer output channel for magnon-scattering reservoir computers than the single scattered modes used in current proposals.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports a computational and theoretical study of self-induced Floquet magnons in a 300 nm × 20 nm vortex-state nanopillar. Under an RF field driving azimuthal spin-wave modes, the vortex core can enter steady gyration through nonlinear scattering; the authors extend the Thiele model with three-particle scattering terms to derive a radial equation dR/dt = -Gamma_g R + f(R) and use f(R) extracted from simulated Floquet spectra to construct bifurcation diagrams. For several drive frequencies they find parameter windows with multiple stable fixed-point radii, and direct simulations with different initial radii confirm coexisting steady states, which yield distinct frequency-comb spectra and imply hysteretic behavior.
Significance. The qualitative central claim—that one RF drive can sustain multiple metastable gyration radii with observable spectral signatures—is well supported by the direct micromagnetic simulations in Figs. 4 and 5, and the use of independent initial radii to test fixed-point predictions is a definite strength. If it holds, the result provides an experimentally accessible hysteresis mechanism and extends the self-induced Floquet magnon effect from micron-sized disks to nanopillar geometries. The paper is clearly written and uses standard open-source simulation tools. However, the quantitative theory is partly phenomenological: f(R) is fitted to simulated Floquet spectra, Eq. (15) is not derived in the manuscript, and the linear-in-brf scaling of f(R) is assumed rather than tested. These issues limit the predictive power of the model but do not undermine the core multistability finding, which stands on the simulations.
major comments (3)
- [Sec. III.D, Eqs. (16)-(17) and Fig. 4] The linear-in-brf scaling of f(R) is a load-bearing assumption that is not justified. The text says "We assume f(R) to scale linearly with brf, based on arguments leading to Eq. (15)," but Eq. (15) only isolates one factor proportional to P0 brf / Gamma0; the Floquet comb amplitudes a_tilde_k in Eq. (16) are themselves produced by the same drive through the nonlinear scattering of Eq. (11) and can carry additional brf dependence. The bifurcation diagrams in Fig. 4(b,d,f,h) are constructed from f(R) curves measured at only two brf values per frequency and linearly rescaled up to 5.5 mT; if the scaling fails, the predicted thresholds and coexistence intervals will shift or disappear. The authors should measure f(R) at several intermediate driving amplitudes, or derive its brf dependence, before the quantitative bifurcation predictions can be accepted. The qualitative multistability is nonetheless independently supported by the direct simulations in Fig. 5.
- [Sec. III.C, Eq. (15)] The derivation of Eq. (15) is not shown and is delegated to the unpublished Ref. [27]; the numerical extraction of f(R) from "the five sidebands" is also not specified. Because Eq. (17) is the central dynamical equation, the manuscript should be self-contained: either provide the derivation of Eq. (15) in the main text or an appendix, or cite a published version of Ref. [27], and give the explicit formula used to convert sideband amplitudes into f(R). In addition, since f(R) is extracted from the same class of micromagnetic simulations used for validation, the "good quantitative agreement" in Fig. 4 demonstrates internal consistency of the reduction rather than an independent prediction; the genuinely predictive element is the bifurcation structure, which is tested via different initial radii. Please state this limitation explicitly.
- [Sec. III.C and Discussion] The assumption that Floquet modes preserve the azimuthal indices of the linear modes is asserted rather than demonstrated, and the scattering coefficients C_k are treated as fitting parameters. The manuscript acknowledges this, but it means the theory is a phenomenological parametrization of the simulations rather than a first-principles derivation. The authors should either verify the azimuthal-index conservation by projecting the simulated Floquet profiles onto the linear eigenmodes, or soften the claim in the Discussion that the model "accurately predict[s]" the fixed-point structure. The central multistability finding does not depend on this assumption alone, because Fig. 5 provides direct simulation evidence.
minor comments (6)
- [Abstract and Sec. I] The word "Di fferent" appears with an extra space in the abstract (and in the typeset text); correct the spacing.
- [Sec. II] The phrase "allowed as to obtain" should read "allowed us to obtain."
- [Sec. III.A] The text says "two broader peaks inTg can also be seen," but the quantity plotted in Fig. 1(c) is Rg, the gyration radius; this appears to be a typo for "in Rg."
- [Sec. III.B] The simulations for Fig. 3 are said to use a finite temperature of 1 K, whereas the rest of the simulations are at zero temperature; clarify why this temperature is introduced and whether it affects the threshold comparison.
- [Sec. III.B, Fig. 3(b)] The comb lines are defined as f = fg(m - m0), but m0 is not defined; specify the offset for the central comb line.
- [Sec. IV] The phrase "we can accurately predict the existence of multiple stable and unstable gyration radii" overstates the role of the fitted f(R); consider "reproduce" or "capture" instead of "predict."
Circularity Check
Bifurcation diagrams are read off from fitted f(R), so the quantitative prediction reduces to a simulation-derived fit; qualitative multistability is independently supported by direct simulations.
-
fitted input called prediction
[Section III.D, Eqs. (16)-(17), Fig. 4; Discussion]
"We assume f (R) to scale linearly with brf, based on arguments leading to Eq. (15). ... By parameterizing the model parameters with simulated Floquet spectra at different steady-state gyration radii, we can accurately predict the existence of multiple stable and unstable gyration radii under azimuthal mode pumping, with good quantitative agreement found with micromagnetics simulations."
The central quantitative output, R*(brf) and the multistability intervals in Fig. 4, is obtained by solving Eq. (17), dR/dt = -Gamma_g R + f(R), with f(R) numerically estimated from simulated Floquet spectra using the five sidebands. The paper explicitly states that scattering coefficients are fitting parameters for the simulated Floquet spectra. Thus the predicted fixed points are the crossings of the damping line with a fitted curve; the threshold and bifurcation diagrams are read off from that same fitted input, linearly rescaled in brf. This is partial circularity, not total: the direct micromagnetics simulations in Fig. 5 and the open-symbol tests in Fig. 4 independently support the qualitative multi-radius effect.
-
self citation load bearing
[Section III.C, Eq. (11)]
"To second order in the fluctuations, it can be shown [27] that additional terms such as L(2)B ~ sum_{n,m} i C_{n,m} (a*_{n,m} a_{n,m+1} V+_g - a_{n,m} a*_{n,m+1} V-_g) appear in the Berry-phase term."
Ref. [27] is an unpublished arXiv preprint whose author list overlaps with the present paper (J.-V. Kim). The three-particle scattering structure that motivates Eq. (16) and the definition of f(R) is imported from that citation rather than derived in this paper, and the text extends it by assuming scattering terms such as Eq. (11) remain valid for Floquet modes. The theoretical skeleton is therefore anchored in a self-citation; this is a justification gap, though the empirical claim does not rest solely on this citation because the coefficients are fitted and direct simulations provide independent support.
full rationale
The model's predictive content is carried by f(R), which is extracted from simulated Floquet spectra of the same physical system and then used in Eq. (17) to generate the fixed-point curves and bifurcation diagrams. The paper honestly calls the scattering coefficients fitting parameters, but it also calls the resulting fixed points predictions. That is a fitted-input-as-prediction pattern: the quantitative thresholds and multistability intervals in Fig. 4 are graphical consequences of a fitted f(R), not independent first-principles results. The untested linear-in-brf scaling of f(R) is an additional correctness risk, but it is not itself circularity. Separately, the derivation of the scattering terms is delegated to Ref. [27], an overlapping-author preprint, making the theoretical chain partly self-citational. The paper does retain independent content: direct LLG simulations (Fig. 5) show three distinct steady-state radii under the same drive, and the open-symbol checks in Fig. 4 test the predictions out-of-sample relative to the spectral fits. These independent simulations keep the central qualitative claim from being merely tautological, but the quantitative bifurcation predictions still reduce largely to the fitted input. Score 6 reflects this partial circularity.
Assumptions & free parameters
free parameters (3)
- Nonlinear interaction contribution f(R) =
Extracted from 5 comb sidebands at each frf; not given as analytic form
- Scattering coefficients C_k (k=1..N) in Eq. (16) =
Not reported numerically
- Linear scaling factor of f(R) with brf =
Implicit
assumptions (6)
- domain assumption Thiele rigid-core ansatz: vortex core position X(t) is the only dynamical variable
- ad hoc to paper Floquet modes preserve azimuthal indices of linear modes, so the three-particle scattering form Eq. (11) still applies
- ad hoc to paper f(R) is linear in brf
- domain assumption Steady-state gyration is reached before Floquet spectra are extracted
- domain assumption Confining potential is quartic, U = kappa1 |X|^2 + kappa2 |X|^4
- domain assumption Micromagnetic Landau-Lifshitz-Gilbert equation accurately captures vortex and spin-wave dynamics at this scale
Cite this review
Pith. "Pith review of Excitation of vortex core gyration in nanopillars through driven Floquet magnons." pith.science (2026). https://pith.science/paper/IOALXWSH
@misc{pith2026250719865,
author = {Pith},
title = {Pith review of: Excitation of vortex core gyration in nanopillars through driven Floquet magnons},
year = {2026},
howpublished = {\url{https://pith.science/paper/IOALXWSH}},
note = {Machine review of arXiv:2507.19865}
}
read the original abstract
The dynamics of vortex states in confined geometries like thin-film disks are characterized by a sub-GHz gyration, representing the damped oscillatory motion of the vortex core about the disk center. It has recently been shown that interactions between the core and azimuthal spin waves, lying in the GHz range and driven by magnetic fields, can result in steady-state core gyration. The gyration in turn provides a time-periodic modulation for the spin waves, resulting in the emergence of Floquet states. Here, we present results of a theoretical and computational study in which we examine how Floquet modes sustain this core gyration. In particular, we find that multiple steady-state gyration radii are possible under certain field conditions, resulting from the nonlinear interactions between the core and Floquet modes. Different gyration radii result in distinct Floquet frequency comb spectra and allow for hysteretic effects, as reported in recent experiments.
Figures
Forward citations
Cited by 2 Pith papers
-
Coupling phase interference effects in a multimode cavity magnonics system
Coupling phases—not just strengths—determine which cavity modes couple to magnons and can produce nonreciprocal transmission at antiresonances in a multimode cavity magnonics system.
-
Time-resolved observation of magnon splitting into vortex gyration and Floquet spin waves
Time-resolved single-shot electrical measurements show the vortex gyration and the first Floquet sideband of a magnon frequency comb emerge synchronously after an incubation delay, identifying three-wave splitting of ...
Reference graph
Works this paper leans on
- [27]
-
[1]
K. Y . Guslienko and K. L. Metlov, Evolution and stability of a magnetic vortex in a small cylindrical ferromagnetic particle under applied field, Phys. Rev. B 63, 100403 (2001)
work page 2001
-
[2]
K. L. Metlov and K. Y . Guslienko, Stability of magnetic vortex in soft magnetic nano-sized circular cylinder, Journal of Mag- netism and Magnetic Materials 242-245, 1015 (2002)
work page 2002
- [3]
-
[4]
K. Y . Guslienko and V . Novosad, V ortex state stability in soft magnetic cylindrical nanodots, Journal of Applied Physics 96, 4451 (2004)
work page 2004
-
[5]
K. Y . Guslienko, B. A. Ivanov, V . Novosad, Y . Otani, H. Shima, and K. Fukamichi, Eigenfrequencies of vortex state excitations in magnetic submicron-size disks, Journal of Applied Physics 91, 8037 (2002)
work page 2002
-
[6]
B. A. Ivanov and C. E. Zaspel, Gyrotropic mode frequency of vortex-state permalloy disks, Journal of Applied Physics 95, 7444 (2004)
work page 2004
-
[7]
J. Ding, G. N. Kakazei, X. Liu, K. Y . Guslienko, and A. O. Adeyeye, Higher order vortex gyrotropic modes in circular fer- romagnetic nanodots, Scientific Reports 4, 4796 (2014)
work page 2014
Show all 41 references
-
[8]
A. V . Bondarenko, S. A. Bunyaev, A. K. Shukla, A. Apolinario, N. Singh, D. Navas, K. Y . Guslienko, A. O. Adeyeye, and G. N. Kakazei, Dominant higher-order vortex gyromodes in circular magnetic nanodots, Nanoscale Horizons 9, 1498 (2024)
2024
-
[9]
Novosad, M
V . Novosad, M. Grimsditch, K. Guslienko, P. Vavassori, Y . Otani, and S. Bader, Spin excitations of magnetic vortices in ferromagnetic nanodots, Physical Review B66, 052407 (2002)
2002
-
[10]
B. A. Ivanov and C. E. Zaspel, Magnon modes for thin circular vortex-state magnetic dots, Applied Physics Letters 81, 1261 (2002)
2002
-
[11]
J. Park, P. Eames, D. Engebretson, J. Berezovsky, and P. Crow- ell, Imaging of spin dynamics in closure domain and vortex structures, Physical Review B 67, 020403 (2003)
2003
-
[12]
Buess, T
M. Buess, T. Knowles, R. H ¨ollinger, T. Haug, U. Krey, D. Weiss, D. Pescia, M. Scheinfein, and C. Back, Excitations with negative dispersion in a spin vortex, Physical Review B 71, 104415 (2005)
2005
-
[13]
Ivanov and C
B. Ivanov and C. Zaspel, High Frequency Modes in V ortex- State Nanomagnets, Physical Review Letters 94, 027205 (2005)
2005
-
[14]
A. A. Awad, K. Y . Guslienko, J. F. Sierra, G. N. Kakazei, V . Metlushko, and F. G. Aliev, Precise probing spin wave mode frequencies in the vortex state of circular magnetic dots, Ap- plied Physics Letters 96, 012503 (2010)
2010
-
[15]
Castel, J
V . Castel, J. B. Youssef, F. Boust, R. Weil, B. Pigeau, G. D. Loubens, V . Naletov, O. Klein, and N. Vukadinovic, Perpendic- ular ferromagnetic resonance in soft cylindrical elements: V or- tex and saturated states, Physical Review B 85, 184419 (2012)
2012
-
[16]
Taurel, T
B. Taurel, T. Valet, V . V . Naletov, N. Vukadinovic, G. D. Loubens, and O. Klein, Complete mapping of the spin-wave spectrum in a vortex-state nanodisk, Physical Review B 93, 184427 (2016)
2016
-
[17]
Schultheiss, R
K. Schultheiss, R. Verba, F. Wehrmann, K. Wagner, L. K¨orber, T. Hula, T. Hache, A. Kakay, A. A. Awad, V . Tiberkevich, A. N. Slavin, J. Fassbender, and H. Schultheiss, Excitation of Whis- pering Gallery Magnons in a Magnetic V ortex, Physical Review Letters 122, 097202 (2019)
2019
-
[18]
J. P. Park and P. A. Crowell, Interactions of Spin Waves with a Magnetic V ortex, Physical Review Letters95, 167201 (2005)
2005
-
[19]
K. Y . Guslienko, A. N. Slavin, V . Tiberkevich, and S.-K. Kim, Dynamic Origin of Azimuthal Modes Splitting in V ortex-State Magnetic Dots, Physical Review Letters 101, 247203 (2008), 0808.3335
2008 arXiv
-
[20]
K ¨orber, I
L. K ¨orber, I. K´ezsm´arki, and A. K ´akay, Mode splitting of spin waves in magnetic nanotubes with discrete symmetries, Physi- cal Review B 105, 184435 (2022)
2022
-
[21]
Uzunova, L
V . Uzunova, L. K ¨orber, A. Kavvadia, G. Quasebarth, H. Schultheiss, A. K´akay, and B. Ivanov, Nontrivial Aharonov- Bohm e ffect and alternating dispersion of magnons in cone- state ferromagnetic rings, Physical Review B 108, 174445 (2023), 2308.08486
2023 arXiv
-
[22]
Salama, J.-V
S. Salama, J.-V . Kim, A. Anane, and J.-P. Adam, Large fre- quency nonreciprocity of azimuthal spin-wave modes in sub- micron vortex state disks, Phys. Rev. B111, 134445 (2025)
2025
-
[23]
K ¨orber, C
L. K ¨orber, C. Heins, I. Soldatov, R. Sch ¨afer, A. K ´akay, H. Schultheiss, and K. Schultheiss, Modification of three- magnon splitting in a flexed magnetic vortex, Applied Physics Letters 122, 092401 (2023)
2023
-
[24]
K. Y . Guslienko, G. R. Aranda, and J. M. Gonzalez, Topolog- ical gauge field in nanomagnets: Spin-wave excitations over a slowly moving magnetization background, Physical Review B 81, 014414 (2010)
2010
-
[25]
Oka and S
T. Oka and S. Kitamura, Floquet Engineering of Quantum Ma- terials, Annual Review of Condensed Matter Physics 10, 387 (2019)
2019
-
[26]
Rodriguez-Vega, M
M. Rodriguez-Vega, M. V ogl, and G. A. Fiete, Low-frequency and Moir ´e–Floquet engineering: A review, Annals of Physics 435, 168434 (2021)
2021
-
[28]
Heins, A
C. Heins, A. K ´akay, J.-V . Kim, G. Hlawacek, J. Fassbender, K. Schultheiss, and H. Schultheiss, Control of magnon fre- quency combs in magnetic rings, (2025), arXiv:2501.05080 [cond-mat.mes-hall]. 9
2025 arXiv
-
[29]
Martins, A
L. Martins, A. S. Jenkins, J. Borme, J. Ventura, P. P. Freitas, and R. Ferreira, Second harmonic injection locking of coupled spin torque vortex oscillators with an individual phase access, Communications Physics 6, 72 (2023)
2023
-
[30]
Bruckner, S
F. Bruckner, S. Koraltan, C. Abert, and D. Suess, magnum.np: a pytorch based gpu enhanced finite di fference micromagnetic simulation framework for high level development and inverse design, Scientific Reports 13, 12054 (2023)
2023
-
[31]
Vansteenkiste, J
A. Vansteenkiste, J. Leliaert, M. Dvornik, M. Helsen, F. Garcia- Sanchez, and B. Van Waeyenberge, The design and verification of MuMax3, AIP Advances 4, 107133 (2014)
2014
-
[32]
Massouras, S
M. Massouras, S. Perna, M. D’Aquino, C. Serpico, and J.-V . Kim, Mode-resolved micromagnetics study of parametric spin wave excitation in thin-film disks, Phys. Rev. B 110, 064435 (2024)
2024
-
[33]
Ivanov, H
B. Ivanov, H. Schnitzer, F. Mertens, and G. Wysin, Magnon modes and magnon-vortex scattering in two-dimensional easy- plane ferromagnets, Physical Review B 58, 8464 (1998)
1998
-
[34]
B. A. Ivanov, G. G. Avanesyan, A. V . Khvalkovskiy, N. E. Kula- gin, C. E. Zaspel, and K. A. Zvezdin, Non-Newtonian dynamics of the fast motion of a magnetic vortex, JETP Letters 91, 178 (2010)
2010
-
[35]
A. A. Thiele, Steady-State Motion of Magnetic Domains, Phys- ical Review Letters 30, 230 (1973)
1973
-
[36]
Ivanov and C
B. Ivanov and C. Zaspel, Excitation of Spin Dynamics by Spin- Polarized Current in V ortex State Magnetic Disks, Physical Re- view Letters 99, 247208 (2007)
2007
-
[37]
Rajaraman, Solitons and Instantons (North Holland, Ams- terdam, 1982)
R. Rajaraman, Solitons and Instantons (North Holland, Ams- terdam, 1982)
1982
-
[38]
Kim, Spin-torque oscillators, inSolid State Physics, edited by R
J.-V . Kim, Spin-torque oscillators, inSolid State Physics, edited by R. E. Camley and R. L. Stamps (Academic Press, 2012) pp. 217–294
2012
-
[39]
K ¨orber, C
L. K ¨orber, C. Heins, T. Hula, J.-V . Kim, S. Thlang, H. Schultheiss, J. Fassbender, and K. Schultheiss, Pattern recognition in reciprocal space with a magnon-scattering reser- voir, Nature Communications 14, 3954 (2023)
2023
-
[40]
Heins, J.-V
C. Heins, J.-V . Kim, L. K ¨orber, J. Fassbender, H. Schultheiss, and K. Schultheiss, Benchmarking a magnon-scattering reser- voir with modal and temporal multiplexing, Physical Review Applied 23, 054087 (2025)
2025
-
[41]
Butschek, A
L. Butschek, A. Akrout, E. Dimitriadou, A. Lupo, M. Hael- terman, and S. Massar, Photonic reservoir computer based on frequency multiplexing, Optics Letters 47, 782 (2022)
2022
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.