{"id":"e9d876ee-5ed0-41f4-84f7-13c024f47802","arxiv_id":"2509.04869","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"Opposite sub-cycle correlation dynamics appear in noble-gas versus alkaline-earth atoms, and the correlation-driven evolution of natural orbitals, rather than occupation numbers, dominates the HHG spectrum.","lead":"This paper uses the one-particle density matrix to describe high-harmonic generation in multi-electron atoms, expressing the emission as a sum over natural-orbital 'antennas' with time-dependent weights. Simulations of He, Ne, Be, and Mg show that electron-correlation measures rise during the laser pulse in noble gases but fall in alkaline-earth atoms.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"MCTDHF convergence is asserted but not documented; sign reversal and occupation-number conclusions could be truncation artifacts.","rationale":"The paper's formal framework—expressing the dipole acceleration as a sum over natural-orbital 'antennas' (Eq. 11) and separating occupation-number dynamics from orbital dynamics via the collision operator—is mathematically sound and clearly presented. The central claim that occupation-number variation has minor impact on HHG is directly supported by the near-coincidence of FFE and SFE spectra with full MCTDHF, while the large TDHF/MCTDHF differences indicate that correlations matter. However, every quantitative trend, especially the surprising sign reversal of δS1 between noble gases and alkaline-earth atoms, rests on the convergence of the MCTDHF propagations. The paper asserts convergence but provides no convergence plots, and the active spaces—particularly Be with M=5—are small enough that field-driven correlation dynamics could be under-resolved. This is exactly the reader's weakest assumption, and it is load-bearing because if the sign reversal or the FFE/SFE near-coincidence changes with active-space size, the central physical conclusion collapses. The concern does not invalidate the paper; it makes the current evidence conditional, consistent with the reader's CONDITIONAL verdict. No other concern appears equally load-bearing: the derivation of Eqs. 18 and 20 is standard, the FFE/SFE decomposition is a valid numerical diagnostic, and the interpretation via the collision operator is plausible. Thus the reader's verdict should stand unchanged.","tokens_in":22471,"tokens_out":6105,"duration_ms":64737,"concrete_test":"Repeat the Be and Mg MCTDHF calculations with increased active spaces (e.g., Be M=7 or 9 adding 3s/3p/3d orbitals; Mg M=15 or 17) and the same pulse parameters, and recompute δS1(t), δ||Δ12(t)||, and the HHG spectra. Check: (1) Does δS1(t) remain negative for Be/Mg and positive for He/Ne? (2) Do the FFE and SFE spectra remain within the stated near-coincidence with MCTDHF (e.g., second-plateau yield changes by less than ~50%)? (3) Do the plateau heights and cutoff shapes stabilize with M? If the sign of δS1 reverses or FFE/MCTDHF differences grow substantially, the central conclusion is a truncation artifact. Also provide explicit convergence plots versus M, lmax, and grid spacing.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative results—the sign reversal of δS1(t) and δ||Δ12(t)|| between noble gases and alkaline-earth atoms (Fig. 4) and the claim that occupation-number dynamics have minor impact on the high-energy plateau (Sec. V.C)—are derived entirely from MCTDHF wavefunctions. Appendix A states 'Convergence of the results for these parameters has been verified' but presents no convergence data. The active spaces are modest (Be: M=5; Ne: M=9; Mg: M=13), and for Be M=5 contains only 1s, 2s, 2p−1, 2p0, 2p+1 orbitals, which captures static 2s–2p correlation but may not be sufficient for field-driven dynamical correlation. If the sign of δS1(t) or the plateau yields change with M, the claimed distinction between noble-gas and alkaline-earth correlation dynamics would be a truncation artifact rather than a physical effect. The paper's 'unambiguous' attribution to natural-orbital dynamics is also slightly undercut by the MF model (Eq. 37), where the initial correlated natural orbitals are not eigenstates of the mean-field propagator; the resulting 'beats' are acknowledged but the removal procedure is not specified. The cleaner FFE/SFE comparison (Eqs. 38–39) does support the minor role of occupation-number variation, but its validity also hinges on the same MCTDHF convergence.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a reduced-density-matrix (1RDM) description of high-harmonic generation in multi-electron atoms. It shows that the dipole acceleration can be written exactly as a sum over natural-orbital 'antennas' weighted by the corresponding occupation numbers (Eq. 11), and derives equations of motion for the natural-orbital occupations (Eq. 18) and orbitals (Eq. 20) that separate mean-field and two-particle-correlation (collision-operator) contributions. The framework is then applied to MCTDHF simulations of He, Ne, Be, and Mg driven by two-cycle pulses with similar Keldysh parameters. The authors report a qualitative distinction: noble gases show a sub-cycle increase of the von Neumann entropy and of the two-particle cumulant norm, whereas alkaline-earth atoms show a decrease. They further compare MCTDHF spectra with TDHF, TDDFT, and controlled models with frozen (FFE) or sharp (SFE) Fermi edges, concluding that the high-energy plateau is strongly influenced by correlations acting on the natural-orbital dynamics, while the time variation of the occupation numbers has only a minor spectral impact.","tokens_in":22807,"tokens_out":5256,"duration_ms":53590,"significance":"The exact decomposition (Eq. 11) and the equations of motion (Eqs. 18 and 20) provide a conceptually clean and potentially widely applicable language for correlation effects in strong-field processes. The FFE/SFE comparison is a principled way to isolate the spectral role of occupation-number dynamics from that of natural-orbital evolution. The predicted sign reversal of the sub-cycle correlation variation between noble gases and alkaline-earth atoms is an interesting, falsifiable claim. If the numerical convergence is properly established, this would be a valuable contribution to attosecond and strong-field correlation physics. The paper is also notable for explicitly placing its results on an exact reduced-density-matrix footing rather than on a phenomenological model.","major_comments":[{"comment":"The central quantitative conclusions—the sign of δS1(t) and δ||Δ12(t)|| (Fig. 4), the relative plateau heights (Figs. 9–12), and the FFE/SFE agreement—are derived entirely from MCTDHF wavefunctions. Appendix A states 'Convergence of the results for these parameters has been verified' but presents no convergence data. The active spaces are modest (Be: M=5; Ne: M=9; Mg: M=13), and for Be the space contains only 1s, 2s, and 2p orbitals, which may not capture all field-driven dynamical correlation in the continuum or high-l partial waves. If the sign of δS1(t) or the plateau yields change with M or grid size, the noble-gas versus alkaline-earth distinction could be a truncation artifact rather than a physical effect. Please provide convergence tests (e.g., δS1 and HHG spectra versus M and against a larger l_max) or explicitly state the estimated uncertainty.","section":"Sec. V.A and Appendix A"},{"comment":"The claim that the HHG spectra are 'strongly and unambiguously influenced' by two-particle correlations through their effect on the natural orbitals relies in part on the MF model. However, with Eq. (37) the initial natural orbitals are not eigenstates of the mean-field propagator; the text acknowledges additional high-frequency 'beats' that 'should be removed' but does not specify the removal procedure. If a spectral filter or smoothing is applied, its effect on the high-energy plateau must be documented. The FFE/SFE comparison robustly supports the minor role of occupation-number variation, but the MF comparison is the main evidence for the strong role of correlations in the orbital evolution. Without a defined beat-removal procedure, the word 'unambiguously' is too strong.","section":"Sec. V.C, Eq. (37)"}],"minor_comments":[{"comment":"Typo: 'cummulant' should be 'cumulant' in the caption.","section":"Fig. 4 caption"},{"comment":"There is a missing closing bracket in the displayed equation; it should read i∂tηj(t) = ⟨ηj(t)| Tr2[Ŵ12, Δ12(t)] |ηj(t)⟩.","section":"Eq. (18)"},{"comment":"MF results are shown for He, Ne, and Be but not for Mg. The text states this is due to numerical-convergence issues; it would be helpful to state this explicitly in the text and to comment on whether the absence of Mg affects the generality of the comparison.","section":"Sec. V.C, Eqs. (37)–(39)"},{"comment":"The recombination probability Prec is used in Eq. (28) but not defined. Please define it explicitly.","section":"Eq. (28)"}],"recommendation":"major_revision","confidential_remarks":"The methodological references include several from the authors' own groups, but these are established formalisms (MCTDHF, TD-CASSCF, 2RDM methods) rather than the target of the paper; I see no circularity. The exact decomposition and the controlled FFE/SFE comparisons are strong points. The main risk to the central empirical claim is the lack of convergence documentation; if the requested convergence tests show instability of the sign reversal or of the plateau-height conclusions, the paper's main message would require substantial qualification. This is addressable within the scope of a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does two things. First, it gives a clean RDM framework for HHG: the dipole acceleration is exactly decomposed into natural-orbital antennas, and the equations of motion for occupations and orbitals are spelled out. That's standard RDM theory, but clearly presented and directly useful. Second, it applies the framework to He, Ne, Be, and Mg and finds a surprising trend: sub-cycle entropy and cumulant variation increase in noble gases, decrease in alkaline earths. The FFE/SFE comparison is a good control; it shows that occupation-number dynamics has little effect on HHG, while the orbital dynamics carries the correlation effect.\n\nThe main soft spot is the one the reader flagged: convergence of the MCTDHF calculations is asserted but not documented. The sign reversal is the central physical claim; if it changes with active space or grid, the narrative collapses. I'm not convinced it's an artifact—the entropy and cumulant norms agree, and the collision-operator maps show plausible spatial patterns—but without a convergence plot I can't rule it out. The word \"unambiguously\" in Sec. V.C is also a bit strong; the MF reference has beats that have to be removed, and the removal procedure isn't specified. That's minor, but it means \"unambiguous\" should be \"consistent with.\"\n\nThe formal derivations check out. Eq. 11 is exact within the Larmor formula; Eqs. 18 and 20 are correct. The citation pattern is fine; the self-citations are to established methods, not to questionable results.\n\nOverall, a solid paper with an interesting but not fully documented numerical result. It deserves a serious referee, but that referee should demand convergence data before the sign-reversal claim is taken as established.","headline":"Solid RDM-based decomposition of HHG with an interesting but convergence-dependent claim about correlation sign reversal.","tokens_in":23305,"tokens_out":1623,"would_cite":true,"duration_ms":17110,"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":"High-harmonic emission from multi-electron atoms decomposes exactly into radiation from natural-orbital antennas, and the correlation-driven re-shaping of those antennas—not their occupation numbers—sets the high-energy plateau.","keywords":["high-harmonic generation","reduced density matrix","natural orbitals","electron correlation","sub-cycle dynamics","strong-field ionization","noble gas atoms","alkaline earth atoms"],"falsifier":"A systematic convergence study varying both the number of active orbitals and the radial grid for one atom of each pair (e.g., Ne and Mg) and checking whether the sign of the entropy change and the full-propagation versus frozen-Fermi-edge plateau difference remain unchanged; alternatively, an experiment that independently measures the sub-cycle ionization-burst enhancement in He and suppression in Be via attosecond transient absorption and compares with the predicted correlation contribution to the ionization rate.","tokens_in":22410,"feed_emoji":"⚛","tokens_out":5706,"duration_ms":54061,"temperature":0.7,"pith_summary":"This paper introduces a way to write the high-harmonic generation (HHG) signal of a correlated multi-electron atom as a sum of independent radiating oscillators, one per natural orbital of the one-particle reduced density matrix. The exact decomposition lets the authors ask which part of the correlation matters: the time evolution of the natural orbitals themselves, or the time evolution of their occupation numbers. Using fully converged many-body propagations for helium, neon, beryllium, and magnesium, they find that noble gases and alkaline earth atoms respond oppositely at the sub-cycle level: correlation measures increase near the field extrema for He and Ne, and decrease for Be and Mg. They trace this sign flip to the shape and parity of the natural orbitals involved in the population transfer, and show that the high-energy plateau is strongly shaped by correlation effects on the orbitals, while occupation-number reshuffling plays only a minor role.","feed_headline":"Correlation dynamics reverse sign between noble gases and alkaline earth atoms","feed_subtitle":"Natural-orbital decomposition traces the HHG plateau to orbital motion, not occupation change.","key_machinery":"The one-particle reduced density matrix and its natural-orbital decomposition, which turns the HHG radiation amplitude into a superposition of independent 'atomic antennas'. The two-particle cumulant (collision operator) enters as the sole source of occupation-number changes and as one of the drivers of orbital motion; comparing full propagation with models that freeze or approximate the Fermi edge isolates which part of the correlation controls the spectrum.","core_discovery":"For an arbitrary N-electron system driven by a strong laser, the Larmor dipole acceleration can be written exactly as a superposition over the natural orbitals of the one-particle reduced density matrix: each term is an independent 'atomic antenna' weighted by a time-dependent occupation number. The equation of motion shows that occupation-number changes are driven exclusively by the two-particle cumulant, while the orbital evolution mixes strong-field, mean-field, and correlation contributions. Comparing noble gases with alkaline earth atoms at matched Keldysh parameters, the paper finds that the sub-cycle variation of both the von Neumann entropy and the norm of the two-particle cumulant i","pith_inferences":["The same natural-orbital antenna decomposition should apply to other one-body observables (dipole moment, momentum density) in molecules and solids, giving a general sub-cycle diagnostic for correlation-driven emission.","For atoms with shells of opposite parity near the Fermi edge, such as xenon, the parity argument suggests the correlation flux may reverse sign within a single driving cycle, offering a sub-cycle interpretation of the giant-dipole-resonance enhancement.","The sign reversal between noble gases and alkaline earths predicts that, in HHG spectroscopy of atoms with tunable ground-state correlation, the sub-cycle modulation of plateau yield should flip sign as the ground-state correlation across the active pair changes.","An experiment that independently measures the sub-cycle ionization-burst enhancement in He and suppression in Be via attosecond transient absorption could test the predicted correlation contribution to the ionization rate without relying on the full HHG spectrum."],"forward_implications":["The radiation amplitude of a correlated atom can be computed exactly from natural orbitals and occupation numbers alone, eliminating the need for the full two-particle wavefunction for this observable.","Because occupation-number evolution is driven solely by two-particle correlations, the von Neumann entropy of the one-particle reduced density matrix is a sub-cycle, in-principle-measurable clock of dynamical correlations.","In noble gases, correlations enhance tunnel ionization and the high-energy plateau; in alkaline earth atoms, correlations suppress the ionization burst, consistently with the reversed direction of collision-operator flux between HONO and LUNO.","The frozen-Fermi-edge model reproduces most of the HHG spectrum, so the orbital dynamics, not the occupation-number dynamics, carries the correlation fingerprint into the plateau."],"supporting_citations":[{"why":"Supplies the 'atomic antenna' picture of HHG that this paper extends to a correlated many-body description.","marker":"[19]"},{"why":"One of the foundational formulations of the multi-configuration time-dependent Hartree-Fock method used to propagate the correlated wavefunction.","marker":"[43]"},{"why":"Provides the time-dependent complete-active-space self-consistent-field implementation and the ionization-time-shift analysis that the present numerical propagations build on.","marker":"[48]"},{"why":"Reports the correlated secondary plateau in helium that the present decomposition aims to explain and quantify via orbital dynamics.","marker":"[30]"},{"why":"Attributes the extended plateau in beryllium to recollision-induced two-electron correlation, the alkaline-earth comparison case for the sign reversal.","marker":"[32]"},{"why":"Defines the two-particle cumulant used to separate mean-field from correlation contributions.","marker":"[52]"}],"fun_headline_variants":["Sub-cycle correlation flips sign in noble vs alkaline earth atoms","Natural orbitals trace HHG: sign reversal in correlation dynamics","Noble gases vs alkaline earth: correlation changes direction","Ultrafast correlation dynamics reverse between atomic families","HHG: correlation sign switch traced to natural orbital motion"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The quantitative conclusions rest on the assertion that the many-body wavefunction propagations are converged at the stated active-space sizes and radial grids; if the observed sign of the correlation measures or the plateau-height differences shift when these parameters are increased, the central claims would not survive.","fun_headline_variants_meta":{"raw":{"variants":["Sub-cycle correlation flips sign in noble vs alkaline earth atoms","Natural orbitals trace HHG: sign reversal in correlation dynamics","Noble gases vs alkaline earth: correlation changes direction","Ultrafast correlation dynamics reverse between atomic families","HHG: correlation sign switch traced to natural orbital motion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000372,"raw_usage":{"total_tokens":1789,"prompt_tokens":672,"completion_tokens":1117,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":416,"completion_tokens_details":{"reasoning_tokens":1038}},"tokens_in":416,"tokens_out":1117,"duration_ms":8263,"temperature":1.0,"reasoning_tokens":1038,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T05:48:27.583666+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A systematic convergence study varying both the number of active orbitals and the radial grid for one atom of each pair (e.g., Ne and Mg) and checking whether the sign of the entropy change and the full-propagation versus frozen-Fermi-edge plateau difference remain unchanged; alternatively, an experiment that independently measures the sub-cycle ionization-burst enhancement in He and suppression in Be via attosecond transient absorption and compares with the predicted correlation contribution to the ionization rate.","supporting_citations":[{"cited_title":"Floss, C","cited_arxiv_id":null,"evidence_quote":"One of the foundational formulations of the multi-configuration time-dependent Hartree-Fock method used to propagate the correlated wavefunction."},{"cited_title":"Kato and H","cited_arxiv_id":null,"evidence_quote":"Provides the time-dependent complete-active-space self-consistent-field implementation and the ionization-time-shift analysis that the present numerical propagations build on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the correlated secondary plateau in helium that the present decomposition aims to explain and quantify via orbital dynamics."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Attributes the extended plateau in beryllium to recollision-induced two-electron correlation, the alkaline-earth comparison case for the sign reversal."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the two-particle cumulant used to separate mean-field from correlation contributions."}],"review_version":1}