H. Норма сферической функции Бесселя для граничных условий Неймана

Для вычисления нормы решения уравнения для геометрии необходимо вычислить норму сферической функции Бесселя. Она отличается от обычной функции Бесселя множителем 1r\frac{\displaystyle 1}{\displaystyle \sqrt r}, а вес для сферических координат ρ(r)=r2\rho(r) = r^2, поэтому мы приходим к тому же интегралу, что и для обычной функции Бесселя

0μ(1rJk+1/2(r))2r2dr=0μrJk+1/22(r)dr,\int_0^{\mu} \left( \frac{\displaystyle 1}{\displaystyle \sqrt r} \cdot J_{k + 1/2}(r) \right)^2 \cdot r^2 \,dr = \int_0^{\mu} r \cdot J_{k + 1/2}^2(r) \,dr,
(H.1)

где μ\mu — одно из решений уравнения dJk+1/2(r)drr=μ=12μJk+1/2(μ)\frac{\displaystyle d J_{k + 1/2}(r)}{\displaystyle dr} \bigg|_{r=\mu} = \frac{\displaystyle 1}{\displaystyle 2 \cdot \mu} \cdot J_{k + 1/2}(\mu), которое получается из условия Неймана для сферической функции Бесселя: ddr(1rJk+1/2(r))r=μ=0\frac{\displaystyle d}{\displaystyle dr} \left( \frac{\displaystyle 1}{\displaystyle \sqrt r} \cdot J_{k + 1/2}(r) \right) \bigg|_{r=\mu} = 0.

Попробуем вычислить норму (H.1) интегрированием по частям, учитывая формулы (C.2) и (C.4):

u=rk1/2Jk+1/2(r),dv=rk+3/2Jk+1/2(r)drdu=(k12)rk3/2Jk+1/2(r)drdu+rk3/2[(k+12)Jk+1/2(r)rJk+3/2(r)]drdu=rk1/2Jk+3/2(r)drv=rk+3/2Jk+3/2(r)\begin{aligned}\left | \begin{aligned} &u = r^{-k - 1/2} \cdot J_{k + 1/2}(r), \quad dv = r^{k+3/2} \cdot J_{k + 1/2}(r) \,dr\\ &du = \left( - k - \frac{1}{2} \right) \cdot r^{-k - 3/2} \cdot J_{k + 1/2}(r) \,dr\\ &\phantom{du} + r^{-k - 3/2} \cdot \left[ \left( k + \frac{1}{2} \right) \cdot J_{k + 1/2}(r) - r \cdot J_{k + 3/2}(r) \right] \,dr\\ &du = - r^{-k - 1/2} \cdot J_{k + 3/2}(r) \,dr\\ &v = r^{k+3/2} \cdot J_{k + 3/2}(r) \end{aligned} \right |\end{aligned}
0μrJk+1/22(r)dr=rJk+1/2(r)Jk+3/2(r)0μ+0μrJk+3/22(r)dr.\int_0^{\mu} r \cdot J_{k+1/2}^2(r) \,dr = r \cdot J_{k+1/2}(r) \cdot J_{k+3/2}(r) \bigg|_0^{\mu} + \int_0^{\mu} r \cdot J_{k+3/2}^2(r) \,dr.

Преобразуем уравнение для собственных значений с помощью формулы (C.4) при r=μr = \mu

μdJk+1/2(r)drr=μ=(k+12)Jk+1/2(μ)μJk+3/2(μ)=12Jk+1/2(μ),\mu \cdot \frac{\displaystyle d J_{k + 1/2}(r)}{\displaystyle dr} \bigg|_{r=\mu} = \left( k + \frac{1}{2} \right) \cdot J_{k + 1/2}(\mu) - \mu \cdot J_{k + 3/2}(\mu) = \frac{\displaystyle 1}{\displaystyle 2} \cdot J_{k + 1/2}(\mu),
kJk+1/2(μ)=μJk+3/2(μ),k \cdot J_{k + 1/2}(\mu) = \mu \cdot J_{k + 3/2}(\mu),

и подставим полученное соотношение во внеинтегральный член

0μrJk+1/22(r)dr=kJk+1/22(μ)+0μrJk+3/22(r)dr.\int_0^{\mu} r \cdot J_{k+1/2}^2(r) \,dr = k \cdot J_{k+1/2}^2(\mu) + \int_0^{\mu} r \cdot J_{k+3/2}^2(r) \,dr.

Полученный интеграл также возьмём по частям:

u=Jk+3/22(r),du=2Jk+3/2(r)1r((k+3/2)Jk+3/2(r)+rJk+1/2(r))drdv=rdr,v=r2/2\begin{aligned}\left | \begin{array}{ll} u = J_{k+3/2}^2(r), &du = 2 \cdot J_{k+3/2}(r) \cdot \frac{\displaystyle 1}{\displaystyle r} \cdot \left( - (k+3/2) \cdot J_{k+3/2}(r) + r \cdot J_{k+1/2}(r) \right) \,dr\\ dv = r \,dr, &v = r^2 / 2 \end{array} \right |\end{aligned}
0μrJk+3/22(r)dr=k22Jk+1/22(μ)+(k+32)0μrJk+3/22(r)dr0μr2Jk+1/2(r)Jk+3/2(r)dr.\begin{split} &\int_0^{\mu} r \cdot J_{k+3/2}^2(r) \,dr = \frac{\displaystyle k^2}{\displaystyle 2} \cdot J_{k+1/2}^2(\mu) + \left( k + \frac{3}{2} \right) \cdot \int_0^{\mu} r \cdot J_{k+3/2}^2(r) \,dr -\\ &\int_0^{\mu} r^2 \cdot J_{k+1/2}(r) \cdot J_{k+3/2}(r) \,dr. \end{split}

Последний интеграл возьмём по частям:

u=rk+5/2Jk+1/2(r),dv=rk1/2Jk+3/2(r)drdu=(k+52)rk+3/2Jk+1/2(r)+rk+3/2[(k+12)Jk+1/2(r)rJk+3/2(r)]drdu=(2k+3)rk+3/2Jk+1/2(r)rk+5/2Jk+3/2(r)drv=rk1/2Jk+1/2(r)\begin{aligned}\left | \begin{aligned} &u = r^{k+5/2} \cdot J_{k+1/2}(r), \quad dv = r^{-k-1/2} \cdot J_{k+3/2}(r) \,dr\\ &du = \left( k + \frac{5}{2} \right) \cdot r^{k+3/2} \cdot J_{k+1/2}(r) + r^{k+3/2} \cdot \left[ \left( k + \frac{1}{2} \right) \cdot J_{k+1/2}(r) - r \cdot J_{k+3/2}(r) \right] \,dr\\ &du = (2 \cdot k + 3) \cdot r^{k+3/2} \cdot J_{k+1/2}(r) - r^{k+5/2} \cdot J_{k+3/2}(r) \,dr\\ &v = - r^{-k-1/2} \cdot J_{k+1/2}(r) \end{aligned} \right |\end{aligned}
0μr2Jk+1/2(r)Jk+3/2(r)dr=μ2Jk+1/22(μ)+(2k+3)0μrJk+1/22(r)dr0μr2Jk+1/2(r)Jk+3/2(r)dr,\begin{split} &\int_0^{\mu} r^2 \cdot J_{k+1/2}(r) \cdot J_{k+3/2}(r) \,dr = - \mu^2 \cdot J_{k+1/2}^2(\mu) +\\ &(2 \cdot k + 3) \cdot \int_0^{\mu} r \cdot J_{k+1/2}^2(r) \,dr - \int_0^{\mu} r^2 \cdot J_{k+1/2}(r) \cdot J_{k+3/2}(r) \,dr, \end{split}
0μr2Jk+1/2(r)Jk+3/2(r)dr=μ22Jk+1/22(μ)+(2k+3)20μrJk+1/22(r)dr.\int_0^{\mu} r^2 \cdot J_{k+1/2}(r) \cdot J_{k+3/2}(r) \,dr = - \frac{\displaystyle \mu^2}{\displaystyle 2} \cdot J_{k+1/2}^2(\mu) + \frac{\displaystyle (2 \cdot k + 3)}{\displaystyle 2} \cdot \int_0^{\mu} r \cdot J_{k+1/2}^2(r) \,dr.

Подставляя найденные равенства последовательно одно в другое, выражаем норму:

0=k22Jk+1/22(μ)+(k+12)0μrJk+3/22(r)dr0μr2Jk+1/2(r)Jk+3/2(r)dr,0 = \frac{\displaystyle k^2}{\displaystyle 2} \cdot J_{k+1/2}^2(\mu) + \left( k + \frac{1}{2} \right) \cdot \int_0^{\mu} r \cdot J_{k+3/2}^2(r) \,dr - \int_0^{\mu} r^2 \cdot J_{k+1/2}(r) \cdot J_{k+3/2}(r) \,dr,
0μrJk+3/22(r)dr=k22k+1Jk+1/22(μ)+22k+10μr2Jk+1/2(r)Jk+3/2(r)dr,\int_0^{\mu} r \cdot J_{k+3/2}^2(r) \,dr = - \frac{\displaystyle k^2}{\displaystyle 2 \cdot k + 1} \cdot J_{k+1/2}^2(\mu) + \frac{\displaystyle 2}{\displaystyle 2 \cdot k + 1} \cdot \int_0^{\mu} r^2 \cdot J_{k+1/2}(r) \cdot J_{k+3/2}(r) \,dr,
0μrJk+1/22(r)dr=k(k+1)2k+1Jk+1/22(μ)+22k+10μr2Jk+1/2(r)Jk+3/2(r)dr,\int_0^{\mu} r \cdot J_{k+1/2}^2(r) \,dr = \frac{\displaystyle k \cdot (k+1)}{\displaystyle 2 \cdot k + 1} \cdot J_{k+1/2}^2(\mu) + \frac{\displaystyle 2}{\displaystyle 2 \cdot k + 1} \cdot \int_0^{\mu} r^2 \cdot J_{k+1/2}(r) \cdot J_{k+3/2}(r) \,dr,
0μrJk+1/22(r)dr=k(k+1)μ22k+1Jk+1/22(μ)+2k+32k+10μrJk+1/22(r)dr,\int_0^{\mu} r \cdot J_{k+1/2}^2(r) \,dr = \frac{\displaystyle k \cdot (k+1) - \mu^2}{\displaystyle 2 \cdot k + 1} \cdot J_{k+1/2}^2(\mu) + \frac{\displaystyle 2 \cdot k + 3}{\displaystyle 2 \cdot k + 1} \cdot \int_0^{\mu} r \cdot J_{k+1/2}^2(r) \,dr,
k(k+1)μ22k+1Jk+1/22(μ)=22k+10μrJk+1/22(r)dr.\frac{\displaystyle k \cdot (k+1) - \mu^2}{\displaystyle 2 \cdot k + 1} \cdot J_{k+1/2}^2(\mu) = - \frac{\displaystyle 2}{\displaystyle 2 \cdot k + 1} \cdot \int_0^{\mu} r \cdot J_{k+1/2}^2(r) \,dr.

В итоге получаем:

0μrJk+1/22(r)dr=μ2k(k+1)2Jk+1/22(μ).\int_0^{\mu} r \cdot J_{k+1/2}^2(r) \,dr = \frac{\displaystyle \mu^2 - k \cdot (k+1)}{\displaystyle 2} \cdot J_{k+1/2}^2(\mu).
(H.2)