pith. sign in

arxiv: 2605.23334 · v1 · pith:4FSHLZDDnew · submitted 2026-05-22 · 🧮 math.NA · cs.NA

Nonconforming Finite Element Approximation and Energy Lower Bound Estimation for the Gross--Pitaevskii Energy Functional

Pith reviewed 2026-05-25 03:52 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Gross-Pitaevskii functionalnonconforming finite elementsenergy lower boundBose-Einstein condensatea priori error estimatesground state approximationEQ1rot element
0
0 comments X

The pith

The EQ_1^rot nonconforming finite element gives a lower bound estimate for the exact Gross-Pitaevskii energy.

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

This paper examines the minimization of the Gross-Pitaevskii energy functional using nonconforming finite element methods to approximate the ground state of Bose-Einstein condensates. It derives a priori error estimates for the discrete energy, eigenvalue, and ground state, with specific rates for the EQ_1^rot element. The key result is that the discrete energy serves as a lower bound to the continuous energy within this framework. A sympathetic reader would care because this provides a way to compute approximations that are guaranteed to be below the true minimum, aiding in reliable numerical simulations of quantum systems.

Core claim

Within the EQ_1^{rot} finite element framework, the discrete ground state energy provides a lower bound estimation to the exact energy. The paper also establishes a priori error estimates for the discrete ground state energy, the discrete eigenvalue, and the discrete ground state, with explicit convergence rates for the EQ_1^{rot} element.

What carries the argument

The nonconforming EQ_1^{rot} finite element space and its approximation and consistency properties that enable the lower bound proof.

If this is right

  • The discrete ground state energy converges to the exact energy from below.
  • Explicit convergence rates are obtained for the error in energy, eigenvalue, and ground state.
  • Numerical experiments confirm the theoretical lower bound property.

Where Pith is reading between the lines

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

  • This lower bound property could be used to certify the accuracy of computations by checking how close the discrete energy is to known upper bounds from other methods.
  • Similar nonconforming elements might yield lower bounds in other constrained variational problems.
  • The approach suggests potential for developing guaranteed error estimators based on the energy difference.

Load-bearing premise

The continuous Gross-Pitaevskii problem has a solution with sufficient regularity, and the nonconforming space satisfies the needed approximation and consistency properties.

What would settle it

Finding a specific Gross-Pitaevskii problem with known exact energy where the computed discrete energy from EQ1rot exceeds that exact value would disprove the lower bound claim.

Figures

Figures reproduced from arXiv: 2605.23334 by Chen Zhang, Heyan Zhu, Wenbin Chen.

Figure 1
Figure 1. Figure 1: illustrates the transformation of the discrete ground state uh as the mesh is refined. On coarser grids, the discontinuities of the EQrot 1 element across element interfaces are clearly visible, providing a direct visual representation of the nonconforming nature of the space. As the mesh size h decreases, the discrete solution gradually converges toward a smooth, exact ground state. (a) 8 × 8 (b) 16 × 16 … view at source ↗
Figure 2
Figure 2. Figure 2: Convergence results for Example 6.3. 6.4. Asymptotic convergence under complicated potentials. Finally, we further evaluate the algo￾rithm using Case II from [8, Example 3]. The domain is D = [−8, 8]2 , with β = 400, and the trapping potential consists of an isotropic harmonic potential combined with a stirrer potential: V (x) = x 2 1 + x 2 2 + 8e −(x1−1)2−x 2 2 . The results are depicted in [PITH_FULL_IM… view at source ↗
Figure 3
Figure 3. Figure 3: Convergence results for Example 6.4. theoretically demonstrated that the EQrot 1 element provides a lower bound approximation to the exact ground state energy when the mesh size is sufficiently small. This lower bound property serves as a complement to the upper bound estimates typically produced by conforming finite element methods. Finally, numerical results have validated our theoretical analysis, demon… view at source ↗
read the original abstract

The ground state of Bose--Einstein condensates can be described as the minimizer of the Gross--Pitaevskii energy functional subject to a mass conservation constraint. In this paper, we study the corresponding discrete optimization problem in nonconforming finite element spaces and establish a priori error estimates for the discrete ground state energy, the discrete eigenvalue, and the discrete ground state. Specifically, we derive explicit convergence rates for the a priori error in the particular case of the $EQ_1^{\mathrm{rot}}$ finite element. Furthermore, we proof that within the $EQ_1^{\mathrm{rot}}$ finite element framework, the discrete ground state energy provides a lower bound estimation to the exact energy. Finally, numerical experiments are presented to validate the theoretical analysis.

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.

Referee Report

2 major / 2 minor

Summary. The paper develops a nonconforming finite element discretization of the Gross-Pitaevskii energy minimization problem using the EQ_1^rot element. It claims a priori error estimates (with explicit rates) for the discrete ground-state energy, the associated eigenvalue, and the ground-state function itself; proves that the discrete energy furnishes a lower bound on the continuous energy; and presents numerical experiments supporting the analysis.

Significance. If the lower-bound result and the error estimates are rigorously established, the work supplies a practical a-priori bound on the ground-state energy that is useful for validation in Bose-Einstein condensate simulations. The explicit rates for the particular EQ_1^rot element and the numerical confirmation add concrete value to the existing literature on nonconforming methods for nonlinear eigenvalue problems.

major comments (2)
  1. [proof of lower bound / abstract] The lower-bound claim for the discrete energy (abstract and the dedicated proof section) is load-bearing and rests on two unverified hypotheses: (i) the continuous minimizer possesses H^2 regularity (or higher) sufficient to control the interpolation error in the broken H^1 seminorm plus the L^2 and quartic terms, and (ii) the specific consistency identity or quadrature rule of the EQ_1^rot space makes the discrete energy of the interpolant strictly less than or equal to the continuous energy. Neither property is automatic for the nonlinear functional; the manuscript must exhibit the precise consistency identity used for the quartic term.
  2. [error estimates section] The a-priori error estimates for the energy, eigenvalue, and eigenfunction (theorems stated after the abstract) are derived under the same regularity and consistency assumptions. If the lower-bound argument is incomplete, the error-analysis section inherits the same gap; explicit constants or rates cannot be claimed without confirming that the broken-norm consistency error for the nonlinear term vanishes at the required order.
minor comments (2)
  1. [abstract] The abstract contains the grammatical error 'we proof that' instead of 'we prove that'.
  2. [preliminaries] Notation for the broken H^1 seminorm and the discrete energy functional should be introduced once and used consistently; several passages mix |·|_{1,h} with the full discrete energy without clear cross-reference.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation of the work's significance and for the detailed comments. We agree that the lower-bound proof requires an explicit consistency identity for the quartic term and will revise the manuscript accordingly to close this gap while preserving the claimed results.

read point-by-point responses
  1. Referee: [proof of lower bound / abstract] The lower-bound claim for the discrete energy (abstract and the dedicated proof section) is load-bearing and rests on two unverified hypotheses: (i) the continuous minimizer possesses H^2 regularity (or higher) sufficient to control the interpolation error in the broken H^1 seminorm plus the L^2 and quartic terms, and (ii) the specific consistency identity or quadrature rule of the EQ_1^rot space makes the discrete energy of the interpolant strictly less than or equal to the continuous energy. Neither property is automatic for the nonlinear functional; the manuscript must exhibit the precise consistency identity used for the quartic term.

    Authors: We accept that the proof section must be strengthened. The H^2 regularity of the continuous minimizer follows from standard elliptic regularity theory for the Gross-Pitaevskii equation under the usual assumptions on the potential (smooth and confining) and domain (convex polygonal); we will add a brief justification citing this. For the second point, the revised manuscript will include an explicit lemma deriving the consistency identity for the quartic term under the EQ_1^rot interpolation operator, showing via direct expansion and the element's quadrature rule that the discrete energy of the interpolant is indeed ≤ the continuous energy. This identity will be stated precisely and used to complete the lower-bound argument. revision: yes

  2. Referee: [error estimates section] The a-priori error estimates for the energy, eigenvalue, and eigenfunction (theorems stated after the abstract) are derived under the same regularity and consistency assumptions. If the lower-bound argument is incomplete, the error-analysis section inherits the same gap; explicit constants or rates cannot be claimed without confirming that the broken-norm consistency error for the nonlinear term vanishes at the required order.

    Authors: The error estimates rely on the same consistency framework as the lower bound. Once the explicit identity for the quartic term is added (as described above), we will verify in the error-analysis section that the broken-norm consistency error for the nonlinear term is of the order needed to support the stated rates. The theorems will be updated to reference this identity, and any intermediate estimates will be adjusted if necessary to maintain the explicit convergence rates for the EQ_1^rot element. revision: yes

Circularity Check

0 steps flagged

No circularity in derivation chain

full rationale

The paper performs a priori analysis and proofs for error estimates and a lower-bound property of the discrete energy in the EQ_1^rot nonconforming space for the Gross-Pitaevskii functional. The abstract and provided excerpts describe standard FEM consistency arguments, regularity hypotheses on the continuous minimizer, and explicit convergence rates derived from approximation properties; none of these steps reduce by construction to self-definitions, fitted parameters renamed as predictions, or load-bearing self-citations. The lower-bound claim is presented as a theorem proved from the element's consistency identity and the continuous problem's properties, which are external to the discrete construction itself. This is a self-contained mathematical derivation with no evident circular reductions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Paper relies on standard Sobolev-space theory and finite-element approximation properties; no free parameters or invented entities are visible from the abstract.

axioms (1)
  • standard math Standard approximation and consistency properties of nonconforming finite-element spaces on Sobolev spaces
    Invoked for a priori error estimates and lower-bound proof.

pith-pipeline@v0.9.0 · 5663 in / 1036 out tokens · 22544 ms · 2026-05-25T03:52:15.525885+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

58 extracted references · 58 canonical work pages

  1. [1]

    Math.148(2021), no

    Robert Altmann, Patrick Henning, and Daniel Peterseim,TheJ-method for the Gross-Pitaevskii eigenvalue problem, Numer. Math.148(2021), no. 3, 575–610. MR4293963

  2. [2]

    Luigi Ambrosio, Nicola Fusco, and Diego Pallara,Functions of bounded variation and free discontinuity problems, Oxford University Press, 2000

  3. [3]

    Xavier Antoine, Antoine Levitt, and Qinglin Tang,Efficient spectral computation of the stationary states of rotating Bose-Einstein condensates by preconditioned nonlinear conjugate gradient methods, J. Comput. Phys.343(2017), 92–109. MR3654050

  4. [4]

    4, 322–333

    Ivo Babuˇ ska,Error-bounds for finite element method, Numerische Mathematik16(1971), no. 4, 322–333

  5. [5]

    Weizhu Bao,Mathematical models and numerical methods for Bose-Einstein condensation, Proceedings of the International Congress of Mathematicians—Seoul 2014. Vol. IV, 2014, pp. 971–996. MR3751161

  6. [6]

    Weizhu Bao and Yongyong Cai,Mathematical theory and numerical methods for Bose-Einstein condensation, Kinet. Relat. Models6(2013), no. 1, 1–135. MR3005624 FINITE ELEMENT APPROXIMATION FOR GROSS–PITAEVSKII ENERGY FUNCTIONAL 21

  7. [7]

    Weizhu Bao, Yongyong Cai, and Hanquan Wang,Efficient numerical methods for computing ground states and dynamics of dipolar Bose-Einstein condensates, J. Comput. Phys.229(2010), no. 20, 7874–7892. MR2674308

  8. [8]

    Weizhu Bao and Qiang Du,Computing the ground state solution of Bose-Einstein condensates by a normalized gradient flow, SIAM J. Sci. Comput.25(2004), no. 5, 1674–1697. MR2087331

  9. [9]

    Weizhu Bao and Jie Shen,A generalized-Laguerre-Hermite pseudospectral method for computing symmetric and central vortex states in Bose-Einstein condensates, J. Comput. Phys.227(2008), no. 23, 9778–9793. MR2469033

  10. [10]

    Markowich,Ground, symmetric and central vortex states in rotating Bose- Einstein condensates, Commun

    Weizhu Bao, Hanquan Wang, and Peter A. Markowich,Ground, symmetric and central vortex states in rotating Bose- Einstein condensates, Commun. Math. Sci.3(2005), no. 1, 57–88. MR2132826

  11. [11]

    1, 178–181

    Bose,Plancks gesetz und lichtquantenhypothese, Zeitschrift f¨ ur Physik26(1924), no. 1, 178–181

  12. [12]

    Eric Canc` es, Rachida Chakir, Lianhua He, and Yvon Maday,Two-grid methods for a class of nonlinear elliptic eigenvalue problems, IMA J. Numer. Anal.38(2018), no. 2, 605–645. MR3800034

  13. [13]

    Eric Canc` es, Rachida Chakir, and Yvon Maday,Numerical analysis of nonlinear eigenvalue problems, J. Sci. Comput.45 (2010), no. 1-3, 90–117. MR2679792

  14. [14]

    Haifan Chen, Guozhi Dong, Wei Liu, and Ziqing Xie,Second-order flows for computing the ground states of rotating Bose-Einstein condensates, J. Comput. Phys.475(2023), Paper No. 111872, 28. MR4530846

  15. [15]

    Simul.12(2014), no

    Huajie Chen, Xiaoying Dai, Xingao Gong, Lianhua He, and Aihui Zhou,Adaptive finite element approximations for Kohn- Sham models, Multiscale Model. Simul.12(2014), no. 4, 1828–1869. MR3284557

  16. [16]

    Huajie Chen, Xingao Gong, Lianhua He, Zhang Yang, and Aihui Zhou,Numerical analysis of finite dimensional approxi- mations of Kohn-Sham models, Adv. Comput. Math.38(2013), no. 2, 225–256. MR3019148

  17. [17]

    Huajie Chen, Xingao Gong, Lianhua He, and Aihui Zhou,Adaptive finite element approximations for a class of nonlinear eigenvalue problems in quantum physics, Adv. Appl. Math. Mech.3(2011), no. 4, 493–518. MR2823345

  18. [18]

    Methods Appl

    Huajie Chen, Lianhua He, and Aihui Zhou,Finite element approximations of nonlinear eigenvalue problems in quantum physics, Comput. Methods Appl. Mech. Engrg.200(2011), no. 21-22, 1846–1865. MR2787541

  19. [19]

    Ziang Chen, Jianfeng Lu, Yulong Lu, and Xiangxiong Zhang,On the convergence of Sobolev gradient flow for the Gross- Pitaevskii eigenvalue problem, SIAM J. Numer. Anal.62(2024), no. 2, 667–691. MR4712406

  20. [20]

    Comp.94(2025), no

    ,Fully discretized Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem, Math. Comp.94(2025), no. 356, 2723–2760. MR4940541

  21. [21]

    Ciarlet,The finite element method for elliptic problems, Society for Industrial and Applied Mathematics, 2002

    Philippe G. Ciarlet,The finite element method for elliptic problems, Society for Industrial and Applied Mathematics, 2002

  22. [22]

    Ionut Danaila and Fr´ ed´ eric Hecht,A finite element method with mesh adaptivity for computing vortex states in fast-rotating Bose-Einstein condensates, J. Comput. Phys.229(2010), no. 19, 6946–6960. MR2677760

  23. [23]

    Ionut Danaila and Parimah Kazemi,A new Sobolev gradient method for direct minimization of the Gross-Pitaevskii energy with rotation, SIAM J. Sci. Comput.32(2010), no. 5, 2447–2467. MR2684722

  24. [24]

    Daniele Antonio Di Pietro and Alexandre Ern,Mathematical aspects of discontinuous galerkin methods, 1st ed., Springer Berlin, Heidelberg, 2012

  25. [25]

    Dion and Eric Canc` es,Ground state of the time-independent Gross-Pitaevskii equation, Comput

    Claude M. Dion and Eric Canc` es,Ground state of the time-independent Gross-Pitaevskii equation, Comput. Phys. Comm. 177(2007), no. 10, 787–798. MR2578134

  26. [26]

    zweite abhandlung, Albert Einstein: Akademie-Vortr¨ age: Sitzungsberichte der Preußischen Akademie der Wissenschaften 1914–1932 (2005), 245–257

    Albert Einstein,Quantentheorie des einatomigen idealen gases. zweite abhandlung, Albert Einstein: Akademie-Vortr¨ age: Sitzungsberichte der Preußischen Akademie der Wissenschaften 1914–1932 (2005), 245–257

  27. [27]

    Engstr¨ om, S

    C. Engstr¨ om, S. Giani, and L. Grubiˇ si´ c,Higher order composite DG approximations of Gross-Pitaevskii ground state: benchmark results and experiments, J. Comput. Appl. Math.400(2022), Paper No. 113652, 15. MR4293703

  28. [28]

    Moritz Hauck and Yizhou Liang,A hybrid high-order method for the gross–pitaevskii eigenvalue problem, arXiv preprint arXiv:2506.19944 (2025)

  29. [29]

    Wihler,Gradient flow finite element discretizations with energy-based adaptivity for the Gross-Pitaevskii equation, J

    Pascal Heid, Benjamin Stamm, and Thomas P. Wihler,Gradient flow finite element discretizations with energy-based adaptivity for the Gross-Pitaevskii equation, J. Comput. Phys.436(2021), Paper No. 110165, 15. MR4234216

  30. [30]

    Patrick Henning, Axel M˚ a lqvist, and Daniel Peterseim,Two-level discretization techniques for ground state computations of Bose-Einstein condensates, SIAM J. Numer. Anal.52(2014), no. 4, 1525–1550. MR3227466

  31. [31]

    Simul.21(2023), no

    Patrick Henning and Anna Persson,On optimal convergence rates for discrete minimizers of the Gross-Pitaevskii energy in localized orthogonal decomposition spaces, Multiscale Model. Simul.21(2023), no. 3, 993–1011. MR4619875

  32. [32]

    Patrick Henning and Daniel Peterseim,Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem: global conver- gence and computational efficiency, SIAM J. Numer. Anal.58(2020), no. 3, 1744–1772. MR4108207

  33. [33]

    Comp.94 (2025), no

    Patrick Henning and Mahima Yadav,On discrete ground states of rotating Bose-Einstein condensates, Math. Comp.94 (2025), no. 351, 1–32. MR4807803

  34. [34]

    Jun Hu, Yunqing Huang, and Qun Lin,Lower bounds for eigenvalues of elliptic operators: by nonconforming finite element methods, J. Sci. Comput.61(2014), no. 1, 196–221. MR3254372 22 CHEN ZHANG, HEYAN ZHU, AND WENBIN CHEN

  35. [35]

    Elias Jarlebring, Simen Kvaal, and Wim Michiels,An inverse iteration method for eigenvalue problems with eigenvector nonlinearities, SIAM J. Sci. Comput.36(2014), no. 4, A1978–A2001. MR3249374

  36. [36]

    Parimah Kazemi and Michael Eckart,Minimizing the Gross-Pitaevskii energy functional with the Sobolev gradient— analytical and numerical results, Int. J. Comput. Methods7(2010), no. 3, 453–475. MR2736464

  37. [37]

    Andris Lasis and Endre S¨ uli,Poincar´ e-type inequalities for broken sobolev spaces, Isaac Newton Institute for Mathematical Science, Cambridge, UK, 2003

  38. [38]

    Youai Li,Lower approximation of eigenvalue by the nonconforming finite element method, Math. Numer. Sin30(2008), no. 2, 195–200

  39. [39]

    Qigang Liang, Xuejun Xu, and Liuyao Yuan,A weak Galerkin finite element method can compute both upper and lower eigenvalue bounds, J. Sci. Comput.93(2022), no. 1, Paper No. 19, 21. MR4476865

  40. [40]

    1, Science Press, Beijing, 2007

    Qun Lin and Jiafu Lin,Finite element methods: accuracy and improvement, Vol. 1, Science Press, Beijing, 2007

  41. [41]

    Qun Lin, Lutz Tobiska, and Aihui Zhou,Superconvergence and extrapolation of non-conforming low order finite elements applied to the Poisson equation, IMA J. Numer. Anal.25(2005), no. 1, 160–181. MR2110239

  42. [42]

    Math.94(2003), no

    Yvon Maday and Gabriel Turinici,Error bars and quadratically convergent methods for the numerical simulation of the Hartree-Fock equations, Numer. Math.94(2003), no. 4, 739–770. MR1990591

  43. [43]

    116, The Clarendon Press, Oxford University Press, Oxford, 2003

    Lev Pitaevskii and Sandro Stringari,Bose-Einstein condensation, International Series of Monographs on Physics, vol. 116, The Clarendon Press, Oxford University Press, Oxford, 2003. MR2012737

  44. [44]

    Rannacher and S

    R. Rannacher and S. Turek,Simple nonconforming quadrilateral Stokes element, Numer. Methods Partial Differential Equations8(1992), no. 2, 97–111. MR1148797

  45. [45]

    Sergey Repin,Poincar´ e type inequalities for vector functions with zero mean normal traces on the boundary and ap- plications to interpolation methods, Contributions to partial differential equations and applications, 2019, pp. 411–432. MR3837916

  46. [46]

    Math.157(2025), no

    Nian Shao and Wenbin Chen,Riemannian acceleration with preconditioning for symmetric eigenvalue problems, Numer. Math.157(2025), no. 1, 307–354. MR4858137

  47. [47]

    Zhongci Shi and Ming Wang,Finite element methods, Science Press, Beijing, 2010

  48. [48]

    3, 449–471, available athttps://doi.org/10.1137/0716037

    Friedrich Stummel,The generalized patch test, SIAM Journal on Numerical Analysis16(1979), no. 3, 449–471, available athttps://doi.org/10.1137/0716037

  49. [49]

    Xinming Wu, Zaiwen Wen, and Weizhu Bao,A regularized Newton method for computing ground states of Bose-Einstein condensates, J. Sci. Comput.73(2017), no. 1, 303–329. MR3704854

  50. [50]

    Hehu Xie and Manting Xie,A multigrid method for ground state solution of Bose-Einstein condensates, Commun. Comput. Phys.19(2016), no. 3, 648–662. MR3480951

  51. [51]

    2, 645–663

    Fei Xu, Hehu Xie, Manting Xie, and Meiling Yue,A multigrid method for the ground state solution of Bose-Einstein condensates based on Newton iteration, BIT61(2021), no. 2, 645–663. MR4258031

  52. [52]

    Yidu Yang,A posteriori error estimates in Adini finite element for eigenvalue problems, J. Comput. Math.18(2000), no. 4, 413–418. MR1773912

  53. [53]

    China Math.53(2010), no

    Yidu Yang, Zhimin Zhang, and Fubiao Lin,Eigenvalue approximation from below using non-conforming finite elements, Sci. China Math.53(2010), no. 1, 137–150. MR2594754

  54. [54]

    Chen Zhang, Patrick Henning, Mahima Yadav, and Wenbin Chen,Convergence analysis of sobolev gradient flows for the rotating gross-pitaevskii energy functional, 2025

  55. [55]

    Ning Zhang, Fei Xu, and Hehu Xie,An efficient multigrid method for ground state solution of Bose-Einstein condensates, Int. J. Numer. Anal. Model.16(2019), no. 5, 789–803. MR3946658

  56. [56]

    Zhimin Zhang, Yidu Yang, and Zhen Chen,Eigenvalue approximation from below by Wilson’s element, Math. Numer. Sin. 29(2007), no. 3, 319–321. MR2370469

  57. [57]

    2, 541–550

    Aihui Zhou,An analysis of finite-dimensional approximations for the ground state solution of Bose-Einstein condensates, Nonlinearity17(2004), no. 2, 541–550. MR2039057

  58. [58]

    Methods Appl

    ,Finite dimensional approximations for the electronic ground state solution of a molecular system, Math. Methods Appl. Sci.30(2007), no. 4, 429–447. MR2293570 AppendixA.Properties of discrete ground states In this section, we present several lemmas concerning the finite element spaceV h and its associated discrete ground statesu h. Lemma A.1 provides an e...