Access the full text.
Sign up today, get DeepDyve free for 14 days.
Loading next page...
/lp/arxiv-cornell-university/an-integral-spectral-representation-of-the-massive-dirac-propagator-in-p4hGn65vD1
References
References for this paper are not available at this time. We will be adding them shortly, thank you for your patience.
- ISSN
- 1095-0761
- eISSN
- ARCH-3337
- DOI
- 10.4310/ATMP.2018.v22.n1.a3
- Publisher site
- See Article on Publisher Site
Abstract
An Integral Spectral Representation of the Massive Dirac Propagator in the Kerr Geometry in Eddington{Finkelstein-type Coordinates 1 1, 2 Felix Finster and Christian R oken Universit at Regensburg, Fakult at fur Mathematik, 93040 Regensburg, Germany Departamento de Geometr a y Topolog a, Facultad de Ciencias - Universidad de Granada, Campus de Fuentenueva s/n, 18071 Granada, Spain (Dated: August 2018) ABSTRACT. We consider the massive Dirac equation in the non-extreme Kerr geometry in horizon-penetrating advanced Eddington{Finkelstein-type coordinates and derive a functional ana- lytic integral representation of the associated propagator using the spectral theorem for unbounded self-adjoint operators, Stone's formula, and quantities arising in the analysis of Chandrasekhar's separation of variables. This integral representation describes the dynamics of Dirac particles out- side and across the event horizon, up to the Cauchy horizon. In the derivation, we rst write the Dirac equation in Hamiltonian form and show the essential self-adjointness of the Hamiltonian. For the latter purpose, as the Dirac Hamiltonian fails to be elliptic at the event and the Cauchy horizon, we cannot use standard elliptic methods of proof. Instead, we employ a new, general method for mixed initial-boundary value problems that combines results from the theory of symmetric hyper- bolic systems with near-boundary elliptic methods. In this regard and since the time evolution may not be unitary because of Dirac particles impinging on the ring singularity, we also impose a suit- able Dirichlet-type boundary condition on a time-like inner hypersurface placed inside the Cauchy horizon, which has no e ect on the dynamics outside the Cauchy horizon. We then compute the resolvent of the Dirac Hamiltonian via the projector onto a nite-dimensional, invariant spectral eigenspace of the angular operator and the radial Green's matrix stemming from Chandrasekhar's separation of variables. Applying Stone's formula to the spectral measure of the Hamiltonian in the spectral decomposition of the Dirac propagator, that is, by expressing the spectral measure in terms of this resolvent, we obtain an explicit integral representation of the propagator. Contents I. Introduction 2 II. Preliminaries 3 III. Hamiltonian Formulation of the Massive Dirac Equation in the Non-extreme Kerr Geometry in Horizon-penetrating Coordinates 7 IV. The Canonical Scalar Product 8 V. Essential Self-adjointness of the Dirac Hamiltonian 11 VI. Resolvent of the Dirac Hamiltonian and Integral Spectral Representation of the Dirac Propagator 15 Acknowledgments 23 e-mail: felix. nster@mathematik.uni-regensburg.de e-mail: christian.roeken@mathematik.uni-regensburg.de arXiv:1606.01509v3 [gr-qc] 31 Aug 2018 2 A. Symmetry of the Dirac Hamiltonian and Dirichlet-type Boundary Condition 23 B. Fundamental Solutions for the Construction of the Radial Green's Matrix 26 References 30 I. INTRODUCTION In [11], an integral spectral representation of the propagator of the massive Dirac equation in the non- extreme Kerr geometry outside the event horizon is derived in Boyer{Lindquist coordinates. It has been used to study the long-time behavior (including decay rates) and the escape probability of Dirac particles in rotating Kerr black hole spacetimes [10]. The shortcoming of this integral spectral representation is, however, that it yields a solution of the associated Cauchy problem only outside the event horizon. In the present paper, we construct a generalized integral spectral representation that describes the complete dynamics of Dirac particles outside, across, and inside the event horizon, up to the Cauchy horizon. The methods used in the derivation of our integral spectral representation are quite di erent from those employed in [11], as is now outlined. We work with horizon-penetrating advanced Eddington{ Finkelstein-type coordinates [22], i.e., with an analytic extension of Boyer{Lindquist coordinates that simultaneously covers both the exterior and the interior black hole region without exhibiting singularities at the horizons and features a proper time function . Furthermore, we employ a regular Carter tetrad, i.e., a symmetric Newman{Penrose null frame, which re ects, on the one hand, the discrete time and angle reversal isometries of the Kerr geometry and, on the other hand, its Petrov type. After computing the corresponding spin coecients by solving the rst Maurer{Cartan equation of structure, we explicitly determine the massive Dirac equation in Hamiltonian form i@ ( ; r; ; ) = H ( ; r; ; ) ; where H denotes the Hamiltonian, is a Dirac 4-spinor, and ( ; r; ; ) are the advanced Eddington{ Finkelstein-type coordinates. Moreover, we introduce a suitable scalar product on the associated space of solutions, and show for smooth and compactly supported Dirac 4-spinors that the Hamiltonian is symmetric with respect to this scalar product. We also establish that it coincides with the canonical scalar product obtained by integrating the normal component (de ned with respect to the level sets of ) of the Dirac current. As we apply the spectral theorem in the derivation of the propagator, we need to establish the essential self-adjointness of the Dirac Hamiltonian. To this end, we rst impose a Dirichlet-type boundary condition on a time-like inner boundary surface placed beyond the Cauchy horizon. This boundary condition prevents Dirac particles from impinging on the curvature singularity without a ecting their dynamics in the region outside the Cauchy horizon, consequently leading to a unitary time evolution. Then, we apply the method of proof for the essential self-adjointness of the Dirac Hamiltonian for mixed initial-boundary value problems that are not uniformly elliptic introduced in [15]. Subsequently, it is possible to derive an integral representation of the Dirac propagator via the spectral theorem for unbounded self-adjoint operators i H i! = e = e dP ; 0 0 ! where dP is the spectral measure of H , ! is the spectral parameter, and := ( = 0; r; ; ) is ! 0 smooth initial data with compact support. We compute the spectral measure of the Dirac Hamiltonian employing Stone's formula [20], which yields an explicit expression in terms of the resolvent of the Hamiltonian Res(H ) := (H ! ) , for which ! 2 CnR. Furthermore, we determine the resolvent c c applying quantities obtained in the analysis of the systems of radial and angular ordinary di erential equations (ODEs) that arise in Chandrasekhar's separation of variables, that is, after factoring out the azimuthal angle modes, we project the Dirac Hamiltonian onto a nite-dimensional, invariant spectral eigenspace of the angular operator, which leaves us with a matrix-valued rst-order ordinary 3 di erential operator in the radial variable. The resolvent of this operator can be calculated by means of the Green's matrix of the radial ODE system. For this purpose, we derive generalized Jost-type equations [21] for this system and study speci c aspects of their solutions, namely the existence, uniqueness, and boundedness. We moreover use the asymptotic radial solutions at in nity, the event horizon, and the Cauchy horizon for guidance in the implicit construction of the fundamental solutions of the radial ODE system required for the computation of the corresponding Green's matrix. Eventually, by summing over all azimuthal angle modes, we obtain the full resolvent of the Dirac Hamiltonian in separated form. The resulting horizon-penetrating generalization of the integral spectral representation of the Dirac propagator describes the complete dynamics of massive Dirac par- ticles outside and across the event horizon of the non-extreme Kerr geometry, up to the Cauchy horizon. Main Theorem. The massive Dirac propagator in the non-extreme Kerr geometry in horizon- penetrating advanced Eddington{Finkelstein-type coordinates can be expressed via the integral spectral representation ik i! 1 1 0 0 0 0 ( ; r; ; ) = e e lim (H ! i) (H ! + i) (r; ; r ; ) (r ; ) d! ; k k 0;k &0 2i k2Z where is the initial data for xed k-modes and (H ! i) are the unique resolvents of the 0;k k Dirac Hamiltonian H for xed k-modes on the upper and lower complex half-planes. In a nal step, we compute the limit & 0 of the di erence of resolvents leading to a rigorously simpli ed form of the propagator. The article is organized as follows. In Section II, we provide the mathematical framework for the Kerr geometry and for the massive Dirac equation. Moreover, we recall required results from the asymptotic analysis of the radial ODE system and from the spectral analysis of the angular ODE system arising in Chandrasekhar's separation of variables without giving proofs. We derive the Hamiltonian formulation and a suitable scalar product for the space of solutions of the associated Cauchy problem in Sections III and IV, respectively. Furthermore, we verify the symmetry of the Hamiltonian with respect to this scalar product in Appendix A. In Section V, we show the essential self-adjointness of the Hamiltonian. Finally, we construct the resolvent of the Hamiltonian and the integral spectral representation of the propagator in Section VI. The fundamental radial solutions required for the computation of the resolvent are determined in Appendix B. II. PRELIMINARIES We recall the necessary basics on the non-extreme Kerr geometry in horizon-penetrating advanced Eddington{Finkelstein-type coordinates, the general relativistic, massive Dirac equation in the Newman{Penrose formalism, and Chandrasekhar's separation of variables. The non-extreme Kerr geometry is a connected, orientable and time-orientable, smooth, asymptot- 2 2 ically at Lorentzian 4-manifold (M;g) with topology S R , for which the metric g is stationary and axisymmetric and given in horizon-penetrating advanced Eddington{Finkelstein-type coordinates ( ; r; ; ) with 2 R; r 2 R ; 2 [0; ], and 2 [0; 2) [22] by >0 2Mr 2Mr 2 2 g = 1 d d dr a sin () d d + d dr a sin () d (1) 2Mr 2 2 2 1 + dr a sin () d dr a sin () d d d sin () d d ; where M is the mass and aM the angular momentum of the black hole, with 0 a < M , and 2 2 2 2 2 = (r; ) := r + a cos (). The event and Cauchy horizons are located at r := M M a , respectively. The advanced Eddington{Finkelstein-type coordinates are an analytic extension of the common Boyer{Lindquist coordinates (t; r; ; ') with t 2 R; r 2 R ; 2 [0; ], and ' 2 [0; 2) [3], >0 covering both the exterior and interior black hole regions while being regular at the horizons. In terms of the Boyer{Lindquist coordinates, the advanced Eddington{Finkelstein-type time and azimuthal angle coordinates read 2 2 2 2 r + a r + a := t + lnjr r j lnjr r j r r r r + + (2) a r r := ' + ln : r r r r This horizon-penetrating coordinate system possesses a proper time function, unlike the original ad- vanced Eddington{Finkelstein (null) coordinates [8, 9]. It is advantageous to describe the Kerr geometry in the Newman{Penrose formalism using a regular Carter tetrad [4, 22] l = p [ + 4Mr] @ + @ + 2a @ 2 r n = p (@ @ ) (3) m = p ia sin () @ + @ + i csc () @ m = ia sin () @ @ + i csc () @ ; 2 2 with = (r) := (r r )(r r ) = r 2Mr + a being the horizon function, because this frame is adapted to the two principal null directions of the Weyl tensor and to the fundamental discrete time and angle reversal isometries. Thus, since the Kerr geometry is algebraically special and of Petrov type D, one has the computational advantage that the four spin coecients ; ; , and as well as the four Weyl scalars ; ; , and vanish [18], and that speci c spin coecients are linearly dependent. 0 1 3 4 Substituting the Carter tetrad (3) into { and solving { the rst Maurer{Cartan equation of structure, we obtain the spin coecients [22] 2 2 r r a 2ia cos () (r M ) p p = = = = 0 ; = ; = 3=2 3=2 2 r ia cos () 2 r r ia cos () ia sin () r p p p (4) = = ; = ; % = 2 r ia cos () 2 r ia cos () 2 r r ia cos () 2 2 = = r + a cot () ira sin () : 3=2 (2) 4 4 Next, introducing a spin bundle SM = MC on M with bers S M ' C , x 2 M, we can formulate the general relativistic, massive Dirac equation (without an external potential) r + im (x ) = 0 ; 2 f0; 1; 2; 3g ; (5) where r is the metric connection on SM, are the Dirac matrices, is the Dirac 4-spinor de ned on the bers S M, and m is the invariant fermion mass. In the Newman{Penrose formalism { by x 5 employing a local dyad spinor frame { (5) becomes the coupled rst-order system of partial di erential equations im (n @ + )G (m @ + )G = p F 1 2 1 im (l @ + " %)G (m @ + )G = F 2 1 2 (6) im (l @ + " %)F + (m @ + )F = p G 1 2 1 im (n @ + )F + (m @ + )F = p G 2 1 2 with = (F ;F ; G ; G ) [5]. Inserting the Carter tetrad (3) and the associated spin coecients 1 2 1 2 (4) into the system (6), and applying the transformation 0 T 0 1 = P = (H ;H ; J ; J ) ; = P P ; (7) 1 2 1 2 where p p p p P := diag r ia cos (); r ia cos (); r + ia cos (); r + ia cos () ; (8) we nd r @ @ J + ia sin () @ @ + i csc () @ 2 cot () J = im r + ia cos () H + r 1 2 1 1 1 r [ + 4Mr] @ + @ + 2a @ + r M J ia sin () @ + @ + i csc () @ + 2 cot () J r 2 1 = im r + ia cos () H 1 1 r [ + 4Mr] @ + @ + 2a @ + r M H ia sin () @ @ + i csc () @ 2 cot () H r 1 2 = im r ia cos () J r @ @ H + ia sin () @ + @ + i csc () @ + 2 cot () H = im r ia cos () J ; + r 2 1 2 (9) which is the starting point for the derivation of the Hamiltonian formulation of the massive Dirac equation on a Kerr background geometry in horizon-penetrating coordinates presented in the next section. We note in passing that the system (9) corresponds to the transformed Dirac equation 0 y 1 0 1 0 P P r + P @ P + im = 0 ; (10) where := diag(1; 1; 1; 1). This will become relevant in the following construction of both the Hamiltonian formulation and the scalar product. Finally, for the explicit computation of the resolvent of the Dirac Hamiltonian, we require speci c results arising from Chandrasekhar's separation of variables of the system (9). More precisely, we need the asymptotic solutions of the radial ODE system at in nity, at the event horizon, and at the Cauchy horizon, as well as certain information about the eigenvalues and eigenfunctions of the angular ODE system. In the following, these results are recalled. For a detailed analysis and proofs see [22]. Substituting the separation ansatz i(! +k ) H = e R (r)T () 1 + + i(! +k ) H = e R (r)T () (11) i(! +k ) J = e R (r)T () 1 + i(! +k ) J = e R (r)T () ; 2 + 6 in which ! 2 R and k 2 Z + 1=2, into (9) yields the rst-order radial and angular ODE systems jj 0 1 e e R(r )R = R (12) 2 2 r + a sign() 0 A()T = T ; where 2 2 2 2 r + a r + a r := r + lnjr r j lnjr r j ? + r r r r + + is the Regge{Wheeler coordinate, 0 p 1 !( + 4Mr) 2ak jj mr @ A R(r ) := 1 1 2 @ + ? C r 2 2 r + a jj sign() mr ! 0 1 ma cos () @ 2 cot () + a! sin () + k csc() @ A A() := (13) @ + 2 cot () + a! sin () + k csc() ma cos () e T are matrix-valued radial and angular operators, R := jjR ; r R and T := (T ;T ) are + + + radial and angular vector-valued functions, and is the constant of separation. The asymptotic solutions and the decay properties of the associated errors of the radial ODE system at in nity, the event horizon, and the Cauchy horizon are speci ed in the lemmas below. Lemma II.1. Every nontrivial solution R of (12) for j!j m is asymptotically as r ! 1 of the oscillatory form (1) i (r ) + ? f e e e R(r ) = R (r ) + E (r ) = D + E (r ) ; ? 1 ? 1 ? 1 1 ? (2) i (r ) f e where 1 ! m cosh ( ) sinh ( D = with := ln ; sinh ( ) cosh ( 4 ! + m the functions p 2 2 2 (r ) = sign(!) ! m r + M 2! p ln (r ) (14) ? ? ? 2 2 ! m (1) (2) are the asymptotic phases, and f = f ; f 6= 0 is a vector-valued constant. The error E has 1 1 1 polynomial decay e e kE (r )k = R(r ) R (r ) 1 ? ? 1 ? for a suitable constant a 2 R . For the case j!j < m, the non-trivial solution R has both contributions >0 p p 2 2 2 2 m ! r m ! r ? ? that show exponential decay e and exponential growth e . Lemma II.2. Every nontrivial solution R of (12) is asymptotically as r & r of the form () (1) 2i !+k r? Kerr g e e e R(r ) = R (r ) + E (r ) = + E (r ) ? r ? r ? r ? (2) (1) (2) () with the constants g = g ; g 6= 0 and := a=(2Mr ), as well as an error with exponential r r Kerr decay q r e e ? kE (r )k = R(r ) R (r ) p e r ? ? r ? for suitable constants p ; q 2 R . >0 The spectral properties of the eigenvalues and eigenfunctions of the angular ODE system are summarized in the following proposition. Proposition II.3. For any ! 2 R and k 2 Z + 1=2, the di erential operator (13) has a complete set of orthonormal eigenfunctions (T ) in L (0; ); sin () d . The corresponding eigenvalues are l l2Z l real-valued and non-degenerate, and can thus be ordered as < . Moreover, the eigenfunctions are l l+1 pointwise bounded and smooth away from the poles, 2 2 1 1 T 2 L (0; ) \ C (0; ) : Both the eigenfunctions T and the eigenvalues depend smoothly on !. l l III. HAMILTONIAN FORMULATION OF THE MASSIVE DIRAC EQUATION IN THE NON-EXTREME KERR GEOMETRY IN HORIZON-PENETRATING COORDINATES For the derivation of the Hamiltonian formulation of the massive Dirac equation in the non-extreme Kerr geometry in horizon-penetrating advanced Eddington{Finkelstein-type coordinates, it is advantageous to rst rewrite the system (9) in the form (R +A) = 0 ; where 0 1 imr 0 D 0 B C 0 imr 0 D B C R := (15) @ A D 0 imr 0 0 D 0 imr and 0 1 ma cos () 0 0 L B C 0 ma cos () L 0 B C A := (16) @ A 0 L ma cos () 0 L 0 0 ma cos () are matrix-valued di erential operators with D := r [ + 4Mr] @ + @ + 2a @ + r M + r D := r (@ @ ) + r L := ia sin () @ + @ + i csc() @ + 2 cot () : We then separate the -derivative and multiply the resulting equation by the inverse of the matrix 0 1 0 0 r ia sin () p B C 0 0 ia sin () r [ + 4Mr] 0 0 y 1 0 B C e := P P = (17) @ A r [ + 4Mr] ia sin () 0 0 ia sin () r 0 0 + 8 (cf. Eq. (10)) as well as by the imaginary unit, which leads to the Schr odinger-type equation 0 0 1 (3) (3) 0 0 i@ = i ( e ) R +A =: H ; (18) (3) (3) where R and A contain the rst-order spatial and all zero-order contributions of the operators (15) and (16), respectively. The Dirac Hamiltonian H may be recast in the more convenient form H = @ + V ; j 2 fr; ; g ; (19) with the matrix-valued coecients 0 1 i r a sin () 0 0 B C 1 r a sin () i ( + 4Mr) 0 0 B C := (20) @ A + 2Mr 0 0 i ( + 4Mr) r a sin () 0 0 r a sin () i 0 1 a sin () ir 0 0 B C 1 ir [ + 4Mr] a sin () 0 0 B C := (21) @ A + 2Mr 0 0 a sin () ir [ + 4Mr] 0 0 ir a sin () 0 1 ia r csc() 0 0 B C 1 r csc() ( 2) ia 0 0 B C := (22) @ A + 2Mr 0 0 ia r csc() ( 2) 0 0 r csc() ia and the potential B B 1 2 V := ; (23) B B + 2Mr 3 4 where the quantities B , k 2 f1; 2; 3; 4g, are the (2 2)-blocks 0 1 2i (r M ) a cos () ir cot () @ A B := r 2a sin () (r M ) + i cot () ( + 4Mr) a cos () 0 1 r ia sin () @ A B := m r ia cos () ia sin () r ( + 4Mr) (24) 0 1 r ( + 4Mr) ia sin () @ A B := m r + ia cos () ia sin () r 0 1 a cos () r 2a sin () (r M ) i cot () ( + 4Mr) @ A B := : ir cot () 2i (r M ) + a cos () IV. THE CANONICAL SCALAR PRODUCT In order to set up a Hilbert space that contains the solutions of (18) 0 2 0 H := Sol(H i@ ); (j ) for which Sol(H i@ ) = 2 L (M; SM)j (H i@ ) = 0 ; and to establish the symmetry property of the Hamiltonian H 0 0 0 0 0 0 ( jH ) = (H j ) with ; 2 Sol(H i@ ) (or its self-adjointness), we require a suitable scalar product (j ). We thus work with the scalar product [11] ( j ) := j= d (25) jN jN de ned on the space-like hypersurface N := f = const.; r; ; g, where j : S M S M ! C ; ( ; ) 7! (26) x x ? y denotes the inde nite spin scalar product of signature (2; 2), := S the adjoint Dirac spinor, = = the Cli ord contraction of the future-directed, time-like normal , and d = jN jdet(g )j d d dr is the invariant measure on N , in which g is the induced Riemannian metric. jN jN The matrix S is de ned via the relation y 1 := S S : (27) We note that this scalar product is independent of the choice of the speci c space-like hypersurface N . This can be easily shown by means of Gauss' theorem and current conservation. In the following, we explicitly compute the above quantities and subsequently derive a more convenient representation for the scalar product. We begin with the calculation of the matrix S . Via (17) and the spinor transformation (7), we nd y 1 0 0 = p (P ) e P and hence, using the de ning equation (27), we obtain the expression 0 1 0 0 1 0 B C 0 0 0 1 B C S = : @ A 1 0 0 0 0 1 0 0 Next, we determine the normal vector eld by means of the conditions hj@ i = hj@ i = hj@ i = 0 and hji = 1 ; jg jg jg jg where hji := g( ; ) = g dx dx ( ; ) is the spacetime inner product on M. Accordingly, employing (1) yields 1=2 1=2 2Mr 2Mr 2Mr = 1 + @ 1 + @ : The corresponding dual co-vector reads 1=2 2Mr = 1 + d : (28) Moreover, the induced metric g on the hypersurface N is simply the restriction of (1) to N and, jN thus, we obtain 2Mr 2 2 2 g = 1 + dr a sin () d dr a sin () d d d sin () d d : jN The associated Jacobian determinant in the volume measure d becomes jN q 1=2 2Mr jdet(g )j = sin () 1 + : (29) jN We now express the scalar product (25) in terms of the primed quantities (7) used in (18), and substitute (28) as well as (29). This results in ZZZ 0 0 0y 0 0 0 ( j ) = S sin () d d dr : (30) jN 0 y 1 0 0 Again employing (17), that is with = P (P ) e = , the scalar product (30) yields ZZZ 0 0 0y 0 y 1 0 0 0 ( j ) = S P (P ) e sin () d d dr jN ZZZ 0y y 0 y 1 0 0 0 = PP S P (P ) e sin () d d dr ZZZ 0y y 1 0 0 0 = PS (P ) e sin () d d dr (31) ZZZ 0y y y 1 0 0 0 = SP (P ) e sin () d d dr ZZZ 0y 0 = sin () d d dr ; where 0 1 r [ + 4Mr] ia sin () 0 0 B C ia sin () r 0 0 0 0 + B C := S e = : (32) @ A 0 0 r ia sin () 0 0 ia sin () r [ + 4Mr] We point out that in the above derivation, we have rst applied the relation 1 1 4 = PP , then the 0 y 0 transformation law for the matrix S , namely S = P S P , and nally we have used the fact that both S and the product PS are self-adjoint, which leads to PS = SP . Besides, the integration limits are suppressed for ease of notation. The eigenvalues ; of the matrix (32) are positive 1 2 1 + 4Mr + 4Mr = r + + r + 4a sin () > 0 1 + + 2 r r + + 1 + 4Mr + 4Mr = r + r + 4( + 2Mr) > 0 2 + + 2 r r + + and with algebraic multiplicities ( ) = ( ) = 2, demonstrating that (31) is indeed positive- A 1 A 2 de nite. The symmetry property of the Dirac Hamiltonian (19) with respect to this scalar product is explicitly proven in Appendix A. 11 V. ESSENTIAL SELF-ADJOINTNESS OF THE DIRAC HAMILTONIAN In this section, we show that the Dirac Hamiltonian in the non-extreme Kerr geometry in horizon- penetrating advanced Eddington{Finkelstein-type coordinates is essentially self-adjoint using the results obtained in [15] (we recently learned that in [19] related results were found with di erent methods). Having an essentially self-adjoint Hamiltonian is necessary for the derivation of the integral representa- tion of the Dirac propagator presented in the subsequent section, as it is based on the spectral theorem for unbounded, self-adjoint operators. The proof of the essential self-adjointness involves the technical diculty that in the Kerr geometry the Dirac Hamiltonian is only almost everywhere elliptic, and hence not uniformly elliptic. More precisely, it fails to be elliptic at the event and the Cauchy horizon. This can be easily seen from the evaluation of the determinant of the principal symbol of the Hamiltonian (19) P (r; ;) = (r; ) ; (33) ? j where 2 T N . In more detail, we rst rewrite the matrices in terms of the original Dirac matrices j 0 1 0j 0 1 0j 1 j 1 = i ( e ) e = i ( ) = iP ( ) P : Substituting this expression into the principal symbol (33) and computing the determinant yields det( ) det P (r; ;) = : det( ) Using the relations 2 i j ij ( ) = g 1 1 and = g 1 1 ; S M i j i j S M x x we obtain ij i j det P (r; ;) = : (34) The Hamiltonian fails to be elliptic if its principal symbol is not invertible, that is, if the determinant (34) vanishes. This is the case for ij g = 0 with 6= 0 ; (35) i j ij where the quantities g are the components of the inverse of the spacetime metric (1) g = [ + 2Mr] @ @ 2Mr (@ @ + @ @ ) @ r r r r a (@ @ + @ @ ) @ @ csc () @ @ : r r Analyzing condition (35), we nd that the Hamiltonian is not elliptic at the event and the Cauchy horizon. We point out that by using the intrinsic Dirac Hamiltonian on the space-like hypersurface N , ellipticity would be conserved even at the horizons, which can be inferred from the analog condition ij g = 0 with 6= 0 ; i j jN ij where the g denote the components of the inverse of the associated induced Riemannian metric jN 2 2 2 2 r + a a sin () g = @ @ + a [@ @ + @ @ ] + @ @ + csc () @ @ : r r r r jN + 2Mr 12 FIG. 1: Carter{Penrose diagram of the region M of the non-extreme Kerr geometry with constant- hyper- surfaces N and N cut-o at the boundary @M . A radial Dirichlet-type boundary condition imposed on @M 1 2 leads to a re ection of Dirac particles away from the singularity without a ecting their dynamics outside the Cauchy horizon. This is represented by Cauchy data propagated in -direction. However, as we work with the Hamiltonian obtained from the Dirac operator in the full Kerr spacetime, ellipticity is broken. Therefore, we cannot employ standard techniques and results from elliptic theory in order to verify the essential self-adjointness of the Dirac Hamiltonian. Instead, we apply the results derived in [15], where near-boundary elliptic methods are combined with results from the theory of symmetric hyperbolic systems (see, e.g., [2, 6, 16, 23]). In the following, we state the geometrical and functional analytic settings for the formulation of the Cauchy problem for the massive Dirac equation in the non-extreme Kerr geometry in horizon-penetrating coordinates in Hamiltonian form, which is used as a technical tool in the proof of the essential self-adjointness of the Dirac Hamiltonian. We let (M;g) be the non-extreme Kerr geometry with the metric (1) in horizon-penetrating advanced Eddington{Finkelstein-type coordinates (2) and consider the subset M := f ; r > r ; ; g M ; where r < r : 0 0 Furthermore, we introduce the time-like inner boundary @M := f ; r = r ; ; g of M and the family of space-like hypersurfaces N = (N ) ; where N := f = const.; r > r ; ; g ; 2R 0 with boundaries @N := @M \ N ' S (see FIG. 1). These hypersurfaces constitute a foliation of M . Moreover, near @M , we de ne a locally time-like Killing vector eld [18] K := @ + @ with = (r ) 2 Rnf0g ; (36) 0 0 0 0 13 where @ and @ are the Killing elds describing the stationarity and axisymmetry of the Kerr geometry, respectively. We note that the speci c proof of existence of unique, global solutions of the Cauchy problem for the massive Dirac equation in Hamiltonian form for a general class of mixed initial-boundary value problems on Lorentzian manifolds presented in [15] makes essential use of a Killing eld K = @ that is tangential to and time-like on the inner boundary @M and represented by a coordinate system describing an observer who is co-moving along the associated ow lines. This Killing eld may also be space-like or null in Mn@M . A direct computation shows that the Killing eld K = @ in advanced Eddington{Finkelstein-type coordinates is not everywhere time-like on @M . To be more precise, the condition for @ not being time-like reads 2Mr g(@ ; @ ) = 1 0 : This inequality is solved by p p 2 2 2 2 2 2 M M a cos () r M + M a cos () ; which is the ergosphere region. But taking K as the linear combination (36), it turns out to be a Killing eld that satis es all the above assumptions. To make the connection between (36) and the Killing eld K = @ , we use the coordinate transformation R R [0; ] [0; 2) ! R R [0; ] [0; 2) ; ( ; r; ; ) 7! (t; r; ; ) >0 >0 with t = and = : (37) This transformation can be easily derived from the condition @ @ K = @ = @ + @ = @ + @ : t 0 @t @t Evaluation of the gradient rt = g (@ t) @ = g @ and subsequently 2Mr t tt g(@ ;rt) = g = 1 > 0 as well as g(rt;rt) = g = 1 + > 0 demonstrates that rt is future-pointing and time-like and, hence, that the coordinate t is a time function as is the original time coordinate [22]. Due to the speci c form of the transformation (37), we nd that the induced metric g on the level sets of t is identical to the induced metric jt=const. 2Mr 2Mr g = g = 1 + dr dr + a sin () 1 + dr d + d dr j =const. jN 2Mr 2 2 2 d sin () + a sin () 1 + d on the level sets of for the advanced Eddington{Finkelstein-type coordinates. As a consequence, all the results obtained for the co-moving coordinate system in [15] will also hold true for the advanced Eddington{Finkelstein-type coordinates. Next, in addition to these geometric structures, we introduce the spin bundle SM of M with bers S M ' C , where x 2 M . We may then consider the Dirac Hamiltonian given in (19) H = @ + V on N ; (38) j 14 which is a symmetric operator with respect to the scalar product speci ed in (31) Z Z 0 0 0y 0 ( j ) = sin () d d dr ; (39) jN r S 0 0 where the Dirac 4-spinors ; 2 S M and the matrix is de ned in (32). We note that for H being symmetric with respect to (39), we have to impose the radial Dirichlet-type boundary condition (see Appendix A) y 1 0 n= iP (P ) = 0 ; j@M in which n is the inner normal to @M and P is determined by (8). This boundary condition has the e ect that Dirac particles are re ected at @M away from the singularity such that, without changing their dynamics outside the Cauchy horizon, we obtain a unitary time evolution. The speci c domain of the Hamiltonian reads 0 1 y 1 0 Dom(H ) = 2 C (N ; SM ) n= iP (P ) = 0 : j@N In this setting, we nd a unique, global solution of the Cauchy problem for the massive Dirac equation in Hamiltonian form in the class C (M; SM ). sc Lemma V.1. The Cauchy problem for the massive Dirac equation in the non-extreme Kerr geometry in horizon-penetrating advanced Eddington{Finkelstein-type coordinates 0 0 i@ = H 0 0 1 =: 2 C (N ; SM ) =0 j =0 0 0 with the radial Dirichlet-type boundary condition at @M given by y 1 0 n iP (P ) = 0 ; (40) j@M where the initial data is smooth, compactly supported outside, across, and inside the event horizon, up to the Cauchy horizon, and is compatible with the boundary condition, i.e., y 1 p 0 n iP (P ) (H ) = 0 8 p 2 N ; j@N 0 has a unique, global solution in the class of smooth functions with spatially compact support 1 0 C (M; SM ). Evaluating this solution at subsequent times and gives rise to a unique unitary sc propagator ; 1 1 U : C (N ; SM ) ! C (N ; SM ) : 0 0 The proof of this lemma is shown in detail for a more general class of non-uniformly elliptic mixed initial-boundary value problems for the Dirac equation in Hamiltonian form on Lorentzian manifolds 0 1 with dimension d 3 in [15]. The existence of a unique, global solution 2 C (M; SM ) is imperative sc for the speci c proof of the essential self-adjointness of the Dirac Hamiltonian presented in the same work. Below, we state the result for the Hamiltonian in the non-extreme Kerr geometry in horizon- penetrating coordinates (38). Theorem V.2. The massive Dirac Hamiltonian H in the non-extreme Kerr geometry in horizon- penetrating advanced Eddington{Finkelstein-type coordinates with domain of de nition 0 1 y 1 p 0 Dom(H ) = 2 C (N ; SM ) n= iP (P ) (H ) = 0 8 p 2 N j@N 0 is essentially self-adjoint. 15 VI. RESOLVENT OF THE DIRAC HAMILTONIAN AND INTEGRAL SPECTRAL REPRESENTATION OF THE DIRAC PROPAGATOR We may now construct an integral spectral representation of the Dirac propagator that yields the dynamics of massive Dirac particles outside, across, and inside the event horizon, up to the Cauchy horizon. More precisely, we derive an explicit expression for the spectral measure dP of the essentially self-adjoint Dirac Hamiltonian H de ned in (38) with the domain speci ed in Theorem V.2, which arises in the formal spectral decomposition of the Dirac propagator 0 i H 0 i! 0 ( ; r; ; ) = e (r; ; ) = e (r; ; ) dP ; (41) 0 0 where is smooth, spatially compact initial data. To this end, we employ Stone's formula and, thus, express the spectral measure in terms of the resolvent Res(H ) = (H ! ) of the Hamiltonian. As the spectrum of the Hamiltonian (H ) R is on the real line, this resolvent exists for all ! 2 CnR with real part Re(! ) = ! 2 (H ) and is given uniquely. In the computation of the resolvent, we make use of quantities obtained in the analysis of Chandrasekhar's separation of variables, namely the angular projector onto a nite-dimensional, invariant eigenspace of the angular operator (13) and the Green's matrix of the radial ODE system (12). Theorem VI.1. The massive Dirac propagator in the non-extreme Kerr geometry in horizon- penetrating advanced Eddington{Finkelstein-type coordinates can be expressed via the integral spectral representation 0 ik i! 1 1 0 0 0 0 0 ( ; r; ; ) = e e lim (H ! i) (H ! + i) (r; ; r ; ) (r ; ) d! ; k k 0;k &0 2i k2Z 0 1 1 where 2 C (r ;1) [0; ]; SM is the initial data for xed k-modes and (H ! i) are the 0 k 0;k resolvents of the Dirac Hamiltonian for xed k-modes H on the upper and lower complex half-planes. The resolvents are unique and of the form Z Z 1 1 G(r; r ) 0 2 1 0 0 0 0 0 0 k;l;!i C (H ! i) (r; ; r ; ) (r ; ) = Q (; ) C k l 0;k 0 0 2 G(r; r ) C k;l;!i 1 r l2Z 0 0 0 0 0 0 0 E (r ; ) (r ; ) dr d cos ( ) 0;k with Q ( ; ) being the integral kernel of the spectral projector onto a nite-dimensional, invariant eigenspace of the angular operator (13) that corresponds to the spectral parameter , G(r; r ) the l k;l;!i two-dimensional Green's matrix of the radial rst-order ODE system (12), and 0 1 0 1 0 0 0 1 0 0 0 i (r ) + 4Mr r a sin ( ) 0 0 B C B C 0 0 0 1 0 0 ir a sin ( ) 0 0 B C B C C = ; E (r ; ) = : 0 0 0 @ A @ A 0 1 0 0 0 0 r a sin ( ) i (r ) + 4Mr 0 0 1 0 a sin ( ) ir 0 0 Proof. We rst compute the resolvents (H ! i) of the Dirac Hamiltonian for xed k-modes H , where ! 2 R and > 0 is suciently small so that it can be considered as a slightly non-self- adjoint perturbation. For this purpose, we begin by substituting Chandrasekhar's mode ansatz with complex-valued frequencies 0 i(! +k ) 2 ( ; r; ; ) = e (r; ) in which 2 L (r ;1) [0; ]; SM into the Dirac equation (18) restricted to M , yielding (H ! ) = 0 ; (42) k 16 r j where H := @ + @ ik + V with the Dirac matrices , j 2 fr; ; g, and the potential k r V given in (20)-(23). We remark that we introduced the abbreviation ! 2 C with Re(! ) = ! and Im(! ) 2 f ; g in order to cover the resolvents in both the upper and lower complex half-planes simultaneously. Next, we de ne the spectral projector 0 0 0 Q := Q (; ) (r; ) d cos ( ) l l onto the nite-dimensional, invariant eigenspace of the matrix-valued angular operator (13), which results from Chandrasekhar's separation of variables, corresponding to the spectral parameter with l 2 Z. This spectral projector is idempotent Q = Q for all n 2 N ; and the family of spectral projectors (Q ) is complete l l2Z Q = 1 1 : (43) l2Z We may therefore express the angular operator (13) by means of the family (Q ) as l l2Z A = Q : l l l2Z Applying this representation of the angular operator and the completeness constraint (43) in Eq. (42), we obtain ( + 2Mr) M(@ ; r; ) Q = 0 ; (44) r k;l;! l l2Z where 0 1 i O a sin () U ir S a sin () S k;! ! + l l B C B C B C a sin () i( + 4Mr) i( + 4Mr) B C O U a sin () S S k;! ! l l B C r r r + + + B C B C M(@ ; r; ) := r k;l;! B C B C i( + 4Mr) i( + 4Mr) a sin () B C S a sin () S U O l l ! k;! B C r r r + + + B C @ A a sin () S ir S a sin () U i O l + l ! k;! with the purely radial di erential operators O , U , and the function S de ned by k;! ! l O := @ + r M i! ( + 4Mr) 2iak k;! r U := r (@ + i! ) ! + r S := + imr : l l In the following, we show that the computation of the resolvent of the operator M(@ ; r; ) in (44) r k;l;! can be reduced to determining the two-dimensional Green's matrix of the radial ODE system (12). Writing the Dirac equation (44) in the factorized form ( + 2Mr) B(r; ) R(@ ; r) Q = 0 ; r k;l;! l l2Z 17 where the matrix B(r; ) and the matrix-valued radial operator R(@ ; r) read r k;l;! 0 1 i a sin () 0 0 B C a sin () i( + 4Mr) B C 0 0 B C r r + + B C B(r; ) := B C i( + 4Mr) a sin () B C 0 0 @ A r r + + 0 0 a sin () i and 0 1 O 0 r S 0 k;! + l B C 0 U 0 S ! l B C R(@ ; r) := ; r k;l;! @ A S 0 U 0 l ! 0 r S 0 O + l k;! we can easily bring it into the block diagonal representation X 22 R (@ ; r) 0 2 1 r k;l;! C 1 E (r; ) C Q = 0 ; (45) 22 0 R (@ ; r) C r k;l;! l2Z in which 22 R (@ ; r) 0 2 r k;l;! C 1 = C R(@ ; r) C (46) r k;l;! 22 0 2 R (@ ; r) C r k;l! with O r S k;! + l 22 R (@ ; r) := ; (47) r k;l;! S U l ! and the matrices E (r; ) and C are de ned by 1 1 E (r; ) := ( + 2Mr) B(r; )C and 0 1 1 0 0 0 B C 0 0 0 1 B C C := : @ A 0 1 0 0 0 0 1 0 From the speci c form of (45), it is obvious that the key quantity in the determination of the resolvent of 1 0 M(@ ; r; ) , and thus of the resolvent (H ! ) , is the solution G(r; r ) of the distributional r k;l;! k k;l;! equation 22 0 0 R (@ ; r) G(r; r ) = (r r ) 1 1 2 : (48) r k;l;! k;l;! C We point out that (47) is identical to the operator (12), but for a complex-valued frequency ! with imaginary part Im(!) 2 f ; g. Hence, G(r; r ) corresponds to the Green's matrix of the radial k;l;! ODE system obtained via Chandrasekhar's separation of variables. However, in the case of a complex- valued frequency, the solution of the radial ODE system has an additional dampening contribution guaranteeing that (r; ) is in L (r ;1) [0; ]; SM , which is in contrast to the original case with ! 2 R. In order to solve the distributional equation (48), we rst introduce the vector-valued functions 0 0 (r; r ) (r; r ) 0 1;1 0 2;1 (r; r ) = and (r; r ) = 1 2 0 0 (r; r ) (r; r ) 1;2 2;2 that 18 for r 6= r are linearly independent solutions of the homogeneous equation 22 0 R (@ ; r) (r; r ) = 0 ; r k;l;! have jump discontinuities at r = r , satisfy the Dirichlet-type boundary condition (40) at r = r , and are square-integrable, that is 2 2 0 0 (r; r ) = (r; r ) dr < 1 : 1=2 1=2 These functions are speci ed in Appendix B. It turns out that their r -dependence can be chosen in such a way that it is solely contained in Heaviside step functions . For clarity, these Heaviside step functions are explicitly stated in what follows, which makes it possible to consider and as functions 1 2 of only the variable r. From the rst and the last of the above properties as well as from Lemma II.1 and Lemma II.2 (but with a complex-valued frequency ! 2 f! + i; ! ig and with the substitution p p 2 2 2 2 ! m ! j! j m ), we can moreover infer that they have the speci c asymptotics ! ( 1;1 Im(! ) < 0 if j! j m i (r (r)) + ? (r) e for r ! 1 and 1=2 Re(! ) 0 if j! j < m 2;1 ! ( 3;1 Im(! ) > 0 if j! j m i (r (r)) (r) e for r ! 1 and 1=2 Re(! ) < 0 if j! j < m 4;1 0 1 () 1;r 2i ! +k r (r) Kerr e r ! r and Im(! ) < 0 @ A (r) for jj 1=2 r ! r and Im(! ) > 0 2;r r ! r and Im(! ) > 0 (r) c for 1=2 3;r r ! r and Im(! ) < 0 ; where c and c , with n 2 f1; 2; 3; 4g as well as m 2 f1; 2; 3g, are scalar constants. We then use n;1 m;r the ansatz 0 0 0 0 0 0 (r r ) (r) P (r ) + (r r) (r) P (r ) for r < r < 1 and r r r 1 1 2 2 + 0 G(r; r ) = k;l;! 0 0 0 0 0 (r r ) (r) P (r ) + (r r ) (r) P (r ) for r < r r 1 1 2 2 + (49) for the radial Green's matrix in case j! j m and Im(! ) < 0 or j! j < m and Re(! ) 0 ; whereas in case j! j m and Im(! ) > 0 or j! j < m and Re(! ) < 0 ; we employ the ansatz 0 0 0 0 0 0 (r r ) (r) P (r ) + (r r) (r) P (r ) for r < r < 1 and r r r 1 1 2 2 + 0 G(r; r ) = k;l;! 0 0 0 0 0 (r r) (r) P (r ) + (r r) (r) P (r ) for r < r r ; 1 1 2 2 + (50) 19 in which P and P are unknowns yet to be determined. Applying the radial operator (47) to (49) and 1 2 (50), we obtain 0 0 0 0 (r )P (r ) (r )P (r ) for (49) 0 1 1 2 2 22 0 0 R (@ ; r) G(r; r ) = (r r ) r k;l;! k;l;! 0 0 0 0 0 r + (r )P (r ) (r )P (r ) for (50) : 1 1 2 2 Comparing these equations with (48) yields the two systems 1 0 0 0 0 (r )P (r ) (r )P (r ) for (49) 1 1 2 2 0 0 0 0 0 r (r )P (r ) (r )P (r ) for (50) : 1 1 2 2 Their solutions P (r ) read 1=2 0 0 0 0 (r ) (r ) (r ) (r ) 2;2 2;1 1;2 1;1 0 0 0 0 P (r ) = ; P (r ) = ; P (r ) = ; P (r ) = 1;1 1;2 2;1 2;2 0 0 0 0 0 0 (r ) W (r ) r W (r ) (r ) W (r ) r W (r ) + + for the ansatz (49) and 0 0 0 0 (r ) (r ) (r ) (r ) 2;2 2;1 1;2 1;1 0 0 0 0 P (r ) = ; P (r ) = ; P (r ) = ; P (r ) = 1;1 1;2 2;1 2;2 0 0 0 0 0 0 (r ) W (r ) r W (r ) (r ) W (r ) r W (r ) + + for ansatz (50), where 0 0 0 0 0 0 W (r ) = W ( ; )(r ) := (r ) (r ) (r ) (r ) 1 2 1;1 2;2 1;2 2;1 is the Wronskian. Substituting these expressions into (49) and (50), respectively, leads in case r < 0 0 r < 1 or r r r to the Green's matrix 2 0 1 0 0 (r) (r ) (r) (r ) 1;1 2;2 1;1 2;1 6 B 0 C (r ) r 6 B C 0 0 6 B C G(r; r ) = (r r ) k;l;! 0 6 B C W (r ) 0 0 4 @ A (r) (r ) (r) (r ) 1;2 2;2 1;2 2;1 (r ) r k;l;! (51) 0 1 3 0 0 (r) (r ) (r) (r ) 2;1 1;2 2;1 1;1 B 0 C 7 (r ) r B C 7 B C 7 +(r r) B C 7 0 0 @ A 5 (r) (r ) (r) (r ) 2;2 1;2 2;2 1;1 (r ) r k;l;! for both (49) and (50), whereas in case r < r r it leads to the Green's matrices 0 1 0 0 0 0 (r) (r ) (r) (r ) (r) (r ) (r) (r ) 1;1 2;2 2;1 1;2 2;1 1;1 1;1 2;1 B 0 C (r ) r B C B C G(r; r ) = k;l;! 0 B C W (r ) 0 0 0 0 @ A (r) (r ) (r) (r ) (r) (r ) (r) (r ) 1;2 2;2 2;2 1;2 2;2 1;1 1;2 2;1 (r ) r k;l;! (r r ) for (49) (r r) for (50) : (52) Subsequently, we may directly read o the resolvent of the Dirac Hamiltonian for xed k-modes from the block-diagonalized representation (45). We thus nd X 0 G(r; r ) 0 2 1 k;l;! C 0 0 0 (H ! ) = Q C E (r ; ) (r ; ) dr ; (53) k l 0 G(r; r ) C k;l;! l2Z 20 where the Green's matrix G(r; r ) is given in (51) and (52). To show that this expression is actually k;l;! the desired resolvent, we verify the identity (H ! )(H ! ) = : (54) k k Accordingly, applying the operator in (45) to (53), we obtain in a rst step X 22 X R (@ ; r) 0 2 1 1 r k;l;! C 1 0 (H ! )(H ! ) = E (r; ) C Q (; ) k k l 22 0 R (@ ; r) C r k;l;! l2Z m2Z Z Z 1 1 G(r; r ) 0 k;m;! C 0 00 0 00 0 00 0 00 0 Q ( ; ) C E (r ; ) (r ; ) dr d cos ( ) d cos ( ) : 0 2 G(r; r ) k;m;! 1 r 0 00 Moving the integral kernel of the spectral projector Q (; ) into the -integral and taking into account that the spectral projectors are idempotent, i.e., their integral kernels satisfy the relation 0 0 00 0 00 Q (; ) Q ( ; ) = cos () cos ( ) Q (; ) ; l m lm m we infer, after evaluating the -integral and the sum over all integers m, that X 22 R (@ ; r) 0 2 r k;l;! C 1 1 1 (H ! )(H ! ) = E (r; ) C k k 22 0 2 R (@ ; r) r k;l;! l2Z Z Z 1 1 G(r; r ) 0 2 k;l;! 00 C 0 00 0 00 0 00 Q (; ) C E (r ; ) (r ; ) dr d cos ( ) : 0 2 G(r; r ) C k;l;! 1 r Next, we can also move the constant matrix C as well as the matrix-valued radial operator (46) into 00 0 the - and the r -integral. Employing (48) yields 1 1 0 0 0 0 (H ! )(H ! ) = E (r; ) Q (r r )E (r ; ) (r ; ) dr : k k l l2Z Solving the integral with respect to the variable r and substituting the completeness constraint for the spectral projectors (43), we immediately obtain the identity (54). Having established the explicit form of the resolvent (H ! ) in (53), we continue deriving the integral spectral representation of the Dirac propagator. To this end, we express the Dirac spinor at ;0 i H 0 time in terms of the propagator U = e applied to smooth, spatially compact initial data at time = 0 and expand the initial data in terms of k-modes 0 i H 0 ik i H 0 = e = e e : (55) 0 0;k k2Z We furthermore introduce the spectral projector of the Dirac Hamiltonian for xed k-modes H onto the interval I R P (H ) := (H ) ; I k I k where denotes the characteristic function 1 for ! 2 I (H ) := I k 0 for ! 2= I with ! 2 (H ). Then, by making use of the identity relation P (H ) = 1 1 ; ( 1;1) 21 we write (55) as X X 0 ik i H 0 ik i H 0 k k = e e lim P (H ) = e lim e P (H )+P (H ) : ( a;a) k ( a;a) k [ a;a] k 0;k 0;k a!1 a!1 k2Z k2Z Employing Stone's formula for the spectral projector of an unbounded, self-adjoint operator [20], which in our framework reads i H 0 i! 1 1 0 e P (H ) + P (H ) = lim e (H ! i) (H ! + i) d! ; ( a;a) k [ a;a] k k k 0;k 0;k &0 i we obtain 0 ik i! 1 1 0 = e lim lim e (H ! i) (H ! + i) d! ; k k 0;k 2i a!1 &0 k2Z where the resolvents are given by (53). Finally, since the fundamental solutions that occur in the resolvents are bounded for all > 0 and all ! 2 R as shown in Appendix B, we can apply Lebesgue's dominated convergence theorem and commute the -limit and the integral with respect to !, yielding 0 ik i! 1 1 0 0 0 0 0 ( ; r; ; ) = e e lim (H ! i) (H ! + i) (r; ; r ; ) (r ; ) d! : k k 0;k 2i &0 k2Z (56) We note that by comparing this expression to formula (41), we may directly identify the spectral measure dP of the Dirac Hamiltonian H . This integral spectral representation can be further simpli ed, on the one hand, by performing the limit r % r and, on the other hand, by computing the di erence of the two resolvents (H ! i) and 0 k (H ! + i) for & 0. In the following, we explicitly work out the case j! j m. The case j! j < m 0 0 may be treated similarly. As the fundamental solutions (r; r ) and (r; r ), which constitute the 1 2 0 1 radial Green's matrix G(r; r ) and therefore the resolvent (H ! ) , are given piecewise for the k;l;! k 0 0 0 domains r < r < 1, r < r r , and r r r (see the second part of Appendix B), we begin + + 0 by splitting the r -integral in the di erence of resolvents in the limit & 0 into the three associated contributions Z Z Z r r 1 1 1 0 lim (H ! i) (H ! + i) = lim Q + + C k k l 0;k &0 &0 r r r 0 + l2Z (57) 0 0 G(r; r ) G(r; r ) 0 2 k;l;! i k;l;!+i C 0 0 0 0 E (r ; ) (r ; ) dr : 0;k 0 0 0 2 G(r; r ) G(r; r ) C k;l;! i k;l;!+i Because the integrands in (57), and hence the lth summand, are bounded for all values of , r , and (see the rst part of Appendix B and keeping in mind that the initial data for xed k-modes has 0;k spatially compact support), we can again employ Lebesgue's dominated convergence theorem, which allows us to commute the limit & 0 with the sum over the integers l, the spectral projector Q , and the integrals with respect to r . Then, applying the limit r % r to the integral spectral representation (56) and commuting this limit with the sum over the integers k, the integral with respect to !, and last the sum over the integers l as well as the spectral projector Q (in the di erence of resolvents) using the same reasoning as before, we obtain the expression Z Z r 1 X + 1 1 0 lim lim (H ! i) (H ! + i) = Q + C k k l 0;k r %r &0 r r l2Z 0 1 0 0 lim G(r; r ) G(r; r ) 0 2 k;l;! i k;l;!+i &0 0 0 0 0 @ A E (r ; ) (r ; ) dr : 0 0 0;k 0 2 lim G(r; r ) G(r; r ) C k;l;! i k;l;!+i &0 (58) 22 In order to compute the limit & 0 of the di erence of the radial Green's matrices G(r; r ) and k;l;! i 0 0 G(r; r ) in the domain r < r < 1, we introduce the auxiliary functions (see the second part of k;l;!+i + Appendix B) (1) (1) q b (r) := lim (r) and (r) := lim (r) ; 1 2 &0 &0 and write the -limits of the fundamental radial solutions and as 1 2 lim = and lim = + for Im(! ) > 0 1 1 2 1 2 &0 &0 (59) lim = and lim = + for Im(! ) < 0 ; 1 2 2 1 2 &0 &0 where , , , and are constants. The corresponding Wronskian yields W ( ; ) for Im(! ) > 0 1 2 lim W ( ; ) = (60) 1 2 &0 W ( ; ) for Im(! ) < 0 : 1 2 Substitution of (59) and (60) into (51) results in 0 0 lim G(r; r ) G(r; r ) k;l;! i k;l;!+i jr <r <1 &0 0 1 0 0 (r) (r ) (r) (r ) u;1 v;2 u;1 v;1 B 0 C (r ) r + (61) B C B C = T u;v 0 B C W ( ; )(r ) 1 2 0 0 @ A (r) (r ) (r) (r ) u;v=1 u;2 v;2 u;2 v;1 (r ) r k;l;! with the coecients T = ; T = T = 1 ; and T = : 1;1 1;2 2;1 2;2 For the domain r < r r on the other hand, we de ne the auxiliary functions ( ) (1) ( ) b b q e (r) := (r r) lim (r) + (r r ) lim (r) and e (r) := (r r ) lim (r) ; 1 + + 2 &0 &0 &0 and express the -limits of the fundamental radial solutions by 0 0 lim = e + e and lim = e for Im(! ) > 0 1 1 2 2 2 &0 &0 (62) 0 0 lim = e and lim = e + e for Im(! ) < 0 ; 1 1 2 1 2 &0 &0 0 0 0 0 in which , , , and are also constants. In this case, the Wronskian becomes W ( e ; e ) for Im(! ) > 0 1 2 lim W ( ; ) = (63) 1 2 &0 W ( e ; e ) for Im(! ) < 0 : 1 2 Using (62) and (63) in (52) to calculate the di erence of Green's matrices for & 0 leads to 0 0 lim G(r; r ) G(r; r ) k;l;! i k;l;!+i 0 jr <r r &0 0 1 0 0 e (r) e (r ) e (r) e (r ) u;1 v;2 u;1 v;1 B C (r ) r X + (64) B C e B C = T u;v B C W ( e ; e )(r ) 0 0 1 2 @ A u;v=1 e (r) e (r ) e (r) e (r ) u;2 v;2 u;2 v;1 (r ) r k;l;! 23 with the coecients e e e e T = T = 0 and T = T = 1 : 1;1 2;2 1;2 2;1 0 0 Abbreviating (61) and (64) by G(r; r < r < 1) and G(r; r < r r ) , respectively, and + k;l;! + k;l;! inserting these quantities into (58), the Dirac propagator (56) yields Z Z + 0 G(r; r < r r ) 0 2 0 ik i! + k;l;! C ( ; r; ; ) = e e Q C 0 2 G(r; r < r r ) 2i C + k;l;! R r k;l2Z G(r; r < r < 1) 0 + k;l;! C 0 0 0 0 + E (r ; ) (r ; ) dr d! 0 0;k 0 2 G(r; r < r < 1) + k;l;! Z Z ik i! 0 = e e Q C 1 1 2 G(r; r < r r ) l C + k;l;! 2i R r k;l2Z 0 0 0 0 0 + G(r; r < r < 1) E (r ; ) (r ; ) dr d! : + k;l;! 0;k Given in this form, the generalized, horizon-penetrating integral spectral representation of the massive Dirac propagator in the non-extreme Kerr geometry resembles the one restricted to the region outside the event horizon derived in [11]. Acknowledgments The authors are grateful to Niky Kamran, Guillaume Idelon{Riton, and Simone Murro for useful discussions and comments. This work was supported by the DFG research grant \Dirac Waves in the Kerr Geometry: Integral Representations, Mass Oscillation Property and the Hawking E ect." Appendix A: Symmetry of the Dirac Hamiltonian and Dirichlet-type Boundary Condition In this appendix, we show the symmetry of the Dirac Hamiltonian H with respect to the canonical scalar product (j ) on the space-like hypersurface N by direct computation. Furthermore, we introduce jN and discuss the relevant radial Dirichlet-type boundary condition imposed on the Dirac spinors. Theorem A.1. The Dirac Hamiltonian (19) is symmetric with respect to the scalar product (31). Proof. To establish the symmetry, namely that 0 0 0 0 ( jH ) = (H j ) ; jN jN we begin by splitting the potential V given in (23) into mass-independent and mass-dependent parts V = V + V ; 0 m where 1 B 0 2 1 0 2 B 1 C C 2 V := and V := : 0 m 2 2 + 2Mr 0 B + 2Mr B 0 C 4 3 C 24 The (2 2)-blocks B , with k 2 f1; 2; 3; 4g, are speci ed in (24). This procedure bears the advantage of obtaining anti-self-adjoint and self-adjoint matrices y y V = V and V = V ; (A1) 0 m 0 m for which = is de ned in (32). We may then write ZZZ 0 0 0y 0 ( jH ) = H sin () d d dr jN ZZZ ZZZ 0y j 0 0y 0 = @ ( ) sin () d d dr + V sin () d d dr j 0 ZZZ 0y 0 + V sin () d d dr : Integration by parts of the rst triple integral in the second line and substitution of the relations (A1) in the remaining two triple integrals results in ZZZ ZZZ 0 0 0y j 0 0y 0 ( jH ) = @ sin () d d dr V sin () d d dr jN 0 ZZZ 0y y 0 + V sin () d d dr ZZZ 0y j 0 = @ ( ) sin () d d dr (A2) ZZZ 0y j j 0 @ ( ) + @ ( ) + cot () sin () d d dr j j ZZZ ZZZ 0 y 0 0 y 0 (V ) sin () d d dr + (V ) sin () d d dr : 0 m We remark that in the integration by parts, the angular derivatives do not give rise to boundary terms 2 2 because the two-dimensional submanifold S in N ' R S is compact without boundary. The >0 radial derivative, on the other hand, yields a boundary term that vanishes as we impose an appropriate Dirichlet-type boundary condition on the Dirac spinors. More precisely, since the computation of the matrix leads to the expression = i diag ; r ; r ; ; + + r r + + the radial boundary term becomes ZZ ZZ r r 2 2 0y r 0 0 0 0 0 0 0 0 0 sin () d d = ir + + sin () d d : 1 1 2 2 3 3 4 4 2 2 r r 2 r 2 r 1 1 S S + + In order for this term to vanish, we impose the radial Dirichlet-type boundary condition i 0 0 0 0 0 0 0 0 ( 1) + + = 0 : (A3) 1 1 2 2 3 3 4 4 2 2 r r + + jr=r i=1 25 In the present work, we consider only Dirac spinors with support from a speci c time-like inner boundary at r = r < r beyond the Cauchy horizon up to in nity, that is 0 2 supp = (r ;1) S : 2 2 Moreover, we require the Dirac spinors to be in L (r ;1)S ; SM , implying proper decay at in nity. Taking this into account, the radial boundary condition (A3) reduces to a condition for the time-like inner boundary at r = r 0 0 0 0 0 0 0 0 + + = 0 ; (A4) 1 1 2 2 3 3 4 4 2 2 r r + + jr=r which can be brought into a more suitable form as follows. By means of the spin scalar product (26) 0 0r r and the relation S = i =, we may represent (A4) as 0 0r 0 j 2 = 0 : (A5) jf gfr gS Now, introducing n as the unit normal to the hypersurfaces f g S , we can write (A5) in the form 0 0 0 jn= 2 = 0 ( (n= iH) = 0 ; (A6) jf gfr gS 0 jf gfr gS where the slash again denotes Cli ord multiplication and H is an arbitrary matrix with the property 0 1 y 0 H = (S ) H S . The implication can be easily veri ed via the calculation 0 0 0 0 0 0 jn 2 = n j 2 = iH j 2 = = jf gfr gS jf gfr gS jf gfr gS 0 0 0 0 0 0 0 2 2 = j iH = jn= : jf gfr gS jf gfr gS 0 0 To guarantee compatibility of the boundary condition on the right hand side of (A6) with a potential product structure of the Dirac 4-spinors, in which the dependences on the variables ; r; , and are separated (such as in Chandrasekhar's separation ansatz (11)), we choose y 1 H = P (P ) ; where P is de ned in (8). We note in passing that this Dirichlet-type boundary condition is a so-called MIT-type boundary condition for Dirac elds [7] that describes a perfect re ection of Dirac particles at the respective boundary surface. Continuing the proof of symmetry, the explicit computation of the square bracket in the fourth line of (A2) yields the result j j @ ( ) + @ ( ) + cot () = 2V : j j Besides, all three matrix products , with j 2 fr; ; g, are anti-self-adjoint j jy y jy = = : Therefore, we immediately nd that ZZZ ZZZ 0 0 0y jy 0 0 y 0 ( jH ) = @ ( ) sin () d d dr + 2 (V ) sin () d d dr jN j 0 ZZZ ZZZ 0 y 0 0 y 0 (V ) sin () d d dr + (V ) sin () d d dr 0 m ZZZ ZZZ j 0 y 0 0 y 0 = ( @ ) sin () d d dr + (V ) sin () d d dr j 0 ZZZ 0 y 0 + (V ) sin () d d dr ZZZ 0 y 0 0 0 = (H ) sin () d d dr = (H j ) : jN Appendix B: Fundamental Solutions for the Construction of the Radial Green's Matrix 0 0 In order to determine the fundamental solutions (r; r ) and (r; r ) of the radial system (12) with 1 2 complex-valued frequencies, which are used for the construction of the Green's matrix de ned via equation (48), we rst study certain aspects of the associated Jost-type solutions [1, 21]. In more detail, we derive radial Jost-type equations that yield solutions with asymptotic behaviors near in nity, the event horizon, and the Cauchy horizon similar to those given in Lemmas II.1 and II.2, and brie y discuss the existence, uniqueness, and boundedness of these solutions. Since we apply Lebesgue's dominated convergence theorem to simplify our integral spectral representation of the Dirac propagator, the latter aspect becomes also relevant for the commutation of speci c limits, sums, and integrals. For the derivation of the Jost-type equations, we rewrite the radial rst-order system (12) as two second-order e e e scalar equations. In terms of the Regge{Wheeler coordinate r and the function R = (R ;R ) = ? + jjR ; r R , these read + + @ + J (r) @ + K (r) R = 0 ; (B1) r r ;k;! r ? ? ? ;k;! where 2 2 1 M (3r a ) im J (r) := r 4i!Mr 2ika ;k;! 2 2 2 2 r + a r + a + imr r M im + 2 2 2 K (r) := i! [ + 4Mr] + 2ika i! + + 2i! (r + M ) m r ;k;! 2 2 2 (r + a ) + imr im 2 2 2 2 K (r) := i! r M 4i!Mr 2ika + + ! m r : ;k;! 2 2 2 (r + a ) imr Employing the ansatzes R (r ) = exp K (r) dr Y (r ) ; ? ? ? ;k;! we may transform (B1) into the Schr odinger-type equations @ + V (r) Y = 0 (B2) r r ? ? ;k;! with the potentials J (r) @ J (r) ;k;! r ;k;! V (r) := K (r) : ;k;! ;k;! 4 2 To obtain Jost-type equations with boundary conditions prescribed at in nity, we split these potentials into an asymptotic contribution e ective at in nity and otherwise regular contributions V = V + V ; (B3) 1 reg. ;k;! where the asymptotic contribution is given by the expression 2Mm 2 2 V = V (r ) := ! m + (B4) 1 1 ? and the regular contributions are on the order of V = O 1=r satisfying the condition reg. ? V (y) dy < 1 for all r > 0 : reg. ? 27 We remark that the asymptotic potential (B4) corresponds to the equation @ + V (r ) Y = 0 ; (B5) r r 1 ? 1 ? ? which has the solution [25] p p 2 2 2 2 1 1 Y = Z W 2i sign(!) ! m r +Z W 2i sign(!) ! m r ; 1 1 ? 2 ? ; + ; 2 2 2 2 where W 1 ( ) are Whittaker functions with := i sign(!) Mm = ! m and Z denote con- 1=2 stants. The asymptotics of this solution at in nity reads in case j!j m Mm 2 2 Y Z exp i sign(!) ! m r + p ln (r ) 1 ? ? 2 2 ! m (B6) Mm 2 2 +Z exp i sign(!) ! m r + p ln (r ) ; ? ? 2 2 ! m whereas for j!j < m it yields Mm 2 2 Y Z exp sign(!) m ! r + ln (r ) 1 ? ? 2 2 m ! Mm 2 2 +Z exp sign(!) m ! r + p ln (r ) ? ? 2 2 m ! with Z also being constants (cf. Lemma II.1). In the following, we restrict our attention to the 1=2 case j!j m. The case j!j < m may be treated accordingly. As in the usual study of Jost equations and their solutions, we complexify the Schr odinger-type equations (B2) via the analytic continuation ! ! ! 2 C of the frequency. Then, by means of the above splittings of the potentials (B3) and the speci c form of the asymptotic solution (B6), we can write the Jost-type equation representation of (B2) as Mm 2 2 Y (r ) = exp i sign Im(! ) sign(! ) j! j m r + ln (r ) ? c c c ? ? 2 2 j! j m (B7) sin V (y) [r y] 1 ? + p V (y)Y (y) dy : reg. V (y) ? 1 We note that the proper complexi cation of the asymptotic Schr odinger-type equation (B5), which is in accordance with the particular representation (B6) of the asymptotic solutions containing signum functions, is obtained by rst rewriting the potential V de ned in (B4) in the form 2Mm 2 2 V = sign (!) j!j m + and subsequently extending the frequency ! to complex values. This is relevant for the derivation of the exponential term in (B7). Applying the series ansatzes Y (r ) = Y (r ) ; ? ;n ? n=0 where the zeroth-order terms are given by Mm 2 2 Y (r ) = exp i sign Im(! ) sign(! ) j! j m r + ln (r ) ; ;0 ? c c c ? ? 2 2 j! j m c 28 in (B7), we nd the recurrence relations sin V (y) [r y] 1 ? Y (r ) = p V (y)Y (r ) dy for n 1 : ;n ? reg. ;n 1 ? V (y) ? 1 In the theorem below, we discuss the relevant points pertaining to the existence, uniqueness, and boundedness of such solutions for the case Im(! ) < 0. Detailed proofs are worked out explicitly in, e.g., [13, 14, 17, 21]. The results and proofs for the case Im(! ) > 0 are in essence identical. Theorem B.1. For each ! 2 C with ! 6= 0 and Im(! ) < 0, the Jost-type equations (B7) have unique c c c solutions Y (r ) obeying Mm 2 2 lim exp i sign(! ) j! j m r + p ln (r ) Y (r ) < 1 : c c ? ? ? 2 2 r !1 j! j m These solutions are moreover continuously di erentiable in r on the interval (0;1) with Mm 2 2 lim exp i sign(! ) j! j m r + p ln (r ) Y (r ) = 1 c c ? ? ? 2 2 r !1 j! j m and p p Mm 2 2 2 2 lim exp i sign(! ) j! j m r + p ln (r ) @ Y (r ) = i sign(! ) j! j m : c c ? ? r ? c c 2 2 r !1 j! j m For each xed value of r , Y (r ) and @ Y (r ) are functions that are analytic in f! j Im(! ) < 0g, ? ? r ? c c continuous in f! j ! 6= 0 and Im(! ) < 0g, and satisfy the bound c c c Mm 2 2 Y (r ) exp i sign(! ) j! j m r + p ln (r ) ? c c ? ? 2 2 j! j m exp Im V (r ) r exp Q (r ) 1 1 ? ? ? as well as Mm 2 2 @ Y (r ) + exp i sign(! ) j! j m r + p ln (r ) r ? c c ? ? 2 2 j! j m p p Mm 2 2 i sign(! ) j! j m + p exp Im V (r ) r + Q (r ) V (y) dy ; c c 1 ? ? ? reg. 2 2 j! j m r c ? ? where p p 4y V (y) reg. Q (r ) := p exp Im V (y) + Im V (y) y dy : ? 1 1 1 + y V (y) ? 1 It remains to determine the Jost-type equations with boundary conditions prescribed at the event horizon and at the Cauchy horizon. This can be done using a similar approach as in the above case. For details, we again refer to [13, 14]. 0 0 We now specify the fundamental solutions (r; r ) and (r; r ) of the radial system (12) with 1 2 complex-valued frequencies. Due to the high degree of complexity of this system, explicit analytical expressions for its fundamental solutions are not known. Thus, we describe them in terms of suitable 29 asymptotic expansions. To this end, we de ne, on the one hand, auxiliary functions that have the proper decay at in nity (cf. Lemma II.1) 1 Im(! ) < 0 if j! j m (1) 1=2 i (r (r)) + ? (r) := jj d e +O for 1;1 r (r) ? Re(! ) 0 if j! j < m Im(! ) > 0 if j! j m 0 1 (1) i (r (r)) (r) := d e +O for 2;1 1 r (r) Re(! ) < 0 if j! j < m ; where the quantities d are scalar constants and the functions are given in (14), but with 1=2;1 p p 2 2 2 2 frequencies ! 2 f! + i; ! ig, for which > 0, and with the substitution ! m ! j! j m . On the other hand, we use auxiliary functions that are nite at the event and the Cauchy horizon and further comply with the associated asymptotics (cf. Lemma II.2) (+) (+) 1=2 2i ! +k r (r) qr (r) ? ? b Kerr (r) := jj d e +O e for Im(! ) < 0 1;r (+) qr (r) (r) := d +O e for Im(! ) > 0 2;r ( ) qr (r) (r) := d +O e for Im(! ) < 0 1;r ( ) ( ) 1=2 2i ! +k r (r) qr (r) ? ? q Kerr (r) := jj d e +O e for Im(! ) > 0 ; 2;r where d are scalar constants as well. To clarify the notation, we point out that the superscripts 1=2;r (1), (+), and ( ) designate asymptotic expansions at in nity, the event horizon, and the Cauchy horizon, respectively. Besides, we remark that the existence and uniqueness of fundamental solutions of the radial system (12) with these particular asymptotic expansions follows from the above study of the radial Jost-type equations and, moreover, that these asymptotic expansions ensure that the fundamental solutions are square-integrable. For example, as the Regge{Wheeler coordinate r tends (+) to minus in nity at the event horizon, the exponential factor in the auxiliary function tends to zero because Im(! ) < 0. However, this exponential factor would not be square-integrable if Im(! ) > 0. Last, we introduce an auxiliary function that satis es the Dirichlet-type boundary condition (40) at r = r (1) (r) := (r) ; @M @M i jj = r (1) with denoting its rst component. Then, in case @M j! j m and Im(! ) < 0 or j! j < m ; Im(! ) < 0 ; and Re(! ) 0 ; the fundamental radial solutions and read 1 2 30 0 0 (1) (r; r < r < 1) = (r r ) (r) 1 + 0 0 (+) (r; r < r < 1) = (r r) (r r ) (r) 2 + + 0 0 ( ) (1) b b (r; r < r r ) = (r r ) (r r) (r) + (r r ) (r) 1 + + + (B8) 0 0 (+) (r; r < r r ) = (r r ) (r r) (r) 2 + + 0 0 ( ) (1) b b (r; r r r ) = (r r ) (r r) (r) + (r r ) (r) 1 0 + + 0 0 (r; r r r ) = (r r) (r) ; 2 0 @M whereas for j! j m and Im(! ) > 0 or j! j < m ; Im(! ) > 0 ; and Re(! ) < 0 ; they are given by 0 0 (1) (r; r < r < 1) = (r r ) (r) 1 + 0 0 (+) (r; r < r < 1) = (r r) (r r ) (r) + (r r) (r) 2 + @M 0 0 (+) (r; r < r r ) = (r r) (r r ) (r) + (r r) (r) 1 + @M (B9) 0 0 ( ) (r; r < r r ) = (r r) (r r ) (r) 2 + 0 0 ( ) (r; r r r ) = (r r ) (r r) (r) 1 0 0 0 (r; r r r ) = (r r) (r) : 2 0 @M In the remaining cases j! j < m ; Im(! ) < 0 ; and Re(! ) < 0 or j! j < m ; Im(! ) > 0 ; and Re(! ) 0 ; we also obtain the fundamental solutions (B8) or (B9), respectively, but with the auxiliary functions (1) (1) b q and interchanged. A case-by-case analysis shows that these solutions are uniquely determined by the conditions and asymptotics listed in the proof of Theorem VI.1. [1] V. de Alfaro and T. Regge, \Potential scattering," North-Holland Publishing Company (1965). [2] R. A. Bartnik and P. T. Chru sciel, \Boundary value problems for Dirac-type equations," arXiv:math/ 0307278 [math.DG], Journal fur die Reine und Angewandte Mathematik 579, 13 (2005). [3] R. H. Boyer and R. W. Lindquist, \Maximal analytic extension of the Kerr metric," Journal of Mathematical Physics 8, 265 (1967). [4] B. Carter, \Black hole equilibrium states," in Black holes/Les astres occlus, Ecole d' et e Phys. Th eor., Les Houches (1972). [5] S. Chandrasekhar, \The mathematical theory of black holes," Oxford University Press (1983). [6] P. R. Cherno , \Essential self-adjointness of powers of generators of hyperbolic equations," Journal of Functional Analysis 12, 401 (1973). [7] A. Chodos, R. L. Ja e, K. Johnson, C. B. Thorn, and V. F. Weisskopf, \New extended model of hadrons," Physical Review D 9, 3471 (1974). [8] A. S. Eddington, \A comparison of Whitehead's and Einstein's formulae," Nature 113, 192 (1924). 31 [9] D. Finkelstein, \Past-future asymmetry of the gravitational eld of a point particle," Physical Review 110, 965 (1958). [10] F. Finster, N. Kamran, J. Smoller, and S. T. Yau, \Decay rates and probability estimates for massive Dirac particles in the Kerr{Newman black hole geometry," arXiv:gr-qc/0107094, Communications in Mathemat- ical Physics 230, 201 (2002). [11] F. Finster, N. Kamran, J. Smoller, and S. T. Yau, \The long-time dynamics of Dirac particles in the Kerr{ Newman black hole geometry," arXiv:gr-qc/0005088, Advances in Theoretical and Mathematical Physics 7, 25 (2003). [12] F. Finster, N. Kamran, J. Smoller, and S. T. Yau, \An integral spectral representation of the propagator for the wave equation in the Kerr geometry," arXiv:gr-qc/0310024, Communications in Mathematical Physics 260, 257 (2005). [13] F. Finster, N. Kamran, J. Smoller, and S. T. Yau, \Decay of solutions of the wave equation in the Kerr geometry," arXiv:gr-qc/0504047, Communications in Mathematical Physics 264, 465 (2006). [14] F. Finster and J. Smoller, \Decay of solutions of the Teukolsky equation for higher spin in the Schwarzschild geometry," arXiv:gr-qc/0607046, Advances in Theoretical and Mathematical Physics 13, 71 (2009). [15] F. Finster and C. R oken, \Self-adjointness of the Dirac Hamiltonian for a class of non-uniformly elliptic boundary value problems," arXiv:1512.00761 [math-ph], Annals of Mathematical Sciences and Applications 1, 301 (2016). [16] F. John, \Partial di erential equations," Springer-Verlag (1991). [17] J. Kronthaler, \The Cauchy problem for the wave equation in the Schwarzschild geometry," arXiv:gr- qc/0601131, Journal of Mathematical Physics 47, id. 042501 (2006). [18] B. O'Neill, \The geometry of Kerr black holes," Dover Publications (2014). [19] J. Rauch and M. Taylor, \Essential self-adjointness of powers of generators of hyperbolic mixed problems," Journal of Functional Analysis 12, 491 (1973). [20] M. Reed and B. Simon, \Methods of modern mathematical physics I: functional analysis," Academic Press (1980). [21] M. Reed and B. Simon, \Methods of modern mathematical physics III: scattering theory," Academic Press (1979). [22] C. R oken, \The massive Dirac equation in Kerr geometry: separability in Eddington{Finkelstein-type coordinates and asymptotics," arXiv:1506.08038 [gr-qc], General Relativity and Gravitation 49, 39 (2017). [23] M. Taylor, \Partial di erential equations I," Springer-Verlag (1996). [24] M. Taylor, \Partial di erential equations III," Springer-Verlag (1997). [25] E. T. Whittaker and G. N. Watson, \A course of modern analysis," Cambridge University Press (1927).
Journal
Mathematical Physics – arXiv (Cornell University)
Published: Jun 5, 2016