{"id":"4790068e-b43b-4e17-bbdb-01af5fc7bbed","arxiv_id":"2508.19091","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Fractal-like families of stable, large-energy multi-mode periodic solutions are found in the cubic wave and beam equations via Galerkin continuation and reducible mode analysis.","lead":"The paper maps a fractal-like web of time-periodic vibration patterns for nonlinear wave and beam equations, including large-energy multi-mode solutions that are linearly stable. It uses a two-mode 'reducible system' approximation to predict where these solution branches appear and tests them against full numerics.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Beam-equation existence claim rests on an unverified N→∞ extrapolation: the reducible-system argument is approximate for the (1,1) branch, and no proof or convergence test shows the Galerkin branches persist for the PDE.","rationale":"Reading in good faith, the paper is transparent about its heuristic status: the N→∞ statement is phrased as 'this behaviour suggests', and the beam proof is explicitly deferred. The reducible-system algebra is clean, the qualitative agreement with Galerkin numerics is plausible, and the wave equation has independent computer-assisted support [24], which counts as real evidence. However, the strongest claim is made for both equations and includes large-energy existence and stability. The beam half lacks any rigorous or residual-convergence evidence, and the non-reducible (1,1) beam branch means the analytic branch predictor is approximate precisely where the numerical example starts. The M=N² vs. m≈2N² discrepancy is a concrete additional check on the completeness of the numerics. These are addressable gaps rather than demonstrated failures, so I do not move the verdict away from CONDITIONAL. The condition should be: prove or certify beam branch persistence in the PDE limit.","tokens_in":7025,"tokens_out":14390,"duration_ms":158142,"concrete_test":"Certify existence for the beam equation by adapting the computer-assisted proof in [24] (arXiv:2506.10839) to ν=2: take the first nontrivial branch (m,n)=(1,1) with the full non-reducible Galerkin system (10) at increasing N, and use interval arithmetic to verify a contraction/Newton–Kantorovich argument with controlled truncation error. If the certification succeeds, the N→∞ extrapolation is supported; if it cannot be completed, the beam-equation claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The strongest claim asserts new large-energy, multi-mode time-periodic solutions for both the wave and beam equations. For the wave equation the cited computer-assisted proof [24] supplies independent support. For the beam equation (ν=2) no such proof is given; the Summary states only that 'an extension of the proof to the beam equation appears straightforward.' That is precisely the missing load-bearing step. The analytical prediction comes from reducible systems (Eq. 6), but for the first nontrivial beam branch (m,n)=(1,1) the full two-mode system (Eq. 10) contains the extra terms −3A²B and −A³, so the reducible prediction is only approximate even at the two-mode level. The finite-Galerkin evidence is also incomplete: for ν=2, the two-mode existence condition (2m+1)<(2n+1)² permits time modes up to m≈2N² for the largest spatial mode n=N−1, while the reported computations use M=N² (Fig. 2); the assertion that increasing M adds no new branches is not demonstrated. Thus the beam-equation half of the central existence claim depends on an unverified limit and an approximate model, with no residual-convergence or computer-assisted certificate supplied.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies time-periodic solutions of the one-dimensional defocusing cubic wave and beam equations (ν=1,2 in Eq. (1)) with Dirichlet/Navier boundary conditions. It uses a Galerkin ansatz with odd temporal and spatial modes, numerically continues a family of solutions from the linearized fundamental mode (the 'trunk'), and observes additional solution branches. These branches are modeled by finite 'reducible' two-mode and three-mode algebraic systems, Eq. (6) and its analogues, which correctly reproduce branch locations and shapes for both equations in the reported examples. The paper argues, by extrapolation from these finite reductions, that in the N→∞ limit the branches populate the trunks densely, implying new large-energy time-periodic solutions with frequencies arbitrarily close to Ω=1. It also performs Floquet stability computations on the numerical solutions and reports regions of linear stability on the trunk and on branches. A computer-assisted proof for the wave equation is cited [24]; for the beam equation the paper only states that a similar proof 'appears straightforward.'","tokens_in":7384,"tokens_out":8090,"duration_ms":79027,"significance":"If the finite-dimensional branch structure and stability computations are taken at face value, the paper provides a useful and systematic map of a complex solution web for two classical PDEs, and the explicit two-mode formulas (7)-(8) are a clean analytical tool. The numerical continuation and Floquet analysis are appropriate methods, and the paper is honest in presenting the beam-equation PDE-level statement as an extrapolation. The main advertised conclusions, however, go beyond the finite-Galerkin evidence, especially for the beam equation, where no rigorous certificate or quantitative convergence study is supplied. The wave-equation half is substantially supported by the cited computer-assisted proof [24], which is a genuine strength. The value of the paper would be much increased if the PDE-level status of the beam claims were clarified and the numerical convergence claims made precise.","major_comments":[{"comment":"The statement 'in the limit N→∞ the branches populate the trunks densely' and the accompanying implication about solutions with arbitrarily large energies and frequencies arbitrarily close to Ω=1 are extrapolated from a finite two-mode algebraic model. Eq. (8) describes solutions of Eq. (6), not of the PDE (1); no estimate of omitted modes or compactness argument is provided. For the wave equation the cited computer-assisted proof [24] supplies independent support, but for the beam equation the Summary's assertion that an extension 'appears straightforward' is not a proof. This is the load-bearing step for the beam-equation existence claim and needs either a proof or an explicit reframing as a numerical conjecture.","section":"Reducible systems, paragraph after Eq. (8)"},{"comment":"The claim that 'further increasing M does not qualitatively alter the solution structure; in particular, no new branches appear' is not demonstrated. For ν=2, the reducible-system existence condition (2m+1)<(2n+1)^2 permits, for the largest spatial mode n=N−1, temporal mode numbers m up to O(N^2), whereas the reported computations use M=N^2. No check is shown that branches with temporal modes beyond N^2 are absent, nor is any residual or convergence measure reported. This affects the completeness of the branch census and the basis for the N→∞ description.","section":"Galerkin scheme, second paragraph and Fig. 2"},{"comment":"For the first nontrivial beam branch (m,n)=(1,1), the exact two-mode Galerkin system contains the nonreducible terms −3A^2B and −A^3, so Eq. (6) is only an approximation. The Supplemental comparison in Fig. 6 is qualitative; no quantitative error bound is given. Since the (1,1) beam branch is the first and most prominent example of the claimed new class, the approximate nature of the reducible model weakens the analytic prediction of branch existence and location for the beam equation unless the omitted terms are shown to be harmless in a controlled sense.","section":"Supplemental Material, Eq. (10)"},{"comment":"The Floquet stability computations are presented without a convergence study in the spatial truncation K or in the temporal mode count M, and without a comparison against an independent numerical method. The sentence 'M large enough for the residue to be small' does not define the residue or give values. Since the existence of linearly stable members is one of the paper's advertised findings, this numerical certificate should be supplied or the claim weakened accordingly.","section":"Linear stability, paragraph starting 'We next investigate' and Figs. 4-5"}],"minor_comments":[{"comment":"The scaling symmetry is stated without explaining how the boundary conditions and the spatial domain are affected. Adding one sentence on why the transformed functions still satisfy the same boundary conditions would improve readability.","section":"Eq. (3)"},{"comment":"The term 'reducible systems' is defined only by reference to [20]. Since it is central to the paper, a short self-contained definition would help readers and would make the algebraic structure of Eq. (6) clearer.","section":"Before Eq. (6)"},{"comment":"The caption does not specify the meaning of the panels (e.g., which N values are shown) or the notation M=N^2. These details are scattered in the text; a clearer caption would avoid ambiguity.","section":"Fig. 2"},{"comment":"The final sentence states that one can show that the only non-reducible two-mode system producing a branch for the beam equation is the (1,1) case, but no proof or reference is given. This is a nontrivial claim and should either be proved or explicitly deferred.","section":"Supplemental Material"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest in presenting the numerical character of most of its evidence, and the explicit algebraic reductions are a strong point. The main concern is the gap between the finite-Galerkin computations and the PDE-level claims for the beam equation: the paper itself contains a self-identified missing proof ('extension appears straightforward'), which is a load-bearing point. This can be fixed within the manuscript's scope by either adding a convergence theorem or computer-assisted proof, or by carefully reformulating the beam claims as conjectures supported by numerics. The wave-equation case, backed by the cited computer-assisted proof, is in better shape."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Jim,\n\nQuick take: the paper is worth reading. It maps a genuinely intricate web of time-periodic solutions for the cubic wave and beam equations on an interval, mostly by numerical continuation plus an explicit two-mode “reducible” analysis, and it adds the first Floquet stability results for these branches. The wave-equation branch structure and reducible framework come from their earlier papers, but the stability analysis and the beam-equation extension are new. If correct, it means that between the known small-amplitude Cantor families there are large-energy multi-mode solutions, some linearly stable, which is a meaningful addition to the picture.\n\nWhat I like. The reducible-system derivation is clean: for most two-mode interactions, the Galerkin equations reduce to a linear system in squared amplitudes, giving closed-form formulas for branch locations. The numerics match those formulas well. They also openly flag the one non-reducible two-mode case, (m,n)=(1,1) for the beam equation, in the supplement — the actual system has extra cubic terms, and the reducible system is only approximate there. That is honest. For the wave equation they point to a computer-assisted proof [24], which is a real anchor. The stability computations are new, and the claim that some branch solutions are linearly stable is plausible.\n\nWhere I’d press. The beam-equation half of the main existence claim is a conjecture, not a theorem. The summary says extending the computer-assisted proof “appears straightforward” — that is not a proof. The N→∞ limit is inferred from finite-Galerkin data and a statement that increasing M does not add branches, but no convergence or residual bounds are shown. This matters because for (m,n)=(1,1) the reducible prediction is only approximate, so the chain from Eq. (6) to the PDE is weakest exactly on the first non-trivial beam branch. The stress-test note lands on this. The Floquet results also lack convergence data; the 1e-12 tolerance is fine, but I’d want to see K-dependence.\n\nThat said, the paper does not oversell: the Fig. 2 caption says “suggests”, the summary says “appears”, and the wave-equation proof is cited. The softness is real but is the standard gap between numerical discovery and rigorous confirmation. The reader’s CONDITIONAL verdict is fair.\n\nMy call: send it out. A serious referee can ask for convergence data and a clearer “what is proven vs. inferred” statement, and the numerics-plus-framework contribution is sturdy enough to be published once those caveats are addressed. I’d cite it for the beam-equation results and the stability findings.","headline":"Solid numerical and semi-analytic map of a multi-mode branch web for the cubic wave and beam equations; the genuinely new value is the Floquet stability analysis and the beam extension, but the beam existence claim is a conjecture, not a proof.","tokens_in":7783,"tokens_out":3239,"would_cite":true,"duration_ms":34125,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B10","35L70","35Q74","35B32","37K55"],"pacs":[],"model":"deepseek-v4-flash","headline":"The defocusing cubic wave and beam equations have a new class of large-energy, multi-mode time-periodic solutions that fill the gaps of previously known Cantor-like families, connect their rescaled copies, and include linearly stable member","keywords":["time-periodic solutions","nonlinear wave equation","beam equation","Galerkin method","reducible systems","Floquet stability","bifurcation structure","fractal pattern"],"falsifier":"Increase the truncation N for a fixed branch of the beam equation and monitor its energy–frequency curve and endpoint: if the branch endpoint does not approach the predicted intersection with the rescaled trunk as N grows, or if a branch labelled stable develops a Floquet multiplier with |λ|>1 once further modes are included, the claimed structure collapses. For the wave equation, the same test is settled by the computer-assisted proof in [24].","tokens_in":6958,"feed_emoji":"🎻","tokens_out":6436,"duration_ms":55227,"temperature":0.7,"pith_summary":"The paper claims that the defocusing cubic wave equation and beam equation on an interval admit far more time-periodic solutions than previously mapped: in addition to the known single-mode, Cantor-like families, there is a hierarchy of large-energy, multi-mode solution branches. These branches emerge from the main solution trunk, exchange dominance between modes, and reconnect at bifurcation points to rescaled copies of the trunk, producing a fractal-like web. The paper proposes a systematic 'reducible systems' framework—small Galerkin truncations solvable in closed form—that predicts the location, shape, and mode content of every branch. A Floquet analysis shows that some members of the new class are linearly stable, so these solutions may shape long-time nonlinear dynamics rather than being transient numerical artefacts. If the picture holds, periodic solutions of these mechanical model equations are organized by a self-similar structure that bridges small-amplitude rigorous results and the finite-amplitude regime.","feed_headline":"New class of stable oscillations found in strings and beams","feed_subtitle":"Large-energy multi-mode solutions fill the gaps of known periodic families and may shape long-time dynamics.","key_machinery":"The argument rests on three pieces: a Galerkin ansatz using only odd spatial and temporal modes, with pseudo-arclength continuation to follow solution curves through folds and bifurcations; the scaling symmetry u(τ,x) → n^ν u(mτ, nx), Ω → n^ν Ω/m, E → n^{4ν} E, which generates rescaled copies of every solution family and forces their connections at branch endpoints; and 'reducible systems'—two- and three-mode truncations whose algebraic equations are linear in squared mode amplitudes, solvable explicitly, giving a systematic catalogue of branch locations and mode compositions, including the non-reducible (m,n)=(1,1) case for ν=2. Floquet theory then tests linear stability of the numerically","core_discovery":"For Eq. (2) with defocusing cubic nonlinearity and Dirichlet (wave) or Navier (beam) boundary conditions, the authors identify a new class of time-periodic solutions whose Galerkin mode composition is dominated by two modes, cos τ sin x and cos(2m+1)τ sin(2n+1)x, with (2m+1)<(2n+1)^ν. Along such a branch the fundamental-mode amplitude decreases while the higher-mode amplitude grows; at the branch end the fundamental mode is absent and the solution coincides with a rescaled copy of the trunk. The two-mode systems are 'reducible'—linear in the squared amplitudes—and their explicit solutions reproduce all numerically observed branches for both equations, with one specially handled non-reducible","pith_inferences":["If the dense-population picture is right, the set of time-periodic solutions may be dense in the energy–frequency plane, which would make exact periodic states a common, not exceptional, feature of these nonlinear PDEs.","The stable branches suggest that Hamiltonian phase space contains islands of periodic or quasi-periodic motion at large energy; numerical long-time evolution of the PDE should show intermittent trapping near these states.","The same reducible-system logic may apply to other odd-power nonlinearities or to nonlinear field theories on bounded domains, where analogue fractal families of periodic solutions could be searched for at large amplitude.","A direct test would be to measure, in a nanomechanical or suspension-bridge setting, whether frequency–response curves show the predicted dense branching pattern at large driving amplitudes."],"forward_implications":["The previously established Cantor-like families of small-amplitude periodic solutions are only the visible part of a much larger, self-similar solution set; the gaps in frequency are filled by higher-energy multi-mode solutions.","Linearly stable large-energy periodic solutions exist, so these states can influence long-time evolution of the wave and beam equations, potentially acting as organized, long-lived oscillations.","The reducible-system framework provides a predictive catalogue of branch structure at any truncation order, making it possible to locate new solutions without solving the full Galerkin system.","The results extend to rescaled copies at arbitrarily large energies and to frequencies arbitrarily close to Ω=1, so periodic solutions with extreme parameters are expected.","For the wave equation a computer-assisted existence proof is already available; the authors state the extension of such a proof to the beam equation appears straightforward."],"supporting_citations":[{"why":"Computer-assisted proof of existence of the new class of time-periodic solutions for the wave equation; grounds the claim that these branches are genuine PDE solutions.","marker":"[24]"},{"why":"Introduces the reducible systems and the 'trees, trunks, and branches' classification used to predict the whole branch structure.","marker":"[20]"},{"why":"Supplies the numerical Galerkin scheme and continuation details for the wave equation that the present paper extends to the beam equation.","marker":"[19]"},{"why":"Establishes the Cantor-like families of periodic solutions for wave equations that the new class is said to fill and connect.","marker":"[9]"},{"why":"Establishes small-amplitude Cantor families for the nonlinear beam equation, the prior result the beam branches populate.","marker":"[16]"},{"why":"Provides the computer-assisted approach for Hamiltonian PDE periodic solutions that inspired the existence proof in [24].","marker":"[14]"},{"why":"Supplies Floquet theory used for the linear stability analysis of the constructed periodic solutions.","marker":"[21]"},{"why":"Pseudo-arclength continuation method used to parametrize solution branches through folds and bifurcations.","marker":"[18]"}],"fun_headline_variants":["Stable high-energy oscillations found in strings and beams","Two-mode solutions explain stable branches in strings and beams","Fractal-like structure reveals stable oscillations in strings and beams","New large-energy solutions linearly stable in strings and beams","Stable periodic branches discovered in nonlinear strings and beams"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The finite-Galerkin branch structure, computed at increasing truncations and modelled by two-mode reductions, survives as genuine solutions of the full partial differential equations in the limit of infinitely many modes—for the beam equation this is asserted rather than proved.","fun_headline_variants_meta":{"raw":{"variants":["Stable high-energy oscillations found in strings and beams","Two-mode solutions explain stable branches in strings and beams","Fractal-like structure reveals stable oscillations in strings and beams","New large-energy solutions linearly stable in strings and beams","Stable periodic branches discovered in nonlinear strings and beams"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000829,"raw_usage":{"total_tokens":3380,"prompt_tokens":590,"completion_tokens":2790,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":334,"completion_tokens_details":{"reasoning_tokens":2712}},"tokens_in":334,"tokens_out":2790,"duration_ms":21013,"temperature":1.0,"reasoning_tokens":2712,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T15:56:06.852934+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Increase the truncation N for a fixed branch of the beam equation and monitor its energy–frequency curve and endpoint: if the branch endpoint does not approach the predicted intersection with the rescaled trunk as N grows, or if a branch labelled stable develops a Floquet multiplier with |λ|>1 once further modes are included, the claimed structure collapses. For the wave equation, the same test is settled by the computer-assisted proof in [24].","supporting_citations":[{"cited_title":"New class of time- periodic solutions to the 1D cubic wave equation","cited_arxiv_id":null,"evidence_quote":"Computer-assisted proof of existence of the new class of time-periodic solutions for the wave equation; grounds the claim that these branches are genuine PDE solutions."},{"cited_title":"Periodic solutions for the 1D cubic wave equation with Dirichlet bound- ary conditions","cited_arxiv_id":null,"evidence_quote":"Supplies the numerical Galerkin scheme and continuation details for the wave equation that the present paper extends to the beam equation."},{"cited_title":"Cantor families of periodic solutions for wave equations via a variational principle","cited_arxiv_id":null,"evidence_quote":"Establishes the Cantor-like families of periodic solutions for wave equations that the new class is said to fill and connect."},{"cited_title":"KAM for the nonlinear beam equation 1: small-amplitude solutions","cited_arxiv_id":"1412.2803","evidence_quote":"Establishes small-amplitude Cantor families for the nonlinear beam equation, the prior result the beam branches populate."},{"cited_title":"Families of Periodic So- lutions for Some Hamiltonian PDEs","cited_arxiv_id":null,"evidence_quote":"Provides the computer-assisted approach for Hamiltonian PDE periodic solutions that inspired the existence proof in [24]."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Pseudo-arclength continuation method used to parametrize solution branches through folds and bifurcations."}],"review_version":1}