エアリー関数
物理科学におけるエアリー函数（エアリーかんすう、英: Airy function）あるいは第一種エアリー函数 (Airy function of the first kind) Ai(x) は、イギリスの天文学者ジョージ・ビドル・エアリー (1801–92) に因んで名づけられた特殊函数である。この函数 Ai(x) および第二種エアリー函数とも呼ばれる関連の函数（A を次の文字 B に変えて、故に冗談めかしてベアリー (Biry) 函数とも）Bi(x) は、エアリー方程式あるいはストークス方程式と呼ばれる微分方程式
- [math]\frac{d^2y}{dx^2} - xy = 0[/math]
の線型独立な解としても言及される。これは転回点（turning point: 方程式の解が振動型から指数型へ変わる特徴点）を持つ最も単純な二階線型微分方程式である。
エアリー函数は三角ポテンシャル井戸に留め置かれた粒子に対する、あるいは一次元定力場における粒子に対するシュレーディンガー方程式の解である。同じ理由により、ポテンシャルが位置の線型函数で局所近似されるときの、転回点の周りでのWKB近似として、エアリー函数は一様半古典近似を与えるのに利用できる。三角ポテンシャル井戸解は、多くの半導体デバイスを理解することに直接的に関係がある。
エアリー函数はまた、虹のようなテンプレート:訳語疑問点範囲。歴史的にはこれがエアリーがこの特殊函数を導入するに至った数学的問題であった。またエアリー函数は顕微鏡学や天文学においても重要である。つまり、エアリー函数は（顕微鏡や望遠鏡の解像限界よりも小さな）点光源によって与えられる回折や干渉のパターンを記述する。
定義
- [math]\begin{align} \operatorname{Ai}(x) & = \frac{1}{\pi} \int_0^\infty \cos \left( \frac{t^3}{3} + xt \right) ~ dt\\ &\equiv \frac{1}{\pi} \lim_{b \to \infty} \int_0^b \cos \left( \frac{t^3}{3} + xt \right) ~ dt \end{align}[/math]
として定義することができる。これが収束することは、激しく振動するグラフの正の成分と負の成分とが互いに打ち消し合う（これは部分積分で確認できる）ことによるものである。
函数 y = Ai(x) はエアリー方程式
- [math]y'' - xy = 0[/math]
を満足する。この方程式は二つの線型独立な解を持つ。スカラー倍の違いを除いて、Ai(x) は x → ∞ で y → 0 なる条件を満たす唯一の解である。もう一つの解として第二種エアリー函数 Bi(x) を取るのが標準的である。第二種エアリー函数は第一種エアリー函数 Ai(x) と同じ振幅を持ち x → −∞ で位相が π/2 だけ異なる解
- [math]\operatorname{Bi}(x) = \frac{1}{\pi} \int_0^\infty \left[ \exp \left( -\frac{t^3}{3} + xt \right) + \sin \left( \frac{t^3}{3} + xt \right) \right] dt[/math]
として定義することができる。
性質
Ai(x), Bi(x) およびそれらの導函数の x = 0 における値は
- [math]\begin{align} \operatorname{Ai}(0) &= \frac{1}{3^{2/3}\Gamma(2/3)}, & \operatorname{Ai}'(0) &= -\frac{1}{3^{1/3}\Gamma(1/3)}, \\ \operatorname{Bi}(0) &= \frac{1}{3^{1/6} \Gamma(2/3)}, & \operatorname{Bi}'(0) &= \frac{3^{1/6}}{\Gamma(1/3)} \end{align}[/math]
で与えられる。ここで Γ はガンマ函数である。これにより Ai(x) および Bi(x) のロンスキアンが 1/π となることが従う。
x が正のとき、Ai(x) は正値凸函数かつ指数函数的に零に減少し、Bi(x) は正値凸函数かつ零の周辺で限りなく小さい振幅を以って無限回振動する。x が負のとき、Ai(x), Bi(x) はともに零の周辺で限りなく小さくなる振幅を以って無限回振動する。これはエアリー函数に対する漸近公式（後述）によって保証される。
任意の変数に対するエアリー函数同士は、広義リーマン積分
- [math] \int_{-\infty}^\infty \operatorname{Ai}(t+x) \operatorname{Ai}(t+y)\,\mathit{dt} = \delta(x-y)[/math]
が成り立つという意味において互いに直交する[1]。
漸近公式
後述のように、エアリー函数は複素数平面上へ延長することができて整函数を与える。そのエアリー函数の偏角 arg(z) が一定値を保ったまま テンプレート:Abs を無限大へ飛ばすときの漸近挙動は arg(z) に依存して決まる（これをストークス現象と言う）。テンプレート:Abs < π のとき、Ai(z) に対して以下の漸近公式
- [math]\operatorname{Ai}(z)\sim \dfrac{e^{-\frac{2}{3}z^{\frac{3}{2}}}}{2\sqrt\pi\,z^{\frac{1}{4}}} \left[ \sum_{n=0}^{\infty} \dfrac{(-1)^n \Gamma(n+\frac{5}{6})\Gamma(n+\frac{1}{6})\left(\frac{3}{4}\right)^n}{2\pi n! z^{3n/2}} \right][/math]
を得る[2]。Bi(z) に対しても同様の漸近公式
- [math]\operatorname{Bi}(z)\sim \frac{e^{\frac{2}{3}z^{\frac{3}{2}}}}{\sqrt\pi\,z^{\frac{1}{4}}}\left[ \sum_{n=0}^{\infty} \dfrac{ \Gamma(n+\frac{5}{6})\Gamma(n+\frac{1}{6})\left(\frac{3}{4}\right)^n}{2\pi n! z^{3n/2}} \right][/math]
存在するが テンプレート:Abs < π/3 でしか適用できない。
よりきつい Ai(z) に対する公式、および π/3 < テンプレート:Abs < π のときの Bi(z) に対する公式、あるいは同じことだが Ai(−z) および テンプレート:Abs < 2π/3 だが零でないときの Bi(−z) の漸近近似が
- [math]\begin{align} \operatorname{Ai}(-z) &{}\sim \frac{\sin \left(\frac23z^{\frac{3}{2}}+\frac{\pi}{4} \right)}{\sqrt\pi\,z^{\frac{1}{4}}} \\[6pt] \operatorname{Bi}(-z) &{}\sim \frac{\cos \left(\frac23z^{\frac{3}{2}}+\frac{\pi}{4} \right)}{\sqrt\pi\,z^{\frac{1}{4}}} \end{align}[/math]
で与えられる[3]。
テンプレート:Abs = 0 のとき、これらは良い近似であるが漸近的でない（Ai(−z) または Bi(−z) と上記の近似式との比は、現れる正弦および余弦の値が零となるところで無限大に発散することによる）。これらの極限に対する漸近展開も可能である。それらについては {{#invoke:Footnotes | harvard_citation }} および {{#invoke:Footnotes | harvard_citation }} にある。
複素エアリー函数
エアリー函数の定義を
- [math]\operatorname{Ai}(z) = \frac{1}{2\pi i} \int_{C} \exp \left( \frac{t^3}{3} - zt \right) dt[/math]
と置くことにより複素数平面まで拡張することができる。ただし、積分路 C は偏角 −π/2 の無限遠点 −∞i から偏角 π/2 の無限遠点 ∞i までとる。あるいは、微分方程式 y′′ − xy = 0 によって Ai(x) および Bi(x) を複素数平面上の整函数に拡張することもできる。
Ai(x) に対する漸近公式は、xテンプレート:Exp の主値をとり、x が負の実軸の近くを除いて有界ならば、複素数平面上でもやはり有効である。Bi(x) に対する漸近公式が有効なのは、x が適当な δ > 0 に対する扇形 {x ∈ C : テンプレート:Abs < テンプレート:Fraction − δ} にあるときである。また、Ai(−x), Bi(−x) に対する公式は x が扇形 {x ∈ C : テンプレート:Abs < テンプレート:Fraction − δ} にあるとき有効である。
エアリー函数の漸近的挙動から分かることは、Ai(x), Bi(x) ともに負の実軸上に無限個の零点を持つことである。函数 Ai(x) は複素数平面上に他の零点を持たないが、Bi(x) は扇形 {z ∈ C : テンプレート:Fraction < テンプレート:Abs < テンプレート:Fraction} にも無限個の零点を持つ。
グラフ
| [math]\Re[\operatorname{Ai}(x + iy)][/math] | [math]\Im[\operatorname{Ai}(x + iy)][/math] | [math]|\operatorname{Ai}(x + iy)|[/math] | [math]\operatorname{arg}[\operatorname{Ai}(x + iy)][/math] |
|---|---|---|---|
| 200px | 200px | 200px | 200px |
| 200px | 200px | 200px | 200px |
| [math]\Re[\operatorname{Bi}(x + iy)][/math] | [math]\Im[\operatorname{Bi}(x + iy)][/math] | [math]|\operatorname{Bi}(x + iy)|[/math] | [math]\operatorname{arg}[\operatorname{Bi}(x + iy)][/math] |
|---|---|---|---|
| 200px | 200px | 200px | 200px |
| 200px | 200px | 200px | 200px |
他の特殊関数との関係
正の引数に対して、エアリー関数は変形ベッセル関数と
- [math]\begin{align} \operatorname{Ai}(x) & = \frac{1}{\pi}\sqrt{\frac{x}{3}} \, K_{1/3}(2x^{3/2}/3), \\[8pt] \operatorname{Bi}(x) & = \sqrt{\frac{x}{3}}\left(I_{1/3}(2x^{3/2}/3) + I_{-1/3}(2x^{3/2}/3)\right) \end{align}[/math]
なる関係を持つ。ここで I±1/3, K1/3 は xテンプレート:Expy" + xyテンプレート:' − (xテンプレート:Exp + テンプレート:Fraction)y = 0 の解である。エアリー函数の一階微分は
- [math]\operatorname{Ai}'(x) = - \frac{x}{\pi\sqrt{3}} \, K_{2/3}(2x^{3/2}/3)[/math]
であり、函数 K1/3 および K2/3 は急速に収束する積分によって表すことができる[4]（変形ベッセル函数の項も参照）。
負の引数に対して、エアリー関数はベッセル関数と
- [math]\begin{align} \operatorname{Ai}(-x) & = \sqrt{\frac{x}{9}} \left(J_{1/3}(2x^{3/2}/3) + J_{-1/3}(2x^{3/2}/3)\right), \\[8pt] \operatorname{Bi}(-x) & = \sqrt{\frac{x}{3}} \left(J_{-1/3}(2x^{3/2}/3) - J_{1/3}(2x^{3/2}/3)\right) \end{align}[/math]
なる関係を持つ。ここに J±1/3 は xテンプレート:Expy" + xyテンプレート:' + (xテンプレート:Exp − テンプレート:Fraction)y = 0 の解。
スコアラーの関数（y" − xy = 1/π の解）もまたエアリー函数を用いて
- [math]\begin{align} \operatorname{Gi}(x) & = \operatorname{Bi}(x) \int_x^\infty \operatorname{Ai}(t) \mathit{dt} + \operatorname{Ai}(x) \int_0^x \operatorname{Bi}(t) \mathit{dt}, \\[8pt] \operatorname{Hi}(x) & = \operatorname{Bi}(x) \int_{-\infty}^x \operatorname{Ai}(t) \mathit{dt} - \operatorname{Ai}(x) \int_{-\infty}^x \operatorname{Bi}(t) \mathit{dt} \end{align}[/math]
と書くこともできる。
フーリエ変換
エアリー函数 Ai(x) の定義から直接に、そのフーリエ変換が
- [math]\mathcal{F}(\operatorname{Ai})(k) := \int_{-\infty}^{\infty} \operatorname{Ai}(x)\ e^{- 2\pi i k x}\,dx = e^{i(2\pi k)^3/3}[/math]
で与えられることが示せる。
歴史
エアリー関数の名の由来は、イギリス人天文学者・物理学者のジョージ・ビドル・エアリー (1801–1892) である。エアリーは光学の初期研究 (Airy 1838) においてこの函数に遭遇した。"Ai(x)" という表記はハロルド・ジェフリースによるものである。エアリーはイギリス王室天文官を1835年から1881年まで務めた人物である。
関連項目
注
- ↑ David E. Aspnes, Physical Review, 147, 554 (1966)
- ↑ Abramowitz & Stegun 1970, p. 448, Eqns 10.4.59 and 10.4.63
- ↑ Abramowitz & Stegun 1970, p. 448, Eqns 10.4.60 and 10.4.64
- ↑ M.Kh.Khokonov. Cascade Processes of Energy Loss by Emission of Hard Photons // JETP, V.99, No.4, pp. 690-707 \ (2004).
参考文献
- テンプレート:AS ref
- Airy (1838), “On the intensity of light in the neighbourhood of a caustic”, Transactions of the Cambridge Philosophical Society (University Press) 6: 379–402
- Frank William John Olver (1974). Asymptotics and Special Functions, Chapter 11. Academic Press, New York.
- Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007), “Section 6.6.3. Airy Functions”, Numerical Recipes: The Art of Scientific Computing (3rd ed.), New York: Cambridge University Press, ISBN 978-0-521-88068-8
- Vallée, Olivier; Soares, Manuel (2004), Airy functions and applications to physics, London: Imperial College Press, ISBN 978-1-86094-478-9, MR 2114198
外部リンク
- テンプレート:Springer
- Weisstein, Eric W. “Airy Functions”. MathWorld（英語）. Template:Cite webの呼び出しエラー：引数 accessdate は必須です。
- Wolfram function pages for Ai and Bi functions. Includes formulas, function evaluator, and plotting calculator.
- テンプレート:Dlmf