{"id":"af6452e5-f1c5-4aa0-a179-a170af99081e","arxiv_id":"2608.11471","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In KTaO3, oxygen-vacancy defect states produce the low-temperature SHG signal, and thermal fluctuations of the soft polar mode suppress that signal at higher temperatures, creating an effective inversion-symmetry crossover without a structural transition.","lead":"This paper explains why KTaO3, a crystal whose structure stays symmetric, still produces a strong nonlinear optical signal that appears only at low temperature. The authors argue that oxygen defects locally break symmetry while soft lattice vibrations scramble that signal at higher temperature, without any structural change.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The phase-only ansatz M_cvd(k;{Q})=M_cvd(k)e^{-iΣgνQν} is the load-bearing step: without it the soft-mode temperature dependence reduces to an uncontrolled fit. A frozen-phonon DFT check would settle it.","rationale":"The reader's weakest assumption identifies the same load-bearing point: the phase-only factorization before Eq. (4) and used in Eq. (5). I agree. The diffuse X-ray scattering and 3D-∆PDF analysis provide credible independent support for the absence of a structural transition and for phonon-dominated correlations, and the HSE06 calculation plausibly establishes a vacancy-mediated nonlinear channel. Those pieces support the qualitative picture of persistent local inversion-symmetry breaking inside a globally centrosymmetric crystal. The quantitative claim—that the SHG crossover is a soft-mode phase-coherence effect with a Debye-Waller-like e^{-W(T)}—rests entirely on the phase-only ansatz and a cumulant expansion whose derivation is deferred to an absent supplement. Since the coupling constants gν are undetermined, the successful fit of I(T) in Fig. 4(d) cannot by itself validate the mechanism; a monotonic suppression factor would also fit. The proposed frozen-phonon calculation is the minimal decisive test: if the matrix-element magnitude or the in-gap-state energy changes appreciably with soft-mode displacement, then the predicted crossover is not controlled by the soft mode alone, and the central claim would need revision. Thus I would keep the reader's CONDITIONAL verdict; no change is needed.","tokens_in":10138,"tokens_out":5536,"duration_ms":53759,"concrete_test":"Using the same HSE06 2×2×4 supercell, freeze in the soft TO phonon eigenvector at amplitudes Q corresponding to the thermal rms at T = 75 K (⟨Q²⟩^{1/2} from the soft-mode frequency), and directly compute the defect-assisted products ⟨c,k|d̂|d[Q]⟩⟨d[Q]|d̂|v,k⟩ and the in-gap-state energy as functions of Q. If the magnitude changes by more than a few percent, or if the defect-state shift is comparable to the detuning εc-εv-2ℏω or to the SHG linewidth, the phase-only parameterization fails and Eq. (5) is not justified. A clean pass would require |M(Q)|≈|M(0)| and |ΔE_d| much smaller than both the detuning and the thermal energy at the crossover.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central mechanism requires that the only temperature dependence of the defect-assisted vertex comes from a phase factor. The paper states, before Eq. (4) and again in reaching Eq. (5), that soft-mode displacements enter M_cvd(k;{Qν}) only as exp(-iΣν gνQν); all magnitude and energy-denominator effects are neglected. This factorization converts the Bose population of the soft mode into a Debye-Waller suppression e^{-W(T)}, and the good fit in Fig. 4(d) depends entirely on it. But the justification given—long-wavelength displacements are locally equivalent to translations—is not demonstrated for a localized oxygen vacancy: a polar TO displacement creates a local electric field at the defect, which should shift the in-gap state and the virtual-state denominator Λ(k), and may also change the dipole magnitudes. If those effects are comparable to the phase modulation, the predicted saturation temperature and the crossover shape are not fixed by the soft mode alone; they would be fit outcomes. The detailed derivation is explicitly deferred to Supplementary Section??, which is absent, so the factorization is currently asserted rather than shown.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a combined optical SHG, diffuse X-ray scattering, and hybrid-functional DFT study of KTaO3. The authors observe a strong SHG enhancement below roughly 75 K with saturation near 50 K, while diffuse scattering shows no static symmetry-lowering transition. They propose that charged oxygen vacancies (V_O^+) provide local inversion-symmetry breaking and a defect-assisted virtual-transition channel, and that thermal fluctuations of the soft polar mode suppress the coherent nonlinear vertex through a Debye-Waller-like factor exp[-W(T)]. The central claim is that the temperature-dependent SHG is a fluctuation-driven coherence crossover rather than an order-parameter transition, and that local inversion-symmetry-breaking defects persist at all temperatures while the crystal remains globally centrosymmetric.","tokens_in":10438,"tokens_out":4872,"duration_ms":42690,"significance":"If the proposed mechanism is validated, the paper would substantially revise the common polar-nano-region interpretation of SHG in quantum paraelectrics and would establish a general framework in which collective lattice fluctuations regulate the optical manifestation of local symmetry breaking. The manuscript has notable strengths: the diffuse-scattering analysis convincingly shows that the lattice stays dynamically disordered with no static transition; the HSE06 calculations identify a localized in-gap defect state with a finite nonlinear vertex; and the observed saturation temperature is plausibly linked to the soft-mode energy scale. These elements make the proposed picture credible, but the theoretical core currently rests on an unproven factorization and on coupling constants that are not derived from first principles.","major_comments":[{"comment":"The factorization M_cvd(k;{Qν}) = M_cvd(k) exp(-i Σν gν Qν), introduced just before Eq. (4) and used to obtain Eq. (5), is the load-bearing step of the manuscript: it is what converts soft-mode displacement fluctuations into the Debye-Waller suppression exp[-W(T)] that produces the agreement in Fig. 4(d). The justification given in the text, namely that long-wavelength lattice displacements are locally equivalent to translations, is not demonstrated for a localized oxygen vacancy; a polar TO displacement creates a local electric field that can shift the in-gap defect state and the virtual-state denominator Λ(k), and may also change the magnitude of M_cvd(k). The detailed derivation is deferred to 'Supplementary Section??', which is absent from the manuscript. I request either a full derivation of the phase-only form or a frozen-phonon DFT calculation of M_cvd(k;{Qν}) as a function of soft-mode amplitude that quantifies the neglected magnitude and energy-denominator effects.","section":"Microscopic mechanism, Eqs. (4)-(5)"},{"comment":"The effective defect-phonon coupling constants gν in Eq. (5) are not computed from the DFT wavefunctions, and the manuscript does not state whether they are fitted parameters. The dephasing rate γk(T) in Eq. (6) is also introduced phenomenologically. The temperature curve in Fig. 4(d) is therefore, as written, a fit with several free parameters rather than a parameter-free prediction of the mechanism. To support the central claim, the authors should report the values of gν and γk(T), specify which parameters are adjusted, and show the sensitivity of the predicted crossover shape and saturation temperature to those choices. Ideally, gν should be computed from the HSE06 electronic states and the soft-mode eigenvectors so that the central prediction is not a disguised fit.","section":"Microscopic mechanism, Eqs. (5)-(6) and Fig. 4(d)"},{"comment":"The claim that the SHG linewidth scales as T^3 is presented as evidence for acoustic-phonon-induced dephasing, but the Voigt analysis is not reported in sufficient detail: the Lorentzian component and the excitation laser bandwidth are not given, and the relation between γk(T) in Eq. (6) and the measured Gaussian FWHM is not stated. This matters because γk(T) in Eq. (6) also contributes to the temperature dependence of the SHG intensity and could be degenerate with exp[-W(T)] in reproducing the data in Fig. 4(d).","section":"Fig. 2(b) and Eq. (6)"}],"minor_comments":[{"comment":"There are unresolved placeholders, including 'Supplementary Section??', 'Supplement Material Sec. ??', and 'Fig. ??', which must be filled before submission.","section":"Throughout"},{"comment":"Reference 17 has a malformed DOI ('10.1103/74d5-4hsw'); please correct it.","section":"References"},{"comment":"The sentence 'and the report to reproduce key nonequilibrium SHG observations' appears to be a typographical error for 'and the failure to reproduce'; please rephrase.","section":"Introduction"},{"comment":"The color scale for the difference map between 26 K and 100 K is not defined; please add a scale bar or color bar so that the 3D-ΔPDF differences can be interpreted quantitatively.","section":"Fig. 3(d)"},{"comment":"The comparison of theory and experiment would benefit from error bars on the experimental data and a clear statement of which parameters enter the theoretical curve.","section":"Fig. 4(d)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a good fit for the journal and the experimental part appears sound. The main risk is that the theory section is currently a fit in disguise: the phase-only ansatz and the coupling constants are not derived. If the authors provide the requested frozen-phonon validation and parameter transparency, I would support publication. I do not see grounds for rejection because the missing derivation is within the scope of a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I think the paper is worth taking seriously, but it is not finished. The genuinely new thing is the mechanism: temperature-dependent SHG in KTaO3 is attributed to a defect-assisted nonlinear channel whose coherence is scrambled by thermal fluctuations of the soft polar mode, not to polar nanoregions or a ferroelectric transition. That is a real departure from the PNR literature, and the experiment supports the absence of a structural transition: the diffuse scattering shows no superlattice peaks, the 3D-ΔPDF stays phonon-like, and the onset-temperature invariance under annealing is a clean constraint. The HSE06 defect channel is also plausible.\n\nThe soft spot is the quantitative model. The Fig. 4(d) agreement relies on writing the defect-assisted matrix element as M(k)e^{-iΣgQ}, a phase-only factorization that turns Bose population into a Debye-Waller factor. That factorization is the load-bearing step, and the justification—long-wavelength displacements are locally equivalent to translations—is not obviously correct for a localized oxygen vacancy. A TO displacement produces a local electric field at the defect, which should shift the in-gap state and the virtual-state denominator. If those magnitude changes are comparable to the phase modulation, the predicted crossover shape and saturation temperature are not fixed by the soft mode alone. The derivation is deferred to Supplementary Section??, which is missing, so the claim is currently asserted rather than shown. A frozen-phonon DFT check of M_cvd as a function of the soft-mode coordinate would settle the point.\n\nThere are also many placeholders (Fig.??, Supplement Sec.??) and the dephasing rate is phenomenological. Those are addressable in revision. The citation pattern looks fine.\n\nBottom line: the mechanism is new and the experiment is coherent, but the central quantitative result is not yet demonstrated. I would send it to peer review—the physics matters enough to spend referee time on it—and would ask for the supplement, the derivation of the phase-only ansatz, and a frozen-phonon validation. Until then, treat Fig. 4(d) as illustrative.","headline":"A genuinely new mechanism for SHG in quantum paraelectrics, backed by good diffuse-scattering evidence, but the central quantitative result rests on an unproven phase-only ansatz and a missing supplement.","tokens_in":10952,"tokens_out":3121,"would_cite":false,"duration_ms":30089,"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":"In KTaO3, the SHG onset is a fluctuation-driven coherence crossover, not a structural transition.","keywords":["second-harmonic generation","quantum paraelectric","KTaO3","soft phonon","inversion symmetry","oxygen vacancy","Debye-Waller factor","diffuse X-ray scattering"],"falsifier":"A first-principles calculation of the full matrix element M_cvd(k;{Q_nu}) as a function of soft-mode displacement amplitude would settle it: if |M_cvd| varies appreciably while the phase wraps by order pi, the Debye-Waller factorization fails. Conversely, an experiment measuring SHG onset in isotopically substituted or strained KTaO3 with fixed defect density should show no onset shift if the mechanism is correct and a shift if soft-mode energy controls the crossover.","tokens_in":9942,"feed_emoji":"🔬","tokens_out":5887,"duration_ms":81516,"temperature":0.7,"pith_summary":"This paper tries to establish that the dramatic rise of second-harmonic generation (SHG) in the quantum paraelectric KTaO3 below 75 K is not the onset of a ferroelectric phase or cooperative polar order. Instead, the crystal stays globally centrosymmetric while oxygen-vacancy defects locally break inversion symmetry, and the temperature dependence is set by thermal fluctuations of the soft polar phonon, which scramble the phase coherence of the defect-assisted virtual transitions that generate SHG. The claim matters because it replaces the polar-nano-region and ferroelectric-ordering explanation for quantum paraelectrics with a purely dynamical, fluctuation-driven crossover, and because it predicts that defect-mediated nonlinear responses can be engineered by controlling low-energy phonons rather than by changing symmetry.","feed_headline":"Soft phonons, not a phase change, switch on KTaO3's SHG","feed_subtitle":"Cooling quenches soft-mode fluctuations, restoring coherent defect-mediated second-harmonic generation without a structural transition.","key_machinery":"The central object is the defect-assisted nonlinear optical vertex Mcvd(k;{Q_nu}) = Mcvd(k) exp(-i sum_nu g_nu Q_nu), the factorized form of the matrix element for a valence-to-defect-to-conduction virtual transition mediated by the oxygen-vacancy in-gap state. Combining this phase-only dependence on soft-mode normal coordinates with a Gaussian average (second-order cumulant expansion) produces the coherence factor $e^{{-W(T)}}$, with W(T) = (1/2) sum_nu |g_nu|^2 <$Q_nu^{2}$>_T; because the soft TO mode has tiny frequency, its Bose-enhanced fluctuations dominate W(T) and set the temperature scale. This factor, multiplied into the semiconductor Bloch-equation polarization, is what converts phonon fluctuations into a suppression of SHG and yields the paper's prediction for I_{2omega}(T).","core_discovery":"The central discovery is an effective inversion-symmetry crossover without a structural transition. In KTaO3, charged oxygen vacancies create a localized in-gap state that mediates a coherent two-step virtual transition between valence and conduction bands, giving a finite second-order nonlinear vertex even though the crystal is centrosymmetric. Because the defect is weakly pinned, soft-mode displacements only imprint phase factors on this matrix element; thermal averaging of harmonic soft-mode fluctuations yields a Debye-Waller-like factor $e^{{-W(T)}}$ that suppresses the nonlinear vertex at elevated temperatures. The measured SHG intensity is reproduced by this factor, with W(T) dominated by the anomalously low-frequency transverse optical soft mode, while diffuse X-ray scattering and 3D-PDF show only phonon fluctuations and no symmetry-lowering order. The paper concludes that symmetry and its manifestation need not be one-to-one: local inversion-symmetry breaking persists at all temperatures, and the crossover is a coherence phenomenon, not an order parameter.","pith_inferences":["Inference: If the phase-scrambling mechanism is generic, the same crossover should appear in SrTiO3 and other quantum paraelectrics, with onset temperature set by each material's soft-mode energy; this is testable by comparing SHG onsets across the family.","Inference: Isotope substitution or epitaxial strain that alters the soft-mode frequency should shift the SHG onset temperature even if the defect concentration is fixed, a signature that distinguishes this mechanism from defect-ordering models.","Inference: Time-resolved THz or mid-IR pumping that transiently heats the soft-mode bath should suppress defect-mediated SHG on picosecond timescales, offering a pump-probe test of the coherence-scrambling picture and a route to fast optical switching of inversion-asymmetry visibility.","Inference: The factorization assumption implies that the sum over soft-mode displacements enters only through the phase; a first-principles check of |Mcvd| versus Q would either confirm the model or reveal where it needs magnitude corrections."],"forward_implications":["No structural phase transition or cooperative polar ordering is required to explain the SHG onset in KTaO3; the 75 K scale is set by soft-mode thermal population, not by an order-parameter transition.","The SHG intensity becomes a quantitative thermometer of soft-mode displacement fluctuations: its saturation temperature tracks the soft-mode energy, and its magnitude is independently set by defect concentration.","Changing oxygen-vacancy concentration should alter the absolute SHG strength but not the onset temperature, exactly the decoupling observed under different annealing atmospheres.","Materials that combine local inversion-symmetry-breaking defects with low-energy polar phonons should show the same effective crossover, and strain, temperature, or coherent phonon excitation can tune the optical manifestation of broken local symmetry without any structural change."],"supporting_citations":[{"why":"Supplies the KTaO3 phonon dispersion and soft TO1 branch frequencies used to model diffuse-scattering intensity and W(T).","marker":"[33]"},{"why":"Documents the soft optic phonon and anomalous acoustic dispersion that establish the low-energy polar fluctuation spectrum.","marker":"[34]"},{"why":"Hybrid-DFT identification of the charged oxygen vacancy in-gap state that carries the defect-mediated SHG channel.","marker":"[36]"},{"why":"Shows oxygen vacancies modify KTaO3 electronic structure, grounding the weak pinning and susceptibility of the defect to polar fluctuations.","marker":"[44]"},{"why":"Provides the phase-modulation parameterization of the defect-phonon coupling used to derive the Debye-Waller factor.","marker":"[45]"},{"why":"Experimental soft-mode dispersion data that fixes the soft-mode energy versus temperature, the scale where SHG saturates.","marker":"[52]"},{"why":"Establishes the diffuse-scattering and 3D-PDF interpretation of phonon displacement correlations used to rule out a structural transition.","marker":"[31]"},{"why":"Demonstrates that nonequilibrium SHG in quantum paraelectrics can be reproduced without cooperative polar correlations, motivating the local-defect and soft-mode picture.","marker":"[17]"}],"fun_headline_variants":["Soft phonons hide local symmetry breaking in KTaO3","Thermal soft phonons turn KTaO3 effectively centrosymmetric","Soft-phonon jitter quenches SHG in KTaO3 without phase change","No symmetry change yet KTaO3 appears centrosymmetric","Soft phonons make KTaO3 look centrosymmetric"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that soft-mode displacements change only the phase of the defect-assisted optical matrix element, never its magnitude or the transition energies, so all temperature dependence enters through the $e^{{-W(T)}}$ coherence factor.","fun_headline_variants_meta":{"raw":{"variants":["Soft phonons hide local symmetry breaking in KTaO3","Thermal soft phonons turn KTaO3 effectively centrosymmetric","Soft-phonon jitter quenches SHG in KTaO3 without phase change","No symmetry change yet KTaO3 appears centrosymmetric","Soft phonons make KTaO3 look centrosymmetric"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001489,"raw_usage":{"total_tokens":5986,"prompt_tokens":963,"completion_tokens":5023,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":579,"completion_tokens_details":{"reasoning_tokens":4931}},"tokens_in":579,"tokens_out":5023,"duration_ms":30328,"temperature":1.0,"reasoning_tokens":4931,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:12:24.640348+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A first-principles calculation of the full matrix element M_cvd(k;{Q_nu}) as a function of soft-mode displacement amplitude would settle it: if |M_cvd| varies appreciably while the phase wraps by order pi, the Debye-Waller factorization fails. Conversely, an experiment measuring SHG onset in isotopically substituted or strained KTaO3 with fixed defect density should show no onset shift if the mechanism is correct and a shift if soft-mode energy controls the crossover.","supporting_citations":[{"cited_title":"H.et al.Phonon dispersion and lattice dynamics of KTaO3 from 4 to 1220 K.Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the KTaO3 phonon dispersion and soft TO1 branch frequencies used to model diffuse-scattering intensity and W(T)."},{"cited_title":"D., Harada, J","cited_arxiv_id":null,"evidence_quote":"Documents the soft optic phonon and anomalous acoustic dispersion that establish the low-energy polar fluctuation spectrum."},{"cited_title":"& Modak, B","cited_arxiv_id":null,"evidence_quote":"Hybrid-DFT identification of the charged oxygen vacancy in-gap state that carries the defect-mediated SHG channel."},{"cited_title":"K.et al.Oxygen vacancy induced electronic structure modification of ktao 3.Phys","cited_arxiv_id":null,"evidence_quote":"Shows oxygen vacancies modify KTaO3 electronic structure, grounding the weak pinning and susceptibility of the defect to polar fluctuations."},{"cited_title":"& Tanaka, K","cited_arxiv_id":null,"evidence_quote":"Experimental soft-mode dispersion data that fixes the soft-mode energy versus temperature, the scale where SHG saturates."},{"cited_title":"J., Talbayev, D","cited_arxiv_id":null,"evidence_quote":"Demonstrates that nonequilibrium SHG in quantum paraelectrics can be reproduced without cooperative polar correlations, motivating the local-defect and soft-mode picture."}],"review_version":1}