{"id":"3f017819-a487-4d8e-b81a-189fc73279a2","arxiv_id":"2411.08366","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In four spatial dimensions, small codimension-1 perturbations of catenoid initial data yield global HVMC solutions that converge modulo translation and boost to a boosted/translated catenoid with explicit decay rates.","lead":"This paper proves that the four-dimensional catenoid, a curved solution of the hyperbolic vanishing mean curvature equation, is stable against a class of small perturbations: solutions exist for all time and settle, up to translation and boost, on a nearby catenoid. It completes the stability picture in dimensions n≥3 and introduces a new weighted commutator method for wave decay in even dimensions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The decisive r^p-commutator estimate for Y in §6.5 is not fully proved: the displayed text asserts the absorption of the ∂τ^2Y and ℘r^{-1/2}-type error terms by integration by parts without verifying the resulting endpoint/flux integrals.","rationale":"I read the paper in good faith. The main theorem is a serious nonlinear PDE result, and the visible parts of the argument are coherent: the decomposition into a modified profile p and q, the orthogonality conditions, the modulation equations, the shooting argument, and the r^p-weighted hierarchy are all set up carefully, and the algebraic identity (6.50) checks out. The paper also contains substantial independent material, including the spectral analysis of the catenoid and the weighted elliptic estimates. My concern is not with the overall strategy but with the precise location of the proof's most fragile step. The reader already identified Section 6.5-6.6 as the weak point, and I agree that the new commutator hierarchy is where the proof is least secure. I sharpen the concern: Proposition 6.22 is the theorem that must deliver the improved decay for Φ, and its proof, as reproduced in full text, stops at the point where the error terms R12 and (6.55) are declared acceptable. The missing estimates are not cosmetic; they concern boundary/flux integrals that are not in the stated energy spaces. Because the bootstrap for ˙℘ and the final asymptotic decay both depend on the τ^{-9/4+κ} decay of Φ, this is load-bearing. However, I do not see an internal contradiction or an obviously false estimate; the gap could well be fillable by a longer calculation. Therefore the appropriate verdict remains CONDITIONAL, unchanged from the reader's verdict, with the condition being an expert check of the displayed Step 1 of the proof of Proposition 6.22 and the analogous higher-order and pointwise-decay steps in Section 6.6.","tokens_in":81315,"tokens_out":11737,"duration_ms":116985,"concrete_test":"Complete the proof of Proposition 6.22 by explicitly writing out and bounding the endpoint and bulk terms generated by the two integrations by parts: (1) integrate R12 by parts in τ and bound ∫_{Σ_{τ_i}} χr^{p-2}∂τY∂rY and ∫_{D_{τ_1}^{τ_2}} χr^{p-3}|∂τY|^2 using only the stated bootstrap assumptions and the definitions of E_p, B_p and LE; (2) integrate the O(℘r^{-1/2}+r^{-3/2})∂r∂θY term by parts in θ and verify that the resulting ∂θ-coefficient term is absorbed at p = 3/2-2α for sufficiently large R and small ε. If either step requires a weighted ∂τY norm not implied by (ILED-∂τΦ) or (BA-ϕ-LE), Proposition 6.22 as stated is unsupported and the bootstrap for Φ collapses.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on Proposition 6.22, the new r^p-weighted hierarchy for Y = r^{3/2}∂r Ũ. The commutator identity (6.50) is algebraically correct, and the improved p-range heuristic is plausible, but the proof of Proposition 6.22 is not actually completed in the text. After multiplying (6.53) by χr^p∂rY, the terms R12 (the O(r^{-2})∂τ^2Y term) and the first-order terms in (6.55), especially O(℘r^{-1/2}+r^{-3/2})∂r∂θY and O(℘r^{-1/2}+r^{-3/2})∂rY, are handled by saying 'integration by parts and Cauchy-Schwarz' without displaying the resulting integrals. For R12, integration by parts in τ produces endpoint terms such as ∫_{Σ_{τ_i}} χr^{p-2}∂τY∂rY and a bulk contribution of ∫χr^{p-3}|∂τY|^2-type. These are not manifestly controlled by the E_p[Y]/B_p[Y] hierarchy, because Y involves r^{3/2}∂rŨ and B_p[Y] as stated controls only r^{p-1}(∂rY)^2, r^{p-3}|∇Y|^2 and r^{p-3}|Y|^2, not ∂τY with weight r^{p-3}. The stated ILED norms control ρ^{-1-α}|∂U|^2 for U, which is not the same as a weighted ∂τY estimate. A similar but more elementary issue arises for (6.55): after integrating by parts in θ, the cross terms must be absorbed using the smallness of ℘ and the largeness of R, and the uniformity of this absorption as p approaches 3/2-2α is not shown. If any of these flux terms is not absorbable, Proposition 6.22 fails, and with it the improved τ^{-9/4+κ} decay for Φ in §6.6, the proof of Proposition 5.4, and ultimately the bootstrap for the parameter derivatives in Theorem 5.7. This is the single most load-bearing step in the paper.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proves asymptotic stability of the Minkowski-space catenoid as a solution of the hyperbolic vanishing mean curvature equation in n = 4 spatial dimensions, for compactly supported codimension-one perturbations and modulo translation and boost modulation. The proof follows the strategy of Lührmann--Oh--Shahshahani for n ≥ 5 and Oh--Shahshahani for n = 3, but replaces the strong Huygens principle used in n = 3 with a new commutator vector field K = r^{3/2}∂r. The announced decay rates are |ℓ̇|, |ξ̇-ℓ| ≲ ε⟨τ⟩^{-9/4+κ} and ∥ψ∥_{L∞(Στ)} ≲ ε⟨τ⟩^{-2+κ}. The main technical novelty is a family of r^p-weighted estimates for Y = r^{3/2}∂r Ũ, culminating in Proposition 6.22, from which an improved late-time tail is extracted in Section 6.6.","tokens_in":81845,"tokens_out":4700,"duration_ms":51231,"significance":"If the central estimate is completed, the result closes the n = 4 case of the catenoid stability problem and, more importantly, introduces a genuinely new commutator-based r^p hierarchy that appears applicable to other even-dimensional quasilinear wave equations and to wave equations with inverse-square potentials. The paper is also valuable for its detailed derivations: Lemma 2.10 provides a clean cancellation giving F0 = O(℘̇ r^{-3}), Section 3.2 gives an explicit spherical-harmonic proof of the zero-mode classification of the catenoid stability operator, and Sections 4.2--4.5 lay out the gauge choice, modified profile, modulation equations, and unstable-mode decomposition carefully. The bootstrap architecture is coherent, with Remark 5.1 correctly identifying and breaking the linear circularity between parameter derivatives and energy via the δ℘ smallness. However, the decisive new estimate is not fully proved in the submitted text: the proof of Proposition 6.22 stops after forming the preliminary identity (6.57), and the absorption of the displayed error terms is asserted rather than shown.","major_comments":[{"comment":"The proof of Proposition 6.22 is not completed: after deriving (6.57) the text announces 'Step 1: Obtain (6.62)' and then breaks off, so the decisive absorption estimates for R12 and for the G2 terms in (6.54)--(6.55) are not displayed. In particular, R12 = ∫ χ r^p ˜m^{ττ}_1 ∂τ^2 Y ∂r Y, after integration by parts in τ, produces endpoint terms of the form ∫_{Σ_{τ_i}} χ r^{p-2} ∂τ Y ∂r Y and bulk terms of type ∫ χ r^{p-3} |∂τ Y|^2. Neither is manifestly controlled by the quantities appearing on the right-hand side of (6.56): B_p[Y] controls r^{p-1}(∂rY)^2, r^{p-3}|∇Y|^2 and r^{p-3}|Y|^2, while the stated ILED norms control ρ^{-1-α}|∂U|^2 for U, not a weighted ∂τY norm. Similarly, the two terms O(℘r^{-1/2}+r^{-3/2})∂r∂θY and O(℘r^{-1/2}+r^{-3/2})∂rY in (6.55) require an integration-by-parts/absorption argument whose uniformity as p approaches 3/2-2α is not shown. Since Proposition 6.22 is the mechanism behind the improved τ^{-9/4+κ} decay in §6.6 and hence behind the bootstrap for Φ in Proposition 5.4 and Theorem 5.7, this is a load-bearing gap that must be filled before the main theorem can be considered established.","section":"§6.5, Proposition 6.22"},{"comment":"The claims that the Y-hierarchy implies pointwise decay for Φ are not verifiable from the submitted text. The passage from weighted-energy bounds on dyadic sequences to the pointwise bounds ∥Φ∥_{L∞(Στ)} ≲ ε⟨τ⟩^{-2+κ} is only summarized, and the fragile endpoint effects as the power p approaches the coercivity threshold 3/2-2α are not addressed. In particular, the final step in Section 6.6 is not displayed. This is not a cosmetic omission because the improved decay for Φ is exactly what feeds back into the parameter-derivative estimates and closes the bootstrap.","section":"§6.6 and proof of Proposition 5.4"}],"minor_comments":[{"comment":"The statement of Proposition 6.22 uses the notation B_p[Y] without defining it in that section; the definition appears only heuristically in Section 1.1.2. A formal definition should be included before the proposition.","section":"§6.5, Proposition 6.22"},{"comment":"The remark stating 'τ is not the same as τ' appears to contain a typo and should be rewritten to distinguish the global coordinate time from the hyperboloidal time coordinate.","section":"§2.6, Remark 2.13"},{"comment":"The notation η′ is used in the metric formulas near (6.1) without an explicit definition at that point; the reader must recall η(τ) = ξ(τ) - γ(τ)R_f ℓ(τ) from Section 2.3. A brief reminder would improve readability.","section":"§6.1.1"},{"comment":"The rough version writes R^{4+1} and R^{1+(4+1)} inconsistently; the notation should be unified with the rest of the paper.","section":"Theorem 1.1"}],"recommendation":"major_revision","confidential_remarks":"The incomplete proof in Section 6.5 is the only substantive obstacle I see to accepting the central claim. The rest of the architecture is careful and coherent, and the gap is local in character: it concerns the explicit absorption of specific flux and error terms in Proposition 6.22. I therefore recommend major revision rather than rejection, giving the author the opportunity to supply the missing estimates."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper delivers the n=4 catenoid stability theorem, which is genuinely the last missing no-symmetry case after n≥5 and n=3. The new tool, the commutator K = r^{3/2}∂r, is a real idea: it extends the r^p hierarchy from p<2 to p<5/2 at the level of U, and it is specifically adapted to even dimensions where the Huygens-based trick of OS24 doesn't work. The visible parts of the paper are careful and honest. The spectral analysis of the zero modes is clean, the modulation setup is handled in detail, and Remark 5.1 explicitly addresses the circularity between parameter derivatives and energy. The derivation of the source term F0 = O(℘̇r^{-3}) is a serious computation, and the paper does not oversell its debt to LOS22 and OS24.\n\nBut the stress-test concern is real and load-bearing. Proposition 6.22, which produces the improved τ^{-9/4+κ} decay for Y, is not proved in the text. After multiplying by χr^p∂rY, the terms R12 and the first-order terms in (6.55) are dismissed with 'integration by parts and Cauchy-Schwarz' without displaying the resulting flux and bulk integrals. For R12, integration by parts in τ gives endpoints like ∫χr^{p-2}∂τY∂rY and a bulk ∫χr^{p-3}|∂τY|^2, neither of which is manifestly controlled by the stated B_p[Y] hierarchy, which only controls r^{p-1}(∂rY)^2, r^{p-3}|∇Y|^2, and r^{p-3}|Y|^2. The same issue appears in the θ-integration of the ℘r^{-1/2}∂r-type terms near the endpoint p→3/2-2α. If these fluxes are not absorbable, the bootstrap for Φ and the parameter derivatives collapses. This is exactly the kind of step a referee must see written out completely.\n\nMinor points: the 'full picture' claim is overstated, and the paper leans heavily on black-box ILED from prior work, but that is not by itself a flaw when the cited results are established.\n\nWho is this for? Specialists in nonlinear wave equations and stability of minimal submanifolds. It deserves a serious referee, but the referee should be instructed to verify Section 6.5 in full before acceptance. I would engage with it myself only after the gap is closed.","headline":"Serious candidate for the final catenoid stability dimension, but the decisive commutator estimate in Section 6.5 is not actually proved, and the paper should be refereed under a requirement to fill that gap.","tokens_in":82418,"tokens_out":1469,"would_cite":false,"duration_ms":18173,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35L72","35B40","53C42","35Q75"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves that sufficiently small codimension-1 perturbations of the catenoid initial data in four spatial dimensions evolve globally and converge, modulo translation and boost, to a boosted and translated catenoid.","keywords":["hyperbolic vanishing mean curvature","catenoid stability","late-time tails","commutator vector field","r^p-weighted estimates","modulation","even-dimensional wave equations","codimension-1 perturbations"],"falsifier":"Directly verify the absorption step in the new hierarchy on the model linear problem in 1+4 Minkowski space: take a solution of Q_0\\tilde U = F with F = O(\\dot\\wp $r^{{-3}}$) and run the Y = $r^{{3/2}}$\\partial_r\\tilde U estimates of Section 6.5 for p arbitrarily close to 3/2. The claim stands only if the integrals generated by the error terms O(\\wp $r^{{-1/2}}$+$r^{{-3/2}}$)\\partial_r in (6.51)-(6.55) can be bounded by a small constant times the coercive energy plus the source term, and if the boundary term \\lim_{r\\to\\infty}\\int_{$S^{3}$} $r^{{p-1}}$$Y^{2}$ d\\$\\theta$ vanishes for p<2; a computed counterexample to either point would falsify the $τ^{{-9/4+\\kappa}}$ decay claim and collapse the bootstrap.","tokens_in":81089,"feed_emoji":"🌊","tokens_out":13866,"duration_ms":133954,"temperature":0.7,"pith_summary":"Sufficiently small, compactly supported, codimension-1 perturbations of the catenoid initial data in four spatial dimensions produce global solutions of the hyperbolic vanishing mean curvature (HVMC) equation, and these solutions converge, modulo a time-dependent translation and boost, to a boosted and translated catenoid. This is the n=4 case of a stability problem already settled for n≥5 and n=3; the case is harder because the catenoid decays only polynomially in space, waves in even dimensions decay slowly in time, and the strong Huygens principle available in three dimensions is absent. The proof achieves the necessary improved late-time tails by introducing a commutator vector field K = $r^{{3/2}}$\\partial_r that extends the standard r^p-weighted decay hierarchy beyond its usual range. The quantitative conclusion is that |\\dot\\ell|, |\\dot\\xi-\\ell| \\lesssim \\epsilon\\langle\\tau\\$rangle^{{-9/4+\\kappa}}$ and \\|\\psi\\|_{L^\\infty(\\Sigma_\\tau)}\\lesssim\\epsilon\\langle\\tau\\$rangle^{{-2+\\kappa}}$.","feed_headline":"Perturbed 4-D catenoids evolve to boosted, shifted catenoids","feed_subtitle":"A new commutator vector field extends the decay hierarchy and closes the stability proof without symmetry assumptions.","key_machinery":"The central object is the commutator vector field K = $r^{{3/2}}$\\partial_r applied to the rescaled variable \\tilde U = $r^{{3/2}}$U in the outgoing hyperboloidal coordinates (τ,r,θ), with Y=K\\tilde U. Acting on the leading operator Q_0 = -\\frac34 $r^{{-2}}$ -2\\partial_r\\partial_\\tau+\\$partial_r^{2}$+$r^{{-2}}$\\Delta_{$S^{3}$}, the commutator identity (K+$2r^{{1/2}}$)Q_0=Q_1K produces a new operator Q_1 from which the inverse-square potential -\\frac34 $r^{{-2}}$\\tilde U has disappeared. The r^p-weighted energies for Y are coercive for p<3/2, and because Y carries an extra weight $r^{{3/2}}$, the effective range for \\tilde U extends from p<2 to p<5/2; this is what yields the improved $τ^{{-9/4+\\kappa}}$ decay. Around this sits the modulation machinery—modified profile, orthogonality conditions, time-smoothing operator, and the shooting and trapping argument—that converts the decay into convergence of the translation and boost parameters.","core_discovery":"On its own terms, the paper claims Theorem 5.7: given compactly supported initial data that are ε-close to the catenoid in suitable weighted Sobolev norms and satisfy one scalar 'codimension-1' constraint, there exists a global embedding Φ solving the HVMC equation, foliated by hyperboloidal leaves Στ, such that the perturbation ψ satisfies the bounds above and the modulation parameters converge. No symmetry assumption is needed. The proof is a bootstrap built on a decomposition ψ=p+q, with a modified profile p absorbing the slowly decaying source term O(\\dot\\wp $r^{{-3}}$), followed by a decomposition q=ϕ+a_+Z_++a_-Z_- separating the exponentially growing mode. Orthogonality conditions keep the perturbation transverse to the translation and boost zero modes and to the unstable neck mode; a shooting argument selects the codimension-1 initial slice so that the unstable coefficient a_+ stays trapped under its decay barrier. Since |\\dot\\ell|+|\\dot\\xi-\\ell| is integrable in τ, the translated and boosted catenoid parameters have limits as τ→∞.","pith_inferences":["Editorial inference: the same commutator mechanism should yield a direct proof of improved late-time tails for the scalar wave equation with inverse-square potential in 1+4 dimensions; the sharp boundary of the p-range for Y could be tested on the model linear equation before committing to the full quasilinear bootstrap.","Editorial inference: the codimension-1 shooting constraint is likely tied to the sign of the initial projection onto the exponentially growing mode; a simpler threshold condition on a_+(0) might replace the abstract shooting choice in applications.","Editorial inference: because the n=3 proof relied on Huygens cancellation and the n=4 proof replaces it with a commutator, the even-dimensional cases n=6,8,... may be approachable by iterating the same commutator or powers r^{k/2}\\partial_r, with the working range of p shrinking each time.","Editorial inference: one could test numerically whether the decay rate τ^{-9/4+\\kappa} for ∂τΦ is sharp by solving the linearized HVMC equation around the catenoid with a source O(r^{-3}) that is compactly supported in time; if the tail is faster, the bootstrap rate is not optimal."],"forward_implications":["The catenoid is asymptotically stable in the borderline even dimension n=4: the two-parameter family of boosted and translated catenoids acts as the attractor for a codimension-1 open set of small perturbations, without any symmetry assumption.","The modulation parameters ξ(τ) and ℓ(τ) converge as τ→∞, because the pointwise decay |\\dot\\ell|, |\\dot\\xi-\\ell| \\lesssim \\epsilon\\langle\\tau\\rangle^{-9/4+\\kappa} is integrable; the limiting solution is exactly a boosted and translated catenoid.","The improved decay mechanism works for even-dimensional wave equations: the commutator K=r^{3/2}\\partial_r extends the standard r^p-weighted hierarchy, heuristically to p<5/2, and is expected to apply to other quasilinear wave equations in even dimensions or with inverse-square potentials.","The spectral analysis justifies that the only instabilities of the linearized catenoid are translations, boosts, and the shrinking of the neck, so the codimension-1 shooting condition is a genuine and sharp restriction.","The proof establishes full bootstrap closure at the level of energy and pointwise norms, including the L^2-type estimates for second derivatives of the parameters needed for the continuation argument."],"supporting_citations":[{"why":"Supplies the overall stability framework for n≥5, including the foliations, first-order formulation, integrated local energy decay estimates, and bootstrap structure that this paper extends to n=4.","marker":"[LOS22]"},{"why":"Supplies the n=3 stability proof and the modified-profile construction that absorbs the slowly decaying source term in the n=4 argument.","marker":"[OS24]"},{"why":"Provides the commutator-weight idea; the paper adapts it with K=r^{3/2}\\partial_r because the r^2\\partial_r version used there loses coercivity on non-spherically-symmetric backgrounds.","marker":"[AAG18]"},{"why":"Provides the r^p-weighted energy hierarchy that the new commutator extends; the improved p-range is the paper's main quantitative gain.","marker":"[DR10]"},{"why":"Establishes that the catenoid's stability operator has Morse index one, namely one positive eigenvalue, which is the spectral input for the unstable-mode decomposition.","marker":"[TZ09]"}],"fun_headline_variants":["4D catenoid stability proven via new commutator field","Catenoid stable in 4D under codim-1 initial data","Commutator vector field yields catenoid stability in 4D","4D HVMC catenoid stable: new decay hierarchy","Catenoid stability in 4D via modulated translation boost"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the new commutator-based hierarchy closes: all error terms in the weighted estimates of Section 6.5, particularly the borderline terms of size (small parameter times $r^{{-1/2}}$ plus $r^{{-3/2}}$) acting on derivatives, can be absorbed at the claimed decay rate $τ^{{-9/4+\\kappa}}$. If that absorption fails, the pointwise decay for the perturbation and the convergence of the modulation parameters collapse.","fun_headline_variants_meta":{"raw":{"variants":["4D catenoid stability proven via new commutator field","Catenoid stable in 4D under codim-1 initial data","Commutator vector field yields catenoid stability in 4D","4D HVMC catenoid stable: new decay hierarchy","Catenoid stability in 4D via modulated translation boost"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000829,"raw_usage":{"total_tokens":3667,"prompt_tokens":1037,"completion_tokens":2630,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":653,"completion_tokens_details":{"reasoning_tokens":2542}},"tokens_in":653,"tokens_out":2630,"duration_ms":17441,"temperature":1.0,"reasoning_tokens":2542,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T21:39:40.429589+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Directly verify the absorption step in the new hierarchy on the model linear problem in 1+4 Minkowski space: take a solution of Q_0\\tilde U = F with F = O(\\dot\\wp $r^{{-3}}$) and run the Y = $r^{{3/2}}$\\partial_r\\tilde U estimates of Section 6.5 for p arbitrarily close to 3/2. The claim stands only if the integrals generated by the error terms O(\\wp $r^{{-1/2}}$+$r^{{-3/2}}$)\\partial_r in (6.51)-(6.55) can be bounded by a small constant times the coercive energy plus the source term, and if the boundary term \\lim_{r\\to\\infty}\\int_{$S^{3}$} $r^{{p-1}}$$Y^{2}$ d\\$\\theta$ vanishes for p<2; a computed counterexample to either point would falsify the $τ^{{-9/4+\\kappa}}$ decay claim and collapse the bootstrap.","supporting_citations":[],"review_version":1}