{"id":"2eb3c081-7416-403e-a0a5-211b2fd9f583","arxiv_id":"2607.18073","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Higher-order parametric resonances let low-frequency pumping excite magnons in YIG, with threshold amplitudes growing as the n-th root of the damping.","lead":"This paper uses Floquet theory to show that magnons in a thin magnetic film (YIG) can be excited by microwave pumping at frequencies well below the magnon energy, via higher-order parametric resonances. It derives analytic thresholds for the first three resonances and confirms them with micromagnetic simulations, pointing toward low-frequency, wave-number-selective magnon generation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eqs. (17) and (21) for the n=2 and n=3 thresholds are dimensionally inconsistent as printed, undermining the analytic derivation of h_th,min ∝ α^{1/n}.","rationale":"The reader's verdict of CONDITIONAL is reasonable, but the weakest assumption identified there — linearization validity near resonance — is not the most load-bearing issue. The printed n=2 and n=3 perturbation formulas (Eqs. 17 and 21) are dimensionally inconsistent, which prevents an independent reader from verifying the derivation of the headline α^{1/n} scaling. The final threshold expressions (Eqs. 19 and 22) are dimensionally homogeneous, so the scaling exponent may survive after correcting the displayed coefficients; nevertheless, the central analytic support needs correction before the paper can be accepted. The numerical Floquet results and MuMax3 simulation provide supporting evidence, but they do not resolve the internal inconsistency of the printed formulas. I therefore recommend CONDITIONAL rather than ACCEPT, and agree with the reader's overall verdict direction, though for a different, more specific reason.","tokens_in":12572,"tokens_out":18763,"duration_ms":177098,"concrete_test":"Independently rederive the n=2 and n=3 quasi-energies from Eq. (13) by symbolic second- and third-order perturbation theory, keeping all factors of ω_k explicitly. Verify whether c0 = -B_k^2/(3ω_k^3) (not /ω_k^2) and Δϵ^(2) ∝ h1^2B_k/ω_k^2 (not /ω_k^3). Then recompute the n=2 threshold at k=2.5 rad/µm, α=0.0002, with the corrected coefficients and compare to the numerical Floquet result in Fig. 2; agreement to within a few percent would show the inconsistency is typographical and the scaling survives.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central scaling h_th,min ∝ α^{1/n} (Eq. 24) rests on degenerate Floquet perturbation theory, but the n=2 and n=3 expressions are not derived and, as written, contain unit errors. In Eq. (17), with A_k, B_k, ω_k, and h1 all having energy/frequency units, c0 = -B_k^2/(3ω_k^2) is dimensionless while c1 = Ā_k^2B_k^2/(4ω_k^6) has units 1/ω^2. The combination Δω2 - c0 h1^2 then adds an energy to an energy^2, so the square-root argument is dimensionally inconsistent. In Eq. (21), Δϵ^(2) ≈ -i h1^2 9B_k/(32ω_k^3) is dimensionless (h1^2B/ω^3), not a quasi-energy correction. If these are merely typographical missing powers of ω, the final α^{1/n} scaling may survive; if not, the claimed thresholds are unsupported. The paper's own Appendix B states the linearized description breaks down at resonance, so the micromagnetic confirmation is only qualitative and cannot resolve this inconsistency.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies parametric magnon excitations in thin-film YIG driven by a microwave field, with emphasis on driving frequencies below the magnon energy. Starting from a spin Hamiltonian with exchange and dipole interactions, the authors use the Holstein-Primakoff transformation, the linear spin-wave approximation, and the uniform-mode approximation to obtain a time-periodic quadratic magnon Hamiltonian. They then set up the linearized Landau-Lifshitz-Gilbert equation and analyze its Floquet quasi-energies. Degeneracies at nωD = 2ωk are identified as parametric resonances, and degenerate Floquet perturbation theory is used to derive threshold amplitudes for n = 1,2,3, culminating in the prediction h_th,min ∝ α^{1/n}. A gauge-transformation argument for strong driving predicts that the instability regions are modulated by Bessel functions J_n(2h1/ωD) and vanish at the zeros of these functions. The analytic instability thresholds are compared with the numerical Floquet solution of the same linearized model, and the n = 3 case is additionally compared with MuMax3 micromagnetic simulations, which show Floquet replicas and resonant magnon growth at drive amplitudes above threshold.","tokens_in":12879,"tokens_out":12802,"duration_ms":179017,"significance":"The paper addresses a topical and technically useful problem: low-frequency parametric pumping of magnons below the magnon band edge. If the analytic predictions are correct, the clean power law h_th,min ∝ α^{1/n} and the Bessel-function structure at strong drive are falsifiable predictions that could guide targeted magnon excitation experiments. The authors are to be credited for deriving the thresholds rather than fitting them, for checking the n = 1 and n = 2/3 thresholds against the exact numerical Floquet solution of the same linearized equation, and for including an independent micromagnetic simulation for the n = 3 resonance. The n = 1 expression is dimensionally sound and is in good agreement with numerics. However, the n = 2 and n = 3 analytic results are stated through final coefficients without derivation, and as printed they contain dimensional inconsistencies; the central α^{1/n} scaling therefore is not yet fully supported. The manuscript has the potential to be a valuable contribution, but the load-bearing n = 2 and n = 3 derivations need to be corrected and supplied.","major_comments":[{"comment":"The n = 2 quasi-energy formula is dimensionally inconsistent as written. Since h1, Bk, ωk, Ak, and Δω2 all have energy/frequency units, c0 = −Bk^2/(3ωk^2) is dimensionless, so c0 h1^2 has units of frequency^2 and cannot be subtracted from Δω2 inside the first square in Eq. (17). Likewise, c1 in Eq. (18) has units of 1/ω^2 but c0^2 is dimensionless, so the denominator c1 − c0^2 mixes incompatible units. If this is a typographical omission of a power of ω in c0, it must be fixed explicitly; as printed, Eq. (18) and the n = 2 threshold hth,min ∝ α^{1/2} are not reproducible.","section":"Theory, Eqs. (17)–(18)"},{"comment":"The n = 3 result has the same class of problem. The second-order quasi-energy correction in Eq. (21), Δϵ^(2) ≈ −i h1^2 9Bk/(32ωk^3), is dimensionless (h1^2 Bk/ω^3) rather than an energy/frequency, which is what Eq. (20) requires. Moreover, Δϵ^(2) and c̃1 are only stated as final expressions; no derivation is provided. Since Eq. (22) and the α^{1/3} law rest on these expressions, the central multi-order scaling claim is unsupported until either the expressions are corrected (e.g., an omitted power of ω) and the second/third-order perturbation calculation is included, or the final formulas are otherwise verified.","section":"Theory, Eqs. (20)–(22)"},{"comment":"The transition to the strong-driving treatment is internally inconsistent as written. The paragraph begins with 'We now consider larger driving h1' and then states 'Progress can be made for ωD/h1 ≫ 1'. That inequality is the opposite of the large-h1 regime and is violated for the large amplitudes shown in Fig. 4. If the actual small parameter is Bk J_n(2h1/ωD)/ωD, this should be stated and justified; otherwise the Bessel-function prediction for the disappearance of resonances at high drive, Eqs. (28)–(30), is not properly grounded.","section":"Strong-drive analysis, before Eq. (25)"},{"comment":"The authors explicitly acknowledge in Appendix B that the linearized description 'quickly breaks down' at resonance. This is acceptable for linear-instability threshold calculations, but it should be made clearer in the main text that the MuMax3 comparison is only a qualitative confirmation: the simulation thresholds are higher than the analytic ones, and the analytic formulas therefore should not be presented as quantitatively predictive for the full nonlinear system without stating this limitation. The sentence in the simulations section already hints at this, but it deserves to be part of the abstract-level claims.","section":"Appendix B / Simulation section"}],"minor_comments":[{"comment":"Typo: 'exciatations' should be 'excitations'.","section":"Abstract"},{"comment":"'complimented by micromagnetic simulations' should be 'complemented by micromagnetic simulations'.","section":"Introduction and Abstract"},{"comment":"The caption refers to 'the analytical predictions for the thresholds in Eqs. (16), (18) and (23)'; Eq. (23) is a width formula, not a threshold expression. The threshold for n = 3 is given in Eq. (22). This should be corrected.","section":"Fig. 2 caption"},{"comment":"The caption says the dashed curve is the 'numerical solution of Eq. (8)', but Eq. (8) is a differential equation; presumably the numerical Floquet solution of Eq. (13) is meant. Please clarify.","section":"Fig. 1 caption"},{"comment":"Typo: 'and and film thickness' should be 'and film thickness'.","section":"Fig. 4 caption"}],"recommendation":"major_revision","confidential_remarks":"The dimensional inconsistencies in Eqs. (17) and (21) look like fixable typos rather than a fundamental flaw in the method, because the n = 1 result is sound and the MuMax3 simulation gives qualitative support for n = 3. However, the paper as submitted states the n = 2 and n = 3 perturbative coefficients without derivation, so a referee cannot verify the central α^{1/n} claim. If the authors can provide a corrected and derived version of these formulas, the paper would likely be acceptable. The strong-driving section also needs a corrected small-parameter statement. I do not see grounds for rejection if these points are properly addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper has a clear central idea: use Floquet degenerate perturbation theory to get threshold amplitudes for the first three parametric magnon resonances in thin-film YIG, including the α^{1/n} scaling and a Bessel-function collapse at strong drive. The n=1 resonance is derived cleanly, matches the numerical solution of the linearized LLG, and the n=3 resonance is qualitatively reproduced in MuMax3. That part is real and useful for magnonics.\n\nBut there is a problem in the analytic core. In Eq. (17), the term c0 h1^2 has units of energy squared while Δω_2 has units of energy, so the square root mixes incompatible quantities. In Eq. (21), Δϵ^(2) is dimensionless, not a quasi-energy correction. These are not cosmetic typos in the final formulas; they sit inside the expressions that produce the n=2 and n=3 thresholds. Since the derivations are not shown beyond stating coefficients, a referee cannot tell whether this is a missing power of ω or a deeper mistake. The n=1 formula is dimensionally sound, so I would bet the underlying scaling survives, but the paper as printed does not support its central claim for n≥2.\n\nThe strong-drive section is interesting, but the predicted Bessel-function zeros are explicitly said to be unobservable as an overall reduction in realistic simulations, and no simulation of that regime is attempted. The MuMax3 comparison is good for n=3, though the paper itself notes in Appendix B that the linear theory breaks down near resonance, so the comparison is only qualitative.\n\nWho is this for: experimentalists and theorists working on parametric pumping of magnons at sub-band frequencies. The topic is timely and the n=1 framework is worth knowing. But the paper needs a proper derivation of the second- and third-order perturbation theory with correct units before the α^{1/n} claim is reliable.\n\nMy recommendation: send it to peer review, not desk reject, but the referees should demand the full calculation and the corrected equations. If the authors fix the units and show the steps, it becomes a credible and citable paper.","headline":"Solid n=1 result but dimensional inconsistencies in the n=2/n=3 threshold formulas undermine the central scaling claim as printed.","tokens_in":13344,"tokens_out":4783,"would_cite":false,"duration_ms":46530,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Parametric resonances at nω_D=2ω_k mean a drive below the magnon energy can still excite magnons, with thresholds that soften as α^{1/n}.","keywords":["parametric resonance","magnons","Floquet theory","spin waves","Landau-Lifshitz-Gilbert","YIG thin films","threshold amplitude","parametric pumping"],"falsifier":"A room-temperature YIG film experiment measuring the n=3 threshold amplitude for films with known damping α, from roughly 0.0002 to 0.05, would test the law h_th∝α^{1/3}; if the measured exponent departs clearly from 1/n, or if no instability appears at ω_D=2ω_k/3 below threshold, the central claim is falsified. Observing that resonance positions shift with damping would also contradict the predicted robustness.","tokens_in":1340,"feed_emoji":"🧲","tokens_out":1303,"duration_ms":53304,"temperature":0.7,"pith_summary":"This paper asks whether magnons in a thin magnetic film can be excited by a microwave drive whose frequency is lower than the magnon's own energy. The authors show that it can: parametric resonances appear whenever n times the drive frequency equals twice the magnon frequency, for n=1,2,3 and beyond. Using Floquet degenerate perturbation theory on the linearized Landau-Lifshitz-Gilbert equation, they derive analytic formulas for the instability regions and for the threshold drive amplitude, which follows a power law in the Gilbert damping that softens as n increases. They confirm the predictions numerically and with micromagnetic simulations at the n=3 resonance, mapping the instability tongues and the off-resonant Floquet replicas. The practical payoff would be the ability to create magnons at a chosen wavenumber with a drive frequency below the magnon band, where ordinary linear excitation is suppressed.","feed_headline":"Drive magnets below magnon energy and they still resonate","feed_subtitle":"Floquet theory yields threshold law h ∝ α^{1/n}, enabling targeted wavenumber excitation with a low-frequency drive.","key_machinery":"The central object is Floquet degenerate perturbation theory applied to the linearized Landau-Lifshitz-Gilbert equation. The driven spin system is expanded in Floquet replicas |σ,m⟩, and degeneracies between replicas at nω_D=2ω_k define the resonances. The resulting two-by-two effective quasi-energy matrix yields instability when the imaginary part of the quasi-energy is positive; exceptional points, where real and imaginary parts both degenerate, separate stable from unstable regions. For strong driving, a gauge transformation re-expresses the coupling through Bessel functions J_n(2h_1/ω_D), predicting that resonance widths oscillate and collapse at Bessel zeros.","core_discovery":"The central claim is that the parametric resonance condition nω_D=2ω_k holds for all n in a driven thin-film ferromagnet, even when the drive frequency ω_D lies below the magnon energy ω_k. The novel quantitative results are the closed-form threshold amplitudes: h_th,min=2αω_k Ā_k/B_k for n=1, roughly ω_k√(2αω_k/B_k) for n=2, and ω_k(2αĀ_k/(B_k√c_k))^{1/3} for n=3, together with the general scaling h_th,min∝α^{1/n}. These thresholds define the boundary of the instability region in drive-amplitude–frequency space, where the imaginary part of the Floquet quasi-energy turns positive. The paper also shows that damping changes the size of the instability regions but not their positions, and that","pith_inferences":["An extension the authors do not spell out: if the α^{1/n} law holds in experiment, a single low-frequency microwave source could address different wavenumbers by tuning amplitude, since higher-order thresholds rise more slowly with damping than the n=1 threshold does.","The linear-theory thresholds are lower than the micromagnetic thresholds, suggesting that quantitative predictions for real YIG films will need a renormalized effective damping or nonlinear corrections; the qualitative resonance map, however, should persist.","The Bessel-function collapse of resonances at high drive amplitudes implies a form of dynamical decoupling at selected wavenumbers, which could be used to selectively suppress excitation of one wavenumber while still pumping others.","The Floquet treatment is transferable to other thin-film magnets and to acoustic or strain-driven parametric pumping, where the effective coupling takes the same mathematical form."],"forward_implications":["At a drive frequency ω_D=2ω_k/3, the n=3 resonance permits magnon creation with an onset around 200 Oe for YIG, so excitations can be created below the magnon energy.","Threshold amplitudes for higher-order resonances scale as α^{1/n}, so low-damping materials make high-order resonances accessible at modest field amplitudes.","Damping shifts the size of the instability regions but not their positions, so the resonance condition nω_D=2ω_k is robust against dissipation.","At large driving amplitudes the resonance width is controlled by Bessel functions and can collapse at special amplitude-frequency combinations, while the dispersion itself is renormalized by the drive.","Off-resonant Floquet replicas appear at intervals of the driving frequency in the magnon spectrum, providing a direct signature of the drive."],"fun_headline_variants":["Low-freq pump excites magnons via parametric resonance","Sub-energy pumping drives magnons into resonance","Magnons resonate under low-frequency drive","Below-energy drive still resonates magnons","Parametric resonance: low-frequency pump excites magnons"],"cache_read_input_tokens":14720,"weakest_assumption_plain":"The threshold formulas assume the linear spin-wave approximation and a constant Gilbert damping hold up to the moment the instability starts; the paper itself notes that at resonance the linear description quickly breaks down because the magnon number grows exponentially.","fun_headline_variants_meta":{"raw":{"variants":["Low-freq pump excites magnons via parametric resonance","Sub-energy pumping drives magnons into resonance","Magnons resonate under low-frequency drive","Below-energy drive still resonates magnons","Parametric resonance: low-frequency pump excites magnons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000271,"raw_usage":{"total_tokens":1439,"prompt_tokens":694,"completion_tokens":745,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":438,"completion_tokens_details":{"reasoning_tokens":675}},"tokens_in":438,"tokens_out":745,"duration_ms":8490,"temperature":1.0,"reasoning_tokens":675,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T16:11:15.170802+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A room-temperature YIG film experiment measuring the n=3 threshold amplitude for films with known damping α, from roughly 0.0002 to 0.05, would test the law h_th∝α^{1/3}; if the measured exponent departs clearly from 1/n, or if no instability appears at ω_D=2ω_k/3 below threshold, the central claim is falsified. Observing that resonance positions shift with damping would also contradict the predicted robustness.","supporting_citations":[],"review_version":1}