{"id":"219fa773-738d-420a-8e9d-9039063e92cc","arxiv_id":"2508.07612","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A sub-cycle optical selection rule links instantaneous light chirality to valley selectivity, enabling independent control of valley currents.","lead":"This paper proposes that the instantaneous rotation direction of a laser field controls which of two electronic valleys in a two-dimensional material gets preferentially excited. Using specially shaped chirality-separated pulses, the authors show in simulations that valley currents can be sent in different directions or canceled to zero net charge flow.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (4)'s symmetric-gap approximation E(K+A(ti))≈E(K'+A(ti)) is not symmetry-protected and may fail for strong fields; the quantitative instantaneous selection rule needs a direct check.","rationale":"The paper's central mechanism is the instantaneous optical valley selection rule in Eq. (4), and the applications depend on this formula. The main-text derivation is not shown; the reader flagged this. I focused on the explicit symmetry approximation Ebar_g = Ecv[K+A(ti)] ≈ Ecv[K'+A(ti)]. This is the load-bearing point because if it fails, the tanh argument in Eq. (4) is wrong even if the sign of c(ti)ζg still gives qualitative valley preference. Time-reversal symmetry relates the valleys as E(K+q)=E(K'-q), not E(K'+q); the equality used in Eq. (4) requires the band dispersion around K to be symmetric under q→-q, i.e., no trigonal warping, which is generically absent in hexagonal 2D materials. At the high intensities and long wavelengths used, A(ti) can be a significant fraction of the Brillouin-zone scale, so the approximation is not obviously perturbative. The numerical demonstrations in Figs. 2–5 are consistent with the sign rule, but they do not separately verify the quantitative tanh formula or scan over intensity where the approximation degrades. A direct computation of the gap asymmetry and a comparison of Eq. (4) with full simulation would settle this. If the check passes, the central claim stands; if not, the paper's quantitative rule needs revision, though the qualitative chirality-control idea may survive. This does not change the reader's conditional verdict.","tokens_in":9436,"tokens_out":5724,"duration_ms":70061,"concrete_test":"Using the same two-band model and parameters as in Figs. 2 and 5 (including peak intensities 5×10^10 and 1×10^12 W/cm^2 and wavelengths 4000/2000 nm), compute δ(t)=|Ecv[K+A(t)]−Ecv[K'+A(t)]|/(average gap) over one optical cycle. If at any ionization time ti where wavelets are launched in Fig. 2 δ exceeds, say, 10%, Eq. (4)'s symmetric-gap approximation is violated in the demonstrated regime. Then re-derive the tanh expression without the Ebar_g equality and compare with the full numerical Nc(K), Nc(K') from the semiconductor Bloch equations; if the full calculation disagrees with Eq. (4) by more than the small-gap error, the instantaneous rule needs modification.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative prediction, Eq. (4), rests on the approximation Ebar_g = Ecv[K+A(ti)] ≈ Ecv[K'+A(ti)]. This is not guaranteed by the symmetries invoked in the paper. Time-reversal symmetry relates the two valleys as Ecv(K+q) = Ecv(K'-q) (up to reciprocal lattice vectors), not Ecv(K'+q). The equality used in Eq. (4) would require the dispersion near K to be invariant under q → -q, i.e., absence of trigonal warping. Such warping is generically present in hexagonal 2D materials, and its effect grows with |A(ti)|. At the strongest intensities considered (up to 1e12 W/cm^2 at a 4000-nm fundamental), A(ti) is not necessarily small on the Brillouin-zone scale, so the gap values at K+A(ti) and K'+A(ti) can differ substantially. If Ebar_g is mis-estimated by a large fraction, the tanh argument in Eq. (4) changes nonlinearly, so the predicted ηvp from a wavelet at ti is quantitatively unreliable. The sign of c(ti) may still control the preferred valley, but the claimed quantitative rule, and any purity values derived from it, remain unverified. The derivation itself is only in the Supplemental Material, so the main text does not establish the limits of Eq. (4).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an instantaneous optical valley selection rule in two-dimensional hexagonal systems: the population imbalance between K and K′ valleys produced by an ionization wavelet at time t_i is given by η_vp = -tanh(2 c(t_i) ζ_g Ebar_g / g0) (Eq. 4), where c(t_i) is the instantaneous optical chirality, ζ_g is the system chirality, and Ebar_g is the (approximately valley-symmetric) gap. This rule is used to design 'chirality-separated' bicircular fields, with opposite subcycle optical chirality and biased vector-potential lobes, claimed to enable independent control of K and K′ valley currents. The authors demonstrate numerically the rotation of valley current directions with phase φ, the complete orthogonal separation of valley currents (P reaching ±1), and the generation of a pure valley current with zero net charge (Q = 1 to four decimal places).","tokens_in":9782,"tokens_out":6404,"duration_ms":75892,"significance":"If established, this is a valuable step beyond cycle-averaged valley selection: it identifies the subcycle ionization time as the control handle and shows how a single field can address the two valleys separately. The explicit closed-form selection rule, the clear field-design rationale, and the numerical demonstrations are strengths. The sign-based mechanism is physically plausible and consistent with prior work on interband excitation. However, the quantitative validity of Eq. (4) is not sufficiently supported in the present manuscript, because its derivation and the symmetric-gap approximation are not presented, and the role of the parameter g0 is unspecified. These gaps prevent full confidence in the claimed 100% purity and pure-valley-current results.","major_comments":[{"comment":"Equation (4) is the quantitative core of the paper, but its derivation is only cited to the Supplemental Material and the approximation Ebar_g = Ecv[K+A(ti)] ≈ Ecv[K′+A(ti)] is stated without proof or validity conditions. Time-reversal symmetry relates Ecv(K+q) to Ecv(K′−q), not to Ecv(K′+q); the equality therefore requires additional conditions, e.g., negligible trigonal warping. At the strongest intensities considered (up to 1e12 W/cm² at 4000 nm, Fig. 5(d)), A(ti) is not parametrically small on the Brillouin-zone scale, so the gap values can differ substantially. The sign of c(ti)ζg may still select the valley, but the quantitative tanh rule and the purity values derived from it are unverified. Please include the derivation or a direct validation of Eq. (4) against the full numerical model for multiple ionization instants and intensities, and state the regime of validity.","section":"Instantaneous optical valley selection rule, Eq. (4)"},{"comment":"The formula contains g0, described only as 'a constant that describes the width of the ionization wavelet.' The main text does not state how g0 is determined. If g0 is fitted to simulations, the agreement in Fig. 2 is not an independent test of the selection rule; if it is computed from first principles, the defining expression should be given. Without this information, the quantitative predictions (including the 100% purity claims in Applications 1 and 2) are not reproducible from the main text alone.","section":"Eq. (4) and Fig. 2"},{"comment":"The pure-valley-current claim rests on the statement that Q = 1 to four decimal places over the intensity range. No convergence tests, numerical grid parameters, or error estimates are provided in the main text. Because the cancellation Jcharge = 0 is partly enforced by mirror symmetry, it would be helpful to state whether Q = 1 is a symmetry-protected result or a numerical outcome, and to specify the numerical accuracy of the computed currents. Similarly, the P = ±1 points in Fig. 4(d) should be clarified: they appear to follow from projection orthogonality of the two valley currents rather than from a valley-resolved carrier-counting purity, which should be stated explicitly to avoid overinterpretation.","section":"Fig. 5(d) and Applications 1–2"}],"minor_comments":[{"comment":"Typo: 'nonvalishing valley currents' should be 'nonvanishing valley currents'.","section":"Introduction"},{"comment":"The caption labels panels as (a) and (c-d), but the text refers to 'Figs. 3(b) and (c)'. Renumber the figure panels and fix the cross-references.","section":"Fig. 3 and accompanying text"},{"comment":"Please define all symbols in one place, especially Rcv and dcv, and clarify the notation n(t) = F(t)/||F(t)|| and ˙n(t). The physical dimensions of the three terms in the integrand should be stated.","section":"Eq. (1)"},{"comment":"The sign convention for εco± is confusing as written: 'εco± = ∓[...]' together with 'co-rotating' should be explained, and the relationship between the handedness of the bases e± and the IOC sign should be made explicit.","section":"Eq. (5)"},{"comment":"The 'dashed line' marking the Lissajous figure is not described in the caption; specify which curve corresponds to the driving field and where t1, t2, t3 are located.","section":"Fig. 2(a)"}],"recommendation":"major_revision","confidential_remarks":"The central quantitative result, Eq. (4), is not established in the main text; the Supplemental Material and the determination of g0 are essential. The editor should ask the authors to provide the full derivation, a direct numerical test of Eq. (4) over the relevant parameter range, and numerical convergence details. The novelty relative to the authors' own previous ionization-wavelet formalism (refs. [46,47]) and bicircular-field valley control (ref. [27]) should also be clarified explicitly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read. The genuinely new thing is the instantaneous optical valley selection rule: use the sign of the instantaneous chirality c(ti) at ionization time, rather than the cycle-averaged helicity, to pick K versus K'. That is a real step beyond the bicircular-field work, and the chirality-separated field construction is elegant. The two applications—complete spatial separation of valley currents and pure valley current with zero net charge—follow from the same picture, and the numerics look consistent.\n\nThe soft spot is exactly the one in the stress-test. Eq. (4) uses E[K+A(ti)] ≈ E[K'+A(ti)]. That is not symmetry-protected. Time reversal relates the valleys as E(K+q)=E(K'-q), so the equality holds only if the dispersion near each valley is locally even in q. Trigonal warping breaks that, and the difference grows with |A|. At the strongest fields quoted (10^12 W/cm^2, 4000 nm), A is not small on the BZ scale, so the tanh argument in Eq. (4) could be off by a large factor. The sign rule likely survives—that is the main claim—but the quantitative population imbalance and the 100% purity numbers are model-dependent. The authors do not claim symmetry protection; they just introduce the approximation without justifying its range. The derivation is in the SM, which I could not see. A referee must check it.\n\nAlso worth saying: the numerical demonstration is not an independent test. The ionization-wavelet picture and the interband phase are the authors' own framework, so the simulations confirm the framework without falsifying it. That is fine for a theory paper, but it means the impressive purity numbers are not external predictions.\n\nCredit: the paper is clearly written, cites the prior work fairly, and flags the approximation rather than hiding it. The field design and the two example applications are new and useful.\n\nRecommendation: send it to a serious referee. Ask specifically whether Eq. (4) holds against a full band-structure calculation with trigonal warping, and where the approximation breaks. If it survives at the quoted intensities, this is a strong PRL-class result. If not, the sign rule stays, but the quantitative claims need revision.","headline":"The instantaneous chirality-based valley selection rule is a real step forward, but Eq. (4) rests on an unverified symmetric-gap approximation that may break at strong fields.","tokens_in":10246,"tokens_out":4620,"would_cite":true,"duration_ms":48916,"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 establishes an instantaneous optical valley selection rule: the sign of the field's optical chirality at the moment of ionization determines whether the K or K' valley is excited, and a single chirality-separated field can then dr","keywords":["valleytronics","instantaneous optical chirality","optical valley selection rule","valley current","chirality-separated field","bicircular field","two-dimensional materials","sub-cycle control"],"falsifier":"Solve the full time-dependent Schrödinger equation for a two-band hexagonal lattice driven by one of the chirality-separated fields of Fig. 3 at an intensity where $|\\mathbf{A}(t_i)|$ is large enough that $E_{cv}[\\mathbf{K}+\\mathbf{A}(t_i)]$ differs from $E_{cv}[\\mathbf{K}'+\\mathbf{A}(t_i)]$ by more than about 10%. Compute $\\eta_{vp}$ at instants like $t_1, t_2, t_3$ of Fig. 2. If the imbalance no longer tracks $-\\tanh(2 c(t_i) \\zeta_g \\bar{E}_g / g_0)$, or if changing $\\varphi$ by $\\pi$ fails to reverse the K/K$'$ current directions as predicted, the instantaneous rule is only an approximatio","tokens_in":9371,"feed_emoji":"🌀","tokens_out":6307,"duration_ms":69089,"temperature":0.7,"pith_summary":"The paper claims that the old optical valley selection rule--which valley gets excited by circularly polarized light--has an instantaneous version: at the exact moment an electron makes an interband transition, the sign of the field's instantaneous optical chirality, the rotation sense of its polarization direction, decides whether the K or K' valley is populated. The central result is a compact formula (Eq. 4) for the valley-population imbalance, $\\eta_{vp} = -\\tanh(2 c(t_i) \\zeta_g \\bar{E}_g / g_0)$, so the rule is quantitative, not just qualitative. The authors then argue that a single synthesized field, made from two co-rotating bicircular components, can carry opposite chirality signs in different sub-cycle lobes while biasing those lobes in different directions. That is enough to push the K and K' electrons along independently chosen paths within one optical cycle, yielding 100%-purity valley-polarized currents and pure valley currents with zero net charge. A sympathetic reader would care because this turns the valley degree of freedom from a synchronized binary switch into a per-cycle, per-valley control dial.","feed_headline":"Instant chirality sets which valley carries current","feed_subtitle":"One chirality-separated pulse drives K and K' electrons in different directions, reaching 100% pure valley currents.","key_machinery":"Equation (4), the tanh valley-polarization formula, is the load-bearing object. It follows from rewriting the third term of the accumulated interband phase $S[k(t)]$ as an integral over $c(t')\\, \\zeta[k(t')]$; the chirality of the two-band system is $\\zeta \\approx +1$ near K and $-1$ near K$'$, and the field chirality is $c(t) = n_y \\partial_t n_x - n_x \\partial_t n_y$. The formula converts the sign of $c$ at ionization time into a population imbalance, and the factor $g_0$ (the ionization-wavelet width) sets the scale over which the rule saturates. On the field-construction side, the co-rotating bicircular basis of Eq. (5) supplies the chirality-separated waveform: the two lobes of the Liss","core_discovery":"On the paper's terms, the central discovery is the instantaneous optical valley selection rule. In the interband transition phase, the term involving the time derivative of the field's polarization direction factorizes into the product of the instantaneous optical chirality $c(t)$ and the valley-dependent chirality $\\zeta(k)$. For an electron born at time $t_i$ via an ionization wavelet, each valley's population is set by a tanh law, $\\eta_{vp} = -\\tanh(2 c(t_i) \\zeta_g \\bar{E}_g / g_0)$, so the sign of $c(t_i)$ alone selects the valley. Because oppositely signed chirality can be concentrated in different halves of an optical cycle of a synthesized field (two co-rotating bicircular beams of","pith_inferences":["The sign-only character of Eq. (4) suggests that any waveform whose Lissajous figure contains sub-cycle segments of opposite rotation sense--not just the two-color co-rotating construction--could serve as a valley router, so harmonic ratios beyond 2:1 or polarization gating may generalize the scheme.","Because the tanh saturates, the purity of each lobe's selection can be tuned continuously by lobe intensity rather than only by waveform geometry; the authors' 100% and zero-charge demonstrations are the two endpoints of a continuous family.","The derivation assumes the ionization-wavelet picture and a symmetric gap at the shifted momenta, so it is likely most faithful for low-frequency fields and moderate intensities; testing the rule with few-cycle mid-infrared pulses where $\\mathbf{A}(t)$ is strong would reveal whether a correction term proportional to $E_{cv}[\\mathbf{K}+\\mathbf{A}] - E_{cv}[\\mathbf{K}'+\\mathbf{A}]$ is needed.","An attosecond pump-probe experiment that measures the K/K$'$ emission asymmetry as a function of carrier-envelope phase could directly read out the instantaneous chirality $c(t_i)$ of a pulse, turning the rule into a diagnostic of sub-cycle field structure."],"forward_implications":["A single chirality-separated field can produce valley-polarized currents with purity $P = \\pm 1$ (100%) by detecting along directions orthogonal to the two valley currents.","Pure valley current with zero net charge transport ($Q = 1$ to at least four decimal places) can be generated robustly against laser intensity fluctuations over two orders of magnitude.","The K and K$'$ currents can be rotated independently by tuning the relative phase $\\varphi$ of the fundamental and second-harmonic components, including the case where one valley's current stays fixed while the other rotates arbitrarily.","The independent control operates on the optical-cycle timescale (a 12-$T_0$ pulse with a 10-$T_0$ plateau), much faster than approaches relying on electrodes or heterostructures."],"supporting_citations":[{"why":"Establishes the valley-dependent optical selection rule based on field helicity that this work extends to the instantaneous level.","marker":"[30]"},{"why":"Supplies the two-band gapped-graphene model and the interband phase-integral approach used in the derivation.","marker":"[35]"},{"why":"Provides the subcycle optical chirality concept and demonstrates its experimental control, which Eq. (2) adapts.","marker":"[44]"},{"why":"Introduces the ionization-wavelet picture in which Eq. (4) is derived.","marker":"[46]"},{"why":"Gives the companion particle-wave perspective on the wavelet picture that supports the derivation.","marker":"[47]"},{"why":"Supplies the co-rotating bicircular field basis used to construct the chirality-separated field.","marker":"[27]"}],"fun_headline_variants":["Instant chirality selects valley for current flow","Chirality-separated light drives independent valley currents","Single pulse achieves 100% pure valley-polarized currents","Zero-net-charge valley current from one optical pulse"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"The tanh formula assumes that the energy gap seen at the two momentarily shifted valleys is the same, $\\bar{E}_g = E_{cv}[\\mathbf{K}+\\mathbf{A}(t_i)] \\approx E_{cv}[\\mathbf{K}'+\\mathbf{A}(t_i)]$, and that a single constant $g_0$ describes the ionization-wavelet width; if a strong vector potential breaks this symmetry, the clean sign-only selection rule can be corrupted.","fun_headline_variants_meta":{"raw":{"variants":["Instant chirality selects valley for current flow","Chirality-separated light drives independent valley currents","Single pulse achieves 100% pure valley-polarized currents","Zero-net-charge valley current from one optical pulse"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00022,"raw_usage":{"total_tokens":1254,"prompt_tokens":684,"completion_tokens":570,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":428,"completion_tokens_details":{"reasoning_tokens":507}},"tokens_in":428,"tokens_out":570,"duration_ms":6557,"temperature":1.0,"reasoning_tokens":507,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T21:58:23.844873+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the full time-dependent Schrödinger equation for a two-band hexagonal lattice driven by one of the chirality-separated fields of Fig. 3 at an intensity where $|\\mathbf{A}(t_i)|$ is large enough that $E_{cv}[\\mathbf{K}+\\mathbf{A}(t_i)]$ differs from $E_{cv}[\\mathbf{K}'+\\mathbf{A}(t_i)]$ by more than about 10%. Compute $\\eta_{vp}$ at instants like $t_1, t_2, t_3$ of Fig. 2. If the imbalance no longer tracks $-\\tanh(2 c(t_i) \\zeta_g \\bar{E}_g / g_0)$, or if changing $\\varphi$ by $\\pi$ fails to reverse the K/K$'$ current directions as predicted, the instantaneous rule is only an approximatio","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the valley-dependent optical selection rule based on field helicity that this work extends to the instantaneous level."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the two-band gapped-graphene model and the interband phase-integral approach used in the derivation."},{"cited_title":"Rozen, A","cited_arxiv_id":null,"evidence_quote":"Provides the subcycle optical chirality concept and demonstrates its experimental control, which Eq. (2) adapts."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the ionization-wavelet picture in which Eq. (4) is derived."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the companion particle-wave perspective on the wavelet picture that supports the derivation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the co-rotating bicircular field basis used to construct the chirality-separated field."}],"review_version":1}