We derive two distinct asymptotic expansions for the zeros $j_{ν,k}^{(n)}$ of the $n$-th derivative of Bessel function $J_ν^{(n)}(x)$. The first is a McMahon-type expansion for the case when $k \to \infty$ with fixed $ν$, for which we also establish an explicit error bound. The second addresses the case when $ν\to \infty$ with fixed $k$ and it involves the zeros of Airy functions and their derivatives. These results extend and refine the classical work of Wong, Lang, and Olver on the zeros of Bessel functions. In the course of obtaining our main results, we also generalize several auxiliary results, which in turn provide a broader framework for the study of zeros of special functions.
論文ID : 2510.12353タイトル : Asymptotic behavior of zeros of Bessel function derivatives著者 : Árpád Baricz, Pranav Kumar, Saminathan Ponnusamy分類 : math.CA (古典解析と常微分方程式)発表日時 : 2025年10月14日論文リンク : https://arxiv.org/abs/2510.12353v1 本論文は、ベッセル関数のn階導関数 J ν ( n ) ( x ) J_\nu^{(n)}(x) J ν ( n ) ( x ) の零点 j ν , k ( n ) j_{\nu,k}^{(n)} j ν , k ( n ) の2つの異なる漸近展開式を導出している。第1の展開式は k → ∞ k \to \infty k → ∞ かつ ν \nu ν が固定の場合のマクマホン型展開であり、明示的な誤差界を確立している。第2の展開式は ν → ∞ \nu \to \infty ν → ∞ かつ k k k が固定の場合を扱い、エアリー関数およびその導関数の零点を含む。これらの結果は、ベッセル関数の零点に関するウォン、ラング、オルバーの古典的業績を拡張・改善している。主要な結果を得る過程で、複数の補助的結果も一般化し、特殊関数の零点研究に対してより広い枠組みを提供している。
ベッセル関数の零点は、応用数学および数学物理学の様々な問題において重要な役割を果たしており、特に量子力学、散乱理論、波動伝播などの関連分野で活躍している。その広範な適用性のため、ベッセル関数およびその零点の漸近挙動を理解することは相当な注目を集めている。
古典的マクマホン展開の制限 : 初期のマクマホンによるベッセル関数 J ν ( x ) J_\nu(x) J ν ( x ) の零点に関する研究は、主に関数そのものを対象としており、その導関数ではなかった。ウォンとラングの業績 : 彼らはこれらの考え方を拡張して J ν ′ ′ ( x ) J''_\nu(x) J ν ′′ ( x ) の零点の漸近展開を得たが、二階導関数に限定されていた。研究上の空白 : n階ベッセル関数導関数の零点の性質を理解する上で大きな進展があったにもかかわらず、これらの零点の漸近展開は文献ではまだ探索されていない。本論文は、マクマホン型漸近展開を提示し、ν → ∞ \nu \to \infty ν → ∞ かつnが固定の場合を分析することで、この空白を埋めることを目指している。研究は、マクマホンが導入した標準的方法を採用して J ν ( n ) ( x ) J_\nu^{(n)}(x) J ν ( n ) ( x ) の漸近展開を導出し、k → ∞ k \to \infty k → ∞ のときの j ν , k ( n ) j_{\nu,k}^{(n)} j ν , k ( n ) の漸近挙動を確立している。
マクマホン型漸近展開の確立 : ベッセル関数n階導関数の零点に対して、k → ∞ k \to \infty k → ∞ かつ ν \nu ν が固定の場合の漸近展開式。明示的誤差界の提供 : マクマホン型展開に対して明確な誤差推定を確立し、近似の精度を定量化している。大パラメータ漸近展開の導出 : ν → ∞ \nu \to \infty ν → ∞ かつ k k k が固定の場合を扱い、エアリー関数およびその導関数の零点を含む。補助的結果の一般化 : 主要な結果を得る過程で、複数の基礎理論的結果を一般化し、特殊関数の零点研究に対してより広い枠組みを提供している。古典理論の拡張 : ウォン、ラング、オルバーの古典的業績をベッセル関数そのものから任意階の導関数へ拡張している。第1種ベッセル関数 J ν ( x ) J_\nu(x) J ν ( x ) のn階導関数 J ν ( n ) ( x ) J_\nu^{(n)}(x) J ν ( n ) ( x ) の正の零点 j ν , k ( n ) j_{\nu,k}^{(n)} j ν , k ( n ) の漸近挙動を研究する。ここで:
ν > 0 \nu > 0 ν > 0 はベッセル関数の位数n ≥ 1 n \geq 1 n ≥ 1 は導関数の位数k k k は零点の指標定理2.1 : n ∈ N 0 n \in \mathbb{N}_0 n ∈ N 0 および大きな x x x に対して、ベッセル関数の2n階導関数は以下のように表現できる:
1 2 π x J ν ( 2 n ) ( x ) = cos ( x − ν π 2 − π 4 ) τ ν ( 2 n ) ( x ) + sin ( x − ν π 2 − π 4 ) θ ν ( 2 n ) ( x ) \sqrt{\frac{1}{2}\pi x} J_\nu^{(2n)}(x) = \cos\left(x - \frac{\nu\pi}{2} - \frac{\pi}{4}\right)\tau_\nu^{(2n)}(x) + \sin\left(x - \frac{\nu\pi}{2} - \frac{\pi}{4}\right)\theta_\nu^{(2n)}(x) 2 1 π x J ν ( 2 n ) ( x ) = cos ( x − 2 ν π − 4 π ) τ ν ( 2 n ) ( x ) + sin ( x − 2 ν π − 4 π ) θ ν ( 2 n ) ( x )
ここで:
τ ν ( 2 n ) ( x ) = ∑ m = 0 ∞ α 2 m , 2 n x 2 m , θ ν ( 2 n ) ( x ) = ∑ m = 0 ∞ α 2 m + 1 , 2 n x 2 m + 1 \tau_\nu^{(2n)}(x) = \sum_{m=0}^{\infty} \frac{\alpha_{2m,2n}}{x^{2m}}, \quad \theta_\nu^{(2n)}(x) = \sum_{m=0}^{\infty} \frac{\alpha_{2m+1,2n}}{x^{2m+1}} τ ν ( 2 n ) ( x ) = ∑ m = 0 ∞ x 2 m α 2 m , 2 n , θ ν ( 2 n ) ( x ) = ∑ m = 0 ∞ x 2 m + 1 α 2 m + 1 , 2 n
係数は以下の漸化式を満たす:
α 2 m , 2 n = α 2 m , 2 n − 1 − 4 m − 1 2 α 2 m − 1 , 2 n − 1 \alpha_{2m,2n} = \alpha_{2m,2n-1} - \frac{4m-1}{2}\alpha_{2m-1,2n-1} α 2 m , 2 n = α 2 m , 2 n − 1 − 2 4 m − 1 α 2 m − 1 , 2 n − 1 α 2 m + 1 , 2 n = − α 2 m − 1 , 2 n − 1 − 4 m + 1 2 α 2 m , 2 n − 1 \alpha_{2m+1,2n} = -\alpha_{2m-1,2n-1} - \frac{4m+1}{2}\alpha_{2m,2n-1} α 2 m + 1 , 2 n = − α 2 m − 1 , 2 n − 1 − 2 4 m + 1 α 2 m , 2 n − 1
反転技法により、零点のマクマホン型展開を得る:
偶数階導関数 (J ν ( 2 m ) J_\nu^{(2m)} J ν ( 2 m ) の零点):
j ν , k + 1 ( 2 m ) = α + α 1 , 2 m α 0 , 2 m ⋅ 1 α + O ( 1 α 3 ) j_{\nu,k+1}^{(2m)} = \alpha + \frac{\alpha_{1,2m}}{\alpha_{0,2m}} \cdot \frac{1}{\alpha} + O\left(\frac{1}{\alpha^3}\right) j ν , k + 1 ( 2 m ) = α + α 0 , 2 m α 1 , 2 m ⋅ α 1 + O ( α 3 1 )
ここで α = k π + ν π 2 − π 4 \alpha = k\pi + \frac{\nu\pi}{2} - \frac{\pi}{4} α = kπ + 2 ν π − 4 π
奇数階導関数 (J ν ( 2 m + 1 ) J_\nu^{(2m+1)} J ν ( 2 m + 1 ) の零点):
j ν , k ( 2 m + 1 ) = β − α 1 , 2 m + 1 α 0 , 2 m + 1 ⋅ 1 β + O ( 1 β 3 ) j_{\nu,k}^{(2m+1)} = \beta - \frac{\alpha_{1,2m+1}}{\alpha_{0,2m+1}} \cdot \frac{1}{\beta} + O\left(\frac{1}{\beta^3}\right) j ν , k ( 2 m + 1 ) = β − α 0 , 2 m + 1 α 1 , 2 m + 1 ⋅ β 1 + O ( β 3 1 )
ここで β = k π + ν π 2 − 3 π 4 \beta = k\pi + \frac{\nu\pi}{2} - \frac{3\pi}{4} β = kπ + 2 ν π − 4 3 π
定理2.2 : 漸近展開の誤差項に対して、以下の有界推定が成立する:
2n階導関数に対して、ν ≥ − 2 n + 1 2 \nu \geq -2n + \frac{1}{2} ν ≥ − 2 n + 2 1 、n ∈ N n \in \mathbb{N} n ∈ N かつ x x x が十分に大きいとき:
∣ δ 2 n ( ν , x ) ∣ ≤ 4 ( ν + 2 n ) 2 − 1 4 x exp { 4 ( ν + 2 n ) 2 − 1 4 x } |\delta_{2n}(\nu,x)| \leq \frac{4(\nu+2n)^2-1}{4x} \exp\left\{\frac{4(\nu+2n)^2-1}{4x}\right\} ∣ δ 2 n ( ν , x ) ∣ ≤ 4 x 4 ( ν + 2 n ) 2 − 1 exp { 4 x 4 ( ν + 2 n ) 2 − 1 }
ν → ∞ \nu \to \infty ν → ∞ かつ k k k が固定のとき、オルバーの一様漸近展開理論とエアリー関数を利用する:
偶数階導関数に対して:
j ν , k ( 2 n ) = ν − a k 2 1 / 3 ν 1 / 3 + 2 1 / 3 a k 2 10 ν − 1 / 3 + O ( ν − 1 ) j_{\nu,k}^{(2n)} = \nu - \frac{a_k}{2^{1/3}}\nu^{1/3} + \frac{2^{1/3}a_k^2}{10}\nu^{-1/3} + O(\nu^{-1}) j ν , k ( 2 n ) = ν − 2 1/3 a k ν 1/3 + 10 2 1/3 a k 2 ν − 1/3 + O ( ν − 1 )
奇数階導関数に対して:
j ν , k ( 2 n + 1 ) = ν − a k ′ 2 1 / 3 ν 1 / 3 + 2 1 / 3 ( a k ′ ) 2 10 ν − 1 / 3 + O ( ν − 1 ) j_{\nu,k}^{(2n+1)} = \nu - \frac{a'_k}{2^{1/3}}\nu^{1/3} + \frac{2^{1/3}(a'_k)^2}{10}\nu^{-1/3} + O(\nu^{-1}) j ν , k ( 2 n + 1 ) = ν − 2 1/3 a k ′ ν 1/3 + 10 2 1/3 ( a k ′ ) 2 ν − 1/3 + O ( ν − 1 )
ここで a k a_k a k と a k ′ a'_k a k ′ はそれぞれエアリー関数 Ai ( x ) \text{Ai}(x) Ai ( x ) と Ai ′ ( x ) \text{Ai}'(x) Ai ′ ( x ) の第k番目の負の零点である。
統一的漸化枠組み : すべての階の導関数係数に対する統一的な漸化関係を確立し、各階の導関数の個別計算を回避している。改善された誤差分析 : ヘスコートの方法を修正することにより、古典的結果より精密な誤差界を得ている。適用範囲の拡張 : 古典理論を特定の低階導関数から任意階導関数へ拡張している。一様漸近理論の応用 : オルバーの一様漸近展開理論を高階導関数の場合に巧妙に適用している。純粋な数学理論研究であるため、「実験」は主に理論検証と数値検証を指す:
数学的帰納法による検証 : 漸化関係の正確性を数学的帰納法で証明漸近解析 : 極限解析により漸近展開の有効性を検証誤差界検証 : 具体的計算により誤差界の厳密性を検証マクマホンの原始的結果(J ν ( x ) J_\nu(x) J ν ( x ) の零点に対して) ウォンとラングの結果(J ν ′ ′ ( x ) J''_\nu(x) J ν ′′ ( x ) の零点に対して) オルバーの一様漸近展開理論 マクマホン型展開の確立 : 任意階ベッセル関数導関数の零点のマクマホン型漸近展開を成功裏に導出し、明確な係数表現を有する。誤差界の改善 : 既存の結果と比較して、新しい誤差界はより精密であり、特に高階導関数の場合に顕著である。大パラメータ展開の完全性 : ν → ∞ \nu \to \infty ν → ∞ の場合に対して、主要項と補正項を含む完全な漸近級数を提供している。三階および四階導関数の零点に対して、以下を得る:
j ν , k ( 3 ) = ν − a k 2 1 / 3 ν 1 / 3 + 2 1 / 3 a k 2 10 ν − 1 / 3 + a k 3 700 ν − 1 + ⋯ j_{\nu,k}^{(3)} = \nu - \frac{a_k}{2^{1/3}}\nu^{1/3} + \frac{2^{1/3}a_k^2}{10}\nu^{-1/3} + \frac{a_k^3}{700}\nu^{-1} + \cdots j ν , k ( 3 ) = ν − 2 1/3 a k ν 1/3 + 10 2 1/3 a k 2 ν − 1/3 + 700 a k 3 ν − 1 + ⋯
j ν , k ( 4 ) = ν − a k ′ 2 1 / 3 ν 1 / 3 + 2 1 / 3 ( a k ′ ) 2 10 ν − 1 / 3 + ( a k ′ ) 3 700 ν − 1 + ⋯ j_{\nu,k}^{(4)} = \nu - \frac{a'_k}{2^{1/3}}\nu^{1/3} + \frac{2^{1/3}(a'_k)^2}{10}\nu^{-1/3} + \frac{(a'_k)^3}{700}\nu^{-1} + \cdots j ν , k ( 4 ) = ν − 2 1/3 a k ′ ν 1/3 + 10 2 1/3 ( a k ′ ) 2 ν − 1/3 + 700 ( a k ′ ) 3 ν − 1 + ⋯
マクマホン型展開に対して、条件を満たすとき:
偶数階導関数:∣ j ν , k + 1 ( 2 n ) − ( k + ν 2 − 1 4 ) π ∣ ≤ 0.26 ( 4 ( ν + 2 n ) 2 − 1 ) π ( k + ν 2 − 0.314 ) |j_{\nu,k+1}^{(2n)} - (k + \frac{\nu}{2} - \frac{1}{4})\pi| \leq \frac{0.26(4(\nu+2n)^2-1)}{\pi(k + \frac{\nu}{2} - 0.314)} ∣ j ν , k + 1 ( 2 n ) − ( k + 2 ν − 4 1 ) π ∣ ≤ π ( k + 2 ν − 0.314 ) 0.26 ( 4 ( ν + 2 n ) 2 − 1 ) 奇数階導関数:∣ j ν , k ( 2 n + 1 ) − ( k + ν 2 − 1 4 ) π ∣ ≤ 0.26 ( 4 ( ν + 2 n + 1 ) 2 − 1 ) π ( k + ν 2 − 0.314 ) |j_{\nu,k}^{(2n+1)} - (k + \frac{\nu}{2} - \frac{1}{4})\pi| \leq \frac{0.26(4(\nu+2n+1)^2-1)}{\pi(k + \frac{\nu}{2} - 0.314)} ∣ j ν , k ( 2 n + 1 ) − ( k + 2 ν − 4 1 ) π ∣ ≤ π ( k + 2 ν − 0.314 ) 0.26 ( 4 ( ν + 2 n + 1 ) 2 − 1 ) マクマホン (1895) : ベッセル関数零点の古典的漸近展開を確立ワトソン (1944) : 理論的基礎を完善し、より厳密な証明を提供オルバー (1954) : 一様漸近展開理論を発展させたウォン・ラング (1990-1991) : 理論をベッセル関数の二階導関数に拡張バリッツ等 (2018) : n階導関数零点の実性と正性を証明本論文はこの発展の流れの中で、理論を特定の低階導関数の場合から完全に任意階導関数へ一般化し、より精密な誤差分析を提供することで貢献している。これは当該分野における重要な理論的完善である。
完全性 : ベッセル関数任意階導関数の零点の完全な漸近理論を確立した精密性 : 既存の結果より精密な誤差界を提供した統一性 : すべての階の導関数に適用可能な統一的理論枠組みを提供した空白の補完 : ベッセル関数導関数零点の漸近理論の最後のピースを完成させた方法論的革新 : 高階導関数の漸近問題を扱うための一般的方法を提供した応用の可能性 : 関連する物理および工学問題に対して理論的ツールを提供した適用条件 : 特定の結果は特定のパラメータ範囲条件を満たす必要がある計算の複雑性 : 高階項の計算は相当に複雑になる可能性がある数値検証 : 理論的結果は実用性を確認するためにさらなる数値検証が必要である修正ベッセル関数への推広 : 類似の方法を修正ベッセル関数およびその導関数に適用する他の特殊関数 : ライト関数、超ベッセル関数などへの拡張数値アルゴリズム : 理論的結果に基づいた効率的な零点計算アルゴリズムの開発理論的完全性 : 当該問題に対する完全な理論的解決策を提供している方法の厳密性 : 成熟した数学解析方法を使用し、証明過程は厳密である結果の一般性 : 特定の問題を解決するだけでなく、一般的な理論枠組みを提供している誤差分析 : 明確な誤差界を提供し、結果の実用性を高めている数値検証の欠如 : 純粋な理論研究として、十分な数値実験による検証が不足している応用例の不足 : 具体的な応用場面と実例の提供が不足している計算複雑性 : 非常に高階の導関数に対して、計算が実用的でなくなる可能性がある学術的影響 : 特殊関数理論の重要な分野を完善している実用的価値 : ベッセル関数導関数の零点を必要とする応用に対して理論的基礎を提供している方法論的価値 : 提供される技術的方法は他の類似問題に適用可能である数学物理 : 波動方程式、熱伝導方程式などの問題の解析解工学応用 : 信号処理、振動解析などの分野数値計算 : 高精度零点計算アルゴリズムの理論的基礎論文は当該分野の重要な文献を引用しており、以下を含む:
マクマホンの開拓的業績 (1895) ワトソンの古典的教科書 (1944) オルバーの漸近理論 (1954, 1962, 1974) ウォン・ラングの拡張業績 (1990, 1991) 最近の関連研究 (バリッツ等, 2018; ディミトロフ・ルン, 2025) これらの引用は研究の歴史的連続性と理論的完全性を体現している。