[論文レビュー] Indefinite Integrals of Spherical Bessel Functions

Jolyon Bloomfield, Stephen H. P. Face+1 more|arXiv (Cornell University)|Mar 6, 2017

Mathematical functions and polynomials参考文献 4被引用数 7

ひとこと要約

本稿では、整数階数 $l$ の第一種ハンケル関数 $j_l(x)$ を含む単項式の積の不定積分について、閉形式の解析解および新規の再帰関係を提示する。この手法により、特に初期宇宙宇宙論において一般的に見られる高周波数振動を示す積分を、初等関数または終了する再帰鎖への還元によって、正確かつ効率的に数値的に評価可能となる。

ABSTRACT

Highly oscillatory integrals, such as those involving Bessel functions, are best evaluated analytically as much as possible, as numerical errors can be difficult to control. We investigate indefinite integrals involving monomials in $x$ multiplying one or two spherical Bessel functions of the first kind $j_l(x)$ with integer order $l$. Closed-form solutions are presented where possible, and recursion relations are developed that are guaranteed to reduce all integrals in this class to closed-form solutions. These results allow for definite integrals over spherical Bessel functions to be computed quickly and accurately. For completeness, we also present our results in terms of ordinary Bessel functions, but in general, the recursion relations do not terminate.

研究の動機と目的

• 整数 $n$, $k,l \geq 0$ および実数 $\alpha, \beta$ を用いた $\int dx\,x^n j_l(\alpha x)$ および $\int dx\,x^n j_k(\alpha x)j_l(\beta x)$ の形の不定積分について、閉形式の解析解を導出すること。
• 高階数の球ベッセル関数を低階数の項に体系的に還元する再帰関係を構築し、既知の閉形式解に到達するように終了を保証すること。
• 一般化超幾何関数を用いず、広く利用可能な特殊関数のみを用いて、標準的な科学計算環境で実装可能な結果を提供すること。
• 特に定積分に対しては十分に文書化されている一方で、不定積分に関しては包括的な解析的ツールが不足している点に注目し、その欠落を補うこと。
• 初期宇宙における量子場の球面調和モードの共分散計算に生じる定積分を、高精度で評価可能にするための基盤を提供すること。

提案手法

• 球ベッセル関数の標準的再帰恒等式、特に $j_l(x) = \frac{l-1}{x}j_{l-1}(x) - \partial_x j_{l-1}(x)$ および $j_l(x) = \frac{2l-1}{x}j_{l-1}(x) - j_{l-2}(x)$ をもとに再帰関係を導出する。
• 部分積分および関数的恒等式を用いて、$\int dx\,x^n j_l(\alpha x)$ の形の積分を、低階数の球ベッセル関数の組み合わせに還元する。
• 積 $j_k(\alpha x)j_l(\beta x)$ の再帰的還元スキームを構築し、対称性と再帰を活用して高階数の積を低階数のものに表現する。
• 終了条件を確立する：すべての再帰鎖は、$j_0(x)$ の sinc 関数としての性質に起因し、解析的に扱える既知の閉形式積分 $\int dx\,x^n j_0(\alpha x)$ で終了する。
• 通常ベッセル関数 $J_l(x)$ に対しても同等の形に変換するが、これらの再帰関係は一般には終了しないことに注意する。
• すべての式が初等関数および科学計算ライブラリ（C や Python におけるものなど）で広く利用可能な標準的特殊関数（例：$J_l(x)$, $\mathrm{sinc}(x)$）で構成されることを保証する。

実験結果

リサーチクエスチョン

• RQ1 単項式と一つまたは二つの球ベッセル関数の積の不定積分は、閉形式で表されるか、または終了する再帰関係で表現可能か？
• RQ2 高階数の球ベッセル関数積分を、解析的に解ける基本ケースに体系的に還元する再帰関係はどのように導出可能か？
• RQ3 一般化超幾何関数に依存せずに、数値コードで効率的かつ高精度に評価可能な方法は何か？
• RQ4 なぜ球ベッセル関数の再帰関係は終了するが、通常ベッセル関数のそれらは終了しないのか？
• RQ5 球ベッセル関数の構造と、振動的積分における解析的可積分性との関係は何か？

主な発見

• 本稿では、$\int dx\,x^n j_l(\alpha x)$ および $\int dx\,x^n j_k(\alpha x)j_l(\beta x)$ の新しい閉形式解を導出し、低階数のケースについて明示的な式を提示している。
• すべての積分が基本ケースに還元される体系的な再帰フレームワークが構築され、$j_l(x)$ の整数階数性に起因して再帰鎖の終了が保証されている。
• 再帰関係は数値実装に最適化されており、科学計算ライブラリで広く利用可能な標準的特殊関数のみを用いている。
• $j_0(x) = \mathrm{sinc}(x)$ の場合、積分 $\int dx\,x^n j_0(\alpha x)$ は解析的に扱えるため、再帰鎖の終端ケースとして機能する。
• 結果は通常ベッセル関数 $J_l(x)$ に拡張可能であるが、対応する再帰関係は一般には終了せず、球ベッセル関数と通常ベッセル関数の間の重要な相違を浮き彫りにしている。
• 本手法により、初期宇宙における量子場の球面調和モードの共分散に生じる高周波数振動を示す積分を、高精度で評価可能となった。