{"id":"418209df-4c13-4df9-8759-bfae2c25e2cc","arxiv_id":"2504.14236","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Elliptically polarized intense laser pulses generate ellipticity-dependent bulk photocurrents in 2D hexagonal materials, with a gap-dependent signal that could be used for spectroscopy.","lead":"This paper uses computer simulations to show that elliptically polarized laser pulses can generate controllable electric currents in two-dimensional honeycomb materials like boron nitride. The current's strength and direction depend on the laser's ellipticity and orientation, and the pattern could reveal the material's electronic band gap.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Gap-spectroscopy claim rests on a single-parameter Gaussian model; the one ab initio check (hBN, 4.2 eV) already disagrees with the model on the orientation of maximum photocurrent, so the small-gap scaling used for spectroscopy is unvalidated.","rationale":"The reader's weakest_assumption correctly identifies the transferability of the Gaussian model for gap spectroscopy as the central vulnerability. I considered whether a more fundamental issue exists, such as the absence of dephasing in the current definition, but that is standard for this class of strong-field simulations and does not single out the paper's claim. The model-versus-TDDFT orientation discrepancy at hBN is the sharpest evidence that multi-band and non-K/K' physics matter. Since the proposed spectroscopy is built on a model that demonstrably misses the correct orientation at the only ab initio point, the load-bearing concern is real. The core BPG effect with elliptical pulses is well supported by symmetry and by the qualitative agreement of the bell-shaped ellipticity curves, so the verdict remains conditional rather than reject. The proposed concrete test - a second ab initio point at a smaller gap - would directly confirm or refute the interpolation, making it the single check that settles whether the concern lands.","tokens_in":13848,"tokens_out":7116,"duration_ms":63990,"concrete_test":"Run time-dependent density-functional theory for a small-gap hexagonal 2D material (e.g., monolayer MoS2, LDA gap ~1.8 eV) with the same 800 nm, 0.3 TW/cm2 elliptical pulse (epsilon around 0.5, theta at 0, pi/6, pi/4, pi/2), and compare the orientation of maximum current, the sign pattern of Ilong(epsilon), and the peak amplitude with the Gaussian-model prediction at the same gap. If the orientation discrepancy seen at hBN (Fig. 7) persists or the gap-vs-current scaling (Fig. 5c) does not transfer, the spectroscopy claim is not supported beyond the model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's primary claim of elliptically controlled BPG currents is supported by symmetry and by both model and TDDFT results (bell-shaped epsilon response, zero current at linear/circular, Fig. 3 and Fig. 7a,c). The load-bearing problem is the gap-spectroscopy extension. All gap-dependent trends (Figs. 5-6) come from a real-space Gaussian-potential model with two bands and no electron interactions. The sole ab initio calibration, hBN at LDA gap ~4.2 eV, shows the model fails to reproduce the orientation of maximum current (TDDFT max at theta=0, model max at theta=pi/2; Fig. 7). The authors attribute this to multi-band contributions away from K/K' - precisely the physics excluded from the model that the gap-scaling extrapolation (0-0.12 eV linear regime and 0.3-1.5 eV double-peak regime) relies on. With no error bars, no data/code release, and no second material point, the 'photocurrent indicative of gap size' claim is an unsupported interpolation. A further internal inconsistency appears in the conclusions: 'amplitude that scales linearly with the laser power' contradicts the Results, which show linear scaling with gap size, not laser power.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript numerically investigates bulk photogalvanic (BPG) currents in two-dimensional hexagonal materials driven by monochromatic, elliptically polarized intense laser pulses. Using a two-band real-space Gaussian-potential model solved by time-dependent Schrödinger equation (TDSE) and ab initio time-dependent density functional theory (TDDFT) for monolayer hBN, the authors show that such pulses generate BPG currents when inversion symmetry is broken, with vanishing current for purely linear and circular polarization and a bell-shaped ellipticity dependence peaking near ellipticity values of about 0.4–0.5. They further report that the current amplitude and direction depend on the orientation of the major elliptical axis, that the current decomposes into longitudinal and transverse (Hall-like) components, and that the current scales linearly with the band gap for small gaps. They propose this behavior as a basis for photocurrent-based gap spectroscopy and benchmark the model against TDDFT for monolayer hBN.","tokens_in":14102,"tokens_out":5951,"duration_ms":50879,"significance":"If the results hold, the paper establishes a new set of control knobs—laser ellipticity and in-plane orientation—for strong-field BPG currents in 2D hexagonal materials, and it suggests a transport-based probe of band gaps. The central qualitative effect is supported by two independent computational methods (model TDSE and TDDFT) and by a symmetry analysis of momentum-resolved conduction-band occupations. The authors are also candid about the model/TDDFT discrepancy in the orientation of maximum current and about the unknown origin of the double-peak structure. However, the gap-spectroscopy extension currently rests on model-only data with a single ab initio comparison at a gap where the model disagrees with TDDFT on a qualitative feature, so the extrapolation to a general gap-spectroscopy scheme is not yet established.","major_comments":[{"comment":"The central gap-spectroscopy claim rests on model-only data. The linear small-gap scaling in Fig. 5c and the double-peak regime in Fig. 6 are generated entirely by the two-band Gaussian-potential model of Methods Eq. (1)–(3), yet the sole ab initio benchmark, monolayer hBN at the LDA gap of about 4.2 eV, already disagrees with the model on a qualitative feature: TDDFT finds the maximum current for θ=0 (Fig. 7a), while the model predicts it for θ=π/2 (Fig. 7c). The authors attribute this discrepancy to multi-band and non-K/K' contributions, which are exactly the physics excluded from the model used to extrapolate the gap dependence. Without a second material point or a TDDFT check at small gaps, the proposed photocurrent-based gap spectroscopy is not validated outside the model.","section":"Gap spectroscopy; Monolayer hBN (Figs. 5–7)"},{"comment":"No convergence or uncertainty assessment is reported for the central quantitative results. The photocurrent amplitudes in Figs. 3, 5, and 6 are quoted without error bars, and the linear fit constants a=0.0087 and b=0.0165 in Fig. 5c and the regime boundaries 0.12 eV, 0.3 eV, and 1.5 eV are presented without tests of sensitivity to the grid spacing (0.28 Bohr), k-grid (100×100), time step (0.2 a.u.), pulse duration, or the averaging window (six cycles for 800 nm versus two cycles for 3000 nm). Because these numbers underpin the spectroscopy claim, a convergence study or numerical uncertainty estimate is required.","section":"Methods: Model calculations; Gap spectroscopy"},{"comment":"The Conclusions state that the photocurrent signal 'scales linearly with the laser power,' but the Results and Discussion section 'Gap spectroscopy' shows scaling with the gap size (Fig. 5c), not with laser power. This is an internal inconsistency in a central claim; if power scaling is intended it is not demonstrated anywhere in the paper, and if gap scaling is meant the sentence should be corrected.","section":"Conclusions"},{"comment":"The double-peak structure in the ellipticity dependence is proposed as a 'fingerprint sign for quantum interference spectroscopy,' yet the same paragraph states that the physical origin of the double peak 'currently remains unclear.' Since the mechanism is unidentified and the feature appears only in the model calculation, the evidential basis for using it as a spectroscopic fingerprint is not established.","section":"Gap spectroscopy (Fig. 6)"}],"minor_comments":[{"comment":"The supersine envelope in Eq. (5) is typeset ambiguously; the exponent is unclear. Please rewrite the formula so that the envelope can be reproduced unambiguously.","section":"Methods, Eq. (5)"},{"comment":"The parameter σ_B is never given a numerical value except implicitly for the gapless case; state the value used for the gapped sublattice B.","section":"Methods, Eq. (1)"},{"comment":"The caption refers to panel (c) for the lineouts, but the figure contains only panels (a) and (b); fix the caption/panel labels.","section":"Fig. 6 caption"},{"comment":"The phrase 'excellent agreement' overstates the level of agreement given the orientation discrepancy identified in the same section; 'qualitative agreement in the ellipticity dependence and the longitudinal-component behavior' would be more precise.","section":"Monolayer hBN"},{"comment":"No data availability statement or code release is included; making the model parameters and TDDFT input files available would strengthen reproducibility of the numerical photocurrent values.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The primary claim of ellipticity-controlled BPG currents is credible and likely publishable after revision, but the gap-spectroscopy claim needs either additional ab initio validation at small gaps or a clear reduction in claim strength. The issues are fixable within the scope of the manuscript, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"For your internal take: this paper is worth a serious referee slot, but the spectroscopy promise is over-sold relative to the evidence.\n\nWhat is actually new and good: the systematic map of how monochromatic elliptical driving controls BPG currents in broken-inversion hexagonal lattices. The zero-current at linear and circular polarization, the bell-shaped ellipticity response peaking around epsilon 0.3–0.6, and the appearance of a transverse Hall-like component are all cleanly demonstrated in both a two-band tight-binding-style model and TDDFT for hBN. The symmetry discussion in Fig. 4 is clear and convincing, and the two-method cross-check is a genuine strength. The gap-dependent sign changes in the longitudinal current for small gaps (Fig. 5) are interesting and, as the authors say, not what you'd guess from semiclassical dynamics. That part is a real, citable contribution for people working on strong-field transport in 2D materials.\n\nWhere it gets soft: the gap-spectroscopy claim. All the gap-dependent trends (Figs. 5–6) come from a Gaussian-potential model tuned through one parameter. The only ab initio calibration is hBN at its LDA gap (~4.2 eV), and there model and TDDFT disagree on which crystal orientation maximizes the current: TDDFT says theta=0, model says theta=pi/2. The authors attribute this to multi-band contributions away from K/K', which is plausible but also exactly the physics their small-gap extrapolation ignores. So the claim that photocurrent ellipticity curves can be used to read off the gap from, say, a Floquet-engineered graphene sample is an interpolation across a regime that has never been tested against a real material. That does not sink the paper, but it should be labeled as a proposal, not a demonstrated technique.\n\nThere is also a real internal inconsistency in the conclusions: it says the photocurrent amplitude 'scales linearly with the laser power' when the results show scaling with gap size. Probably a typo, but it needs fixing. And the paper has no convergence plots, no error bars, and no code/data release, which matters for a purely numerical study where the quantitative claims are the product.\n\nVerdict: engage with it. The central qualitative result—elliptically controlled BPG currents with a transverse component—is solid and well-explained. The spectroscopy extension needs more evidence, ideally a second material or a controlled comparison with a known time-dependent gap. I would tell the editor: send it to review, with a request for the convergence data and a rewrite of the gap-spectroscopy language so it matches what is actually shown.","headline":"A solid numerical study of ellipticity-controlled bulk photocurrents in 2D hexagonal materials, but the gap-spectroscopy claim is an extrapolation supported by only one ab initio point that already disagrees with the model on a key detail.","tokens_in":14654,"tokens_out":1568,"would_cite":true,"duration_ms":15798,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Elliptical light pulses drive controllable photocurrents in 2D hexagonal materials, and the current's ellipticity pattern reflects the band gap.","keywords":["bulk photogalvanic effect","2D hexagonal materials","elliptically polarized light","Hall-like photocurrent","gap spectroscopy","strong-field dynamics","time-dependent density functional theory","monolayer hBN"],"falsifier":"Measure the ellipticity-dependent photocurrent in a monolayer hBN sample under 800 nm driving, comparing the current for the major axis along zigzag and along armchair. The model predicts the armchair orientation carries the maximum current, while the ab initio calculation predicts zigzag; if experiment follows the ab initio result, the model-based extrapolation to other gap sizes is not validated. A stronger test is to tune the gap continuously by strain or alloying and check whether the small-gap linear scaling and sign-change pattern predicted by the model appear in the measured current.","tokens_in":13610,"feed_emoji":"⚡","tokens_out":7511,"duration_ms":64163,"temperature":0.7,"pith_summary":"This paper argues that a single-color, elliptically polarized laser pulse is enough to create bulk photogalvanic currents in two-dimensional hexagonal materials that lack inversion symmetry. The current has both a longitudinal and a Hall-like transverse component, and both can be steered by the laser's ellipticity, polarization angle, wavelength, and intensity. The paper's central proposal is that the photocurrent's dependence on ellipticity is a material fingerprint: as a band gap opens at the $K$ and $K'$ valleys, the ellipticity curve changes sign repeatedly and its amplitude grows linearly with the gap for small gaps. If correct, this makes photocurrent detection a compact, all-optical way to measure band gaps in 2D materials and to probe how gaps change under ultrafast or Floquet driving.","feed_headline":"Elliptical pulses drive controllable photocurrents in 2D materials","feed_subtitle":"The current's ellipticity pattern scales with the band gap, making photocurrent readout a compact gap spectroscopy.","key_machinery":"The central object is the bulk photogalvanic current, a steady current generated by light absorption in a non-centrosymmetric material without an applied voltage. The mechanism is the symmetry of the light-induced conduction-band occupation in momentum space: the paper computes the k-resolved occupation after the pulse and shows that elliptical polarization leaves only one mirror plane, producing an occupation imbalance that flows as a net current, while linear and circular polarization leave symmetries that cancel it. The modeling machinery is a real-space honeycomb-lattice Hamiltonian with Gaussian on-site potentials; varying the A/B sublattice potential tunes the $K$/$K'$ gap from 0 to 2 eV, and the time-dependent Schrödinger equation is solved in the velocity gauge to obtain the current. Ab initio validation uses time-dependent density functional theory in the adiabatic approximation for monolayer hBN.","core_discovery":"Monochromatic elliptical pulses generically drive bulk photogalvanic currents in broken-inversion-symmetric 2D hexagonal lattices, because elliptical polarization is the minimal symmetry-breaking drive: linear polarization leaves mirror symmetries in the conduction-band occupation that cancel the current, and circular polarization leaves a threefold valley symmetry, while elliptical polarization leaves only a single mirror plane and therefore a net current. The current vanishes for linear and circular light, peaks for ellipticities around 0.3-0.6, and its direction is set by the major axis of the ellipse, with zigzag-oriented driving producing a longitudinal current and armchair-oriented driving producing a transverse Hall-like current. As the sublattice asymmetry grows from zero to a 2 eV gap, the ellipticity dependence of the longitudinal current develops sign changes and double-peak structures, and for small gaps (below about 0.12 eV) both components scale linearly with the gap. The same bell-shaped ellipticity response is found in ab initio simulations of monolayer hBN, supporting the proposal that ellipticity-dependent photocurrents can serve as a gap spectroscopy.","pith_inferences":["If the linear small-gap scaling survives in real materials, ellipticity-dependent photocurrents could be used to time-resolve a dynamically opened gap, for instance in a Floquet-engineered graphene-like system, by reading the gap from the current amplitude during the pump.","The transverse Hall-like component may offer a transport-based route to Berry-curvature spectroscopy, since the current direction is tied to valley occupation imbalance; this connection is implicit in the paper but not directly measured.","The unexplained double-peak structure at intermediate gaps is a natural target for a two-photon resonance test: if the peak positions shift with laser wavelength according to multi-photon energies, the interference picture would be confirmed.","Combining elliptical driving with a weak second harmonic could expose interference between the two symmetry-breaking mechanisms, potentially increasing current magnitude or adding new control axes beyond monochromatic ellipticity."],"forward_implications":["Ellipticity, polarization angle, wavelength, and intensity become independent knobs for directing both the longitudinal and the transverse (Hall-like) photocurrent in 2D hexagonal materials.","A photocurrent measurement at fixed driving conditions can report the $K$/$K'$ gap size, with small gaps producing a linear current-vs-gap scaling and larger gaps producing structured multi-peak ellipticity curves.","The vanishing of the current for linear and circular polarization, and its bell-shaped maximum near intermediate ellipticity, provides a symmetry-based check that a measured signal is a bulk photogalvanic current rather than heating or injection noise.","The sign reversals of the longitudinal current as ellipticity is varied offer a fingerprint of broken inversion symmetry and of interfering multi-photon pathways that open as the Dirac cone gaps out.","Because the effect needs only a monochromatic source, it extends BPG-based probing to simpler laser setups than the few-cycle or two-color schemes used previously."],"supporting_citations":[{"why":"Establishes the baseline that gapped non-centrosymmetric crystals support bulk photogalvanic currents under monochromatic illumination.","marker":"[62]"},{"why":"Provides the experimental precedent for photogalvanic currents in non-centrosymmetric systems.","marker":"[63]"},{"why":"Supplies the Floquet-dressed-state picture used to interpret the ring-shaped conduction-band occupations that determine current direction.","marker":"[21]"},{"why":"Sets the symmetry conditions for photocurrents, used to explain why linear and circular driving yield zero net current while elliptical driving does not.","marker":"[61]"},{"why":"Earlier demonstration of BPG currents with tailored polychromatic pulses; this paper's monochromatic-elliptical result extends that mechanism.","marker":"[38]"},{"why":"Connects photocurrent and Hall response to Berry curvature and valley occupations, supporting the transverse-current interpretation.","marker":"[57]"},{"why":"Real-space time-dependent density functional theory implementation used for the ab initio hBN validation.","marker":"[77]"},{"why":"Continued description of the same TDDFT implementation, supporting reproducibility of the ab initio results.","marker":"[78]"}],"fun_headline_variants":["Elliptical light unlocks tunable photocurrents in 2D crystals","Photocurrent control via light ellipticity in 2D materials","Gap sensing through ellipticity-driven photocurrents","Elliptical pulses steer currents and probe gaps in 2D","Light shape governs photocurrents in 2D hexagonal lattices"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proposed gap spectroscopy assumes that the simplified Gaussian-potential model and the independent-particle/adiabatic density-functional dynamics faithfully reproduce the real material's nonlinear photocurrent at every gap size, even though the model and ab initio calculations already disagree about which crystal orientation gives the maximum current in monolayer hBN.","fun_headline_variants_meta":{"raw":{"variants":["Elliptical light unlocks tunable photocurrents in 2D crystals","Photocurrent control via light ellipticity in 2D materials","Gap sensing through ellipticity-driven photocurrents","Elliptical pulses steer currents and probe gaps in 2D","Light shape governs photocurrents in 2D hexagonal lattices"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00086,"raw_usage":{"total_tokens":3762,"prompt_tokens":1003,"completion_tokens":2759,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":619,"completion_tokens_details":{"reasoning_tokens":2666}},"tokens_in":619,"tokens_out":2759,"duration_ms":17534,"temperature":1.0,"reasoning_tokens":2666,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:53:06.748267+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the ellipticity-dependent photocurrent in a monolayer hBN sample under 800 nm driving, comparing the current for the major axis along zigzag and along armchair. The model predicts the armchair orientation carries the maximum current, while the ab initio calculation predicts zigzag; if experiment follows the ab initio result, the model-based extrapolation to other gap sizes is not validated. A stronger test is to tune the gap continuously by strain or alloying and check whether the small-gap linear scaling and sign-change pattern predicted by the model appear in the measured current.","supporting_citations":[{"cited_title":"Nonperturbative Nonlinear Transport in a Floquet-Weyl Semimetal","cited_arxiv_id":"2409.04531","evidence_quote":"Provides the experimental precedent for photogalvanic currents in non-centrosymmetric systems."},{"cited_title":"Han, APL Materials 4, 032401 (2016), https://pubs.aip.org/aip/apm/article- pdf/doi/10.1063/1.4941712/19947471/032401 1 1.4941712.pdf","cited_arxiv_id":null,"evidence_quote":"Supplies the Floquet-dressed-state picture used to interpret the ring-shaped conduction-band occupations that determine current direction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Sets the symmetry conditions for photocurrents, used to explain why linear and circular driving yield zero net current while elliptical driving does not."},{"cited_title":"Quantum interference and occupation control in high harmonic generation from monolayer $WS_2$","cited_arxiv_id":"2503.04335","evidence_quote":"Earlier demonstration of BPG currents with tailored polychromatic pulses; this paper's monochromatic-elliptical result extends that mechanism."},{"cited_title":"Bauer and K","cited_arxiv_id":null,"evidence_quote":"Connects photocurrent and Hall response to Berry curvature and valley occupations, supporting the transverse-current interpretation."},{"cited_title":"Broers and L","cited_arxiv_id":null,"evidence_quote":"Continued description of the same TDDFT implementation, supporting reproducibility of the ab initio results."}],"review_version":1}