{"id":"79c1e4f7-2e70-4206-ad34-a06a6bd85863","arxiv_id":"2501.04743","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The paper presents several equivalent first-order symmetric forms of linear elastodynamics in 1D and 2D by adding compatibility, momentum, and time-differentiated constitutive equations, with a correctness issue in the last 2D set.","lead":"This paper rewrites the classical equations of linear elasticity as first-order systems with symmetric coefficient matrices, using several choices of variables (velocity, strain or displacement gradient, stress, momentum). The new forms are mostly equivalent combinations of known formulations and could be useful for numerical analysis, but the final two-dimensional version contains matrix errors as printed.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"§3.4's four-variable 2D matrices are not symmetric and do not match the equations from which they are assembled, so the most complex advertised symmetric writing is unsupported as printed.","rationale":"The reader's explicit weakest_assumption was the time-differentiated constitutive constraint, which is a real but minor omission (the constraint is preserved for admissible data). However, the reader's rationale already identifies the decisive issue as the §3.4 matrix asymmetry. My stress-test confirms that issue concretely: the printed B1 has asymmetric entries and mismatches with the very equations it is built from, so the claim for the (ρ0ut, ut, e, σ) system is not supported as submitted. This is not a conceptual failure; the 1D systems and simpler 2D systems appear correct, and the correction of the earlier [8] matrix issue is a useful contribution. The problems are typographical and localized, so a revised version with corrected matrices and a symbolic verification script would make the claims acceptable. Since the reader's verdict was already CONDITIONAL and my concern reinforces it rather than changing it, the verdict stays UNCHANGED.","tokens_in":22431,"tokens_out":14019,"duration_ms":97076,"concrete_test":"Use a computer algebra system (e.g., SymPy) to assemble the system from the ten printed equations (155)–(182) by coefficient matching into A q_t + B1 q_x1 + B2 q_x2 = 0 for a generic anisotropic stiffness tensor with C1122 ≠ C1222 (for instance C1111=5, C1122=2, C1222=3, C2111=4, C1112=0.7, C2122=1.1). Then (i) check whether the extracted B1 and B2 are symmetric and (ii) compare them with the printed eqs. (185)–(186). If B1 is not symmetric or differs from the printed matrix, the §3.4 claim fails as written; a minimal correction would also need to be re-verified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of the paper is that each listed variable set admits a legitimate symmetric first-order system. The most complex 2D writing, §3.4 (q = (ρ0ut, ut, e, σ)), fails this test as printed. The displayed matrix B1 (eq. 185) is not symmetric: entry (1,6) is −C1122 while entry (6,1) is −C1222, and these are not equal under the stated minor/major symmetries. Several entries also do not reproduce the equations (155)–(182) they are claimed to come from: (4,7) and (7,4) are −2C2122, but momentum equations (164) and (173) require −2C2112; (5,6) and (6,5) are −C1222, but (167) and a correct reading of (170) require −C1122. Equation (170) itself appears to contain typos, printing −C1222 v1,1 and −C1222 v5,1 where the derivation from the §3.3 compatibility equation (138) requires −C1122. Thus the advertised symmetric system for (ρ0ut, ut, e, σ) is not established. The constraint-preservation issue with time-differentiated constitutive laws is a smaller omission: for admissible initial data satisfying σ = C:e, the constraint is preserved by a linear ODE (e.g., in 1D wt = (αt/α)w, so w=0 persists). The decisive, checkable problem is the §3.4 matrix assembly.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an inverse construction of first-order symmetric systems for classical linear elastodynamics. In 1D it presents symmetric writings for seven variable sets, including (u_t,u_x), (u_t,σ), (u_t,u_x,σ), momentum-based versions, and a four-variable set (ρ0u_t,u_t,u_x,σ); in 2D it presents four writings, the most complex being (ρ0u_t,u_t,e,σ). The construction forms linear combinations of the momentum equation, compatibility equations, and the time-differentiated constitutive law, then verifies symmetry of the resulting matrices A and B_i. The paper also contains a remark correcting a mistake in the author's earlier work [8].","tokens_in":22675,"tokens_out":29458,"duration_ms":207748,"significance":"If the construction is correct, the paper offers a useful catalog of symmetric first-order formulations of anisotropic linear elastodynamics, with potential applications to finite element methods, Hamiltonian formulations, and semigroup theory. The approach is fully algebraic and parameter-free, and the 1D constructions together with the simpler 2D systems are largely plausible. The main weakness is that the most complex 2D writing, §3.4, is not correct as printed, due to numerous inconsistencies between the displayed equations and the assembled matrices; these appear to be repairable typographical errors rather than a fundamental flaw, but they must be fixed before the central claim can be accepted.","major_comments":[{"comment":"Equation (158) is missing the factor 2 multiplying v2,t. Adding Eqs. (156) and (157) yields 2v2,t − ... = 0, and the (2,2) entry of A in Eq. (184) is 2. As printed, the equation does not match the assembled system.","section":"§3.4, Eq. (158)"},{"comment":"The coefficient of v1,1 in Eq. (168) and of v5,1 in Eq. (169) is printed as −C1222, but the compatibility equation (138) from §3.3 requires −C1122. The error is carried into the sum (170) and then into the B1 matrix.","section":"§3.4, Eqs. (168)–(170)"},{"comment":"Matrix B1 is not symmetric as printed and does not reproduce the stated field equations. Row 5, column 6 should be −C1122 (from Eq. (167)) but is printed as −C1222; row 6, columns 1 and 5 should be −C1122 (from the corrected Eq. (170)) but are printed as −C1222; row 3, column 7 should be −C2111 (from Eq. (161)) but is printed as −C1211; row 4, column 7 should be −2C2112 (from Eq. (164)) but is printed as −2C2122. With these corrections B1 becomes symmetric.","section":"§3.4, Eq. (185)"},{"comment":"The terms multiplying v7,t and v8,t are interchanged. The equation should read C−1_2211 v6,t + C−1_2222 v7,t + 2C−1_2212 v8,t − v5,2 = 0 to agree with the A matrix (144) and with the structure of the time-differentiated constitutive law.","section":"§3.3, Eq. (141)"},{"comment":"The terms C−1_1122 v4,2, C−1_2222 v4,2, and C−1_1222 v4,2 must be time derivatives v4,t. As printed, these equations are inconsistent with the A matrix (131) and with the fact that v4 = σ22.","section":"§3.2, Eqs. (127)–(129)"},{"comment":"The systems that use the time-differentiated constitutive law as a field equation are equivalent to classical elastodynamics only for initial data satisfying σ = C:e (or σ = α ux in 1D). The paper does not state or prove that this constraint is preserved. For time-independent coefficients preservation is immediate (wt = 0); for the 1D case with α = α(t) one obtains wt = (αt/α)w. This gap should be closed explicitly. Relatedly, Eq. (119) in §3.2 omits the (C−1)t σ term that appears in the 1D treatment (Eqs. (15)–(16)); the paper should state whether C is assumed time-independent.","section":"§2.2–§3.4, constraint preservation"}],"minor_comments":[{"comment":"The last term should be −2C2212 v5,2, not −C2212 v5,2, to match the A and B matrices.","section":"§3.3, Eq. (136)"},{"comment":"The expressions 'C11122' and 'C12122' in Eqs. (153)–(158), (167), and (173) should be written as '2C1112' and '2C1212' to avoid ambiguity.","section":"§3.4, notation"},{"comment":"Equation (178) is garbled: '2C−1_12122ρ0v10,t' should be '2ρ0 C−1_2212 v10,t' (or equivalently '2ρ0 C−1_2122 v10,t') to match Eq. (177) multiplied by ρ0.","section":"§3.4, Eq. (178)"},{"comment":"The paper verifies only that det A is nonzero, but for a symmetric hyperbolic system in the sense of Friedrichs the matrix A must be positive definite. The relevant blocks are Hessians of the (positive definite) elastic energy and complementary energy, so positive definiteness holds, but this should be stated explicitly.","section":"§3.1 and §3.3"},{"comment":"The time-differentiated constitutive law is written as C−1 σ̇ = ė, which silently assumes time-independent C; if time-dependent coefficients are allowed, as the 1D treatment does, the term (C−1)t σ must be included.","section":"§3.2, Eq. (119)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript would benefit from a careful proofreading pass; the number of typographical errors in the displayed matrices is high, and the most complex claimed construction is not reproducible as printed. The author's correction of a mistake in their own earlier work [8] is a positive sign of scholarly care. The errors in §3.4 appear correctable, and the core method is sound, so major revision rather than rejection is appropriate."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful part of this paper is the collection of equivalent symmetric first-order writings for linear elastodynamics, and the concrete correction it makes to the author's earlier work: the positive semi-definiteness issue in [8] is real, and the 5-equation (u_t,e) system here fixes it. The 1D derivations are clean and the simpler 2D systems—(u_t,e), (u_t,σ), and (u_t,e,σ)—are assembled correctly, at least as far as I can tell by hand. The three-variable 2D system (eqs. 149–151) is genuinely symmetric and consistent with the equations it comes from. The momentum-based writings are formal but legitimate substitutions; no free parameters, no fitting, and the self-citation to [8] is used as a starting point rather than a crutch.\n\nThe soft spot is §3.4, and it is not minor. The advertised four-variable 2D system (ρ0u_t, u_t, e, σ) fails as printed. The displayed B1 in eq. (185) is not symmetric: entry (1,6) is −C1122 while entry (6,1) is −C1222, which are not equal under the stated symmetries. Entry (4,7) is −2C2122 when the momentum equations (164) and (173) require −2C2112. Entry (5,6) is −C1222 when (167) and a consistent reading of (170) require −C1122. Equation (170) itself contains typos, printing −C1222 where the compatibility equation (138) demands −C1122. These are specific and checkable, and they mean the central claim for that system is unsupported as submitted. I think these are typographical errors rather than a conceptual flaw—the assembly pattern is the same as in the three-variable system—but the burden is on the author to fix them.\n\nOne smaller omission: the time-differentiated constitutive law means each system is equivalent to classical elastodynamics only on the constraint manifold σ = C:e. The paper never says this. It is a minor issue because in the linear case the constraint is preserved by the evolution (e.g., in 1D, w_t = (α_t/α)w), so admissible initial data stay admissible. Still, a sentence or two would clear it up.\n\nWho is this for? People working on symmetric hyperbolic formulations of elasticity, especially anyone classifying implicit theories by hyperbolicity or looking for first-order systems for numerics. It deserves a serious referee: the idea is sound, the correction of [8] is useful, and the problems in §3.4 are local and fixable. I would not cite it in its current state, but I would happily revisit a revised version with corrected matrices and, ideally, a symbolic verification script. Send it to review, with the request that the author fix §3.4 and state the constraint condition.","headline":"The 1D systems and the simpler 2D systems check out, but the four-variable 2D matrices in §3.4 have concrete symmetry and assembly errors, so the paper is a conditionally useful working draft rather than a finished result.","tokens_in":903,"tokens_out":856,"would_cite":false,"duration_ms":25931,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["74B05","74B99","35L02","35Q74"],"pacs":[],"model":"deepseek-v4-flash","headline":"Classical elastodynamics can be written as a first-order symmetric system for eleven different choices of variables.","keywords":["Linear theory","First order symmetric system","Elastodynamics","Anisotropic elasticity","Compatibility equations","Time-differentiated constitutive law","Symmetric hyperbolic system","Momentum-velocity conjugacy"],"falsifier":"Solve the 1D three-variable system of eqs. (30)-(32) with constant $\\alpha$ and initial data $u_t(0,x)=0$, $u_x(0,x)=0$, $\\sigma(0,x)=1$. If the solution develops $\\sigma-\\alpha u_x\\neq 0$ while satisfying the system, the symmetric writing admits non-elastodynamic motions and full equivalence fails; if the constraint is preserved, equivalence is confirmed.","tokens_in":22137,"feed_emoji":"📐","tokens_out":7861,"duration_ms":70987,"temperature":0.7,"pith_summary":"The paper tries to establish that the equations of classical linear anisotropic elastodynamics, in one and two space dimensions, admit first-order symmetric writings of the form $A \\frac{\\partial q}{\\partial t} + \\sum_i B_i \\frac{\\partial q}{\\partial x_i}=0$ for a menu of variable sets: velocity with displacement gradient or strain, velocity with stress, all three together, and variants in which momentum $\\rho_0 u_t$ replaces or accompanies velocity. The construction is inverse: the author writes the desired system, forces the spatial-derivative matrices $B_i$ to be symmetric by combining compatibility equations, the momentum equation, and the time-differentiated constitutive law, and then checks that the time-derivative matrix $A$ is symmetric and positive definite. If the central claim is right, every listed writing is a legitimate symmetric hyperbolic formulation of classical elastodynamics, and that matters because symmetric form is the standard route to energy estimates, well-posedness arguments, and finite-element treatments of wave propagation.","feed_headline":"Elastodynamics takes symmetric first-order form in 11 variable sets","feed_subtitle":"New writings unify compatibility, momentum, and differentiated material laws for anisotropic wave equations.","key_machinery":"The machinery is an inverse symmetrization procedure. Starting from the first-order form, the author fixes a variable vector $q$ and writes candidate field equations drawn from three sources: the momentum equation, the compatibility equations linking strain components to velocity gradients, and the time-differentiated constitutive law $\\dot\\sigma = C:\\dot e$ (written through the compliance $C^{-1}$ for stress-based systems). These equations are multiplied and added so that every spatial-derivative matrix $B_i$ becomes symmetric; the same additions alter the time-derivative matrix $A$, and the writing is accepted when $A$ comes out symmetric and, for hyperbolicity, positive definite. The compatibility equations carry the strain- or displacement-gradient-based writings, while the time-differentiated constitutive law carries the stress-based writings.","core_discovery":"The paper's central claim is a catalog: for the one-dimensional case it produces symmetric writings with respect to $(u_t,u_x)$, $(u_t,\\sigma)$, $(u_t,u_x,\\sigma)$, $(\\rho_0 u_t,u_x)$, $(\\rho_0 u_t,\\sigma)$, $(\\rho_0 u_t,u_x,\\sigma)$, and $(\\rho_0 u_t,u_t,u_x,\\sigma)$; for the two-dimensional case it produces symmetric writings with respect to $(u_t,\\mathbf e)$, $(u_t,\\sigma)$, $(u_t,\\mathbf e,\\sigma)$, and $(\\rho_0 u_t,u_t,\\mathbf e,\\sigma)$. In each writing the matrices $A$ and $B_i$ are explicitly symmetric, and $A$ is positive definite when the material parameters satisfy the stated conditions ($\\alpha>0$ in 1D; the usual major and minor symmetries plus the displayed invertibility and determinant conditions in 2D). The paper also reports a correction to an earlier two-dimensional $(u_t,\\mathbf e)$ writing: the earlier version had dependent compatibility equations and only positive semidefinite $A$, while the new five-equation version has independent equations and an invertible $A$.","pith_inferences":["Because the stress- and combined-variable systems use the time-differentiated constitutive law, a reader should test numerically whether the algebraic constraint $\\sigma=C:\\mathbf e$ is preserved along solutions; if it drifts, the equivalence to classical elastodynamics holds only on the constrained initial-data submanifold.","The same inverse symmetrization recipe should extend to three dimensions and to nonlinear or implicit constitutive laws, at the cost of heavier bookkeeping; the paper stops at 2D and describes 3D as a straightforward generalization.","The four-variable writings with both momentum and velocity are deliberately redundant, so they are better read as templates for variational or port-Hamiltonian formulations than as minimal evolution systems."],"forward_implications":["Each of the eleven writings is a symmetric first-order system, so Friedrichs' theory of symmetric hyperbolic systems applies and energy estimates can be obtained in every variable set.","The combined writings that use velocity, strain or displacement gradient, and stress tie the compatibility equations and the differentiated constitutive law into one system, which the paper suggests as a setting for studying hyperbolicity of implicit constitutive theories.","The momentum-based writings make the velocity-momentum conjugacy explicit, connecting the symmetric forms directly to kinetic energy and to candidate Hamiltonian structures.","The corrected two-dimensional $(u_t,\\mathbf e)$ writing has five independent equations with an invertible symmetric $A$, replacing the earlier six-equation version whose redundant compatibility equations left $A$ only positive semidefinite."],"supporting_citations":[{"why":"Defines the symmetric hyperbolic systems that are the target format for all writings in the paper.","marker":"[2]"},{"why":"Supplies the standard stress-velocity symmetric writing built from the time-differentiated constitutive law, the baseline extended here.","marker":"[3]"},{"why":"Introduces the inverse symmetrization method and the compatibility-equation route for writing elastodynamics with $(u_x,u_t)$.","marker":"[8]"},{"why":"Another classical source for using the time-differentiated constitutive law to obtain first-order symmetric systems.","marker":"[10]"},{"why":"Gives the anisotropic initial-value formulation in stress and velocity that the paper follows for the 2D $(u_t,\\sigma)$ writing.","marker":"[12]"},{"why":"Provides the implicit-constitutive-theory hyperbolicity question that motivates the combined writings with all three variables.","marker":"[9]"}],"fun_headline_variants":["Elastodynamics: 11 symmetric first-order forms","Symmetric rewritings of elastodynamic equations for 11 variable sets","Corrected and extended symmetric first-order elastodynamics","Unified symmetric system for elastodynamics from momentum to stress","Elastodynamic equations: symmetric first-order versions for all variable choices"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The rewrite replaces the material law (stress equals elasticity times strain) by its time-derivative form, so the new systems agree with classical elastodynamics only for initial data that already satisfy the material law, and the paper does not prove that the evolution keeps those data on it.","fun_headline_variants_meta":{"raw":{"variants":["Elastodynamics: 11 symmetric first-order forms","Symmetric rewritings of elastodynamic equations for 11 variable sets","Corrected and extended symmetric first-order elastodynamics","Unified symmetric system for elastodynamics from momentum to stress","Elastodynamic equations: symmetric first-order versions for all variable choices"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001138,"raw_usage":{"total_tokens":4851,"prompt_tokens":1194,"completion_tokens":3657,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":810,"completion_tokens_details":{"reasoning_tokens":3574}},"tokens_in":810,"tokens_out":3657,"duration_ms":26343,"temperature":1.0,"reasoning_tokens":3574,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:37:19.780241+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the 1D three-variable system of eqs. (30)-(32) with constant $\\alpha$ and initial data $u_t(0,x)=0$, $u_x(0,x)=0$, $\\sigma(0,x)=1$. If the solution develops $\\sigma-\\alpha u_x\\neq 0$ while satisfying the system, the symmetric writing admits non-elastodynamic motions and full equivalence fails; if the constraint is preserved, equivalence is confirmed.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the symmetric hyperbolic systems that are the target format for all writings in the paper."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the standard stress-velocity symmetric writing built from the time-differentiated constitutive law, the baseline extended here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the inverse symmetrization method and the compatibility-equation route for writing elastodynamics with $(u_x,u_t)$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Another classical source for using the time-differentiated constitutive law to obtain first-order symmetric systems."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the anisotropic initial-value formulation in stress and velocity that the paper follows for the 2D $(u_t,\\sigma)$ writing."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the implicit-constitutive-theory hyperbolicity question that motivates the combined writings with all three variables."}],"review_version":1}