エアリー関数

エアリー関数（エアリーかんすう、英: Airy function）あるいは第一種エアリー関数 (Airy function of the first kind) Ai(x) は、イギリスの天文学者ジョージ・ビドル・エアリー (1801–92) に因んで名づけられた特殊関数である。この関数 Ai(x) および第二種エアリー関数とも呼ばれる関連の関数（A を次の文字 B に変えて、故に冗談めかしてベアリー (Biry) 関数とも）Bi(x) は、エアリー方程式あるいはストークス方程式と呼ばれる微分方程式 d 2 y d x 2 − x y = 0 {\displaystyle {\frac {d^{2}y}{dx^{2}}}-xy=0} の線型独立な解としても言及される。これは転回点（turning point: 方程式の解が振動型から指数型へ変わる特徴点）を持つ最も単純な二階線型微分方程式である。 エアリー関数は三角ポテンシャル井戸に留め置かれた粒子に対する、あるいは一次元定力場における粒子に対するシュレーディンガー方程式の解である。同じ理由により、ポテンシャルが位置の線型関数で局所近似されるときの、転回点の周りでのWKB近似として、エアリー関数は一様半古典近似を与えるのに利用できる。三角ポテンシャル井戸解は、多くの半導体デバイスを理解することに直接的に関係がある。 エアリー関数はまた、虹のような方向性集光（英語版）の周辺強度の形でも根底にある[訳語疑問点]。歴史的にはこれがエアリーがこの特殊関数を導入するに至った数学的問題であった。またエアリー関数は顕微鏡学（英語版）や天文学においても重要である。つまり、エアリー関数は（顕微鏡や望遠鏡の解像限界よりも小さな）点光源（英語版）によって与えられる回折や干渉のパターンを記述する。 実変数 x に対する第一種エアリー関数は広義リーマン積分 Ai ⁡ ( x ) = 1 π ∫ 0 ∞ cos ⁡ ( t 3 3 + x t )   d t ≡ 1 π lim b → ∞ ∫ 0 b cos ⁡ ( t 3 3 + x t )   d t {\displaystyle {\begin{aligned}\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{aligned}}} として定義することができる。これが収束することは、激しく振動するグラフの正の成分と負の成分とが互いに打ち消し合う（これは部分積分で確認できる）ことによるものである。 関数 y = Ai(x) はエアリー方程式 y ″ − x y = 0 {\displaystyle y''-xy=0} を満足する。この方程式は二つの線型独立な解を持つ。スカラー倍の違いを除いて、Ai(x) は x → ∞ で y → 0 なる条件を満たす唯一の解である。もう一つの解として第二種エアリー関数 Bi(x) を取るのが標準的である。第二種エアリー関数は第一種エアリー関数 Ai(x) と同じ振幅を持ち x → −∞ で位相が π/2 だけ異なる解 Bi ⁡ ( x ) = 1 π ∫ 0 ∞ [ exp ⁡ ( − t 3 3 + x t ) + sin ⁡ ( t 3 3 + x t ) ] d t {\displaystyle \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} として定義することができる。 Ai(x), Bi(x) およびそれらの導関数の x = 0 における値は Ai ⁡ ( 0 ) = 1 3 2 / 3 Γ ( 2 / 3 ) , Ai ′ ⁡ ( 0 ) = − 1 3 1 / 3 Γ ( 1 / 3 ) , Bi ⁡ ( 0 ) = 1 3 1 / 6 Γ ( 2 / 3 ) , Bi ′ ⁡ ( 0 ) = 3 1 / 6 Γ ( 1 / 3 ) {\displaystyle {\begin{aligned}\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{aligned}}} で与えられる。ここで Γ はガンマ関数である。これにより Ai(x) および Bi(x) のロンスキアンが 1/π となることが従う。 x が正のとき、Ai(x) は正値凸関数かつ指数関数的に零に減少し、Bi(x) は正値凸関数かつ零の周辺で限りなく小さい振幅を以って無限回振動する。x が負のとき、Ai(x), Bi(x) はともに零の周辺で限りなく小さくなる振幅を以って無限回振動する。これはエアリー関数に対する漸近公式（後述）によって保証される。 任意の変数に対するエアリー関数同士は、広義リーマン積分 ∫ − ∞ ∞ Ai ⁡ ( t + x ) Ai ⁡ ( t + y ) d t = δ ( x − y ) {\displaystyle \int _{-\infty }^{\infty }\operatorname {Ai} (t+x)\operatorname {Ai} (t+y)\,{\mathit {dt}}=\delta (x-y)} が成り立つという意味において互いに直交する。 後述のように、エアリー関数は複素数平面上へ延長することができて整関数を与える。そのエアリー関数の偏角 arg(z) が一定値を保ったまま |z| を無限大へ飛ばすときの漸近挙動は arg(z) に依存して決まる（これをストークス現象（英語版）と言う）。|arg(z)| < π のとき、Ai(z) に対して以下の漸近公式 Ai ⁡ ( z ) ∼ e − 2 3 z 3 2 2 π z 1 4 [ ∑ n = 0 ∞ ( − 1 ) n Γ ( n + 5 6 ) Γ ( n + 1 6 ) ( 3 4 ) n 2 π n ! z 3 n / 2 ] {\displaystyle \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]} を得る。Bi(z) に対しても同様の漸近公式 Bi ⁡ ( z ) ∼ e 2 3 z 3 2 π z 1 4 [ ∑ n = 0 ∞ Γ ( n + 5 6 ) Γ ( n + 1 6 ) ( 3 4 ) n 2 π n ! z 3 n / 2 ] {\displaystyle \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]} 存在するが |arg(z)| < π/3 でしか適用できない。 π/3 < |arg(z)| < π のときのより正確な Ai(z) および Bi(z) に対する公式、あるいは同じことだが |arg(z)| < 2π/3 だが零でないときの Ai(−z) および Bi(−z) の漸近近似が A i ( − z ) ∼ sin ⁡ ( 2 3 z 3 2 + π 4 ) π z 1 4 [ ∑ n = 0 ∞ ( − 1 ) n Γ ( 2 n + 5 6 ) Γ ( 2 n + 1 6 ) ( 3 4 ) 2 n 2 π ( 2 n ) ! z 3 n ] − cos ⁡ ( 2 3 z 3 2 + π 4 ) π z 1 4 [ ∑ n = 0 ∞ ( − 1 ) n Γ ( 2 n + 11 6 ) Γ ( 2 n + 7 6 ) ( 3 4 ) 2 n + 1 2 π ( 2 n + 1 ) ! z 3 n + 3 / 2 ] B i ( − z ) ∼ cos ⁡ ( 2 3 z 3 2 + π 4 ) π z 1 4 [ ∑ n = 0 ∞ ( − 1 ) n Γ ( 2 n + 5 6 ) Γ ( 2 n + 1 6 ) ( 3 4 ) 2 n 2 π ( 2 n ) ! z 3 n ] + sin ⁡ ( 2 3 z 3 2 + π 4 ) π z 1 4 [ ∑ n = 0 ∞ ( − 1 ) n Γ ( 2 n + 11 6 ) Γ ( 2 n + 7 6 ) ( 3 4 ) 2 n + 1 2 π ( 2 n + 1 ) ! z 3 n + 3 / 2 ] {\displaystyle {\begin{aligned}\mathrm {Ai} (-z)\sim &{}{\frac {\sin \left({\frac {2}{3}}z^{\frac {3}{2}}+{\frac {\pi }{4}}\right)}{{\sqrt {\pi }}\,z^{\frac {1}{4}}}}\left[\sum _{n=0}^{\infty }{\dfrac {(-1)^{n}\Gamma (2n+{\frac {5}{6}})\Gamma (2n+{\frac {1}{6}})\left({\frac {3}{4}}\right)^{2n}}{2\pi (2n)!z^{3n}}}\right]\\[6pt]&{}-{\frac {\cos \left({\frac {2}{3}}z^{\frac {3}{2}}+{\frac {\pi }{4}}\right)}{{\sqrt {\pi }}\,z^{\frac {1}{4}}}}\left[\sum _{n=0}^{\infty }{\dfrac {(-1)^{n}\Gamma (2n+{\frac {11}{6}})\Gamma (2n+{\frac {7}{6}})\left({\frac {3}{4}}\right)^{2n+1}}{2\pi (2n+1)!z^{3n+3/2}}}\right]\\[6pt]\mathrm {Bi} (-z)\sim &{}{\frac {\cos \left({\frac {2}{3}}z^{\frac {3}{2}}+{\frac {\pi }{4}}\right)}{{\sqrt {\pi }}\,z^{\frac {1}{4}}}}\left[\sum _{n=0}^{\infty }{\dfrac {(-1)^{n}\Gamma (2n+{\frac {5}{6}})\Gamma (2n+{\frac {1}{6}})\left({\frac {3}{4}}\right)^{2n}}{2\pi (2n)!z^{3n}}}\right]\\[6pt]&{}+{\frac {\sin \left({\frac {2}{3}}z^{\frac {3}{2}}+{\frac {\pi }{4}}\right)}{{\sqrt {\pi }}\,z^{\frac {1}{4}}}}\left[\sum _{n=0}^{\infty }{\dfrac {(-1)^{n}\Gamma (2n+{\frac {11}{6}})\Gamma (2n+{\frac {7}{6}})\left({\frac {3}{4}}\right)^{2n+1}}{2\pi (2n+1)!z^{3n+3/2}}}\right]\end{aligned}}} で与えられる。 |arg(z)| = 0 のとき、これらは良い近似であるが漸近的でない（Ai(−z) または Bi(−z) と上記の近似式との比は、現れる正弦および余弦の値が零となるところで無限大に発散することによる）。これらの極限に対する漸近展開も可能である。それらについては (Abramowitz & Stegun 1954) および (Olver 1974) にある。 エアリー関数の定義を Ai ⁡ ( z ) = 1 2 π i ∫ C exp ⁡ ( t 3 3 − z t ) d t {\displaystyle \operatorname {Ai} (z)={\frac {1}{2\pi i}}\int _{C}\exp \left({\frac {t^{3}}{3}}-zt\right)dt} と置くことにより複素数平面まで拡張することができる。ただし、積分路 C は偏角 −π/2 の無限遠点 −∞i から偏角 π/2 の無限遠点 ∞i までとる。あるいは、微分方程式 y′′ − xy = 0 によって Ai(x) および Bi(x) を複素数平面上の整関数に拡張することもできる。 Ai(x) に対する漸近公式は、x2/3 の主値をとり、x が負の実軸の近くを除いて有界ならば、複素数平面上でもやはり有効である。Bi(x) に対する漸近公式が有効なのは、x が適当な δ > 0 に対する扇形 {x ∈ C : |arg(x)| < π⁄3 − δ} にあるときである。また、Ai(−x), Bi(−x) に対する公式は x が扇形 {x ∈ C : |arg(x)| < 2π⁄3 − δ} にあるとき有効である。 エアリー関数の漸近的挙動から分かることは、Ai(x), Bi(x) ともに負の実軸上に無限個の零点を持つことである。関数 Ai(x) は複素数平面上に他の零点を持たないが、Bi(x) は扇形 {z ∈ C : π⁄3 < |arg(z)| < π⁄2} にも無限個の零点を持つ。 正の引数に対して、エアリー関数は変形ベッセル関数と Ai ⁡ ( x ) = 1 π x 3 K 1 / 3 ( 2 x 3 / 2 / 3 ) , Bi ⁡ ( x ) = x 3 ( I 1 / 3 ( 2 x 3 / 2 / 3 ) + I − 1 / 3 ( 2 x 3 / 2 / 3 ) ) {\displaystyle {\begin{aligned}\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{aligned}}} なる関係を持つ。ここで I±1/3, K1/3 は x2y" + xy' − (x2 + 1⁄9)y = 0 の解である。エアリー関数の一階微分は Ai ′ ⁡ ( x ) = − x π 3 K 2 / 3 ( 2 x 3 / 2 / 3 ) {\displaystyle \operatorname {Ai} '(x)=-{\frac {x}{\pi {\sqrt {3}}}}\,K_{2/3}(2x^{3/2}/3)} であり、関数 K1/3 および K2/3 は急速に収束する積分によって表すことができる（変形ベッセル関数の項も参照）。 負の引数に対して、エアリー関数はベッセル関数と Ai ⁡ ( − x ) = x 9 ( J 1 / 3 ( 2 x 3 / 2 / 3 ) + J − 1 / 3 ( 2 x 3 / 2 / 3 ) ) , Bi ⁡ ( − x ) = x 3 ( J − 1 / 3 ( 2 x 3 / 2 / 3 ) − J 1 / 3 ( 2 x 3 / 2 / 3 ) ) {\displaystyle {\begin{aligned}\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{aligned}}} なる関係を持つ。ここに J±1/3 は x2y" + xy' + (x2 − 1⁄9)y = 0 の解。 スコアラーの関数（英語版）（y" − xy = 1/π の解）もまたエアリー関数を用いて Gi ⁡ ( x ) = Bi ⁡ ( x ) ∫ x ∞ Ai ⁡ ( t ) d t + Ai ⁡ ( x ) ∫ 0 x Bi ⁡ ( t ) d t , Hi ⁡ ( x ) = Bi ⁡ ( x ) ∫ − ∞ x Ai ⁡ ( t ) d t − Ai ⁡ ( x ) ∫ − ∞ x Bi ⁡ ( t ) d t {\displaystyle {\begin{aligned}\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{aligned}}} と書くこともできる。 エアリー関数 Ai(x) の定義から直接に、そのフーリエ変換が F ( Ai ) ( k ) := ∫ − ∞ ∞ Ai ⁡ ( x )   e − 2 π i k x d x = e i ( 2 π k ) 3 / 3 {\displaystyle {\mathcal {F}}(\operatorname {Ai} )(k):=\int _{-\infty }^{\infty }\operatorname {Ai} (x)\ e^{-2\pi ikx}\,dx=e^{i(2\pi k)^{3}/3}} で与えられることが示せる。 エアリー関数の名の由来は、イギリス人天文学者・物理学者のジョージ・ビドル・エアリー (1801–1892) である。エアリーは光学の初期研究 (Airy 1838) においてこの関数に遭遇した。"Ai(x)" という表記はハロルド・ジェフリースによるものである。エアリーはイギリス王室天文官を1835年から1881年まで務めた人物である。 ウィッテン予想: コンツェヴィッチによる証明にエアリー関数の行列値への一般化が使われている。 エアリーゼータ関数（英語版）