{"id":"67b7db1c-9bdf-4f7e-af93-27246b655385","arxiv_id":"1908.02863","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For slightly deformed triangles and for triangles with a small added potential, the L2 Neumann data mass on each side is the side length divided by the triangle area, up to an error of order epsilon.","lead":"This paper proves that a special balance in how wave-like solutions of a quantum triangle touch its boundary survives small changes to the shape or to a background potential. For each side, the boundary activity stays close to a simple fraction given by side length over area, with error set by the perturbation size.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lemma 2.1's obtuse-case proof has a sign/cancellation gap: with X=(y-a2/l x)∂y, the B-boundary contribution is negative when a1>a2, so the O(1) total does not control the positive C' integral.","rationale":"The reader's weakest-assumption analysis correctly identified Lemma 2.1 as the load-bearing step and noted sign issues in the obtuse case. My stress-test confirms this and isolates a specific mechanism: in the obtuse part of Section 4, the vector field X=(y-a2/l x)∂y produces a boundary contribution on B whose sign depends on a1-a2. For a1>a2, the B-term is negative and can cancel the positive C'-term in the boundary identity, so the O(1) bound on the total does not imply the needed bound on the C' integral. This is not merely a cosmetic issue: equations (2.1)-(2.3) use the boundedness of IC' to convert the exact commutator identity into the O(ε) estimates for the three sides, and Theorem 1 states the result for all obtuse triangles. A likely correction is to use X=(y-a1/l x)∂y, which vanishes on B and restores the acute-style argument; the manuscript's 'a2' is inconsistent with that. Because this is a concrete gap in the proof as written, rather than a demonstrated counterexample, the conditional verdict remains appropriate. The potential-theorem existence issue and the missing uniform-bound lemma for the potential case are secondary; they are also presentation gaps but do not change the verdict. The paper's explicit commutator identities and the recovery of the exact triangle identity at ε=0 give independent support for the claimed result, but the obtuse lemma needs to be fixed before the proof is complete.","tokens_in":11527,"tokens_out":25717,"duration_ms":302255,"concrete_test":"Recompute the obtuse part of Lemma 2.1 with an explicit obtuse triangle, e.g. l=1, a1=2, a2=1, and g=0. For X=(y-a2 x)∂y, write out the boundary integral ∫_{∂D} (hXu)(h∂νu) dS: on B it equals -((a1-a2)/b)∫_{B} x|h∂νu|² dS, which is negative, while on C' it is positive. Then check whether the displayed bound ∫_{C'∩{x≥δ}} |h∂νu|² dS = O(1) follows from the commutator identity; it does not, because the total can be small due to cancellation. Then repeat with X=(y-a1 x)∂y; since X vanishes on B and gives a positive C' contribution, the acute-style bound is recovered. This determines whether the published a2 is a load-bearing error or a harmless typo.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of Theorem 1 depends on Lemma 2.1, which supplies the uniform bound ∫_{C'} |h∂νu|² dS ≤ Γ used in (2.1)-(2.3) to control the O(ε) errors for IA, IB, and IC'. The acute-case proof of Lemma 2.1 works because the B-term in the boundary identity has the same sign as the C'-term. The obtuse-case proof as written does not have this property. In Section 4, the vector field X=(y-a2/l x)∂y is used. On B, y=a1/l x and h∂y=-(l/b)h∂ν, so Xu·h∂νu = -((a1-a2)x/b)|h∂νu|². When a1>a2, this B-term is negative, while the C'-term is positive. The commutator identity only gives an O(1) bound on the total boundary integral; if the negative B-term cancels a large positive C'-term, the conclusion ∫_{C'∩{x≥δ}} |h∂νu|² = O(1) does not follow. The fix is likely to replace a2 by a1 in X, which makes X vanish on B and leaves a positive C'-term, but this correction is not present in the manuscript. Since Theorem 1 states the result for all obtuse triangles, and the a1>a2 case is not excluded, the proof of the load-bearing bound is incomplete as written.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper claims two stability results for Christianson's exact Neumann data mass equidistribution on triangles. Theorem 1 replaces one side of a triangle by a smooth curve whose graph distance from the original side is at most epsilon and asserts that, for every normalized Dirichlet eigenfunction of -h^2 Delta - 1, the integral of |h d_nu u|^2 over each side equals the side length divided by the area of the perturbed domain, up to O(epsilon), with constants uniform in h and in the perturbation. Theorem 2 proves the same form of equidistribution, up to O(epsilon), when an epsilon-small potential is added to the Laplacian on a fixed triangle. The proof is based on commutators with affine vector fields and on an a priori bound, Lemma 2.1, for the Neumann data on the perturbed side.","tokens_in":11812,"tokens_out":16220,"duration_ms":175450,"significance":"If the proof is completed, this is a genuinely useful stability statement: the exact triangle identity from [Chr17] is shown to be robust under small boundary and potential perturbations, uniformly in the spectral parameter h and in the perturbation. The argument is elementary, self-contained apart from the zero-perturbation baseline, and no quantity in the conclusion is used as an input, so the derivation is not circular. The main value is as a benchmark for how rigid the triangle equidistribution law is, and the methods could plausibly extend to higher dimensions. The significance is moderate but solid for a specialized analysis paper.","major_comments":[{"comment":"In the obtuse case of Lemma 2.1, the vector field X=(y-(a2/l)x)d_y does not have the positivity property used in the acute case. On the side B, parametrized by y=(a1/l)x, the boundary contribution in (4.1) is ((a2-a1)x/b)|h d_nu u|^2 dS, which is negative when a1>a2, while the C' contribution f(x) gamma^{-1}|h d_nu u|^2 is positive for x>=delta. The commutator identity bounds only the signed sum of the boundary terms, so the conclusion int_{C' cap {x>=delta}} |h d_nu u|^2 = O(1) does not follow for all obtuse triangles as written. A likely correction is X=(y-(a1/l)x)d_y, which vanishes on B and makes the C' term positive, but this is not the vector field used in the manuscript.","section":"Section 4, obtuse case"},{"comment":"The vector field Y=psi(x)(d_x+(a2/l)d_y) is asserted to be tangential to the side B, but B has slope a1/l, so a tangential field is d_x+(a1/l)d_y. As written, Yu does not vanish on B. The subsequent computation on C' uses hYu=psi(x)(-f'/gamma + a1/(l gamma))h d_nu u, which corresponds to the field d_x+(a1/l)d_y, not the displayed one. This inconsistency must be corrected and the boundary estimates rechecked, since the near-corner control of int_{C' cap {0<=x<=delta}} |h d_nu u|^2 depends on it.","section":"Section 4, obtuse case"},{"comment":"After equation (3.2), the identity has the form 2+O(epsilon)(1+|m|+|n|) equal to the boundary sum, and the text says 'the rest of the proof proceeds exactly as the proof of Theorem 1.' However, deriving the analogues of (2.2) and (2.3) requires differentiating the identity with respect to m and n. Since the error term depends on m and n, one must justify that differentiation produces only O(epsilon) corrections. This is likely fixable by a direct first-order Taylor expansion in m and n, but it is not written and is load-bearing for Theorem 2.","section":"Section 3, proof of Theorem 2"}],"minor_comments":[{"comment":"The parametrization of C' in the obtuse case is written as f(x)=a2+a1/l x+g(x), which is inconsistent with the later derivative f'(x)=(a1+a2)/l+g'(x). Presumably f(x)=((a1+a2)/l)x+g(x), and the displayed formula should be corrected.","section":"Section 2.3"},{"comment":"The notation in the proof of Lemma 2.1 is inconsistent: the slope of C' is a2/l+g'(x), but the text uses a/b in the lower-bound estimates, and the B-boundary coefficient is a1/b, not a2/b. Please define all symbols consistently so that the constants in (4.2) and the near-corner estimate can be verified.","section":"Section 4, acute case"},{"comment":"The remark says that g does not have to be of the form epsilon g_tilde, but the domain construction still requires g(0)=g(l)=0 for the perturbed side to close the domain; this endpoint condition should be stated explicitly for the more general g.","section":"Remark 1.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is likely salvageable with a focused revision. The main obstruction is the obtuse-case proof of Lemma 2.1, where the displayed vector fields appear to be mistyped; replacing a2 by a1 in the two vector fields and rechecking the boundary terms would probably repair the proof. The remainder of the argument is coherent and the paper is short, so I expect a revised version to be acceptable. There are no concerns about novelty or citation practice."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper takes Christianson's exact equidistribution identity for Dirichlet eigenfunctions on triangles and asks whether it survives two small perturbations: deforming one side by O(ε) and adding an O(ε) potential. The answer is the expected one—the per-side Neumann data mass stays within O(ε) of length/area—and the proof is mostly a direct commutator computation in the spirit of the original. That is a genuine, if modest, new result; I do not see it in the cited papers.\n\nThe acute-triangle case is worked out in detail and looks correct. The bookkeeping with the vector field X=(x+m)∂x+(y+n)∂y, the derivation of equations (2.1)–(2.3), and the O(ε) control on the C′ side all hold together. The explicit independence of the error from h is a nice touch.\n\nThe trouble is in the obtuse case of Lemma 2.1, which supplies the a priori bound ∫_{C′} |h∂νu|² dS ≤ Γ that everything else leans on. The proof for obtuse triangles uses X=(y−a2/l x)∂y. On side B, y=a1/l x, so Xu·h∂νu = −(a1−a2)x/b |h∂νu|². If a1>a2, that term is negative, while the C′ term is only O(ε) (since Xu on C′ is g(x)∂y, with |g|≤ε). The commutator identity only gives an O(1) bound on the whole boundary sum, so a large positive C′ integral could be cancelled by the negative B term. The remark in the text that \"the same argument as in the acute case\" works does not address this. A likely fix is to choose X=(y−a1/l x)∂y, which kills the B term and leaves a positive C′ term, but that choice is not made. So the proof of the load-bearing lemma is incomplete as written for at least some obtuse triangles.\n\nA second, smaller issue is Theorem 2. The statement assumes there are normalized eigenfunctions of −h²Δ+w_ε with eigenvalue exactly 1. For a fixed h and generic potential this is not automatic; one would need to work with spectral windows or quasi-modes. The proof uses only |w|<ε and |∇w|<ε, so the gap is probably easy to repair, but it is not just notational.\n\nThere are also a few typos and garbled display equations (e.g., in Section 4 the expression for f′ looks off), but nothing that suggests dishonesty or circular reasoning. The paper is an honest, competent extension of an existing program.\n\nWho is it for: spectral geometers working on eigenfunction restrictions and boundary traces. With a corrected Lemma 2.1 and a clarified spectral assumption, it would be a solid contribution. I would send it to peer review rather than desk reject, but I would insist on the obtuse case being fixed. For my own work, I would not cite the current version yet.","headline":"Solid extension of an exact equidistribution result; acute case works, but the obtuse case has a real gap in the key lemma and Theorem 2 has a spectral assumption issue.","tokens_in":12337,"tokens_out":6340,"would_cite":false,"duration_ms":60367,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35P20","58J51","35J25"],"pacs":[],"model":"deepseek-v4-flash","headline":"The exact triangle identity for Neumann data mass is stable: an epsilon-small bump on one side, or an epsilon-small potential inside, changes each side's boundary L2 mass from side-length over area by at most O(epsilon), uniformly in…","keywords":["Neumann data mass","Dirichlet eigenfunctions","semiclassical analysis","perturbed triangles","equidistribution on triangles","boundary traces","commutator method","small potential perturbations"],"falsifier":"Fix a right triangle and set $g(x)=\\epsilon\\sin(\\pi x/l)$ on one side, then compute $I_A,I_B,I_{C'}$ for a sequence of normalized Dirichlet eigenfunctions as $h\\to 0$ at one small $\\epsilon$. If for some eigenfunction $I_{C'}$ grows without bound in $h$, or if any of the three integrals deviates from the corresponding side length divided by the area by more than a fixed multiple of $\\epsilon$, then Lemma 2.1 or the uniformity in Theorem 1 is false; a numerical search over low-lying eigenfunctions would settle it.","tokens_in":11297,"feed_emoji":"📐","tokens_out":17506,"duration_ms":154533,"temperature":0.7,"pith_summary":"The paper proves that the exact equidistribution of Neumann data on triangles - each side's boundary $L^2$ mass of the semiclassical normal derivative equals the side's length divided by the area - survives small perturbations of both the domain and the operator. If one side of a triangle is replaced by a smooth curve that is $\\epsilon$-close to the original side in $C^1$ norm, then every normalized Dirichlet eigenfunction of $-h^2\\Delta-1$ on the new domain satisfies $\\int_A |h\\partial_\\nu u|^2\\,dS = a/\\operatorname{Area}(D)+O(\\epsilon)$, with the same statement on the other sides. The same conclusion holds when an $\\epsilon$-small potential is added to the Laplacian on a fixed triangle. The error constant is independent of $h$ and of the particular perturbation, so the law is not an artifact of straight edges.","feed_headline":"Perturbed triangles still split boundary mass by side length","feed_subtitle":"A smooth epsilon-bump on one side changes each side's boundary L2 mass from length/area by at most O(epsilon).","key_machinery":"The central mechanism is the commutator identity $[-h^2\\Delta-1,X]=-2h^2\\Delta$ for the affine vector field $X=(x+m)\\partial_x+(y+n)\\partial_y$. Because the eigenfunction is normalized, this identity converts the interior norm into the boundary identity $2=\\int_{\\partial D}(hXu)(h\\partial_\\nu u)\\,dS$. The Dirichlet condition makes the tangential part of $hXu$ vanish, so on each side $hXu$ becomes a known multiple of $h\\partial_\\nu u$; matching the three sides yields a $3\\times 3$ linear system for the three masses $I_A,I_B,I_{C'}$. The smallness assumptions on $g$ and $g'$ enter as $O(\\epsilon)$ remainder terms, and Lemma 2.1, proved by separate vector fields $y\\partial_y$ and $\\psi(x)\\partial_x$, provides the uniform bound on $I_{C'}$ needed to close the estimates.","core_discovery":"On a triangle with one side $C$ replaced by the graph $f(x)=a_2 x/l+g(x)$, where $|g|,|g'|\\le\\epsilon$ and $g(0)=g(l)=0$, the paper proves that for the Dirichlet eigenfunction problem $(-h^2\\Delta-1)u=0$ with $\\|u\\|_{L^2}=1$, the Neumann data masses obey $I_A=a/\\operatorname{Area}(D)+O(\\epsilon)$, $I_B=b/\\operatorname{Area}(D)+O(\\epsilon)$, and $I_{C'}=\\ell(C')/\\operatorname{Area}(D)+O(\\epsilon)$, where $I_\\cdot=\\int_\\cdot|h\\partial_\\nu u|^2\\,dS$ and the constants are uniform in $h$ and in the perturbation. For a fixed triangle, the same equalities hold up to $O(\\epsilon)$ when the Laplacian is replaced by $-h^2\\Delta+w_\\epsilon$ with $|w_\\epsilon|,|\\nabla w_\\epsilon|\\le\\epsilon$. Thus the exact triangle identity is stable: small geometric or operator perturbations change each side's boundary mass only by the size of the perturbation, never by a factor that grows as the frequency goes up.","pith_inferences":["A direct consequence the authors leave implicit: because the three boundary masses are measurable, one can recover the area and the three side lengths of an unknown near-triangular domain up to $O(\\epsilon)$ from Neumann data alone, a stable inverse problem.","The mechanism is the commutator identity for a linear scaling vector field, so the result should extend to any domain with three boundary pieces on which the Dirichlet condition turns $Xu$ into a tangential multiple; a polygon with more sides would likely give the analogous statement with side lengths weighted by the $C^1$ mismatch.","Only $|g|,|g'|\\le\\epsilon$ are used outside the cut-off argument in Lemma 2.1, so perturbations with large second derivative but small height and slope may probe whether curvature is hidden in the boundedness lemma.","If the $O(\\epsilon)$ error is not just an estimate artifact, the boundary mass distribution on a slightly non-triangular domain could serve as a sensitive probe of the domain's area and side geometry, with no need to know individual eigenfunctions."],"forward_implications":["The exact triangle identity extends to a neighborhood of triangles in the $C^1$ topology on the boundary, so the formula is stable under small deformations, not just for polygons.","Because the error is uniform in $h$, the result applies to every eigenfunction in the whole spectrum, including high-frequency modes where boundary concentration might otherwise be expected.","Adding an $\\epsilon$-small potential does not break the law, so the same equidistribution holds for perturbed Hamiltonians $-h^2\\Delta+w_\\epsilon$; the formula is robust for operator perturbations as well as domain perturbations.","The proof supplies explicit bounds from which the $O(\\epsilon)$ constants could in principle be extracted from the side lengths and the a priori constant $\\Gamma$ of Lemma 2.1.","The authors expect the analogous statement in all dimensions $n\\ge3$ for simplices with perturbed faces, following the same commutator scheme."],"supporting_citations":[{"why":"Supplies the exact Neumann-data equidistribution identity on triangles and the commutator/vector-field proof template that this paper perturbs.","marker":"[Chr17]"},{"why":"Gives the same exact identity on simplices in dimensions $n\\ge3$, the result whose perturbed analogue the authors expect to extend.","marker":"[Chr18]"},{"why":"Supplies the $O(1)$ scale for Neumann data along hypersurfaces that Lemma 2.1 must preserve under perturbation.","marker":"[CHT15]"},{"why":"Gives the boundary-value quantum ergodicity asymptotic that the exact triangle identity matches, framing why a whole-sequence exact formula is significant.","marker":"[HZ04]"}],"fun_headline_variants":["Epsilon-bump on triangle side leaves Neumann mass stable","Boundary L2 mass stays length/area under tiny perturbations","Small perturbations don't shift Neumann mass on triangles","Triangle Neumann mass robust to epsilon perturbations","Boundary mass split survives side bumps on triangles"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole proof rests on Lemma 2.1: the integral over the perturbed side $C'$ of $|h\\partial_\\nu u|^2$ must stay bounded by a constant independent of $h$ and $\\epsilon$; if that bound fails, the $O(\\epsilon)$ errors in the linear system cannot be controlled and the equidistribution formulas collapse.","fun_headline_variants_meta":{"raw":{"variants":["Epsilon-bump on triangle side leaves Neumann mass stable","Boundary L2 mass stays length/area under tiny perturbations","Small perturbations don't shift Neumann mass on triangles","Triangle Neumann mass robust to epsilon perturbations","Boundary mass split survives side bumps on triangles"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000599,"raw_usage":{"total_tokens":2795,"prompt_tokens":938,"completion_tokens":1857,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":554,"completion_tokens_details":{"reasoning_tokens":1783}},"tokens_in":554,"tokens_out":1857,"duration_ms":12178,"temperature":1.0,"reasoning_tokens":1783,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:32:59.715146+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fix a right triangle and set $g(x)=\\epsilon\\sin(\\pi x/l)$ on one side, then compute $I_A,I_B,I_{C'}$ for a sequence of normalized Dirichlet eigenfunctions as $h\\to 0$ at one small $\\epsilon$. If for some eigenfunction $I_{C'}$ grows without bound in $h$, or if any of the three integrals deviates from the corresponding side length divided by the area by more than a fixed multiple of $\\epsilon$, then Lemma 2.1 or the uniformity in Theorem 1 is false; a numerical search over low-lying eigenfunctions would settle it.","supporting_citations":[],"review_version":1}