arxiv: v1 [gr-qc] 31 Mar 2009

Transcription

1 axiv: v1 [g-qc] 31 Ma 2009 Hamonic Gauge Petubations of the Schwazschild Metic by Mak V. Bendtson B.A., Univesity of Coloado, 1996 A thesis submitted to the Faculty of the Gaduate School of the Univesity of Coloado in patial fulfillment of the equiements fo the degee of Docto of Philosophy Depatment of Physics 2007

2 This thesis entitled: Hamonic Gauge Petubations of the Schwazschild Metic witten by Mak V. Bendtson has been appoved fo the Depatment of Physics Neil Ashby David Batlett Date The final copy of this thesis has been examined by the signatoies, and we find that both the content and the fom meet acceptable pesentation standads of scholaly wok in the above mentioned discipline.

3 iii Bendtson, Mak V. Ph.D., Physics) Hamonic Gauge Petubations of the Schwazschild Metic Thesis diected by Pof. Neil Ashby The satellite obsevatoy LISA will be capable of detecting gavitational waves fom exteme mass atio inspials EMRIs), such as a small black hole obiting a supemassive black hole. The gavitational effects of the much smalle mass can be teated as the petubation of a known backgound metic, hee the Schwazschild metic. The petubed Einstein field equations fom a system of ten coupled patial diffeential equations. We solve the equations in the hamonic gauge, also called the Loentz gauge o Loenz gauge. Using sepaation of vaiables and Fouie tansfoms, we wite the solutions in tems of six adial functions which satisfy decoupled odinay diffeential equations. We then use the solutions to calculate the gavitational self-foce fo cicula obits. The self-foce gives the fist ode petubative coections to the equations of motion.

4 Dedication To my family

5 v Acknowledgements I thank my adviso, Neil Ashby, fo his advice and patience, and the othe committee membes fo thei sevice.

6 Contents Chapte 1 Intoduction Black Hole Petubation Theoy Summay of Thesis Odd Paity Solutions Non-Zeo Fequency Solutions Solutions fo l Solutions fo l = Zeo Fequency Solutions Solutions fo l Solutions fo l = Homogeneous Solutions Even Paity Solutions Non-Zeo Fequency Solutions Solutions fo l Solutions fo l = 0, Zeo Fequency Solutions Solutions fo l Solution of Two Systems of Equations

7 vii Solutions fo l = Solutions fo l = Homogeneous Solutions Inteim Summay of Odd and Even Paity Solutions Equations of Motion Backgound Geodesic Equations Equations fo Gavitational Self-Foce Calculation of the Stess Enegy Tenso fo a Point Mass Multipole Decomposition Fouie Tansfoms Solution of Genealized Regge-Wheele Equations Non-Zeo Fequency Homogeneous Solutions Zeo Fequency Homogeneous Solutions Inhomogeneous Solutions Radiation Wavefoms Enegy and Angula Momentum Flux Numeical Results Conclusion 205

8 Bibliogaphy 207 viii Appendix A Non-Zeo Fequency Even Paity Solutions fo l B Five Zeo Fequency Even Paity Solutions fo l C Zeo Fequency Even Paity Solution fo l = D Zeo Fequency Even Paity Solutions fo l = 0 233

9 ix Tables Table 3.1 Numeical Values of the Constant c 3 in Equation 3.148) Numeical Values of the Constant c 3 in Equation 3.193) Radial Component of Self-Foce fo Cicula Obits Compaison of Radiation Reaction Self-Foce with Gavitational Wave Enegy Flux Cicula Obit Enegy Flux Compaison fo R = 10M Radial Component of Self-Foce fo Cicula Obits Using Detweile and Poisson Fomula Radial Component of Self-Foce Compaison

10 x Figues Figue 4.1 Schematic of Self-Foce fo Cicula Obit Plot of Radial Self-Foce fo Cicula Obits Sample Plot of Least Squaes Fit fo R = 10M Compaison Plot of Radial Self-Foce fo Cicula Obits

11 Chapte 1 Intoduction The satellite obsevatoy LISA will be capable of detecting gavitational waves fom exteme mass atio inspials EMRIs). These occu when a compact sta, such as a black hole, neuton sta, o white dwaf, is captued by a supemassive black hole [80]. The evolution of EMRI obits and thei gavitational wavefoms can be calculated using black hole petubation theoy. An intoduction to black hole petubation theoy, including the hamonic gauge, is in section 1.1. An outline of the emainde of the thesis is in section Black Hole Petubation Theoy The Schwazschild metic is ds 2 = g µν dx µ dx ν = 1 2M ) dt M )d 2 2 dθ 2 2 sin 2 θdφ 2, 1.1) whee the standad coodinates have been used. It is a solution of the Einstein field equations, which ae R µν 1 2 g µνr = 8π G c 4 T µν. 1.2) The Schwazschild metic is a vacuum solution, meaning T µν = 0. In 1.2), the gavitational constant G and speed of light c ae shown explicitly. Geneally though, we will use geometized units, fo which G = c = 1. We will use this and othe notational

12 2 conventions as descibed by Misne, Thone and Wheele [71], including the metic signatue of 1.1) and fom of the field equations 1.2). Black hole petubation theoy fo the Schwazschild metic was fomulated by Regge and Wheele [98] and extended by Zeilli [115]. A summay of thei method follows, which is taken mainly fom thei aticles. A small petubation h µν is added to the backgound Schwazschild metic g µν. Ou physical poblem involves a small mass m 0 obiting a much lage black hole M, so h µν is popotional to the mass atio m 0 /M. Accodingly, the total petubed metic g µν is g µν = g µν h µν Om 0 /M) 2 ). 1.3) The invese petubed metic is g µν = g µν h µν Om 0 /M) 2 ). 1.4) The petubed field equations ae linea in h µν : ] [h ;α µν;α 2R α µ β ;α ;α ν h αβ h µα ;ν h να ;µ ) h ;µ;ν R α νh µα R α µh να g µν h λα ;α;λ h ;λ ;λ ) h µν R g µν h αβ R αβ = 16πT µν. 1.5) The semicolons epesent covaiant diffeentiation with espect to the backgound metic g µν. Fo the Schwazschild metic, the Ricci tenso R µν and Ricci scala R ae zeo, so the field equations simplify to ] [h ;α µν;α 2R α µ β ;α ;α ν h αβ h µα ;ν h να ;µ ) h ;µ;ν g µν h λα ;α;λ h ;λ ;λ ) = 16πT µν. 1.6) As discussed in [71], the left side of the petubed equations epesents the popagation of a wave inteacting with the backgound spacetime cuvatue. The stess enegy tenso T µν is the covaiant fom of T µν = m 0 δ 4 x zτ)) g dz dτ µ ν dz dτ, 1.7) dτ

13 3 whee the delta function epesents a point mass o test paticle. The stess enegy tenso divegence equation is T µν ;ν = 0, 1.8) which, when applied to 1.7), gives the backgound geodesic equations of motion [32]. Regge, Wheele and Zeilli solved the petubed field equations using sepaation of vaiables. The angula dependence is contained in tenso hamonics, which ae obtained fom the familia spheical hamonics. A Fouie tansfom sepaates the time dependence, leaving a set of adial odinay diffeential equations to be solved. The solution was simplified by making a paticula choice of gauge, the so-called Regge- Wheele gauge. A gauge is a choice of coodinates. A change of gauge is a small change of coodinates x µ new = xµ old ξµ, 1.9) which causes the metic petubation to change as h new µν = h old µν ξ µ;ν ξ ν;µ. 1.10) Hee, small means Om 0 /M). Although the petubation is gauge dependent, the petubed field equations in 1.6) ae gauge invaiant [32]. We will use the notation of Ashby [4], which is diffeent fom, and simple than, Zeilli s notation [115]. Because the backgound metic is spheically symmetic, the petubation can be split into odd and even paity pats. This decomposition gives h µν t,,θ,φ) = l l=0 m= l ) e iωt h o,lm µν ω,,θ,φ) h e,lm µν ω,,θ,φ) dω. 1.11) The odd paity tems ae 0 0 h lm 0 ω,)csc θ Y lmθ,φ) φ h lm 0 ω,)sin θ Y lmθ,φ) θ h o,lm µν ω,,θ,φ) = 0 h lm 1 ω,)csc θ Y lmθ,φ) φ h lm 1 ω,)sin θ Y lmθ,φ) θ h lm 2 ω,)x lmθ,φ) h lm 2 ω,)sin θw, 1.12) lmθ,φ) h lm 2 ω,)sin2 θx lm θ,φ)

14 4 whee W lm θ,φ) = 2 Y lm θ,φ) θ 2 The even paity pat is X lm θ,φ) = 2 sin θ 1 2 Y lm θ,φ) sin 2 θ ) cot θ Y lmθ,φ) θ Ylm θ,φ) cot θ Y lm θ,φ) φ θ φ 2, 1.13). 1.14) h e,lm µν ω,,θ,φ) = 1 2M ) H lm 0 ω,)y lm H1 lmω,)y lm h lm 0 ω,) Y lm θ h lm 0 ω,) Y lm φ H2 lmω,)y lm h 1 2M ) lm 1 ω,) Y lm θ h lm 1 ω,) Y lm φ 2 K lm ω,)y lm ) 2 sinθ G lm ω,)x lm G lm ω,)w lm 2 sin 2 θ K lm ω,)y lm G lm ω,)w lm ). 1.15) Asteisks epesent symmetic components. The angula functions ae the tenso hamonics. The Regge-Wheele gauge is defined by setting fou adial factos equal to zeo: the odd paity h 2 and the even paity h 0, h 1 and G. The tace h of the metic petubation is ht,,θ,φ) = l l=0 m= l Y lm θ,φ) e iωt h lm ω,)dω, 1.16) whee h lm ω,) = H lm 0 ω,) H lm 2 ω,) 2K lm ω,). 1.17) Similaly, the stess enegy tenso may be decomposed in tems of tenso hamonics and Fouie tansfoms. The covaiant components ae T µν t,,θ,φ) = l l=0 m= l ) e iωt Tµν o,lm ω,,θ,φ) Tµν e,lm ω,,θ,φ) dω, 1.18)

15 5 whee 0 0 So lm 02 ω,)csc θ Y lm φ So lm 02 ω,)sin θ Y lm θ Tµν o,lm 0 So lm 12 ω,,θ,φ) = ω,)csc θ Y lm φ So lm 12 ω,)sin θ Y lm θ So lm 22 ω,)x lm So lm 22 ω,)sin θw lm So lm 22 ω,)sin2 θx lm, 1.19) and T e,lm µν ω,,θ,φ) = Se lm 00 ω,)y lm Se lm 01 ω,)y lm Se lm 02 ω,) Y lm θ Se lm 11 ω,)y lm Se lm 12 ω,) Y lm θ Ue lm 22 ω,)y lm Se lm 22 ω,)w lm The tace of the stess enegy tenso is Se lm 02 ω,) Y lm φ Se lm 12 ω,) Y lm φ Se lm 22 ω,)sin θx lm sin 2 θ Ue lm 22 ω,)y lm) Se lm 22 ω,)w lm. 1.20) T = g µν T µν, 1.21) and its multipole decomposition is Tt,,θ,φ) = l l=0 m= l Y lm θ,φ) e iωt T lm ω,)dω, 1.22) whee T lm ω,) = 2M Selm 00 ω,) 1 2M ) Se lm 11 ω,) 2 2Uelm 22 ω,). 1.23) The gauge change vecto ξ µ is ξ µ t,,θ,φ) = l l=0 m= l {[ ξ e,lm µ ω,,θ,φ) e iωt ξµ o,lm ω,,θ,φ) ) ] dω δ l0 δ µt C 0 1 2M ) t Y 00 θ,φ) [ δ l1 C 1 t 2 δ µθ csc θ Y 1mθ,φ) δ µφ sinθ Y 1mθ,φ) φ θ ]}, 1.24)

16 6 whee δ ab is the Konecke delta and whee ) ξµ o,lm ω,,θ,φ) = 0,0,Z lm ω,)csc θ Y lm φ, Zlm ω,)sin θ Y lm θ, 1.25) ) ξµ e,lm ω,,θ,φ) = M0 lmω,)y lm,m1 lmω,)y lm,m2 lmω,) Y lm θ,mlm 2 ω,) Y lm φ. 1.26) Because a gauge change is a small change of coodinates, the quantities Z lm ω,) and M lm i ω,) ae of ode m 0 M, just as the petubation h µν is. The emaining two tems in ξ µ 1.24) ae also ode m 0 M and ae discussed in subsections and 2.2.2, espectively. Regge and Wheele showed that the odd paity components could be expessed in tems of a single adial scala function, the Regge-Wheele function, which satisfies a second ode odinay diffeential equation. Zeilli did the same fo the even paity components; his function is called the Zeilli function. Late, Moncief showed that these two functions ae gauge invaiant [72]. The Regge-Wheele gauge is compaatively simple, but we will wok in a diffeent gauge, the hamonic gauge. The eason fo this choice is that the equations fo the gavitational self-foce wee deived in the hamonic gauge [70], [96]. The basic hamonic gauge field equations ae given in [71], which uses the tem Loentz gauge. The spelling Loenz is also used in the liteatue [86]. We fist define the tace-evesed metic h µν = h µν 1 2 g µνh, 1.27) whee the tace h is h = g αβ h αβ. 1.28) The hamonic gauge is defined by the condition h µν ;ν = ) Using 1.27) and 1.29), the field equations 1.6) simplify to h µν;α ;α 2R α µ β ν h αβ = 16πT µν. 1.30)

17 7 Adding 1.6) and 1.30) gives h µα ;α ;ν h να ;α ;µ g µν h λα ;α;λ = 0, 1.31) which ae the tems eliminated in going fom 1.6) to 1.30). The hamonic gauge condition 1.29) and field equations 1.30) ae peseved afte a gauge change which satisfies ξ ;ν µ;ν = ) The hamonic gauge equations above apply when the backgound metic descibes a cuved spacetime, as the Schwazschild metic does. Befoe discussing futhe the hamonic gauge fo the Schwazschild metic, it is helpful to eview a much simple example, the plane wave. The analysis below of the plane wave is taken fom Weinbeg s chapte on gavitational adiation [110], supplemented by [104]. The plane wave is a petubation h µν of the flat space metic η µν. The petubed metic g µν is g µν = η µν h µν. 1.33) We can choose the gauge so that x µhµ ν = 1 2 x ν hµ µ. 1.34) In this gauge, the homogeneous field equations ae h µν = 0, 1.35) whee is the flat spacetime D Alembetian opeato. The gauge 1.34) and field equations 1.35) ae peseved by gauge changes ξ ν which satisfy ξ ν = ) Equations 1.34)-1.36) ae the flat spacetime backgound metic equivalents of the cuved backgound spacetime equations 1.29)-1.30), 1.32). Moeove, we can de-

18 8 ive equations 1.34)-1.36) fom the cuved spacetime expessions, although Weinbeg deives his equations ab initio. Weinbeg shows that the solution to the field equations 1.35) is a wave h µν x) = e µν expik λ x λ ) e µν exp ik λ x λ ), 1.37) whee e µν is the symmetic polaization tenso. Hee and elsewhee in this thesis, an oveba usually epesents complex conjugation; howeve, this notational ule does not apply to h µν 1.27) and its tace, h. Substituting h µν 1.37) into the field equations 1.35) and gauge definition 1.34) yields k µ k µ = 0, k µ e µ ν = 1 2 k νe µ µ, 1.38) espectively. A symmetic 4 4 matix has at most ten independent components. The gauge definition 1.34) epesents fou constaints, which educe the numbe of independent components, o polaizations, to six. The gauge tansfomation vecto appopiate to h µν 1.37) is ξ µ x) = ib µ expik λ x λ ) ib µ exp ik λ x λ ), 1.39) whee B µ is a constant vecto equation 9.14) of [104]). A change of gauge modifies the polaization tensos as The gauge tansfomation vecto 1.39) satisfies 1.36). vecto k λ is e new µν = e old µν k µ ξ ν k ν ξ µ. 1.40) The solution 1.37) is a plane wave taveling in the positive z-diection if the wave k x = k y = 0, k z = k t > ) Fo such a plane wave, the six independent polaizations can be witten as the following linea combinations of e µν : e ± = e xx ie xy, f ± = e zx ie zy, e tt, e zz. 1.42)

19 The emaining components ae not independent and depend on the six above, by symmety and by the following equations: 9 e tx = e zx, e ty = e zy, e tz = 1 2 e tt e zz ), e yy = e xx. 1.43) Weinbeg shows the components e ± ae gauge invaiant, because they can not be emoved by a coodinate tansfomation. Howeve, f ±, e tt and e zz can be eliminated by a coodinate tansfomation which satisfies 1.36) and peseves the gauge condition 1.34). The six independent components 1.42) behave diffeently when a otation is made about the z-axis, the diection of wave popagation. Specifically, Weinbeg shows e ± = exp±2iθ)e ±, f ± = exp±iθ)f ±, e tt = e tt, e zz = e zz. 1.44) whee the pimes denote polaizations afte otation though an angle θ. A plane wave ψ has helicity h if ψ = expihθ)ψ. 1.45) The otation esults 1.44) show that the gavitational plane wave h µν 1.37) can be decomposed into six pieces: two of helicity ±2, two of helicity ±1 and two of helicity 0. The tace of the polaization tenso is e zz e tt, so the zeo helicity pieces ae elated to the tace of h µν. Only the helicity ±2 components ae physically meaningful, because they alone can not be emoved by a coodinate tansfomation. Schwazschild metic petubation theoy has an analogue to the diffeent helicity functions of the plane wave example. The genealized Regge-Wheele equation is d 2 ψ s ) d 2 ω 2 ψ s ) V sl )ψ s ) = S slm ω, ), s = 0,1,2, 1.46) whee the potential V is V sl ) = 1 2M ) 2λ 1) 2 1 s 2 ) 2M 3 ). 1.47)

20 10 The genealized Regge-Wheele equation is discussed in [55], [63], [64]. It epesents a wave inteacting with an effective potential that esults fom the backgound spacetime cuvatue. The paamete s is the spin o spin weight, and it coesponds to the diffeent helicities of the plane wave example. The case ψ 2 is equal to the Regge-Wheele function that was deived in the Regge-Wheele gauge. Descibed in [98] and [115], the coodinate is so that d = d 1 2M ) ), = 2M ln 2M 1 d = 1 2M d ) d d, 1.48) 1.49) and d 2 d 2 = 1 2M ) 2 d 2 d 2 1 2M ) 2M 2 d d. 1.50) Because 2M < <, we have < <. Due to this elation, is called the totoise coodinate [71]. Thoughout this thesis, the paamete λ is defined as [115] λ = 1 l 1)l 2), 1.51) 2 which implies 2λ 1) = ll 1). 1.52) The genealized Regge-Wheele equation will often be abbeviated as L s ψ s = S s, 1.53) whee the opeato L s is defined by the left side of 1.46). The souce S s is constucted fom the adial coefficients of the stess enegy tenso 1.18)-1.20). The solutions fo the hamonic gauge ae descibed in tems of genealized Regge- Wheele functions. Specifically, the odd paity solutions ae witten in tems of two genealized Regge-Wheele functions, one with s = 2 and the othe with s = 1. The even paity solutions ae witten in tems of two functions with s = 0, one with s = 1

21 11 and one with s = 2. The even paity spin 2 function is actually the Zeilli function, but can be elated to the spin 2 Regge-Wheele function by diffeential opeatos. The even and odd paity spin 2 functions ae gauge invaiant. These six functions coespond to the six diffeent helicity states of the plane wave example. 1.2 Summay of Thesis Most of this thesis discusses the solution of the field equations in the hamonic gauge 1.30). Chaptes 2 and 3 summaize the deivation of the odd and even paity solutions, espectively. The solutions ae expessions fo the adial coefficients such as h 0 and H 0 ) of the angula functions in the odd 1.12) and even 1.15) paity metic petubations h µν. The cases of non-zeo and zeo angula fequency ae handled sepaately. Both inhomogeneous and homogeneous solutions ae coveed. Some of the even paity solutions ae listed sepaately in the appendices. With a few cited exceptions, the solutions deived in these two chaptes ae new and constitute the main eseach esults of this thesis. Chapte 4 coves the equations of motion. Both the backgound geodesic and gavitational self-foce equations ae pesented. Chapte 5 discusses the stess enegy tenso fo a point mass, its multipole decomposition and its Fouie tansfoms fo cicula and elliptic obits. The chapte includes fomulae fo the adial coefficients such as So 22 and Se 22 ) of the angula functions of the stess enegy components T µν in equations 1.19) and 1.20). Chapte 6 has homogeneous and inhomogeneous solutions fo the genealized Regge-Wheele and Zeilli equations. Chapte 7 explains how to use the hamonic gauge solutions to calculate gavitational wavefoms and enegy and angula momentum fluxes. Chapte 8 contains esults of numeical calculations, mainly of the gavitational self-foce fo cicula obits. Chapte 9 is a bief conclusion. Mathematica was used extensively fo the deivations. Some esults ae attibuted to unpublished wok of Neil Ashby, which is cited as efeence [4]. Among othe things,

22 12 he edeived and coected the published solutions fo the Regge-Wheele gauge. He also povided Mathematica tools fo simplifying expessions including deivatives of angula functions) and fo deiving ecusion elations fo infinite seies. Quantities elated to the backgound Schwazschild metic, such as Chistoffel symbols and Riemann cuvatue tenso components, wee calculated with Mathematica-based softwae witten by him.

23 Chapte 2 Odd Paity Solutions This chapte contains the deivation of the odd paity field equations and thei solutions in the hamonic gauge. We will use sepaation of vaiables to educe the petubed field equations to a system of thee coupled odinay diffeential equations, with independent vaiable. The esulting equations ae solved in tems of genealized Regge-Wheele functions, with s = 2 and s = 1. Section 2.1 descibes the non-zeo fequency solutions, with the cases of spheical hamonic index l 2 and l = 1 handled sepaately. Section 2.2 does the same fo zeo fequency. Section 2.3 concludes with a discussion of homogeneous solutions. 2.1 Non-Zeo Fequency Solutions Subsection coves the case l 2. Subsection shows how the solutions fo l = 1 can be obtained fom the solutions fo l Solutions fo l 2 The fist step is to deive the adial field equations, using sepaation of vaiables. The method of deivation is simila to that used by Regge, Wheele and Zeilli fo the Regge-Wheele gauge [98], [115]. We substitute the odd paity metic petubations fom 1.11)-1.12) and stess enegy tenso components fom 1.18)-1.19) into the hamonic gauge field equations 1.30), which ae a system of coupled patial diffeential

24 equations. Because the equations ae Fouie tansfoms, the time dependence is in the exponentials e iωt, which divide off. Patial deivatives with espect to time become 14 factos of angula fequency, using the ule t iω. The angula vaiables θ and φ ae contained entiely in the tenso hamonics and thei deivatives. Afte simplifying the angula deivatives, the tenso hamonics also sepaate off, and we ae left with the following adial factos to the odd paity field equations: 1 2M ) 2 h 0 8M2 42 λ)m 2 2 2λ iω) 2 2 )) 4 h 0 2iωM2M ) 3 h 1 = 16π 1 2M ) So 02, 2.1) 1 2M ) 2 h 1 4M 2 1 2M ) h 1 2iωM λ8m 4) 2M 2h 0 16M2 45 λ)m 2 6 2λ iω) 2 2 ) 4 h 1 = 16π 4 h 2 1 2M ) So 12, 2.2) 1 2M ) 2 h 22M )3M ) 2 3 h 2 2M)2 2 3 h 1 16M2 4 3 λ)m 2 2 2λ iω) 2 2 ) 4 h 2 = 16π 1 2M ) So ) Pimes signify diffeentiation with espect to. Geneally, functional dependence on l,m,ω and will be suppessed, so that, fo example, h 0 abbeviates h lm 0 ω,) and So 02 epesents So lm 02 ω,). Howeve, the system of equations 2.1)-2.3) must be solved fo each combination of indices l, m and ω. Equation 2.1) comes fom the tφ component of the field equations multiplied by a facto of 1 2M/)). Similaly, 2.2) and 2.3) ae obtained fom the φ and φφ components, espectively. Because thee ae only thee odd paity adial functions h 0, h 1 and h 2 ), thee ae only thee odd paity adial field equations. The emaining odd paity components of the field equations ae eithe zeo o have the same adial factos as 2.1)-2.3).

25 Each of the thee adial field equations can be witten in the fom 1 2M 15 ) 2 h i iω)2 h i othe tems, i = 0,1,2, 2.4) o, altenatively, d 2 h i d 2 iω) 2 h i othe tems, 2.5) whee the othe tems ae at most fist ode deivatives. The invese Fouie tansfom of 2.5) is 2 h i 2 2 h i othe tems, 2.6) t2 which is a wave equation. In the time domain, the field equations ae a hypebolic system of patial diffeential equations. The system is a well-posed initial value poblem and can be solved numeically in the time domain [7], [70]. The field equations above agee with those deived by Baack and Lousto [7], who used diffeent notation and woked in the time domain. Related equations ae also sepaable. The hamonic gauge condition h µν ;ν = 0 leads to the adial equation 1 2M ) h 1 iω 2M h 0 The stess enegy tenso divegence equation, T µν ;ν = 0, is [4] 1 2M ) So 12 iω 2M So 02 2M ) 2 h 1 2λ 2 h 2 = ) 2M ) 2 So 12 2λ 2 So 22 = ) Fo a change of gauge descibed by equations 1.10) and 1.24)-1.25), the adial petubation factos tansfom as h new 0 = h old 0 iωz, 2.9) h new 1 = h old 1 2 Z Z, 2.10) h new 2 = h old 2 Z, 2.11)

26 16 whee the function Z is ode m 0 M [4], [98], [115]. To peseve the hamonic gauge, a change of gauge must satisfy ξ µ;ν ;ν = 0, and the associated adial equation is d 2 Z d 2 ω 2 Z 1 2M ) 2λ 1) 2 Z = ) This is the homogeneous genealized Regge-Wheele equation, with s = 1. Odd paity gauge changes which peseve the hamonic gauge ae implemented by adding homogeneous spin 1 solutions to the metic petubations. To deive anothe fist ode equation, diffeentiate 2.7) with espect to and use 2.2) and 2.7) to eliminate h 1 and h 1. This gives iωh 0 2λ 2 1 2M ) h 2 2iω h 0 4λ 3 1 2M iω) 2 λ4m 2) 3 ) h 2 = 16π ) h 1 1 2M ) So ) Unlike the hamonic gauge condition 2.7), which applies only in the hamonic gauge, equation 2.13) is gauge invaiant and applies in any gauge. As shown by 2.9)-2.11), the adial functions ae gauge dependent; howeve, a gauge change to one say h 0 ) is canceled by changes to the emaining functions, leaving equation 2.13) the same in the new gauge. The stess enegy tenso tem So 12 is coodinate dependent; howeve, it is ode m 0 M, so any changes to it ae ode m 0 ) 2. M Thus, 2.13) is gauge invaiant to linea ode in m 0 M. In contast, a gauge change applied to 2.7) esults in additional tems involving Z and its deivatives, unless the change satisfies 2.12) and peseves the hamonic gauge. Equation 2.7) is invaiant only unde changes which peseve the hamonic gauge, but equation 2.13) is invaiant unde abitay gauge changes. Anothe way of deiving 2.13) is to substitute the odd paity metic petubations into 1.6), the geneal petubation field equations applicable to any gauge. Equation 2.13) is obtained fom the φ component, multiplied by 1 2M ). Because 1.6) is gauge invaiant, so is 2.13). Regge and Wheele showed that, in the Regge-Wheele gauge, the odd paity

27 petubations can be witten in tems of a single scala function which satisfies a wave equation called the Regge-Wheele equation [98]. In the notation of this thesis, this function is the odd paity ψ 2, and the Regge-Wheele equation is d 2 ψ 2 d 2 ω 2 ψ 2 1 2M ) 2λ 1) 2 6M 3 17 ) ψ 2 = S ) Equation 2.14) is the genealized Regge-Wheele equation fom 1.46), with s = 2. Subsequently, Moncief poved that the Regge-Wheele function ψ 2 is gauge invaiant [72]-[73]. He deduced a fomula fo ψ 2 in tems of the metic petubations. Adjusting fo diffeences in notation between his pape and this thesis, his expession is ψ 2 = 1 2M ) h1 2h ) 2 2 h ) Although the adial petubation functions ae gauge dependent, ψ 2 is gauge invaiant. Moncief showed that any gauge change in h 1 would be offset by changes to h 2 and h 2, as follows. Fom equation 2.10), the tem containing h 1 changes by 2Z 2 Z. Using 2.11), the tems with h 2 and h 2 change by 2Z Z 2. The changes cancel each othe, leaving ψ 2 invaiant. Moncief did not deive his esult in the hamonic gauge, but his wok can be used hee, because 2.15) is gauge invaiant. In equation 2.14), the quantity S 2 is a souce tem constucted fom the adial components of the stess enegy tenso. To find S 2, substitute 2.15) into the Regge- Wheele equation and simplify with the othe hamonic gauge equations, obtaining S 2 = 16π 1 2M ) 2 So 12 32π6M2 5M 2 ) 4 So 22 16π 1 2M ) 2 So ) Taking into account diffeences in notation, equation 2.16) agees with Zeilli s esult in the Regge-Wheele gauge [115], as coected by othes [4], [102]. h 1, we find Solutions fo h 1 and h 0 can be witten tems of h 2 and ψ 2. By solving 2.15) fo h 1 = 2 h 2 2 2M ψ 2 h )

28 We then use 2.17) and its adial deivative to eliminate h 1 and h 1 fom 2.7) and solve 18 fo h 0 to get h 0 = 2M ) ψ 2 iω 1 2M ) ψ2 iω iωh 2 1 2M ) 16π iω So ) Substituting fo h 1 in the field equation 2.3), we have L 1 h 2 = 2 1 2M ) ψ 2 16π 1 2M ) So 22, 2.19) whee the opeato L 1 was defined in 1.53). The left side is the genealized Regge- Wheele equation, with s = 1. To complete the odd paity solutions, we will solve 2.19) fo h 2 and substitute the esult into 2.18) and 2.17). The fom of 2.19) suggests the following tial solution h ty 2 = e d C ψ 2 ψ 1 C e C ψ 2 16π f 02 So 02 f 12 So 12 f 22 So 22 ), 2.20) whee C is a constant and the othe quantities ae functions of. To find the unknowns, we inset h ty 2 into the left side of 2.19) and obtain L 1 ψ 1 = ψ 2,ψ 2 tems souce tems. 2.21) The souce tems ae complicated expessions involving So 02, So 12 and So 22, and thei fist and second adial deivatives. The ψ 2,ψ 2 tems ae tems popotional to eithe ψ 2 o ψ 2. In ode that ψ 1 does not couple to ψ 2, we set the coefficients of ψ 2 and ψ 2 equal to zeo, which poduces two second ode diffeential equations fo e and e d. The coefficient of ψ 2 gives 2M ) 2 2 e 2M 2M ) 3 e 6M2M ) 4 e 48M λ)M 2 41 λ) 3 2M λ 2iω) 2 2)) 2M ) 5 e d 2 12M λ)M 2 2 2λ iω) 2 2)) 4 e 2C 2M) d =, 2.22) and the coefficient of ψ 2 yields 2M ) 2 2 e d 2M2M ) 3 e d 2 2M )2 2 e 2M4M ) 4 e d = )

29 Equations 2.22) and 2.23) can be solved by substituting seies tial solutions, 19 namely, Doing so leads to e ty = n=4 n= 3 a n n, ety d = n=4 n= 3 b n n. 2.24) and e = C 2CM iω) 2 3C d M 2 3CM2 iω) 2 2 e d = C iω) 2 C d 2C dm 4CM2 iω) 2 The seies teminate if iω)2 2C d 1 λ) 3C M) CM 2λM) 2iω) 4 3 O 4 ) 2.25) iω)2 2C d 1 λ) 3C M) CM 2λM) 2iω) 4 2 M 4C4 3λ)M iω) 2 C d 11 8λ) 6C M) ) 2iω) 4 3 O 4 ). 2.26) C = C3 2λ) 3iω) 2, C d = 2CM iω) 2, 2.27) which gives e = C 6M 3 2λ)) 3iω) 2, e d = C 2M ) iω) ) These solutions can be veified by substitution into 2.22) and 2.23). If we choose C = iω) 2, we have h ty 2M ) 2 = iω) 2 ψ 2 ψ 1 6M 3 2λ)) iω) 2 3iω) 2 ψ 2 16π f 02 So 02 f 12 So 12 f 22 So 22 ). 2.29) As mentioned peviously, the souce tems in equation 2.21) include second adial deivatives of So 02, So 02 and So 02. The second deivatives ae eliminated if f 02 = f 12 = 0 and f 22 = 2M iω) )

30 Substituting these esults into h ty 2, we finally obtain the solution to equation 2.19), ] )So 22 h 2 = 1 [ iω) 2 2M)ψ 2 ψ 1 and its adial deivative, 6M 3 2λ) ψ 2 16π 1 2M 3, ) h 2 = 1 [ 2 6M 3 λ)) iω) 2 ψ 2 ψ 1 3 8M2 42 λ)m 2 2 2λ iω) 2 2 ) 2M ) 2 ψ 2 16π 1 2M ) So 12 32π 1 2M )So 22 ]. 2.32) Equation 2.15) defines ψ 2 in tems of odd paity adial functions. An analogous expession fo ψ 1 is deived by solving 2.31) fo ψ 1 and using 2.18) and 2.17) to eliminate ψ 2 and ψ 2, espectively. The esult is ψ 1 = iωh 0 1 2M ) 2λ h 1 2 ) 3 h 2 h ) Substituting 2.33) into the spin 1 genealized Regge-Wheele equation gives the diffeential equation fo ψ 1, L 1 ψ 1 = 323 λ)π 2M)2 32λπ 2M) 3 3 So So 22 16π 2M)3 3 So 32λπ 2M) So ) Equation 2.8) may be used to eliminate So 12. Although 2.15) is gauge invaiant, equation 2.33) is not: it is valid only in the hamonic gauge. If a gauge change is made which satisfies 2.12) and theeby peseves the hamonic gauge, ψ 1 changes by a homogeneous spin 1 solution in othe wods, a homogeneous solution of 2.34)). Stated diffeently, ψ1 new = ψ1 old ψ1 hom = iωh new 0 1 2M ) 2λ h new hnew 2 dhnew 2 d ), 2.35) whee ψ hom 1 efes to the homogeneous solution. Even though ψ 1 has changed by ψ hom 1, the ight side of 2.35) has the same fom as 2.33). In this limited sense, the fomula

31 fo ψ 1 in 2.33) is invaiant unde gauge changes which peseve the hamonic gauge. Nevetheless, the behavio of ψ 1 is diffeent fom that of ψ 2, because ψ new 2 = ψ old 2 afte any change of gauge. The above solutions fo h 2 and h 2 ae substituted into 2.18) and 2.17) to obtain h 0 = 1 ψ 1 2λ ) iω 3 ψ ) and h 1 = 1 [ iω) 2 2λ 3 ψ 2 2 ψ 1 2λ 3 ψ 2 16π 1 2M ) ] So 12 ψ 1, 2.37) and thei adial deivatives, and h 0 = 1 ψ 1 2λ ) iω 3 ψ ) h 1 = 1 [ 2λ 4M ) iω) 2 32M ) ψ 2 8M2 43 λ)m 2 4 2λ iω 2 ) 2 ) 2 2M ) 2 ψ 1 2λ8M2 23 2λ)M 2 1 2λ iω) 2 2 ) 3 2 2M ) 2 ψ 2 32πM ) So 12 32πλ 2 M ) 2 So 22 2M ) ψ 1 ]. 2.39) One may veify by substitution that these solutions satisfy the field equations, as well as the hamonic gauge condition, the fist ode equation 2.13), and the definitions of ψ 2 and ψ 1. The solutions also can be checked by tansfoming fom the hamonic gauge H ) to the Regge-Wheele gauge RW ). By definition, h RW 2 = 0. Applying equation 2.11), we set Z = h H 2 h RW 0 = and substitute into equations 2.9) and 2.10) to obtain 2M ) ψ 2 iω 1 2M ) ψ2 iω 1 2M ) 16π iω So 22, 2.40) and h RW 1 = 2 2M ψ )

32 22 Adjusting fo diffeences in notation, the Regge-Wheele gauge solutions agee with those obtained by Zeilli and othes [4], [102], [115]. Equation 2.41) can be solved fo ψ 2 to give the Regge-Wheele expession fo ψ 2. In the limit h 2 0, Moncief s fomula, equation 2.15), educes to theis. The peceding paagaph suggests anothe way of deiving the hamonic gauge solutions. Instead of woking thoughout in the hamonic gauge, we could have stated in the Regge-Wheele gauge whee the solutions ae known) and looked fo the gauge tansfomation vecto that would take us fom the Regge-Wheele gauge to the hamonic gauge. To do so, we would need to deive the adial function Z in 1.25) that we would then substitute into 2.9)-2.11), with the supescipt old efeing to the Regge- Wheele gauge and new being the hamonic gauge. This appoach was begun in [102], whee a diffeential equation fo the gauge tansfomation vecto was deived. The equation is simila to 2.19) above, although the deivation in [102] is vey diffeent. Howeve, the authos of [102] did not solve thei gauge tansfomation equation they put it aside fo futue study ) and did not complete the deivation of the hamonic gauge solutions. Because equation 2.13) is gauge invaiant, the fomula fo ψ 2 in 2.15) is not unique. We can solve 2.13) fo h 2 and substitute fo h 2 in 2.15). The esult is ψ 2 = iωψ JT 2 8π 2M) So 12, 2.42) λ whee ψ2 JT = 1 h 0 iω λ 2 h 1 ) 2 h ) The supescipt JT stands fo Jhingan and Tanaka, who deived ψ JT 2 in the Regge- Wheele gauge using a diffeent method and notation and showed that it is gauge invaiant [61]. We will use the supescipt JT because they discuss this fom extensively, even though they acknowledge it was deived ealie by othes. In 2.43), the facto of iω indicates that ψ 2 is the time deivative of ψ2 JT, and Jhingan and Tanaka

33 constucted ψ JT 2 so that it would be the Fouie tansfom of a time integal of ψ 2. Applying equations 2.9) and 2.10), a gauge change in h 1 is canceled by coesponding gauge changes to h 0 and h 0, leaving ψjt 2 invaiant. Afte substituting ψ2 JT into the Regge-Wheele equation and simplifying, we obtain a souce tem given by L 2 ψ2 JT 8π 2M) = iωso12 So ) ) λ Adjusting fo diffeences in notation, the souce tem agees with that in [61]. Using 2.42), the field equation solutions can be ewitten in tems of ψ2 JT, if desied. This completes the deivation of the odd paity non-zeo fequency solutions, fo l 2. We stated with the hamonic gauge field equations given by 1.30), which ae patial diffeential equations. Using sepaation of vaiables, that system was educed to the thee adial field equations in 2.1)-2.3), a system of coupled odinay diffeential equations. The solutions to the adial equations ae witten in tems of ψ 2 and ψ 1, each of which satisfies its own decoupled second ode diffeential equation. To calculate h 0, h 1 and h 2, we would fist solve the decoupled equations fo ψ 2 and ψ 1 and then substitute the esults into the solutions 2.36), 2.37) and 2.31). In this manne, we have simplified the poblem fom a system of coupled patial diffeential equations to two decoupled odinay diffeential equations Solutions fo l = 1 Although subsection assumed that l 2, most of the esults deived thee can be used fo l = 1. One diffeence is that h 2 is no longe pesent. In equation 1.12) fo the odd paity metic petubations, the angula functions W and X ae zeo [115], so h 2 does not exist fo l = 1. Similaly, So 22 is non-existent. Futhe, thee ae only two field equations, 2.1) and 2.2), because 2.3) is the adial coefficient of W o X. To show that W and X ae zeo, substitute the spheical hamonics fo l = 1 [3] into the definitions of W 1.13) and X 1.14).

34 24 If l = 1, then λ = 0, fom the definition of λ 1.51). In the emaining field equations, any tems involving h 2 have a facto of λ, so such tems ae zeo. Othe, elated equations 2.7)-2.10) and 2.12)-2.13) also still apply, with the substitution λ 0 to ensue that tems with h 2 and So 22 ae zeo. Finally, the solutions fo h 0 and h 1 apply, again with λ 0. This means the solutions do not depend on ψ 2, so ψ 2 is not defined fo l = 1. The solutions now depend on ψ 1, as given by 2.33) and 2.34). 2.2 Zeo Fequency Solutions The zeo-fequency equations and solutions have ω = 0. Factos of ω ae due to time deivatives of the facto e iωt. Because this exponential contains the time dependence of the metic petubation, solutions with ω = 0 ae time independent solutions. Below, the cases l 2 and l = 1 ae discussed sepaately Solutions fo l 2 We can use the field and elated equations fom subsection 2.1.1, afte substituting ω = 0. Howeve, h 0 no longe couples to h 1 and h 2 though the adial field equations. Instead, h 0 has a sepaate field equation, 2.1), while h 1 and h 2 solve a coupled system, 2.2) and 2.3). This decoupling of h 0 also appeas in othe equations, namely, the hamonic gauge condition 2.7) and the gauge invaiant fist ode equation 2.13). Neithe contains h 0 if ω = 0. Futhe, the souce fo 2.1), So 02, no longe couples to So 12 and So 22 though the stess enegy tenso divegence equation 2.8). We will solve fo h 1 and h 2 fist. Using 2.13) to eliminate h 1 fom 2.3), we have 1 2M ) 2 h 2 1 2M ) 2M 21 λ)2m ) 2 h 2 3 h 2 = 16π 1 2M ) [ ] 2M) So 12 So ) λ

35 25 This is the genealized Regge-Wheele equation, with s = 1. Accodingly, h 2 = ψ 1, h 2 = ψ 1, 2.46) whee L 1 ψ 1 = 16π 1 2M ) [ ] 2M) So 12 So 22 λ. 2.47) Hee, ψ 1 is diffeent fom ψ 1 in subsection 2.1.1, even though the same notation is used fo both functions. Substituting 2.46) into 2.13) gives and the adial deivative is h 1 h 1 = 2 ψ 1 ψ 1 8π2 λ So 12, 2.48) 2 λ)) 2 M) 16π M) = 2 2M 2 ψ 1 2M) 2M) ψ 1 λ2m ) So ) By substitution, the eade may veify that 2.46) and 2.48) ae solutions to the zeo fequency field and elated equations that involve h 1 and h 2. The two solutions depend on ψ 1, but not ψ 2. This means the definition of ψ 2 in 2.15) is not valid fo zeo fequency, as that definition is in tems of h 1 and h 2. The next step is to solve the emaining field equation, 2.1), fo h 0. Fom 2.9), h 0 is gauge invaiant when ω = 0, and the solution will be in tems of the spin 2 Regge- Wheele function, which is also gauge invaiant. Using 2.43), we set ψ 2 = ψ JT 2, so that and ψ 2 = 1 λ L 2 ψ 2 = h 0 ) 2 h 0, 2.50) 8π 2M) So 02 λ. 2.51) Afte diffeentiating 2.50) and using 2.1) to eliminate h 0, we find that ψ 2 λ) 8π 2 = 2M λ2m 2 h 0 ) λ 2M ) So λ h )

36 Equations 2.50) and 2.52) can be solved fo h 0 and h 0, giving h 0 = 2M)ψ 2 1 2M 26 )ψ 2 8π2 λ So 02, 2.53) and h 0 = 2 1 2M ) ψ 2 2 2M λ) 2 ψ 2 16π λ So ) Fo zeo fequency and l 2, the metic petubation can be witten in tems of spin 2 and spin 1 genealized Regge-Wheele functions, just as fo non-zeo fequency modes Solutions fo l = 1 Fo l = 1, the function h 2 is not pesent, as explained in subsection Fom 2.1) and 2.2), the field equations become 1 2M ) 2 h 0 2 2M)2 4 h 0 = 16π 1 2M ) So 02, 2.55) 1 2M ) 2 h 1 4M 2 1 2M ) h 1 2 8M 2 10M 3 2) 4 h 1 = 16π 1 2M ) So ) The hamonic gauge condition simplifies to 1 2M ) h 2M ) 1 2 h 1 = ) Fo the non-zeo fequency mode having l = 1, we used the solutions fo l 2, but we can not do that hee. The zeo fequency solutions fo l 2 have factos of 1 λ, but λ = 0 if l = 1. The ules fo a change of gauge ae somewhat diffeent fom 2.9)-2.10) and ae discussed in [108]. Rewitten in tems of ξ1 o t,) instead of Zω,), the ules become h new 0 = h old 0 ξo 1 t,), 2.58) t h new 1 = h old 1 2 ξo 1 t,) ξo 1 t,), 2.59)

37 27 whee the eplacement iω t has been used in 2.9). The symbol ξo 1 is shot fo ξl=1 odd. Refeing to 1.24)-1.25), the odd paity gauge change vecto fo l = 1 is 0,0,ξ1 ot,)csc θ Y 1m φ, ξo 1 t,)sin θ Y 1m θ The ules ae modified in ode to allow a gauge change of the fom ). 2.60) ξ o 1t,) = C 1 t 2, 2.61) whee C 1 is a constant. Although witten in tems of t, this is actually a change in the coodinate φ [108]. The modification is necessay in ode to allow a gauge change in h 0. If we substitute 2.61) into 2.58) and 2.59), then h 0 changes by C 1 2, but h 1 is unalteed [108]. Because ξ1 ot,) is only linea in t, h 0 emains time independent. This fom of gauge change could not have been used fo othe modes. Fo non-zeo fequency modes, the gauge vecto time dependence is in the facto e iωt. Fo zeo fequency modes with l 2, the new fom would cause h 2 to gow linealy with time, as 2.11) shows. As a esult, the zeo fequency h 0 is effectively gauge invaiant fo l 2 by 2.9)), but gauge dependent fo l = 1. Fom 2.12), a gauge change which peseves the hamonic gauge must be a solution of 2 ξ o 1 t,) 2 2 ξ o 1 t,) t 2 1 2M ) 2 2ξo 1t,) = 0, 2.62) which is the genealized Regge-Wheele equation with s = 1, witten in the time domain. The most geneal zeo fequency gauge change which satisfies 2.62) and leaves the petubation time independent is ξ1t,) o = C 1 t 2 2M 2 2M 2 ln [ 1 2M C 2 8M 3 C ) If ξ1 o t,) wee non-linea in t, o had some othe time dependence, hnew 0 and h new 1 would be time dependent. The fist tem of 2.63) changes only h 0 ; the second, only h 1. The last tem, C 3 2, changes neithe h 0 no h 1, so we will disegad it. ]

38 to give Fo l = 1, the zeo and non-zeo fequency gauge changes ules may be combined ξ1t,) o = C 1 t 2 e iωt Zω,)dω. 2.64) This satisfies 2.62), povided that Zω,) is a solution of 2.12) with λ = 0. Fo nonzeo fequency, Zω,) is the function Z efeed to in section 2.1. Fo zeo fequency, Zω = 0,) is equal to the last two tems of 2.63). We will solve the field equations fist fo h 1 and then fo h 0. If we use 2.57) and its adial deivative to eliminate h 1 and then h 1 fom 2.56), we find that the left side of 2.56) is educed to zeo. This means So 12 = 0. Moe diectly, the same esult can be obtained fom 2.13), by substituting λ = 0 and ω = 0. Because So 12 = 0, the function h 1 is a homogeneous solution of 2.56) and 2.57). Accodingly, h 1 = 28 C 2M), 2.65) whee C is a constant. If C is non-zeo, we can nevetheless zeo out h 1 by a gauge tansfomation descibed in 2.63), with C 1 = C 3 = 0 and C 2 = C. Because h 1 is entiely gauge dependent, we can set C = 0 in 2.65), so h 1 = 0 fo l = 1 and ω = 0. The othe adial function, h 0, epesents the obital angula momentum of the small mass [115]. Equation 2.55) has homogeneous solutions of 2M and 2M )2. Fom them, we can constuct an inhomogeneous solution using vaiation of paametes [67], o a Geen s function [60]. Eithe method poduces h 0 = 1 16π 3 So 02 ) 2M 3 d 2 16πSo 02 ) 2M) 3 2M) d, 2.66) h 0 = 1 16π 3 So 02 ) 2 3 d 16πSo 02 ) 2 2M) 3 2M) d, 2.67) 2M whee So 02 ) is shot fo So 02 ω = 0, ). Altenatively, the solution fo h 0 may be witten in tems of a spin 1 Regge-Wheele function and its adial deivative; howeve, the esulting expessions ae moe complicated than those given above. The limits 2M and ae geneic and should be eplaced by limits appopiate to the obital motion.

39 Fo example, an elliptic obit would have limits of min peiaston) instead of 2M and max apaston) instead of. We can analytically evaluate the integals in h ) fo a cicula obit of adius R, with the stess enegy tenso expessions in Chapte 5. Doing so gives h 0 = m 0 4 L 2 R 3 29 π 3 θr ) m 4 L π 0 θ R), 2.68) 3 whee we have used [3] Y 10 θ,φ) θ 3 = 4π sin θ. 2.69) The facto L is the obital angula momentum pe unit mass 4.10), 4.26). To find the petubation, we substitute h ) and h 1 = 0 into h o,lm µν 1.12), getting ) ) h 10 tφ = h10 φt = 2 L 2 2 L m 0 R 3 sin2 θ θr ) m 0 sin2 θ θ R). 2.70) The othe petubation components ae zeo fo this mode, and the supescipt 10 is shot fo l = 1,m = 0. The solution 2.70) was also calculated by Detweile and Poisson [34], following Zeilli [115]. They note the solution fo > R esembles the lineaized Ke metic. Specifically, the 1 tem has the same fom as the g tφ component of the Ke metic lineaized to Oa), with m 0 L/M taking the place of the Ke angula momentum paamete a. In contast, the solution fo < R diffes fom the backgound metic only by a gauge tansfomation [34]. As discussed in [115] and [34], the solution fo this multipole epesents the obital angula momentum of the mass m 0 and the change in angula momentum that occus at the obital adius. 2.3 Homogeneous Solutions This section discusses homogeneous solutions of the odd paity field equations. We will stat with the non-zeo fequency case, with l 2. The adial field equations in 2.1)-2.3) fom a system of thee second ode linea odinay diffeential equations. Fom the theoy of diffeential equations, such a

40 30 system may be efomulated as an equivalent system of six fist ode equations, and the new system will have at most six linealy independent homogeneous solutions [117]. Accodingly, one might expect that the field equations could also have six homogeneous solutions. Howeve, conside the following system of fou fist ode homogeneous linea diffeential equations: h 2 d 2 = 0, 2.71) 1 2M ) 2 d 2 22M )3M ) 2 2M)2 3 d 2 3 h 1 1 2M 16M2 4 3 λ)m 2 2 2λ iω) 2 2 ) 4 h 2 = 0, 2.72) ) h 1 iω 2M h 0 2M ) 2 h 1 2λ 2 h 2 = 0, 2.73) iωh 0 2λ 2 1 2M ) d 2 2iω h 0 iω) 2 λ4m 2) 4λ 3 3 ) h 1 1 2M ) h 2 = ) Equation 2.72) is the homogeneous fom of 2.3), ewitten in tems of d ). Similaly, 2.73) is the hamonic gauge condition, 2.7), and 2.74) is fom 2.13). Thus, the system of fou is deived fom the field equations and the hamonic gauge condition. Moeove, it is possible to show that the homogeneous fom of the field equations can be deived fom 2.71)-2.74). Each system implies the othe, so the system of fou and the homogeneous field equations ae equivalent. The system of fou has only fou linealy independent solutions. Because these fou equations ae equivalent to the homogeneous fom of the field equations, the odd paity field equations also have only fou linealy independent homogeneous solutions. This eduction fom six to fou evidently occus because of the hamonic gauge condition, which is a diffeentiable constaint. Solutions to a system of linea diffeential equations ae witten in the fom of column vectos. If a matix fomed fom the solution vectos has a non-zeo deteminant, the solutions ae linealy independent [117].

41 We need to find fou linealy independent solution vectos fo the system 2.71)- 2.74). As explained in Chapte 6, the non-zeo fequency genealized Regge-Wheele equation has two linealy independent homogeneous solutions: ψ in s, which epesents an ingoing wave nea the event hoizon, and ψs out, which epesents an outgoing wave at lage. This suggests that we fom one solution vecto fom each of ψ1 in, ψout 1, ψ2 in and 31 ψ out 2. The fist solution vecto is [ ] ψ X1 T in = 1 iω, 1 2 iω) 2 ψin 1 ψin 1 ), ψ1 in iω) 2, ψin 1 ) ) iω) ) The T is fo tanspose, because X 1 is a column vecto, athe than the ow vecto displayed above. The fist component of 2.75) comes fom the ψ 1 tem of h ), the second component fom the ψ 1 tems of h ), the thid fom h ) and the fouth fom d 2 = h ). The second solution, X 2, is obtained fom X 1 by eplacing ψ in 1 with ψout 1. The thid vecto is 2λ X3 T = 3iω ψin 2, 2λ [ 3iω) 2 ψ in 2 ) ψin 2 ], 1 iω) 2 ], 6M 3 2λ) ψ2 in 3 [ 2M)ψ in 2 ) 1 iω) 2 [ Aψ in 2 ) B ψ2 in ] ), 2.76) whee A = 2 6M 3 λ)) 3, 2.77) B = 8M2 42 λ)m 2 2 2λ iω) 2 2 ) 2M ) ) The fouth solution, X 4, is obtained fom X 3 by eplacing ψ in 2 with ψout 2. The function 4 W 1 W 2 iω) 5 2M) ) is the deteminant of the matix whose fou columns ae X 1 though X 4. In equation 2.79), the factos W 1 and W 2 ae the Wonskians W s 6.6) of the genealized Regge-Wheele homogeneous solutions. Because the deteminant 2.79) is non-zeo, the solution vectos X 1 though X 4 ae linealy independent. Fom the theoy of odinay

42 32 diffeential equations, any homogeneous solution to the system 2.71)-2.74) and, by extension, the field equations 2.1)-2.3)) can be witten as a linea combination of the solution vectos descibed above. Accodingly, homogeneous solutions of the odd paity field equations ae fomed fom combinations of the genealized Regge-Wheele functions ψ 2 and ψ 1. This esult is to be expected, given the fom of the inhomogeneous solutions deived in subsection Gauge tansfomations which peseve the hamonic gauge ae implemented by adding homogeneous spin 1 solutions 2.12). Because of this gauge feedom, the ψ 1 contibutions to the homogeneous solutions above can be emoved by a gauge tansfomation which peseves the hamonic gauge. Howeve, the ψ 2 solutions can not be so emoved. This is because the spin 2 Regge-Wheele functions ae gauge invaiant. To modify the two solution vectos attibutable to ψ 2 X 3, X 4 ), we would have to tansfom to a diffeent gauge, such as the Regge-Wheele gauge. Also, a hamonic gauge peseving change adds only homogeneous solutions, so it can not emove the inhomogeneous ψ 1 contibutions deived in subsection The homogeneous solutions above assumed l 2. Fo l = 1, the adial functions h 2 and h 2 ae not pesent, as explained in subsection We ae left with a two equation system that is composed of 2.73)-2.74), modified by the substitution λ = 0. Using simila aguments to the l 2 case, we can show that the l = 1 homogeneous solutions can be witten solely in tems of ψ in 1 and ψout 1. Specifically, the solution vectos ae the fist two components of X 1 and X 2. Moeove, these homogeneous solutions can be emoved by a gauge tansfomation which peseves the hamonic gauge. The homogeneous solutions above ae fo non-zeo fequency. Fo zeo fequency and l 2, we use the following fist ode system instead: 1 2M h 0 d 0 = 0, 2.80) ) 2 d 0 8M2 42 λ)m 2 2 2λ)) 4 h 0 = 0, 2.81)

43 2λ 2 1 2M 1 2M ) h 2 ) h 1 2M ) 2 h 1 2λ λ4m 2) 3 h 1 4λ h 2 = 0, 2.82) 1 2M ) h 2 = ) We use these fou fist ode equations because of the decoupling of the zeo fequency h 0 fom h 1 and h 2. If we had wished, we could have used the non-zeo fequency equivalents of 2.80)-2.83) instead of 2.71)-2.74) in ou discussion of the non-zeo fequency case. Equation 2.81) is the homogeneous fom of the field equation 2.1), ewitten in tems of d ). Equation 2.82) is the zeo fequency hamonic gauge condition 2.7), and equation 2.83) is the zeo fequency fom of 2.13). It is possible to show that the system 2.80)-2.83) is equivalent to the zeo fequency field equations, in the same way that 2.71)-2.74) ae equivalent to the non-zeo fequency field equations. Accodingly, the zeo fequency field equations have only fou linealy independent homogeneous solutions. As discussed in Chapte 6, thee ae two linealy independent zeo fequency homogeneous solutions of the genealized Regge-Wheele equation. One solution, ψ in s, is finite as 2M, but diveges like l1 as. The othe solution, ψs out, is bounded fo lage, but diveges logaithmically nea the event hoizon. With aguments simila to the non-zeo fequency case, we can show that all zeo fequency homogeneous solutions can be witten as combinations of ψ1 in, ψout 1, ψ2 in and ψout 2. Futhe, the ψ in 1 and ψ out 1 solutions can be emoved by means of a gauge tansfomation which peseves the hamonic gauge. Although the spin 2 solutions emain, they ae divegent, eithe at the hoizon o fo lage. To pevent an unphysical divegence, we set them equal to zeo by choice of integation constants. A simila conclusion was eached by Vishveshwaa [108], although he did not wok in the hamonic gauge. Fo l = 1, h 2 is no longe pesent, and the homogeneous solution fo h 1 is zeo, as discussed in subsection This leaves the single second ode diffeential equation fo h 0, which is witten as two fist ode equations in 2.80)-2.81). Fo efeence, the