{"id":"dd4b87c1-2d29-405d-bbe7-3d2e0aee26f4","arxiv_id":"2506.10839","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"For Ω=69/40, three distinct 2π-periodic solutions of u_tt - u_xx + u^3 = 0 on [0,π] are rigorously shown to exist close to explicitly given rational approximations.","lead":"This paper proves, with computer-assisted exact arithmetic, that three distinct time-periodic solutions of the 1D cubic wave equation exist at frequency 69/40. Two of them belong to a new branch structure beyond the standard family, supporting the authors' earlier numerical conjectures.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Computational certificate is the sole load-bearing risk; tight margins in (2.6) make a single script or formula error fatal to Theorem 1.","rationale":"The reader's weakest assumption correctly identifies the computer-assisted verification as the only meaningful risk. I read the analytic sections carefully: the functional setting, the contraction argument, and the operator norm estimates (Lemmas 1-4 and formula (5.2)) are coherent, and the proof of Theorem 1 from Theorem 2 is straightforward. The remaining uncertainty is purely computational: whether the supplied Mathematica scripts, running in exact rational arithmetic, actually verify inequalities (2.6) and (2.7) with the rational numbers stated in Section 7. Nothing in the manuscript provides an independent audit of those scripts or a frozen certificate; the absence of a commit hash or checksum is a genuine reproducibility weakness. The margins are tight enough (on the order of 10^-5 for the first two solutions) that even a small error in formula (B.6), in the g-coefficients, or in the data files would destroy the proof. The proposed concrete test — an independent exact-arithmetic recomputation of ||H0|| from the provided data — would settle this concern: if it reproduces the stated bound, the certificate is valid and Theorem 1 follows; if not, the proof collapses. I therefore agree with the reader's conditional verdict and recommend no change.","tokens_in":15178,"tokens_out":17713,"duration_ms":175158,"concrete_test":"Recompute the operator-norm bound ||H0|| for u0^(1) using an independent exact-arithmetic implementation (e.g., Python fractions or Sage) from the supplied data files u0hat_1st.m and Acal_1st.m, implementing (B.6)-(B.7) with rho_tau=rho_x=1+10^-20 as in CAPCode.m, and compare the result exactly to 10426318695/10439379622. If the values differ, or if the calculation cannot be completed exactly, the certificate fails and the proof of Theorem 1 is incomplete.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The analytic framework is sound: Theorem 1 follows from Theorem 2 by a standard contraction argument, and the operator bounds (5.2), (B.6)-(B.7) are plausible. However, the proof of Theorem 2 is entirely delegated to Mathematica scripts using exact rational arithmetic, and the paper does not freeze the certificate (no commit hash, no checksums). The published rational bounds in Section 7 are extremely tight: for the first two solutions, K0 - ||H0|| is about 5.2e-5 and 4.1e-5 respectively, comparable in magnitude to the delta-dependent terms in (2.6). Any undetected error in the scripts — especially in the evaluation of (B.6)/(B.7), in the g-coefficients from (B.4), or in the assembly of the rational matrix A — could make (2.6) fail while the printed bounds appear to hold. The paper also notes that A is constructed by approximating a floating-point inverse, yet the run scripts do not explicitly verify that A is invertible, a hypothesis of Theorem 2. The text offers no independent check, so the existence claim rests on unverified code and data.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proves the existence of three nontrivial 2π-periodic solutions to the defocusing cubic wave equation on an interval with Dirichlet boundary conditions, all at frequency Ω = 69/40. The proof combines a weighted ℓ1 fixed-point framework with a computer-assisted verification of operator bounds using exact rational arithmetic. One solution is close to a low-mode Galerkin approximation; the other two are close to higher-energy approximations, supporting the authors' earlier numerical conjectures about 'branches' beyond the classical 'trunk' family. The analytic part is self-contained and reduces the existence claim to Theorem 2, whose hypotheses are verified by supplied Mathematica scripts and data files.","tokens_in":15440,"tokens_out":16686,"duration_ms":174940,"significance":"If the computational certificate is correct, the result is a valuable contribution to the rigorous theory of time-periodic solutions of resonant 1D nonlinear wave equations, providing the first proof of solutions outside the classical Cantor-family trunk. The analytic framework is clean: Lemmas 1–4 give transparent bounds on multiplication, the inverse of L_Ω, and tail estimates, and the reduction to the finite verification (2.6) is rigorous. The use of exact rational arithmetic is a sound choice for computer-assisted proofs and avoids rounding-error concerns. The authors are to be commended for supplying the code and data. The main risk is that the existence claim rests entirely on the correctness of unverified scripts with very tight margins.","major_comments":[{"comment":"The proof of Theorem 2 is delegated to the Mathematica scripts, but the manuscript does not include the output logs or an explicit evaluation of the inequalities (2.6). The margins are extremely tight: for the first and second solutions, K0 − ||H0|| is approximately 5.2×10^-5 and 4.1×10^-5, respectively, which is of the same order as the δ-dependent terms in (2.6). A small error in the scripts' evaluation of (B.6)–(B.7) or in the g-coefficients from (B.4) could flip the inequality and invalidate Theorem 1. Please provide a machine-readable certificate of all rational bounds and the final inequalities, including exact values of ||N(0)|| and of the left-hand side of the first inequality in (2.6), and freeze the exact code and data versions (e.g., with checksums).","section":"Section 7 and Appendix A"},{"comment":"Theorem 2 assumes that A is a linear isomorphism, but the construction of A as a rationalized approximate floating-point inverse of à does not include a proof that the resulting matrix A (and hence the block-diagonal operator A) is invertible. The text says the matrix is 'sufficiently close' to an inverse, but no exact determinant or explicit inverse is provided. Please supply an exact rational determinant or an explicit inverse for each A^(i), or alternatively modify Theorem 2 and the proof of Theorem 1 to remove the invertibility hypothesis, since the contraction argument itself only requires A to be a bounded linear operator.","section":"Section 4"},{"comment":"The convolution formulas (B.2)–(B.4) are central to assembling the matrix A and computing the bound on ||H0||, yet they are stated without derivation. An off-by-one error or a sign error in these formulas would propagate directly into the claimed bounds. Please add a derivation or a reference to one, and include a validation script that checks (B.2)–(B.4) against direct symbolic trigonometric products for small random inputs, so that the correctness of these formulas is independently verifiable.","section":"Appendix B"}],"minor_comments":[{"comment":"The contraction argument is applied on the open ball Bδ(0), which is not a complete metric space. Since the inequalities in (2.6) are strict, the map sends the closed ball into itself; the proof should be formulated on the closed ball to apply the Banach contraction principle.","section":"Section 2.3"},{"comment":"In the final line of the proof, the last sum is written with |\\hat u_{m1,n1}|; it should be |\\hat w_{m3,n3}| to match the product of the three norms.","section":"Lemma 1 proof"},{"comment":"In the display after the first inequality, the supremum is written with P_{m1,n1} and the denominator ρ^{2m+1}ρ^{2n+1}; the subscripts m1,n1 should be m,n for consistency.","section":"Lemma 2 proof"},{"comment":"The ε^(i) values are given only as decimal expansions in (7.1). Since the paper emphasizes exact rational arithmetic, please provide rational upper bounds for these quantities or state explicitly that the decimals are rigorous upper bounds computed from the rational data.","section":"Section 7"},{"comment":"Please specify the exact version of Wolfram Mathematica used, the operating system, and the hardware environment, and provide checksums for the data files and scripts to allow exact reproduction.","section":"Appendix A"},{"comment":"The formula for φ(m,n) is typeset ambiguously in the text due to line breaking; write it as a single fraction, \\frac{4q^2}{2\\max(2q(2n+1),(2p+1)(2m+1))-1}, to avoid confusion.","section":"Lemma 3"}],"recommendation":"major_revision","confidential_remarks":"The central claim is plausible and the analytic framework is sound, but the computational verification is the sole load-bearing point and it is not yet presented in a form that allows a referee to trust it without re-running the code. The tight margins make this a genuine risk. The issues are fixable: provide output logs, exact intermediate values, a check of invertibility of A, and validation of the convolution formulas. I would not recommend rejection, because the analytic reduction is rigorous and the computational steps are, in principle, verifiable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new thing here is that two solutions on the conjectured branch-family of the 1D cubic wave equation are now proven to exist, not just numerically indicated. The analytic setup is genuinely solid: the weighted l1 space, the contraction argument, and the lemmas bounding the nonlinearity and the inverse of the linear operator are all self-contained and correct. Lemma 1 (the product bound) and Lemma 3 (the rational-frequency inverse bound) are particularly clean. The paper also does something rare: it is honest about what is and is not proved. Remarks 1 and 2 put the single-frequency result in context and do not overclaim. The real load-bearing part is the computer verification of Theorem 2. The supplied Mathematica scripts and data files make the exact rational arithmetic reproducible in principle, and the formulas in Appendix B are detailed enough to reimplement. That is exactly what a computer-assisted proof should provide. Where I would push back is the tightness of the margins. For the two branch solutions, K0 - ||H0|| is about 5e-5. That is not a lot of room. If any script contains an error - especially in the g-coefficients from (B.4) or the norm evaluation in (B.6)-(B.7) - the printed rational bounds could be wrong and the inequality (2.6) could fail. The paper would be stronger if the certificate were frozen (checksums or a commit hash) and if the scripts explicitly verified that A is invertible, as Theorem 2 assumes. I suspect invertibility is not actually needed for the fixed-point argument to produce a solution, since the fixed point of N still gives u = u0 + Ah, but the paper asserts it and should check it. I have not rerun the scripts. If I were the referee, I would ask the authors to confirm the run environment, provide a frozen version of the code, and ideally get an independent reimplementation of the norm computation. That said, I see no red flag that suggests the computation is wrong; the risk is the usual one for computer-assisted proofs, and it is properly disclosed. This paper deserves a serious referee. It is an honest, technically sound advance in a small but active subfield. The contribution is modest but real, and the methods are worth having on the record. Send it to review, but make sure the referee actually engages with the computer-assisted part.","headline":"A clean analytic framework with a tight, code-delivered computational certificate; the existence proof is believable but the margins are thin enough that the scripts deserve a careful referee.","tokens_in":731,"tokens_out":1204,"would_cite":true,"duration_ms":78050,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B10","68V05","35B32","35L71"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper rigorously constructs three distinct time-periodic solutions of the defocusing cubic wave equation on an interval, two of them belonging to a new 'branch' class beyond the classical trunk family.","keywords":["Time-periodic solutions","Nonlinear wave equation","Bifurcations","Computer-assisted proof","Rational arithmetic","Fixed-point argument","Cubic wave equation","Dirichlet boundary conditions"],"falsifier":"Rerun the supplied scripts with the stated data (or independently recompute the maxima in formula (5.2) using interval arithmetic) and check whether, for each $i$, the quantity $\\|H_0\\|+6\\|L_\\Omega^{-1}\\|\\|u_0^{(i)}\\|\\|A\\|^2\\delta+3\\|L_\\Omega^{-1}\\|\\|A\\|^3\\delta^2$ is smaller than $K_0$ and $\\|N_\\Omega(0)\\|$ is smaller than $(1-K_0)\\delta$. If any of these inequalities fails, the contraction argument does not close and the claimed solutions are not established.","tokens_in":14966,"feed_emoji":"🌊","tokens_out":14656,"duration_ms":149207,"temperature":0.7,"pith_summary":"The paper proves that the one-dimensional defocusing cubic wave equation $u_{tt}-u_{xx}+u^3=0$ with Dirichlet boundary conditions has at least three pairwise distinct $2\\pi$-periodic solutions, all with the same frequency $\\Omega=69/40$ after time rescaling. One of these sits on the classical 'trunk' family of solutions known since the 1980s, while the other two have larger energy and stronger high-mode content, and belong to a conjectured 'branch' structure. The proof is a computer-assisted fixed-point argument: it starts from approximate solutions $u_0^{(i)}$ with rational Fourier coefficients, constructs a finite-dimensional correction operator $A$, and verifies with exact rational arithmetic that a contraction mapping on a weighted $\\ell^1$ space has a fixed point close to each $u_0^{(i)}$. If correct, this is the first rigorous confirmation that the solution set of this equation is richer than the classical Cantor-like family.","feed_headline":"New class of periodic wave solutions proven to exist","feed_subtitle":"Rigorous proof confirms branch solutions beyond the classical trunk family.","key_machinery":"The argument works in the Banach space $X$ of $2\\pi$-periodic functions spanned by $P_{m,n}(\\tau,x)=\\cos((2m+1)\\tau)\\sin((2n+1)x)$, with norm $\\|v\\|=\\sum \\rho_\\tau^{2m+1}\\rho_x^{2n+1}|\\hat v_{m,n}|$ and weights $\\rho_\\tau=\\rho_x=1+10^{-20}$. The load-bearing object is the linear part $H_0(h)=-3L_\\Omega^{-1}(u_0^2 A h)+h-Ah$ of the fixed-point map $N_\\Omega(h)=F_\\Omega(u_0+Ah)-u_0+(I-A)h$, where $F_\\Omega(u)=-L_\\Omega^{-1}u^3$ and $A$ is a finite-dimensional rational matrix approximating the inverse of $I+3L_\\Omega^{-1}\\Lambda_{u_0^2}$. Lemma 3 gives the explicit bound $\\|L_\\Omega^{-1}v\\|\\le \\phi(m,n)\\|v\\|$ with $\\phi(m,n)=4q^2/(2\\max(2q(2n+1),(2p+1)(2m+1))-1)$, and the choice $\\Omega=(2p+1)/(2q)$ makes the denominators differences of an even and an odd integer, so the small-divisor problem---the near-resonances that plague generic perturbative constructions---does not arise. Lemma 4 bounds the action of $L_\\Omega^{-1}(u_0^2P_{m,n})$ on high modes, and formula (5.2) reduces $\\|H_0\\|$ to a finite maximum that the supplied scripts evaluate exactly in rational arithmetic.","core_discovery":"For frequency $\\Omega=69/40$, the paper's Theorem 1 asserts that for each of three explicitly given rational-coefficient approximate solutions $u_0^{(i)}$, $i=1,2,3$, there is an exact solution $u^{(i)}$ of the rescaled equation satisfying $\\|u^{(i)}-u_0^{(i)}\\|<\\varepsilon^{(i)}$ in a weighted $\\ell^1$ norm, with $\\varepsilon^{(1)}\\approx 1.79\\times 10^{-8}$, $\\varepsilon^{(2)}\\approx 1.40\\times 10^{-8}$, and $\\varepsilon^{(3)}\\approx 2.18\\times 10^{-7}$; the six functions $\\pm u^{(i)}$ are pairwise distinct. One solution is dominated by the lowest mode and lies on the known trunk family, whereas the other two carry substantial higher-mode content and are the first rigorously constructed members of the new branch class conjectured from numerical Galerkin computations. A symmetry argument then yields time-periodic solutions of the focusing equation $u_{tt}-u_{xx}-u^3=0$ with frequency $40/96$.","pith_inferences":["Editorial extension: if the branch pairs persist at other rational frequencies of the form (1.3), the solution set would contain infinitely many such pairs accumulating along the trunk, matching the fractal-like picture suggested by numerics.","Editorial extension: the same exact-rational-arithmetic verification scheme should transfer to other 1D semilinear wave equations with polynomial nonlinearities, since it only needs product-to-sum identities and an explicit $\\phi(m,n)$ bound for the linear resolvent.","Editorial extension: a continuation in $\\Omega$ from $69/40$ could test whether the two branch solutions remain close to their numerical approximations, and whether new branch pairs appear at nearby rational frequencies.","Editorial extension: the tiny radii $\\varepsilon^{(i)}$ suggest the approximate Galerkin solutions are very accurate; a reader could rerun the supplied scripts at higher truncation to look for additional branch solutions at the same frequency."],"forward_implications":["There exist at least three pairwise distinct $2\\pi$-periodic solutions of the rescaled equation (1.2) at frequency $\\Omega=69/40$, and hence of the original equation (1.1) after undoing $\\tau=\\Omega t$.","Two of these solutions are the first rigorously confirmed members of the branch family that had previously been seen only numerically.","The same symmetry gives time-periodic solutions of the focusing cubic wave equation with frequency $40/96$.","The rational-frequency condition $\\Omega=(2p+1)/(2q)$ makes the denominators in $L_\\Omega^{-1}$ odd, so the same finite rational verification scheme is not blocked by small divisors at any such frequency.","The explicit bounds $\\varepsilon^{(i)}$ locate actual solutions inside tiny balls around rational approximations, giving quantitative control of the construction."],"supporting_citations":[{"why":"Supplies the numerical Galerkin approximations $u_0^{(i)}$ and the trunk/branch conjecture that this paper makes rigorous.","marker":"[FM25]"},{"why":"Provides the frequency diagram and further numerical evidence of the branch structure used to select $\\Omega=69/40$.","marker":"[FM24]"},{"why":"Contributes the computer-assisted fixed-point and operator-norm-bound strategy that the present proof adapts to rational arithmetic.","marker":"[AK17]"},{"why":"Establishes the classical trunk family of time-periodic solutions, the baseline that one of the three solutions belongs to.","marker":"[LS88]"},{"why":"Establishes Cantor families of periodic solutions for completely resonant nonlinear wave equations, the classical existence theory that the new branch solutions extend and go beyond.","marker":"[BB06]"},{"why":"Provides the software environment in which the supplied exact-rational-arithmetic verification scripts are written and run; the verification of (2.6) depends on these scripts.","marker":"[MMA]"}],"fun_headline_variants":["New periodic wave solutions are rigorously proven to exist","A new branch of wave solutions is confirmed beyond the trunk","Time-periodic solutions are rigorously constructed for a new class","Rigorous proof establishes existence of new periodic wave solutions","Existence of a new class of time-periodic solutions is proven"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof rests on the correctness of the supplied computer scripts, which check the inequalities (2.6) by exact rational arithmetic; if a script or the formula (B.6)/(B.7) for the norm of $H_0$ contains an implementation error, the bounds and hence the existence conclusions could fail.","fun_headline_variants_meta":{"raw":{"variants":["New periodic wave solutions are rigorously proven to exist","A new branch of wave solutions is confirmed beyond the trunk","Time-periodic solutions are rigorously constructed for a new class","Rigorous proof establishes existence of new periodic wave solutions","Existence of a new class of time-periodic solutions is proven"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00054,"raw_usage":{"total_tokens":2525,"prompt_tokens":820,"completion_tokens":1705,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":436,"completion_tokens_details":{"reasoning_tokens":1624}},"tokens_in":436,"tokens_out":1705,"duration_ms":15431,"temperature":1.0,"reasoning_tokens":1624,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:17:51.932104+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Rerun the supplied scripts with the stated data (or independently recompute the maxima in formula (5.2) using interval arithmetic) and check whether, for each $i$, the quantity $\\|H_0\\|+6\\|L_\\Omega^{-1}\\|\\|u_0^{(i)}\\|\\|A\\|^2\\delta+3\\|L_\\Omega^{-1}\\|\\|A\\|^3\\delta^2$ is smaller than $K_0$ and $\\|N_\\Omega(0)\\|$ is smaller than $(1-K_0)\\delta$. If any of these inequalities fails, the contraction argument does not close and the claimed solutions are not established.","supporting_citations":[],"review_version":1}