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. Incidentally, while solving this problem, Chandrasekhar discovered a separable solution to the Dirac equation in flat space-time in oblate spheroidal coordinates for the first time. By assuming a normal mode decomposition of the form e i ( m ϕ − ω t ) {\displaystyle e^{i(m\phi -\omega t)}} (with m {\displaystyle m} being the azimuthal component of the particle angular momentum and takes half integer values and with ω {\displaystyle \omega } being the frequency) for the time and the azimuthal component of the spherical polar coordinates ( r , θ , ϕ ) {\displaystyle (r,\theta ,\phi )} , Chandrasekhar showed that the four Dirac spinor components of the wave function,
[ F 1 ( r , θ ) F 2 ( r , θ ) G 1 ( r , θ ) G 2 ( r , θ ) ] e i ( m ϕ − ω t ) {\displaystyle {\begin{bmatrix}F_{1}(r,\theta )\\F_{2}(r,\theta )\\G_{1}(r,\theta )\\G_{2}(r,\theta )\end{bmatrix}}e^{i(m\phi -\omega t)}}
can be expressed as product of radial and angular functions. The separation of variables is effected for the functions f 1 = ( r − i a cos θ ) F 1 {\displaystyle f_{1}=(r-ia\cos \theta )F_{1}} , f 2 = ( r − i a cos θ ) F 2 {\displaystyle f_{2}=(r-ia\cos \theta )F_{2}} , g 1 = ( r + i a cos θ ) G 1 {\displaystyle g_{1}=(r+ia\cos \theta )G_{1}} and g 2 = ( r + i a cos θ ) G 2 {\displaystyle g_{2}=(r+ia\cos \theta )G_{2}} (with a {\displaystyle a} being the angular momentum per unit mass of the black hole) as in
f 1 ( r , θ ) = R − ( r ) S − ( θ ) , f 2 ( r , θ ) = R + ( r ) S + ( θ ) , {\displaystyle f_{1}(r,\theta )=R_{-}(r)S_{-}(\theta ),\quad f_{2}(r,\theta )=R_{+}(r)S_{+}(\theta ),}
g 1 ( r , θ ) = R + ( r ) S − ( θ ) , g 2 ( r , θ ) = R − ( r ) S + ( θ ) . {\displaystyle g_{1}(r,\theta )=R_{+}(r)S_{-}(\theta ),\quad g_{2}(r,\theta )=R_{-}(r)S_{+}(\theta ).}
Chandrasekhar–Page angular equations The angular functions satisfy the coupled eigenvalue equations,
L 1 2 S + = − ( λ − a μ cos θ ) S − , L 1 2 † S − = + ( λ + a μ cos θ ) S + , {\displaystyle {\begin{aligned}{\mathcal {L}}_{\frac {1}{2}}S_{+}&=-(\lambda -a\mu \cos \theta )S_{-},\\{\mathcal {L}}_{\frac {1}{2}}^{\dagger }S_{-}&=+(\lambda +a\mu \cos \theta )S_{+},\end{aligned}}}
where μ {\displaystyle \mu } is the particle's rest mass (measured in units so that it is the inverse of the Compton wavelength), and
L n = d d θ + ( m csc θ − a ω sin θ ) + n cot θ , {\displaystyle {\mathcal {L}}_{n}={\frac {d}{{d}\theta }}+(m\csc \theta -a\omega \sin \theta )+n\cot \theta ,}
L n † = d d θ − ( m csc θ − a ω sin θ ) + n cot θ . {\displaystyle {\mathcal {L}}_{n}^{\dagger }={\frac {d}{{d}\theta }}-(m\csc \theta -a\omega \sin \theta )+n\cot \theta .}
Eliminating S + ( θ ) {\displaystyle S_{+}(\theta )} between the two equations, one obtains
( L 1 2 L 1 2 † + a μ sin θ λ + a μ cos θ L 1 2 † + λ 2 − a 2 μ 2 cos 2 θ ) S − = 0. {\displaystyle \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_{-}=0.}
The function S + {\displaystyle S_{+}} satisfies the adjoint equation, that can be obtained from the above equation by replacing θ {\displaystyle \theta } with π − θ {\displaystyle \pi -\theta } . The boundary conditions for these second-order differential equations are that S − {\displaystyle S_{-}} (and S + {\displaystyle S_{+}} ) be regular at θ = 0 {\displaystyle \theta =0} and θ = π {\displaystyle \theta =\pi } . The eigenvalue problem presented here in general requires numerical integrations for it to be solved.
Properties and limiting solutions The eigenvalue problem depends on two continuous parameters, namely a ω {\displaystyle a\omega } and a μ {\displaystyle a\mu } . For a given a ω {\displaystyle a\omega } and a μ {\displaystyle a\mu } , the eigenstates are characterised by three discrete numbers: the particle angular momentum, j = 1 / 2 , 3 / 2 , … {\displaystyle j=1/2,3/2,\dots } , its azimuthal component m = − j , − j + 1 , … , j − 1 , j {\displaystyle m=-j,-j+1,\dots ,j-1,j} and the parity P = ± 1. {\displaystyle {\mathcal {P}}=\pm 1.} The spectrum elements may then be explicitly labelled as
S ± =
s = ± 1 2 S j , m , P ( a ω , a μ ) , λ = λ j , m , P ( a ω , a μ ) . {\displaystyle S_{\pm }={}_{s=\pm {\frac {1}{2}}}S_{j,m,{\mathcal {P}}}^{(a\omega ,a\mu )},\lambda =\lambda _{j,m,{\mathcal {P}}}^{(a\omega ,a\mu )}.}
The eigenvalue λ {\displaystyle \lambda } has the physical interpretation of being the square root of the generalised total anagular momentum squared. The knowledge of the spectrum in the positive quadrant a ω > 0 , a μ > 0 {\displaystyle a\omega >0,a\mu >0} is sufficient to determine the full spectrum, as implied by the symmetry:
λ j , m , P ( a ω , a μ ) = − λ j , − m , − P ( − a ω , a μ ) = − λ j , m , − P ( a ω , − a μ ) = λ j , − m , P ( − a ω , − a μ ) . {\displaystyle \lambda _{j,m,{\mathcal {P}}}^{(a\omega ,a\mu )}=-\lambda _{j,-m,-{\mathcal {P}}}^{(-a\omega ,a\mu )}=-\lambda _{j,m,-{\mathcal {P}}}^{(a\omega ,-a\mu )}=\lambda _{j,-m,{\mathcal {P}}}^{(-a\omega ,-a\mu )}.}
Furthermore
s S j , m , P ( a ω , a μ ) ( θ ) = ( − 1 ) s − 1 / 2
s S j , m , − P ( a ω , − a μ ) ( Θ ) = P ( − 1 ) m − 1 / 2
− s S j , − m , − P ( − a ω , a μ ) ( θ ) = P ( − 1 ) j + m
− s S j , m , P ( a ω , a μ ) ( π − θ ) {\displaystyle {}_{s}S_{j,m,{\mathcal {P}}}^{(a\omega ,a\mu )}(\theta )=(-1)^{s-1/2}{}_{s}S_{j,m,-{\mathcal {P}}}^{(a\omega ,-a\mu )}(\Theta )={\mathcal {P}}(-1)^{m-1/2}{}_{-s}S_{j,-m,-{\mathcal {P}}}^{(-a\omega ,a\mu )}(\theta )={\mathcal {P}}(-1)^{j+m}{}_{-s}S_{j,m,{\mathcal {P}}}^{(a\omega ,a\mu )}(\pi -\theta )}
and the combinations thereof. Non-rotating black hole (Schwarzschild black hole) a = 0 {\displaystyle a=0} : The problem can be solved explicitly. The eigenvalues and eigenfunctions are given by
λ j , m , P ( 0 , 0 ) = P ( j + 1 / 2 ) , {\displaystyle \lambda _{j,m,{\mathcal {P}}}^{(0,0)}={\mathcal {P}}(j+1/2),}
[
+ 1 / 2 S j , m , P ( 0 , 0 ) ( θ )
− 1 / 2 S j , m , P ( 0 , 0 ) ( θ ) ] = A [ cos θ 2 sin θ 2 − sin θ 2 cos θ 2 ] [ P j + P / 2 m + 1 / 2 ( cos θ ) c j , m , P ( 0 , 0 ) P j + P / 2 m − 1 / 2 ( cos θ ) ] , {\displaystyle {\begin{bmatrix}{}_{+1/2}S_{j,m,{\mathcal {P}}}^{(0,0)}(\theta )\\{}_{-1/2}S_{j,m,{\mathcal {P}}}^{(0,0)}(\theta )\end{bmatrix}}=A{\begin{bmatrix}\cos {\frac {\theta }{2}}&\sin {\frac {\theta }{2}}\\-\sin {\frac {\theta }{2}}&\cos {\frac {\theta }{2}}\end{bmatrix}}{\begin{bmatrix}P_{j+{\mathcal {P}}/2}^{m+1/2}(\cos \theta )\\c_{j,m,{\mathcal {P}}}^{(0,0)}P_{j+{\mathcal {P}}/2}^{m-1/2}(\cos \theta )\end{bmatrix}},}
where P {\displaystyle P} is the associated Legendre polynomials and
c j , m , P ( 0 , 0 ) = P ( j + P / 2 + 1 / 2 ) − m , A = ( j − m ) ! 2 π ( j + m ) ! . {\displaystyle c_{j,m,{\mathcal {P}}}^{(0,0)}={\mathcal {P}}(j+{\mathcal {P}}/2+1/2)-m,\quad A={\sqrt {\frac {(j-m)!}{2\pi (j+m)!}}}.}
Special case ω = ± μ {\displaystyle \omega =\pm \mu } : For the special case where ω = + μ {\displaystyle \omega =+\mu } , the solutions are given by
λ j , m , P ( a ω , a ω ) = − 1 2 + P ( j + P / 2 + 1 / 2 ) 2 − 2 m a ω + a 2 ω 2 , {\displaystyle \lambda _{j,m,{\mathcal {P}}}^{(a\omega ,a\omega )}=-{\frac {1}{2}}+{\mathcal {P}}{\sqrt {(j+{\mathcal {P}}/2+1/2)^{2}-2ma\omega +a^{2}\omega ^{2}}},}
[
+ 1 / 2 S j , m , P ( 0 , 0 ) ( θ )
− 1 / 2 S j , m , P ( 0 , 0 ) ( θ ) ] = A [ cos θ 2 sin θ 2 − sin θ 2 cos θ 2 ] [ P j + P / 2 m + 1 / 2 ( cos θ ) c j , m , P ( a ω , a ω ) P j + P / 2 m − 1 / 2 ( cos θ ) ] , {\displaystyle {\begin{bmatrix}{}_{+1/2}S_{j,m,{\mathcal {P}}}^{(0,0)}(\theta )\\{}_{-1/2}S_{j,m,{\mathcal {P}}}^{(0,0)}(\theta )\end{bmatrix}}=A{\begin{bmatrix}\cos {\frac {\theta }{2}}&\sin {\frac {\theta }{2}}\\-\sin {\frac {\theta }{2}}&\cos {\frac {\theta }{2}}\end{bmatrix}}{\begin{bmatrix}P_{j+{\mathcal {P}}/2}^{m+1/2}(\cos \theta )\\c_{j,m,{\mathcal {P}}}^{(a\omega ,a\omega )}P_{j+{\mathcal {P}}/2}^{m-1/2}(\cos \theta )\end{bmatrix}},}
where
c j , m , P ( a ω , a ω ) = [ ( j + P / 2 + 1 / 2 ) 2 − m 2 ] / ( λ j , m , P ( a ω , a ω ) + m + 1 / 2 − a ω ) , {\displaystyle c_{j,m,{\mathcal {P}}}^{(a\omega ,a\omega )}=[(j+{\mathcal {P}}/2+1/2)^{2}-m^{2}]/(\lambda _{j,m,{\mathcal {P}}}^{(a\omega ,a\omega )}+m+1/2-a\omega ),}
A = ( j + P / 2 + 1 / 2 ) 2 π ( j + P / 2 − m − 1 / 2 ) ! ( j + P / 2 + m + 1 / 2 ) ! ( 1 + P ( m − a ω ) ( j + P / 2 + 1 / 2 ) 2 − 2 m a ω + a 2 ω 2 ) . {\displaystyle A={\sqrt {{\frac {(j+{\mathcal {P}}/2+1/2)}{2\pi }}{\frac {(j+{\mathcal {P}}/2-m-1/2)!}{(j+{\mathcal {P}}/2+m+1/2)!}}\left(1+{\frac {{\mathcal {P}}(m-a\omega )}{\sqrt {(j+{\mathcal {P}}/2+1/2)^{2}-2ma\omega +a^{2}\omega ^{2}}}}\right)}}\,.}
When ω = − μ {\displaystyle \omega =-\mu } , we have
λ j , m , P (