球贝塞尔函数

                     

贡献者: addis

预备知识 贝塞尔函数

  1球贝塞尔方程(spherical Bessel's equation)

\begin{equation} x^2 \frac{\mathrm{d}^{2}{y}}{\mathrm{d}{x}^{2}} + 2x \frac{\mathrm{d}{y}}{\mathrm{d}{x}} + [x^2 - l(l + 1)]y = 0~. \end{equation}

图
图 1:球贝塞尔函数和球汉克尔函数(来自 Wikipedia)

   两个线性无关的解分别为第一类球贝塞尔函数 $j_l(x)$ 和第二类球贝塞尔函数 $y_l(x)$,见图 1 。它们可以通过贝塞尔函数 $J_l$,$Y_l$ 来定义

\begin{equation} j_l(x) = \sqrt{\frac{\pi}{2x}} J_{l+1/2}(x)~, \qquad y_l(x) = \sqrt{\frac{\pi}{2x}} Y_{l+1/2}(x)~. \end{equation}
同样也可以定义两类球汉克尔函数(spherical Hankel's function)
\begin{equation} h_l^{(1)}(x) = \sqrt {\frac{\pi }{2x}} H_{l+1/2}^{(1)}(x) = j_l(x) + \mathrm{i} y_l(x)~, \end{equation}
\begin{equation} h_l^{(2)}(x) = \sqrt{\frac{\pi }{2x}} H_{l+1/2}^{(2)}(x) = j_l(x) - \mathrm{i} y_l(x)~. \end{equation}
另一种等效的定义方式使用 Rayleigh's 方程
\begin{equation} j_l(x) = (-x)^l \left(\frac{1}{x} \frac{\mathrm{d}}{\mathrm{d}{x}} \right) ^l \frac{\sin x}{x}~, \end{equation}
\begin{equation} y_l(x) = -(-x)^l \left(\frac{1}{x} \frac{\mathrm{d}}{\mathrm{d}{x}} \right) ^l \frac{\cos x}{x}~. \end{equation}
$l$ 为整数的球贝塞尔函数都可以用正弦余弦函数表示,例如
\begin{align} &j_0(x) = \frac{\sin x}{x}~,\\ &j_1(x) = \frac{\sin x}{x^2} - \frac{\cos x}{x}~,\\ &j_2(x) = \left(\frac{3}{x^2} - 1 \right) \frac{\sin x}{x} - \frac{3\cos x}{x^2}~,\\ &j_3(x) = \left(\frac{15}{x^3} - \frac{6}{x} \right) \frac{\sin x}{x} - \left(\frac{15}{x^2}-1 \right) \frac{\cos x}{x}~. \end{align}
\begin{align} &y_0(x) = -\frac{\cos x}{x}~,\\ &y_1(x) = -\frac{\cos x}{x^2} - \frac{\sin x}{x}~,\\ &y_2(x) = \left(-\frac{3}{x^2}+1 \right) \frac{\cos x}{x} - \frac{3\sin x}{x^2}~,\\ &y_3(x) = \left(-\frac{15}{x^3}+\frac{6}{x} \right) \frac{\cos x}{x} - \left(\frac{15}{x^2} - 1 \right) \frac{\sin x}{x}~. \end{align}
更多项可以用 Mathematica 命令如2
l = 4; Series[SphericalBesselJ[l, x], {x, \[Infinity], 1000}] // 
  Normal // Simplify

1. 性质

   原点值

\begin{equation} j_l(0) = \delta_{l,0} ~,\qquad y_l(0) = -\infty~. \end{equation}

   奇偶性:当 $l$ 为偶数,$j_l, y_l$ 分别是偶函数和奇函数,而 $l$ 为奇数时则相反。

\begin{equation} j_l(-x) = (-1)^l j_l(x)~, \qquad y_l(-y) = -(-1)^l y_l(x)~. \end{equation}

   一阶导数满足($f$ 是 $j, y, h^{(1)}, h^{(2)}$ 中的任意一种)

\begin{equation} f'_l(z) = f_{l-1}(z) - \frac{l+1}{z} f_l(z)~. \end{equation}
渐进形式 当 $x \gg 1$ 时,球贝塞尔函数的渐进表达式为
\begin{equation} j_l(x) \to \sin\left(x - l\pi /2\right) /x~, \qquad y_l(x) \to - \cos\left(x - l\pi /2\right) /x~, \end{equation}
\begin{equation} h_l^{(1)}(x) \to ( - \mathrm{i} )^{l+1} \mathrm{e} ^{ \mathrm{i} x}/x~, \qquad h_l^{(2)}(x) \to \mathrm{i} ^{l + 1} \mathrm{e} ^{- \mathrm{i} x}/x~. \end{equation}
与勒让德多项式的关系(式 9
\begin{equation} j_l(\rho) = \frac{(- \mathrm{i} )^l}{2} \int_{-1}^1 P_l(x) \mathrm{e} ^{ \mathrm{i} \rho x} \,\mathrm{d}{x} ~. \end{equation}
分解($[\dots]$ 表示 Clebsch–Gordan 系数
\begin{equation} j_l(x_1 \pm x_2) = \sum_{l_1,l_2} \mathrm{i} ^{l_1\pm l_2-l} \frac{(2l_1+1)(2l_2+1)}{2l+1} \begin{bmatrix}l_1 & l_2 & l\\ 0 & 0 & 0\end{bmatrix} ^2 j_{l_1}(x_1) j_{l_2}(x_2)~. \end{equation}

2. 积分性质

   由渐进形式可得径向归一化积分为($\delta$ 是狄拉克 $\delta$ 函数,推导见式 20

\begin{equation} \begin{aligned} \int_0^\infty k'r j_l(k'r) \cdot kr j_l(kr) \,\mathrm{d}{r} &= \int_0^\infty \sin\left(k'r - l\pi/2\right) \sin\left(kr - l\pi/2\right) \,\mathrm{d}{r} \\ & = \frac{\pi}{2}\delta(k'-k) \qquad (k, k' > 0)~, \end{aligned} \end{equation}
\begin{equation} \int_{0}^{\infty} j_l(x) \,\mathrm{d}{x} = \frac{\pi l!}{2^{l+1} [(l/2)!]^2}~. \end{equation}
使用奇偶性式 16 易得
\begin{equation} \int_{-\infty}^{+\infty} j_l(x) \,\mathrm{d}{x} = \left\{\begin{aligned} &\frac{\pi l!}{2^l [(l/2)!]^2} &\quad (l = \text{偶数})\\ &0 &\quad (l = \text{奇数}) \end{aligned}\right. ~.\end{equation}
正交性
\begin{equation} \int_{-\infty}^{+\infty} j_l(x) j_{l'}(x) \,\mathrm{d}{x} = \frac{\pi}{2l+1}\delta_{l,l'}~, \end{equation}
\begin{equation} \int_0^\infty j_l(k' r) j_l(kr) \,\mathrm{d}{r} = \frac{\pi}{2(2l+1)} \frac{k_<^l}{k_>^{l+1}}~. \end{equation}
傅里叶变换
\begin{equation} \int_{-\infty}^{+\infty} \mathrm{e} ^{ \mathrm{i} k x} j_l(x) \,\mathrm{d}{x} = \pi \mathrm{i} ^l P_l(k)~. \end{equation}

3. 修正球贝塞尔函数

   修正球贝塞尔方程

\begin{equation} x^2 \frac{\mathrm{d}^{2}{y}}{\mathrm{d}{x}^{2}} + 2x \frac{\mathrm{d}{y}}{\mathrm{d}{x}} - [x^2 + l(l + 1)]y = 0~. \end{equation}
两个线性无关解称为第一类修正球贝塞尔函数(modified spherical Bessel function of the first kind)第二类修正球贝塞尔函数
\begin{equation} i_l(x) = \sqrt{\frac{\pi }{2x}} I_{l+1/2}(x) = \mathrm{i} ^{-l} j_l( \mathrm{i} x)~, \end{equation}
\begin{equation} k_l(x) = \sqrt{\frac{2}{\pi x}} K_{l+1/2}(x) = \frac{\pi }{2} \mathrm{i} ^{l + 2} h_l^{(1)}( \mathrm{i} x)~, \end{equation}
渐进形式为
\begin{equation} i_l(x) \to \frac{ \mathrm{e} ^x}{2x}~, \qquad k_l(x) \to \frac{\pi}{2} \frac{ \mathrm{e} ^{-x}}{x}~. \end{equation}


1. ^ 参考 Wikipedia 相关页面。以及 R. Mehrem, arXiv:0909.0494 (2009)。
2. ^ 笔者也不明白为什么要这么做。


致读者: 小时百科一直以来坚持所有内容免费,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 10 元,我们一个星期内就能脱离亏损, 并保证在接下来的一整年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

                     

友情链接: 超理论坛 | ©小时科技 保留一切权利