[论文解读] A New Effective Asymptotic Formula for the Stieltjes Constants

Lazhar Fekih‐Ahmed|arXiv (Cornell University)|Jul 21, 2014

Analytic Number Theory Research参考文献 13被引用 4

一句话总结

本文提出了一种基于积分的新型渐近公式，用于计算Stieltjes常数γₙ，该方法通过函数方程法与鞍点分析推导得出。该方法能准确捕捉符号振荡与主导阶增长行为，且在n=137等数值病态点上表现优于以往公式——通过三重渐近展开，使n ≥ 8816时的误差降低至5%以下。

ABSTRACT

We derive a new integral formula for the Stieltjes constants. The new formula permits easy computations as well as an exact approximate asymptotic formula. Both the sign oscillations and the leading order of growth are provided. The formula can also be easily extended to generalized Euler constants.

研究动机与目标

• 推导Stieltjes常数γₙ的新型有效渐近公式，γₙ是黎曼ζ函数在s=1处洛朗展开的系数。
• 提供对γₙ的精确近似，准确反映其符号振荡与主导阶增长行为。
• 将该方法推广至与Hurwitz zeta函数ζ(s,a)相关的广义Euler常数γₙ(a)。
• 改进现有渐近公式，尤其在处理数值病态情形（如n=137时出现符号错误）方面表现更优。
• 通过高精度计算与已知结果对比，验证新公式的鲁棒性，覆盖广泛的n值范围。

提出的方法

• 通过涉及t/(eᵗ−1)的生成函数的二阶导数μ(t)，推导出s(s−1)ζ(s)的新型积分表示，适用于Re(s) > −1。
• 应用黎曼ζ函数的函数方程与复指数恒等式，将积分重写为包含复共轭指数形式的表达式，其中a = log(2π) + iπ/2。
• 在s=1处对所得表达式进行泰勒级数展开，提取系数μₙ，再通过幂级数求逆将μₙ与Stieltjes常数γₙ关联。
• 利用截断指数级数与余项估计，将γₙ表示为包含Re{(log t − a)ⁿ⁺² / (n+2)!}及高阶项的积分形式。
• 对积分I(n) = ∫₀^∞ μ(t)(log t − a)ⁿ dt应用鞍点近似，导出γₙ的渐近展开。
• 采用三重渐近公式（M=3）显著降低相对误差，尤其在n=137等关键点上，这些点在其他方法中常出现符号错误。

实验结果

研究问题

• RQ1 能否推导出一种基于积分的新型渐近公式，准确捕捉Stieltjes常数γₙ的符号振荡与主导阶增长行为？
• RQ2 所提出的三重渐近公式在精度与稳定性方面，相较于Knessl-Coffey公式与单重近似公式，尤其在n=137时表现如何？
• RQ3 为何Knessl-Coffey公式与新公式均无法正确预测γ₁₃₇的符号？这反映了该问题的何种数值病态特性？
• RQ4 该方法能否推广至计算Hurwitz zeta函数ζ(s,a)的广义Euler常数γₙ(a)？
• RQ5 该新渐近公式的收敛行为与误差分布如何，尤其在n > 8000的大n值范围内？

主要发现

• 新型三重渐近公式（定理3.2，M=3）在所有n ≥ 8816时将相对误差降至5%以下，仅在n=71158与n=84589两处例外，误差分别为−2.4%与4.5%。
• 当n ≥ 8816时，新公式的相对误差始终低于1.6%，表现出高度的鲁棒性与精度。
• 该公式在n=137时正确预测了γₙ的符号，而Knessl-Coffey公式与定理3.1的单重近似公式均失败，表明其数值稳定性显著提升。
• 利用新公式计算得γ₁₃₇ ≈ −0.00079×10²⁹，而Knessl-Coffey公式则出现符号错误。
• 该方法揭示，n=137为数值病态点，原因在于鞍点方程涉及W(ni/(2π))，导致微小扰动（如n=137.017与n=137.018）即引发符号快速变化。
• 尽管未发现公式间的直接关联，新公式在符号准确性与误差降低方面仍显著优于Knessl-Coffey公式与单重渐近公式。