{"id":"3db3b87f-68f3-4b05-b95b-482a311eeffe","arxiv_id":"2412.19612","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A TDDFT simulation study showing that circular dichroism in attosecond transient absorption between the first and second conduction bands can quantitatively map valley polarization dynamics in h-BN and MoS2.","lead":"This paper proposes an all-optical method to read out valley polarization in 2D materials with sub-femtosecond time resolution. It measures the difference in absorption of left and right circularly polarized attosecond probe pulses, which tracks the electron population imbalance between valleys, and demonstrates the idea in simulations of h-BN and MoS2.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The advertised 250 as time resolution is not established: with Eq. (3)'s 1 fs (0.5 fs FWHM) probe, Fig. 2(e) only shows a 250 as delay shift, not the ability to resolve VP changes on that timescale.","rationale":"Good-faith reading: the core mapping Delta mu proportional to Delta N_CB1 is supported by TDDFT for two materials (Fig. 2(f), Fig. 4(b)), and the matrix-element checks in Fig. 2(a,b) address valley selectivity. Those are real evidence. But the central advertised advance is an unprecedented 250 as time resolution. The reasoning in Fig. 2(e) is not a resolution criterion; it only shows the delay stage can be stepped by 250 as. Given the 1 fs probe envelope, the retrieved signal is inherently time-integrated over about 0.5 fs, so sub-250 as resolution would require an explicit deconvolution or a sharper effective gate, neither of which is described. The paper's own protocol sweeps tau_p in 250 as steps, but sampling rate is not temporal resolution. This is an internal gap in the argument, not a disagreement with consensus. It does not invalidate the underlying proposal, but it means the strongest advertised claim should be conditional on a proper resolution test. The reader's identified weakest assumption (matrix-element asymmetry) is a different, though related, structural risk; I judge the unproven resolution claim to be more immediately load-bearing because it is a headline quantitative assertion and is checkable from the equations in the main text. Since the reader already issued CONDITIONAL and my concern strengthens the conditions, no change to the verdict is needed.","tokens_in":8718,"tokens_out":8697,"duration_ms":98998,"concrete_test":"Run the same retrieval protocol on a controlled step-function VP: initialize Delta N_CB1(t) as a Heaviside step at a known t0 (e.g., by artificially injecting a narrow-band CB1 population in the TDDFT run, or by using a delta-like switching perturbation), then apply the Eq. (3) probe and compute Delta mu(tau_p) at 7.5 eV. Measure the 10-90% rise time of the retrieved signal. If the rise time is not below 250 as (or, more realistically, is close to the 0.5 fs probe FWHM), the 250 as resolution claim fails. A complementary check: place two step switches 250 as apart and test whether the retrieved trace resolves two distinct events.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (3) defines the probe envelope as sin^2(pi(t-tau_p)/T_probe) with T_probe = 1 fs, i.e., a pulse of about 0.5 fs FWHM. The paper's evidence for 250 as resolution (Fig. 2(e)) is that responses for tau_p = 7.5 fs and 7.75 fs are clearly separated. That is a trivial consequence of delaying the probe by 250 as; it does not show that two VP features separated by 250 as can be distinguished. Transient absorption with a weak probe is, to first order, a convolution of the pump-induced change with the probe envelope, so the instrument response is at best the probe duration (~0.5 fs) unless a deconvolution or nonlinear correlation scheme is provided; none is. The statement that 'the large energy separation ... supports subfemtosecond time resolution' conflates the CB1-CB2 transition energy (7.5 eV) with temporal resolution: the relevant quantity is the duration/bandwidth of the probe and the time-integrated nature of the measured absorption. The quantitative retrieval may still be valid, but the headline 'unprecedented ... 250 as resolution' is not supported by the evidence presented.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an all-optical scheme for retrieving valley polarization (VP) dynamics in two-dimensional hexagonal materials using attosecond circular dichroism (CD) in transient absorption. The central idea, expressed in Eq. (1), is that the difference in absorption coefficients for left- and right-circularly polarized probe pulses, Δμ = μ⁺ − μ⁻, is proportional to the CB1 population imbalance between the K1 and K2 valleys, ΔN_CB1. The scheme is demonstrated numerically with TDDFT (Octopus) for h-BN and MoS₂: a circularly polarized pump creates VP, and a 1-fs sin²-envelope probe with photon energy matching the CB1–CB2 gap produces a CD signal whose magnitude follows the VP dynamics. The manuscript reports retrieval of transient VP switching events with a claimed 250-as time resolution.","tokens_in":8985,"tokens_out":3071,"duration_ms":34984,"significance":"If the central mapping of Eq. (1) is valid, the scheme offers a direct, quantitative, all-optical readout of valley polarization on sub-femtosecond timescales, complementing or surpassing existing photoluminescence, Kerr rotation, and high-harmonic spectroscopy approaches. The paper's strengths include the explicit linear relationship between an observable and the valley population imbalance, the use of ab initio TDDFT simulations for two prototypical materials, a clear calibration procedure via one auxiliary ARPES measurement, and a visually convincing demonstration that the retrieved signal tracks the simulated VP in both h-BN and MoS₂. These features make the core proposal attractive. However, the manuscript's headline claim of 250-as temporal resolution is not supported by the evidence presented, and the derivation of the central proportionality is relegated to the Supplemental Material.","major_comments":[{"comment":"The claim of '250 as time resolution' is not established. The curves in Fig. 2(e) for probe delays of 7.5 fs and 7.75 fs show that two probe pulses with a 250-as delay difference produce different time-integrated CB2 populations, which is a trivial consequence of the delay. This does not demonstrate the ability to resolve two VP features that are separated by 250 as in time. For a weak probe, transient absorption is to first approximation a convolution of the pump-induced change with the probe envelope, so the instrument response is limited by the probe duration, here T_probe = 1 fs (FWHM ≈ 0.5 fs), unless a deconvolution or nonlinear correlation scheme is provided. The statement that the large CB1–CB2 energy separation 'supports subfemtosecond time resolution' conflates spectral separation with temporal resolution. I recommend either providing a two-feature test that actually resolves two VP events separated by 250 as, or replacing '250 as time resolution' in the abstract and conclusions with claims about a 250-as delay sampling step and explicitly stating that the temporal resolution is limited by the probe duration.","section":"Fig. 2(e), Fig. 3(e)–(f), and the concluding paragraph"},{"comment":"The central proportionality Δμ ∝ N_CB1^K1 − N_CB1^K2 is stated in Eq. (1) but its derivation is deferred to the Supplemental Material. As submitted, the main text does not contain the assumptions and steps leading to this relation, making the central claim difficult to verify. In particular, the derivation must establish that the CB1–CB2 transition is perfectly valley-selective with equal matrix-element magnitudes at K1 and K2 for the two helicities, and that any helicity-dependent background or off-resonant contribution is absent or negligible. Please include the derivation in the manuscript (or ensure the Supplemental Material is part of the submission) and state explicitly the conditions under which Eq. (1) holds.","section":"Eq. (1) and Sec. I of the Supplemental Material"},{"comment":"The quantitative retrieval relies on the proportionality constant in Eq. (1) being time-independent and identical for the pump-driven dynamics. The linear fits in Fig. 2(f) and Fig. 4(b) are performed for VP states prepared with different pump strengths, which is not the same as verifying that the proportionality holds dynamically during the switching process, where the pump pulse is present, band populations change, and field dressing may modify transition matrix elements. The manuscript should address whether the calibration constant obtained in a static or slowly varying VP condition remains valid during the sub-femtosecond switching transients, or alternatively restrict the quantitative retrieval claims to the post-switching regime.","section":"Sec. IV and Fig. 2(f), Fig. 4(b)"}],"minor_comments":[{"comment":"The phrase 'unprecedented quantitative retrieval ... with a time resolution of 250 as' should be softened unless the temporal-resolution concern raised above is resolved; at present it overstates what the simulations show.","section":"Abstract and conclusion"},{"comment":"The label 'Populaion' in Fig. 3(b) contains a typo; it should read 'Population'.","section":"Fig. 3(b) label"},{"comment":"The units of the vector potential in Fig. 4(a) are not specified in the caption; please add the appropriate units or state that atomic units are used.","section":"Fig. 4(a) caption"},{"comment":"The terms N_CB1^K1 and N_CB1^K2 are defined in the text, but the superscripts are used inconsistently in places (e.g., 'K1(K2)' and 'K1' vs 'K1'); please harmonize the notation.","section":"Notation near Eq. (1)"}],"recommendation":"major_revision","confidential_remarks":"The paper has a potentially interesting and useful idea, and the numerical demonstrations are suggestive. My main concern is the temporal-resolution claim, which is currently unsupported and is likely to draw criticism in the community. This is fixable with a more careful statement of what the simulations show and, ideally, an additional numerical test with two well-separated short-duration VP features. The derivation of Eq. (1) should also be moved into the main text or the SM made available. If these points are addressed, the manuscript could be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi — wanted your read on 2412.19612. The core idea is genuinely new and worth taking seriously: instead of using the valence-to-CB1 transition, they use circular dichroism between CB1 and higher empty bands, so the measured absorption difference is a direct linear proxy for the valley population imbalance in CB1. That is a real step past PL, Kerr, SHG, and harmonic schemes, which give indirect or time-averaged signals. The numerical demonstration is also legitimate: TDDFT simulations in h-BN and MoS2, matrix-element checks, a linear fit between Δμ and ΔN, and retrieval of a three-pulse VP switching sequence. This is not a curve-fit of the target; the same simulation produces both the populations and the absorption.\n\nSoft spots, in order of seriousness.\n\nFirst, the 250 as resolution claim is not supported. The probe is a 1 fs sin^2 envelope, about 0.5 fs FWHM. A weak-probe transient absorption is, to first order, the pump-induced change convolved with the probe envelope; you cannot resolve features separated by 250 as with a 500 as probe unless you do deconvolution or a nonlinear correlation scheme, and none is provided. Figure 2(e) merely shows that shifting the probe delay by 250 as shifts the response—that is trivial. The argument that 'the large energy separation supports subfemtosecond time resolution' is a category error: transition energy is not temporal resolution.\n\nSecond, Eq. (1), the central relation, is deferred to the Supplemental Material. The main text should contain the derivation or at least the assumptions behind the proportionality—specifically, equal CB1->CB2 matrix elements at the two valleys and no helicity-dependent background. Fig. 2(f) checks this numerically for h-BN, but a structural claim needs the analytic statement.\n\nThird, the 'all-optical' label is slightly oversold: extracting the absolute VP requires one auxiliary ARPES measurement to calibrate the proportionality constant. That is a reasonable one-time calibration, not a fatal flaw, but it should be stated honestly.\n\nMinor: the valley integration radius R and the proportionality constant are free parameters; there are no error bars on the retrieval. For a theory proposal this is acceptable, but it limits how strongly one can claim 'unprecedented quantitative retrieval.' The citation pattern looks solid; they engage the recent HHG and valleytronics literature, including the HHS-concerns reference.\n\nOverall: the central physics is plausible and the numerical package is coherent. The paper deserves a serious referee; the referee should push for the SM derivation in the main text, a real resolution criterion, and a softened or better-evidenced 250 as claim. I'd bring it to reading group—it will spark a useful argument about what 'time resolution' means in attosecond pump-probe.","headline":"A genuinely new all-optical valley-polarization readout idea, backed by clean TDDFT numerics, but the 250 as resolution claim is overstated and the central derivation lives in the SM.","tokens_in":9527,"tokens_out":2686,"would_cite":false,"duration_ms":25999,"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":"The paper claims that the difference in absorption of two circularly polarized attosecond probe pulses is proportional to the valley population imbalance in a 2D hexagonal material, enabling quantitative retrieval of valley polarization…","keywords":["valley polarization","attosecond transient absorption","circular dichroism","two-dimensional hexagonal materials","ultrafast valley dynamics","time-dependent density functional theory","h-BN","MoS2"],"falsifier":"Perform the same pump-probe simulation on a material with slightly broken K/K' symmetry, such as a strained h-BN monolayer, and check whether the proportionality constant in Eq. (1) stays time-independent when the directly counted $\\Delta N_{\\mathrm{CB1}}$ is compared with $\\Delta\\mu$; a drift would show that the linear mapping is not universal.","tokens_in":8528,"feed_emoji":"⚛️","tokens_out":7748,"duration_ms":67163,"temperature":0.7,"pith_summary":"This paper proposes an all-optical scheme that measures valley polarization in two-dimensional hexagonal materials with sub-femtosecond time resolution. The central idea is that two circularly polarized attosecond probe pulses with opposite helicities absorb differently by an amount proportional to the valley population imbalance in the first conduction band. That proportionality turns transient absorption spectroscopy into a direct, quantitative readout of valley polarization dynamics, including the brief switching events between valleys. The authors demonstrate the scheme in time-dependent density functional theory (TDDFT) simulations for h-BN and MoS2, retrieving valley polarization dynamics with a 250 attosecond time step. If the claim holds, it gives quantum materials researchers an all-optical, quantitative probe of valley dynamics at timescales that were previously out of reach.","feed_headline":"Circular dichroism reads valley polarization in 250 attoseconds","feed_subtitle":"Opposite-helicity attosecond probes turn valley population imbalance into an absorption difference for quantitative readout.","key_machinery":"The load-bearing identity is Eq. (1), $\\Delta\\mu = \\mu_+ - \\mu_- \\propto N^{K_1}_{\\mathrm{CB1}} - N^{K_2}_{\\mathrm{CB1}}$, which maps the valley population imbalance onto a difference in absorption coefficients for the two probe helicities. The mechanism that makes the identity hold is the valley-selective circular dichroism transition from the first conduction band (CB1) to a higher conduction band (CB2), whose momentum matrix elements for the two helicities are equal in magnitude but couple to opposite valleys. This transition converts the population imbalance into an optical asymmetry, while the large CB1–CB2 energy separation supports the attosecond probe duration needed for sub-femtosecond time resolution.","core_discovery":"The central claim, expressed as Eq. (1), is that the circular dichroism $\\Delta\\mu = \\mu_+ - \\mu_-$ of two opposite-helicity probe pulses is proportional to the instantaneous valley population imbalance $\\Delta N_{\\mathrm{CB1}} = N^{K_1}_{\\mathrm{CB1}} - N^{K_2}_{\\mathrm{CB1}}$ in the first conduction band. The proportionality arises because the probe couples the first and second conduction bands in a valley-selective way, so a prepared valley-polarization state makes one helicity excite more electrons than the other. The authors verify the linear mapping numerically for h-BN and MoS2 by varying the pump strength and fitting $\\Delta\\mu$ at the CB1–CB2 resonance against the directly counted $\\Delta N_{\\mathrm{CB1}}$. They then use the mapping to retrieve the time-dependent valley polarization during multi-pulse switching sequences, resolving the switching steps with 250 as resolution.","pith_inferences":["The same circular-dichroism observable at the CB1–CB2 resonance could be used to monitor spin-valley locking in transition-metal dichalcogenides, provided the spin-split matrix elements still obey the equality assumed in Eq. (1); this extension would need explicit verification.","A residual $\\Delta\\mu$ in an unpolarized or symmetrically pumped sample would expose intrinsic asymmetry in the CB1–CB2 transition strengths, offering a clean experimental test of the load-bearing assumption.","Because the probe photon energies lie in the ultraviolet/XUV range (7.5 eV for h-BN, 35 eV for MoS2), excitonic and many-body effects may alter the CB1–CB2 oscillator strengths; including them could change the proportionality constant and would need to be checked in a full many-body treatment."],"forward_implications":["Valley polarization dynamics can be retrieved quantitatively from a single differential absorption measurement, without requiring a separate electron-spectroscopy probe during the ultrafast evolution.","The 250 as delay step resolves the transient switching of valley polarization between the K1 and K2 valleys, not merely the polarization state after switching has completed.","The linear mapping is confirmed for two different materials, h-BN and MoS2, indicating the scheme transfers to other 2D hexagonal materials with valley-contrasting band structure.","With one auxiliary absolute measurement (e.g., ARPES) under static or slowly varying valley polarization, the proportionality constant can be calibrated, making the ultrafast valley-polarization readout absolute."],"supporting_citations":[{"why":"Supplies the optical valley selection rules that make both the pump excitation and the probe transitions valley-selective.","marker":"[3]"},{"why":"Gives the valley-contrasting band structure and selection rules for 2D hexagonal materials, grounding the proportionality in Eq. (1).","marker":"[4]"},{"why":"Provides the transient absorption spectroscopy formalism used to compute the absorption coefficients $\\mu_+$ and $\\mu_-$.","marker":"[25]"},{"why":"Demonstrates experimental feasibility of circularly polarized attosecond probe pulses, which the scheme requires.","marker":"[28]"},{"why":"Supports the claim that valley-selective transitions in higher conduction bands are universal for valley-active materials.","marker":"[30]"},{"why":"Provides the time-dependent density functional theory implementation used for all numerical pump-probe experiments.","marker":"[39]"}],"fun_headline_variants":["Attosecond CD maps valley polarization in 250 as","All-optical valley readout with 250-attosecond resolution","Circular dichroism retrieves valley dynamics at 250 as","Sub-femtosecond valley polarization read via CD","Ultrafast valley polarization measured in 250 attoseconds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The retrieval assumes that the CB1-to-CB2 transition has equal matrix-element magnitudes for the two probe helicities at the K1 and K2 valleys, so that any absorption asymmetry comes purely from the valley population imbalance and not from material-specific differences in transition strength.","fun_headline_variants_meta":{"raw":{"variants":["Attosecond CD maps valley polarization in 250 as","All-optical valley readout with 250-attosecond resolution","Circular dichroism retrieves valley dynamics at 250 as","Sub-femtosecond valley polarization read via CD","Ultrafast valley polarization measured in 250 attoseconds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000637,"raw_usage":{"total_tokens":2900,"prompt_tokens":871,"completion_tokens":2029,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":487,"completion_tokens_details":{"reasoning_tokens":1948}},"tokens_in":487,"tokens_out":2029,"duration_ms":12754,"temperature":1.0,"reasoning_tokens":1948,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:09:23.072311+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Perform the same pump-probe simulation on a material with slightly broken K/K' symmetry, such as a strained h-BN monolayer, and check whether the proportionality constant in Eq. (1) stays time-independent when the directly counted $\\Delta N_{\\mathrm{CB1}}$ is compared with $\\Delta\\mu$; a drift would show that the linear mapping is not universal.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the transient absorption spectroscopy formalism used to compute the absorption coefficients $\\mu_+$ and $\\mu_-$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Demonstrates experimental feasibility of circularly polarized attosecond probe pulses, which the scheme requires."},{"cited_title":"Geondzhian, A","cited_arxiv_id":null,"evidence_quote":"Supports the claim that valley-selective transitions in higher conduction bands are universal for valley-active materials."},{"cited_title":"Tancogne-Dejean, M","cited_arxiv_id":null,"evidence_quote":"Provides the time-dependent density functional theory implementation used for all numerical pump-probe experiments."}],"review_version":1}