Quantitative Kröger inequalities for Neumann eigenvalues of convex domains

这篇论文《凸域 Neumann 特征值的定量 Kröger 不等式》（Quantitative Kröger inequalities for Neumann eigenvalues of convex domains）由 Dorin Bucur、Andrea Gentile 和 Antoine Henrot 撰写。文章主要研究了凸域上 Neumann-Laplacian 算子特征值的定量稳定性问题，特别是针对 Kröger 提出的 sharp 上界进行了改进，给出了包含几何“平坦度”项的定量不等式。

以下是该论文的详细技术总结：

1. 研究问题 (Problem)

• 背景：对于 $\mathbb{R}^d$ 中的有界开凸集 $\Omega$，Neumann-Laplacian 的特征值问题定义为 $-\Delta u = \mu u$ 且 $\frac{\partial u}{\partial \nu} = 0$。特征值按重数排序为 $0 = \mu_0 < \mu_1 \le \mu_2 \le \dots$。
• 已知结果：Kröger (1999) 证明了对于直径为 $D_\Omega$ 的凸域，第 $k$ 个特征值满足 sharp 上界：
  $$D_\Omega^2 \mu_k(\Omega) < \mu^*_{k,d}$$
  其中 $\mu^*_{k,d}$ 是仅依赖于 $k$ 和 $d$ 的常数。
• 核心问题：上述不等式是严格的（即等号在 $d \ge 2$ 时无法取到），等号仅在一系列“坍缩”（collapsing）到线段的域序列中渐近达到。自然的问题是：如果一个凸域的特征值非常接近这个上界，那么该域的几何形状必须多么“扁平”？即，能否给出一个定量的不等式，其中包含一个衡量域与线段（或一维结构）距离的项？
• 目标：建立形如 $D_\Omega^2 \mu_k(\Omega) \le \mu^*_{k,d} - C(k, d) \frac{a_2(\Omega)^2}{D_\Omega^2}$ 的不等式，其中 $a_2(\Omega)$ 是 $\Omega$ 的 John 椭球体的第二长半轴（衡量平坦度）。

2. 方法论 (Methodology)

论文采用了反证法结合变分分析和渐近分析的方法：

1. 反证假设：假设存在一个序列 $\{\Omega_n\}$，其直径固定为 1，且特征值 $\mu_k(\Omega_n)$ 趋近于上界 $\mu^*_{k,d}$ 的速度快于 $a_2(\Omega_n)^2$（即假设不等式不成立）。
2. 坍缩分析 (Collapsing Domains)：
   • 利用 Hausdorff 距离的紧性，分析序列 $\Omega_n$ 的极限行为。
   • 证明如果特征值接近上界，$\Omega_n$ 必须坍缩到一个低维流形（对于 $k=1$ 是线段）。
   • 引入 John 椭球体的半轴 $a_i(\Omega_n)$，通过重缩放变换 $T_n$ 将 $\Omega_n$ 映射为非退化的序列 $\tilde{\Omega}_n$，使其收敛到一个极限集 $\tilde{\Omega}$。
3. Sturm-Liouville 问题与加权特征值：
   • 利用 Brunn-Minkowski 不等式，将高维特征值问题转化为加权 Sturm-Liouville 问题（权重函数 $p(x)$ 与截面体积有关）。
   • 证明在极限情况下，加权特征值 $\mu_k(\omega, p)$ 必须等于 Kröger 的上界 $\mu^*_{k,d}$。
   • 利用 Kröger 极值问题的解的结构（$p^{1/(d-1)}$ 是分段线性的，对应圆锥截断结构），推导出极限几何形状必须是特定的锥形结构。
4. 变分扰动 (Variational Perturbation)：
   • 构造一个特殊的测试函数 $w_n$，基于极限域的第一特征函数 $u$ 进行扰动，形式为 $w_n \approx u(x_1)(1 + \tau (x_2 - c(x_1))^2)$。
   • 将 $w_n$ 代入 Rayleigh 商，展开关于小参数 $a_2(\Omega_n)$ 的级数。
   • 通过精细的积分估计，证明如果假设成立，则会导致矛盾（即导出的不等式方向与假设相反），从而证明原不等式成立。
5. 显式常数计算：
   • 对于 $d=2, k=1$ 且域具有对称轴的情况，通过求解一个具体的变分问题（最小化涉及 Bessel 函数的泛函），给出了常数 $C(1, 2)$ 的显式下界。

3. 主要贡献与结果 (Key Contributions & Results)

定理 1.1 (定量 Kröger 不等式)

对于任意 $k, d \in \mathbb{N}$ ($d \ge 2$)，存在仅依赖于 $k, d$ 的常数 $C(k, d) > 0$，使得对于任意有界凸集 $\Omega \subset \mathbb{R}^d$：
$$D_\Omega^2 \mu_k(\Omega) \le \mu^*_{k,d} - C(k, d) \frac{a_2(\Omega)^2}{D_\Omega^2}$$
其中 $a_2(\Omega)$ 是 John 椭球的第二长半轴。

• 意义：这是首次将 Kröger 不等式推广到包含几何平坦度项的定量形式，且适用于任意阶特征值 $k$。

定理 1.2 (指数最优性)

在 $d=2, k=1$ 的情况下，上述不等式中平坦度项 $a_2(\Omega)$ 的指数 2 是最优的（sharp）。

• 证明：通过构造超等边三角形（superequilateral triangle）序列，当三角形变扁时，特征值的衰减率恰好与 $w^2$（宽度的平方）同阶。

定理 1.3 (平面情形显式常数)

对于关于直径中垂线对称的平面凸域 $\Omega$，给出了显式常数：
$$D_\Omega^2 \mu_1(\Omega) \le 4j_{0,1}^2 - 0.432 \frac{w_\Omega^2}{D_\Omega^2}$$
其中 $j_{0,1}$ 是第一类 Bessel 函数 $J_0$ 的第一个正零点，$w_\Omega$ 是垂直于直径方向的宽度。

• 对比：虽然 0.432 不是最优常数，但相比于之前文献中关于类似问题的常数（如 $10^{-16}$ 量级），这是一个显著改进且可计算的数值。

4. 技术细节与难点

• 几何刚性：证明的关键在于利用 Brunn-Minkowski 不等式的等号成立条件，推导出当特征值接近上界时，域的截面必须具有特定的线性比例关系（即域必须是截断的圆锥）。
• 测试函数的构造：构造的测试函数不仅包含一维特征函数，还引入了垂直于主轴的二次项扰动。这一项捕捉到了域在“扁平”方向上的几何变化对特征值的二阶影响。
• 积分估计：在展开 Rayleigh 商时，需要处理分子和分母中关于 $a_2$ 的高阶项，并证明交叉项的符号和大小足以产生矛盾。
• 变分问题求解：在计算显式常数时，作者解决了一个涉及 Bessel 函数和凹函数约束的变分问题，证明了最优的轮廓函数是线性的（$h(x)=2x$），这对应于三角形域。

5. 意义与影响 (Significance)

1. 填补理论空白：在等周不等式和特征值稳定性研究中，大多数结果关注的是域与球体的距离（Fraenkel 不对称性）。本文处理的是非紧的极值情况（域坍缩到线段），这类问题在文献中较少见，且通常只涉及第一特征值。本文将其推广到了任意 $k$。
2. 几何与谱的联系：文章深刻揭示了凸域的几何形状（特别是其“扁平度”或各向异性）如何定量地控制 Neumann 特征值偏离理论上限的程度。
3. 方法论创新：结合了凸几何（John 椭球、Brunn-Minkowski）、谱理论（Sturm-Liouville 问题）和变分法，为处理此类“退化”极值问题提供了新的分析框架。
4. 应用潜力：这些定量不等式在形状优化、材料科学（如各向异性材料的热传导或振动分析）以及数值分析中（误差估计）具有潜在的应用价值。

综上所述，该论文通过严谨的数学分析，成功地将经典的 Kröger 不等式提升为定量的稳定性估计，不仅证明了误差项的存在和形式，还给出了具体情形下的显式常数，是该领域的重要进展。