{"id":"e438fe93-1dfc-4501-b289-93660eb4e517","arxiv_id":"2507.18019","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Electron occupation dynamics, especially delayed recombination from symmetry-forbidden transitions, controls and can broaden attosecond pulses emitted by bulk silicon under intense laser light.","lead":"Using simulations of bulk silicon driven by intense mid-infrared laser pulses, the authors show that the timing of attosecond pulses emitted from solids is set by how electrons populate different energy bands over time. They find that symmetry-forbidden optical transitions delay part of the emission and can broaden the attosecond pulse train, which matters for building solid-state attosecond sources.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Appendix C's diagonal-density approximation excludes interband coherences, so the central 'occupation-driven' mechanism is not actually derived; the SBE's own emission (Eq. 19) is coherence-based, not occupation-based.","rationale":"The paper's title and abstract make a strong causal claim: the temporal structure of solid-state attosecond pulses is 'occupation-driven,' and the analytical theory 'captures' this. I trace the derivation in Appendix C. The key step (Eq. 11) replaces the full one-body density n(r,t)=|Ψ|^2 by its diagonal part in an adiabatic basis, Σ O |φ|^2, discarding all off-diagonal terms. These off-diagonal terms are exactly the interband coherences that appear in the SBE as p (Eqs. 14-16) and that generate the macroscopic polarization P (Eq. 19). Equation (12), derived under this approximation, yields an APT that depends only on band-diagonal forces and the trivial N_e E(t) term; it cannot describe interband recombination. The SBE used for validation, however, computes emission from P and J (Eq. 21), so the delayed emission comes from coherences, not occupations. The occupation delay shown in Fig. 6(d) is a downstream effect of coherence-driven population transfer (Eqs. 17-18). Thus the central theoretical foundation of the paper is not merely incomplete—it is in tension with the paper's own numerical model. This is an internal inconsistency, not a matter of consensus. If the test confirms that the FWHM broadening disappears without P(t), the paper's interpretation must be revised from 'occupation-driven' to 'coherence-driven emission, with occupation as a proxy.' Given that the title and abstract assert the former, the paper needs a major revision, but the underlying numerical findings (nonmonotonic FWHM) may still be valid; therefore I keep the reader's CONDITIONAL verdict, i.e., UNCHANGED. The reader's weakest assumption (TDM k-dependence) is a different and secondary issue; it concerns the accuracy of the SBE's input, not the logical status of the theory.","tokens_in":15895,"tokens_out":13606,"duration_ms":148841,"concrete_test":"In the Appendix E SBE, suppress the interband polarization p_{k}^{he}(t) in the emission formula while retaining the full population dynamics, and recompute the APT FWHM versus intensity from the intraband current J(t) alone. If the nonmonotonic broadening and the peak delay vanish, the effect is coherence-driven, contradicting the occupation-only theory of Eq. (12). If they persist, occupation-driven intraband emission is sufficient and the concern does not land.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Appendix C's central step (Eq. 11) expands the time-dependent density as n(r,t) = Σ_{l,k} O_{l,k}(t)|φ_{l,k}(r,t)|^2, dropping all off-diagonal components of the density matrix in the adiabatic basis. Under this diagonal approximation, the derived APT expression (Eq. 12) contains only band-diagonal matrix elements of ∇v_ion plus the N_e E(t) term. This is an intraband-only expression; it cannot describe interband recombination. Yet the paper's central mechanism—the delayed 'anomalous interband' emission that broadens the APT—is an off-diagonal (coherence) contribution. In the SBE validation (Appendix E), the emitted field is P(t) = Σ μ p (Eq. 19), with p the interband coherence; the occupation delay in Fig. 6(d) is a consequence of coherence-driven population dynamics (Eqs. 17-18), not an independent driver of emission timing. Thus the analytic derivation does not support the abstract's claim that APT temporal structure is 'characterized by occupation' or that the emission is 'occupation-driven.' The causal arrow in the title is not established by the paper's own equations; at best, occupation dynamics are correlated with the coherence-mediated delay.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript studies the temporal structure of attosecond pulse trains (APTs) from high-order harmonic generation (HHG) in bulk silicon using real-time time-dependent density functional theory (rt-TDDFT), supported by an analytical expression for APT emission decomposed in an adiabatic basis and a one-dimensional semiconductor Bloch equation (SBE) validation. The authors report that as laser intensity increases, the maximum emission peak shifts later and the FWHM of the low-pass-filtered APT varies nonmonotonically (Fig. 1(d,e)). They attribute this to symmetry-forbidden interband transitions (transition dipole moment between VBM and CBM+3 existing only at the Gamma point) causing delayed 'anomalous interband' recombination, and argue that electron occupation in different energy bands characterizes APT timing. The paper claims this establishes a fundamental time-domain constraint for solid-state attosecond sources.","tokens_in":16117,"tokens_out":6293,"duration_ms":57034,"significance":"Understanding the temporal profile of solid-state HHG is important for attosecond pulse synthesis, and the paper identifies a concrete mechanism—symmetry-restricted recombination at the Gamma point—that could explain intensity-dependent pulse broadening. The rt-TDDFT calculations reproduce the main features of experimental HHG spectra of silicon (Fig. 2), lending credibility to the numerical setup. The analytical APT expression is explicitly shown to reduce to a known result of Ref. [25] (Appendix C, Eq. (10)), which is a useful consistency check. However, the paper's central causal claim is not established by its own derivation: the diagonal-density approximation in Appendix C drops interband coherences, and the SBE emission is coherence-mediated, so occupation is at best correlated with, not a driver of, emission timing. No code or data availability statement is provided, so reproducibility rests primarily on the standard nature of the methods.","major_comments":[{"comment":"The derivation of the central occupation-resolved APT expression (Eq. (12)) assumes the time-dependent density can be written as n(r,t) = sum_{l,k} O_{l,k}(t) |phi_{l,k}(t)|^2, dropping all off-diagonal density-matrix elements in the adiabatic basis. Under this diagonal approximation the APT expression contains only band-diagonal matrix elements of the ionic potential gradient plus the N_e E(t) term; it is an intraband-only expression and cannot represent interband recombination emission. This matters because the delayed anomalous interband emission, which is the paper's explanation for the FWHM broadening, is a coherence (off-diagonal) contribution: in the SBE validation in Appendix E the emission is P(t) = sum mu p (Eq. (19)), with p the interband coherence, not the occupation. The analytic derivation therefore does not support the abstract's claim that APT temporal structure is 'characterized by occupation' or that emission is occupation-driven. At best, the equations show occupation and emission are correlated; the causal arrow in the title is not established by the paper's own formalism.","section":"Appendix C, Eqs. (11)-(13)"},{"comment":"The central nonmonotonic FWHM-versus-intensity trend and the peak-center shift are obtained by low-pass filtering all signals with a single cutoff energy chosen from the lowest-intensity signal (4.3 eV) and then tracking only the maximum peak. The manuscript reports no sensitivity analysis with respect to the filter cutoff, no uncertainty estimates (e.g., from the k-point sampling or the pulse envelope), and no discussion of how many peaks exist in the filtered signal at each intensity. Since the quantitative claim of the paper rests on this post-processing, the robustness of the nonmonotonic trend to the cutoff choice and to the peak-selection criterion must be demonstrated before the conclusion can be accepted.","section":"Section III.B and Appendix B (Fig. 1(d,e))"},{"comment":"The 1D SBE validation extends the equilibrium transition dipole moment between VBM and CBM+3, which is nonzero only at the Gamma point, to the entire 1D Brillouin zone using k.p perturbation theory, and employs a single dephasing time T2 = 484 as. The k-dependence of this TDM is precisely what controls whether the anomalous recombination is delayed; if the k.p extension artificially suppresses the off-Gamma TDM, the SBE reproduces the occupation delay by construction. An independent check is needed, for example by comparing SBE occupations directly with the 3D TDDFT band-resolved occupations at matched parameters, or by varying T2 and showing the delay is robust.","section":"Appendix E, Eqs. (14)-(19)"}],"minor_comments":[{"comment":"The sentence 'The mechanism underlying high harmonic generation (HHG) in gases has been well clarified, characterizing attosecond pulse trains (APT) in the time domain, significantly advances the synthesis of isolated attosecond pulse (IAP)' is grammatically broken and should be rewritten.","section":"Abstract"},{"comment":"The caption refers to panels (c), (d) and (e) where the figure appears to contain six panels (a)-(f); the main text also refers to 'Figure 3 (d), (e)' for chirp panels. Please correct the cross-references.","section":"Figure 3 caption"},{"comment":"The caption lists panels (a), (c), (e) for the low-intensity case but the text refers to (c), (e), (g); please unify the panel numbering.","section":"Figure 4 caption"},{"comment":"The phrase 'only can the anomalous temporal phenomenon be controlled by the first term of equation (8)' is unclear: the first term is an occupation-weighted potential integral, not a quantity that is 'controlled by' occupation. Please rephrase to state that the first term depends on occupation.","section":"Section IV.A"},{"comment":"Reference [66] is cited in Section III.B but the reference list contains only 60 entries; either add the missing reference or remove the citation.","section":"References"},{"comment":"There are numerous typos and grammatical issues, e.g., 'diffierent' in the Figure 3 caption, 'preform' in Appendix E, 'can could' in Section IV.A, and 'obstacles IAP separation' in the Abstract. A thorough language edit is needed.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The manuscript reads like a rough draft: the arXiv title differs from the full-text title, and the abstract differs from the arXiv abstract. More importantly, the analytical derivation's diagonal-density approximation removes the interband coherences that the SBE validation identifies as the emission source, so the occupation-driven causality claim is not supported. I recommend major revision with a request to either reframe the conclusion as a correlation or extend the derivation to include off-diagonal contributions. The TDDFT results themselves are plausible and the reproduction of experimental HHG spectra is a strength."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things. First, the paper reports a plausible mechanism for why attosecond pulse trains from bulk silicon broaden non-monotonically with laser intensity: certain interband recombination channels are symmetry-forbidden except at the Gamma point, so those carriers recombine half a cycle later, and the mixture of normal and 'anomalous' interband emission broadens the pulse. The TDDFT calculation looks standard, reproduces the experimental HHG spectra, and the 1D Bloch equations qualitatively reproduce the delay. Second, the paper's headline claim that 'APT can be characterized by the occupation of electrons' is not backed by its own equations. The analytic derivation in Appendix C expands n(r,t) in the adiabatic basis and keeps only diagonal terms (Eq. 11), which is an intraband-only expression. Interband emission, including the delayed recombination the paper cares about, is mediated by off-diagonal coherence—the SBE's own emission formula (Eq. 19) is exactly that. So the derivation does not show occupation drives emission timing; at best, occupations are correlated with the coherence dynamics. The stress-test note is right about this.\n\nWhat is new and good: the symmetry analysis of the transition dipoles (VBM-CBM+3 TDM only at Gamma) is a nice observation, and using it to explain the delayed 'anomalous' emission is a legitimate new contribution. The paper also does the right thing in acknowledging that its analytical expression reduces to ref [25]. The TDDFT setup is solid and matches experimental silicon HHG. The authors are transparent about approximations: no dephasing in TDDFT, a single T2 in the SBE, and a k.p extension of the TDM.\n\nWhere it is soft: (1) The causal claim is overstated. The title and abstract say 'occupation-driven' and 'fundamental constraint', but the paper demonstrates a correlation, not a derivation, and only in silicon. A single material with a specific TDM structure is not a fundamental constraint. (2) The central nonmonotonic FWHM result is produced by low-pass filtering with cutoff 4.3 eV (chosen as the upper limit of the lowest-intensity signal) and then selecting the maximum peak. There are no error bars, no sensitivity to the cutoff, and no discussion of how peak selection affects the trend. That is a fragile basis for the paper's main quantitative claim. (3) The 1D SBE uses a fixed dephasing T2=484 as and a k.p extrapolation of the Gamma-point TDM to the whole Brillouin zone; these are reasonable but unvalidated, so the mechanism is plausible rather than proven.\n\nBottom line: this paper will interest people working on solid-state HHG and attosecond source design, and it deserves a serious referee. But it should not pass as-is. I would recommend sending it to review, with expectations of major revision: the filter sensitivity needs to be shown, and the causal language needs to be toned down to match what the equations actually support.","headline":"A plausible symmetry-based mechanism for delayed interband emission in solid-state HHG, but the 'occupation-driven' claim is not supported by the paper's own analytic derivation, and the headline trend rests on fragile filtering.","tokens_in":16708,"tokens_out":2971,"would_cite":false,"duration_ms":31962,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["42.65.Ky"],"model":"deepseek-v4-flash","headline":"This paper argues that electron band occupation, not the laser field alone, sets the timing of attosecond pulses emitted by solids.","keywords":["attosecond pulse trains","solid-state high-harmonic generation","electron occupation dynamics","transition dipole moment symmetry","time-dependent density functional theory","semiconductor Bloch equations","bulk silicon","isolated attosecond pulses"],"falsifier":"An experimental or independent computational test would be to measure, or recalculate with fully k-resolved transition dipoles, the time-resolved harmonic emission of silicon from 0.1 to 1.0 TW/cm$^2$; if no emission component appears roughly half a laser cycle after the driving-field peak at high intensity, or if the filtered pulse width changes monotonically instead of nonmonotonically, the asynchronous-occupation mechanism is falsified.","tokens_in":15661,"feed_emoji":"⚡","tokens_out":9635,"duration_ms":89975,"temperature":0.7,"pith_summary":"This paper argues that the temporal shape of attosecond pulse trains emitted by a solid under intense mid-infrared light is set by the time-dependent occupation of electrons in each energy band, not just by the laser field. Using real-time time-dependent density functional theory on bulk silicon, the authors show that when the field is strong enough to populate higher conduction bands, a symmetry-forbidden interband channel recombines only at the $\\Gamma$ point half a laser cycle later, so part of the emission lags the field-synchronous part and the pulse broadens. They derive an occupation-resolved expression for the attosecond pulse train and reproduce the delayed occupation with a three-band semiconductor Bloch model. Because the delayed channel turns on and then loses to normal interband emission, the pulse width varies nonmonotonically with intensity. If correct, this identifies a time-domain constraint on solid-state attosecond sources that acts independently of material damage thresholds.","feed_headline":"Electron band occupancy controls attosecond pulse timing in solids","feed_subtitle":"A symmetry-forbidden transition delays part of the emission, so pulse width varies nonmonotonically with laser intensity.","key_machinery":"The central object is the time-dependent band occupation $O^{(e)}_{l,k}(t)$ in the adiabatic Kohn-Sham basis, because it appears directly inside the attosecond pulse train expression. The mechanism is carried by the transition dipole moment between the valence band maximum and the conduction band CBM+3, which is nonzero only at $\\Gamma$ under the laser polarization; recombination through it is therefore delayed by half a laser cycle. The authors support the TDDFT result with a one-dimensional three-band semiconductor Bloch equation, using one valence band and two conduction bands with a single dephasing time $T_2 = 484\\ \\mathrm{as}$, which reproduces the occupation delay.","core_discovery":"The central claim is that the time-resolved emission of a solid can be written as a sum of a field-synchronous term $N_e \\mathbf{E}(t)$ and an occupation-weighted ionic-potential term $\\int_\\Omega d^3r\\, \\sum_{l,k} O^{(e)}_{l,k}(t)\\, \\varphi^*_{l,k}(t)\\varphi_{l,k}(t)\\,\\nabla v_{\\mathrm{ion}}(\\mathbf{r})$, so the attosecond pulse train is governed by which Kohn-Sham orbitals are occupied at each instant. In bulk silicon at low intensity, the occupied states follow the field half-cycle and the pulse tracks the laser; at high intensity, occupation of higher bands lags the field. The lag is traced to the transition dipole moment between the valence band maximum and the band called CBM+3, which is nonzero only at the $\\Gamma$ point, so recombination through that channel can happen only after the electron is carried back to $\\Gamma$ half a cycle later. This asynchronous occupation delays part of the emission and broadens the pulse, and the competition among intraband, normal interband, and anomalous interband emission produces a nonmonotonic dependence of the pulse width on intensity.","pith_inferences":["The authors leave implicit that the predicted band-resolved occupation delay could be observed directly with time-resolved photoemission or transient absorption, giving a band-selective probe of emission timing.","Extending to nearby problems, strain or alloying that shifts the CBM+1/CBM+3 band crossing should tune the intensity at which the pulse width maximum occurs.","A practical corollary the authors do not draw is that a single intensity value cannot define the temporal quality of a solid-state attosecond source; pulse width, center, and yield must be mapped together.","In a neighboring geometry, the same mechanism should appear in two-dimensional semiconductors with multiple conduction bands and symmetry-forbidden dipoles, where the delay may be tunable by stacking or twist angle."],"forward_implications":["Isolated attosecond pulse extraction from solid-state high-harmonic generation is harder than the intraband-only picture suggests, because the delayed anomalous channel adds emission after the main burst.","The attosecond pulse width is not a simple saturation curve: it first grows as the anomalous interband channel turns on, then shrinks as normal interband emission dominates.","The same mechanism predicts that time-domain attosecond pulse characterization must report intensity together with pulse center and width, since all three change with field strength.","Materials with no symmetry-forbidden high-lying conduction-band recombination channels should preserve shorter pulses at high intensity.","The occupation-resolved derivation gives a quantitative criterion for searching laser parameters that minimize the delay when generating isolated attosecond pulses."],"supporting_citations":[{"why":"Supplies the silicon prototype and the picture of strong-field excitation at the valence-band maximum.","marker":"[10]"},{"why":"Provides the theoretical separation of intraband and interband mechanisms that the paper is built on.","marker":"[20]"},{"why":"Supplies the generalized electron-hole recollision picture used for the normal interband channel.","marker":"[21]"},{"why":"Gives the analytic high-harmonic generation expression in the adiabatic basis that the occupation-resolved derivation extends.","marker":"[25]"},{"why":"Provides the reciprocal-space trajectory and pre-acceleration picture used to describe field-driven excitation.","marker":"[32]"},{"why":"Provides the real-time time-dependent density functional theory method used for all ab initio pulse and spectrum calculations.","marker":"[35]"},{"why":"Supplies the experimental silicon high-harmonic generation results used to validate the computational setup.","marker":"[40]"},{"why":"Establishes the central role of transition dipole amplitude and phase in interband emission, motivating the transition dipole moment analysis.","marker":"[42]"},{"why":"Used to interpret the intensity-dependent shift of the emission center through dipole-phase interferometry.","marker":"[43]"}],"fun_headline_variants":["Solid-state attosecond pulse width hinges on band occupation lag","Symmetry-forbidden transition delays emission, widens attosecond pulses","Asynchronous band occupation sets fundamental limit on attosecond pulses","Nonmonotonic pulse width in solids traced to delayed band emission","Occupation-driven asynchrony governs solid-state attosecond pulse shape"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The explanation rests on the assumption that the transition from the top valence band to the conduction band labeled CBM+3 can really occur only at the $\\Gamma$ point in reciprocal space; if that restriction is wrong, the delayed recombination and the predicted pulse broadening would not occur.","fun_headline_variants_meta":{"raw":{"variants":["Solid-state attosecond pulse width hinges on band occupation lag","Symmetry-forbidden transition delays emission, widens attosecond pulses","Asynchronous band occupation sets fundamental limit on attosecond pulses","Nonmonotonic pulse width in solids traced to delayed band emission","Occupation-driven asynchrony governs solid-state attosecond pulse shape"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000279,"raw_usage":{"total_tokens":1637,"prompt_tokens":906,"completion_tokens":731,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":522,"completion_tokens_details":{"reasoning_tokens":645}},"tokens_in":522,"tokens_out":731,"duration_ms":7620,"temperature":1.0,"reasoning_tokens":645,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T14:38:40.604641+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An experimental or independent computational test would be to measure, or recalculate with fully k-resolved transition dipoles, the time-resolved harmonic emission of silicon from 0.1 to 1.0 TW/cm$^2$; if no emission component appears roughly half a laser cycle after the driving-field peak at high intensity, or if the filtered pulse width changes monotonically instead of nonmonotonically, the asynchronous-occupation mechanism is falsified.","supporting_citations":[{"cited_title":"Ghimire, A","cited_arxiv_id":null,"evidence_quote":"Supplies the silicon prototype and the picture of strong-field excitation at the valence-band maximum."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the theoretical separation of intraband and interband mechanisms that the paper is built on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the generalized electron-hole recollision picture used for the normal interband channel."},{"cited_title":"Ikemachi, Y","cited_arxiv_id":null,"evidence_quote":"Gives the analytic high-harmonic generation expression in the adiabatic basis that the occupation-resolved derivation extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the reciprocal-space trajectory and pre-acceleration picture used to describe field-driven excitation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the real-time time-dependent density functional theory method used for all ab initio pulse and spectrum calculations."},{"cited_title":"Chini, K","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental silicon high-harmonic generation results used to validate the computational setup."},{"cited_title":"Suthar, F","cited_arxiv_id":null,"evidence_quote":"Establishes the central role of transition dipole amplitude and phase in interband emission, motivating the transition dipole moment analysis."},{"cited_title":"Floss, C","cited_arxiv_id":null,"evidence_quote":"Used to interpret the intensity-dependent shift of the emission center through dipole-phase interferometry."}],"review_version":1}