{"id":"f6598dbc-2a7c-4c1c-b015-dbe38da010a1","arxiv_id":"2412.13329","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A mean-field model with a tiny anisotropy of the CDW order parameter reproduces the measured polarization-dependent high harmonic spectra of TiSe2, proposing an all-optical probe of CDW order.","lead":"This paper builds a simplified quantum model of the charge-density-wave material TiSe2 and shows that the angle-dependent high harmonic light emission seen in a recent experiment can be reproduced when the ordered state is slightly anisotropic. It suggests harmonic polarimetry could become an all-optical probe of internal order in such materials.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Low-temperature inference rests on unmeasured 0.1 meV CDW anisotropy plus an undisclosed dephasing time; the same asymmetry may be reproducible by other symmetry-breaking terms, so the claim is underdetermined.","rationale":"The reader's weakest assumption is the same load-bearing concern: the anisotropic CDW order parameters are chosen to match the target data and are never independently measured. My read sharpens this into a parameter-uniqueness problem: the fit is underdetermined not only by the unmeasured DeltaQi values but also by the undisclosed dephasing time and the absence of any computed alternative symmetry-breaking mechanism. The high-temperature mechanism and the 45-degree projection argument are credible and well supported by the figures, so a rejection is not warranted. The requested changes from the conditional verdict are correct: disclose tau, provide code and data, quantify the fit, and either measure the CDW anisotropy independently or reframe the low-temperature claim as a consistency check rather than a unique revelation.","tokens_in":14762,"tokens_out":8443,"duration_ms":89302,"concrete_test":"Run the same mean-field simulation and reproduce Fig. 4 with isotropic DeltaQ but add one alternative symmetry-breaking term, such as a uniaxial strain correction delta_t to the tight-binding hopping t(pd_sigma) along the Q1 direction, or a two-domain population weight p, plus a scanned dephasing time tau from 5 fs to 100 fs. Optimize the extra parameters against the experimental H3/H5/H7 angular curves including error bars and report normalized residuals. If any strain-only or domain-only model fits the low-temperature pattern as well as the DeltaQ-anisotropy model, the central attribution fails; if the anisotropy model is uniquely best across the scanned range, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section V and Fig. 4 carry the paper's central claim: the low-temperature H5/H7 asymmetry is attributed to static anisotropy DeltaQ1=0.1148 eV, DeltaQ2=0.1147 eV, DeltaQ3=0.1146 eV. These values are not independently measured, and their smallness (about 0.1 meV, roughly 0.09%) makes them difficult to distinguish from other symmetry-breaking perturbations. The paper does not report the dephasing time tau introduced after Eq. (4), a free parameter known to control harmonic intensities and angular distributions. The theoretical curves are normalized per harmonic and no goodness-of-fit is reported, so the claim of an excellent match is not quantified. The authors assert that crystal-axis offset alone cannot explain the sign of the low-temperature asymmetry and that disorder, residual strain, or surface inhomogeneity would behave similarly in both phases, but they do not compute these alternatives. They also leave underspecified whether DeltaQi(t) in Eq. (4) is updated self-consistently during the pulse or frozen at the initial values; if frozen, the model assumes the anisotropy rather than revealing it. A uniaxial strain correction to the tight-binding hoppings, a domain population imbalance among the three Q states, or an intensity-dependent alignment shift could in principle generate the same left-peak-higher H5 pattern. The experiment may reveal low-temperature symmetry breaking with the Q1>Q3 orientation, but the specific conclusion that the CDW order parameters themselves are anisotropic is not uniquely established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a phenomenological mean-field tight-binding model for the CDW phase of TiSe2 and computes high-harmonic generation (HHG) spectra by solving the time-dependent Schrödinger equation in the velocity gauge. The authors show that the experimentally observed polarization dependence of H3, H5, and H7 in the high-temperature (semimetal) phase is reproduced when the actual 45-degree incidence geometry is used, with the H5 double peak emerging from the projection of the laser field onto the material plane. For the low-temperature CDW phase, they reproduce the asymmetric H5/H7 polarization traces by introducing a small static anisotropy among the three CDW order parameters, ΔQ1 = 0.1148 eV, ΔQ2 = 0.1147 eV, and ΔQ3 = 0.1146 eV, together with a crystal-axis offset. The central claim is that this agreement reveals an anisotropic CDW order in the low-temperature phase.","tokens_in":15037,"tokens_out":3982,"duration_ms":41474,"significance":"The high-temperature part of the paper is a valuable and robust result: the simplification to a two-dimensional tight-binding model together with the 45-degree incidence geometry explains the qualitative peak structure of H3, H5, and H7 without ad hoc symmetry breaking, and Fig. 5 shows that this conclusion is stable over a range of laser frequencies and amplitudes. If the low-temperature inference could be strengthened, the paper would establish HHG polarimetry as a genuinely informative probe of CDW anisotropy. However, as it stands, the low-temperature claim is not yet established: the specific anisotropy of the CDW order parameters is inserted by hand to match the data rather than extracted from it, and the dephasing time that controls the harmonic intensities is undisclosed. The paper is therefore a useful contribution with a promising framework, but its flagship conclusion requires further work before it can be accepted as stated.","major_comments":[{"comment":"The central low-temperature conclusion is obtained by choosing the CDW order-parameter anisotropy to reproduce the experiment. The green curves in Fig. 4 are generated with ΔQ1 = 0.1148 eV, ΔQ2 = 0.1147 eV, and ΔQ3 = 0.1146 eV, and the corresponding interaction strengths UQ1 = -1.80965 eV, UQ2 = -1.80935 eV, UQ3 = -1.80905 eV; the text states that this choice is made so that the left H5 peak becomes higher than the right peak. Since these values are not independently measured and no inversion or fitting residual is reported, the statement that the model 'reveals' an anisotropic CDW order is circular: the asymmetry is an input, not an output. Please either perform a quantitative parameter scan and show that the experimental data uniquely select this small anisotropy within the model, or reframe the claim as 'can be reproduced by' and discuss the degeneracy with other symmetry-breaking mechanisms.","section":"Section V, Fig. 4"},{"comment":"The dephasing time tau, introduced in the phenomenological dephasing step after Eq. (4), is never specified anywhere in the manuscript. Harmonic intensities and their angular distributions are strongly sensitive to tau, so omitting its value makes the theoretical curves in Figs. 3, 4, 5, and 6 unreproducible. Please report the value of tau used for every calculation, and ideally show a short convergence check with respect to tau.","section":"Section III, after Eq. (4)"},{"comment":"It is unclear whether the CDW order parameters ΔQi(t) in Eq. (4) are updated self-consistently during the laser pulse or frozen at their initial self-consistent values. The text says the Hamiltonian depends on the 'time-evolved' order parameter and gives the self-consistency expression, but the numerical implementation is not described. If the order parameters are frozen, then the low-temperature anisotropy is imposed at t=0 and the calculation does not demonstrate that the CDW order responds to or is revealed by the laser field. Please specify the update procedure and state explicitly whether ΔQi(t) is evolved in the simulations shown.","section":"Section III, Eq. (4) and Section V"},{"comment":"The paper dismisses alternative symmetry-breaking sources for the low-temperature H5 asymmetry—such as disorder, residual strain, surface inhomogeneity, domain population imbalance, or multiband effects—without computing them. Since the fitted anisotropy is only about 0.1 meV, many small perturbations could in principle produce the same H5 peak-height reversal, and the authors' argument that the high-temperature data constrain these effects is qualitative rather than quantitative. To support the uniqueness of the CDW-anisotropy interpretation, please provide at least one concrete counter-check, such as a calculation with a uniaxial strain correction to the tight-binding hoppings or a domain-population imbalance, showing that these alternatives do not reproduce the observed low-temperature asymmetry.","section":"Section V, low-temperature discussion"}],"minor_comments":[{"comment":"There is a typo in the text: 'paramter' should be 'parameter' in the sentence introducing the mean-field order parameter.","section":"Section II"},{"comment":"The caption does not identify which shaded region corresponds to experimental error bars; please label the experimental data and the error bars explicitly in the figure.","section":"Section V, Fig. 3 caption"},{"comment":"The experimental offset angle is quoted as 7° ± 2°, but the simulations use α = 5°. A brief comment on why 5° was chosen and how the result depends on α would clarify the fit.","section":"Section V, α discussion"},{"comment":"The conclusion says the anisotropic CDW 'could be induced by the possible onset of strain,' but the body text offers no calculation or estimate for the strain magnitude needed to produce the 0.1 meV anisotropy. Adding such an estimate would help the reader judge the plausibility of the proposed mechanism.","section":"Section VI"},{"comment":"The word 'stability' in the appendix is used loosely: Fig. 5 demonstrates robustness with respect to laser parameters, and Fig. 6 demonstrates persistence of the asymmetry when the anisotropy is increased, but no quantitative criterion for 'stable' is given. Please define what level of variation is considered acceptable.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the good news. The paper's high-temperature story is convincing and useful. The authors show that the measured HHG polarization pattern—single-peaked H3 and H7, double-peaked H5—comes from the 45-degree incidence geometry projecting the laser field onto the sample plane, not from any exotic CDW physics. Fig. 5 backs this up by showing the peak structure is stable across laser frequencies and amplitudes. That is a genuine insight and will be referenced.\n\nThe low-temperature conclusion is another story. The claim that HHG reveals an anisotropic CDW order rests on choosing the CDW order parameters to be slightly anisotropic (DeltaQ1=0.1148, DeltaQ2=0.1147, DeltaQ3=0.1146 eV). The authors are transparent about this—they say they 'chose' these values to reproduce the left/right peak behavior—but that means the agreement in Fig. 4 is a fit, not a prediction. The differences are ~0.1 meV, about 0.09% of the order parameter. At that level, any symmetry-breaking perturbation (uniaxial strain, domain imbalance, laser-induced alignment shift) could generate the same asymmetry. The paper does not compute those alternatives, so the specific attribution to static CDW anisotropy is underdetermined.\n\nThere is also a reproducibility issue: the dephasing time tau introduced after Eq. (4) is never specified. Dephasing is known to change harmonic intensities and angular distributions; omitting it leaves the calculation incomplete. And there is no goodness-of-fit statistic; the curves are normalized per harmonic, so 'excellent match' is visual.\n\nThe paper also notes its own limitations (the paragraph after Eq. 3 acknowledges the tight-binding caveats and multi-band effects), which is honest. I am not saying the low-T mechanism is wrong—it is plausible. But it is a consistency check, not a revelation.\n\nWho is this for? People working on HHG as a probe of phase transitions and the CDW community. They will get value from the high-T geometry mechanism and from the cautionary example of how easy it is to fit HHG angular data.\n\nRecommendation: This deserves peer review. It is technically sound, well-structured, and the high-T result is a real contribution. A referee should push for (1) reporting tau, (2) quantifying the fit, (3) either measuring the anisotropy independently or reframing the low-T claim as a consistency check that motivates future experiments, and (4) testing at least one alternative symmetry-breaking mechanism. With those revisions, the paper would be publishable.","headline":"Solid high-temperature geometry mechanism; the low-temperature anisotropy claim is a fitted input, not a revealed property.","tokens_in":15673,"tokens_out":2773,"would_cite":true,"duration_ms":23332,"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":"This paper claims that a tiny anisotropy in the three charge-density-wave order parameters of TiSe2 explains the asymmetric high-harmonic spectra observed at low temperature.","keywords":["TiSe2","charge-density-wave","high harmonic generation","polarization-resolved HHG","mean-field model","time-dependent Schrödinger equation","triple-Q CDW","anisotropy"],"falsifier":"Directly measure the three CDW order parameters at low temperature, for example by X-ray diffraction or scanning tunneling microscopy, and check whether $\\Delta Q_1$, $\\Delta Q_2$, and $\\Delta Q_3$ differ by roughly 0.1 percent in the direction the model requires; alternatively, apply a controlled uniaxial strain and test whether the H5 peak-height asymmetry follows the predicted relation between strain and harmonic polarization curves.","tokens_in":14516,"feed_emoji":"💡","tokens_out":6461,"duration_ms":55125,"temperature":0.7,"pith_summary":"This paper argues that the polarization-dependent high-harmonic generation (HHG) spectra of the layered material TiSe2 are explained by a simplified phenomenological mean-field model of its charge-density-wave (CDW) phase. Solving the time-dependent Schrödinger equation for an 8×8 Bloch Hamiltonian, the model reproduces the measured single-peaked third- and seventh-harmonic and double-peaked fifth-harmonic polarization curves. At low temperature, the strongly asymmetric spectra are captured only when the three CDW order parameters differ by about 0.1 percent, a weak anisotropy attributed to strain or disorder. If this is right, HHG polarimetry becomes a sensitive all-optical probe of CDW anisotropy.","feed_headline":"0.1% CDW anisotropy explains TiSe2 harmonic asymmetry","feed_subtitle":"A 0.1% spread in the three CDW order parameters explains the asymmetric H3, H5, H7 polarization curves.","key_machinery":"The machinery is an 8×8 momentum-dependent mean-field Hamiltonian in the reduced Brillouin zone, built from a two-band tight-binding model of TiSe2 and three CDW order parameters $\\Delta Q_1$, $\\Delta Q_2$, $\\Delta Q_3$. The harmonic spectrum is computed by integrating the velocity-gauge density-matrix equation $i\\hbar\\, d\\rho/dt = [H(k + eA(t)/\\hbar; \\{\\Delta Q_i(t)\\}), \\rho]$, with a dephasing step in the adiabatic basis. Two geometric elements do the explanatory work: the projection of the 45°-incident laser vector potential onto the crystal plane makes the in-plane field strength vary with polarization angle, and a small crystal-axis offset $\\alpha = 5^\\circ$ accounts for the weak high-temperature asymmetry. Adding the tiny CDW anisotropy then produces the strong low-temperature asymmetry.","core_discovery":"The central claim is that three ingredients together shape the HHG response of TiSe2: the hexagonal band structure, the 45-degree incidence geometry of the driving laser, and a weak anisotropy among the three CDW order parameters. In the low-temperature phase, the measured intensity distributions of H3, H5, and H7 as functions of polarization angle are matched by a mean-field solution with $\\Delta Q_1 = 0.1148\\,\\mathrm{eV}$, $\\Delta Q_2 = 0.1147\\,\\mathrm{eV}$, and $\\Delta Q_3 = 0.1146\\,\\mathrm{eV}$. This one-part-in-a-thousand breaking of the three-fold CDW symmetry tilts the harmonic peak heights in just the way the experiment shows, while the single- versus double-peaked structures are produced by projecting the 45°-incident field onto the material plane. The paper therefore attributes the low-temperature harmonic asymmetry to anisotropic CDW order.","pith_inferences":["Controlled uniaxial strain experiments on TiSe2 would provide a direct test: if strain is the source of the anisotropy, the H5 peak-height asymmetry should vary systematically with applied strain.","The apparent sensitivity to order-parameter differences of about 0.1% suggests that HHG polarimetry might also detect CDW domain patterns or inhomogeneous strain, since different regions would contribute different effective anisotropies.","An alternative mechanism, such as laser-induced dynamics of the order parameters or multi-band effects, would weaken the static-anisotropy conclusion; measuring the harmonic response as a function of pulse duration or intensity could distinguish static from dynamic symmetry breaking.","The mechanism of amplifying a tiny ground-state symmetry breaking through the nonlinear optical response may generalize to other correlated phases, where HHG asymmetry could act as a fingerprint of hidden order."],"forward_implications":["If the model is correct, polarization-resolved HHG can reveal anisotropy in CDW order parameters at the 0.1 percent level, a sensitivity that linear optical probes do not obviously provide.","The 45-degree incidence geometry is an essential part of the explanation: it converts the intrinsic three-peaked harmonic pattern into the single-peaked H3/H7 and double-peaked H5 curves seen experimentally.","The high-temperature H5 asymmetry arises from a crystal-axis offset of about 5 degrees, whereas the low-temperature asymmetry requires CDW anisotropy and cannot be explained by the offset alone.","The same simulation scheme could be applied to other triple-Q CDW materials, predicting harmonic polarization curves that would test whether their order parameters are similarly anisotropic."],"supporting_citations":[{"why":"Provides the experimental HHG polarization data that the model must reproduce.","marker":"[54]"},{"why":"Supplies the tight-binding band structure of TiSe2 used to construct the two-band dispersions.","marker":"[57]"},{"why":"Gives the mean-field temperature dependence of the CDW order parameters used in the model.","marker":"[58]"},{"why":"Supplies the velocity-gauge equation of motion for the density matrix in solids.","marker":"[59]"},{"why":"Provides the dephasing-in-adiabatic-basis scheme used in the HHG simulation.","marker":"[32]"},{"why":"Is the Slater-Koster method used to derive the tight-binding transfer integrals.","marker":"[60]"}],"fun_headline_variants":["Tiny CDW anisotropy flips TiSe2 harmonic asymmetry","0.1% order-parameter spread drives TiSe2 HHG asymmetry","Anisotropic CDW order revealed by TiSe2 high harmonics","TiSe2's 0.1% CDW anisotropy shapes harmonic emission","Mean-field model matches TiSe2 harmonic polarization asymmetry"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central claim rests on the assumption that the real sample's three CDW order parameters are statically anisotropic by the specific tiny amounts chosen, with the anisotropy attributed to strain or disorder but never independently measured.","fun_headline_variants_meta":{"raw":{"variants":["Tiny CDW anisotropy flips TiSe2 harmonic asymmetry","0.1% order-parameter spread drives TiSe2 HHG asymmetry","Anisotropic CDW order revealed by TiSe2 high harmonics","TiSe2's 0.1% CDW anisotropy shapes harmonic emission","Mean-field model matches TiSe2 harmonic polarization asymmetry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000474,"raw_usage":{"total_tokens":2362,"prompt_tokens":960,"completion_tokens":1402,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":576,"completion_tokens_details":{"reasoning_tokens":1311}},"tokens_in":576,"tokens_out":1402,"duration_ms":10135,"temperature":1.0,"reasoning_tokens":1311,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T13:14:18.401644+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Directly measure the three CDW order parameters at low temperature, for example by X-ray diffraction or scanning tunneling microscopy, and check whether $\\Delta Q_1$, $\\Delta Q_2$, and $\\Delta Q_3$ differ by roughly 0.1 percent in the direction the model requires; alternatively, apply a controlled uniaxial strain and test whether the H5 peak-height asymmetry follows the predicted relation between strain and harmonic polarization curves.","supporting_citations":[{"cited_title":"Tyulnev, L","cited_arxiv_id":null,"evidence_quote":"Provides the experimental HHG polarization data that the model must reproduce."},{"cited_title":"Kaneko, Y","cited_arxiv_id":null,"evidence_quote":"Supplies the tight-binding band structure of TiSe2 used to construct the two-band dispersions."},{"cited_title":"Monney, C","cited_arxiv_id":null,"evidence_quote":"Gives the mean-field temperature dependence of the CDW order parameters used in the model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the velocity-gauge equation of motion for the density matrix in solids."}],"review_version":1}