{"id":"a7023915-87bf-4c28-9f46-3642865cf0c4","arxiv_id":"2502.00849","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Simulations of a two-orbital Hubbard model show that NiO's Néel vector changes length rather than rotating under ultrafast laser drive, and that proximity spin-orbit coupling, not an interfacial spin current, explains the observed terahertz emission.","lead":"A numerical study of a small model for the antiferromagnet NiO finds that ultrafast laser pulses shrink the Néel vector and the magnetization without rotating them, and that terahertz emission appears only when spin-orbit coupling from an adjacent heavy metal is included. The work offers a microscopic alternative to the usual spin-current picture of terahertz emission from NiO/Pt bilayers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The magnitude of the Rashba SOC tSO is never specified, making the central 'SOC turns on THz emission' claim untestable; a tSO scan is needed.","rationale":"The reader identified the simplified model (no O, no Pt, unspecified tSO, impurity field) as the weakest assumption. Among these, the unspecified tSO is the single most load-bearing for the THz mechanism claim because it is the only parameter that distinguishes the 'plain' and 'SOC' cases that the paper compares to experiments. Without its value, the reader cannot judge whether the computation used a realistic proximity SOC or a contrived strength. The missing O and Pt layers are more model-level approximations whose impacts are partially mitigated by the effective parameters (U, JH, J) and by the monolayer geometry; they do not directly control the on/off comparison. The impurity field B_imp is a concern but it is at least specified (0.1 eV) and serves as a symmetry-breaking proxy for magnetic anisotropy; tSO is entirely suppressed. The proposed tSO scan would turn the qualitative 'significant only with SOC' into a quantitative, testable statement. If the scan passes, the central claim is substantially strengthened; if it fails, the claim is refuted. Thus the appropriate verdict remains CONDITIONAL: accept only after the tSO value and scan are provided.","tokens_in":24,"tokens_out":9290,"duration_ms":165937,"concrete_test":"Rerun the subgap-pump ED calculation of Fig. 3 for tSO = 0.02, 0.05, 0.1, 0.2 eV (and state the value used originally), keeping all other parameters fixed. Plot the integrated THz spectral weight (0.1–3 THz) of E_FF from both current and magnetization sources as a function of tSO, and also plot max|Mz(t)| over the simulation window. If the THz enhancement and nonzero Mz appear for all tSO ≳ 0.01 eV with monotonically increasing weight, the qualitative claim is robust; if they only emerge above a threshold tSO > 0.1 eV inconsistent with NiO/Pt proximity, the central mechanism fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Central claim: in subgap-driven NiO, THz radiation becomes significant only when proximity Rashba SOC is present, and this explains NiO/Pt experiments without interlayer spin current. The only control parameter distinguishing the two cases is tSO in Eq. (4), yet its value is never given—no number, no units, no range, no derivation. The simulated spectra in Fig. 3(h) are therefore not reproducible, and the comparison between the plain and SOC cases is not quantitative. More seriously, the strength of tSO relative to t0 ≈ 1 eV, JH ≈ 1 eV, and J = 0.1 eV determines whether the observed Mz(t) and THz fields are realistic for NiO/Pt proximity. If the chosen tSO is far above the actual interfacial spin-orbit coupling (which DFT puts in the tens of meV range), the effect could be artificially inflated; if far below, the claimed enhancement might vanish. The absence of a scan over tSO leaves the central mechanism unfalsified. The paper's claim 'in full accord with experiments' thus rests on an unreported parameter.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies femtosecond-laser-driven NiO within a two-orbital Hubbard-Hund-Heisenberg model on a 4×2 ladder, with and without a Rashba-type spin-orbit coupling term intended to describe proximity to a Pt layer. Using exact diagonalization and tDMRG, it computes the time-dependent Néel vector, nonequilibrium magnetization, and far-field THz/high-harmonic radiation from both bond charge currents and the magnetic dipole. The central claims are: (i) subgap excitation produces nonclassical, length-only, non-rotating dynamics of the Néel vector and magnetization; (ii) THz emission from local charge currents and the magnetic dipole becomes significant only when Rashba SOC is switched on, which the authors argue explains NiO/Pt experiments without invoking an interlayer spin current; and (iii) above-gap excitation gives a sizable Néel-vector reduction and noninteger harmonics, while subgap excitation gives odd integer harmonics from currents and even integer harmonics from the magnetic dipole.","tokens_in":14630,"tokens_out":8042,"duration_ms":87498,"significance":"If the central mechanism is robust, the paper offers a conceptually useful alternative to the FM/HM picture of spintronic THz emission from antiferromagnetic insulators: THz radiation from NiO/Pt would originate from SOC-enabled local spin and charge dynamics inside NiO, rather than from an assumed interlayer spin current. The work also proposes a concrete, falsifiable prediction (even-order magnetic-dipole harmonics) and uses numerically exact time evolution for the chosen Hamiltonian, with parameters for U, t0, JH, and J inherited from prior DFT/DMFT work. The radiation calculation through the Jefimenko formula is a strength. However, the model is minimal—a single 8-site ladder, no explicit oxygen or Pt layers, an impurity field used to pin Néel order, and, most importantly, a Rashba hopping amplitude tSO whose numerical value is never specified. These gaps currently make the comparison between the plain and SOC cases qualitative and not fully reproducible.","major_comments":[{"comment":"The Rashba hopping amplitude tSO is never assigned a numerical value in eV or as a ratio to t0. The text states only that tSO is a 2×2 matrix hopping 'with values −itSO σ_y (itSO σ_x)' on horizontal (vertical) bonds; the scalar amplitude is absent. Since the central claim is that THz emission becomes significant only when this term is present, the entire plain-NiO vs SOC-NiO comparison rests on an unreported parameter. I request the value actually used, and a scan over at least the tens-of-meV range expected from DFT for AFI/HM interfaces (e.g., Ref. [42]), with the THz electric-field amplitude or power plotted versus tSO. Without this, Fig. 3(h) cannot be reproduced, and the statement 'in full accord with experiments' is not quantitatively grounded.","section":"Model and Methods, Eq. (4), and Fig. 3(h)"},{"comment":"The study uses one 4×2 ladder with two orbitals per site, an impurity field B_imp^z = 0.1 eV at site 1, no explicit oxygen atoms, and no explicit Pt layer. These simplifications are acknowledged, but the footnote in Ref. [75] argues that omitting the HM layer is acceptable because experiments see the strongest THz signal for the thinnest NiO. That argument does not establish that an interlayer spin current, oxygen-mediated exchange, or three-dimensional geometry leaves the spin-charge dynamics unchanged. Because the paper's central claim is to explain the experiments without an interlayer spin current, a finite-size and model-robustness check is needed. At minimum, the authors should show that length-only dynamics, SOC-induced Mz(t), and the THz enhancement survive in a larger ladder (or a different lattice shape) simulated with tDMRG, and that the results do not depend strongly on the impurity-field strength or position.","section":"Model and Methods, Eqs. (1)–(4), Fig. 1, and footnote [75]"},{"comment":"The central qualitative claim that the Néel vector and magnetization 'are changing length along the z-axis while not rotating at all' is supported only by the sentence that the x and y components are 'vanishingly small.' No time traces, Fourier amplitudes, or numerical upper bounds for Nx, Ny, Mx, and My are shown. Since this is the basis for labeling the dynamics 'nonclassical,' please provide these components for the subgap case and state the maximum ratios |N⊥|/|Nz| and |M⊥|/|Mz|. If the components are nonzero but small, the phrase 'not rotating at all' should be softened accordingly.","section":"Results and Discussion, Fig. 3(a)–(b), (e)–(f)"},{"comment":"The phrase 'in full accord with experiments [11,12]' is supported only by the qualitative presence of THz radiation in the SOC case. There is no overlay with the experimental THz spectra of Refs. [11,12], no clear statement of the simulated emission band within the 0.1–3 THz range, and no intensity ratio between plain NiO and NiO/HM. Because the SOC strength is unconstrained, this comparison is not falsifiable as written. Please specify which experimental observable is being matched and show the simulated E-field or power spectrum over the relevant THz window, rather than only an FFT cut whose normalization is not described.","section":"Figs. 3(d),(h) and 4(d),(h), and Refs. [11,12]"}],"minor_comments":[{"comment":"There are several typos: 'Ie case' should be 'In the case', 'lectrons' should be 'electrons', 'efects' should be 'effects', and 'only only when' should be 'only when'.","section":"Introduction and Fig. 3 caption"},{"comment":"The impurity field is described as a magnetic field on site i = 1, but Eq. (1) writes it as g μ_B B_imp^z s^z_{1α}, which is an orbital-resolved operator. Please clarify whether the field couples equally to both orbitals and whether it is meant to represent a local symmetry-breaking field rather than a physical magnetic impurity.","section":"Eq. (1) and surrounding text"},{"comment":"The dimensionless pump intensity zmax = e a0 A_max/ℏ is used, but no corresponding peak electric field in V/Å is given. Stating the peak field would help experimentalists gauge the intensity regime of the calculation.","section":"Model and Methods, after Eq. (3)"},{"comment":"The noninteger harmonics for above-gap pumping are reported as an unusual result, but the explanation is limited to a brief reference to multiple Floquet-state populations. A sentence connecting the present observation to the mechanism of Ref. [29] would make the claim more useful to readers.","section":"Results and Discussion, noninteger harmonics"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the missing numerical value of tSO is a substantive reproducibility gap, not a presentation issue, and it directly affects the paper's central claim. The 'nonclassical' language relies on the authors' prior work (Refs. [71,76]), but the present simulation is self-contained in showing length-only dynamics, so this is not circular. The agreement with experiments is qualitative and post hoc; the even-HHG prediction is the most falsifiable element and should be highlighted. There is no indication of a novelty disclosure problem, but the lack of any system-size check makes robustness hard to judge."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a plausible and well-executed small-system numerical experiment, but the central quantitative claim is underdetermined by a missing parameter value, and the paper oversells the agreement. Worth refereeing with major revision.\n\nWhat is actually new: the combination of the established 2HHH model with ED/tDMRG to compute THz and high-harmonic radiation from a NiO-like ladder with Rashba SOC, plus the concrete proposal that THz emission from NiO/HM does not require an interlayer spin current. The even high harmonics from the magnetic dipole channel and the noninteger harmonics for above-gap pumping are new spectral predictions. The dynamics itself – Neel vector and magnetization changing length along z without rotating – is largely an extension of the authors' prior work (refs [71,76]) to this setup, so the conceptual novelty is modest.\n\nWhat is done well: the ED time evolution is numerically exact for the model, parameters are mostly inherited from prior DFT/DMFT work, and the radiation calculation via Jefimenko formulas is concrete. The qualitative contrast – THz appears when Rashba SOC is switched on and not otherwise – is suggestive and aligns with the experimental phenomenology.\n\nThe soft spots are real. Most important: the Rashba hopping tSO in Eq. (4) is never given a number, unit, or range. I checked; the parameter list ends at J and B_imp. The entire plain-NiO vs NiO/HM comparison is controlled by this one parameter, and the label \"full accord with experiments\" rests on an unreported number. A simple scan over tSO, or at least a stated value with a sensitivity check, is needed before the mechanism can be trusted. Also, the system is a single 8-site ladder with no finite-size or convergence analysis; the impurity field B_imp = 0.1 eV at site 1 is an ad hoc way to pin Neel order and could affect the dynamics. No code or data is shipped. Some of these are addressable, but the missing tSO is not cosmetic: it is load-bearing.\n\nThe comparison to experiment is also qualitative. The paper says \"only when proximity SOC is present... in full accord with experiments,\" but the theoretical spectrum is not directly compared to the measured waveform or amplitude. That language should be tempered.\n\nBottom line: this is a serious proposal, and the radiation signatures (even HHG, THz only with SOC) give experimentalists something to look for. But in its current form the central claim is underdetermined. I would send it to a referee, but with a clear request to specify tSO, add a scan or justify the value, and soften the \"full accord\" phrasing.\n\nWho should read it: anyone working on THz spintronics from antiferromagnets, or on nonclassical spin dynamics in pumped Mott insulators. I would bring it to reading group and cite it if the tSO issue is fixed; as is, I would treat it as a promising preprint.","headline":"Plausible mechanism, but the missing Rashba hopping value makes the central quantitative claim untestable as written.","tokens_in":15161,"tokens_out":2876,"would_cite":true,"duration_ms":30371,"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":"Subgap-laser-driven NiO emits THz only when heavy-metal proximity adds spin-orbit coupling, and its Neel vector changes length without rotating.","keywords":["antiferromagnetic insulator","THz spintronics","Neel vector","nonequilibrium magnetization","Hubbard-Hund-Heisenberg model","Rashba spin-orbit coupling","high-harmonic generation","femtosecond laser pulse"],"falsifier":"A direct test would be to measure THz emission from NiO/Pt bilayers while inserting an atomically thin insulating spacer between NiO and Pt: if emission persists when proximity spin-orbit coupling is suppressed but spin currents might still flow, or vanishes while spin currents would still be expected, the mechanism would be contradicted. A simpler check is to scan NiO thickness and compare the THz intensity with the predicted dominance of the thinnest proximitized layers.","tokens_in":14185,"feed_emoji":"🧲","tokens_out":6530,"duration_ms":61800,"temperature":0.7,"pith_summary":"This paper tries to establish a microscopic mechanism for THz emission from NiO/Pt bilayers driven by subgap femtosecond laser pulses, replacing the borrowed ferromagnet/heavy-metal picture of an interlayer spin current. Using a two-orbital Hubbard-Hund-Heisenberg model on an 8-site ladder, with Rashba spin-orbit coupling added to mimic the platinum layer's proximity effect, it finds that both the Neel vector and the nonequilibrium magnetization change only in length along the z-axis and do not rotate, a dynamics no Landau-Lifshitz equation can describe. It further finds that THz radiation from both bond charge currents and the time-dependent magnetization is significant only when the proximity spin-orbit coupling is present, matching experiments on NiO/Pt. A sympathetic reader would care because this removes the need for speculative interlayer spin currents and offers a calculable route from strongly correlated many-body dynamics to radiated electromagnetic fields.","feed_headline":"NiO emits THz only with heavy-metal proximity, simulations show","feed_subtitle":"Exact simulation traces emission to spin-orbit-coupled dynamics inside NiO, not a spin current.","key_machinery":"The carrying object is a two-orbital Hubbard-Hund-Heisenberg (2HHH) model on a 4x2 ladder of Ni sites, with oxygen omitted but its mediated interactions encoded in realistic parameters ($U \\approx 8$ eV, $t_0 \\approx 1$ eV, $J_H \\approx 1$ eV). A Rashba-type spin-orbit hopping $t_{SO}$ is added to represent proximity to a heavy metal; the laser enters as a Peierls phase multiplying the hoppings, and a static impurity field at site 1 selects Neel order. Time evolution is performed with numerically exact exact-diagonalization and tensor-network (tDMRG) methods, and the emitted electric field is obtained by feeding the bond current and magnetization into the Jefimenko far-field formula. This machinery lets the authors compute both the many-body spin-charge dynamics and the radiation it produces, and it is the source of the claim that length-changing, nonrotating spin dynamics generates the THz signal.","core_discovery":"The central claim is that subgap femtosecond pumping of NiO produces a highly nonclassical response: the Neel vector (the difference of the two sublattice magnetizations) and the nonequilibrium magnetization (their sum) both shrink or grow along the out-of-plane direction while never rotating. The magnetization stays exactly zero in plain NiO but becomes nonzero once Rashba spin-orbit coupling, modeling the proximity of a heavy-metal layer, is switched on. The paper computes far-field radiation from both time-dependent bond charge currents and the second time derivative of the magnetization, and finds that THz-frequency radiation from either source is significant only in the proximitized case, in full accord with NiO/Pt experiments. It therefore concludes that THz emission from antiferromagnetic-insulator/heavy-metal bilayers can arise from spin-orbit-enabled local charge and magnetization dynamics inside the antiferromagnet, with no interlayer spin current or inverse spin Hall conversion required. Above the THz range, the spectra show odd integer harmonics from the current channel, even integer harmonics from the magnetic-dipole channel, and unusual noninteger harmonics for above-gap pumping.","pith_inferences":["Beyond the paper, replacing the platinum layer with a different heavy metal should change the THz intensity through the strength of the induced Rashba coupling rather than through the metal's spin Hall angle, a testable distinction from the spin-current picture.","The model implies a strong thickness dependence: THz emission should be largest for the thinnest NiO films, where the proximity effect reaches the whole antiferromagnet; inserting a few atomic layers of a spacer should suppress emission even if spin currents could still pass.","The appearance of noninteger harmonics under above-gap pumping could be used as a probe of multi-Floquet-state population; a pump-intensity scan should show whether their spectral positions are tunable."],"forward_implications":["THz emission from NiO/heavy-metal bilayers can be explained without any interlayer spin current, so experiments should focus on spin-orbit-induced charge and magnetization dynamics inside the antiferromagnet rather than on inverse spin Hall conversion.","Because magnetic-dipole radiation can rival or exceed bond-current radiation in the proximitized case, the common assumption that charge currents always dominate THz emission from magnetic bilayers does not carry over to antiferromagnetic insulators.","Even-integer high harmonics in the magnetic-dipole channel and odd-integer harmonics in the current channel give a symmetry-based way to identify which source is radiating.","Above-gap pumping partially destroys the Neel order (up to 15 percent) while subgap pumping leaves it almost intact, meaning the two pumping regimes drive qualitatively different magnetic responses."],"supporting_citations":[{"why":"Supplies the two-orbital Hubbard-Hund-Heisenberg model and the above-gap NiO pump experiment that this paper extends to subgap pumping.","marker":"[2]"},{"why":"The subgap NiO/HM THz experiment whose observed emission the paper aims to explain without invoking an interlayer spin current.","marker":"[11]"},{"why":"The NiO/Pt THz emission experiment that motivates adding Rashba SOC and reports the strongest emission for the thinnest NiO layer.","marker":"[12]"},{"why":"Establishes the TDDFT-plus-Jefimenko approach for computing electromagnetic radiation from spintronic THz emitters, adapted here to an antiferromagnetic insulator.","marker":"[24]"},{"why":"First-principles evidence for Rashba-type spin splitting at heavy-metal interfaces, used to justify the proximity SOC term in the model.","marker":"[10]"},{"why":"First-principles study of spin pumping from an antiferromagnetic insulator spin-orbit-proximitized by a heavy metal, supporting the proximity SOC picture for NiO/HM.","marker":"[42]"},{"why":"Reports sub-gap optical excitation effects in NiO, such as midgap states and gap reduction, motivating coupled spin-charge dynamics.","marker":"[23]"}],"fun_headline_variants":["NiO emits THz only with heavy-metal proximity","No spin current: NiO's THz from local magnetization dynamics","Neel vector shrinks, doesn't rotate, under ultrafast light","Heavy-metal layer unlocks NiO's THz emission","THz from NiO: spin-orbit coupling does it locally"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole mechanism rests on an 8-site two-orbital ladder with no oxygen atoms, no explicit platinum layer, and a single impurity field reproducing the essential physics of a NiO/Pt bilayer; if omitted interfacial spin currents, oxygen-mediated exchange, or three-dimensional geometry materially changes the spin-charge dynamics, the mechanism as stated would not survive.","fun_headline_variants_meta":{"raw":{"variants":["NiO emits THz only with heavy-metal proximity","No spin current: NiO's THz from local magnetization dynamics","Neel vector shrinks, doesn't rotate, under ultrafast light","Heavy-metal layer unlocks NiO's THz emission","THz from NiO: spin-orbit coupling does it locally"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000916,"raw_usage":{"total_tokens":4036,"prompt_tokens":1155,"completion_tokens":2881,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":771,"completion_tokens_details":{"reasoning_tokens":2794}},"tokens_in":771,"tokens_out":2881,"duration_ms":23078,"temperature":1.0,"reasoning_tokens":2794,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T17:30:01.380316+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would be to measure THz emission from NiO/Pt bilayers while inserting an atomically thin insulating spacer between NiO and Pt: if emission persists when proximity spin-orbit coupling is suppressed but spin currents might still flow, or vanishes while spin currents would still be expected, the mechanism would be contradicted. A simpler check is to scan NiO thickness and compare the THz intensity with the predicted dominance of the thinnest proximitized layers.","supporting_citations":[{"cited_title":"Kefayati and B","cited_arxiv_id":null,"evidence_quote":"Establishes the TDDFT-plus-Jefimenko approach for computing electromagnetic radiation from spintronic THz emitters, adapted here to an antiferromagnetic insulator."},{"cited_title":"Dolui, A","cited_arxiv_id":null,"evidence_quote":"First-principles study of spin pumping from an antiferromagnetic insulator spin-orbit-proximitized by a heavy metal, supporting the proximity SOC picture for NiO/HM."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports sub-gap optical excitation effects in NiO, such as midgap states and gap reduction, motivating coupled spin-charge dynamics."}],"review_version":1}