{"id":"020efe4a-72ff-419a-8d5c-c212289f3e97","arxiv_id":"2507.19865","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"RF-driven azimuthal spin waves in a 300 nm vortex nanopillar support multiple steady-state gyration radii, each producing a distinct Floquet frequency comb, so the device can be hysteretic.","lead":"Micromagnetic simulations and a collective-coordinate model show that radio-frequency driving of azimuthal spin waves can hold a vortex core in several different steady gyration orbits in a 300 nm nanopillar. Since each orbit emits a different frequency comb, the output spectrum depends on the drive's history, which could be useful in magnonic signal processing.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The bifurcation diagrams of Fig. 4 are built on an untested linear-in-brf scaling of f(R) (Sec. III.D, Eqs. 16-17), so the predicted thresholds and multistability intervals are not yet quantitatively established; the qualitative effect is independently supported by direct simulations.","rationale":"The reader's weakest-assumption analysis identifies the same load-bearing concern: the assertion that f(R) scales linearly with brf, used to construct all bifurcation diagrams. I agree with that identification and with the moderate conditionality of the verdict. The concern is genuine because Eqs. (15)-(16) contain comb amplitudes that are field-driven, so no first-principles argument establishes pure linear scaling; the only offered basis is an analogy to Eq. (15), which does not by itself exclude higher-order brf dependence. The two measured f(R) curves per frequency are insufficient to validate rescaling over the factor-of-ten amplitude range explored in Fig. 4. At the same time, the central qualitative claim of multiple coexisting gyration radii is not solely a model artifact: Fig. 5 shows direct simulations with three different initial radii relaxing to three distinct steady-state radii under identical drive, and the predicted stable branches in Fig. 4 are checked against independent simulations with different initial conditions. Those checks give the claim independent support even if the quantitative scaling assumption later proves to be incomplete. I therefore would not move the verdict to accept or reject; the appropriate action is to keep the conditional verdict pending a direct test of the linear-scaling assumption and a fuller description of the f(R) extraction procedure.","tokens_in":13753,"tokens_out":7443,"duration_ms":93265,"concrete_test":"Use the two f(R) curves already shown in Fig. 4(a) for frf = 11.5 GHz: for each R in the overlap region, compute q(R) = f(R, brf = 1.3 mT)/1.3 - f(R, brf = 0.5 mT)/0.5. If max|q(R)| exceeds the estimated numerical noise (for example, 10% of max|f|/mT), the linearity assumption is falsified. Then extend the same collapse test to frf = 12.29 GHz by extracting f(R) at brf = 2.0, 3.2, and 4.7 mT; if f(R, brf)/brf does not collapse onto a single curve, Fig. 4(d,f,h) must be recomputed with the measured brf-dependent f(R) and the multistability intervals re-determined.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative step is Eq. (17), dR/dt = -Gamma_g R + f(R), with f(R) defined in Eq. (16) as a sum of products of Floquet comb amplitudes. In Section III.D the authors state: 'We assume f(R) to scale linearly with brf, based on arguments leading to Eq. (15).' Eq. (15), however, only isolates one factor proportional to P0 brf/Gamma_0; the adjacent comb amplitudes a_tilde_k entering Eq. (16) are themselves driven by the same rf field through nonlinear scattering, so the true brf dependence of f(R) may be quadratic or higher. The bifurcation diagrams in Fig. 4(b,d,f,h) are generated by rescaling f(R) curves measured at only two brf values across the full range of brf, up to 5.5 mT. If linear scaling fails, the predicted thresholds and, more importantly, the predicted coexistence intervals of multiple stable R*_g will move or disappear. The extraction procedure is also not specified analytically ('we obtain numerical estimates of f(R) by using the five sidebands'), and the underlying justification is delegated to unpublished Ref. [27]. This combination makes the linear-scaling assumption the most load-bearing element of the theoretical analysis. It should be stressed that the qualitative multistability claim does not rest on this assumption alone: the direct micromagnetic simulations in Fig. 5 show three distinct steady-state radii for the same drive, which independently supports the central claim. The concern therefore weakens the quantitative predictive power of the model and the specific bifurcation curves, not the existence of the effect.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":14130,"tokens_out":8744,"duration_ms":86571,"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":[{"comment":"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.","section":"Sec. III.D, Eqs. (16)-(17) and Fig. 4"},{"comment":"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.","section":"Sec. III.C, Eq. (15)"},{"comment":"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.","section":"Sec. III.C and Discussion"}],"minor_comments":[{"comment":"The word \"Di fferent\" appears with an extra space in the abstract (and in the typeset text); correct the spacing.","section":"Abstract and Sec. I"},{"comment":"The phrase \"allowed as to obtain\" should read \"allowed us to obtain.\"","section":"Sec. II"},{"comment":"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.\"","section":"Sec. III.A"},{"comment":"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.","section":"Sec. III.B"},{"comment":"The comb lines are defined as f = fg(m - m0), but m0 is not defined; specify the offset for the central comb line.","section":"Sec. III.B, Fig. 3(b)"},{"comment":"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.\"","section":"Sec. IV"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's core result is a useful but incremental extension of the authors' previous work in Ref. [27], and the heavy reliance on this unpublished preprint for the central derivation is a concern for the review process. I would encourage the editor to require the authors to make the manuscript self-contained, either by deriving Eq. (15) or by updating the reference to a published version. The paper is otherwise within the journal's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: this paper convincingly shows that a 300 nm vortex nanopillar driven by RF fields can settle into different metastable gyration orbits under identical drive, with distinct Floquet frequency combs. That is a real extension of the earlier larger-disk work (Ref. [27]), and it is supported by direct micromagnetic simulations, not just a model.\n\nThe fixed-point analysis is the strongest part. The authors extend the Thiele model to include magnon-core scattering, derive a radial equation dR/dt = -Gamma_g R + f(R), and show that the shape of f(R) determines multistability. They extract f(R) from simulated Floquet spectra at two field amplitudes, then predict bifurcation diagrams. The direct simulations with different initial radii agree well with the predicted stable and unstable branches, including the three-state case in Fig. 5. That is good practice: model predictions checked against independent runs.\n\nThe soft spots are real but not fatal. The biggest is the linear-in-brf scaling of f(R), assumed in Sec. III.D and never tested. The arguments leading to Eq. (15) isolate one factor proportional to brf, but the adjacent comb amplitudes in Eq. (16) are themselves drive-dependent, so f(R) could easily have higher-order terms. The bifurcation diagrams in Fig. 4(b,d,f,h) are generated by rescaling f(R) measured at only two brf values, so the quantitative positions of thresholds and coexistence windows could shift. The authors note the extraction uses five sidebands but do not show the analytic formula. The key derivation justifying the scattering terms is also delegated to unpublished Ref. [27]. None of this invalidates the qualitative claim, because Fig. 5 shows three distinct steady orbits from direct simulation, but it does cap the model's quantitative reliability.\n\nA lesser point: no code or data is shipped. Given the paper relies on a semi-empirical f(R) and a fitting procedure, that is a missed opportunity.\n\nWho is this for? People working on vortex dynamics, magnon frequency combs, or unconventional computing schemes based on hysteresis. If the linear-scaling assumption is verified or the model is made fully self-contained, the quantitative predictions will be worth building on. As is, the paper deserves serious refereeing; the core result is sound enough to warrant publication after revision.\n\nRecommendation: send it to peer review. The authors should be asked to test the linear-scaling assumption explicitly or soften the quantitative claims, and to provide the code and data.","headline":"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.","tokens_in":14647,"tokens_out":2173,"would_cite":true,"duration_ms":21425,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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…","keywords":["vortex gyration","Floquet magnons","frequency combs","nonlinear magnon-core scattering","Thiele model","hysteresis","nanopillar magnonics","micromagnetic simulation"],"falsifier":"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.","tokens_in":13533,"feed_emoji":"🧲","tokens_out":4839,"duration_ms":55884,"temperature":0.7,"pith_summary":"This paper tries to establish that the self-induced Floquet mechanism—where GHz spin-wave pumping drives sub-GHz vortex-core gyration—can produce multiple metastable gyration radii under a fixed drive condition. The authors extend the Thiele model to include a nonlinear interaction contribution f(R) arising from three-particle core–magnon scattering, and show that the radial equation dR/dt = -ΓgR + f(R) can have several coexisting stable fixed points. This would mean vortex nanopillars exhibit hysteresis in the gyration radius and in their frequency-comb spectra, making the magnetic response history-dependent. The work is significant because it predicts observable, testable spectral signatures of nonlinear magnon–core coupling in vortex-based nanomagnets.","feed_headline":"One RF drive can lock a vortex core into several orbits","feed_subtitle":"Nonlinear core–magnon scattering creates multiple stable gyration radii and hysteretic frequency combs.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Reports the self-induced Floquet magnon phenomenon and its hysteresis in vortex disks, which this paper extends to nanopillars and the multi-stability analysis.","marker":"[27]"},{"why":"Supplies the Thiele equation for steady-state magnetic domain motion that is the base model extended here.","marker":"[35]"},{"why":"Provides the collective-coordinate expansion used to derive the magnon–core scattering terms in the Lagrangian.","marker":"[37]"},{"why":"Gives the vortex confinement potential and the spin-torque vortex oscillator analogy used to interpret the Hopf bifurcation and auto-oscillation behavior.","marker":"[36]"},{"why":"The mode-resolved supercell method used to compute power spectral densities and extract f(R) from the simulated Floquet spectra.","marker":"[32]"},{"why":"The mumax3 code performs the full time integration of the Landau–Lifshitz–Gilbert equation in the simulations.","marker":"[31]"},{"why":"The magnum.np code computes the linear spin-wave eigenmodes about the static vortex state, fixing the mode indices and frequencies used in the analysis.","marker":"[30]"}],"fun_headline_variants":["One RF drive locks vortex core into multiple orbits","Single drive produces multiple stable vortex gyrations","Vortex core gyration: one field, several stable radii","Hysteretic frequency combs from one RF drive on vortex core","Multiple vortex core orbits from a single RF frequency"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["One RF drive locks vortex core into multiple orbits","Single drive produces multiple stable vortex gyrations","Vortex core gyration: one field, several stable radii","Hysteretic frequency combs from one RF drive on vortex core","Multiple vortex core orbits from a single RF frequency"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000258,"raw_usage":{"total_tokens":1557,"prompt_tokens":892,"completion_tokens":665,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":508,"completion_tokens_details":{"reasoning_tokens":598}},"tokens_in":508,"tokens_out":665,"duration_ms":7391,"temperature":1.0,"reasoning_tokens":598,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T13:56:18.851632+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Thiele equation for steady-state magnetic domain motion that is the base model extended here."},{"cited_title":"Rajaraman, Solitons and Instantons (North Holland, Ams- terdam, 1982)","cited_arxiv_id":null,"evidence_quote":"Provides the collective-coordinate expansion used to derive the magnon–core scattering terms in the Lagrangian."},{"cited_title":"Ivanov and C","cited_arxiv_id":null,"evidence_quote":"Gives the vortex confinement potential and the spin-torque vortex oscillator analogy used to interpret the Hopf bifurcation and auto-oscillation behavior."},{"cited_title":"Massouras, S","cited_arxiv_id":null,"evidence_quote":"The mode-resolved supercell method used to compute power spectral densities and extract f(R) from the simulated Floquet spectra."},{"cited_title":"Vansteenkiste, J","cited_arxiv_id":null,"evidence_quote":"The mumax3 code performs the full time integration of the Landau–Lifshitz–Gilbert equation in the simulations."},{"cited_title":"Bruckner, S","cited_arxiv_id":null,"evidence_quote":"The magnum.np code computes the linear spin-wave eigenmodes about the static vortex state, fixing the mode indices and frequencies used in the analysis."}],"review_version":1}