Relación de recursión de la función de Bessel

Question

Estoy leyendo un artículo y se deriva el siguiente conjunto de ecuaciones radiales:

$ -i \lbrack \partial_r + \frac{1}{r} \left( \frac{1}{2} - \nu \right) \rbrack u(r) = \pm k v(r) $

$ -i \lbrack \partial_r + \frac{1}{r} \left( \frac{1}{2} + \nu \right) \rbrack v(r) = \pm k u(r) $

Donde $\nu$ es una constante cuyo significado no importa aquí. Los autores afirman que los lados izquierdos de estas ecuaciones son una versión de las relaciones de recursión para Funciones de Bessel , por lo que las soluciones tienen la forma:

$ \left( \begin{array} &u(r) \\ v(r) \end{array} \right) = \left( \begin{array} &(i \epsilon)^{-1/2} J_{\epsilon (\nu - 1/2)}(kr) \\ \pm (i \epsilon)^{-1/2} J_{\epsilon (\nu + 1/2)}(kr) \end{array} \right)$

y $\epsilon = \pm 1$ .

No se proporciona ninguna derivación y no entiendo este paso. ¿Podría alguien describir brevemente el método utilizado aquí o proporcionar algún material de lectura útil? No he tenido suerte hasta ahora.

Answer 1

Puedes verificar la corrección de las soluciones utilizando las dos relaciones de recursión $$ J_n(z) = \frac{z}{2n}(J_{n-1}(z)+J_{n+1}(z)) $$ y $$ \frac{\mathrm d}{\mathrm dz}J_n(z) = \frac{J_{n-1}(z)-J_{n+1}(z)}{2} $$ que puede encontrar, por ejemplo, en Wikipedia .

Answer 2

De hecho, es sólo un sistema de EDOs lineales de primer orden.

$\begin{cases}-i\left[\partial_r+\dfrac{1}{r}\left(\dfrac{1}{2}-\nu\right)\right]u(r)=\pm kv(r)\\-i\left[\partial_r+\dfrac{1}{r}\left(\dfrac{1}{2}+\nu\right)\right]v(r)=\pm ku(r)\end{cases}$

$\begin{cases}\partial_ru(r)+\dfrac{1}{r}\left(\dfrac{1}{2}-\nu\right)u(r)=\pm kiv(r)\\\partial_rv(r)+\dfrac{1}{r}\left(\dfrac{1}{2}+\nu\right)v(r)=\pm kiu(r)\end{cases}$

$\therefore\partial_{rr}u(r)+\dfrac{1}{r}\left(\dfrac{1}{2}-\nu\right)\partial_ru(r)-\dfrac{1}{r^2}\left(\dfrac{1}{2}-\nu\right)u(r)=\pm ki\partial_rv(r)$

$\partial_{rr}u(r)+\dfrac{1}{r}\left(\dfrac{1}{2}-\nu\right)\partial_ru(r)-\dfrac{1}{r^2}\left(\dfrac{1}{2}-\nu\right)u(r)=\pm ki\left(-\dfrac{1}{r}\left(\dfrac{1}{2}+\nu\right)v(r)\pm kiu(r)\right)$

$\partial_{rr}u(r)+\dfrac{1}{r}\left(\dfrac{1}{2}-\nu\right)\partial_ru(r)-\dfrac{1}{r^2}\left(\dfrac{1}{2}-\nu\right)u(r)=\mp\dfrac{ki}{r}\left(\dfrac{1}{2}+\nu\right)v(r)-k^2u(r)$

$\partial_{rr}u(r)+\dfrac{1}{r}\left(\dfrac{1}{2}-\nu\right)\partial_ru(r)-\dfrac{1}{r^2}\left(\dfrac{1}{2}-\nu\right)u(r)=-\dfrac{1}{r}\left(\dfrac{1}{2}+\nu\right)\left(\partial_ru(r)+\dfrac{1}{r}\left(\dfrac{1}{2}-\nu\right)u(r)\right)-k^2u(r)$

$\partial_{rr}u(r)+\dfrac{1}{r}\left(\dfrac{1}{2}-\nu\right)\partial_ru(r)-\dfrac{1}{r^2}\left(\dfrac{1}{2}-\nu\right)u(r)=-\dfrac{1}{r}\left(\dfrac{1}{2}+\nu\right)\partial_ru(r)-\dfrac{1}{r^2}\left(\dfrac{1}{4}-\nu^2\right)u(r)-k^2u(r)$

$\partial_{rr}u(r)+\dfrac{1}{r}\partial_ru(r)+\left(k^2-\dfrac{1}{r^2}\left(\nu^2-\nu+\dfrac{1}{4}\right)\right)u(r)=0$

$r^2\partial_{rr}u(r)+r\partial_ru(r)+\biggl(k^2r^2-\left(\nu-\dfrac{1}{2}\right)^2\biggr)u(r)=0$

$u(r)=\begin{cases}C_1J_{\nu-\frac{1}{2}}(kr)+C_2Y_{\nu-\frac{1}{2}}(kr)&\text{when}~\nu-\dfrac{1}{2}\text{is an integer}\\C_1J_{\nu-\frac{1}{2}}(kr)+C_2J_{\frac{1}{2}-\nu}(kr)&\text{when}~\nu-\dfrac{1}{2}\text{is not an integer}\end{cases}$

Según http://people.math.sfu.ca/~cbm/aands/page_361.htm ,

$\partial_ru(r)=\begin{cases}-C_1\left(J_{\nu+\frac{1}{2}}(kr)-\dfrac{1}{r}\left(\nu-\dfrac{1}{2}\right)J_{\nu-\frac{1}{2}}(kr)\right)-C_2\left(Y_{\nu+\frac{1}{2}}(kr)-\dfrac{1}{r}\left(\nu-\dfrac{1}{2}\right)Y_{\nu-\frac{1}{2}}(kr)\right)&\text{when}~\nu-\dfrac{1}{2}\text{is an integer}\\-C_1\left(J_{\nu+\frac{1}{2}}(kr)-\dfrac{1}{r}\left(\nu-\dfrac{1}{2}\right)J_{\nu-\frac{1}{2}}(kr)\right)+C_2\left(J_{-\nu-\frac{1}{2}}(kr)-\dfrac{1}{r}\left(\dfrac{1}{2}-\nu\right)J_{\frac{1}{2}-\nu}(kr)\right)&\text{when}~\nu-\dfrac{1}{2}\text{is not an integer}\end{cases}$

$\therefore v(r)=\begin{cases}\mp\dfrac{i}{k}\left(-C_1\left(J_{\nu+\frac{1}{2}}(kr)-\dfrac{1}{r}\left(\nu-\dfrac{1}{2}\right)J_{\nu-\frac{1}{2}}(kr)\right)-C_2\left(Y_{\nu+\frac{1}{2}}(kr)-\dfrac{1}{r}\left(\nu-\dfrac{1}{2}\right)Y_{\nu-\frac{1}{2}}(kr)\right)+\dfrac{1}{r}\left(\dfrac{1}{2}-\nu\right)(C_1J_{\nu-\frac{1}{2}}(kr)+C_2Y_{\nu-\frac{1}{2}}(kr))\right)&\text{when}~\nu-\dfrac{1}{2}\text{is an integer}\\\mp\dfrac{i}{k}\left(-C_1\left(J_{\nu+\frac{1}{2}}(kr)-\dfrac{1}{r}\left(\nu-\dfrac{1}{2}\right)J_{\nu-\frac{1}{2}}(kr)\right)+C_2\left(J_{-\nu-\frac{1}{2}}(kr)-\dfrac{1}{r}\left(\dfrac{1}{2}-\nu\right)J_{\frac{1}{2}-\nu}(kr)\right)+\dfrac{1}{r}\left(\dfrac{1}{2}-\nu\right)(C_1J_{\nu-\frac{1}{2}}(kr)+C_2J_{\frac{1}{2}-\nu}(kr))\right)&\text{when}~\nu-\dfrac{1}{2}\text{is not an integer}\end{cases}$

$v(r)=\begin{cases}\pm\dfrac{i}{k}\left(C_1J_{\nu+\frac{1}{2}}(kr)+C_2Y_{\nu+\frac{1}{2}}(kr)\right)&\text{when}~\nu-\dfrac{1}{2}\text{is an integer}\\\pm\dfrac{i}{k}\left(C_1J_{\nu+\frac{1}{2}}(kr)-C_2J_{-\nu-\frac{1}{2}}(kr)\right)&\text{when}~\nu-\dfrac{1}{2}\text{is not an integer}\end{cases}$