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
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.
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
- 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
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.
Referee Report
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)
- [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.
- [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)
- [abstract] The abstract contains the grammatical error 'we proof that' instead of 'we prove that'.
- [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
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
-
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
-
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
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
axioms (1)
- standard math Standard approximation and consistency properties of nonconforming finite-element spaces on Sobolev spaces
Reference graph
Works this paper leans on
-
[1]
Robert Altmann, Patrick Henning, and Daniel Peterseim,TheJ-method for the Gross-Pitaevskii eigenvalue problem, Numer. Math.148(2021), no. 3, 575–610. MR4293963
work page 2021
-
[2]
Luigi Ambrosio, Nicola Fusco, and Diego Pallara,Functions of bounded variation and free discontinuity problems, Oxford University Press, 2000
work page 2000
-
[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
work page 2017
-
[4]
Ivo Babuˇ ska,Error-bounds for finite element method, Numerische Mathematik16(1971), no. 4, 322–333
work page 1971
-
[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
work page 2014
-
[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
work page 2013
-
[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
work page 2010
-
[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
work page 2004
-
[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
work page 2008
-
[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
work page 2005
-
[11]
Bose,Plancks gesetz und lichtquantenhypothese, Zeitschrift f¨ ur Physik26(1924), no. 1, 178–181
work page 1924
-
[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
work page 2018
-
[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
work page 2010
-
[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
work page 2023
-
[15]
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
work page 2014
-
[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
work page 2013
-
[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
work page 2011
-
[18]
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
work page 2011
-
[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
work page 2024
-
[20]
,Fully discretized Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem, Math. Comp.94(2025), no. 356, 2723–2760. MR4940541
work page 2025
-
[21]
Philippe G. Ciarlet,The finite element method for elliptic problems, Society for Industrial and Applied Mathematics, 2002
work page 2002
-
[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
work page 2010
-
[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
work page 2010
-
[24]
Daniele Antonio Di Pietro and Alexandre Ern,Mathematical aspects of discontinuous galerkin methods, 1st ed., Springer Berlin, Heidelberg, 2012
work page 2012
-
[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
work page 2007
-
[26]
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
work page 1914
-
[27]
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
work page 2022
- [28]
-
[29]
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
work page 2021
-
[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
work page 2014
-
[31]
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
work page 2023
-
[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
work page 2020
-
[33]
Patrick Henning and Mahima Yadav,On discrete ground states of rotating Bose-Einstein condensates, Math. Comp.94 (2025), no. 351, 1–32. MR4807803
work page 2025
-
[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
work page 2014
-
[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
work page 2014
-
[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
work page 2010
-
[37]
Andris Lasis and Endre S¨ uli,Poincar´ e-type inequalities for broken sobolev spaces, Isaac Newton Institute for Mathematical Science, Cambridge, UK, 2003
work page 2003
-
[38]
Youai Li,Lower approximation of eigenvalue by the nonconforming finite element method, Math. Numer. Sin30(2008), no. 2, 195–200
work page 2008
-
[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
work page 2022
-
[40]
1, Science Press, Beijing, 2007
Qun Lin and Jiafu Lin,Finite element methods: accuracy and improvement, Vol. 1, Science Press, Beijing, 2007
work page 2007
-
[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
work page 2005
-
[42]
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
work page 2003
-
[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
work page 2003
-
[44]
R. Rannacher and S. Turek,Simple nonconforming quadrilateral Stokes element, Numer. Methods Partial Differential Equations8(1992), no. 2, 97–111. MR1148797
work page 1992
-
[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
work page 2019
-
[46]
Nian Shao and Wenbin Chen,Riemannian acceleration with preconditioning for symmetric eigenvalue problems, Numer. Math.157(2025), no. 1, 307–354. MR4858137
work page 2025
-
[47]
Zhongci Shi and Ming Wang,Finite element methods, Science Press, Beijing, 2010
work page 2010
-
[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]
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
work page 2017
-
[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
work page 2016
-
[51]
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
work page 2021
-
[52]
Yidu Yang,A posteriori error estimates in Adini finite element for eigenvalue problems, J. Comput. Math.18(2000), no. 4, 413–418. MR1773912
work page 2000
-
[53]
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
work page 2010
-
[54]
Chen Zhang, Patrick Henning, Mahima Yadav, and Wenbin Chen,Convergence analysis of sobolev gradient flows for the rotating gross-pitaevskii energy functional, 2025
work page 2025
-
[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
work page 2019
-
[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
work page 2007
-
[57]
Aihui Zhou,An analysis of finite-dimensional approximations for the ground state solution of Bose-Einstein condensates, Nonlinearity17(2004), no. 2, 541–550. MR2039057
work page 2004
-
[58]
,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...
work page 2007
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.