{"id":"30387f5a-5630-4d0b-ab5e-ab9a0f5b0fc7","arxiv_id":"2608.06129","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In MgO, tuning pulse duration and intensity selects between multicycle band-climbing and subcycle multiband pathways, extending high-harmonic emission to 50 eV.","lead":"Experiments and simulations show that changing laser pulse duration from 29 to 5 femtoseconds switches how electrons in magnesium oxide climb conduction bands, and can push emitted extreme-ultraviolet light to 25-50 eV. This makes pulse duration a practical control knob for compact solid-state XUV light sources.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Threshold-model delta-windows set both It1 and It2; the 5-fs onset is stated as ~30 TW/cm2 vs predicted 22, so the quantitative pathway-selective claim is not yet secure.","rationale":"Read in good faith: the paper combines a clean pulse-duration/intensity experiment with first-principles TDDFT and a deliberately simple threshold model. The experimental trends, 29-fs high plateau, 17-fs loss of spectral clarity, 5-fs reappearance, are reproduced in simulations using the measured pulse shapes, which is independent support. The weakest link is the quantitative mapping between the threshold model and the data. Both It1 and It2 are derived from delta-window estimates; these windows are not measured or computed to a stated tolerance. The SI dipole plot is visual only. Moreover, the experimental 5-fs onset is quoted as ~30 TW/cm2, not the predicted 22. This discrepancy is either a detection-threshold effect or a model error; without a decision, the pathway-selective claim is conditional. The suggested intensity scan in TDDFT would settle whether the 22 TW/cm2 subcycle threshold is real in the microscopic model. If it is, the paper's central interpretation survives; if not, the two-threshold picture needs revision. This does not call into question the authors' honesty or the qualitative pulse-duration control; it asks for one quantitative check on the load-bearing assumption.","tokens_in":15766,"tokens_out":13908,"duration_ms":135587,"concrete_test":"Run a TDDFT (Octopus) intensity scan using the experimental 5-fs pulse shape at peak intensities 18, 22, 26, 30, and 34 TW/cm2 and the experimental 29-fs pulse at 4, 5.5, 7, and 10 TW/cm2, using the same numerical settings as in SI Sec. III. For each run, extract the yield in the 25-50 eV window and the Gabor time-frequency emission profile around the pulse center. If the 5-fs high-energy plateau turns on near 22 TW/cm2, the It2 formula is confirmed and the experimental ~30 onset is a detection or signal-level effect; if it requires ~30 TW/cm2, the threshold model overestimates the subcycle coupling and must be revised. The same runs should also output k-resolved interband population transfer to identify the actual transition windows and compare them quantitatively with delta_Gamma and delta_X.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that pulse duration selects between a multicycle (It1 ~ 5.5 TW/cm2) and a subcycle (It2 = 4It1 ~ 22 TW/cm2) pathway depends on the threshold vector potentials A_t1 = (2*pi/a - delta_Gamma - delta_X)/2 and A_t2 = 2*A_t1. These are not fitted; they use the ad hoc windows delta_j ~ sqrt(m*_j omega_L) around Gamma and X. The SI (Fig. 3) presents transition dipoles only visually as corroboration, with no quantitative comparison. The windows matter because delta_Gamma + delta_X ~ 0.38 a.u.-1 is a large fraction of the 2*pi/a ~ 0.79 a.u.-1 Gamma-X distance: a 30% change in the sum shifts A_t1 by ~15% and the intensity threshold by ~30%, and a factor-of-two change in the windows moves It1 by roughly a factor of four. A concrete symptom appears in the experimental section: the 5-fs high-energy plateau is said to reappear only above ~30 TW/cm2 (Fig. 2c), whereas the model gives It2 ~ 22 TW/cm2. The later text writes 'once this intensity threshold is exceeded' as if 30 and 22 coincide, but they do not. If the actual k-resolved transition regions are broader or narrower than delta_j, or if the true subcycle onset is 30 rather than 22, the quantitative match that anchors the pathway assignment is degraded. The qualitative pulse-duration trends would survive, since the TDDFT simulation reproduces the spectra, but the specific mechanism, the central novelty, would need revision.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports an experimental and theoretical study of high-harmonic generation (HHG) in MgO driven by 29 fs, 17 fs, and 5 fs pulses at intensities of 0.8–74 TW/cm^2. The authors observe that the 25–50 eV plateau appearing for 29 fs pulses above ~6 TW/cm^2 disappears at 17 fs and reappears for 5 fs pulses only above ~30 TW/cm^2. TDDFT simulations using the experimentally retrieved pulse shapes reproduce the main spectral trends. The paper attributes the long-pulse plateau to cumulative multicycle carrier transfer between Γ and X with threshold I_t1 ≈ 5.5 TW/cm^2, and the short-pulse reappearance to a subcycle Γ–X transfer channel with threshold I_t2 = 4 I_t1 ≈ 22 TW/cm^2. Pulse duration is proposed as a control knob that selects between these pathways.","tokens_in":16180,"tokens_out":6622,"duration_ms":40837,"significance":"Confidence in the experimental trends is increased by the TDDFT simulations that use measured pulse shapes rather than idealized pulses, and by the time-frequency analysis supporting different emission times for the two plateaus. The comparison with a Fourier-limited pulse is a useful falsifiable prediction. If the threshold model is correct, the work would establish pulse duration as a decisive parameter in solid-state HHG and would connect the 50 eV cutoff to band-structure geometry. The paper is less strong on quantitative experimental characterization: no error bars are given for the spectra or intensity calibration, and the predicted subcycle threshold (22 TW/cm^2) is not matched by the observed onset (~30 TW/cm^2). The conceptual framework is appealing, but the quantitative pathway-selective claim currently rests on an ad hoc momentum-window estimate.","major_comments":[{"comment":"The predicted subcycle threshold I_t2 = 4 I_t1 ≈ 22 TW/cm^2 is not equal to the experimentally reported reappearance threshold of ~30 TW/cm^2. The text states that \"once this intensity threshold is exceeded, the high-energy plateau reappears\", implying that the observed onset is I_t2, but 30 TW/cm^2 is ~36% larger than 22 TW/cm^2. Since the central novelty is the assignment of the reappearing plateau to the subcycle A_t2 channel, this discrepancy must be resolved: either the experiment should measure the 5 fs intensity dependence with smaller steps to locate the onset, or the model should be revised to explain why the onset occurs at ~30 TW/cm^2. The 5 fs TDDFT simulation at 41.4 TW/cm^2 does not discriminate, as it is well above both values.","section":"Figs. 2(c) and 5; paragraph defining A_t2"},{"comment":"The thresholds A_t1 = (2π/a − δ_Γ − δ_X)/2 and A_t2 = 2A_t1 depend on the ad hoc momentum-space windows δ_j ≈ sqrt(m*_j ω_L). These windows are estimated, not derived or fitted, and the SI (Sec. III, Fig. 3) offers only a visual \"corroboration\" from transition dipoles without a quantitative comparison. Given that δ_Γ + δ_X ≈ 0.38 a.u.^{-1} is comparable to 2π/a ≈ 0.79 a.u.^{-1}, the threshold intensities are highly sensitive to the window choice: a 30% change in the sum alters I_t1 by several tens of percent, and a factor-of-two change can move I_t1 by more than a factor of two. The claimed match between I_t1 ≈ 5.5 TW/cm^2 and the 29 fs onset is therefore not an independent confirmation. Please provide a quantitative, dipole-weighted estimate of the effective transition regions or a robustness study over a range of δ_j.","section":"Fig. 4 and SI Sec. III, Fig. 3"},{"comment":"No error bars, repeated measurements, or shot-to-shot statistics are reported for the harmonic spectra, and no uncertainty is given for the intensity calibration (pulse energy, focal spot size, pulse duration). The central quantitative claims are threshold intensities (~6 and ~30 TW/cm^2), so without an uncertainty estimate it is not possible to judge whether the 22 vs 30 TW/cm^2 discrepancy is statistically significant. The authors should report at least the spread of measured spectra and the estimated systematic uncertainty of the intensity scale.","section":"Fig. 2 and Table I"}],"minor_comments":[{"comment":"The formula is attributed to \"Lamor\" rather than \"Larmor\" in both the main text and the supplementary material.","section":"Main text and SI Eq. (1)"},{"comment":"The periodicity assumption j(t) ≈ j(t + 2π/ω_L) is an approximation, and the SI acknowledges that long-trajectory accumulation and ground-state depletion break it; this caveat should appear in the main text where the intercycle analysis is used to interpret Fig. 2.","section":"SI Sec. II"},{"comment":"The transition dipole plots are normalized per panel, which prevents cross-panel comparison of absolute dipole strengths; an absolute scale would make the corroboration of δ_Γ and δ_X more convincing.","section":"SI Sec. III, Fig. 3"},{"comment":"The statement that all conclusions are robust to the choice of exchange-correlation functional is not accompanied by any LDA results; either show the supporting calculations or soften the claim.","section":"SI Sec. III"},{"comment":"The phrase \"before decoherence can suppress coherent emission\" is presented as an explanation, but the TDDFT calculations do not include electron–electron scattering decoherence; this is an inference from timing rather than a computed result.","section":"Abstract and conclusion"}],"recommendation":"major_revision","confidential_remarks":"The manuscript relies heavily on self-citations (e.g., Refs. 5, 32, 34, 36, 38, and related SI entries), which may be defensible given the close topical connection, but the editor may wish to check that the novelty relative to Refs. 32 and 5 is sufficiently demarcated. The paper would also benefit from a data availability statement for the experimental spectra and simulation outputs."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the pulse-duration control of the MgO HHG cutoff is real and new, and the TDDFT reproduction with measured pulse shapes is the strongest part of the paper. The two-threshold mechanism is plausible but less secure than the abstract implies, mainly because the predicted subcycle threshold (22 TW/cm2) doesn't match the observed onset (~30 TW/cm2) and the delta-window estimates are too rough to settle it.\n\nWhat's genuinely new: their own Ref [32] already reported a 50 eV cutoff in MgO, but the systematic pulse-duration scan—plateau present at 29 fs, gone at 17 fs, back at 5 fs above a higher intensity—is new. The reappearance at 5 fs is the key observation, and it's reproduced by TDDFT using the actual measured pulse shapes. That's strong evidence that the trends come from electronic dynamics, not propagation artifacts. The threshold model is not fitted; it derives thresholds from lattice constant and effective masses, and the time-frequency analysis gives independent support for the multicycle vs subcycle distinction. Credit where due: this is a well-executed experiment-theory package.\n\nSoft spots. The quantitative match for the subcycle threshold is loose: the model gives It2 ~ 22 TW/cm2, but the 5-fs plateau is observed only above ~30 TW/cm2. The text says 'once this intensity threshold is exceeded' as if 22 and 30 were the same. They are not, and the paper never addresses the gap. The threshold calculation rests on ad hoc momentum windows delta_j ~ sqrt(m* omega_L) around Gamma and X; the SI shows transition dipoles visually but doesn't quantify how well those windows match the actual k-resolved transition regions. A 30% change in the sum of the windows shifts the intensity threshold by ~30%, so the claimed It2 is not a sharp prediction. Also, there are no error bars or shot-to-shot statistics, so the observed onset intensities are approximate. Finally, the 'before decoherence can suppress coherent emission' phrasing is an assertion; no dephasing measurement or timescale is given.\n\nMy take: the qualitative pathway-selective picture—pulse duration selects between cumulative multicycle transfer and subcycle multiband excitation—is credible and well supported by the experiment-theory agreement. The specific quantitative thresholds and the clean separation of the two channels are not yet proven. That's a fixable weakness, not a fatal one.\n\nThis paper deserves a serious referee. I'd send it out, and I'd ask the authors to (1) show error bars or at least a conservative intensity uncertainty, (2) address the 22 vs 30 discrepancy directly, and (3) either strengthen the delta-window justification with quantitative dipole comparisons or soften the claim. If those are addressed, it's a nice contribution to solid-state HHG. For my own work, I'd cite the experimental pulse-duration dependence, though not the specific threshold numbers.","headline":"Solid experiment-theory package on pulse-duration control of the HHG cutoff in MgO; the two-threshold mechanism is plausible but the quantitative subcycle threshold is loose, so the pathway assignment is suggestive rather than proven.","tokens_in":16700,"tokens_out":2681,"would_cite":true,"duration_ms":24050,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.65.Ky","78.47.J-"],"model":"deepseek-v4-flash","headline":"This paper demonstrates that varying laser pulse duration and intensity can select distinct multiband electron pathways in MgO, extending high-harmonic emission to photon energies of 25–50 eV before decoherence suppresses it.","keywords":["high-harmonic generation","solid-state HHG","MgO","pulse duration control","subcycle electron dynamics","multiband dynamics","extreme-ultraviolet emission","threshold vector potential"],"falsifier":"Measure the time-frequency map of MgO harmonic emission for a 5 fs pulse at intensities just above and below $I_{t2} \\approx 22\\ \\mathrm{TW/cm^2}$: the subcycle model predicts the 20–30 eV burst appears and disappears within a half-cycle immediately at the intensity threshold, whereas the multicycle model predicts a delayed build-up over several cycles; a build-up delay comparable to the pulse duration would refute the subcycle channel.","tokens_in":15605,"feed_emoji":"⚡","tokens_out":13319,"duration_ms":100124,"temperature":0.7,"pith_summary":"High-order harmonic generation in solids usually treats intensity, wavelength, and crystal orientation as the adjustable parameters. This paper adds pulse duration as an equally decisive knob: by compressing the driving pulse from 29 fs to 5 fs while raising peak intensity, it separates two distinct excitation pathways in the insulator MgO. Long pulses above about $6\\ \\mathrm{TW/cm^2}$ accumulate carrier transfers over many optical cycles into higher conduction bands, pushing emission to 25–50 eV. Short pulses above about $22\\ \\mathrm{TW/cm^2}$ instead complete the same high-energy transfer within a single half-cycle, before decoherence can erase the coherent emission. The result matters because it identifies pulse duration as a practical control for reaching extreme-ultraviolet photon energies in compact solid-state sources.","feed_headline":"Pulse duration pushes solid-state harmonics to 50 eV","feed_subtitle":"Tuning pulse length and intensity selects fast multiband electron paths before decoherence kills the emission.","key_machinery":"The load-bearing machinery is a momentum-space threshold model: efficient carrier transfer into higher conduction bands is assumed to occur only within narrow windows of radius $\\delta_\\Gamma \\approx 0.05\\pi\\ \\mathrm{a.u.}^{-1}$ around the $\\Gamma$ point (valence-to-first-band excitation) and $\\delta_X \\approx 0.07\\pi\\ \\mathrm{a.u.}^{-1}$ around the $X$ point (transfer into the second and third conduction bands), with widths set by $\\delta \\approx \\sqrt{m^{*}\\omega_L}$. The acceleration theorem $k(t)=k_0+A(t)$ turns the laser vector potential into a momentum displacement, so the existence of a half-cycle $\\Gamma$-to-$X$ trajectory is fixed by the threshold $A_{t1} = (2\\pi/a - \\delta_\\Gamma - \\delta_X)/2$, and the full-cycle trajectory by $A_{t2} = 2A_{t1}$. These thresholds map to intensities $I_{t1} \\approx 5.5\\ \\mathrm{TW/cm^2}$ and $I_{t2} \\approx 22\\ \\mathrm{TW/cm^2}$, which the paper identifies with the observed onset of the high-energy plateau for 29 fs and 5 fs pulses. A complementary intra- and intercycle factorization of the harmonic yield explains why shorter pulses broaden and weaken harmonic peaks via an intercycle interference factor that scales with $N_c^2$ and widens as $(N_c^2-1)^{-1/2}$.","core_discovery":"The paper demonstrates that the harmonic cutoff in MgO is not simply set by peak intensity: pulse duration selects which multiband pathway carries electrons to high recombination energies. Below the first threshold, electrons climb the lowest conduction band and recombine with the usual low-order spectrum. Once the vector potential reaches $A_{t1} = (2\\pi/a - \\delta_\\Gamma - \\delta_X)/2$, electrons can traverse from the $\\Gamma$ point to the $X$ point within a half-cycle and, over successive cycles, transfer into higher conduction bands, creating a broadened high-energy plateau from 25 eV toward 50 eV. With few-cycle pulses that reach intensities above $I_{t2} \\approx 22\\ \\mathrm{TW/cm^2}$ (the intensity corresponding to $A_{t2}=2A_{t1}$), the same high-band transfer occurs within a single subcycle excursion, so the high-energy plateau reappears with narrow harmonic peaks. The intermediate 17 fs pulse exceeds the first threshold but its shorter interaction time broadens the harmonic peaks so much that the plateau is not observable, explaining the non-monotonic reappearance.","pith_inferences":["A testable extension the paper leaves open is carrier-envelope phase control: for a 5 fs pulse, the CEP should shift the subcycle $\\Gamma$-to-$X$ trajectory and modulate the high-energy plateau on a half-cycle timescale, an effect the CEP-averaged simulations do not address.","If the threshold picture holds for other crystals, it gives a material-search rule: reducing the reciprocal-space distance between a strong excitation point and the first high-band transfer point (e.g., via strain or different lattices) should proportionally lower $I_{t1}$ and $I_{t2}$, making 50 eV-class emission easier to reach.","The subcycle pathway's speed implies it should be far less sensitive to temperature and phonon scattering than the multicycle pathway, so a temperature-dependent HHG measurement could separate the two channels without changing the pulse shape."],"forward_implications":["At a fixed intensity, compressing the pulse from 29 fs to 17 fs broadens harmonic peaks by decreasing the number of coherent cycles $N_c$, which is why the 25–50 eV plateau disappears in the 17 fs data even though the field exceeds the first threshold.","For 29 fs pulses, exceeding $I_{t1} \\approx 5.5\\ \\mathrm{TW/cm^2}$ opens a multicycle pathway that progressively transfers carriers into higher conduction bands, producing the 25–50 eV plateau at the cost of suppressing low-order harmonic yield.","For 5 fs pulses, exceeding $I_{t2} \\approx 22\\ \\mathrm{TW/cm^2}$ opens a subcycle pathway: a single accelerated wavepacket visits both $\\Gamma$ and $X$ within one half-cycle, re-creating the high-energy plateau with harmonic peaks as narrow as the low-order ones.","Pulse duration thus becomes a design parameter for solid-state XUV sources, alongside intensity, wavelength, and crystal orientation.","The same threshold logic predicts that along the $\\Gamma$–$K$ direction the onset intensities rise to about $7.3$ and $29\\ \\mathrm{TW/cm^2}$, consistent with the much weaker high-energy emission measured there."],"supporting_citations":[{"why":"Supplies the crystal-momentum-resolved multiple-plateau model that the paper adapts to MgO with Gamma- and X-point transfer windows.","marker":"[52]"},{"why":"Reports the earlier observation of a 50 eV cutoff in MgO that the present work explains through pulse-duration-dependent pathways.","marker":"[32]"},{"why":"Provides the intra- and intercycle factorization of intraband harmonic generation used to interpret peak narrowing and broadening.","marker":"[36]"},{"why":"Introduces the intra- and intercycle interference analysis of strong-field ionization that the solid-state factorization generalizes.","marker":"[34]"},{"why":"Established the first observation of high-harmonic generation in a bulk crystal, the foundational phenomenon studied here.","marker":"[1]"},{"why":"Supports the acceleration-theorem picture of intraband motion and recombination used to assign harmonic energies.","marker":"[15]"},{"why":"Provides the semiclassical analysis of high-harmonic generation in bulk crystals underlying the trajectory and cutoff arguments.","marker":"[53]"},{"why":"Supplies the first-principles simulation approach from which the computed harmonic spectra and transition dipole elements are obtained.","marker":"[38]"},{"why":"Contains the pulse reconstructions, time-frequency analysis, numerical methods, and Gamma-K comparison on which the paper's claims rely.","marker":"[33]"}],"fun_headline_variants":["Pulse-length knob selects electron paths to push harmonics to 50 eV","Subcycle multiband dynamics extend solid-state harmonic cutoff","Pathway-selective pulse control extends harmonic cutoff to 50 eV","Few-cycle pulses exploit subcycle paths to reach 50 eV harmonics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes that efficient transfer of electrons into higher conduction bands happens only within narrow momentum windows of width $\\delta \\approx \\sqrt{m^{*}\\omega_L}$ around the $\\Gamma$ and $X$ points, with the high-energy plateau threshold set by those windows; if the real transfer regions are wider or the transition dipoles differ, the match between $I_{t1} \\approx 5.5\\ \\mathrm{TW/cm^2}$ and the observed onset is largely coincidental.","fun_headline_variants_meta":{"raw":{"variants":["Pulse-length knob selects electron paths to push harmonics to 50 eV","Subcycle multiband dynamics extend solid-state harmonic cutoff","Pathway-selective pulse control extends harmonic cutoff to 50 eV","Few-cycle pulses exploit subcycle paths to reach 50 eV harmonics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0009,"raw_usage":{"total_tokens":3867,"prompt_tokens":930,"completion_tokens":2937,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":546,"completion_tokens_details":{"reasoning_tokens":2872}},"tokens_in":546,"tokens_out":2937,"duration_ms":16343,"temperature":1.0,"reasoning_tokens":2872,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:22:11.861283+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the time-frequency map of MgO harmonic emission for a 5 fs pulse at intensities just above and below $I_{t2} \\approx 22\\ \\mathrm{TW/cm^2}$: the subcycle model predicts the 20–30 eV burst appears and disappears within a half-cycle immediately at the intensity threshold, whereas the multicycle model predicts a delayed build-up over several cycles; a build-up delay comparable to the pulse duration would refute the subcycle channel.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the crystal-momentum-resolved multiple-plateau model that the paper adapts to MgO with Gamma- and X-point transfer windows."},{"cited_title":"Allegre, J","cited_arxiv_id":null,"evidence_quote":"Reports the earlier observation of a 50 eV cutoff in MgO that the present work explains through pulse-duration-dependent pathways."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the intra- and intercycle factorization of intraband harmonic generation used to interpret peak narrowing and broadening."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the intra- and intercycle interference analysis of strong-field ionization that the solid-state factorization generalizes."},{"cited_title":"Vampa, C","cited_arxiv_id":null,"evidence_quote":"Provides the semiclassical analysis of high-harmonic generation in bulk crystals underlying the trajectory and cutoff arguments."},{"cited_title":"Tancogne-Dejean, M","cited_arxiv_id":null,"evidence_quote":"Supplies the first-principles simulation approach from which the computed harmonic spectra and transition dipole elements are obtained."},{"cited_title":"[5,34-48], for details on (I) pulse measurement, (II) intra- and intercycle analysis, (III) numerical methods, (IV) time-frequency analysis, and (V) comparison with theΓ−Kcase","cited_arxiv_id":null,"evidence_quote":"Contains the pulse reconstructions, time-frequency analysis, numerical methods, and Gamma-K comparison on which the paper's claims rely."}],"review_version":1}