{"id":"ac9d8fb5-1cf1-40b2-b352-8bb6069beb33","arxiv_id":"2411.14178","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A ray-theoretic framework using vertical acoustic modes predicts signal parameters and space-time caustics in shallow water.","lead":"The paper develops a mathematical method, called space-time horizontal rays, to describe how sound pulses travel in shallow water. It derives equations for the path, amplitude, and frequency changes of a pulse, including the locations of caustics where the signal becomes focused.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The inner product (2.22) is incompatible with the interface condition (2.19), so the identity (2.32) used to close the transport equation has a missing boundary term and fails for realistic density contrasts.","rationale":"The paper's central claim is that the method predicts signal parameters. The derivation of the amplitude prediction passes through Eq (2.36), which is obtained from Eq (2.35) by substituting Eq (2.32). If Eq (2.32) fails because of the density weighting, the amplitude law and the rest of Section 4 lose their foundation. This is more specific than the reader's 'uncontrolled approximation': it is an identifiable missing boundary term, not just an absent error estimate. I checked the possible rescue: replacing rho by 1/rho in (2.22) makes the operator self-adjoint and Eq (2.32) exact via the Hellmann-Feynman formula, so the issue is likely fixable. Therefore I leave the verdict at CONDITIONAL, but the revision must correct the inner product or the interface condition and re-derive the transport equation; the paper should not be accepted in its current form.","tokens_in":14008,"tokens_out":18329,"duration_ms":167018,"concrete_test":"In a two-layer Pekeris waveguide with water density 1 and bottom density 2, solve (2.17)-(2.19) for the first two modes, compute the integral (2.22) for the first mode, and compare q/k0 * dq/dk0 with ⟨n^2ψ,ψ⟩. If the relative difference is of order the density contrast and nonzero, Eq (2.32) is not a harmless uncontrolled approximation but a false identity; the transport equation must be re-derived with the boundary term. An even simpler analytic check is to keep the boundary term in ⟨ψ'',ψ⟩ and verify that it does not vanish.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The transport equation (2.36) and the amplitude formula (2.53) rest on the identity (2.32), ⟨n^2ψ,ψ⟩ = (q/k0) ∂q/∂k0. This identity is exact for a Sturm-Liouville operator that is self-adjoint in the inner product used. But the weighted inner product (2.22) uses rho_+ and rho_- as weights, while the interface condition (2.19) imposes continuity of (1/rho) ψ_z. With these two choices, L = d^2/dz^2 + n^2 k0^2 is not self-adjoint: integrating by parts gives ⟨ψ'',ψ⟩ = -⟨ψ',ψ'⟩ + (rho_+^2/rho_- - rho_-) ψ'(h-)ψ(h), and the boundary term vanishes only if rho_+ = rho_- or ψ'(h)=0. Hence the displayed expression for ⟨n^2ψ,ψ⟩ in (2.32) is not valid, and the subsequent divergence-form transport equation and the predicted signal amplitudes inherit the error. If (2.22) was intended to use 1/rho as weight, the identity is restored; as written, the derivation is inconsistent.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a space-time horizontal ray method for acoustic wave propagation in a shallow-water waveguide. It starts from the wave equation (2.10) and an ansatz (2.11) separating a vertical mode ψ from a two-dimensional phase φ. Substitution yields an eikonal equation (2.23) in which the vertical-mode eigenvalue q plays the role of a refractive index, and a transport equation (2.36) for the amplitude A. The associated Hamiltonian system (2.29) defines space-time rays; the paper then derives a linearized system (3.30)-(3.33) for ray-coordinate variations, uses it to compute the Jacobian of the transition to ray coordinates, and characterizes space-time caustics as zeros of that Jacobian. It also treats phase fronts, τ-fronts, s-fronts, and multiple-ray interference. The derivation is formal: no numerical examples, no error estimates, and no comparison with full-wave solutions are given.","tokens_in":14299,"tokens_out":17473,"duration_ms":148774,"significance":"If the derived equations are correct, the method provides a compact description of pulsed signal propagation in a dispersive waveguide, with amplitude, frequency modulation, and front geometry expressed through the single function q(⃗r,k0). The derivation is self-contained and does not introduce fitted parameters; I also find no circularity, and the potential inconsistency between the inner product (2.22) and the interface condition (2.19) is resolved by the interface condition itself. However, the transport equation rests on several unquantified approximations, and the absence of numerical validation means the claimed predictive power is not demonstrated. The paper is a useful formal framework that could become a practical tool after further justification and testing.","major_comments":[{"comment":"The approximation ⟨n^2ψ,ψ⟩ ≈ (q/k0)∂q/∂k0 is used without derivation or error estimate. Under the normalization (2.22), this identity is exact if one differentiates the eigenvalue equation Lψ = q^2ψ with respect to k0 and uses self-adjointness; the '≈' is therefore misleading. More importantly, Eq. (2.34) omits the term (1/2)⟨∂n^2/∂τ ψ,ψ⟩, which is generally nonzero and contributes to the transport equation. The authors should prove the simplified transport equation (2.36) under explicit scale-separation assumptions or provide a remainder estimate.","section":"§2.4, Eqs. (2.32)-(2.36)"},{"comment":"The equality ⟨(∇φ,∇ψ),ψ⟩ = ½(∇φ,∇⟨ψ,ψ⟩) is not generally valid, because the inner product (2.22) depends on the horizontal coordinates through ρ± and h(⃗r). Differentiation with respect to x introduces boundary terms from ∇h and terms from ∇ρ that are not negligible in the asymptotic ordering of §2.1. These missing terms affect the divergence form (2.36) and the amplitude formula (2.53). The transport equation should be derived from the exact modal relations, keeping all parametric derivatives of the inner product.","section":"§2.4, Eq. (2.33)"},{"comment":"The paper contains no numerical example or comparison with a direct solution. The abstract promises a simple method for predicting signal form and parameters, but no demonstration is provided. I request at least one benchmark (for instance, a range-dependent Pekeris waveguide) in which the amplitude predicted by (2.53) is compared with a full-wave solution, to assess the quantitative accuracy of the approximations and to establish the method's practical range of validity.","section":"Throughout (no numerical section)"},{"comment":"The coherence condition on the initial data is stated as a first-order PDE for φ0, but the existence of solutions is not analyzed. For a prescribed family of rays (ρ0,r0,k0,α0), the compatibility of (3.59) is not obvious. The paper should give conditions on the initial functions that guarantee solvability, or explain how φ0 is constructed in practice, since this condition is essential for the ray-coordinate framework.","section":"§3.4, Eqs. (3.58)-(3.59)"}],"minor_comments":[{"comment":"The definition of Δα(τ) subtracts the same expression from itself, making it identically zero; this is a typo and should read α(τ,μ,ν) - α(τ,μ',ν') or the reverse. The same issue appears in the initial-condition formula that follows.","section":"§3.1, Eq. (3.13)"},{"comment":"The notation 'Δ_|||' (with three vertical bars) appears in equations (3.7), (3.11), (3.18), (3.30), and elsewhere; it should be 'Δ_||' (two bars).","section":"§3.1-§3.2"},{"comment":"There are minor typos: 'vertical mode s' in the abstract should be 'vertical modes', and 'ommit' in footnote 1 should be 'omit'.","section":"Abstract and footnote 1"},{"comment":"The bare '≈' in (2.32) and (2.34) should be replaced by an explicit statement of the approximation regime (e.g., weak frequency dependence of the vertical modes) and an indication of the order of the neglected terms.","section":"§2.4, Eq. (2.32)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a formal derivation with plausible structure but no numerical validation. The main results depend on a transport equation whose derivation contains unquantified approximations; these must be addressed before the predictive claims can be accepted. I would encourage the editor to request a revision that includes a concrete example and a rigorous derivation of the transport equation under stated assumptions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The variation system (3.30)-(3.33) and the Jacobian-based caustic criterion are the genuine additions here, and they are worked out carefully. The broader framework — vertical modes plus space-time rays — is known from Babich and from Burridge-Weinberg, so the novelty is incremental but legitimate. The eikonal and ray equations are standard but cleanly derived, and the front analysis in Section 4 is clear. Credit where due: the paper is honest about omitting mode coupling, and the bibliography is appropriate.\n\nThe problem is the transport equation. The inner product (2.22) uses rho as a weight, but the interface condition (2.19) makes L = d^2/dz^2 + n^2 k0^2 non-self-adjoint in that inner product. A direct integration by parts gives a boundary term at z = h proportional to (rho_+ - rho_-^2/rho_+) psi'(h-) psi(h). That term is missing from eq. (2.31). Consequently eq. (2.32) is not the identity they claim, and the simplified transport equation (2.36) — and the amplitude formula (2.53) — do not follow. If the inner product were weighted with 1/rho, the identity would hold; as written it does not. The authors might be able to fix this, but the paper does not.\n\nThis is a load-bearing flaw, not a cosmetic one. The abstract promises a simple method for predicting signal amplitude and frequency modulation; that promise depends on (2.36). There is also no numerical test or worked example, and the approximation in (2.32) is uncontrolled even apart from the missing boundary term. So the paper currently supports the ray geometry and the caustic criterion, but not the amplitude predictions.\n\nWho gets value from this: a specialist in underwater acoustics or asymptotic methods who is interested in the variation system and the caustic Jacobian, and who is willing to do some algebra to repair the transport equation. The paper is not ready for use as a predictive tool as written.\n\nRecommendation: send it to peer review, but flag the self-adjointness issue prominently. It is the sort of specific, fixable mathematical error that referees can work through. With the transport equation corrected and at least one illustrative example or error estimate, the paper would be a solid contribution. Without that, it remains a promising sketch.","headline":"The variation system and Jacobian caustic criterion are real extensions of the vertical-mode horizontal-ray framework, but the transport equation rests on a self-adjointness identity that fails for realistic density contrasts, so the amplitude formula and signal parameter predictions are not justified.","tokens_in":14774,"tokens_out":5255,"would_cite":false,"duration_ms":44519,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35A17","35L05"],"pacs":[],"model":"deepseek-v4-flash","headline":"Predict a pulse's shape from one waveguide eigenvalue","keywords":["space-time rays","vertical modes","horizontal rays","underwater acoustics","eikonal equation","transport equation","space-time caustics","frequency modulation"],"falsifier":"A decisive check would be a numerical experiment in a range-dependent shallow-water waveguide: compute the exact vertical-mode inner product $\\langle n^2\\psi,\\psi\\rangle$ as a function of $k_0$, compare it with $(q/k_0)\\partial q/\\partial k_0$, and then propagate a broadband pulse with the ray formulas and with a full-wave solver. If the relative error in the inner product is not small, or if the ray-predicted amplitude and frequency modulation diverge from the full-wave result near a known caustic, the transport equation (2.36) would be falsified.","tokens_in":13823,"feed_emoji":"🌊","tokens_out":8367,"duration_ms":72703,"temperature":0.7,"pith_summary":"This paper develops the space-time ray version of the standard 'vertical modes and horizontal rays' method for sound propagation in shallow water. It works out the full set of equations—an eikonal equation for the phase, a transport equation for the amplitude, and the corresponding Hamiltonian ray equations—using only the squared vertical-mode eigenvalue $q^2(\\mathbf r,k_0)$. The central claim is that once this function is known, one can predict a signal's amplitude, frequency modulation, phase and pulse fronts, and space-time caustic positions at an observation point without solving the wave field for many frequencies separately. The reason this matters is practical: it turns a computationally heavy multi-frequency mode calculation into a ray-geometry reading of one dispersion-related function.","feed_headline":"Predict a pulse's shape from one waveguide eigenvalue","feed_subtitle":"A ray-based derivation turns the vertical-mode eigenvalue into amplitude, fronts, and frequency modulation.","key_machinery":"The central object is the eigenvalue function $q(\\mathbf r,k_0)$, the positive square root of the eigenvalue of the vertical Sturm-Liouville problem $\\psi''+n^2(\\mathbf r,z)k_0^2\\psi=q^2\\psi$ with boundary conditions at the surface and bottom. This one function carries all information about dispersion and horizontal variability: it defines the eikonal equation, the Hamiltonian $H=|\\mathbf k|^2-q^2$, the space-time ray velocity $\\hat{\\boldsymbol\\kappa}=(-\\partial q/\\partial k_0,\\mathbf k/q)$, and the group velocity $v=-(\\partial q/\\partial k_0)^{-1}$. The paper's distinctive machinery is the transport equation written as a three-dimensional divergence $\\widehat{\\mathrm{div}}(qA^2\\hat{\\boldsymbol\\kappa})=0$, together with the fundamental variation matrix $M$ that supplies the Jacobian $D(\\tau,\\mu,\\nu)$ used for amplitudes and for the caustic condition $D=0$.","core_discovery":"The paper establishes that a time-dependent acoustic field in a shallow-water waveguide, written as a sum of adiabatic vertical modes times a slowly varying amplitude and a fast phase $e^{i\\varepsilon^{-1}\\varphi}$, obeys a Hamiltonian ray system in the four-dimensional phase space $(\\tau,\\mathbf r,k_0,\\mathbf k)$. The phase satisfies the eikonal equation $q^2(\\mathbf r,\\partial\\varphi/\\partial\\tau) - |\\nabla\\varphi|^2 = 0$, where $q(\\mathbf r,k_0)$ is the square root of the vertical-mode eigenvalue, and the leading amplitude satisfies the transport equation $\\widehat{\\mathrm{div}}(q A^2 \\hat{\\boldsymbol\\kappa}) = 0$ in the three-dimensional space-time of $\\tau,\\mathbf r$. From these two equations the paper derives explicit formulas: the amplitude along a ray as $A(\\tau)=A(0)\\sqrt{g(0)D(0)/(g(\\tau)D(\\tau))}$, the phase advance $\\varphi(\\tau)=\\varphi(0)+k_0\\tau+\\int_0^\\tau q v\\,d\\tau'$, and the condition $D(\\tau,\\mu,\\nu)=0$ that locates space-time caustics. The paper's central claim is that the observable signal parameters at a given point—frequency modulation, amplitude, phase-front and pulse-front tilts—are all read off these ray objects, so a single function $q$ controls the propagation.","pith_inferences":["This formulation suggests a natural extension beyond non-caustic points: a Maslov-type or uniform asymptotic continuation of the amplitude formula would remove the $D^{-1/2}$ singularities at space-time caustics, which the paper does not attempt.","Because the derivation uses only the dispersion relation $q(\\mathbf r,k_0)$, the same space-time ray apparatus should transfer to electromagnetic pulses in dispersive media, where pulse-front tilt and angular dispersion are already known effects.","A direct numerical check would be feasible: compare the predicted $k_0(\\rho)$ and amplitude $A(\\tau)$ against a full-wave solution for a range-dependent wedge with a strongly frequency-dependent mode; the paper does not report such a comparison.","The weakest numerical step, the approximation (2.32), could be tested in isolation by comparing exact eigenfunction inner products with $(q/k_0)\\partial q/\\partial k_0$ over the frequency band of interest; if the discrepancy is small, the amplitude and frequency predictions are on solid ground."],"forward_implications":["In any waveguide where the eigenvalue function $q(\\mathbf r,k_0)$ is available, the method gives closed-form predictions for the frequency modulation $k_0(\\rho)$ and the phase-front normals at a specified observation point.","The vanishing of the Jacobian $D(\\tau,\\mu,\\nu)$ marks space-time caustics, so caustic surfaces and their projections can be located directly from the variation matrix $M$.","When several rays reach the observation point, the total field is a coherent sum of their contributions; instantaneous frequency and wave vector are no longer unique, and the model accounts for the resulting interference pattern.","In the constant-group-velocity limit the $\\tau$-front and the $s$-front coincide, and the amplitude along a ray decays as $1/\\sqrt{s}$, giving a simple characterization of the amplitude front.","All of these predictions depend on the model only through $q$, so the same equations apply to any shallow-water environment once the vertical eigenvalue problem is solved."],"supporting_citations":[{"why":"Introduces horizontal ray theory for ocean acoustics, the vertical-mode/horizontal-ray decomposition this paper extends to time-dependent space-time rays.","marker":"[21]"},{"why":"Provides the 'vertical modes and horizontal rays' formulation of waveguide propagation used throughout the derivation.","marker":"[6]"},{"why":"Defines the space-time ray method and the canonical Hamiltonian ray system on which the paper builds.","marker":"[1]"},{"why":"Supplies the caustics and catastrophe-theory framework for interpreting the vanishing Jacobian as space-time caustics.","marker":"[17]"},{"why":"Gives the geometrical-optics transport-equation and amplitude-divergence framework adapted here to dispersive waveguides.","marker":"[18]"},{"why":"Provides general properties of vertical modes and the physical setting of shallow-water acoustics used in Section 2.","marker":"[13]"}],"fun_headline_variants":["One eigenvalue, entire pulse forecast","Space-time rays: eigenvalue shapes the pulse","From one mode to full signal prediction","Eigenvalue to signal: a ray-based shortcut","Waveguide eigenvalue dictates pulse shape"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a certain weighted average of the squared refractive index over a vertical mode can be replaced by the simple expression $(q/k_0)\\partial q/\\partial k_0$ with no error bound. That approximation is what converts the transport equation into its final divergence form, and therefore everything the paper says about amplitude and frequency modulation inherits it. If vertical modes change noticeably with frequency, the replacement is not controlled and the predictions can be wrong.","fun_headline_variants_meta":{"raw":{"variants":["One eigenvalue, entire pulse forecast","Space-time rays: eigenvalue shapes the pulse","From one mode to full signal prediction","Eigenvalue to signal: a ray-based shortcut","Waveguide eigenvalue dictates pulse shape"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000914,"raw_usage":{"total_tokens":3922,"prompt_tokens":941,"completion_tokens":2981,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":557,"completion_tokens_details":{"reasoning_tokens":2918}},"tokens_in":557,"tokens_out":2981,"duration_ms":20452,"temperature":1.0,"reasoning_tokens":2918,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:26:06.254103+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check would be a numerical experiment in a range-dependent shallow-water waveguide: compute the exact vertical-mode inner product $\\langle n^2\\psi,\\psi\\rangle$ as a function of $k_0$, compare it with $(q/k_0)\\partial q/\\partial k_0$, and then propagate a broadband pulse with the ray formulas and with a full-wave solver. If the relative error in the inner product is not small, or if the ray-predicted amplitude and frequency modulation diverge from the full-wave result near a known caustic, the transport equation (2.36) would be falsified.","supporting_citations":[{"cited_title":"Horizontal ray the ory for ocean acoustics","cited_arxiv_id":null,"evidence_quote":"Introduces horizontal ray theory for ocean acoustics, the vertical-mode/horizontal-ray decomposition this paper extends to time-dependent space-time rays."},{"cited_title":"Horizontal rays and vertical modes","cited_arxiv_id":null,"evidence_quote":"Provides the 'vertical modes and horizontal rays' formulation of waveguide propagation used throughout the derivation."},{"cited_title":"Babich, Ivan A","cited_arxiv_id":null,"evidence_quote":"Defines the space-time ray method and the canonical Hamiltonian ray system on which the paper builds."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the caustics and catastrophe-theory framework for interpreting the vanishing Jacobian as space-time caustics."},{"cited_title":"Kravtsov and Yuri Ilich Orlov","cited_arxiv_id":null,"evidence_quote":"Gives the geometrical-optics transport-equation and amplitude-divergence framework adapted here to dispersive waveguides."},{"cited_title":"Fundamentals of shallow water acoustics","cited_arxiv_id":null,"evidence_quote":"Provides general properties of vertical modes and the physical setting of shallow-water acoustics used in Section 2."}],"review_version":1}