In this paper we revisit a two-level discretization based on the Localized Orthogonal Decomposition (LOD). It was originally proposed in [P.Henning, A.M{\aa}lqvist, D.Peterseim. SIAM J. Numer. Anal.52-4:1525-1550, 2014] to compute ground states of Bose-Einstein condensates by finding discrete minimizers of the Gross-Pitaevskii energy functional. The established convergence rates for the method appeared however suboptimal compared to numerical observations and a proof of optimal rates in this setting remained open. In this paper we shall close this gap by proving optimal order error estimates for the $L^2$- and $H^1$-error between the exact ground state and discrete minimizers, as well as error estimates for the ground state energy and the ground state eigenvalue. In particular, the achieved convergence rates for the energy and the eigenvalue are of $6$th order with respect to the mesh size on which the discrete LOD space is based, without making any additional regularity assumptions. These high rates justify the use of very coarse meshes, which significantly reduces the computational effort for finding accurate approximations of ground states. In addition, we include numerical experiments that confirm the optimality of the new theoretical convergence rates, both for smooth and discontinuous potentials.
翻译:在本文中,我们重新审视了一种基于局部正交分解(LOD)的双层离散化。它最初是由[P.Henning, A.M{\aa}lqvist, D.Peterseim. SIAM J. Numer. Anal.52-4:1525-1550, 2014]提出的,用于通过寻找Gross-Pitaevskii能量泛函的离散极小化来计算玻色-爱因斯坦冷凝体的基态。然而,该方法的建立的收敛速率与数值观察相比似乎是亚最优的,并且在这种情况下获得最优速率的证明仍然未能解决。在本文中,我们将通过证明拟合光滑势和不光滑势的精确基态和离散极小化之间的$L^2$误差和$H^1$误差的最优误差估计,以及基态能量和基态本征值的误差估计,来弥补这个差距。特别的,能量和本征值的实现收敛速率都是6阶,关于离散LOD空间的网格尺寸,而不会进行任何额外的规则性假设。这些高速率证明了使用非常粗糙的网格是有合理性的,这显著降低了寻找准确的基态近似的计算工作量。此外,我们包括对光滑和不连续势的数值实验证实了新的理论收敛速率的最优性。