{"id":"9e433364-eafc-49cc-97f4-23dc3a4e6990","arxiv_id":"1908.07472","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Quadratic observables of Gaussian beam wave solutions have wavelength-independent stochastic regularity, including time-averaged two-mode observables.","lead":"This paper proves that certain quadratic observables of high-frequency wave solutions, such as energy and Arias intensity, have random-parameter derivatives bounded independently of the wavelength when computed with Gaussian beam approximations. This enables fast spectral uncertainty quantification for seismic, acoustic, and optical wave propagation.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: Proposition 5.5 and the non-stationary phase step are sound; only minor typos in Eq. (38) and a support-radius statement need correction.","rationale":"The reader's conditional verdict is driven by delegation to prior work and by typos. My independent read found the two-mode cross-term proof internally sound. The sign typo in Eq. (38) is confusing but does not alter the lower bound, because the absolute value of the leading term is unchanged. The support-radius typo after Eq. (41) does not affect the estimates, since Proposition 5.5 holds on the larger set Σ_μ. The reliance on [23] for Theorem 4.6 is a presentation and verifiability issue rather than a detected correctness gap; the cited lemmas appear to cover the needed estimates with t-dependent ψ and compact Γ_c. I therefore would not change the reader's verdict: the paper still merits conditional acceptance with minor corrections.","tokens_in":25357,"tokens_out":28176,"duration_ms":300972,"concrete_test":"Re-derive Eq. (38) from (11) and Proposition 3.1 with the correct sign qdot^- = −c(q^-)p^-/|p^-|; if the leading terms are +c^-|p^-| + c^+|p^+|, recompute the R_k bound and confirm |∂_t ϑ_k| ≥ ν with ν > 0 for the chosen μ. This settles the only step that could undermine the non-stationary phase argument.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Theorem 5.2) appears to be supported by the proof. I examined the two places where the argument could fail. First, Proposition 5.5 requires |∂_t ϑ_k| ≥ ν on Σ_μ. Up to the sign typo in Eq. (38) (with (11) one gets +c(q^-)|p^-| + c(q^+)|p^+|, not the displayed negative signs), the lower bound 2γ − C_k μ is valid: both Hamiltonians are conserved and at least γ > 0 by (A1),(A3), and the remainder R_k is O(μ) on Σ_μ. Second, the partition-of-unity application of Lemma 5.3 is valid: on the g1 part the phase derivative is bounded away from zero, and on the g2 part Im ϑ_k ≥ δ μ^2 gives exponential decay. The one-mode input (Theorem 4.2/4.6) is delegated to [23] and [22], which is a verifiability weakness for a standalone paper, but I did not find an internal inconsistency; the appendix supplies the needed derivative expansion. The false support statement 'supp g... ⊂ Σ_{μ/2}' after Eq. (41) should read Σ_μ, but the subsequent estimate only needs the larger set, so it is harmless. No load-bearing flaw identified.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the stochastic regularity of quadratic observables of high-frequency solutions to the scalar wave equation with random coefficients and initial data, in the regime where the wavelength ε is small. The quantities of interest have the form ε^{2(p+|α|)}∫∫ g |∂_t^p ∂_x^α u|^2 ψ dx dt, and the solution u is approximated by a sum of two Gaussian-beam modes. The main result, Theorem 5.2, states that under hypotheses (A1)–(A5) and with an admissible beam cutoff width η, every mixed derivative of the Gaussian-beam QoI Q^{GB}_{p,α}(y) with respect to the parameters y is bounded uniformly in ε on compact parameter sets. The proof splits the QoI into same-mode contributions, handled by the one-mode theorem (Theorem 4.2), and a cross term. The cross term is controlled by a non-stationary phase argument in time: Proposition 5.5 shows that the phase derivative |∂_t ϑ_k| is bounded below by a positive constant on the relevant spatial region, and Lemma 5.3 converts this into arbitrarily high powers of ε. Numerical examples in Sections 5.1 and 5.4 illustrate the difference between the oscillatory space-only QoI and the non-oscillatory time-averaged QoI.","tokens_in":25639,"tokens_out":12013,"duration_ms":122993,"significance":"The result is a substantive extension of the authors' earlier one-mode result [23] to two-mode Gaussian-beam solutions and to observables involving higher derivatives of the wave field. If correct, it provides a theoretical justification for stochastic collocation and stochastic Galerkin methods applied to Gaussian-beam approximations in the high-frequency regime, which is the regime where direct numerical simulation is prohibitively expensive. The proof is careful and largely explicit: the ε-power bookkeeping is tracked through the non-stationary phase lemma, the constants are kept visible, and the assumptions (A1)–(A5) are used in an essential way. The numerical examples are helpful and confirm the qualitative difference between space-only and space-time averaged observables. The paper relies on published prior results for the one-mode input, which is acceptable but makes the paper not fully self-contained.","major_comments":[],"minor_comments":[{"comment":"There is a sign error in the displayed formula for ∂_t ϑ_k. Using (11), ∂_t q^± = ± c(q^±)p^±/|p^±|, so the right-hand side should read +c(q^-,y)|p^-| + c(q^+,y)|p^+| + R_k, not −c(q^-,y)|p^-| − c(q^+,y)|p^+| + R_k. The subsequent lower bound in Proposition 5.5 uses the positive sum, so the display should be corrected to match the argument.","section":"Section 5.3, Eq. (38)"},{"comment":"The support statement for g_{ℓmjσ} := f_{ℓmjσ} ψ g_1 says supp g_{ℓmjσ}(t,·) ⊂ Σ_{μ/2}, but g_1 = ϱ_μ(x−q^+)ϱ_μ(x−q^-) is supported only in Σ_μ (since ϱ_μ has support in B_{2μ}), not in Σ_{μ/2}. The correct statement is supp g_{ℓmjσ}(t,·) ⊂ Σ_μ, and the argument still works because Proposition 5.5 provides the lower bound on Σ_μ.","section":"Section 5.3, after Eq. (41)"},{"comment":"Theorem 4.6, which is the load-bearing input for the one-mode result, is not proved in this paper: the appendix invokes [23, Lemma 5] and [23, Lemma 6] and says the rest of the proof of [23, Theorem 1] can be used as is. Proposition 4.1 is also delegated to [22] and [23]. Since Theorem 5.2 depends on these results, the paper would be easier to verify if the authors stated explicitly which lemmas from [23] are reused and confirmed that they apply under the present slightly generalized hypotheses (t-dependent ψ, compact subsets of an open set Γ).","section":"Section 4.2 and Appendix"},{"comment":"There are small notation inconsistencies in the statements: the spatial multi-index is written as α ∈ N_0^N in places where it should be α ∈ N_0^n, and the parameter derivative multi-index is σ ∈ N_0^N. Also, Eq. (5) in the introduction writes ∂^σ Q/∂y^σ twice.","section":"Theorems 4.2 and 5.2"},{"comment":"The illustrative examples use Gaussian initial data with unbounded support and, in Section 5.4, a non-smooth phase containing |x_1|, so they do not satisfy assumptions (A2) and (A3). This is acceptable for heuristic illustration, but the text should note explicitly that these examples lie outside the theorem's hypotheses.","section":"Sections 5.1 and 5.4"}],"recommendation":"minor_revision","confidential_remarks":"The paper is technically sound in my reading: the main estimates are plausible, the non-stationary phase step is valid, and the ε-balance at the end of Theorem 5.2 works. The only issues are local typos and the heavy reliance on the authors' prior published results for the one-mode input. If the editor values self-containedness, I would encourage asking the authors to expand the appendix; otherwise the paper is publishable after minor revisions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nRead it. The genuinely new content is in Section 5: for a two-mode Gaussian beam superposition, the quadratic observable becomes non-oscillatory when averaged in time, because the cross-term phase has no stationary point in t. The paper proves epsilon-uniform bounds on all parameter derivatives for arbitrary derivative orders p, alpha, and it handles both the one-mode case (extending the earlier p=alpha=0 result) and the two-mode case. The main mechanism is clean: scaled beam derivative classes plus non-stationary phase with tracked constants. I checked the places where a proof like this usually fails, and the argument holds up.\n\nWhat the paper does well: it identifies a real obstacle (oscillatory cross terms for opposite-propagating modes), gives a simple 1D example showing the failure, and then shows the time-averaging cure is not just heuristic but provable under assumptions (A1)-(A5). The numerical experiments are consistent with the theorems. The self-citations to [22,23,24] are substantial but legitimate: those results are published with proofs, and the new theorems do not assume their own conclusion.\n\nSoft spots, in proportion. The biggest one for a standalone reader is delegation: Theorem 4.6, the one-mode workhorse, is proved by saying \"same as Theorem 1 in [23]\" plus an appendix that still calls on lemmas from [23], and Proposition 4.1 is also delegated to [22]. The estimates are probably right, but an independent referee will have to go fetch a thesis and a prior paper to verify the foundation. That should be fixed or at least acknowledged by expanding the appendix. There are also two harmless typos in the main proof: Eq. (38) has the wrong sign for the leading term (with the sign convention in (11) it should be +c(q^-)|p^-| + c(q^+)|p^+|), and the support statement after Eq. (41) says supp g2 subset Sigma_{mu/2} when it should be Sigma_mu; the subsequent estimate only needs the larger set, so no damage. The one place I could imagine a real objection is Proposition 5.5: the lower bound |partial_t theta_k| >= nu relies on both Hamiltonians staying above gamma and on the beam cutoff width mu being small. If a mode had a stationary phase point in t or the two phase speeds canceled, the non-stationary phase argument would fail. But that is not a hidden loophole; the paper's 1D example shows such failure can genuinely happen, and the assumptions are chosen to exclude it.\n\nWho is this for: people working on stochastic collocation or Galerkin for high-frequency waves, especially with Gaussian beams. A specialist in UQ for wave propagation will get real value; a general numerical analyst will probably not care much. Verdict: serious referee, yes. Send one. Ask for the typos to be fixed and for the one-mode estimates to be made less dependent on external sources; neither should require a change to the main results.","headline":"Solid, honest generalization of the authors' earlier one-mode regularity result to the two-mode time-averaged case; the core estimates check out and it deserves a serious referee with minor revisions.","tokens_in":26173,"tokens_out":1395,"would_cite":true,"duration_ms":17963,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M99","35L05","60H35"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves that time-averaged quadratic observables of high-frequency waves, computed by Gaussian beams, have parameter derivatives uniformly bounded as the wavelength tends to zero.","keywords":["stochastic regularity","Gaussian beams","high-frequency waves","uncertainty quantification","quadratic observables","sparse grid collocation","wave equation","non-stationary phase"],"falsifier":"Take the two-mode example in Section 5.4 with $\\phi_0=x_1+(x_2-y_1)^2$ and compute $\\partial_y^2 Q^{\\mathrm{GB}}_{0,0}$ for $\\varepsilon = 1/80$, $1/160$, and $1/320$ at the same parameters. Theorem 5.2 predicts the plotted values remain $O(1)$ as $\\varepsilon$ shrinks; if the peak values grow like $\\varepsilon^{-2}$, the $\\varepsilon$-uniformity claim is false.","tokens_in":25155,"feed_emoji":"🌊","tokens_out":9820,"duration_ms":102300,"temperature":0.7,"pith_summary":"High-frequency wave fields oscillate on a scale $\\varepsilon$, and quantities computed from them typically inherit that oscillation, making their dependence on uncertain parameters increasingly rough as $\\varepsilon$ shrinks, which is the opposite of what fast stochastic methods require. The paper identifies a class of quadratic observables that escape this problem: weighted space-time integrals of $|\\partial_t^p\\partial_x^\\alpha u|^2$, evaluated on a Gaussian-beam approximation of the wave field. For these observables, every derivative with respect to the stochastic parameters is bounded by a constant independent of $\\varepsilon$ (Theorem 5.2), and the same holds for spatial-only one-mode observables uniformly in time (Theorem 4.2). The reason is that the interference term between the two counter-propagating beam modes is killed by time integration through non-stationary phase, so the parameter-to-observable map stays uniformly smooth in the high-frequency limit.","feed_headline":"Time averaging restores smoothness in high-frequency wave outputs","feed_subtitle":"For Gaussian-beam approximations, all parameter derivatives of these observables stay bounded independently of the wavelength.","key_machinery":"The argument is carried by the $k$-th order Gaussian beam superposition $u_k=u_k^++u_k^-$, whose phase and amplitude coefficients are solutions of a set of $\\varepsilon$-independent ODEs. When the observable is expanded, the difficulty is concentrated in the cross term whose integrand contains the phase $\\vartheta_k(t,x,y,z,z')=\\Phi_k^-(t,x-q^-,y,z')-(\\Phi_k^+)^*(t,x-q^+,y,z)$. Proposition 5.5 shows $|\\partial_t\\vartheta_k|\\ge\\nu>0$ on the set where both beams overlap the test function, because the conserved ray Hamiltonians $c(q^\\pm)|p^\\pm|$ are bounded below by $\\gamma>0$ and the beam cutoff width $\\mu$ is chosen small enough. The non-stationary phase lemma then bounds the time integral by $\\varepsilon^K$ times constants involving derivatives of the amplitude, and the positivity of $\\operatorname{Im}\\vartheta_k$ from the admissible-cutoff condition gives exponential decay $e^{-\\delta|x-q^\\pm|^2/\\varepsilon}$ outside the overlap. This combination converts every oscillatory integral into an $\\varepsilon$-independent bound after taking finitely many derivatives in $y$.","core_discovery":"The central result is the $\\varepsilon$-uniform stochastic regularity of a general class of quadratic observables. Under assumptions (A1)--(A5), for a fixed final time $T$, beam order $k$, compact parameter set $\\Gamma_c$, and an admissible cutoff width $\\eta$, the Gaussian-beam observable $$$Q^{{\\mathrm{GB}}$}_{p,\\$\\alpha$}(y)=\\$varepsilon^{{2(p+|\\alpha|)}}$\\int\\!\\int g(t,x,y)|\\partial_t^p\\partial_x^\\$\\alpha$ u_k(t,x,y)|^2\\psi(t,x)\\,dx\\,dt$$ satisfies $\\sup_{y\\in\\Gamma_c}|\\partial_y^\\sigma Q^{\\mathrm{GB}}_{p,\\alpha}(y)|\\le C_\\sigma$, where $C_\\sigma$ is independent of $\\varepsilon$, for every $p$, multi-index $\\alpha$, and derivative multi-index $\\sigma$. The one-mode version, with $u_k$ replaced by a single beam family, satisfies the stronger statement that the same bound holds for each fixed time $t\\in[0,T]$. The new work beyond the earlier one-mode result is the treatment of the two-mode cross term: the relative phase of the two beams has a time derivative uniformly bounded away from zero wherever both beams reach the measurement window, so the non-stationary phase lemma makes the cross contribution of order $\\varepsilon^K$ for any $K$, while Gaussian decay and the cutoffs control everything else.","pith_inferences":["Editor's extension: the same cross-term mechanism should cover bilinear observables such as $\\varepsilon^{2m}\\int\\!\\int g\\,\\overline{D^m u}\\,D^m v$ for two different solutions, since only the smooth amplitude changes and the phase difference has the same structure.","Editor's extension: the uniform positivity of $\\partial_t\\vartheta_k$ is a no-resonance condition; in media with periodic structure or with beams that turn around, stationary phase points can appear, and one would expect $\\varepsilon$-dependent observables unless additional averaging in space or parameters is introduced.","Editor's extension: because the proof only uses the sign and size of the conserved ray Hamiltonians, the result should extend to systems of hyperbolic equations and to Schrödinger-type equations whose characteristics have conserved group speeds; the analogous phase derivative would be a difference of group velocities."],"forward_implications":["Sparse-grid stochastic collocation for these observables keeps its fast convergence as $\\varepsilon\\to0$, because the polynomial and spline interpolation errors in $y$ are controlled by the same $\\varepsilon$-independent derivative bounds.","Physical outputs of the form (3), including acoustic potential energy, total energy, and Arias intensity, are covered, so the result applies to seismic and acoustical quantities of interest used in practice.","When one mode dominates, the spatial-only observable is regular at every fixed time, allowing time-dependent statistics without temporal averaging.","The estimates hold for arbitrarily high derivatives in the stochastic parameters, so refinement algorithms that estimate error from higher derivatives will not see $\\varepsilon$-dependent growth for these outputs."],"supporting_citations":[{"why":"Supplies the one-mode regularity theorem and the proof lemmas (symmetrization, derivative expansion) that the two-mode proof directly extends.","marker":"[23]"},{"why":"Proves the imaginary-part lower bound for the Gaussian-beam phase, which defines the admissible cutoff width used in both theorems.","marker":"[22]"},{"why":"Provides the Gaussian-beam initial data and coefficient ODEs, giving the smoothness and compact-support properties of the beam coefficients used throughout.","marker":"[19]"},{"why":"Establishes the Riccati ODE solvability and the higher-order beam construction that make the phase and amplitude smooth functions of time and parameters.","marker":"[28]"},{"why":"Supplies the non-stationary phase lemma with explicit constants, which is the mechanism that damps the two-mode cross term.","marker":"[13]"},{"why":"Provides the Sobolev error estimate between the beam superposition and the exact solution, connecting the Gaussian-beam observable regularity to the original quantity of interest.","marker":"[20]"}],"fun_headline_variants":["High-frequency waves get smooth statistics via time averaging","Gaussian beams yield wavelength-independent derivative bounds","Quadratic wave observables attain epsilon-uniform regularity","Time averaging restores smoothness in wave-field statistics","Stochastic regularity for high-frequency wave outputs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The two-mode estimate rests on the claim that, whenever both beams reach the measurement region, the relative phase between them changes at a rate that stays strictly positive; this requires the wave speed and the initial phase gradient to be bounded below by positive constants and the beam width to be chosen small enough.","fun_headline_variants_meta":{"raw":{"variants":["High-frequency waves get smooth statistics via time averaging","Gaussian beams yield wavelength-independent derivative bounds","Quadratic wave observables attain epsilon-uniform regularity","Time averaging restores smoothness in wave-field statistics","Stochastic regularity for high-frequency wave outputs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000978,"raw_usage":{"total_tokens":4170,"prompt_tokens":975,"completion_tokens":3195,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":591,"completion_tokens_details":{"reasoning_tokens":3125}},"tokens_in":591,"tokens_out":3195,"duration_ms":25057,"temperature":1.0,"reasoning_tokens":3125,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:19:40.544751+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the two-mode example in Section 5.4 with $\\phi_0=x_1+(x_2-y_1)^2$ and compute $\\partial_y^2 Q^{\\mathrm{GB}}_{0,0}$ for $\\varepsilon = 1/80$, $1/160$, and $1/320$ at the same parameters. Theorem 5.2 predicts the plotted values remain $O(1)$ as $\\varepsilon$ shrinks; if the peak values grow like $\\varepsilon^{-2}$, the $\\varepsilon$-uniformity claim is false.","supporting_citations":[{"cited_title":"Malenová, M","cited_arxiv_id":null,"evidence_quote":"Supplies the one-mode regularity theorem and the proof lemmas (symmetrization, derivative expansion) that the two-mode proof directly extends."},{"cited_title":"Malenová.Uncertainty quantification for high frequency waves","cited_arxiv_id":null,"evidence_quote":"Proves the imaginary-part lower bound for the Gaussian-beam phase, which defines the admissible cutoff width used in both theorems."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Gaussian-beam initial data and coefficient ODEs, giving the smoothness and compact-support properties of the beam coefficients used throughout."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the Riccati ODE solvability and the higher-order beam construction that make the phase and amplitude smooth functions of time and parameters."},{"cited_title":"Hörmander","cited_arxiv_id":null,"evidence_quote":"Supplies the non-stationary phase lemma with explicit constants, which is the mechanism that damps the two-mode cross term."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the Sobolev error estimate between the beam superposition and the exact solution, connecting the Gaussian-beam observable regularity to the original quantity of interest."}],"review_version":1}