{"id":"38262f65-a4ef-4b5d-8d57-c31df1d1a03b","arxiv_id":"2506.10472","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In twinned LaAlO3, the 1.1 THz Eg Raman phonon is excited by a two-photon THz process, while the 0.86 and 0.36 THz oscillations are propagation artifacts, not phonons.","lead":"This paper combines two-dimensional terahertz Kerr effect spectroscopy with four-wave mixing simulations to show that, in LaAlO3 crystals, the 1.1 THz signal is a real Raman phonon while the 0.86 and 0.36 THz signals are light-propagation artifacts in a birefringent sample. The results give researchers a framework for separating genuine material response from spurious beats in nonlinear THz experiments.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The FWM simulation's spurious-peak assignment relies on unmeasured birefringence values (Δn_THz=0.06, Δn_pr=0.007) and a 20 fs probe; if these are tuned, the 0.86/0.36 THz attribution is not established.","rationale":"The paper's headline contribution is that two of the three observed THz Kerr features are propagation artifacts rather than material excitations. The direct experimental evidence for this—abrupt temporal cutoffs and the delay-dependent onset of the 0.86 and 0.36 THz signals—is suggestive but not by itself conclusive; the unambiguous identification comes from the FWM simulation claiming quantitative agreement. That simulation is the load-bearing pillar of the central claim. Its input parameters are not tightly constrained: the birefringence values are quoted only as 'on the order of previous estimates,' and the probe is modeled as a 20 fs pulse while the experimental probe is likely the 170 fs pulse used elsewhere in the setup. Without a sensitivity analysis or independent measurement, one cannot exclude the possibility that these values were chosen because they reproduce the target frequencies. The proposed independent birefringence measurement followed by a no-free-parameter rerun of the simulation would settle this point directly. The TPA mechanism for the 1.1 THz Eg mode is more robust: the field scaling in Fig. 2(c) and the 2D-TKE τ-dependence support it independently of the FWM parameters, and the absence of an IR-active mode in the pulse bandwidth strengthens it. Therefore the conditional verdict remains appropriate: the concern requires parameter anchoring but does not invalidate the core experimental observations. I agree with the reader that this is the weakest assumption, and I would not change the verdict.","tokens_in":8953,"tokens_out":9752,"duration_ms":125161,"concrete_test":"Rerun the FWM simulation of Sec. III.A using independently measured birefringence values for the specific 0.5 mm LAO crystal at 5 K and at 295 K (e.g., polarization-resolved THz time-domain spectroscopy for Δn_THz and 800 nm ellipsometry for Δn_pr) together with the actual experimental probe pulse duration, keeping all other model choices fixed. If the simulated 0.85/0.39 THz features shift by more than the ~0.05 THz experimental linewidth or disappear, the agreement in Fig. 3(f-g) is parameter-dependent rather than a robust prediction of the propagation-effect model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central attribution of the 0.86 THz and 0.36 THz features to propagation-induced spurious signals rests entirely on the FWM simulation in Sec. III.A and Appendix B. The simulation uses an instantaneous delta-function electronic response and equalizes all χ(3) tensor elements, but its output frequencies are set by the assumed birefringence values Δn_THz=0.06 and Δn_pr=0.007, which are only described as 'on the order of previous estimates' without quoting the actual prior values or their uncertainties. The reported agreement of 0.85 vs 0.86 THz and 0.39 vs 0.36 THz is close, but if these parameters were effectively fitted, the simulation cannot establish that the features are not genuine phonon modes. The same concern applies to the use of a 20 fs Gaussian probe pulse in the simulation while the experimental probe pulse is not explicitly stated and is likely the 170 fs pulse used for THz generation; probe duration changes the phase-matching convolution and could materially affect the simulated spectra. The paper provides no sensitivity scan over birefringence, no error bars on these parameters, and no comparison to an independent birefringence measurement of the same crystal, so the quantitative anchor for the central claim is missing.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports 2D-TKE measurements on a 0.5 mm [100]-cut LaAlO3 crystal at 5 K and in a companion room-temperature run. Three spectral features are observed at 1.1, 0.86, and 0.36 THz. The authors assign the 1.1 THz feature to the Eg Raman phonon driven quasi-instantaneously by two-photon THz excitation, supported by field scaling, the τ≈0 concentration of the 2D signal, and an equation-of-motion model. The 0.86 and 0.36 THz features are assigned to propagation-induced artifacts from co- and counter-propagating THz and 800 nm probe pulses in birefringent twin domains. The evidence includes short-time Fourier analysis showing that the 0.86 THz feature exists only for probe delays below t≈5 ps and the 0.36 THz feature between t≈5 and 17 ps, coinciding with calculated pulse-encounter times; Fabry-Perot sidebands; and an FWM simulation with instantaneous electronic hyperpolarizability that produces short-lived features at 0.85 and 0.39 THz. The Introduction states that the simulation quantitatively reproduces all novel features, while Section III.A and the Fig. 3 caption describe the agreement as qualitative.","tokens_in":9354,"tokens_out":7436,"duration_ms":94984,"significance":"If the central assignment survives scrutiny, the paper is valuable: it provides a cautionary and generalizable framework for separating genuine coherent phonons from propagation artifacts in birefringent crystals, and it offers a resolution to conflicting interpretations of the 1.1 THz excitation mechanism in LaAlO3. The strongest parts are the time-window and world-line arguments, which are largely independent of the FWM simulation parameters: an intrinsic phonon should persist across the full probe-delay window, whereas the observed features start and stop at precisely the calculated pulse-encounter times. The field-scaling data and the Fabry-Perot spacing provide additional support. The weaker part is the quantitative anchoring of the FWM simulation, which relies on unmeasured birefringence values and an idealized 20 fs probe pulse. With a sensitivity analysis or an independent birefringence measurement, the paper would be a solid contribution to nonlinear THz spectroscopy.","major_comments":[{"comment":"The FWM simulation's artifact frequencies (0.85 and 0.39 THz) are set by the assumed static birefringence values Δn_THz=0.06 and Δn_pr=0.007, which are described only as 'on the order of previous estimates' without quoting the prior values, their uncertainties, or an independent measurement on the same crystal. Since the simulation is invoked as confirmation that the 0.86 and 0.36 THz features are spurious, the agreement could be partly by construction. Please add a sensitivity analysis showing how the simulated peak frequencies and lifetimes depend on these two parameters, or provide a direct birefringence measurement; otherwise the quantitative anchor for the assignment is missing.","section":"III.A and Appendix B"},{"comment":"The simulation assumes a Fourier-limited 20 fs Gaussian probe pulse, but the experimental probe pulse duration is not stated anywhere. The only pulse duration given is the 170 fs, 800 nm pulse used for THz generation; if the probe is similar, the convolution and phase-matching effects in the FWM simulation could change materially. Please specify the measured probe duration and either justify the 20 fs assumption or repeat the key simulations with the actual probe envelope.","section":"Appendix B and Section II.A"},{"comment":"The Introduction claims that the FWM simulations 'quantitatively reproduce all of the novel features', but Section III.A states that 'the simulation results agree qualitatively with the experimental measurements' and the Fig. 3 caption uses 'qualitatively' as well. This distinction matters because the quantitative claim depends on the unvalidated birefringence and probe-duration assumptions discussed above. Please reconcile the wording and, if only qualitative agreement is claimed, soften the Introduction and abstract accordingly.","section":"Introduction and Section III.A"},{"comment":"The TPA model for the Eg phonon sets δ/M=1 and γ=0.005 and excludes back-reflection and counter-propagation, yet the paper concludes that TPA, rather than anharmonic coupling, is the excitation mechanism. The 2D-TKE data (strongest signal at τ=0 and linear scaling in both fields) are consistent with TPA, but a quantitative uniqueness test would require comparing the measured 2D pattern with at least one alternative anharmonic-coupling model, or showing that the τ≠0 echo-line delayed onset is reproduced by including counter-propagating fields in the Eq. (4) simulation. Please report the parameter sensitivity of the TPA simulations and provide the counter-propagation simulation or explicitly state that the echo-line attribution is only a qualitative argument.","section":"III.B and Eq. (4)"}],"minor_comments":[{"comment":"The material is described as 'twinned LaAlO3', but no twinning characterization (e.g., optical or diffraction evidence, twin-domain orientation) is presented anywhere in the text. Please clarify what is meant by twinned and how the twin domains are identified or modeled.","section":"Abstract and Introduction"},{"comment":"The derivation of the echo-line relation, Eq. (2), is referenced to Ref. [20], which is listed as 'To be submitted' and is therefore not currently accessible. Please provide the derivation in the main text or an appendix, or cite a published source.","section":"Appendix A"},{"comment":"The assumption that all allowed tensor elements of the third-order susceptibility have the same magnitude is a strong simplification for R-3c symmetry. Please state whether this simplification affects only amplitudes or could also affect the simulated frequencies or lifetimes.","section":"Appendix B"},{"comment":"The caption says the signal is sampled along the τ=0 line, but the text says the field scaling of ηNL is measured with the other THz pulse set to maximum. Please clarify how the τ=0 sampling is performed experimentally and whether the extracted scaling is along the pump-pump delay axis or along the probe axis.","section":"Fig. 2(c)"},{"comment":"There is a typo, 'THz-infared mixing', that should read 'THz-infrared mixing'.","section":"Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope and the time-window evidence is strong. The main revision is to provide sensitivity or independent measurements for the birefringence values and the probe-duration assumption, and to calibrate the 'quantitative' claims. I do not see grounds for rejection if these points are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short take: this paper deserves a serious referee. The experimental core is persuasive, and the simulation, despite rough parameters, is supportive rather than load-bearing.\n\nWhat's new: it reassigns the 0.86 and 0.36 THz peaks in LAO's THz Kerr signal from possible lattice modes to propagation artifacts, and it identifies the 1.1 THz Eg phonon as driven by two-photon absorption. The 2D-TKE data are clean, and the time-window behaviour is strong evidence on its own: the 0.86 THz feature disappears abruptly at exactly the delay when pump and probe meet at the rear surface, and the 0.36 THz feature appears at that moment and ends when the probe meets the back-reflected pump at the front surface. That's not something an intrinsic phonon would do. The echo-line structure for the 1.1 THz mode at nonzero pump delays is also nicely explained.\n\nThe weakest part is the FWM simulation. The birefringence values Δn_THz=0.06 and Δn_pr=0.007 are said to be 'on the order of previous estimates', not directly measured for this crystal, and no sensitivity scan is given. The simulation uses a 20 fs probe pulse while the experimental probe is presumably much longer. The paper calls the agreement 'quantitative' but Fig. 3 shows qualitative agreement at best.\n\nHowever, these issues do not vitiate the central claim. The cutoff times are set by the sample thickness and known refractive indices, not by the assumed birefringence. The simulation is a consistency check, not the primary proof. I would want the authors to add a sensitivity analysis, explain the probe duration choice, and either include the methods from the unpublished Ref. [20] or derive the echo-line equation in the appendix.\n\nThe paper will be useful to anyone doing nonlinear THz spectroscopy on birefringent crystals, and it provides a framework for spotting propagation artifacts that are easy to misread as phonons. I'd send it to peer review and ask for those revisions.","headline":"A solid 2D-TKE reassignment of LAO's low-frequency features as propagation artifacts; the temporal evidence carries the paper, and the simulation's rough parameters are a fixable soft spot.","tokens_in":9839,"tokens_out":3982,"would_cite":true,"duration_ms":47768,"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 that the 0.86 and 0.36 THz features in the THz Kerr response of twinned LaAlO3 are propagation artifacts, while the 1.1 THz Eg mode is excited by two-photon THz absorption.","keywords":["terahertz Kerr effect","LaAlO3","two-dimensional THz spectroscopy","four-wave mixing","birefringence","Raman phonon","two-photon absorption","propagation effects"],"falsifier":"Cut or obtain LAO samples of a different thickness, say 1.0 mm instead of 0.5 mm, and repeat the single-pulse TKE measurement: the propagation model predicts the 0.86 THz feature should end at the probe delay where THz and probe meet at the rear surface, which shifts with thickness, and the 0.36 THz feature should start at that same shifted time; if the onset and end times stay at the original roughly 5 ps and 17 ps, the propagation assignment fails. A second check is to rotate the probe polarization so the effective birefringence changes sign, which should alter the simulated artifact frequencies if the model is right.","tokens_in":8777,"feed_emoji":"⚡","tokens_out":6577,"duration_ms":65816,"temperature":0.7,"pith_summary":"The paper is trying to establish that apparent coherent phonon oscillations in the THz Kerr response of LaAlO3 must be interpreted with propagation effects included, and that doing so resolves earlier puzzles without invoking exotic couplings. Concretely, it claims the 1.1 THz Eg mode is driven by two-photon THz absorption, while the 0.86 THz and 0.36 THz features arise from co- and counter-propagation of THz and 800 nm probe pulses in birefringent twin domains, not from lattice dynamics. A reader should care because this gives a general way to distinguish genuine ionic oscillations from instantaneous electronic polarization in birefringent crystals, preventing false phonon and anharmonic-coupling assignments.","feed_headline":"Two THz 'phonon' peaks are propagation artifacts","feed_subtitle":"2D-TKE and FWM show the 0.86 and 0.36 THz signals come from pulse walk-off in birefringent twin domains.","key_machinery":"The central mechanism is a four-wave mixing simulation of the THz-THz-optical probe interaction in an extended birefringent medium: it propagates the two THz pump pulses and the 800 nm probe on a time-space grid, computes the local third-order polarization $P^{(3)}(t,z)$ from an instantaneous electronic response with all allowed $\\chi^{(3)}$ tensor elements set equal, and emits a new field that interferes with the probe in balanced detection. Static birefringence values $\\Delta n_{\\mathrm{THz}}=0.06$ and $\\Delta n_{\\mathrm{pr}}=0.007$ produce phase-matching oscillations near 0.85 and 0.39 THz whose lifetimes are set by pulse walk-off, matching the measured 0.86 and 0.36 THz features. The Eg phonon is modeled separately with a driven oscillator equation, $d^2 Q_R/dt^2 + 2\\gamma\\Omega_R\\, dQ_R/dt + \\Omega_R^2 Q_R = (\\delta/M)\\, E(t)^2$, representing two-photon absorption, which reproduces the 1.1 THz 2D-TKE pattern.","core_discovery":"On the paper's own terms, the central discovery is that the three features in the low-temperature THz Kerr response of twinned LaAlO3 come from two distinct mechanisms. The 1.1 THz mode is the Eg Raman phonon, excited quasi-instantaneously by two-photon THz absorption: 2D-TKE shows it peaks when the two THz pump pulses overlap in time, and its nonlinear amplitude scales linearly in each pump field. The 0.86 THz and 0.36 THz features are not lattice modes at all. They appear only over probe-delay windows bounded by the times at which the THz and probe pulses meet at the rear or front sample surfaces, and a four-wave mixing simulation that propagates both pulses through a birefringent medium with an instantaneous electronic hyperpolarizability reproduces them with no anharmonic phonon coupling and no unknown modes. Surface-reflected THz pulses produce additional 1.1 THz oscillations that can mimic anharmonic coupling, but they begin on a calculable echo line rather than at the arrival of the latest pump pulse.","pith_inferences":["Beyond the paper: the same four-wave mixing propagation model could be applied to previously reported THz Kerr peaks in other twinned oxide crystals to test whether some assigned low-energy modes are artifacts.","Beyond the paper: because the artifact frequencies are set by sample thickness and refractive indices, varying the sample thickness or probe wavelength should shift the 0.86 and 0.36 THz features in a predictable way, a testable prediction the paper does not spell out.","Beyond the paper: a polarization-resolved measurement that varies the angle between the THz polarization and the crystal c-axis could probe the equal-$\\chi^{(3)}$ assumption; if the tensor ratios differ, the simulated artifact frequencies and lifetimes would change, allowing the model's simplifying assumption to be checked."],"forward_implications":["In any birefringent crystal, apparent low-frequency coherent-phonon modes in THz Kerr data must be checked against co- and counter-propagation artifacts before an intrinsic lattice assignment is made.","The previously proposed anharmonic coupling between IR-active acoustic phonons and the Raman mode is not needed to explain the LAO data; two-photon absorption plus propagation effects account for all three observed features.","The 1.1 THz Eg mode can be treated as driven by $E(t)^2$, giving a simple route to modeling 2D-TKE signals in centrosymmetric insulators without invoking phonon-phonon coupling.","Time-domain cuts of 2D-TKE data can identify surface-reflection echoes by the echo-line condition $t + \\tau/2\\,(1 \\pm v/c) = L(1/v - 1/c)$, providing a diagnostic for spurious oscillations."],"supporting_citations":[{"why":"Provides the THz-THz-VIS four-wave mixing simulation scheme used to model propagation and birefringence artifacts in bulk samples.","marker":"[10]"},{"why":"The prior 2D-TKE study that attributed the 0.34 THz feature to parametric acoustic-Raman interaction; this paper's reinterpretation must beat that baseline.","marker":"[12]"},{"why":"Source of the birefringence estimates on the order of $\\Delta n_{\\mathrm{THz}}=0.06$ and $\\Delta n_{\\mathrm{pr}}=0.007$ used in the simulation.","marker":"[11]"},{"why":"Supplies the optical refractive index $n_{\\mathrm{pr}}\\approx 2$ at 800 nm used to compute pulse walk-off timing.","marker":"[16]"},{"why":"Explains the equidistant Fabry-Perot peaks observed on top of the 0.86 THz feature in the single-pulse spectrum.","marker":"[17]"},{"why":"Introduces the two-dimensional THz Kerr effect spectroscopy method used to separate excitation pathways by pump-pump delay.","marker":"[6]"},{"why":"Provides the two-photon absorption excitation model for Raman phonons in centrosymmetric media used in the equation of motion for the Eg mode.","marker":"[23]"},{"why":"Shows there is no IR-active phonon within the THz pump bandwidth, ruling out anharmonic IR-Raman coupling as the Eg excitation channel.","marker":"[25]"}],"fun_headline_variants":["Two low-THz peaks in LaAlO3 are propagation artifacts","Only 1.1 THz mode is real in LaAlO3 Kerr response","THz Kerr artifacts: 0.86 and 0.36 THz are false","Surface reflections mimic anharmonic THz coupling","LaAlO3 THz Kerr: two modes are pulse walk-off"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The spurious-signal assignment rests on the four-wave mixing simulation's assumed refractive-index difference values, about 0.06 for the THz pulse and 0.007 for the probe, and on taking all allowed nonlinear tensor components equal; if those numbers were effectively tuned to reproduce the 0.86 and 0.36 THz peaks, the agreement is partly by construction.","fun_headline_variants_meta":{"raw":{"variants":["Two low-THz peaks in LaAlO3 are propagation artifacts","Only 1.1 THz mode is real in LaAlO3 Kerr response","THz Kerr artifacts: 0.86 and 0.36 THz are false","Surface reflections mimic anharmonic THz coupling","LaAlO3 THz Kerr: two modes are pulse walk-off"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000305,"raw_usage":{"total_tokens":1776,"prompt_tokens":997,"completion_tokens":779,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":613,"completion_tokens_details":{"reasoning_tokens":683}},"tokens_in":613,"tokens_out":779,"duration_ms":8266,"temperature":1.0,"reasoning_tokens":683,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:25:13.257197+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Cut or obtain LAO samples of a different thickness, say 1.0 mm instead of 0.5 mm, and repeat the single-pulse TKE measurement: the propagation model predicts the 0.86 THz feature should end at the probe delay where THz and probe meet at the rear surface, which shifts with thickness, and the 0.36 THz feature should start at that same shifted time; if the onset and end times stay at the original roughly 5 ps and 17 ps, the propagation assignment fails. A second check is to rotate the probe polarization so the effective birefringence changes sign, which should alter the simulated artifact frequencies if the model is right.","supporting_citations":[{"cited_title":"Frenzel, M","cited_arxiv_id":null,"evidence_quote":"Provides the THz-THz-VIS four-wave mixing simulation scheme used to model propagation and birefringence artifacts in bulk samples."},{"cited_title":"Terahertz-driven parametric excitation of Raman-active phonons in LaAlO$_{3}$","cited_arxiv_id":"2410.06748","evidence_quote":"The prior 2D-TKE study that attributed the 0.34 THz feature to parametric acoustic-Raman interaction; this paper's reinterpretation must beat that baseline."},{"cited_title":"Kovalev, C","cited_arxiv_id":null,"evidence_quote":"Source of the birefringence estimates on the order of $\\Delta n_{\\mathrm{THz}}=0.06$ and $\\Delta n_{\\mathrm{pr}}=0.007$ used in the simulation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the optical refractive index $n_{\\mathrm{pr}}\\approx 2$ at 800 nm used to compute pulse walk-off timing."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Explains the equidistant Fabry-Perot peaks observed on top of the 0.86 THz feature in the single-pulse spectrum."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the two-dimensional THz Kerr effect spectroscopy method used to separate excitation pathways by pump-pump delay."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the two-photon absorption excitation model for Raman phonons in centrosymmetric media used in the equation of motion for the Eg mode."},{"cited_title":"Abrashev, A","cited_arxiv_id":null,"evidence_quote":"Shows there is no IR-active phonon within the THz pump bandwidth, ruling out anharmonic IR-Raman coupling as the Eg excitation channel."}],"review_version":1}