{"id":"f2e1c625-66fe-4946-92a5-35d4c066bbdb","arxiv_id":"2507.22557","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Direct simulation of plane-symmetric deep-water gravity waves shows a strongly nonlinear regime with omega^-4 and k^-4 spectra, shock fronts, and -7/2 steepness tails, with no weak-turbulence spectrum detected.","lead":"Direct numerical simulations of one-dimensional deep-water waves with random forcing and viscous damping find that the system always enters a strongly nonlinear regime controlled by wave breaking and sharp crests, with energy spectra falling as frequency to the minus fourth and wavenumber to the minus fourth.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No resolution/dissipation convergence study is given, so the reported k^-4 and omega^-4 spectra may be set by the ad hoc small-scale damping and N/3 filter rather than by genuine wave-breaking cusps.","rationale":"The paper's central positive claim is plausible and consistent with prior theory, but the proof depends on the measured spectrum being independent of numerical parameters. The authors explicitly state that enhanced dissipation and filtering are introduced to suppress singularity formation, so a convergence study is the minimum evidence needed to rule out a regularization-induced slope. The reader's CONDITIONAL verdict remains appropriate; this concern sharpens the condition but does not move the verdict. The absence of multiple realizations and error bars reinforces the conditionality but is not the single most load-bearing issue, which is the possible numerical origin of the -4 spectra.","tokens_in":6941,"tokens_out":6620,"duration_ms":89427,"concrete_test":"Repeat the main run at N = 65,536 with kd = 10,000 and the same gamma0; then vary gamma0 to 10^-6 and 10^-4 at fixed N and kd, and separately move the filter to N/4. For each run, fit the exponent of E_k on [20, 0.1 kd] and of E_omega on the corresponding omega interval. If the fitted exponent changes by more than about 0.2, or if the upper end of the -4 plateau tracks kd or the filter cutoff rather than staying fixed, the observed slope is not established as the physical inertial-range exponent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equations (7)-(8) include artificial dissipation gamma_k = gamma0 k^2 with gamma0 = 10^-5 and kd = 5000, and every time step applies a low-pass filter removing k >= N/3 ~ 1.1 x 10^4. The cusp singularities that are claimed to generate the predicted -4 spectra are curvature singularities; the stated purpose of the enhanced damping and the filter is to prevent exactly those singularities from forming overturning, multivalued shapes. If the numerical scheme smooths each incipient cusp on the scale set by kd or the filter, then the Fourier tail of a smoothed slope discontinuity will show a -4 law up to the smoothing scale and steepen beyond it. The observed plateau in Figs. 5-6 extends nearly to the filter cutoff, so the claimed 'inertial range' may be the regularization range. No run with different N, gamma0, kd, or filter location is reported, and no error bars are given. The negative weak-turbulence claim also rests on an unshown 'wide range of pumping amplitudes' scan.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports direct numerical simulations of plane-symmetric deep-water gravity waves in conformal variables, with random-phase low-wavenumber pumping and high-wavenumber viscous damping. The authors find that, rather than a weak-turbulence Kolmogorov-Zakharov spectrum, the system enters a strongly nonlinear quasi-steady state dominated by wave breaking and cusp formation. They report spatial and frequency spectra consistent with k^-4 and omega^-4 over about three decades, a linear omega ~ k branch in the spatiotemporal spectrum, and power-law tails |eta_x|^-7/2 in the steepness PDF. They interpret these as the first direct numerical confirmation of the strong-turbulence spectra predicted in [2].","tokens_in":7194,"tokens_out":4749,"duration_ms":55489,"significance":"If the reported spectra are genuine, this is a notable result: it gives a numerical test of a parameter-free strong-turbulence theory, links the -4 spectrum to wave-breaking cusps, and shows intermittency via a -7/2 steepness tail. The paper's central strength is that the observed exponents are not fitted and arise from direct integration of the standard potential-flow equations. However, the evidence is currently based on one run with strong artificial regularization, with no convergence study or amplitude scan, so the confirmation claim is not yet established. The result would be more convincing with a systematic study, but the core idea is defensible and the missing tests are within the scope of a revision.","major_comments":[{"comment":"The central claim that the measured k^-4 and omega^-4 spectra are the strong-turbulence spectra of [2] requires ruling out that they are imposed by the numerical regularization. The run uses N=32768, gamma0=1e-5, kd=5000, and a low-pass filter removing k>=N/3 ~ 1.1e4 at every time step, with gamma0 enhanced by three orders of magnitude for k>=kd; the apparent spectral range in Figs. 5-6 extends close to this cutoff. Because the curvature singularities responsible for the predicted -4 tail are exactly what the enhanced damping and filter prevent from forming, the observed slope may simply be the Fourier signature of events smoothed at the regularization scale. I request convergence runs with at least two higher resolutions, different kd and gamma0 values, and different filter cutoffs, plus a demonstration that the spectral exponent and its plateau width are unchanged outside the dissipation/filter range.","section":"Basic equations and numerical scheme (Eqs. (7)-(8); Figs. 5-6)"},{"comment":"The abstract and the text state that for a wide range of pumping amplitudes the weak turbulence regime was not detected, but the paper reports a single simulation configuration (F0=0.6e5, kf=3, L0=0.64). No amplitude scan, ensemble averaging, or error bars are shown, so the negative claim about weak turbulence and the claim of a wide amplitude range are not supported by the presented data. Either the scan should be shown (even in summary form), or the claims should be restricted to the single studied amplitude.","section":"Abstract and Simulation Results"},{"comment":"The paper states that the inertial interval spans almost three decades in both k and omega, but the figures do not include error bars, confidence bands, or a statement of how the fitting interval was chosen. Given that the spectra are time averages from one run and the PDF tail in Fig. 2 is also from a single trajectory, quantitative uncertainties are needed before the exponent -4 and the exponent -7/2 can be asserted as robust. At a minimum, the authors should report the time span used, the number of independent data points in the fitting ranges, and the sensitivity of the fitted exponents to the chosen intervals.","section":"Turbulence spectra (Figs. 5-6) and PDF (Fig. 2)"}],"minor_comments":[{"comment":"The phrase 'leas to form mutually ambiguous regions of the boundary shape' contains a typo; it should read 'leads to form multivalued regions of the boundary shape' or 'overturning regions.'","section":"Basic equations and numerical scheme"},{"comment":"The text uses 'power-low exponent' in the introduction; this should be 'power-law exponent.'","section":"Introduction and Concluding remarks"},{"comment":"The abscissa label appears truncated or unclear; please specify that the plotted variable is the steepness eta_x normalized by its standard deviation sigma_{eta_x}.","section":"Fig. 2"},{"comment":"The vertical axis label should be defined explicitly, for example as |eta_omega|^2 or the frequency spectrum E_omega, so that the reader can connect it to Eq. (2).","section":"Fig. 6"},{"comment":"The term 'shock fronts' is used for slope discontinuities; a sentence defining the terminology and distinguishing it from compressive shocks would improve clarity.","section":"Spatiotemporal analysis"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. Short numerical paper claiming that plane-symmetric deep-water gravity wave turbulence in a conformal-variable Euler solver settles into a strongly nonlinear regime dominated by wave breaking, with E_omega ~ omega^-4 and E_k ~ k^-4, linear omega~k dispersion, and -7/2 steepness PDF tails, and that weak-turbulence spectra are not observed. The genuinely new element is the direct numerical realization: shock fronts moving at nearly constant velocity, and spectra and tails that match the earlier theory in [2] and the experimental papers [16-18]. That is a real and useful observation, and the equations and numerical method are standard.\n\nThe soft spots are real but not fatal. The central one is the absence of any convergence study. The inertial range in Figs. 5-6 extends almost to the low-pass filter at k = N/3, and the enhanced damping for k >= kd = 5000 is three orders of magnitude larger than the background. The stated purpose of that regularization is to stop the cusps from overturning. Since the -4 law is generated by cusp-like slope discontinuities, the regularization could be setting the very scaling it is supposed to reveal. No run with a different N, gamma0, kd, or filter location is shown, so the stress-test concern is not answered by the paper.\n\nSecond, the paper says a 'wide range of pumping amplitudes' was probed and weak turbulence was never seen, but only one set of parameters is presented. The negative weak-turbulence result therefore rests on an unshown scan. No error bars or ensemble averages are reported anywhere, which is common for a Letter but still limits how much weight the exponents can carry.\n\nThird, the claim 'first direct numerical confirmation' is a little aggressive given the previous numerical literature on the same model, but the self-citation to [2] is legitimate: the theory is parameter-free and the simulation is not fitted to it.\n\nFor all that, the qualitative picture is credible and the paper is doing something real. The right next step is to give the authors a chance to run the convergence tests and show a few pumping amplitudes. I would send it to review, with a request for those additions. If the -4 slope survives resolution and damping changes, this becomes a solid confirmation; as it stands, it is a well-posed observation with a load-bearing missing control.","headline":"Plausible numerical evidence for strong-turbulence -4 spectra, but the missing resolution/damping convergence study leaves the scaling possibly set by the filter.","tokens_in":7712,"tokens_out":3094,"would_cite":true,"duration_ms":36119,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["47.35.-i","47.27.-i"],"model":"deepseek-v4-flash","headline":"Wave breaking sets the water-wave turbulence spectrum to slope -4","keywords":["water wave turbulence","wave breaking","cusp singularity","strong turbulence","conformal variables","power-law spectrum","intermittency","surface gravity waves"],"falsifier":"Run the same simulation with the dissipation threshold $k_d$ moved well below 5000 and well above it, and with the low-pass cutoff changed, and see whether the measured $k^{-4}$ and $\\omega^{-4}$ slopes and the $|\\eta_x|^{-7/2}$ tail survive; if the exponents move with the damping scale or filter, the spectrum is shaped by the numerics rather than by wave breaking itself. Alternatively, a wave-tank experiment with breaking-dominated waves that resolves both spatial and temporal spectra would directly test the linear $\\omega\\sim k$ ridge and the equal -4 exponents.","tokens_in":6720,"feed_emoji":"🌊","tokens_out":7358,"duration_ms":82031,"temperature":0.7,"pith_summary":"The paper claims that once plane-symmetric deep-water gravity waves are driven beyond a weakly nonlinear regime, the energy spectrum of the turbulent surface is controlled by wave breaking: sharp cusps form on crests and their sudden, random appearances imprint a universal $\\omega^{-4}$ law on the frequency spectrum and a $k^{-4}$ law on the wavenumber spectrum. Because both spectra share the same exponent, the effective dispersion relation in this state is linear, $\\omega \\sim k$, instead of the linear gravity-wave law $\\omega = \\sqrt{gk}$. These numerical results are put forward as the first direct confirmation of the strong-turbulence spectra predicted in [2], and they come with strongly non-Gaussian, power-law tails in the probability density of surface steepness that signal intermittency. The message matters because it separates the weak-turbulence picture, in which many resonant wave interactions transfer energy, from a strongly nonlinear picture in which actual breaking events set the statistical law.","feed_headline":"Wave breaking sets the water-wave turbulence spectrum to slope -4","feed_subtitle":"Simulations show breaking crest cusps dominate both spectra and make surface steepness strongly intermittent.","key_machinery":"The carrying object is the cusp singularity that forms on a wave crest immediately before breaking. At a fixed spatial point the surface height $\\eta(t)$ stays smooth until the breaking time $t_0$, when $d^2\\eta/dt^2$ acquires a delta-function spike; after averaging over many such random events with jump amplitudes $\\Gamma_i$ and occurrence rate $\\nu$, the singular part of the Fourier transform yields exactly $E_\\omega = g\\nu \\langle\\Gamma^2\\rangle/(2\\pi\\omega^4)$. Because the same cusp produces a jump in the slope $\\eta_x$ as a function of $x$, the same power law shows up in $k$-space, and the agreement of the two exponents implies $\\omega\\sim k$. The numerical scheme uses a conformal mapping of the fluid region onto a half-plane, reducing the free-boundary problem to two coupled equations for the surface coordinate and velocity potential, with random low-frequency forcing, viscous damping $\\gamma_0 k^2$, extra strong damping for $k\\ge k_d=5000$, and a low-pass filter removing $k\\ge N/3$ each step.","core_discovery":"In a one-dimensional, plane-symmetric model of deep-water surface gravity waves, written in conformal variables so the free surface is tracked exactly, the authors drive the system with random low-wavenumber forcing and damp high-wavenumber motion. For a wide range of pumping amplitudes they never see the weak-turbulence regime; instead the surface develops narrow shock-like fronts where the slope $\\eta_x$ jumps. At these fronts the second time derivative of the surface elevation behaves as a sum of random delta functions, and Fourier transforming that singular part gives the spectrum $E_\\omega \\propto \\omega^{-4}$; the same algebraic exponent $-4$ appears in the wavenumber spectrum $E_k$, and the spatiotemporal spectrum shows that disturbances move along straight lines $\\omega \\sim k$, not along the linear gravity wave dispersion curve. The authors claim these results are the first direct numerical confirmation of the strong-turbulence spectra predicted in [2], and that the measured probability density of steepness, with tails $\\propto |\\eta_x|^{-7/2}$, matches the intermittency expected for such cusp-dominated turbulence.","pith_inferences":["If the -4 slope is a genuine physical state, then the spectral collapse to $\\omega \\sim k$ suggests that in the inertial range the wave field behaves as a set of nearly non-dispersive shock fronts, so theories of strong turbulence that start from weak-wave interactions are missing the relevant degrees of freedom.","The numerical regularization could, in principle, set the cusp shape and therefore the high-frequency tail; a natural test is to vary the damping threshold $k_d$ and the low-pass cutoff and check whether the $-4$ exponents and the $-7/2$ PDF tails remain unchanged.","A laboratory experiment with breaking-dominated waves that measures both frequency and wavenumber spectra simultaneously could test the predicted $\\omega \\sim k$ ridge and the equality of the two -4 exponents, which would corroborate this mechanism beyond numerics."],"forward_implications":["If these spectra are the genuine strong-turbulence state, then in the plane-symmetric setup the classical Kolmogorov-Zakharov weak-turbulence cascade does not appear even at small pumping amplitudes, so weak-turbulence theory is not the right starting point for one-dimensional gravity waves.","An $\\omega^{-4}$ frequency law and a $k^{-4}$ wavenumber law with linear dispersion mean that measuring only the spectral slope cannot distinguish weak from strong turbulence; the linear $\\omega\\sim k$ signature is the distinguishing observable.","The same cusp mechanism should make strong water-wave turbulence highly intermittent, with rare large steepness events following $\\sim |\\eta_x|^{-7/2}$ tails, so wave-breaking statistics, not Gaussian fluctuations, dominate extreme events.","Because this behavior is shared with Burgers turbulence and with sound turbulence, the $\\omega^{-4}$/$k^{-4}$ result connects surface gravity waves to a broader family of nondispersive shock-dominated turbulent systems."],"supporting_citations":[{"why":"Supplies the predicted strong-turbulence spectra $E_\\omega\\propto\\omega^{-4}$ and $E_k\\propto k^{-4}$ that this paper's numerics confirm.","marker":"[2]"},{"why":"Gives the conformal-variable equations of motion used for the direct numerical simulation.","marker":"[5]"},{"why":"Companion derivation of the conformal formulation that underlies the numerical scheme.","marker":"[6]"},{"why":"Provides the weak-turbulence Zakharov-Filonenko spectrum whose $\\omega^{-4}$ frequency law coincides numerically with the strong-turbulence result.","marker":"[3]"},{"why":"Predicts the five-wave weak-turbulence spectrum for plane-symmetric gravity waves that the present simulations fail to observe.","marker":"[7]"},{"why":"Describes the procedure for reconstructing $\\eta(x,t)$ from conformal coordinates, which is used to compute the displayed spectra.","marker":"[9]"},{"why":"Documents the finite-time collapse and curvature singularity formation that motivates the cusp mechanism behind the spectra.","marker":"[13]"},{"why":"Reports laboratory measurements of the strong-turbulence spectrum that the present numerical results corroborate.","marker":"[16]"},{"why":"Supplies the Burgers-turbulence analogy for the $|\\eta_x|^{-7/2}$ power-law tails seen in the steepness distribution.","marker":"[24]"}],"fun_headline_variants":["Water-wave turbulence skips weak regime, cusps set -4 slope","Breaking crest cusps dominate water-wave spectra: -4 exponent","Simulations: strong turbulence governs water waves, no weak regime","Water-wave spectra follow -4 power law from cusp formation","Plane-symmetric water waves show strong turbulence, no weak scaling"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central claim assumes that the numerical damping applied to short waves and the low-pass filter used to prevent overturning wave shapes do not influence the measured spectral exponents; the paper does not report tests with different damping strengths or resolutions.","fun_headline_variants_meta":{"raw":{"variants":["Water-wave turbulence skips weak regime, cusps set -4 slope","Breaking crest cusps dominate water-wave spectra: -4 exponent","Simulations: strong turbulence governs water waves, no weak regime","Water-wave spectra follow -4 power law from cusp formation","Plane-symmetric water waves show strong turbulence, no weak scaling"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000624,"raw_usage":{"total_tokens":2855,"prompt_tokens":875,"completion_tokens":1980,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":491,"completion_tokens_details":{"reasoning_tokens":1889}},"tokens_in":491,"tokens_out":1980,"duration_ms":18255,"temperature":1.0,"reasoning_tokens":1889,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T11:32:13.305454+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same simulation with the dissipation threshold $k_d$ moved well below 5000 and well above it, and with the low-pass cutoff changed, and see whether the measured $k^{-4}$ and $\\omega^{-4}$ slopes and the $|\\eta_x|^{-7/2}$ tail survive; if the exponents move with the damping scale or filter, the spectrum is shaped by the numerics rather than by wave breaking itself. Alternatively, a wave-tank experiment with breaking-dominated waves that resolves both spatial and temporal spectra would directly test the linear $\\omega\\sim k$ ridge and the equal -4 exponents.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the predicted strong-turbulence spectra $E_\\omega\\propto\\omega^{-4}$ and $E_k\\propto k^{-4}$ that this paper's numerics confirm."},{"cited_title":"Dyachenko, E.A","cited_arxiv_id":null,"evidence_quote":"Gives the conformal-variable equations of motion used for the direct numerical simulation."},{"cited_title":"Dyachenko, V.E","cited_arxiv_id":null,"evidence_quote":"Companion derivation of the conformal formulation that underlies the numerical scheme."},{"cited_title":"Zakharov, N.N","cited_arxiv_id":null,"evidence_quote":"Provides the weak-turbulence Zakharov-Filonenko spectrum whose $\\omega^{-4}$ frequency law coincides numerically with the strong-turbulence result."},{"cited_title":"Dyachenko, Y.V","cited_arxiv_id":null,"evidence_quote":"Predicts the five-wave weak-turbulence spectrum for plane-symmetric gravity waves that the present simulations fail to observe."},{"cited_title":"Kochurin, JETP Lett., 118(12), 893-898 (2023)","cited_arxiv_id":null,"evidence_quote":"Describes the procedure for reconstructing $\\eta(x,t)$ from conformal coordinates, which is used to compute the displayed spectra."},{"cited_title":"Dyachenko, P.M","cited_arxiv_id":null,"evidence_quote":"Documents the finite-time collapse and curvature singularity formation that motivates the cusp mechanism behind the spectra."},{"cited_title":"Denissenko, S","cited_arxiv_id":null,"evidence_quote":"Reports laboratory measurements of the strong-turbulence spectrum that the present numerical results corroborate."},{"cited_title":"Gurbatov , A.N","cited_arxiv_id":null,"evidence_quote":"Supplies the Burgers-turbulence analogy for the $|\\eta_x|^{-7/2}$ power-law tails seen in the steepness distribution."}],"review_version":1}