{"id":"aa567667-cb32-470e-884d-878365095782","arxiv_id":"2411.10160","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A unified many-body and Maxwell propagation theory explains how phonon-polaritons enter four-wave versus three-wave mixing experiments, predicting that broadband terahertz pumps access a broad region of the polariton dispersion in one shot.","lead":"This paper builds a theoretical framework for how light and lattice vibrations combine into phonon-polaritons during ultrafast laser experiments. It shows that broadband terahertz pumps can map these hybrid modes over a wider range of momenta than previously possible.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"TWM factorization (bare D0 kernel + dressed propagation) is asserted, not derived; if the correct kernel contains D_Q, the central dispersion-via-propagation claim fails.","rationale":"The reader's weakest_assumption identifies the same load-bearing concern I find: the paper's central claim depends on the factorization K(2)∝D_0 with all polariton dressing moved into the linear propagation of the THz field, and this factorization is not derived in the main text. The distinction between FWM and TWM rests entirely on which fields are integrated out, but the main text only states the result and defers the derivation to the missing Supplementary Material [47]. Because the physical phonon-photon coupling is the same in both protocols, it is not obvious why the TWM cubic vertex carries the bare phonon propagator rather than the dressed one; a fully contracted calculation would naturally produce the mixed propagator ⟨A Q⟩. This is precisely the kind of step that could be internally consistent but is not demonstrated. The proposed test settles it by comparing the factorized route with the fully contracted route; if they coincide, the concern is resolved and the paper's conditional status is appropriate. If they do not, the central physical conclusion requires revision. I therefore keep the reader's CONDITIONAL verdict unchanged rather than moving to accept or reject, because the defect is a missing derivation rather than an established contradiction.","tokens_in":9855,"tokens_out":15033,"duration_ms":162628,"concrete_test":"Perform an independent path-integral derivation of the TWM transmitted field Atr(ωpp) in two ways: (A) the paper's factorized route, with the bare D0 kernel contracted with the local (Fresnel-propagated, polariton-dressed) THz field as in Eq. (4); (B) the fully contracted route, with the external THz source connected through the exact mixed photon-phonon propagator ⟨A Q⟩ and a bare propagation factor for the generated signal. Compare the analytic expressions to leading order in the nonlinear coupling. If (A) and (B) are identical after incorporating the linear polariton self-energy, the factorization is justified and the central claim stands. If they differ, recompute Fig. 3(d) with K(2) replaced by the dressed D_Q kernel; if peak positions shift by more than a linewidth, the conclusion that TWM sees the dispersion only via propagation is invalid.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equations (2) and (3) carry the central contrast of the paper: FWM kernel K(3)∝R^2D_Q with the dressed phonon propagator, TWM kernel K(2)∝Ω(ZR)D_0 with the bare propagator. The justification, deferred to the Supplementary Material [47], is that in TWM the THz field is external and is not integrated out, so the phonon line in the cubic vertex is D0, while the polariton enters later through the refractive index n(ω) in Eq. (4). This factorization is the load-bearing assumption for the paper's main conclusion that TWM accesses the polariton dispersion only via propagation and phase-matching. It is not self-evident: the same phonon-photon coupling that dresses D0 into D_Q in FWM should, in a fully contracted calculation, dress the TWM vertex through the mixed photon-phonon propagator ⟨A Q⟩. If K(2) must actually contain D_Q (or the mixed propagator) when the external THz field is used undressed, then the relative heights and peak positions predicted by Eq. (4) change, and the statement that dispersion enters 'mainly via propagation' is not established. The main text gives no derivation of this factorization; it is an ansatz whose correctness is delegated to the missing supplement.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a many-body effective-action and perturbative Maxwell-equation framework for phonon-polariton nonlinear spectroscopy in non-centrosymmetric cubic crystals. In four-wave mixing (FWM), the third-order nonlinear kernel is argued to be proportional to the photon-dressed phonon propagator, K^(3) ∝ R^2 D_Q (Eq. 2), with D_Q given in Eq. (1). In three-wave mixing (TWM), the second-order kernel is instead argued to be proportional to the bare phonon propagator, K^(2) ∝ Ω(ZR)D_0 (Eq. 3), with all polariton effects entering through propagation and phase matching of the THz pump, encoded in Eq. (4). The paper then analyzes narrowband and broadband THz pump-optical probe signals in transmission and reflection, and claims that broadband THz pumps can access the full polariton phase space in a single measurement. Numerical illustrations are provided for representative parameters.","tokens_in":10197,"tokens_out":16337,"duration_ms":145163,"significance":"If the central contrast between Eqs. (2) and (3) holds, this framework is valuable: it unifies spontaneous and stimulated Raman-like scattering with THz-driven three-wave mixing, accounts for finite momentum mismatch, Fabry-Perot effects, and interface screening, and it makes falsifiable predictions, e.g., the phase-matched peak in transmission and the ω_TO-dominated reflection spectrum with a secondary ω_LO feature, which the authors connect to published experiments. The main strengths are the closed-form dressed phonon propagator, the analytic perturbative solution of the nonlinear Maxwell equation, and the explicit symmetry-based discussion of the TWM kernel. However, the paper does not include machine-checkable code, and the derivations of the most load-bearing equations are delegated to a Supplementary Material that is not part of the arXiv submission, so the present assessment of those steps relies on plausibility and on the consistency of the displayed expressions.","major_comments":[{"comment":"There appears to be an error in the printed form of the dressed phonon propagator. If the denominator of the last term is Ω_m^2 + c^2 ε∞ |k|^2 as written, then after analytic continuation the pole condition becomes ε(ω)ω^2 = ε∞ c^2 k^2, giving ω ≈ c k in the high-frequency limit, whereas the equivalent expression quoted in the same paragraph, D_Q = D_0(ε∞Ω_m^2 + c^2|k|^2)/(ε(iΩ_m)Ω_m^2 + c^2|k|^2), gives ε(ω)ω^2 = c^2 k^2 and is consistent with the axes of Fig. 1, where momenta are measured in units of c k/(ω_TO √ε∞). The denominator should presumably read Ω_m^2 + c^2|k|^2/ε∞. Since Eq. (1) controls all subsequent numerical results, this must be corrected and the figures rechecked.","section":"Eq. (1)"},{"comment":"The central claim that the TWM kernel is proportional to the bare propagator D_0, Eq. (3), while all polariton effects enter only through the refractive index n(ω) in Eq. (4), is asserted but not derived in the main text. In the FWM case the same linear Z-coupling is integrated out to produce the dressed propagator D_Q (Eq. (1) and Fig. 1(b)); in the TWM case one must show that the fully contracted nonlinear current genuinely contains no dressed phonon or mixed photon-phonon propagator. The text only states that one integrates out the transverse bare phonon and refers to the Supplementary Material [47], which is not included in this submission. If the exact TWM kernel contained D_Q or an ⟨A Q⟩ propagator, the relative peak heights and the paper's conclusion that TWM accesses the polariton dispersion 'mainly via propagation' would be invalid. I request that the explicit path-integral contraction for the TWM current be presented in the main text, at least at leading order, or that the relevant part of the supplement be reproduced.","section":"Three-wave mixing interaction and Propagation effects"},{"comment":"All numerical results for the TWM pump-probe signal are based on Eq. (4), but the perturbative solution of the nonlinear Maxwell equation is not derived. The reader is told that the procedure is detailed in [47], and the definitions of the internal fields A_t, A_r, the reflection coefficients r_±, and the Fabry-Perot factor f(ω) are only sketched. This is a load-bearing gap because the phase-matching factor and the interface terms are what convert the local kernel K^(2) into the predicted transmission and reflection spectra. The main text should at least state the Green's-function solution, the boundary conditions at the sample interfaces, and the order of the perturbative expansion; otherwise the figures cannot be independently checked.","section":"Propagation effects, Eq. (4)"},{"comment":"The abstract and the Summary state that broadband THz pumps allow one to access the full polariton phase space in a single transmission measurement, but the supporting evidence is not shown. Fig. 3(d) presents three curves with different values of n_eV/√ε∞, and the text says the phase-matched frequency is tuned by changing n_eV. No single-measurement protocol that maps the entire dispersion is described or demonstrated. Please clarify whether 'single measurement' means one pump spectrum with a fixed probe refractive index (in which case only one phase-matched point is selected, albeit with a finite phase-matching width) or whether the full dispersion is reconstructed from the entire transmitted spectrum. If the latter, an illustrative inversion or reconstruction should be shown.","section":"Pump-probe experiments and Abstract"}],"minor_comments":[{"comment":"The diagram labels in the submitted Figure 1 text appear garbled (e.g., 'DQ D0 D0 DQ Z Z = + R R!1,k1...'). Please provide a clean version with all propagator and field labels correctly placed.","section":"Figure 1"},{"comment":"In Eq. (4), the fields A^{σ1}(ω_pp) and A^{σ2}(ω-ω_pp) and the index σ are not fully defined before the equation is used; the subsequent definitions help, but the summation convention should be made explicit at the point of first use.","section":"Eq. (4) notation"},{"comment":"The paper relies heavily on Supplementary Material [47] for the derivations of Eqs. (1), (3), and (4), as well as for the tensorial structure of the TWM kernel and the sample-thickness discussion. Since the arXiv submission does not include this material, the manuscript is not self-contained; at minimum, the main text should state which results are essential and which are peripheral.","section":"Supplementary material"},{"comment":"The phrase 'non-centrosymmetric cubic systems' is used throughout, but the TWM selection rules are only fully meaningful once the crystal point group is specified; the sentence mentioning that the signal can disappear in specific polarization geometries according to the structure of the crystal (ref. [29]) would benefit from one concrete example, even if only in words.","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":"The paper should be suitable for publication in a condensed-matter physics journal after the listed revisions. The main technical risk is the TWM factorization of Eq. (3) and Eq. (4), which is the conceptual novelty of the paper, and the apparent typo in Eq. (1). The authors' use of their earlier work (Ref. [46]) for the Raman coupling is reasonable and not circular. In its current form, the manuscript is not fully self-contained because the Supplementary Material is absent from the arXiv posting; I would encourage the editor to request that the derivation material be made available or be moved into an appendix."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know: this is a solid, useful theory paper, not a breakthrough, and the main worry in the stress test does not land. What is new is that the authors put FWM and TWM phonon-polariton signals in one many-body language, showing that the dressed polariton propagator controls the FWM kernel while in TWM the polariton enters mainly through propagation and phase matching of the THz pump. That contrast makes the broadband-pump prediction concrete: a single transmission measurement can sweep a wide region of the polariton dispersion. The derivations look standard (path integrals plus perturbative Maxwell equation) and I see no obvious errors. The TWM factorization concern: I think it is a red herring. If the THz field is treated as an external source, the phonon line in the cubic vertex is indeed the bare D0; the polariton dressing shows up in how that same field propagates through the material via n(omega). This is the usual separation between a local nonlinear susceptibility and linear propagation. The stress-test worry about a mixed propagator <A Q> would only apply if you integrated out the THz fluctuations, which is not what the paper does. So the central claim is on solid ground. The weakness is that the main text states this factorization without deriving it, leaving the bookkeeping to the supplementary material, which is not included in the arXiv posting. A referee cannot verify the details without asking for the supplement. Soft spots: the supplementary material is missing from this submission, and the comparisons to experiments are qualitative. The phrase \"excellent agreement with experimental data [29]\" overshoots; no data are shown side by side. For a theory paper that is acceptable, but the language should be toned down. The figures are clear, the phase-matching intuition is pedagogically nice, and the parameter dependence (sample thickness, damping, refractive index) is sensible. Bottom line: this deserves a serious referee. I would send it to review with a request that the supplement be included or the TWM factorization justified in the main text. The experimental claims should be softened or shown quantitatively. For researchers working on THz pump-probe spectroscopy or phonon-polariton theory, this is a worthwhile read and a likely citation, though I would not cite it for my own work in the next year.","headline":"A credible and useful unified theory of phonon-polariton FWM and TWM responses, with the real derivation buried in a missing supplement; the TWM factorization concern is a red herring, and the paper deserves review.","tokens_in":663,"tokens_out":856,"would_cite":false,"duration_ms":42922,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["71.36.+c","78.20.Bh","78.47.jh"],"model":"deepseek-v4-flash","headline":"The paper shows that in THz pump–optical probe spectroscopy the phonon-polariton dispersion enters through propagation and phase matching rather than through the nonlinear kernel, so a broadband pump can access the full dispersion in one…","keywords":["phonon-polaritons","THz pump-probe spectroscopy","three-wave mixing","four-wave mixing","nonlinear optical kernel","phase matching","non-centrosymmetric crystals","polariton dispersion"],"falsifier":"Measure the transmitted TWM signal in a sample with known polariton dispersion while varying the optical probe frequency (i.e., $n_{eV}$). If the phase-matched peak does not shift according to Eq. (4) but instead stays pinned to the bare phonon frequency $\\omega_{TO}$ for all probe frequencies, the claim that polariton effects enter only through propagation would be falsified. Conversely, a clear shift of the peak with $n_{eV}$ in a thick sample supports the propagation-dominated picture.","tokens_in":9661,"feed_emoji":"🔬","tokens_out":5896,"duration_ms":52548,"temperature":0.7,"pith_summary":"The paper develops a many-body theory of how phonon-polaritons—hybrid lattice-light modes formed when infrared phonons couple to light—appear in two types of time-resolved nonlinear spectroscopy in non-centrosymmetric crystals. It shows that in four-wave mixing (e.g., impulsive stimulated Raman scattering) the measured nonlinear kernel is proportional to the dressed phonon-polariton propagator, allowing direct mapping of the polariton dispersion. In three-wave mixing (THz pump–optical probe), by contrast, the nonlinear kernel is proportional to the bare phonon propagator, and the polariton enters only through the propagation and phase-matching of the THz pump field. The paper predicts that broadband THz pumps can access the entire polariton phase space in a single transmission measurement, a practical advantage for studying light-matter hybridization.","feed_headline":"Terahertz pumps map the full phonon-polariton dispersion in one shot","feed_subtitle":"Because the polariton enters only through pump propagation, a single broadband pulse covers the whole phase space.","key_machinery":"The central object is the dressed phonon (phonon-polariton) propagator $D_Q(i\\Omega_m,k)=2/(\\omega_{TO}^2+\\Omega_m^2+\\Omega_P^2 \\Omega_m^2/(\\Omega_m^2+c^2|k|^2/\\epsilon_\\infty))$, obtained by integrating out the THz electromagnetic field to dress the bare phonon propagator $D_0=2/(\\omega_{TO}^2+\\Omega_m^2)$. Its poles give the two polariton branches $\\epsilon(\\omega)\\omega^2=c^2|k|^2$ weighted by a factor measuring the phonon character. The nonlinear current in each protocol is derived by functional differentiation of an effective action, and propagation is treated by a perturbative solution of the nonlinear Maxwell equation with momentum mismatch $\\Delta k$ and Fabry–Perot factors. This combination—many-body kernel plus perturbative propagation—is what allows the paper to separate the roles of hybridization from phase matching.","core_discovery":"The central result is a pair of proportionality relations for the nonlinear optical kernels. In four-wave mixing, the kernel is $K^{(3)}(\\Omega,k)\\propto R^2 D_Q(\\Omega,k)$, where $D_Q$ is the phonon propagator dressed by its coupling to the electromagnetic field (Eq. 1); the response therefore directly measures the phonon component of the polariton. In three-wave mixing, the kernel is $K^{(2)}(\\Omega)\\propto \\Omega (ZR) D_0(\\Omega)$, proportional to the bare phonon propagator $D_0$ rather than the dressed one. The polariton still shapes the measured signal, but only through the linear propagation of the THz pump inside the sample, as encoded in the refractive index $n(\\omega)$ and the phase-matching factor in Eq. (4). The consequence is that in THz pump–optical probe experiments the polariton dispersion is probed primarily via propagation and phase matching, and a broadband pump naturally covers many momenta at once.","pith_inferences":["If the factorization holds, then a natural extension is to use the THz-pump–optical-probe signal as a direct measure of the polariton refractive index rather than of the bare phonon: the phase-matched peak position is a clean proxy for $n(\\omega)$ at the pump frequency.","The same framework could be adapted to extract the density of states of the polariton from broadband pump spectra, since the phase-matching factor acts as a filter along the dispersion.","A testable extension: in a thin sample ($d \\lesssim c/(\\omega_{TO}\\sqrt{\\epsilon_\\infty})$) the phase-matching condition loses meaning, so the TWM signal should become nearly independent of $n_{eV}$; this crossover is a direct fingerprint of the propagation-dominated mechanism.","The many-body effective-action derivation suggests that coupling the phonon to an optical cavity should modify the TWM kernel through the dressed propagator even when the bare phonon is probed, providing a route to cavity-enhanced nonlinear spectroscopy."],"forward_implications":["In four-wave mixing (ISRS), the measured signal directly tracks the phonon component of the polariton; tuning the optical pulse frequency moves the phase-matched point along the dispersion.","In three-wave mixing with a narrowband THz pump, a sizable signal is predicted even when the momentum mismatch $\\Delta k\\neq 0$, explaining why stimulated IR experiments see response on the upper polariton branch far from perfect phase matching.","In three-wave mixing with a broadband pump, the transmitted signal in a single measurement covers a wide range of polariton momenta, so the full dispersion can be reconstructed by varying the optical refractive index $n_{eV}$.","In reflection geometry the phase-matched point sits at high momentum where the polariton is nearly bare phonon, so the reflected signal peaks near $\\omega_{TO}$, with a secondary feature at $\\omega_{LO}$ from interface screening.","The formalism extends to anisotropic crystals, anharmonic couplings, and magnon-phonon-polaritons, and can describe cavity polariton pump-probe experiments by treating light and phonons on equal footing."],"supporting_citations":[{"why":"Establishes the ISRS protocol whose four-wave-mixing signal the theory explains as proportional to the dressed phonon propagator.","marker":"[3]"},{"why":"Supplies the phase-matching condition (intersection of the polariton dispersion with ck/n_eV) that the paper uses to locate signal peaks.","marker":"[18]"},{"why":"Provides the experimental THz pump-optical probe reflection data that the theory reproduces qualitatively, including the peak near ω_TO.","marker":"[29]"},{"why":"Gives the stimulated-IR Raman experiment against which the narrowband TWM results are compared, validating the finite-momentum-mismatch response.","marker":"[30]"},{"why":"Introduces the Raman-like quadratic coupling term (QT RAA) that is the detection mechanism in both FWM and TWM protocols.","marker":"[46]"},{"why":"Contains the detailed derivations of the dressed-phonon propagator, the nonlinear currents, and the perturbative propagation solution used throughout the main text.","marker":"[47]"}],"fun_headline_variants":["One THz pulse maps entire phonon-polariton phase space","Polariton dispersion from a single broadband THz pump","THz pump-optical probe: full polariton dispersion without tuning","Single broadband THz pulse accesses full polariton dispersion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The factorization assumes that in three-wave mixing the nonlinear current can be computed from the bare phonon propagator and the undressed THz field, with all polariton effects confined to linear propagation; a perturbative ansatz that is not derived from the full theory.","fun_headline_variants_meta":{"raw":{"variants":["One THz pulse maps entire phonon-polariton phase space","Polariton dispersion from a single broadband THz pump","THz pump-optical probe: full polariton dispersion without tuning","Single broadband THz pulse accesses full polariton dispersion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000469,"raw_usage":{"total_tokens":2320,"prompt_tokens":911,"completion_tokens":1409,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":527,"completion_tokens_details":{"reasoning_tokens":1340}},"tokens_in":527,"tokens_out":1409,"duration_ms":10369,"temperature":1.0,"reasoning_tokens":1340,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T19:54:22.331785+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the transmitted TWM signal in a sample with known polariton dispersion while varying the optical probe frequency (i.e., $n_{eV}$). If the phase-matched peak does not shift according to Eq. (4) but instead stays pinned to the bare phonon frequency $\\omega_{TO}$ for all probe frequencies, the claim that polariton effects enter only through propagation would be falsified. Conversely, a clear shift of the peak with $n_{eV}$ in a thick sample supports the propagation-dominated picture.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the ISRS protocol whose four-wave-mixing signal the theory explains as proportional to the dressed phonon propagator."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the phase-matching condition (intersection of the polariton dispersion with ck/n_eV) that the paper uses to locate signal peaks."},{"cited_title":"Fu and M","cited_arxiv_id":null,"evidence_quote":"Provides the experimental THz pump-optical probe reflection data that the theory reproduces qualitatively, including the peak near ω_TO."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the stimulated-IR Raman experiment against which the narrowband TWM results are compared, validating the finite-momentum-mismatch response."},{"cited_title":"Udina, T","cited_arxiv_id":null,"evidence_quote":"Introduces the Raman-like quadratic coupling term (QT RAA) that is the detection mechanism in both FWM and TWM protocols."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Contains the detailed derivations of the dressed-phonon propagator, the nonlinear currents, and the perturbative propagation solution used throughout the main text."}],"review_version":1}