{"id":"40f1d610-6163-4f44-a77d-ba015b750b85","arxiv_id":"2511.11493","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Bulk PdTe2 shows strong SHG, four-wave mixing, and THz harmonic generation despite its centrosymmetric structure, with visible SHG peaking near a 2.9 eV resonance.","lead":"PdTe2 crystals show strong second- and third-order nonlinear optical responses in both visible and terahertz light, including harmonic generation and frequency mixing. The work suggests PdTe2 could be used for THz rectification, frequency mixing, and beam focusing, though the THz analysis relies on fits.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"THz χ2/χ3 extraction is ill-posed per paper's own collinearity admission; abstract overclaims fingerprints.","rationale":"The reader's weakest assumption precisely pinpoints the ill-posed fit in Eq. (4). The paper's own admission that quadratic and cubic terms become nearly collinear means the fit cannot separate χ2 and χ3, so the abstract's claim to have extracted fingerprints of both is not supported. However, the visible SHG/FWM data are compelling, and the observation of THz spectral gain itself is a real effect. Thus the conditional verdict—requiring revision of THz claims or additional data—is appropriate. My proposed re-fitting test would decisively reveal whether the existing data can actually distinguish the two orders.","tokens_in":13300,"tokens_out":5017,"duration_ms":41689,"concrete_test":"Re-analyze the raw power-dependence data at a representative frequency in the gain band (e.g., 2.1 THz) by fitting three nested models: linear+quadratic, linear+cubic, and full linear+quadratic+cubic. Compute AIC/BIC and parameter covariance. If the nested single-higher-order models fit comparably and the full model yields near-degenerate, large-uncertainty χ2 and χ3, then the data cannot discriminate the two orders; the extraction claim must be withdrawn. If the raw data are not available, repeat the THz power scan over at least a decade of pulse energies and/or use a narrowband (multi-cycle) source to spectrally resolve 2f and 3f components and verify their power slopes are 2 and 3 respectively.","verdict_should_be":"UNCHANGED","load_bearing_attack":"We agree with the reader's weakest assumption. The central THz claim—that the power-dependence analysis 'extracts fingerprints of both second- and third-order processes'—rests on Eq. (4), |E_rad|^2 = |χ1 P^{1/2} + χ2 P + χ3 P^{3/2}|^2. The paper itself concedes (Sec. IV) that 'although quadratic and cubic terms are mathematically independent, over the limited available power range with dominant linear response they become nearly collinear. Consequently, the fit converges to the dominant higher-order term rather than distributing weight across both.' This is an identifiability failure: within the measured power range, the basis functions P and P^{3/2} (and their cross terms with P^{1/2}) are nearly linearly dependent, so χ2 and χ3 cannot be uniquely determined. No error bars or confidence intervals are reported for the fitted coefficients, and the subsequent comparison to the photocurrent theory (Eq. 6, Appendix B) is a qualitative look-up after selecting the injection-current lifetime (τ_inj ≈ 100 fs) to match the data. Therefore the attribution of the 1.5–2.7 THz gain to THG and the near-DC excess to a second-order process is not established by the fitting; it is an assumption. The visible SHG/FWM results are well-supported, but the abstract's claim to have extracted fingerprints of both THz orders is overreaching.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports nonlinear optical measurements on single-crystal PdTe2 in two spectral regimes. In the visible, it presents SHG with quadratic power dependence and C3v-compatible six-fold polarization symmetry, with a resonant enhancement near 2.9 eV attributed to topological surface states; it also presents degenerate four-wave mixing with the expected quadratic-in-pump/linear-in-probe scaling. In the THz regime, it reports reflected-pulse spectra showing spectral gain in certain bands and fits the power dependence to a sum of linear, second-order, and third-order terms (Eq. 4). The authors then interpret the THz response within a radiative-photocurrent framework, using injection/shift-current kernels with adjustable lifetimes, and conclude that they have extracted fingerprints of both second- and third-order THz processes, with an injection-current lifetime of about 100 fs.","tokens_in":13731,"tokens_out":3575,"duration_ms":33864,"significance":"If the THz claims were fully supported, the paper would be a valuable demonstration of broadband second- and third-order nonlinear activity in a centrosymmetric type-II Dirac semimetal, with potential relevance for THz frequency mixing and rectification. The visible SHG/FWM measurements are a solid contribution: the power scalings and polarization patterns are internally consistent, the JDOS comparison is plausible, and the resonant-enhancement interpretation is well motivated. The paper also benefits from a clear presentation of the photocurrent framework. However, the central abstract claim to have 'extracted fingerprints of both second- and third-order processes in the THz regime' is not established by the current analysis. The paper itself admits that the quadratic and cubic terms are nearly collinear over the available power range, and no identifiability analysis is provided. The THz interpretation therefore needs substantial strengthening before the central claim can be accepted.","major_comments":[{"comment":"The central THz claim rests on the fit of |E_rad(f,P)|^2 = |chi1 P^{1/2} + chi2 P + chi3 P^{3/2}|^2. The authors explicitly concede in Sec. IV that 'over the limited available power range with dominant linear response they become nearly collinear. Consequently, the fit converges to the dominant higher-order term rather than distributing weight across both.' This is a model-identifiability failure. With only real-valued intensity data and three complex coefficients per frequency, the near-collinearity means chi2 and chi3 cannot be uniquely extracted. No error bars, confidence ellipses, or covariance estimates are reported, so the apparent separation in Fig. 5(b) may be an artifact of the fitting procedure. To support the claim of extracting both second- and third-order fingerprints, the authors need to demonstrate simultaneous identifiability, e.g., by reporting parameter uncertainties, p","section":"Sec. IV, Eq. (4), Fig. 5"},{"comment":"The interpretation of the near-DC second-order response as an injection current with tau_inj ~ 100 fs is circular. The model contains free parameters (tau_inj, tau_sh, sigma_inj, sigma_sh), and the authors select tau ~ 100 fs because it 'more closely matches the data in Fig. 5' — the same power-dependence data used to identify the nonlinear orders. This does not independently validate the injection-current mechanism; it only shows that the model family is flexible enough to reproduce the measured spectra. An independent test is required, such as helicity-dependent photocurrent measurements, a different pump-pulse duration, or direct time-resolved current detection. The conclusion's statement that the authors 'extracted both second and third order NLO processes' is therefore overstated.","section":"Sec. IV, Fig. 6 and Appendix B/C"},{"comment":"Eq. (6) assumes the nonlinear response kernels are constant over the incoming THz spectrum. However, the later analysis in Fig. 6 introduces strongly frequency- and lifetime-dependent kernels for injection and shift currents. Using the constant-kernel approximation to compute 'theoretical' second- and third-order spectra, then comparing them qualitatively with the power-dependence results, conflates two different models. In addition, the atmospheric-transmission correction is applied to the theoretical spectra but not consistently to the experimental spectra, making the comparison of peak positions (e.g., 1.75 THz vs. 1.5 THz) ambiguous. The authors should state explicitly which model is being tested and should apply the same transfer function to both experiment and theory.","section":"Sec. IV, Eq. (6)"}],"minor_comments":[{"comment":"The heading 'F our-W ave Mixing' contains a formatting typo; it should read 'Four-Wave Mixing'.","section":"Sec. III.B heading"},{"comment":"The caption says 'Se atoms shown in blue' but the compound is PdTe2; this should be Te atoms.","section":"Fig. 1 caption"},{"comment":"The sentence 'Figure 2(c) next presents the the spectra' has a duplicated 'the'.","section":"Sec. III.A"},{"comment":"The text refers to a 'SHG peak at ~1.75 eV' in the THz context; these are THz frequencies, so the label should be 'second-order response peak at ~1.75 THz' rather than 'SHG peak at ~1.75 eV'.","section":"Sec. IV, Fig. 6 discussion"},{"comment":"The phrase 'the reected and incoming pulses' is missing an 'f' in 'reflected'.","section":"Appendix C"}],"recommendation":"major_revision","confidential_remarks":"The visible SHG/FWM work is solid and likely publishable. The THz portion needs a major revision: either a rigorous identifiability analysis of Eq. (4) with uncertainties, or a substantial toning down of the claims to what the data actually support. A nested model comparison or an independent experimental control (e.g., changing pulse bandwidth/duration) would help. The citation of Faizanuddin et al. for surface SHG is appropriate, but the novel THz claim needs stronger statistical support."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The visible-laser part of this paper is genuinely good. The SHG power scaling is quadratic, the FWM scales as I_pump^2 * I_probe, and the polarization patterns match C3v surface symmetry. The new resonant SHG scan around 2.8–2.95 eV and its rough agreement with a JDOS calculation is a nice qualitative observation. The FWM result in PdTe2 is also new and cleanly executed. You can trust these data.\n\nThe THz reflection-geometry experiment is more interesting and more problematic. Measuring harmonic generation from a bulk crystal in reflection with a broadband pulse, in air, with water-vapor lines cut out, is not trivial. The raw observation—spectral gain above the reference in the 1.5–2.7 THz band and near DC—is credible and worth reporting.\n\nThe problem is the decomposition. Equation (4) fits |E_rad|^2 as a mixture of chi1 P^1/2, chi2 P, and chi3 P^3/2. The authors themselves say the quadratic and cubic terms are nearly collinear over the available power range, so the fit converges to one higher-order term rather than distributing weight across both. That is an identifiability failure. No error bars or confidence intervals are given for the fitted coefficients, so the separate chi2 and chi3 values in Fig. 5(b) are effectively placeholders. The later comparison to injection/shift-current models does not rescue this: they scan tau values, pick the one that looks like the data, and call it a lifetime of ~100 fs. That is a look-up, not a test.\n\nSo the abstract's claim to \"extract fingerprints of both second- and third-order processes in the THz regime\" is too strong. What is actually established is that there is nonlinear spectral gain, and that it is consistent with the presence of some combination of even- and odd-order nonlinearities. The specific assignment of the gain to THG and to a second-order injection current is not proven by this analysis.\n\nThat said, the paper does not hide the problem—the collinearity admission is right there in Section IV. The experimental observations are likely real, and the framework for analyzing stimulated THz emission via photocurrent kernels is reasonable. This is a paper that deserves serious refereeing, because the raw results are worth publishing and the analysis can be fixed by either narrowing the claims, adding narrowband THz data, or reporting proper uncertainties. I would send it to review, but the referee should insist on a rewrite of the THz sections and the abstract.\n\nFor a reading group, I'd probably skip it unless you want a concrete example of why power-dependence fits alone cannot separate nearly collinear nonlinear terms. I would not cite the THz extraction in my own work, but the visible FWM and resonant SHG data might be worth citing as a materials characterization point.","headline":"Solid visible NLO characterization of PdTe2, but the THz chi2/chi3 extraction is underdetermined by the authors' own admission—abstract overclaims.","tokens_in":860,"tokens_out":1850,"would_cite":true,"duration_ms":31595,"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":"The centrosymmetric semimetal PdTe2 shows strong second- and third-order nonlinear optical responses in both visible and terahertz light: second-harmonic generation peaks near 2.9 eV via topological surface states, while terahertz gain come","keywords":["PdTe2","topological Dirac semimetal","second-harmonic generation","terahertz spectroscopy","four-wave mixing","surface states","harmonic generation","photocurrent"],"falsifier":"Re-measure the THz power dependence with a narrowband or frequency-tunable source at a few discrete frequencies; if the radiated intensity in the 'gain' bands does not scale separately as P^2 and P^3, the claimed fingerprints of both second- and third-order processes would be falsified. Alternatively, a helicity-dependent photocurrent measurement that fails to show the predicted injection-current response would refute the near-DC mechanism.","tokens_in":13216,"feed_emoji":"⚡","tokens_out":9482,"duration_ms":71945,"temperature":0.7,"pith_summary":"PdTe2, a type-II Dirac semimetal with an inversion-symmetric crystal lattice, is shown to produce strong second- and third-order nonlinear optical responses in two very different frequency ranges: visible light and terahertz radiation. Visible second-harmonic generation follows the surface symmetry of the crystal and is resonantly enhanced near 2.9 eV, which the authors attribute to optical transitions between topological surface-state Dirac points. In the terahertz range, reflected spectra show gain at low and high frequencies, and a power-dependence analysis based on a radiative photocurrent model assigns the high-frequency gain to third-harmonic generation and the near-DC gain to a second-order process. If these claims hold, PdTe2 could function as a single material for terahertz frequency mixing, rectification, and beam focusing, offering a pathway to next-generation sensors and detectors.","feed_headline":"PdTe2 makes visible and terahertz harmonics despite symmetry","feed_subtitle":"Resonant visible harmonics and terahertz gain point to RF mixing and rectification in one material.","key_machinery":"The key analytical objects are (i) the radiative photocurrent model, which writes the emitted terahertz field as proportional to the time derivative of a current j = σ(1)*E + σ(2)*E^2 + σ(3)*E^3, and through the power-dependence relation |E_rad|^2 = |χ1 P^{1/2} + χ2 P + χ3 P^{3/2}|^2 yields linear, second-order, and third-order coefficients at each frequency; (ii) injection- and shift-current response kernels with finite lifetimes, used to reproduce the near-DC THz peak and identify the second-order mechanism as an injection current; and (iii) for the visible data, the computed joint density of states of PdTe2, which locates the SHG resonance near 2.9 eV, complemented by C3v surface point-gr","core_discovery":"On its own terms, the paper establishes that PdTe2 hosts measurable second- and third-order nonlinear optical responses in two widely separated spectral bands, despite its centrosymmetric 1T structure. Visible second-harmonic generation obeys the C3v surface point group, scales quadratically with pump power, and reaches maximum efficiency near 2.8–2.95 eV, matching a joint-density-of-states calculation and the energy separation between the topological surface-state Dirac points. Degenerate four-wave mixing confirms a third-order process with cubic power scaling and the expected polarization symmetry. In the terahertz regime, reflection measurements show the radiated spectrum exceeding the re","pith_inferences":["The acknowledged collinearity between the quadratic and cubic terms in the THz power fit makes the quantitative split between SHG and THG coefficients fragile; a narrowband THz source would test whether the gain bands independently scale as P^2 and P^3.","The JDOS match is a phenomenological proxy; a first-principles calculation of the full χ(2) tensor would determine whether the 2.9 eV resonance is genuinely driven by surface states or by bulk interband contributions.","If the near-DC response is an injection current, helicity-dependent photocurrent measurements—which the authors propose—could confirm the mechanism and measure its lifetime directly.","Since the THz spot is much larger than the visible one, the THz second-order signal may include contributions from defects and edges; comparing exfoliated monolayer or patterned samples would isolate the intrinsic surface response."],"forward_implications":["A single material could perform frequency upconversion and rectification across the visible-to-THz range, enabling compact THz sources and detectors without cryogenic or transmission-geometry constraints.","The SHG resonance tied to topological surface-state Dirac points suggests that electrostatic gating or thickness control could tune the nonlinear response, since the surface states can be modified.","THz reflection geometry, open to air, works for nonlinear spectroscopy of materials where transmission is impossible, expanding the range of candidate materials.","The near-DC second-order signal implies low-energy photocurrent generation, making PdTe2 a candidate for THz photodetection and nonlinear Hall effect studies.","Because THG is the sister process of fundamental enhancement, PdTe2 may be usable for self-focusing or refocusing of THz pulses."],"fun_headline_variants":["PdTe2 doubles and triples light from visible to terahertz","Topological metal PdTe2 mixes frequencies across two bands","PdTe2 emits harmonics despite centrosymmetry","Surface states give PdTe2 nonlinear response in visible and THz"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The THz assignment of spectral gain to third-harmonic generation versus a second-order process assumes the quadratic and cubic terms in the power-dependence fit are separable, but the paper notes they become nearly collinear with the dominant linear response over the available power range, so the fit can converge to one higher-order term rather than distributing weight across both.","fun_headline_variants_meta":{"raw":{"variants":["PdTe2 doubles and triples light from visible to terahertz","Topological metal PdTe2 mixes frequencies across two bands","PdTe2 emits harmonics despite centrosymmetry","Surface states give PdTe2 nonlinear response in visible and THz"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000242,"raw_usage":{"total_tokens":1341,"prompt_tokens":703,"completion_tokens":638,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":447,"completion_tokens_details":{"reasoning_tokens":568}},"tokens_in":447,"tokens_out":638,"duration_ms":6432,"temperature":1.0,"reasoning_tokens":568,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T22:09:15.581490+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-measure the THz power dependence with a narrowband or frequency-tunable source at a few discrete frequencies; if the radiated intensity in the 'gain' bands does not scale separately as P^2 and P^3, the claimed fingerprints of both second- and third-order processes would be falsified. Alternatively, a helicity-dependent photocurrent measurement that fails to show the predicted injection-current response would refute the near-DC mechanism.","supporting_citations":[],"review_version":1}