{"id":"03e9f9eb-ef9d-46e1-927a-ea14e41f16f9","arxiv_id":"1908.05096","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"The elastic Dirichlet-to-Neumann map determines the real-analytic metric up to isometry, and its heat trace expansion gives explicit spectral invariants such as boundary volume and total mean curvature.","lead":"An analysis of the boundary measurement map for the equations of elasticity shows that, for real-analytic curved bodies, the surface map determines the interior geometry up to isometry. It also computes new spectral coefficients that encode boundary volume and curvature from elastic vibration frequencies.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The leading heat coefficient a0 in Eqs. (5.22)/(1.9) is arithmetically wrong: the e^{-2μ|ξ|t} term contributes a factor 1/2^{n-1}, not 1, so the stated Weyl-law constant and boundary-volume invariant in Theorem 1.3 are incorrect.","rationale":"I read the paper as making two central claims: the elastic analogue of Lee–Uhlmann boundary determination (Theorem 1.1) and explicit elastic Steklov heat invariants (Theorem 1.3). The reader's weakest-assumption analysis points to the Lee–Uhlmann–Myers extension lemma. That is a real dependency, but it is also an explicit hypothesis of Theorem 1.1 and is quoted from a standard source, so I would not call it the single most load-bearing flaw. The sharper, verifiable problem is the arithmetic in the leading heat coefficient a0: the ξ-integral in Eq. (5.22) is written inconsistently with the standard radial integral formula, and the same error is repeated in Eq. (1.9). This affects the claimed boundary-volume spectral invariant and the derived Weyl law, which are presented as answers to Problem B. The metric-determination theorem is not refuted by this check, but the spectral half of the paper is not correct as stated. Since the overall verdict CONDITIONAL already captures the need for correction, I recommend leaving the verdict unchanged.","tokens_in":69246,"tokens_out":26484,"duration_ms":245336,"concrete_test":"Recompute the t→0+ trace coefficient for n=2, λ=μ=1. Eq. (5.20) then gives the contour trace e^{-2|ξ|t} + e^{-|ξ|t} (since (λ+μ)/(λ+3μ)=1/2 and n−2=0). The exact integral is (1/2π)∫_R(e^{-2|ξ|t}+e^{-|ξ|t}) dξ = (1/2π)(1/t + 2/t) = 3/(2π t), so a0 = 3/(2π). The paper's Eq. (1.9)/(5.22) yields a0 = 2/π. If the independent evaluation confirms 3/(2π), then the leading coefficient in Theorem 1.3 and the Weyl constant in Remark 5.2 must be corrected.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's strongest spectral claim, Theorem 1.3 and Remark 5.2, rests on the evaluation of the trace integral around Eq. (5.20). In Eq. (5.22), the three exponentials e^{-2μ|ξ'|t}, e^{-2μ(λ+μ)|ξ'|t/(λ+3μ)}, and (n−2)e^{-μ|ξ'|t} are integrated over R^{n−1}. The first integral is vol(S^{n−2}) Γ(n−1) (2μ t)^{1−n} = (n−2)!/(2^{n−1} μ^{n−1} t^{n−1}). After factoring t^{1−n} (n−2)! vol(S^{n−2})/((2π)^{n−1} μ^{n−1}), the first exponential contributes 1/2^{n−1}, not 1. Equation (5.22), and the identical Eq. (1.9), instead put 1 in the first slot of the curly bracket, so a0 is too large by a factor 2^{n−1}. For n=2 this makes the claimed boundary-volume coefficient wrong by a factor 2. This is an internal arithmetic inconsistency, independent of any geometric hypothesis, and it directly affects the advertised spectral invariants. It does not by itself invalidate Theorem 1.1, whose global step is governed by the imported Lee–Uhlmann–Myers extension lemma; that lemma is a stated hypothesis rather than the weakest point of the proof.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the elastic Dirichlet-to-Neumann map for the Lamé system on a compact Riemannian manifold with boundary. It derives an explicit full symbol of the map by a matrix factorization of the Lamé operator, proves that the map determines the real-analytic metric up to isometry under the Lee–Uhlmann hypotheses (Theorem 1.1), and computes the first n coefficients of the heat-trace expansion of the associated semigroup (Theorem 1.3). The boundary-recovery argument is constructive: it recovers the Taylor series of the metric in boundary normal coordinates from the symbol of the Dirichlet-to-Neumann map. The spectral part evaluates contour and radial integrals to extract boundary-volume and curvature invariants.","tokens_in":69531,"tokens_out":8011,"duration_ms":80707,"significance":"If the massive algebraic derivation is correct, Theorem 1.1 would settle the real-analytic anisotropic analogue of the Lee–Uhlmann boundary determination problem for elasticity, including the Euclidean-domain corollary, and would be the first result of this type for the elastic Dirichlet-to-Neumann map. The spectral calculation is also conceptually novel: it gives explicit formulas for all coefficients a_0, ..., a_{n-1} rather than only the leading term. A genuine strength is that the construction is parameter-free: the constants s_1,...,s_5 and the tilded constants are fixed by solving internal algebraic equations, not fitted to any target. However, the leading heat coefficient in Theorem 1.3 contains an arithmetic inconsistency, so the spectral claims as stated are not yet reliable; the length and density of the symbolic computations make independent verification essential.","major_comments":[{"comment":"The first term in the curly bracket of Eq. (5.22), and the identical expression in Eq. (1.9), is inconsistent with the preceding integral evaluation. Since ∫_{R^{n-1}} e^{-2μ|ξ'|t} dξ' = vol(S^{n-2}) Γ(n-1) (2μt)^{-(n-1)}, this contribution equals t^{1-n} (n-2)! vol(S^{n-2}) / ((2π)^{n-1} μ^{n-1}) · 2^{-(n-1)}, not the printed factor 1. The exponential with rate 2μ(λ+μ)/(λ+3μ) also carries a factor 2^{-(n-1)} in this normalization; the factor printed in Eq. (1.9), if read as 2^{n-1}, is also inconsistent with the integrand. Consequently the coefficient a_0 and the Weyl constant in Eq. (5.31) are too large by a factor 2^{n-1}; for n=2 the boundary-volume invariant is wrong by a factor of 2. This is an internal arithmetic inconsistency that invalidates Theorem 1.3 and Remark 5.2 as printed.","section":"Eq. (5.22) and Eq. (1.9), Remark 5.2"},{"comment":"The global step of Theorem 1.1 is not self-contained: Lemma 4.3 (Lee–Uhlmann–Myers) is quoted without proof, and the hypotheses of strong convexity or extendability together with the fundamental-group condition are imported from [42] rather than proved in the elastic setting. This is a legitimate use of the literature, but the paper should state precisely where Lemma 4.3 is proved and confirm that the topological condition π(Ω,∂Ω)=0 is exactly the condition π_1(Ω,∂Ω)=0 used by Lee and Uhlmann. Because Lemma 4.2 only constructs an isometry on a neighborhood of the boundary, Theorem 1.1 cannot be considered established without this imported result.","section":"Proof of Theorem 1.1 and Lemma 4.3"},{"comment":"The simplification leading from Eq. (5.28) to Eq. (5.29) is the computational core of Theorem 1.3 and depends on dozens of intermediate displayed formulas and the constants s_i and s̃_i. The paper announces the result but does not provide a verifiable derivation, a computer-algebra file, or an independent check of the simplifications. Given the demonstrated arithmetic slip in the leading coefficient a_0, a verified symbolic computation or a detailed derivation of (5.29) and (5.31) should be required before the spectral invariants are accepted.","section":"Eqs. (5.28)–(5.29) and the symbolic reduction"}],"minor_comments":[{"comment":"The reference to “Theorem 5.3” appears to be a typo; the statement being discussed is Theorem 1.1 (or possibly Theorem 1.3), not a theorem in Section 5.","section":"Remark 4.4(ii)"},{"comment":"Remark 4.5 claims the results hold for smooth Lamé parameters with μ<0, although ellipticity elsewhere requires μ>0; this is likely a typo, and the claimed extension to variable parameters needs either a proof or deletion.","section":"Remark 4.5"},{"comment":"The heat-trace sum in Theorem 1.3 begins at k=0 while the eigenvalues are indexed from k=1; the indexing should be harmonized.","section":"Section 5, Theorem 1.3"},{"comment":"Several displayed formulas contain garbled inequalities such as “n /greaterorequalslant2”, and the notation for the sphere volume and Gamma function is crowded; the manuscript would benefit from a careful typesetting pass.","section":"Section 5, displayed formulas"},{"comment":"The paper alternates between π(Ω,∂Ω)=0 and π_1(Ω,∂Ω)=0; the notation should be made consistent and the meaning of the topological condition stated explicitly in both places.","section":"Lemma 4.3 and Theorem 1.1"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the arithmetic error in Eq. (5.22)/(1.9) is unambiguous and affects a headline result (Theorem 1.3 and Remark 5.2). It is fixable, but the author should recheck the entire spectral section and provide a computer-algebra verification of the reductions leading to (5.29) and (5.31). The boundary-determination part (Theorem 1.1) appears substantially more reliable, although the computations in Sections 3–4 are very long and merit a careful independent reading. I recommend major revision rather than rejection because the central boundary-recovery claim is defensible and the spectral error can be corrected within the scope of the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: Theorem 1.1 is a genuine extension of Lee-Uhlmann to the Lamé system and looks like the strongest part of the paper. The construction of the matrix-valued full symbol of the elastic Dirichlet-to-Neumann map is substantial: the factorization of the Lamé operator, the ring-theoretic solution of the quadratic symbol equation, and the Sylvester-equation recursions are all real work. On reading, the boundary recovery argument is coherent; the algebraic inequality on λ and μ is not a cosmetic condition. The global step relies on the Lee-Uhlmann-Myers extension lemma, which is stated and imported, so the theorem is only as strong as those hypotheses, but that is not a hidden flaw.\n\nThe spectral part is in trouble. There is a concrete arithmetic error in the leading heat coefficient. Equation (5.22) integrates three exponentials over R^{n-1}. Using the paper's own formula (5.21), ∫ e^{-2μ|ξ'|t} dξ' = vol(S^{n-2}) Γ(n-1) (2μt)^{-(n-1)} = vol(S^{n-2})(n-2)! / (2^{n-1} μ^{n-1} t^{n-1}). After factoring the common prefactor, the first curly-bracket term should be 1/2^{n-1}, not 1. The same error appears in Eq. (1.9), so a0, the advertised boundary-volume invariant, and the Weyl law in Remark 5.2 are off by a factor of 2^{n-1} in that term. This is not a cosmetic typo; it undermines Theorem 1.3 as stated. It does not by itself damage Theorem 1.1.\n\nTwo smaller caveats. The abstract says “explicitly obtain all coefficients a0,...,a_{n-1},” but the text gives closed forms for a0 and a1 only; for m ≥ 2 it gives a procedure plus a structural description of a2. That is a mismatch between advertising and content. And given the size of the algebra, the absence of machine-checked calculation is felt; the author thanks someone for MATLAB checks, but no code or data is provided. I would want the corrected spectral computation independently verified before relying on any of the invariants.\n\nThe boundary-determination theorem is citable once the spectral section is corrected. It deserves a serious referee; a good referee can verify Theorem 1.1 and force the arithmetic fix.","headline":"The elastic boundary determination result is real and worth refereeing; the spectral invariants have a concrete arithmetic error and need correction.","tokens_in":70084,"tokens_out":3016,"would_cite":false,"duration_ms":31459,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["53C24","74B05","35K50","35P20","35S05"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that on a compact real-analytic Riemannian manifold with boundary, the elastic Dirichlet-to-Neumann map determines the metric up to a boundary-fixing isometry under Lamé and extension hypotheses.","keywords":["elastic Lamé operator","elastic Dirichlet-to-Neumann map","isometric uniqueness","real-analytic metric","boundary determination","elastic Steklov eigenvalues","heat trace asymptotics","spectral invariants"],"falsifier":"Find or compute two real-analytic metrics on a compact real-analytic manifold with real-analytic boundary that have the same elastic Dirichlet-to-Neumann map but admit no global boundary-fixing isometry, while the extension hypotheses fail; or, at Lamé parameters satisfying $(n-1)\\lambda^3+(4n-2)\\lambda^2\\mu+(n+5)\\lambda\\mu^2+(14-8n)\\mu^3=0$, exhibit two real-analytic metrics whose DtN symbols agree but whose first normal derivatives at the boundary differ.","tokens_in":68997,"feed_emoji":"⚙️","tokens_out":8524,"duration_ms":80154,"temperature":0.7,"pith_summary":"The paper studies the elastic Dirichlet-to-Neumann map $\\Xi_g$, the operator sending boundary displacements to boundary tractions for the stationary Lamé system on a compact Riemannian manifold. It proves that under real-analyticity plus a strong-convexity or extendability condition, the full map $\\Xi_g$ determines the metric $g$ up to a real-analytic diffeomorphism that fixes the boundary, for all dimensions $n\\ge 2$ and Lamé constants outside one explicit algebraic exceptional set. It also computes the first $n-1$ heat-trace coefficients of the elastic Steklov eigenvalues, showing that boundary volume, total mean curvature, and other total boundary curvatures are spectral invariants measurable from the elastic Steklov spectrum. The results transfer the scalar Laplace boundary-determination theorem and its spectral consequences to a non-Laplace-type operator, answering two open problems posed about the elastic boundary map.","feed_headline":"Elastic boundary data fix the metric up to isometry","feed_subtitle":"Displacement-to-traction data determine the metric; its eigenvalues yield boundary volume and curvature.","key_machinery":"The load-bearing mechanism is the factorization of the elastic Lamé operator in boundary normal coordinates as $$\\partial_{x_n}^2 I + B\\partial_{x_n}+C = (\\partial_{x_n}I+B-Q)(\\partial_{x_n}I+Q),$$ where $Q$ is a pseudodifferential operator with symbol $q\\sim q_1+q_0+q_{-1}+\\cdots$. The principal symbol $q_1$ is found by solving the matrix quadratic $q_1^2-b_1q_1+c_2=0$ inside an invariant subring generated by three geometric matrices over a coefficient ring of diagonal matrices; the lower-order symbols solve Sylvester equations $(q_1-b_1)q_{j-1}+q_{j-1}q_1=E_j$, with explicit solutions obtained by inverting the Kronecker-sum matrix $U=I_n\\otimes L+M^t\\otimes I_n$. This full symbol then defines $\\Xi_g$, and its subprincipal terms recover $g$ and its normal derivatives inductively. The same symbol feeds the resolvent $(\\Xi_g-\\tau I)^{-1}$, whose trace over a contour gives the heat-trace coefficients.","core_discovery":"The central discovery is that the matrix-valued full symbol of $\\Xi_g$ carries every order of normal-derivative information about the boundary metric. Explicitly, the paper proves that for $\\mu>0$, $\\lambda+\\mu\\ge 0$, and the nondegeneracy condition $(n-1)\\lambda^3+(4n-2)\\lambda^2\\mu+(n+5)\\lambda\\mu^2+(14-8n)\\mu^3\\ne 0$, the symbol determines the boundary value of $g$ and all its normal derivatives at every boundary point; with real-analyticity this gives a local boundary-fixing isometry, and the extension lemma upgrades it to a global isometry, yielding $g=\\varrho^*\\tilde g$. On the spectral side, the paper derives the asymptotic expansion $\\sum_{k=1}^\\infty e^{-t\\tau_k}\\sim \\sum_{m=0}^{n-1}a_m t^{m+1-n}+o(1)$ and gives explicit formulas for $a_0$ and $a_1$, the former proportional to boundary volume and the latter to the integral of total mean curvature, with all $a_m$ obtainable by an explicit recursive procedure.","pith_inferences":["A testable consequence the author leaves implicit: because the recovery of the metric is by explicit Taylor series in boundary normal coordinates, the theorem implies a practical local reconstruction algorithm for real-analytic elastic media that measures displacements and tractions on the boundary and then extends the recovered collar metric.","The exceptional Lamé polynomial that must not vanish is a natural place to look for genuine counterexamples: at those parameters the subprincipal symbol may fail to determine the first normal derivative of the metric, and two metrics with the same full DtN symbol might differ to first order across the boundary.","The spectral-invariant method is transferable in principle to other non-Laplace-type boundary operators, such as anisotropic elasticity, thermoelasticity, or poroelasticity, where no explicit heat invariants are known; the same subring-and-Sylvester recipe could yield their first $n-1$ invariants.","Because $a_2$ is written in terms of scalar and Ricci curvatures of both the domain and boundary, the elastic Steklov spectrum may give access to interior curvature once sufficiently many heat coefficients are measured, not just boundary geometry."],"forward_implications":["Equality of elastic DtN maps $\\Xi_g=\\Xi_{\\tilde g}$ forces $g=\\varrho^*\\tilde g$ for a boundary-fixing real-analytic diffeomorphism, so all isometry-invariant interior geometry is encoded in the elastic boundary response.","In the special case of a simply connected bounded real-analytic domain in Euclidean space, any real-analytic metric with the same displacement-to-traction map as the Euclidean metric is isometric to it.","The leading heat coefficient $a_0$ determines the boundary volume $\\mathrm{vol}(\\partial\\Omega)$ from the elastic Steklov spectrum, and $a_1$ determines the total mean curvature and related boundary curvature integrals.","The leading-term Tauberian argument yields the Weyl law $N(\\tau)\\sim C(n,\\mu,\\lambda)\\,\\mathrm{vol}(\\partial\\Omega)\\,\\tau^{n-1}$ for the counting function of elastic Steklov eigenvalues.","The same explicit recovery procedure is stated to work for variable Lamé parameters $\\mu>0$, $\\lambda+\\mu\\ge0$, so the uniqueness statement is not limited to constant coefficients."],"supporting_citations":[{"why":"provides the scalar boundary-determination prototype and the boundary-to-global extension lemma used to complete the isometric uniqueness proof.","marker":"[42]"},{"why":"supplies the integration-by-parts identity and elastic-energy variational formulation used to derive the Lamé operator on a Riemannian manifold.","marker":"[15]"},{"why":"gives the heat-trace expansion method and the sphere integrals used to evaluate the spectral coefficients.","marker":"[43]"},{"why":"supplies the Sylvester-equation theory and Kronecker-product formulation used to solve for the lower-order symbols of Q.","marker":"[8]"},{"why":"provides the pseudodifferential symbol calculus and contour-resolvent formalism used in both the symbol recovery and spectral trace computations.","marker":"[67]"},{"why":"underlies the geometric extension principle that turns a collar isometry into a global boundary-fixing diffeomorphism.","marker":"[53]"},{"why":"supplies the classical factorization method for scalar boundary value problems that the elastic factorization adapts to the matrix Lamé system.","marker":"[69]"}],"fun_headline_variants":["Elastic boundary map fixes metric up to isometry","Elastic Steklov eigenvalues expose boundary geometry","Boundary data reveal metric and spectral invariants","Elastic boundary response determines manifold shape","Uniqueness in elastic boundary inverse problem"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof stands on a global extension lemma: after the DtN symbol recovers the metric only near the boundary, that local isometry must extend across the whole manifold, which requires either strong convexity of the manifold with respect to both metrics or extendability of one metric, together with $\\pi(\\Omega,\\partial\\Omega)=0$.","fun_headline_variants_meta":{"raw":{"variants":["Elastic boundary map fixes metric up to isometry","Elastic Steklov eigenvalues expose boundary geometry","Boundary data reveal metric and spectral invariants","Elastic boundary response determines manifold shape","Uniqueness in elastic boundary inverse problem"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00028,"raw_usage":{"total_tokens":1727,"prompt_tokens":1080,"completion_tokens":647,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":696,"completion_tokens_details":{"reasoning_tokens":578}},"tokens_in":696,"tokens_out":647,"duration_ms":6056,"temperature":1.0,"reasoning_tokens":578,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:26:30.411823+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Find or compute two real-analytic metrics on a compact real-analytic manifold with real-analytic boundary that have the same elastic Dirichlet-to-Neumann map but admit no global boundary-fixing isometry, while the extension hypotheses fail; or, at Lamé parameters satisfying $(n-1)\\lambda^3+(4n-2)\\lambda^2\\mu+(n+5)\\lambda\\mu^2+(14-8n)\\mu^3=0$, exhibit two real-analytic metrics whose DtN symbols agree but whose first normal derivatives at the boundary differ.","supporting_citations":[{"cited_title":"Lee and G","cited_arxiv_id":null,"evidence_quote":"provides the scalar boundary-determination prototype and the boundary-to-global extension lemma used to complete the isometric uniqueness proof."},{"cited_title":"Duduchava, D","cited_arxiv_id":null,"evidence_quote":"supplies the integration-by-parts identity and elastic-energy variational formulation used to derive the Lamé operator on a Riemannian manifold."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the heat-trace expansion method and the sphere integrals used to evaluate the spectral coefficients."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the Sylvester-equation theory and Kronecker-product formulation used to solve for the lower-order symbols of Q."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the pseudodifferential symbol calculus and contour-resolvent formalism used in both the symbol recovery and spectral trace computations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"underlies the geometric extension principle that turns a collar isometry into a global boundary-fixing diffeomorphism."},{"cited_title":"Treves, Introduction to pseudodifferential and Fourier integral op erator, Plenum Press, New York, 1980","cited_arxiv_id":null,"evidence_quote":"supplies the classical factorization method for scalar boundary value problems that the elastic factorization adapts to the matrix Lamé system."}],"review_version":1}