Sharp $L^2$-stability and gradient stability hierarchies of the $L^2$-Poincar\'e inequalities on Euclidean balls and Gaussian Poincar\'e inequality

We study the sharp $L^2$-Poincar\'e inequalities on Euclidean balls and with respect to the Gaussian measure, and we focus on their stability. In both settings, the sharp constant is the first nonzero eigenvalue of a self-adjoint operator and the set of optimizers is a finite dimensional linear space. Using the spectral decomposition, we establish sharp $L^2$-stability estimates and sharp gradient stability estimates, together with the stability of the stability inequalities in both norms, with explicit optimal constants and with the complete characterization of the equality cases. We then show that this process can be continued. The Poincar\'e deficit is equal to an infinite sum of $L^2$-distances to the increasing spaces of optimizers of the successive stability inequalities. It is also equal to an infinite sum of $L^2$-gradient distances to the same spaces. The coefficients of the first sum are the spectral gaps ${\mu }_j-{\mu }_{j-1}$, and the coefficients of the second sum are ${\mu }_1$ times the gaps of the reciprocal spectrum, ${\mu }_1\left({\mu }_{j-1}^{-1}-{\mu }_j^{-1}\right)$. On the Euclidean ball, all the constants are expressed through the roots of Bessel functions, and we prove that the first two spectral levels are given by the degree one and the degree two Neumann modes in all dimensions, while the third level is given by the first radial mode when $N=2,3$ and by the degree three modes when $N\ge 4$. In the Gaussian case, all the constants are rational numbers and the $L^2$-sum has all coefficients equal to $1$; we prove that this property characterizes the arithmetic spectrum.