[論文レビュー] Locating the zeros of partial sums of exp(z) with Riemann-Hilbert methods

Thomas Kriecherbauer, Arno B. J. Kuijlaars+2 more|arXiv (Cornell University)|Sep 8, 2007

Mathematical functions and polynomials参考文献 8被引用数 3

ひとこと要約

本稿では、リーマン・ヒルベルト法を用いて、$ e^z $ の部分和の零点の均一漸近展開を構築し、誤差項を明示的に得る。特に、臨界点 $ z = 1 $ の近傍で、補完誤差関数の零点を用いた均一展開が得られ、古典的手法に比べて収縮する近傍内での精度が向上する。

ABSTRACT

In this paper we derive uniform asymptotic expansions for the partial sums of the exponential series. We indicate how this information will be used in a later publication to obtain full and explicitly computable asymptotic expansions with error bounds for all zeros of the Taylor polynomials $p_{n-1}(z) = \sum_{k=0}^{n-1} z^k/ k!$. Our proof is based on a representation of $p_{n-1}(nz)$ in terms of an integral of the form $\int_γ \frac{e^{nϕ(s)}}{s-z}ds$. We demonstrate how to derive uniform expansions for such integrals using a Riemann-Hilbert approach. A comparison with classical steepest descent analysis shows the advantages of the Riemann-Hilbert analysis in particular for points $z$ that are close to the critical points of $ϕ$.

研究の動機と目的

• Taylor多項式 $ p_{n-1}(z) = \frac{1}{k!} \sum_{k=0}^{n-1} z^k $ の零点に対する均一漸近展開を、明示的な誤差項を伴って導出すること。
• 臨界点 $ z = 1 $ の近傍で漸近的挙動が変化するため、古典的勾配降下法の限界を克服すること。
• 複素平面上の臨界点近傍を含めた一様漸近的挙動を扱うにあたり、リーマン・ヒルベルト法の利点を示すこと。
• 後続の明示的零点計算の基盤を提供することを目的とし、パラメトリクス構成と漸近的マッチングを用いる。
• 臨界点 $ z = 1 $ 近傍の零点が、補完誤差関数 $ \text{erfc}(z) $ の零点を用いて漸近的に展開可能であることを示すこと。

提案手法

• 部分和のスケーリング版 $ p_{n-1}(nz) $ を、$ \phi(s) $ が $ s = 1 $ で臨界点を持つ位相関数を用いた contour 積分 $ \int_{\gamma} \frac{e^{n\phi(s)}}{s - z} ds $ で表現する。
• 積分表現の解に対するリーマン・ヒルベルト問題を定式化し、$ \Sigma $ に沿ったジャンプ行列を定義し、Deift-Zhou の非線形勾配降下法を用いる。
• 臨界点 $ z = 1 $ の近傍に、補完誤差関数を用いて局所的パラメトリクスを構築し、グローバルモデルとマッチさせ、一様な誤差項を保証する。
• 変換 $ \zeta = w_k / \sqrt{n} $ を用いて、$ z = 1 $ 近傍の局所的挙動を $ \text{erfc}(\zeta) $ の零点に写像し、$ w_k $ を用いた漸近展開を可能にする。
• 近似式 $ A_n(z) = \tilde{\phi}(z) - \frac{1}{n} \ln P_n(z) $ を導出し、$ P_n(z) $ が局所的パラメトリクスであることを示し、$ A_n(\alpha_{k,n}) = -\frac{2\pi i k}{n} $ が成り立つことを示す。
• Deift と Zhou の局所化技術により、複素平面全体（特に $ z = 1 $ 近傍を含む）で一様な誤差項を確立する。

実験結果

リサーチクエスチョン

• RQ1 リーマン・ヒルベルト法を用いて、特に臨界点 $ z = 1 $ 近傍で、$ p_{n-1}(z) $ の零点に対する均一漸近展開をどのように導出できるか？
• RQ2 位相関数に臨界点を持つ積分を解析するにあたり、古典的勾配降下法に比べてリーマン・ヒルベルト法が持つ利点は何か？
• RQ3 $ p_{n-1}(nz) $ の零点が、$ z = 1 $ 近傍で補完誤差関数 $ \text{erfc}(z) $ の零点を用いて漸近的に展開可能であるか、明示的な誤差項を伴って示せるか？
• RQ4 収縮する $ z = 1 $ 近傍において、Szegő曲線に基づく展開と補完誤差関数の零点に基づく展開の精度を比較すると、どちらが優れているか？
• RQ5 局所的パラメトリクスの正確な構造は何か？ また、それが複素平面の異なる領域にまたがる漸近的マッチングを可能にする仕組みは何か？

主な発見

• 本稿では、$ p_{n-1}(nz) $ の零点 $ z_{k,n} $ について、補完誤差関数の零点 $ w_k $ を用いた均一漸近展開を確立した。$ 1 < k < n^\beta $ かつ $ 0 < \beta < 1 $ の範囲で有効であり、誤差項は $ O\left(\left(\frac{k}{n}\right)^{r/2}\right) $ である。
• 展開式は $ z_{k,n} = 1 + \sum_{j=1}^{r-1} \frac{q_j(\sqrt{2}w_k)}{n^{j/2}} + O\left(\left(\frac{k}{n}\right)^{r/2}\right) $ で与えられ、$ q_j $ は次数 $ j $ の明示的な多項式であり、$ q_1(x) = x $、$ q_2(x) = \frac{x^2 - 1}{3} $、$ q_3(x) = \frac{x^3 - 7x}{36} $ である。
• 補完誤差関数の零点に基づく近似は、$ O(n^{-1/3}) $ のサイズの近傍でより高い精度を示すが、Szegő曲線に基づく展開はこの範囲外でより有効である。
• 本手法により、複素平面全体（特に $ z = 1 $ 近傍を含む）で一様な誤差項が保証され、古典的勾配降下法が臨界点近傍で示す限界を克服する。
• 補完誤差関数を用いた局所的パラメトリクスの構築は、$ z = 1 $ 近傍での正しい漸近的挙動を捉えるために不可欠であり、恒等式 $ v(\zeta) = 2h(i\zeta)e^{-\zeta^2} $ によりパラメトリクスと補完誤差関数が結びつけられる。
• 結果として、距離 $ |z_{k,n} - 1| $ が $ \sqrt{k/n} $ のオーダーであることが確認され、展開が半径 $ O(n^{(\beta-1)/2}) $ の収縮する円板内で一様に成立することが示された。