{"id":"18bcf595-f972-44e1-b2d4-12170b9fcdc6","arxiv_id":"2607.06890","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A pure-pseudostress mixed FEM for the elasticity eigenproblem is locking-free, needs no symmetry constraint, and comes with a priori and residual a-posteriori error estimates that remain valid as Poisson ratio approaches 1/2.","lead":"The paper gives a locking-free mixed finite-element scheme for the linear elasticity eigenvalue problem that uses only the pseudostress tensor, without enforcing symmetry. Displacement is recovered by post-processing, and both a priori and residual a-posteriori estimates are proved and tested in 2-D and 3-D, including the incompressible limit.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified beyond the already-flagged Assumption 2.6","rationale":"The manuscript’s central claim—that the pure-pseudostress RT method is locking-free for the elasticity eigenproblem—rests on the spectral convergence theory for non-compact operators and on the residual estimator remaining robust as λ → ∞. Both pieces are rigorously derived once the eigenfunction regularity (Lemma 2.4) is independent of λ. The authors themselves flag that this independence “is not completely evident” and adopt it as Assumption 2.6 solely on the basis of numerical observation. That is precisely the weakest link the reader identified; no stronger or more fundamental gap (e.g., failure of the discrete inf-sup, loss of self-adjointness, or λ-dependent constants inside the a-posteriori analysis) is present. Consequently the CONDITIONAL verdict with high confidence is already the correct one, and no adjustment is required. The concrete test simply quantifies how far the numerical evidence can be pushed before the assumption would have to be revisited.","tokens_in":25526,"tokens_out":545,"duration_ms":5852,"concrete_test":"Re-run the unit-square experiment of Table 1 with a sequence of Poisson ratios ν = 0.499, 0.4999, 0.49999 (i.e., λ → ∞) on the same family of uniform meshes and extract the observed order of |κ − κ_h|; if the order remains 2 (within fitting tolerance) down to the largest λ, the numerical support for Assumption 2.6 is strengthened; any systematic drop would confirm that the constants deteriorate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader's identification of Assumption 2.6 is accurate and is the only soft spot that could make the locking-free rates λ-dependent. The remainder of the argument is internally consistent: the pure-pseudostress formulation (2.5)–(2.10), the non-compact spectral theory via Descloux–Nassif–Rappaz (Lemmas 3.1, 4.1–4.3, Theorems 4.4–4.5), the limit operator T_∞ and spectrum separation (Lemma 2.8, Theorem 2.9), and the residual estimator with λ-independent weights (Theorems 5.1, 5.4) all follow standard mixed-FEM reasoning once the regularity constants are granted. No hidden inconsistency, circularity, or missing estimate appears that would undermine the central claim more severely than the already-noted assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper proposes a pure-pseudostress mixed variational formulation of the linear elasticity eigenvalue problem (2.5)–(2.10) that avoids any symmetry constraint on the tensor. The associated solution operator T_λ is non-compact; spectral characterization (Theorem 2.3), convergence of the nearly-incompressible spectrum to the Stokes spectrum (Lemma 2.8, Theorem 2.9), and discrete analysis via Raviart–Thomas elements are carried out with the Descloux–Nassif–Rappaz theory. Optimal a-priori rates for eigenfunctions and double-order rates for eigenvalues (Theorems 4.4–4.5) are proved under an explicit regularity assumption independent of λ. A residual a-posteriori estimator η (and its incompressible counterpart) is shown reliable and efficient with λ-independent weights (Theorems 5.1, 5.4). Numerical experiments on convex and non-convex domains in 2D/3D confirm locking-free convergence and adaptive recovery of optimal rates.","tokens_in":25732,"tokens_out":936,"duration_ms":9341,"significance":"A locking-free pure-pseudostress eigenvalue formulation that never enforces symmetry is a useful addition to the mixed-FEM literature for elasticity eigenproblems. The analysis is self-contained once standard Sobolev regularity is granted, the residual estimator is new for this setting, and the numerical tests (including adaptive refinement near re-entrant corners and the incompressible limit) are reproducible and support the claims. If Assumption 2.6 holds, the method supplies a practical, symmetry-free alternative to existing mixed schemes that remains robust as ν→1/2.","major_comments":[{"comment":"Assumption 2.6 (independence of the regularity exponent r and constant Ĉ of Lemma 2.4 with respect to λ) is load-bearing for the locking-free rates of Theorems 4.4–4.5 and for the claim that the estimator constants are independent of λ. The manuscript itself notes that the independence “is not completely evident” and is adopted only because numerics still show the expected orders at λ=∞. Either a reference establishing λ-uniform regularity for the pure-pseudostress source problem, or a short remark quantifying possible deterioration of the constants, is needed before the locking-free claim can be regarded as fully rigorous.","section":null},{"comment":"Reliability (Theorem 5.1) is stated with an extra term ∥κρ−κ_h ρ_h∥_0 + ∥ρ−ρ_h∥_0 on the right-hand side; efficiency (Theorem 5.4) absorbs a generic higher-order remainder Θ. For the estimator to be fully practical one needs either a proof that these terms are of higher order (or controlled by η itself) or a clear statement that they are neglected only after the discrete eigenpair has already converged. The present wording leaves a small gap between the proved bounds and the quantity that is actually used for marking.","section":null}],"minor_comments":[{"comment":"Several typographical slips appear: “the discrete the eigenvalue problem” (heading 3.2), “let su assume” (p. 16), “con the perfectly incompressible case” (abstract), and inconsistent boldface for tensor/vector fields.","section":null},{"comment":"The definition of the weights ρ_R and ρ_E (and their incompressible analogues) is dense; a short table or displayed list would improve readability of Section 5.1.","section":null},{"comment":"In the numerical section the effectivity index is defined as err(κ_i)/η²; a brief justification why the square appears (consistent with the double-order eigenvalue estimate) would help the reader.","section":null},{"comment":"References [13,15,20,22] are heavily used for background lemmas; a one-sentence pointer to the precise statements that are imported would make the paper more self-contained for non-specialists.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The paper is a solid, incremental contribution that fits a numerical-analysis journal. The only soft spot is the explicit regularity assumption; once the authors address it (even by a short literature citation or a caveat), the manuscript is ready for acceptance. No concerns about novelty disclosure or citation pattern."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This is a clean, complete analysis of a pure-pseudostress mixed formulation for the linear elasticity eigenvalue problem. The only unknown is the pseudostress; displacement is recovered by post-processing, and symmetry is never enforced. That is the real novelty relative to the Stokes pseudostress work and Meddahi’s DG stress paper: a conforming RT discretization that stays locking-free all the way to the incompressible limit, with both a-priori rates under the Descloux–Nassif–Rappaz theory and a residual estimator that is reliable and efficient independently of λ.\n\nThe continuous spectral theory is carefully done: the solution operator is non-compact, they characterize the spectrum (including the infinite-multiplicity eigenvalue 1 and the kernel of the deviatoric part), prove that T_λ converges to the Stokes operator T_∞ as λ→∞, and obtain the usual gap and double-order eigenvalue estimates once the discrete operators satisfy the standard P1/P2 conditions. The residual indicator is new for this pure-stress setting and the bubble-function efficiency argument is standard and clean. Numerics on squares, re-entrant sectors and a notched cube confirm the rates and the adaptive recovery of optimal dof-decay even at ν=1/2.\n\nThe single soft spot is Assumption 2.6: the Sobolev regularity exponent and constant for the eigenfunctions are taken independent of λ. The authors themselves say this “is not completely evident” and adopt it because the experiments still give the expected orders at λ=∞. If the constants deteriorate, the locking-free claim becomes λ-dependent. Everything else is self-contained once that is granted; there is no circularity and the citation trail is appropriate.\n\nThis is for people who work on mixed methods for elasticity or Stokes eigenvalues and want a locking-free pure-stress scheme with a posteriori control. It deserves a serious referee. I would send it out.","headline":"Solid pure-pseudostress mixed FEM for the elasticity eigenproblem; locking-free rates rest on one flagged regularity assumption.","tokens_in":26348,"tokens_out":481,"would_cite":true,"duration_ms":5603,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35P15","65N15","65N25","65N30","74B05"],"pacs":[],"model":"grok-4.5","headline":"A pure-pseudostress mixed finite-element method for elasticity eigenvalues is locking-free for any Poisson ratio, recovers displacement by post-processing, and comes with a residual estimator that is reliable and efficient independently of","keywords":["elasticity eigenvalue problem","pseudostress formulation","locking-free mixed FEM","Raviart–Thomas elements","non-compact operators","a-posteriori error estimation","nearly incompressible materials"],"falsifier":"Compute the first few eigenvalues on a sequence of uniformly refined meshes for Poisson ratios successively closer to 1/2 (or for the formal limit λ=∞) and check whether the observed orders of convergence for both eigenvalues and eigenfunctions remain exactly those predicted by the a-priori theory; any systematic degradation would falsify the locking-free claim.","tokens_in":26388,"feed_emoji":"⚙️","tokens_out":1019,"duration_ms":11533,"temperature":0.7,"pith_summary":"The paper introduces a mixed variational formulation of the linear-elasticity eigenvalue problem whose only unknown is the pseudostress tensor. Because the formulation never enforces symmetry of the stress and never involves the displacement as a primary variable, standard Raviart–Thomas elements can be used without locking when the material becomes nearly or perfectly incompressible. The continuous solution operator is non-compact, so spectral convergence and optimal a-priori rates for eigenvalues and eigenfunctions are proved by the theory of non-compact operators. As the Lamé parameter tends to infinity the spectrum converges to that of the Stokes eigenvalue problem, recovering the incompressible limit. A residual-based a-posteriori indicator is shown to be both reliable and efficient uniformly in the compressibility parameter; adaptive meshes driven by this indicator restore optimal rates on non-convex domains in two and three dimensions. Displacement and true stress are recovered by elementary post-processing of the computed pseudostress.","feed_headline":"Locking-free elasticity eigenvalues from pure pseudostress","feed_subtitle":"Raviart–Thomas elements stay accurate for any Poisson ratio; residual estimator drives optimal adaptive meshes","key_machinery":"The self-adjoint non-compact solution operator T_λ that maps a load to the unique pseudostress solving the shifted source problem; its spectral decomposition and the discrete counterpart T_λ,h allow the application of the Descloux–Nassif–Rappaz theory to obtain locking-free a-priori estimates.","core_discovery":"A pure-pseudostress formulation of the elasticity eigenproblem, discretized by tensor-valued Raviart–Thomas elements, is locking-free: eigenvalues and eigenfunctions converge at the rates predicted by the regularity of the eigenfunctions, the spectrum of the nearly-incompressible operator converges to the Stokes spectrum, and a residual estimator remains reliable and efficient independently of the Lamé parameter.","pith_inferences":["The same pure-pseudostress idea should extend, with only minor changes, to related non-self-adjoint or damped eigenvalue problems in viscoelasticity.","Because the formulation never requires symmetry, it is a natural candidate for hybridization or static condensation that would further reduce the algebraic cost of three-dimensional computations.","The residual estimator’s independence of λ suggests it could serve as a reliable stopping criterion in iterative solvers that themselves become ill-conditioned near the incompressible limit."],"forward_implications":["Standard H(div)-conforming elements can be used for vibration analysis of nearly-incompressible elastic bodies without artificial stiffening.","The spectrum of the discrete elasticity operator converges to the Stokes spectrum as the Poisson ratio approaches 1/2, giving a practical route to incompressible eigencomputations.","The residual estimator can drive adaptive refinement that recovers optimal rates on domains with re-entrant corners or edges, uniformly in compressibility.","Displacement and true Cauchy stress are obtained by simple post-processing of the pseudostress, so no additional mixed system needs to be solved."],"fun_headline_variants":["Locking-free elasticity eigenpairs from pure pseudostress FEM","Tensor RT elements solve elasticity eigenvalues without locking","Pure-pseudostress mixed method recovers displacement via postprocess","Residual estimator stays reliable for any Lamé parameter","Non-compact theory guarantees locking-free spectrum convergence"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"The Sobolev regularity of the eigenfunctions and the constant that bounds it are assumed independent of the Lamé parameter, even though the paper itself notes that this independence is not completely evident from the analysis.","fun_headline_variants_meta":{"raw":{"variants":["Locking-free elasticity eigenpairs from pure pseudostress FEM","Tensor RT elements solve elasticity eigenvalues without locking","Pure-pseudostress mixed method recovers displacement via postprocess","Residual estimator stays reliable for any Lamé parameter","Non-compact theory guarantees locking-free spectrum convergence"]},"model":"grok-4.5","effort":"low","cost_usd":0.003658,"raw_usage":{"total_tokens":1073,"prompt_tokens":695,"num_sources_used":0,"completion_tokens":80,"cost_in_usd_ticks":36580000,"prompt_tokens_details":{"text_tokens":695,"audio_tokens":0,"image_tokens":0,"cached_tokens":0},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":298,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":695,"tokens_out":80,"duration_ms":4320,"temperature":1.0,"reasoning_tokens":298,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-10T20:07:14.350142+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Compute the first few eigenvalues on a sequence of uniformly refined meshes for Poisson ratios successively closer to 1/2 (or for the formal limit λ=∞) and check whether the observed orders of convergence for both eigenvalues and eigenfunctions remain exactly those predicted by the a-priori theory; any systematic degradation would falsify the locking-free claim.","supporting_citations":[],"review_version":2}