Chandrasekhar–Page equations

Chandrasekhar–Page equations describe the wave function of the spin-1/2 massive particles, that resulted by seeking a separable solution to the Dirac equation in Kerr metric or Kerr–Newman metric. In 1976, Subrahmanyan Chandrasekhar showed that a separable solution can be obtained from the Dirac equation in Kerr metric. Later, Don Page extended this work to Kerr–Newman metric, that is applicable to charged black holes. In his paper, Page notices that N. Toop also derived his results independently, as informed to him by Chandrasekhar.

By assuming a normal mode decomposition of the form $e^{i(\sigma t+m\phi )}$ (with $m$ being a half integer and with the convention $\mathrm {Re} \{\sigma \}>0$) for the time and the azimuthal component of the spherical polar coordinates $(r,\theta ,\phi )$, Chandrasekhar showed that the four bispinor components of the wave function,

${\begin{bmatrix}F_{1}(r,\theta )\\F_{2}(r,\theta )\\G_{1}(r,\theta )\\G_{2}(r,\theta )\end{bmatrix}}e^{i(\sigma t+m\phi )}$

can be expressed as product of radial and angular functions. The separation of variables is effected for the functions $f_{1}=(r-ia\cos \theta )F_{1}$, $f_{2}=(r-ia\cos \theta )F_{2}$, $g_{1}=(r+ia\cos \theta )G_{1}$ and $g_{2}=(r+ia\cos \theta )G_{2}$ (with $a$ being the angular momentum per unit mass of the black hole) as in

$f_{1}(r,\theta )=R_{-{\frac {1}{2}}}(r)S_{-{\frac {1}{2}}}(\theta ),\quad f_{2}(r,\theta )=R_{+{\frac {1}{2}}}(r)S_{+{\frac {1}{2}}}(\theta ),$

$g_{1}(r,\theta )=R_{+{\frac {1}{2}}}(r)S_{-{\frac {1}{2}}}(\theta ),\quad g_{2}(r,\theta )=R_{-{\frac {1}{2}}}(r)S_{+{\frac {1}{2}}}(\theta ).$

Chandrasekhar–Page angular equations

The angular functions satisfy the coupled eigenvalue equations,

${\begin{aligned}{\mathcal {L}}_{\frac {1}{2}}S_{+{\frac {1}{2}}}&=-(\lambda -a\mu \cos \theta )S_{-{\frac {1}{2}}},\\{\mathcal {L}}_{\frac {1}{2}}^{\dagger }S_{-{\frac {1}{2}}}&=+(\lambda +a\mu \cos \theta )S_{+{\frac {1}{2}}},\end{aligned}}$

where $\mu$ is the particle's rest mass (measured in units so that it is the inverse of the Compton wavelength),

${\mathcal {L}}_{n}={\frac {\mathrm {d} }{\mathrm {d} \theta }}+Q+n\cot \theta ,\quad {\mathcal {L}}_{n}^{\dagger }={\frac {\mathrm {d} }{\mathrm {d} \theta }}-Q+n\cot \theta$

and $Q=a\sigma \sin \theta +m\csc \theta$. Eliminating $S_{+1/2}(\theta )$ between the foregoing two equations, one obtains

$\left({\mathcal {L}}_{\frac {1}{2}}{\mathcal {L}}_{\frac {1}{2}}^{\dagger }+{\frac {a\mu \sin \theta }{\lambda +a\mu \cos \theta }}{\mathcal {L}}_{\frac {1}{2}}^{\dagger }+\lambda ^{2}-a^{2}\mu ^{2}\cos ^{2}\theta \right)S_{-{\frac {1}{2}}}=0.$

The function $S_{+{\frac {1}{2}}}$ satisfies the adjoint equation, that can be obtained from the above equation by replacing $\theta$ with $\pi -\theta$. The boundary conditions for these second-order differential equations are that $S_{-{\frac {1}{2}}}$(and $S_{+{\frac {1}{2}}}$) be regular at $\theta =0$ and $\theta =\pi$. The eigenvalue problem presented here in general requires numerical integrations for it to be solved. Explicit solutions are available for the case where $\sigma =\mu$.

Chandrasekhar–Page radial equations

The corresponding radial equations are given by

${\begin{aligned}\Delta ^{\frac {1}{2}}{\mathcal {D}}_{0}R_{-{\frac {1}{2}}}&=(\lambda +i\mu r)\Delta ^{\frac {1}{2}}R_{+{\frac {1}{2}}},\\\Delta ^{\frac {1}{2}}{\mathcal {D}}_{0}^{\dagger }R_{+{\frac {1}{2}}}&=(\lambda -i\mu r)R_{-{\frac {1}{2}}},\end{aligned}}$

where $\Delta =r^{2}-2Mr+a^{2},$ $M$ is the black hole mass,

${\mathcal {D}}_{n}={\frac {\mathrm {d} }{\mathrm {d} r}}+{\frac {iK}{\Delta }}+2n{\frac {r-M}{\Delta }},\quad {\mathcal {D}}_{n}^{\dagger }={\frac {\mathrm {d} }{\mathrm {d} r}}-{\frac {iK}{\Delta }}+2n{\frac {r-M}{\Delta }},$

and $K=(r^{2}+a^{2})\sigma +am.$ Eliminating $\Delta ^{\frac {1}{2}}R_{+{\frac {1}{2}}}$ from the two equations, we obtain

$\left(\Delta {\mathcal {D}}_{\frac {1}{2}}^{\dagger }{\mathcal {D}}_{0}-{\frac {i\mu \Delta }{\lambda +i\mu r}}{\mathcal {D}}_{0}-\lambda ^{2}-\mu ^{2}r^{2}\right)R_{-{\frac {1}{2}}}=0.$

The function $\Delta ^{\frac {1}{2}}R_{+{\frac {1}{2}}}$ satisfies the corresponding complex-conjugate equation.

Reduction to one-dimensional scattering problem

The problem of solving the radial functions for a particular eigenvalue of $\lambda$ of the angular functions can be reduced to a problem of reflection and transmission as in one-dimensional Schrödinger equation; see also Regge–Wheeler–Zerilli equations. Particularly, we end up with the equations

$\left({\frac {d^{2}}{d{\hat {r}}_{*}^{2}}}+\sigma ^{2}\right)Z^{\pm }=V^{\pm }Z^{\pm },$

where the Chandrasekhar–Page potentials $V^{\pm }$ are defined by

$V^{\pm }=W^{2}\pm {\frac {dW}{d{\hat {r}}_{*}}},\quad W={\frac {\Delta ^{\frac {1}{2}}(\lambda +\mu ^{2}r^{2})^{3/2}}{\varpi ^{2}(\lambda ^{2}+\mu ^{2}r^{2})+\lambda \mu \Delta /2\sigma }},$

and ${\hat {r}}_{*}=r_{*}+\tan ^{-1}(\mu r/\lambda )/2\sigma$, $r_{*}=r+2M\ln(r/2M-1)$ is the tortoise coordinate and $\varpi ^{2}=r^{2}+a^{2}+am/\sigma$. The functions $Z^{\pm }({\hat {r}}_{*})$ are defined by $Z^{\pm }=\psi ^{+}\pm \psi ^{-}$, where

$\psi ^{+}=\Delta ^{\frac {1}{2}}R_{+{\frac {1}{2}}}\mathrm {exp} \left(+{\frac {i}{2}}\tan ^{-1}{\frac {\mu r}{\lambda }}\right),\quad \psi ^{-}=R_{-{\frac {1}{2}}}\mathrm {exp} \left(-{\frac {i}{2}}\tan ^{-1}{\frac {\mu r}{\lambda }}\right).$

Unlike the Regge–Wheeler–Zerilli potentials, the Chandrasekhar–Page potentials do not vanish for $r\to \infty$, but has the behaviour

$V^{\pm }=\mu ^{2}\left(1-{\frac {2M}{r}}+\cdots \right).$

As a result, the corresponding asymptotic behaviours for $Z^{\pm }$ as $r\to \infty$ becomes

$Z^{\pm }=\mathrm {exp} \left\{\pm i\left[(\sigma ^{2}-\mu ^{2})^{1/2}r+{\frac {M\mu ^{2}}{(\sigma ^{2}-\mu ^{2})^{1/2}}}\ln {\frac {r}{2M}}\right]\right\}.$