{"id":"9b01aa9c-b3ed-40ec-a978-a8b0906039f3","arxiv_id":"2608.18826","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For the damped linear wave equation on a path graph with Dirichlet boundary conditions, the paper provides explicit solution formulas and a complete classification of self-similar solutions.","lead":"This paper studies the semi-discrete hyperbolic curvature flow, a damped wave equation on a chain of curve vertices with prescribed endpoints. It writes down explicit formulas for all solutions and classifies self-similar, periodic, and curve-to-curve evolutions of this linear system.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 2's explicit diagonal solution formula is invalid at critical damping values where M has a Jordan block, so the advertised 'general solution for all beta' is only conditional.","rationale":"The reader's weakest_assumption correctly identifies the point where the central claim is least secure. The paper's advertised contribution is explicit finite-Fourier formulas for general boundary data, and every one of those formulas—Theorem 2, Corollary 3, Theorem 3, Corollary 4, and Proposition 11—passes through P diag(...) P^{-1}. At beta = 4 sin(k pi / (2(m+1))) the quadratic eigenvalue problem has a repeated root and Proposition 2 provides only one eigenvector for that mode; hence P is singular and the displayed formulas are undefined. This is not an external disagreement with consensus but an internal inconsistency in the proof of Theorem 2. The proposed check is a small explicit computation that exhibits the singular P and verifies the Jordan-form correction. Because the underlying system is linear with constant coefficients, existence and uniqueness survive, the diagonal form can be repaired locally, and the self-similar classification and difference-flow results do not appear to collapse. The reader's CONDITIONAL verdict with high confidence is therefore appropriate, and my pass does not move it.","tokens_in":23920,"tokens_out":20972,"duration_ms":188137,"concrete_test":"Set n = 4 (m = 2) and beta = 2, so mode k = 1 is critical. Form A = [[-2, 1], [1, -2]] and M = [[0, I], [A, -2I]]. Build the eigenvector matrix P by the recipe in Proposition 2 and verify that the two k = 1 eigenvectors coincide, so det(P) = 0. Then compare the literal expression (16) with direct numerical integration of dW/dt = M W + F for one initial condition and constant forcing; the formula cannot be evaluated because P^{-1} does not exist, and any pseudo-inverse substitute will not reproduce the direct solution. Finally, replace the critical 2x2 block by the Jordan exponential e^{-t} [[1, t], [0, 1]] and confirm that the corrected expression matches direct integration to round-off. This check separates the fixable formula gap from any failure of existence or uniqueness.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central explicit-solution theorem, Theorem 2, states that the first-order matrix M = [[0, I], [A, -beta I]] is diagonalisable by the eigenvector matrix P from Proposition 2, and writes the solution as P diag(e^{lambda_j t}) P^{-1}, equation (16). This fails for every beta >= 0 such that beta^2 + 4 mu_k = 0 for some mode k, where mu_k = -4 sin^2(k pi / (2(m+1))) are the eigenvalues of A. At such a critical value the two quadratic roots coincide at lambda = -beta/2, and Proposition 2 supplies two identical eigenvectors (u_k, -beta/2 u_k)^T for that mode. The geometric multiplicity is one, so P is singular and the displayed diagonal formula is undefined. The exceptional set is nonempty: for n = 4, m = 2, mode k = 1 gives beta = 2. The defect propagates to Corollary 3, Theorem 3, Corollary 4, and Proposition 11 whenever they invoke P diag(...) P^{-1}. Existence and uniqueness of the linear ODE do survive through the matrix exponential, and a Jordan-form correction repairs the explicit formula, but as stated the theorem overclaims that formula (16) is valid for all beta >= 0.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the semi-discrete hyperbolic curvature flow (SDHF) of piecewise linear curves with prescribed boundary trajectories, defined by U¨ + βU˙ = A U + f(t). The authors derive an energy dissipation law, compute the spectra of the spatial matrix A and the first-order block matrix M, classify all self-similar solutions into four classes (pure scaling-rotation, pure translation, scaling-rotation with translation, and static equilibrium), construct time-periodic solutions, give an explicit variation-of-parameters solution formula for general initial data and arbitrary boundary forcing, introduce a hyperbolic curvature difference flow that evolves one polygonal curve to another, and present anisotropic and matrix-coefficient variants. The main theoretical tools are the quadratic eigenvalue problem, Sylvester's theorem, and the matrix exponential.","tokens_in":24150,"tokens_out":28903,"duration_ms":216677,"significance":"If the claims are correct, the paper provides a complete explicit theory of the linear semi-discrete hyperbolic flow with boundary: a closed-form solution for all times and boundary data, a full classification of self-similar solutions (including translations, which are impossible in the closed-polygon setting), and a Yau-type flow between arbitrary polygonal curves. The reduction of the self-similar classification to a Sylvester equation and the treatment of the damped quadratic eigenvalue problem are elegant and likely useful for the wider discrete geometric flow community. The paper is largely self-contained and its spectral formulas are explicit and directly verifiable. However, the main representation theorem is invalid at critical damping values where the first-order matrix is defective, and one anisotropic solution formula is incorrect; these issues are local and repairable but currently limit the advertised generality.","major_comments":[{"comment":"The proof of Theorem 2 asserts that the matrix M = [[0,I],[A,-βI]] is diagonalisable for every β ≥ 0, with eigenvector matrix P from Proposition 2. This fails at critical damping: when β² + 4µ_k = 0 for some mode k, the two quadratic roots coincide at λ = -β/2 and Proposition 2 supplies two identical eigenvectors (u_k, -β/2 u_k)^⊤ for that mode, so P is singular and the displayed diagonal formula (16) is undefined. The exceptional values are nonempty (e.g., n=4, m=2, k=1 gives β=2). The defect propagates to Corollary 3, Theorem 3, Corollary 4, and Proposition 11, all of which invoke P diag(...)P^{-1}. Existence and uniqueness of the linear ODE do survive through the matrix exponential, and a Jordan-form correction or an explicit restriction on β would repair the statement; as written, the theorem overclaims validity for all β ≥ 0.","section":"§3.3, Theorem 2 (with Proposition 2)"},{"comment":"The solution formula for the coordinate-wise anisotropic flow (23) uses the eigenvector matrix P of the block matrix for the isotropic operator A and the same eigenvalues λ_i, rescaled by c_j: P diag(e^{c_j λ_i t})P^{-1}. But for the coordinate equation with coefficient matrix c_j A, the relevant eigenvalues are the roots of λ² + βλ - c_j µ_k = 0, not c_j λ_i, and the corresponding eigenvector matrix P_j has columns (u_k, λ u_k) with these new λ, so P_j ≠ P unless c_j = 1. A simple scalar check (m=1, β=0, c_j=2: the true eigenvalues are ±2, while c_j λ_i gives ±2√2) confirms the formula is incorrect. The proof (\"follows in the same manner as Theorem 2\") does not justify the rescaling, and the displayed solution should be corrected or the proposition removed.","section":"§5.1, Proposition 9"}],"minor_comments":[{"comment":"The discrete Dirichlet energy W sums i=2 to n-1, omitting the edge (X_{n-1},X_n). With this definition the summation-by-parts identity in the proof leaves the boundary term <X_n - X_{n-1}, dX_{n-1}/dt>, so the claimed identity (5) and the dissipation law (4) do not follow. The correct statement holds if the sum is taken to i=n; please correct the index in (3).","section":"§2.1, Eq. (3)"},{"comment":"The four classes are not mutually exclusive as stated, since the static equilibrium (g≡1, R≡I, h≡0) is a special case of the pure scaling-rotation class (i) with γ=0, ω0=0 and satisfies the Sylvester equation (13) with M0=0. Recommend phrasing (i) as the non-trivial scaling-rotation case or explicitly noting the overlap.","section":"§3.1, Theorem 1"},{"comment":"In the displayed evaluation of the integral, the terms (e^{λ_i} - 1)/λ_i are missing the time variable; they should read (e^{λ_i t} - 1)/λ_i.","section":"§3.3, Corollary 3 proof"},{"comment":"The condition \"fn = Xn(0)\" should be \"fn(0) = Xn(0)\" to be consistent with the boundary condition notation.","section":"§4, Theorem 3 statement"},{"comment":"The index \"j=1,2,...,2m\" in the eigenvalue listing is misleading; for each k there are exactly two eigenvalues, so the set should be written as λ_k^±.","section":"§2.2, Proposition 2"}],"recommendation":"major_revision","confidential_remarks":"The errors identified are substantive but local: the diagonalisability overclaim can be fixed by a Jordan-form treatment of the critical damping cases, and the coordinate-wise anisotropy formula needs the correct eigenvalues of c_j A. With those corrections the paper would be a solid contribution to the semi-discrete geometric flows literature. I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Core issue: the paper's headline explicit solution formula (16) is not valid for all β ≥ 0. At values where β² = -4μ_k for some mode, the first-order matrix M becomes defective; the two eigenvalues coincide at -β/2 but there is only one eigenvector for that mode, so the eigenvector matrix P is singular and P diag(...) P^{-1} is undefined. This is a genuine overclaim, though a repairable one: existence and uniqueness still follow from the matrix exponential, and a Jordan block or a limiting argument fixes the formula. The exceptional β's are isolated, so the advertised 'general solution' is correct on a dense open set.\n\nWhat is actually new: the self-similar classification (Theorem 1) is a real contribution. The reduction to a Sylvester equation via the Riccati/ODE lemma is clean and gives the four mutually exclusive classes. The energy dissipation law, the undamped standing waves, the forced periodic stability, and the hyperbolic curvature difference flow are all argued well and appear correct. The numerical experiments are modest but support the qualitative claims. The paper is self-contained and cites the prior parabolic boundary work appropriately; the self-citations are to the direct predecessor, not padding.\n\nThe main soft spot beyond the defective case is Section 5.1. In the coordinate-wise anisotropic flow, the eigenvalues of the first-order matrix do not simply scale as c_j times the original eigenvalues when β > 0; the displayed e^{c_j λ t} formula looks off unless the λ's are redefined. This is a minor but genuine issue in a variant section. The vertex-wise and coefficient variants are fine.\n\nBottom line: this paper is a solid, clearly written extension, useful for people who want an explicit solvable semi-discrete hyperbolic flow with boundary. It needs a moderate revision to patch the defective damping case and clean up Proposition 9. I would send it to a serious referee, but I would not accept the current version as-is.","headline":"Competent extension with explicit solution formulas and a correct self-similar classification, but the main diagonal formula fails at critical damping where the system matrix becomes defective.","tokens_in":24697,"tokens_out":6743,"would_cite":false,"duration_ms":57473,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":["34A26","37C60"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that the semi-discrete hyperbolic curvature flow with prescribed endpoints is completely solvable: every self-similar motion falls into one of four classes, and for arbitrary boundary trajectories the unique solution is…","keywords":["semi-discrete curvature flow","hyperbolic geometric flow","self-similar solutions","quadratic eigenvalue problem","piecewise linear curves","damped wave equation","curve shortening flow","curvature difference flow"],"falsifier":"Choose $n = 5$ and set $\\beta = 2\\sqrt{2 - 2\\cos(\\pi/4)} \\approx 1.531$, so that the discriminant $\\beta^2 + 4\\mu_1$ vanishes for the first mode. Compute the null space of $M - (-\\beta/2)I$; if it has dimension one while the eigenvalue has algebraic multiplicity two, the matrix is not diagonalisable and formula (16) is undefined at a value of $\\beta$ the paper claims to cover.","tokens_in":23675,"feed_emoji":"📐","tokens_out":5638,"duration_ms":54136,"temperature":0.7,"pith_summary":"The paper studies a discrete model of inertia in curvature-driven motion: each interior vertex of a piecewise linear curve accelerates in the direction of its discrete curvature, with damping proportional to velocity, while the two endpoints move along prescribed trajectories. The central claim is that this second-order geometric evolution is completely solvable. All self-similar evolutions — scaling, rotation, translation, their combinations, and the static straight-line interpolant — are classified, and for any initial curve and any endpoint motion the unique solution is written down explicitly using the eigenmodes of the discrete Laplacian. The same machinery produces a curvature-difference flow that evolves any piecewise linear curve to any other with the same number of vertices. A sympathetic reader would care because this makes a nonlinear-looking geometric flow exactly solvable, exposing the same spectral structure that governs damped vibrating chains.","feed_headline":"Open-curve hyperbolic flow solved exactly","feed_subtitle":"All self-similar motions fall into four classes; any curve can evolve to any other.","key_machinery":"The central object is the block matrix $M = \\begin{pmatrix} 0 & I \\\\ A & -\\beta I \\end{pmatrix}$ obtained by reducing the second-order flow to a first-order system; its eigenvalues solve the quadratic eigenvalue problem $\\lambda^2 + \\beta\\lambda - \\mu_k = 0$, where $\\mu_k = -2 + 2\\cos(k\\pi/(n-1))$ are the eigenvalues of the discrete Laplacian $A = \\operatorname{tridiag}(1,-2,1)$ with Dirichlet boundary conditions. This spectral data produces the matrix exponential $e^{Mt}$, whose diagonalised form $P \\operatorname{diag}(e^{\\lambda_j t}) P^{-1}$ carries the variation-of-parameters solution (16) and the closed-form constant-forcing version (18). The self-similar classification instead runs through the Sylvester equation $A\\mathbf{U}_0 - \\mathbf{U}_0 M_0 = -\\alpha$, which decides when a prescribed scaling-rotation admits a unique interior profile. The discrete energy $E = K + W$ and its dissipation law $dE/dt = -2\\beta K$ control long-time behaviour, including convergence to the straight-line interpolant.","core_discovery":"The load-bearing discovery is Theorem 2 combined with Theorem 1. Theorem 2 states that the initial-value problem for the semi-discrete hyperbolic curvature flow with arbitrary boundary trajectories has a unique solution, given by the variation-of-parameters formula (16), in which the first-order system matrix $M = \\begin{pmatrix} 0 & I \\\\ A & -\\beta I \\end{pmatrix}$ is diagonalised through the eigenvalues of the tridiagonal Toeplitz matrix $A = \\operatorname{tridiag}(1,-2,1)$. Theorem 1 states that every self-similar solution, for a non-degenerate initial curve and self-similar boundary data in the plane, falls into exactly one of four mutually exclusive classes: pure scaling-rotation, pure translation, scaling-rotation with translation, or static equilibrium. In the scaling-rotation classes, the interior profile is the unique solution of a Sylvester equation $A\\mathbf{U}_0 - \\mathbf{U}_0 M_0 = -\\alpha$, with existence controlled by disjointness of the spectra of $A$ and $M_0$. The paper also proves a hyperbolic curvature-difference flow that carries any piecewise linear curve to any other with the same vertex count, with exponential convergence for $\\beta > 0$.","pith_inferences":["Beyond the paper's own claims, the critical-damping values where $\\beta^2 + 4\\mu_k = 0$ fall inside the stated domain $\\beta \\ge 0$; a Jordan-form correction, adding factors like $t e^{-\\beta t/2}$, would likely extend the explicit solution formula to those parameter values.","Because the coefficient matrix is the path graph Laplacian, the open-curve results suggest that wave-diffusion processes on graphs with Dirichlet boundaries are explicitly solvable, and that periodic boundary forcing creates a boundary-driven resonance that has no analogue on closed graphs.","The coordinate-wise anisotropic variant decouples by ambient coordinate, so the predicted evolutions could be checked directly with a physical mass-spring chain whose endpoint displacements follow the prescribed trajectories."],"forward_implications":["All self-similar evolutions of the flow fall into exactly four classes — pure scaling-rotation, pure translation, combined scaling-rotation-translation, and static equilibrium — with no other shapes possible.","For any initial curve and any prescribed endpoint trajectories there is a unique solution, given explicitly by a finite spectral sum; with constant boundary data the integral term is evaluated in closed form.","With fixed boundary points, every solution converges to the straight line interpolant between the endpoints, and the discrete energy decreases monotonically at rate $-2\\beta K$.","For $\\beta > 0$ and $T$-periodic boundary forcing, a unique $T$-periodic solution exists and is globally attracting.","The hyperbolic curvature-difference flow evolves any piecewise linear curve to any other with the same number of vertices, converging exponentially when the damping is positive."],"supporting_citations":[{"why":"Chow and Glickenstein's semidiscrete geometric flows of polygons provide the original flow and the definition of discrete curvature used here.","marker":"[1]"},{"why":"McCoy and Meyer's closed-polygon hyperbolic polyharmonic flows supply the closed-curve analogue and the breather and periodic-solution ideas that the boundary case modifies.","marker":"[10]"},{"why":"The authors' earlier parabolic flow with boundary is the direct predecessor, contributing the boundary forcing setup, spectral properties, and the Yau-type curvature difference problem.","marker":"[11]"},{"why":"Tisseur and Meerbergen's quadratic eigenvalue problem survey frames the spectral analysis of the damped second-order system.","marker":"[12]"},{"why":"Smith's numerical PDE text is the source for the tridiagonal Toeplitz eigenvalue and eigenvector formulas used throughout the spectral lemmas.","marker":"[17]"},{"why":"Horn and Johnson's Matrix Analysis provides the Sylvester equation solvability theorem used to prove existence and uniqueness of self-similar profiles.","marker":"[23]"},{"why":"Hartman's ODE text supplies the variation-of-parameters formula on which the general solution Theorem 2 and its corollaries rest.","marker":"[24]"},{"why":"Lin and Tsai's length-preserving linear flow frames the problem of evolving one curve to another, which the curvature-difference flow addresses in the semi-discrete hyperbolic setting.","marker":"[25]"}],"fun_headline_variants":["Exact solution for semi-discrete hyperbolic curvature flow","Four classes of self-similar curves, all solved exactly","Any curve evolves to any other: exact hyperbolic flow","Hyperbolic flow carries any polygon to any other"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole argument leans on the assumption that the system matrix can always be diagonalised; at special damping values two modes merge and only one eigenvector remains, so that assumption fails and the formula's eigenvector matrix does not exist.","fun_headline_variants_meta":{"raw":{"variants":["Exact solution for semi-discrete hyperbolic curvature flow","Four classes of self-similar curves, all solved exactly","Any curve evolves to any other: exact hyperbolic flow","Hyperbolic flow carries any polygon to any other"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000469,"raw_usage":{"total_tokens":2278,"prompt_tokens":829,"completion_tokens":1449,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":445,"completion_tokens_details":{"reasoning_tokens":1386}},"tokens_in":445,"tokens_out":1449,"duration_ms":10459,"temperature":1.0,"reasoning_tokens":1386,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-27T19:51:29.347693+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Choose $n = 5$ and set $\\beta = 2\\sqrt{2 - 2\\cos(\\pi/4)} \\approx 1.531$, so that the discriminant $\\beta^2 + 4\\mu_1$ vanishes for the first mode. Compute the null space of $M - (-\\beta/2)I$; if it has dimension one while the eigenvalue has algebraic multiplicity two, the matrix is not diagonalisable and formula (16) is undefined at a value of $\\beta$ the paper claims to cover.","supporting_citations":[{"cited_title":"The American Mathematical Monthly 114(4), 316–328 (2007)","cited_arxiv_id":null,"evidence_quote":"Chow and Glickenstein's semidiscrete geometric flows of polygons provide the original flow and the definition of discrete curvature used here."},{"cited_title":"Mathematics in Engineering 7(3), 281–315 (2025) https://doi.org/10","cited_arxiv_id":null,"evidence_quote":"McCoy and Meyer's closed-polygon hyperbolic polyharmonic flows supply the closed-curve analogue and the breather and periodic-solution ideas that the boundary case modifies."},{"cited_title":"Geometriae Dedicata 220(3) (2026) https://doi.org/10.1007/ s10711-026-01080-3","cited_arxiv_id":null,"evidence_quote":"The authors' earlier parabolic flow with boundary is the direct predecessor, contributing the boundary forcing setup, spectral properties, and the Yau-type curvature difference problem."},{"cited_title":"SIAM Review 43(2), 235–286 (2001)","cited_arxiv_id":null,"evidence_quote":"Tisseur and Meerbergen's quadratic eigenvalue problem survey frames the spectral analysis of the damped second-order system."},{"cited_title":"Oxford Press, ??? (1979)","cited_arxiv_id":null,"evidence_quote":"Smith's numerical PDE text is the source for the tridiagonal Toeplitz eigenvalue and eigenvector formulas used throughout the spectral lemmas."},{"cited_title":"Cambridge University Press, ??? (2012)","cited_arxiv_id":null,"evidence_quote":"Horn and Johnson's Matrix Analysis provides the Sylvester equation solvability theorem used to prove existence and uniqueness of self-similar profiles."},{"cited_title":"Journal of Differential Equations 247(9), 2620–2636 (2009) https://doi.org/10.1016/j.jde.2009.07.024","cited_arxiv_id":null,"evidence_quote":"Lin and Tsai's length-preserving linear flow frames the problem of evolving one curve to another, which the curvature-difference flow addresses in the semi-discrete hyperbolic setting."}],"review_version":1}