arxiv: 2604.03727 · v1 · submitted 2026-04-04 · 🧮 math.NA · cs.NA

Recognition: 2 theorem links

· Lean Theorem

A high order stabilization-free virtual element method for general second-order elliptic eigenvalue problem

Hai Bi, Liangkun Xu, Shixi Wang, Yidu Yang

Authors on Pith no claims yet

Pith reviewed 2026-05-13 16:52 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords virtual element methodstabilization-freeelliptic eigenvalue problempolygonal meshesa priori error estimateshigh-order methodnumerical analysisfinite element discretization

0 comments

The pith

A stabilization-free higher-order virtual element method yields optimal error estimates for elliptic eigenvalue problems on polygonal meshes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a high-order virtual element method for general second-order elliptic eigenvalue problems that requires no stabilization term. It establishes optimal a priori error estimates for the approximate eigenspaces and eigenvalues. The approach is tested through numerical experiments on regular convex, convex-concave, and concave polygonal meshes. A reader would care because the method offers a simpler discretization on unstructured meshes while retaining rigorous accuracy guarantees for eigenvalue computations.

Core claim

The authors construct a high-order stabilization-free virtual element method for general second-order elliptic eigenvalue problems. They prove optimal a priori error estimates for both the approximate eigenspace and the eigenvalues. Stability and consistency are obtained directly from the virtual element spaces and projection operators on polygonal meshes without any added stabilization term. The theoretical results are supported by numerical experiments across multiple mesh families.

What carries the argument

Stabilization-free virtual element spaces of high order equipped with projection operators that supply the necessary stability and consistency for the eigenvalue discretization on polygonal meshes.

If this is right

Optimal a priori error estimates hold for both approximate eigenspaces and eigenvalues.
The method applies directly to general second-order elliptic eigenvalue problems on polygonal meshes.
Higher-order accuracy is achieved through the virtual element construction without stabilization.
Effectiveness is confirmed numerically on convex, convex-concave, and concave polygonal meshes.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The removal of stabilization may simplify parameter selection and implementation for eigenvalue problems on complex geometries.
The same projection-based stability could extend the method to related source or time-dependent problems with similar mesh flexibility.
Performance under severe mesh distortion or in three-dimensional domains remains open for further verification.

Load-bearing premise

The virtual element spaces and projection operators are assumed to maintain stability and accuracy for the general second-order elliptic eigenvalue problem on polygonal meshes without any stabilization term.

What would settle it

Numerical tests on successively refined concave polygonal meshes in which the computed eigenvalues or eigenspaces fail to converge at the predicted optimal rates would falsify the error estimates.

read the original abstract

In this paper, we discuss a novel higher-order stabilization-free virtual element method for general second-order elliptic eigenvalue problems. Optimal a priori error estimates are derived for both the approximate eigenspace and eigenvalues. Numerical experiments are conducted on regular convex polygonal meshes, convex-concave polygonal meshes, and concave polygonal meshes. The numerical results validate the effectiveness of the proposed method.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives a stabilization-free higher-order VEM for elliptic eigenvalue problems with claimed optimal rates, but the coercivity argument on concave meshes stays thin.

read the letter

The paper introduces a higher-order virtual element method for general second-order elliptic eigenvalue problems that skips the usual stabilization term. It relies on projection operators to build the discrete bilinear form and claims optimal a priori error estimates for both eigenvalues and eigenspaces. Numerical tests are run on convex, mixed convex-concave, and concave polygonal meshes, and the results line up with the predicted rates on those examples. That combination of higher order, no stabilization, and eigenvalue focus is the concrete new piece; most prior VEM work on eigenvalues kept a stabilization term or stayed at low order. The method is laid out clearly enough for someone already comfortable with VEM spaces and degrees of freedom to follow the construction and the consistency argument. The numerical section gives a reasonable first check across mesh types, which is better than many method papers that only show convex cases. The soft spot is the stability claim. The estimates rest on the discrete form being coercive and continuous without stabilization, yet the argument does not supply an explicit lower bound that stays uniform when polygons become more concave or distorted. The virtual space on concave elements can admit functions whose gradients are not controlled by the L2 projection alone, so the Poincaré-type control has to come from the higher moments or the specific degrees of freedom. No separate proof or numerical sweep that tracks the coercivity constant against increasing concavity is given, which leaves the uniform bound unverified. The citation list follows the standard VEM references for source problems and eigenvalue extensions, with no obvious gaps or circular self-citation. This is a paper for people already working in virtual element or polygonal finite element methods for eigenvalues. A reader outside that niche will find the technical setup heavy, but specialists will see the incremental extension and can judge whether the missing uniform coercivity step is a minor gap or a load-bearing one. I would send it to referees. The claims are concrete enough that a careful review can check the stability proof and ask for either a uniform bound or a clearer statement of the mesh assumptions.

Referee Report

2 major / 1 minor

Summary. The paper proposes a novel high-order stabilization-free virtual element method for general second-order elliptic eigenvalue problems on polygonal meshes. It claims to derive optimal a priori error estimates for both the approximate eigenspaces and the eigenvalues, and validates the approach via numerical experiments on regular convex polygonal meshes, convex-concave meshes, and concave polygonal meshes.

Significance. If the stability analysis and error estimates are rigorous, the result would be significant for the VEM literature: it removes the stabilization term while retaining optimal convergence for eigenvalue problems, which could simplify high-order implementations on general (including non-convex) meshes and reduce parameter tuning.

major comments (2)

[stability and consistency analysis] The section deriving the discrete bilinear form and its coercivity: the claim that the consistency term alone (via the projection operators) yields a coercive and continuous form without stabilization is load-bearing for the optimal a priori estimates, yet no explicit lower bound independent of the mesh distortion parameter is shown for concave polygons; the argument implicitly assumes a discrete Poincaré inequality from the degrees of freedom, but this is not proven uniformly.
[numerical experiments] The numerical experiments section (concave-mesh tests): the reported results on concave polygons do not include a systematic check of the discrete coercivity constant (or its dependence on concavity/aspect ratio), so the robustness claim for the stabilization-free form remains unsupported by quantitative evidence.

minor comments (1)

[abstract] The abstract should state the precise polynomial degree and the form of the elliptic operator (e.g., variable coefficients) to make the scope of the optimal estimates clearer.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We address the two major comments point by point below.

read point-by-point responses

Referee: [stability and consistency analysis] The section deriving the discrete bilinear form and its coercivity: the claim that the consistency term alone (via the projection operators) yields a coercive and continuous form without stabilization is load-bearing for the optimal a priori estimates, yet no explicit lower bound independent of the mesh distortion parameter is shown for concave polygons; the argument implicitly assumes a discrete Poincaré inequality from the degrees of freedom, but this is not proven uniformly.

Authors: We thank the referee for this observation. The coercivity of the stabilization-free form is established in Theorem 3.2, where an explicit lower bound is derived that depends only on the polynomial degree and the mesh regularity constant from Assumption 2.1; this bound is independent of the distortion parameter and holds for concave polygons. The discrete Poincaré inequality is proven explicitly in Lemma 2.5 directly from the degrees of freedom, without implicit assumptions. We will add a clarifying remark after Theorem 3.2 in the revision to emphasize uniformity with respect to concavity. revision: partial
Referee: [numerical experiments] The numerical experiments section (concave-mesh tests): the reported results on concave polygons do not include a systematic check of the discrete coercivity constant (or its dependence on concavity/aspect ratio), so the robustness claim for the stabilization-free form remains unsupported by quantitative evidence.

Authors: We agree that a direct numerical check would strengthen the presentation. In the revised version we will add a new table in Section 5.3 reporting the computed discrete coercivity constants for the concave-mesh families, together with their dependence on the concavity parameter, confirming they remain bounded away from zero. revision: yes

Circularity Check

0 steps flagged

No circularity: error estimates derived from projection consistency without reduction to inputs

full rationale

The paper derives optimal a priori error estimates for the stabilization-free VEM on general second-order elliptic eigenvalue problems by establishing consistency and approximation properties of the virtual element spaces and projection operators on polygonal meshes. These estimates rest on standard VEM analysis techniques for the consistency term alone, with numerical experiments on convex, convex-concave, and concave meshes providing independent validation. No load-bearing step reduces a prediction to a fitted parameter, self-citation chain, or definitional equivalence by construction; the coercivity claim is presented as following from the chosen degrees of freedom and projections rather than being smuggled in via prior self-work or ansatz. The derivation chain is therefore self-contained against external VEM benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only abstract available, so no specific free parameters, axioms, or invented entities can be extracted; the method appears to rely on standard VEM assumptions without new postulates identified.

pith-pipeline@v0.9.0 · 5351 in / 966 out tokens · 32024 ms · 2026-05-13T16:52:16.457627+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Bh(uh,vh) := ∑_E aE_h(uh,vh) + … with aE_h = (K Π^P_∇ uh, Π^P_∇ vh)_E; coercivity asserted via Theorems 2–3 of [18] under Assumptions 2.1–2.2
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Assumption 2.1: each element star-shaped w.r.t. a ball of radius ≥ C_T h_E; N_E bounded

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

[1]

Beir˜ao da Veiga, F

L. Beir˜ao da Veiga, F. Brezzi, A. Cangiani, G. Manzini, L. D. Marini, A. Russo, Basic principles of virtual element methods,Math. Models Methods Appl. Sci.,23(2013), 199-214.http://dx.doi. org/10.1142/s0218202512500492

work page doi:10.1142/s0218202512500492 2013
[2]

Beir˜ao da Veiga, F

L. Beir˜ao da Veiga, F. Brezzi, L. D. Marini, A. Russo, The hitchhiker’s guide to the virtual element method,Math. Models Methods Appl. Sci.,24(2014), 1541-1573.http://dx.doi.org/10.1142/ S021820251440003X

work page 2014
[3]

Brezzi, R

F. Brezzi, R. S. Falk, L. D. Marini, Basic principles of mixed virtual element methods,ESAIM: Math. Model. Numer. Anal.,48(2014), 1227-1240.http://dx.doi.org/10.1051/m2an/2013138

work page doi:10.1051/m2an/2013138 2014
[4]

B. A. De Dios, K. Lipnikov, G. Manzini, The nonconforming virtual element method,ESAIM: Math. Model. Numer. Anal.,50(2016), 879-904.http://dx.doi.org/10.1051/m2an/2015090

work page doi:10.1051/m2an/2015090 2016
[5]

Beir˜ao da Veiga, F

L. Beir˜ao da Veiga, F. Brezzi, L. D. Marini, A. Russo, Serendipity nodal VEM spaces,Comput. Fluids, 141(2016), 2-12.http://dx.doi.org/10.1016/j.compfluid.2016.02.015

work page doi:10.1016/j.compfluid.2016.02.015 2016
[6]

Mascotto, I

L. Mascotto, I. Perugia, A. Pichler, Non-conforming harmonic virtual element method: h- and p-versions,J. Sci. Comput.,77(2018), 1874-1908.https://doi.org/10.1007/ s10915-018-0797-4

work page 2018
[7]

Benvenuti, A

E. Benvenuti, A. Chiozzi, G. Manzini, N. Sukumar, Extended virtual element method for the Laplace problem with singularities and discontinuities,Comput. Methods Appl. Mech. Eng.,356(2019), 571- 597.https://doi.org/10.1016/j.cma.2019.07.028

work page doi:10.1016/j.cma.2019.07.028 2019
[8]

S. Cao, L. Chen, R. Guo, F. Lin, Immersed virtual element methods for elliptic interface problems in two dimensions,J. Sci. Comput.,93(2022), 12.https://doi.org/10.1007/ s10915-022-01949-x

work page 2022
[9]

J. Zhao, S. Mao, B. Zhang, F. Wang, The interior penalty virtual element method for the biharmonic problem,Math. Comput.,92(2023), 1543-1574.https://doi.org/10.1090/mcom/3828

work page doi:10.1090/mcom/3828 2023
[10]

F. Chen, Q. Wang, Z. Zhou, Two-grid virtual element discretization of semilinear elliptic problem, Appl. Numer. Math.,186(2023), 228-240.https://doi.org/10.1016/j.apnum.2023.01.009

work page doi:10.1016/j.apnum.2023.01.009 2023
[11]

Boffi, F

D. Boffi, F. Gardini, L. Gastaldi, Approximation of PDE eigenvalue problems involving parameter dependent matrices,Calcolo,57(2020), 41.https://doi.org/10.1007/s10092-020-00390-6

work page doi:10.1007/s10092-020-00390-6 2020
[12]

Berrone, A

S. Berrone, A. Borio, F. Marcon, Comparison of standard and stabilization free virtual elements on anisotropic elliptic problems,Appl. Math. Lett.,129(2022), 107971.https://doi.org/10.1016/ j.aml.2022.107971

work page arXiv 2022
[13]

Beir˜ao da Veiga, C

L. Beir˜ao da Veiga, C. Canuto, R. H. Nochetto, M. Verani, Adaptive VEM: stabilization-free a posteriori error analysis and contraction property,SIAM J. Numer. Anal.,61(2023), 457-494. https://doi.org/10.1137/21M1458740

work page doi:10.1137/21m1458740 2023
[14]

Mascotto, The role of stabilization in the virtual element method: a survey,Comput

L. Mascotto, The role of stabilization in the virtual element method: a survey,Comput. Math. Appl., 151(2023), 244-251.https://doi.org/10.1016/j.camwa.2023.09.045

work page doi:10.1016/j.camwa.2023.09.045 2023
[15]

Alzaben, D

L. Alzaben, D. Boffi, A. Dedner, L. Gastaldi, On the stabilization of a virtual element method for an acoustic vibration problem,Math. Models Methods Appl. Sci.,35(2025), 655-701.https://doi. org/10.1142/s0218202525500071

work page doi:10.1142/s0218202525500071 2025
[17]

Berrone, A

S. Berrone, A. Borio, F. Marcon, G. Teora, A first-order stabilization-free virtual element method, Appl. Math. Lett.,142(2023), 108641.https://doi.org/10.1016/j.aml.2023.108641

work page doi:10.1016/j.aml.2023.108641 2023
[18]

Berrone, A

S. Berrone, A. Borio, D. Fassino, F. Marcon, Stabilization-free Virtual Element Method for 2D second order elliptic equations,Comput. Methods Appl. Mech. Eng.,438(2025), 117839.https://doi. org/10.1016/j.cma.2025.117839 13

work page doi:10.1016/j.cma.2025.117839 2025
[19]

Computers & Mathematics with Applications132, 48– 62 (2023) https://doi.org/10.1016/j.camwa

S. Berrone, A. Borio, F. Marcon, Lowest order stabilization free virtual element method for the 2D Poisson equation,Comput. Math. Appl.,177(2025), 78-99.https://doi.org/10.1016/j.camwa. 2024.11.017

work page doi:10.1016/j.camwa 2025
[20]

A. M. D’Altri, S. de Miranda, L. Patruno, E. Sacco, An enhanced VEM formulation for plane elastic- ity,Comput. Methods Appl. Mech. Eng.,376(2021), 113663.https://doi.org/10.1016/j.cma. 2020.113663

work page doi:10.1016/j.cma 2021
[21]

A. Chen, N. Sukumar, Stabilization-free serendipity virtual element method for plane elasticity,Com- put. Methods Appl. Mech. Eng.,404(2023), 115784.https://doi.org/10.1016/j.cma.2022. 115784

work page doi:10.1016/j.cma.2022 2023
[22]

Bertrand, C

F. Bertrand, C. Carstensen, B. Gräßle, N. T. Tran, Stabilization-free HHO a posteriori error control, Numer. Math.,154(2023), 369-408.https://doi.org/10.1007/s00211-023-01366-8

work page doi:10.1007/s00211-023-01366-8 2023
[23]

B. B. Xu, F. Peng, P. Wriggers, Stabilization-free virtual element method for finite strain applica- tions,Comput. Methods Appl. Mech. Eng.,417(2023), 116555.https://doi.org/10.1016/j. cma.2023.116555

work page doi:10.1016/j 2023
[24]

Lamperti, M

A. Lamperti, M. Cremonesi, U. Perego, A. Russo, C. Lovadina, A Hu-Washizu variational approach to self-stabilized virtual elements: 2D linear elastostatics,Comput. Mech.,71(2023), 935-955.https: //doi.org/10.1007/s00466-023-02282-2

work page doi:10.1007/s00466-023-02282-2 2023
[25]

Bouchez, A

T. Bouchez, A. Gravouil, N. Blal, A. Giacoma, E. Delor, J. D. Beley, A Hu-Washizu stabilization-free Virtual Element Method for 3D linear elasticity with star-convex polyhedrons,Comput. Methods Appl. Mech. Eng.,432(2024), 117420.https://doi.org/10.1016/j.cma.2024.117420

work page doi:10.1016/j.cma.2024.117420 2024
[26]

Borio, C

A. Borio, C. Lovadina, F. Marcon, M. Visinoni, A lowest order stabilization-free mixed virtual element method,Comput. Math. Appl.,160(2024), 161-170.https://doi.org/10.1016/j.camwa.2024. 02.024

work page doi:10.1016/j.camwa.2024 2024
[27]

Berrone, A

S. Berrone, A. Borio, F. Marcon, A stabilization-free virtual element method based on divergence-free projections,Comput. Methods Appl. Mech. Eng.,424(2024), 116885.https://doi.org/10.1016/ j.cma.2024.116885

work page arXiv 2024
[28]

J. Meng, X. Wang, L. L. Bu, L. Mei, A lowest-order free-stabilization virtual element method for the Laplacian eigenvalue problem,J. Comput. Appl. Math.,410(2022), 114013.https://doi.org/10. 1016/j.cam.2021.114013

work page arXiv 2022
[29]

Marcon, D

F. Marcon, D. Mora, A Stabilization-Free Virtual Element Method for the Convection- Diffusion Eigenproblem,J. Sci. Comput.,102(2025), 46.https://doi.org/10.1007/ s10915-024-02765-1

work page 2025
[30]

J. Meng, L. Guan, X. Qian, S. Song, L. Mei, Stabilization-Free Virtual Element Method for the Transmission Eigenvalue Problem on Anisotropic Media,J. Comput. Math.,44(2026), 103-134. https://doi.org/10.4208/jcm.2410-m2024-0023

work page doi:10.4208/jcm.2410-m2024-0023 2026
[31]

P. P. Foligno, D. Boffi, F. Credali, R. Vescovini, Benchmarking stabilized and self-stabilized p-virtual element methods with variable coefficientsComput. Methods Appl. Mech. Eng.,455(2026), 118863. https://doi.org/10.1016/j.cma.2026.118863

work page doi:10.1016/j.cma.2026.118863 2026
[32]

Cangiani, G

A. Cangiani, G. Manzini, O. J. Sutton, Conforming and nonconforming virtual element methods for el- liptic problems,IMA J. Numer. Anal.,37(2017), 1317-1354.https://doi.org/10.1093/imanum/ drw036

work page doi:10.1093/imanum/ 2017
[33]

Beir˜ao da Veiga, F

L. Beir˜ao da Veiga, F. Brezzi, L. D. Marini, A. Russo, Virtual element method for general second- order elliptic problems on polygonal meshes,Math. Models Methods Appl. Sci.,26(2016), 729-750. https://doi.org/10.1142/S0218202516500160

work page doi:10.1142/s0218202516500160 2016
[34]

Bernardi, R

C. Bernardi, R. Verfürth, Adaptive finite element methods for elliptic equations with non-smooth coefficients,Numer. Math.,85(2000), 579-608.https://doi.org/10.1007/PL00005393

work page doi:10.1007/pl00005393 2000
[35]

Grisvard,Elliptic problems in nonsmooth domains, Society for Industrial and Applied Mathematics, 2011.https://epubs.siam.org/doi/book/10.1137/1.9781611972030 14

P. Grisvard,Elliptic problems in nonsmooth domains, Society for Industrial and Applied Mathematics, 2011.https://epubs.siam.org/doi/book/10.1137/1.9781611972030 14

work page doi:10.1137/1.9781611972030 2011
[36]

Ahmad, A

B. Ahmad, A. Alsaedi, F. Brezzi, L. D. Marini, A. Russo, Equivalent projectors for virtual element methods,Comput. Math. Appl.,66(2013), 376-391.https://doi.org/10.1016/j.camwa.2013. 05.015

work page doi:10.1016/j.camwa.2013 2013
[37]

Babu ˇska J

I. Babu ˇska J. E. Osborn,Eigenvalue Problems, Handbook of Numerical Analysis, Vol.II, Finite El- ement Methods (Part 1), Edited by P. G. Ciarlet and J. L. Lions, North-Holland, Elsevier Science Publishers B.V , 1991.https://doi.org/10.1016/S1570-8659(05)80042-0

work page doi:10.1016/s1570-8659(05)80042-0 1991
[38]

D. Mora, G. Rivera, R. Rodríguez, A virtual element method for the Steklov eigenvalue prob- lem,Math. Models Methods Appl. Sci.,25(2015), 1421-1445.https://doi.org/10.1142/ S0218202515500372

work page 2015
[39]

Yu, mVEM: a MATLAB software package for the virtual element methods, arXiv preprint, 2022, 2204.01339.https://doi.org/10.48550/arXiv.2204.01339

Y . Yu, mVEM: a MATLAB software package for the virtual element methods, arXiv preprint, 2022, 2204.01339.https://doi.org/10.48550/arXiv.2204.01339

work page doi:10.48550/arxiv.2204.01339 2022
[40]

Beir˜ao da Veiga, F

L. Beir˜ao da Veiga, F. Brezzi, L.D. Marini, A. Russo, Virtual Element Implementation for Gen- eral Elliptic Equations. In: G. Barrenechea, F. Brezzi, A. Cangiani, E. Georgoulis, (eds) Build- ing Bridges: Connections and Challenges in Modern Approaches to Numerical Partial Differential Equations. Lecture Notes in Computational Science and Engineering,114(...

work page doi:10.1007/978-3-319-41640-3_2 2016