{"id":"6b000abe-08dd-4ff1-b77d-a35f41d3daf6","arxiv_id":"2507.08147","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Amplitude-modulated optical light can parametrically excite THz antiferromagnetic magnons, generating spin currents, squeezed magnon pairs, and ordered spin patterns.","lead":"The paper proposes using amplitude-modulated light to drive terahertz-frequency spin waves in antiferromagnets without needing terahertz sources. The mechanism could enable spin pumping, squeezed magnon pairs, and ordered spin patterns in future antiferromagnetic spintronic devices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Supplement B's definition of the parametric coupling C is inconsistent with Eqs. (8)/(10); plugging the supplement's C into the threshold gives δJ/J ≈ 3.6%, about 1800× larger than the claimed 2×10^-5, so the central low-power estimate is unsupported as written.","rationale":"The reader's weakest-assumption, that the steady-state amplitudes in Eq. (12) ignore Gilbert damping, is a valid quantitative concern: at δJ near threshold the true damped amplitude is parametrically smaller than the undamped sqrt(δJ) result, so the spin-pumping and ISH estimates in Sec. VI need revision. However, I find a more load-bearing internal inconsistency in the definition of the parametric coupling C, which directly controls the threshold (Eq. (10)) and hence the central practical claim of low-power THz magnon excitation. Supplement B defines C via derivatives of the eigenmode normalization factors D1, D4; for the small-anisotropy regime quoted in the paper, this C is of order Dx/(8J^2) and would push the threshold to δJ/J ≈ 3.6%, more than three orders of magnitude above the quoted value. A direct derivation from Eq. (A5) instead gives an effective C ≈ 1/2, which would restore the order of magnitude of the quoted threshold. The manuscript as written does not reconcile these two definitions; the main text labels C dimensionless while the supplement's expression has units of inverse energy. Because the threshold is the foundation of the proposed experimental protocol, the paper cannot be accepted without correcting this inconsistency and recomputing the threshold and the derived estimates. The qualitative mechanism—parametric excitation via modulated exchange—remains plausible, so the appropriate action is to keep the reader's CONDITIONAL verdict rather than rejecting outright; the flaws are fixable but require substantive revision.","tokens_in":14243,"tokens_out":26812,"duration_ms":270400,"concrete_test":"Take k = 0 in the linearized LLG, Eq. (A5), substitute J(t) = J̄ + δJ cos(2ωdt), and solve the resulting 4×4 system by Floquet analysis to extract the exponential growth rate for the ω1 resonance. Compare this rate with the main text's Eq. (9) to determine the effective C; equivalently, project Eq. (A5) explicitly onto the instantaneous eigenmodes, keeping both the δM matrix element and the time dependence of D1/D4, and check whether the resulting C is ≈ 1/2 (as the direct Mathieu reduction suggests) or ≈ Dx/(8J) ≈ 2.5×10^-4 (as Supplement B's formula gives).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that MFPD can excite THz AF magnons at modest optical powers rests on the threshold in Eq. (10), which depends on the coupling coefficient C. Supplement B defines C in Eq. (A7) as C = ∂logD1/∂J − ∂logD4/∂J. For k = 0 with Dz,Dx ≪ J, D1 ≈ D4 ≈ 2, so C ≈ Dx/(8J^2); with the paper's parameters (J = 50 meV, Dx = 0.1 meV), the dimensionless combination J·C ≈ 2.5×10^-4. Inserting this into Eq. (10) gives δJ/J ≈ 8αωd/(DxS) ≈ 3.6×10^-2 for α = 10^-4, ωd ≈ 4.5 meV, S = 1, roughly 1800× larger than the δJ/J ≈ 2×10^-5 quoted in Sec. VI. In contrast, direct linearization of Eq. (A5) at k = 0 yields the reduced equations u˙ = Dxv and v˙ = −(4J + Dz)u − 4δJ cos(2ωdt)u, a Mathieu-type system with effective dimensionless coupling C ≈ 2J/(4J + Dz) ≈ 1/2. Thus the supplement's C is not the coefficient that appears in the main text's Eqs. (8) and (10); the main text calls C dimensionless while the supplement's definition has units of 1/J. This internal inconsistency makes the threshold derivation unreproducible and undermines the quantitative power-reduction argument.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes Modulated Floquet Parametric Driving (MFPD) as a way to excite THz antiferromagnetic magnons using amplitude-modulated optical light. The modulation of the exchange coupling J(t) = J̄ + δJ cos(2ω_d t) is argued to parametrically drive a magnon mode when the modulation frequency matches the magnon frequency, with a damping-controlled threshold given by Eq. (10). Above threshold, nonlinearities saturate the instability into steady states, which the paper uses to predict spin pumping into an adjacent metal, dynamical pattern formation in one and two dimensions, and the generation of two-mode squeezed magnon states. The analytic results are supplemented by classical LLG simulations with Dedalus and tensor-network simulations with TenPy.","tokens_in":14649,"tokens_out":5496,"duration_ms":64311,"significance":"If correct, the mechanism is significant: it offers a route to drive THz magnons without THz sources, with potential consequences for antiferromagnetic spintronics and quantum magnonics. The paper includes a clean linear parametric instability analysis, a quantum two-mode squeezing derivation, and numerical simulations of pattern formation. However, the central quantitative claims—the low-power threshold and the spin-pumping voltage estimate—currently rest on an inconsistent definition of the coupling coefficient C and on steady-state amplitudes derived without damping. These issues are load-bearing and need to be resolved before the proposal's practical viability can be assessed.","major_comments":[{"comment":"The definition of the coupling coefficient C is inconsistent with the threshold formula and with the text. In the main text, immediately after Eq. (8), C is called dimensionless, and Eq. (10) uses it in δJ/J > (1/(JSC))·4αω_dℏ/(4J + D_x + 2D_z). In Supplement B, however, C is defined as C = ∂log D_1/∂J − ∂log D_4/∂J, which has dimensions of inverse energy. For k = 0 and D_z, D_x ≪ J, D_1 ≈ D_4 ≈ 2, so J·C ≈ D_x/(8J) ≈ 2.5×10^{-4} for J = 50 meV, D_x = 0.1 meV. Inserting this into Eq. (10) yields δJ/J ≈ 3.6×10^{-2}, about 1800 times larger than the quoted δJ/J ≈ 2×10^{-5}. If instead C is treated as the dimensionless O(1) coefficient suggested by a direct Mathieu-type linearization of Eq. (A5), then Eq. (10) is dimensionally inconsistent as written. The supplement's C therefore does not appear to be the coefficient governing the coupling in Eqs. (8) and (9). This must be reconciled, and the threshold estimate and the field-reduction factor in Sec. VI must be recomputed with the correct coefficient.","section":"Supplement B and Eq. (10)"},{"comment":"The steady-state amplitudes in Eq. (12) are derived from the amplitude equations (A11) without including Gilbert damping. As a result they scale as √δJ and vanish only as δJ → 0, not at the finite threshold δJ_th of Eq. (10). The spin-pumping formulas (13)-(15) and the inverse spin Hall voltage estimate in Sec. VI are then evaluated at 'δJ of the order of the threshold'. At threshold, the saturated amplitude in a damped parametrically driven system should instead vanish, scaling roughly as √(δJ − δJ_th) above threshold. The quantitative spin-pumping estimate therefore needs revision: either the amplitude equations must be solved with damping included, or the estimate must be restricted to δJ sufficiently above threshold with the appropriate suppression factor.","section":"Supplement C and Sec. VI, Eqs. (12)-(15)"},{"comment":"The dimensions/prefactors in the quantum treatment should be checked carefully. The Hamiltonian in Eq. (16) has a pair-generation term with coefficient δJ·C/2, while the interaction-picture Hamiltonian is later written as H_int = δJ/(4J)·C(α†_{k*}α†_{−k*} + α_{k*}α_{−k*}), and the squeezing parameter is r = tδJ C/(4ℏ). Depending on whether C is dimensionless or has units of 1/J, and on how δJ is normalized relative to J, these expressions change by factors of J and ω_d. The relation between the classical coupling in Eq. (8), which contains ω_1 ≈ ω_d and δJ/J, and the quantum coupling in Eq. (16) should be stated explicitly so that the squeezing parameter is unambiguous.","section":"Sec. V, Eq. (16) and the squeezing parameter r"}],"minor_comments":[{"comment":"The caption says 'according to Eq. (2)' but the dispersions are given by Eq. (5); please correct the reference.","section":"Figure 2 caption"},{"comment":"The simulation parameters list 'D_z = 0.1, D_z = −0.2' twice. The second should presumably be D_x, since the model requires D_x > 0; please fix this typo and confirm the sign convention used in the simulation.","section":"Supplement D"},{"comment":"The sentence 'Since δJ ∼ √E' appears to be a typo: from Eq. (2), δJ ∝ E^2, so a threshold δJ/J = 2×10^{-5} corresponds to a field-amplitude reduction factor of about 1/√(2×10^{-5}) ≈ 200. Please clarify the proportionality.","section":"Sec. VI, Discussion"},{"comment":"The symbols in Eq. (A11) mix δJ, ϵ, and δJ̃; the sentence 'we wrote ϵδJ̃ = δJ' should define the dimensionless and dimensional quantities more explicitly to avoid confusion about the expansion parameter.","section":"Supplement C, Eq. (A11)"}],"recommendation":"major_revision","confidential_remarks":"The paper's central idea is interesting and worth pursuing, but the current manuscript has a load-bearing internal inconsistency in the definition and numerical value of C that makes the headline threshold and power-reduction estimates unreproducible, and the steady-state amplitudes used for spin pumping are derived without damping. These are fixable with a careful revision, but the quantitative conclusions should not be taken at face value until they are corrected. I do not see grounds for rejection, provided the authors can reconcile the linearized coupling with the threshold formula and re-derive the saturated amplitudes with damping."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a serious paper, worth refereeing, but it needs a round of revision on two specific points. The MFPD idea is theirs from the plasmon papers, so the new thing here is the application to AF magnons: threshold, steady states, spin pumping, two-mode squeezing, and stripe pattern formation. The linear instability and the squeezing derivation are clean; the pattern-formation numerics (Dedalus) and the TenPy quantum simulation are concrete evidence.\n\nI checked the stress-test note about C. It doesn't hold up. The supplement defines C = ∂logD1/∂J − ∂logD4/∂J; at k=0 with Dx,Dz ≪ J, D1≈2 and D4≈sqrt((Dx+Dz)/J), so C≈1/(2J), not Dx/(8J^2). Plugging C≈1/(2J) into Eq. (10) gives δJ/J roughly 2–3×10^-5 for their parameters, matching the paper's estimate. So the central low-power claim is not off by 1800x.\n\nThat said, there is a real units inconsistency: the supplement's C has units 1/energy, while the main text Eq. (8) uses it as dimensionless. The threshold formula Eq. (10) only makes dimensional sense with C in 1/energy. So the paper is internally consistent in one place and wrong in another; the authors need to clarify what C is and fix Eq. (8) accordingly. This is a typo-level fix, not a load-bearing flaw.\n\nThe bigger soft spot is the steady state. Supplement C derives Eq. (12) from the LLG with no Gilbert damping. So the amplitudes scale as sqrt(δJ) and are finite for any δJ>0, including below the threshold. In the damped system the saturated amplitude should vanish at threshold and grow as sqrt(δJ−δJ_th). The paper then uses Eq. (12) at 'δJ of the order of the threshold' for the spin-pumping and inverse-spin-Hall estimates, which will overestimate the effect near threshold. That needs a quantitative redo. The mechanism survives, but the numbers in Sec. VI are not trustworthy as written.\n\nThere are also minor typos (e.g., 'eigenvelue', Dz=−0.2 in Supplement D) that should be caught.\n\nBottom line: the qualitative mechanism is sound, the derivations are mostly clean, and the paper opens a practical route to THz magnons without THz sources. With the damping issue fixed and C clarified, it would be a solid contribution. A serious referee should engage with it; I'd send it out with a request for revision.","headline":"Good physics with a units slip in C and an undamped steady-state; the stress-test's 1800x complaint doesn't survive contact with the supplement.","tokens_in":15143,"tokens_out":8755,"would_cite":true,"duration_ms":82742,"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":"Amplitude-modulated light can parametrically drive terahertz magnons in antiferromagnets, through a mechanism the authors call Modulated Floquet Parametric Driving.","keywords":["antiferromagnetic magnons","terahertz spin dynamics","modulated Floquet parametric driving","parametric instability","spin pumping","two-mode squeezing","magnon pattern formation","inverse spin Hall effect"],"falsifier":"Drive a high-quality antiferromagnet capped with a normal metal using an amplitude-modulated optical beam, and record the dc inverse spin Hall voltage while sweeping modulation frequency and power. The paper predicts a sharp onset at the damping-controlled threshold, a microvolt-scale voltage at resonance for realistic hematite/Pt parameters, and no signal for an unmodulated beam of the same intensity; observing a signal that grows smoothly from arbitrarily small drive would falsify the quantitative picture.","tokens_in":14078,"feed_emoji":"🧲","tokens_out":14092,"duration_ms":142832,"temperature":0.7,"pith_summary":"The paper proposes that a coherent optical beam with a modulated amplitude can resonantly excite antiferromagnetic magnons in the terahertz range without requiring a terahertz source. The mechanism works because the light makes the exchange coupling oscillate, $J(t)=\\bar{J}+\\delta J\\cos(2\\omega_d t)$; the $2\\omega_d$ component acts as a parametric pump on magnon modes with $\\omega_1(k^*)=\\omega_d$, and once the modulation depth exceeds a damping-controlled threshold the magnon amplitude grows until nonlinearities saturate it. In the saturated state the authors derive concrete experimental signatures: dc spin pumping into a neighboring normal metal, two-mode squeezed and entangled magnon pairs, and, at finite wavevectors, standing-wave or stripe spin patterns that break the symmetries of the lattice. This matters because antiferromagnetic resonances are naturally in the THz range but are hard to address directly, and modulated optical frequencies up to 10 THz are already achievable in the lab.","feed_headline":"Optical beat notes can parametrically drive THz magnons","feed_subtitle":"Amplitude-modulated lasers can replace THz sources, enabling spin pumping, squeezing, and spin-pattern formation.","key_machinery":"The central object is the modulated Floquet parametric drive (MFPD): a high-frequency optical carrier whose amplitude is modulated at a lower frequency, so the exchange constant acquires the oscillating term $\\delta J\\cos(2\\omega_d t)$. The argument is carried by the linearized Landau-Lifshitz-Gilbert dynamics, where this term couples the magnon eigenmodes $\\delta S_{\\mathrm{eig},1}$ and $\\delta S_{\\mathrm{eig},2}$; a slowly varying envelope approximation near $\\omega_1(k^*)=\\omega_d$ gives exponential growth with rate proportional to $\\delta J C\\omega_d/(2J)$, and including Gilbert damping turns that growth into the threshold condition of Eq. (10). On the quantum side the same coupling becomes the two-mode squeezing Hamiltonian $\\frac{\\delta J C}{2}\\cos(2\\omega_d t)(\\alpha_k^\\dagger \\alpha_{-k}^\\dagger+\\alpha_k\\alpha_{-k})$, which directly produces the entangled magnon pairs of Eq. (17).","core_discovery":"As the authors state it, the central claim is that the time-dependent exchange coupling created by amplitude-modulated light acts as a parametric pump for antiferromagnetic magnons: $J(t)=\\bar{J}+\\delta J\\cos(2\\omega_d t)$, with the resonance condition $\\omega_1(k^*)=\\omega_d$ for the magnon mode that is pumped. Above the damping-controlled threshold of Eq. (10) the Néel state becomes linearly unstable, and the nonlinear terms in the Landau-Lifshitz-Gilbert equation saturate the growth into a steady state with the magnon amplitudes of Eq. (12). From that steady state the paper obtains three observable consequences: a dc spin current $I_{s,z}=2\\omega_d G(a\\,\\delta b_{H,A}+b\\,\\delta a_{H,A})$ pumped into an adjacent normal metal (Eq. (15)), a two-mode squeezed magnon state with squeezing parameter $r=t\\,\\delta J C/(4\\hbar)$ (Eq. (17)), and, for $k^*\\neq 0$, symmetry-breaking standing-wave and stripe spin patterns in one and two dimensions. All of these follow from the same parametric instability, without the need for direct THz driving.","pith_inferences":["The authors do not spell this out, but the sharp threshold in Eq. (10) gives an experimental discriminant between MFPD and laser heating: an unmodulated beam of identical average intensity should produce no dc spin-pumping signal, while a modulated beam should show an abrupt onset.","Since the resonant wavevector $k^*$ is set by where $\\omega_1(k)=\\omega_d$, sweeping the modulation frequency should continuously tune the stripe wavelength in two-dimensional systems, turning the predicted pattern formation into a frequency-controlled texture.","The same parametric-pair Hamiltonian should apply to any bosonic collective mode whose coupling responds quadratically to light, such as optical phonons or exciton-polaritons; the squeezing and pattern-formation predictions may therefore transfer beyond magnons."],"forward_implications":["Terahertz magnon resonances can be excited with amplitude-modulated optical light, bypassing the need for dedicated THz sources.","A driven antiferromagnet in contact with a normal metal should inject a dc spin current, measurable as a microvolt-scale inverse spin Hall voltage for realistic hematite/Pt parameters.","The parametrically driven magnon state is a two-mode squeezed vacuum in which magnons of opposite momenta are entangled, offering a route to squeezed and entangled states for quantum magnonics.","Driving at finite wavevector selects standing-wave and stripe spin textures that break translational and rotational symmetry, so the drive acts as a switch into a symmetry-broken dynamical phase.","Because the high quality factor of the magnon mode lets it accumulate drive energy over many cycles, the required optical powers are roughly two orders of magnitude lower than in direct exchange-modulation schemes."],"supporting_citations":[{"why":"Supplies the microscopic basis for light-induced modification of the exchange constant, the starting point for the parametric drive.","marker":"[35–38]"},{"why":"Shows that amplitude-modulated laser signals with modulation frequencies up to 10 THz can be produced by interfering detuned sources, making the proposed drive experimentally feasible.","marker":"[33, 34]"},{"why":"Reports the ultra-low damping ($\\alpha \\le 10^{-5}$) in hematite that puts the threshold drive into an accessible power range.","marker":"[10, 11]"},{"why":"Established Modulated Floquet Parametric Driving for plasmonic modes, the mechanism this paper extends to antiferromagnetic magnons.","marker":"[31, 32]"},{"why":"Provides the interfacial spin-pumping formula used to predict the dc spin current and inverse spin Hall voltage.","marker":"[43, 44]"},{"why":"Specifies the spectral solver used for the 1D and 2D Landau-Lifshitz-Gilbert simulations that demonstrate the predicted spin patterns.","marker":"[62]"},{"why":"Supports the expectation that parametrically driven nonlinear systems select discrete wavevector sets, the basis for the predicted pattern formation.","marker":"[59–61]"}],"fun_headline_variants":["Modulated light parametrically drives THz magnons","Amplitude-modulated light pumps THz magnons","Optical modulation excites THz magnon pairs","THz magnons from parametric light drive","Light beat notes parametrically drive THz magnons"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the saturated magnon amplitudes used to predict spin pumping are computed without Gilbert damping, so the predicted signal size near threshold ignores the very damping that creates the threshold.","fun_headline_variants_meta":{"raw":{"variants":["Modulated light parametrically drives THz magnons","Amplitude-modulated light pumps THz magnons","Optical modulation excites THz magnon pairs","THz magnons from parametric light drive","Light beat notes parametrically drive THz magnons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000748,"raw_usage":{"total_tokens":3315,"prompt_tokens":908,"completion_tokens":2407,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":524,"completion_tokens_details":{"reasoning_tokens":2332}},"tokens_in":524,"tokens_out":2407,"duration_ms":18200,"temperature":1.0,"reasoning_tokens":2332,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T18:28:10.431020+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Drive a high-quality antiferromagnet capped with a normal metal using an amplitude-modulated optical beam, and record the dc inverse spin Hall voltage while sweeping modulation frequency and power. The paper predicts a sharp onset at the damping-controlled threshold, a microvolt-scale voltage at resonance for realistic hematite/Pt parameters, and no signal for an unmodulated beam of the same intensity; observing a signal that grows smoothly from arbitrarily small drive would falsify the quantitative picture.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Specifies the spectral solver used for the 1D and 2D Landau-Lifshitz-Gilbert simulations that demonstrate the predicted spin patterns."}],"review_version":1}