arxiv: v1 [gr-qc] 24 Mar 2014 The Self-Force Problem: Local Behaviour of the Detweiler-Whiting Singular Field

Transcription

1 arxiv: v1 [gr-qc] 4 Mar 14 The Self-Force Problem: Local Behaviour of the Detweiler-Whiting Singular Field Anna Heffernan The thesis is submitted to University College Dublin in fulfillment of the requirements for the degree of PhilosophiæDoctor PhD) School of Mathematical Sciences and Complex and Adaptive Systems Laboratory Head of School: Dr. Patrick Murphy Principal Supervisor: Prof. Adrian Ottewill 1 October

2 Contents List of Figures vi List of Tables viii Abstract ix 1 Introduction Einstein s Theory of General Relativity Testing the Theory Gravitational Wave Astronomy The Two-Body Problem Black Hole Binary Sources Modelling Techniques The Self Force Problem - Thesis Outline History of Theory Main Approaches Thesis Outline Background 18.1 Bitensors and Basics Geodesics Penrose Diagrams Synge s World Function Bivector of Parallel Transport ii

3 CONTENTS.1.5 The Van Vleck Determinant Black Hole Space-Times Spherically Symmetric Space-times Schwarzschild Space-Time Reissner-Nordström Space-Time A General Static, Spherically Symmetric Space-Time 48.. Axially Symmetric, Stationary Space-Time Kerr Space-Time Kerr-Newman Space-Time The Detweiler-Whiting Singular Field The wave equation Singular Field in Flat Space-time Singular Green s Function in Curved Space-time Restrictions on U A B x, x ) and V A B x, x ) The Singular Field The Self-Force and the Singular Field High-Order Expansions of the Singular Field Coordinate and Covariant Expansions of Fundamental Bitensors Coordinate Approach Covariant Approach * Expansions of the Singular Field* Scalar singular field Electromagnetic singular field Gravitational singular field Mode-Sum Decomposition of the Singular Field Mode Sum Concept Rotated Coordinates Mode decomposition Schwarzschild Space-time Scalar case Electromagnetic case iii

4 CONTENTS Gravitational case Self-force regularization huu regularization Example - Scalar Self-Force Impact of Regularisation Parameters in Schwarzschild Space- Time Kerr Space-time Scalar Case Electromagnetic Case Gravitational Case Self-Force Regularisation huu Regularisation Results Regularisation Parameters to Date Effective Source Methods World Tube Method Window Function Method Effective Source from the Singular Field m-mode Scheme Derivation of m-mode regularization parameters m-mode regularization parameters Example - Kerr Scalar Self-Force Further Extensions Non-geodesic Motion Cosmic Censorship Conjecture Singular Field Mode Sum Decomposition fr) Regularisation Parameters Reissner-Nordström Regularisation Parameters Second Order Self-force iv

5 CONTENTS 6..1 Mode-sums for the Second Derivative of the Singular Field Regularisation Parameters Example Discussion 1 A Coordinate Expansions in Schwarzschild Space-Time 7 A.1 Synge World function A. Van Vleck determinant A.3 Expansions of an arbitrary point on the world line about xā A.4 Expansions of retarded and advanced points A.5 Bivector of Parallel Transport A.6 Scalar Singular Field A.7 Electromagnetic Singular Field A.8 Gravitational Singular Field B Covariant Bitensor Expansions* 19 B.1 Advanced and retarded points B. Advanced and retarded distances B.3 Van Vleck Determinant B.4 Bivector of parallel transport B.5 Scalar tail B.6 Electromagnetic tail B.7 Gravitational tail References 9 v

6 List of Figures 1.1 Gravitational Wave Stretching Space Gravitational Wave Detectors Sensitivity Curves Orbits Around Supermassive Black Hole Sagittarious A* Black Hole Parameter Space Light-cone of Events Minkowski Space-Time Inverse Tangent Penrose Diagram for Minkowski Space-Time Derivaties of Synge s World Function Kruskal Diagram for Schwarzschild Space-Time Kruskal Diagram for Infalling Particle in Schwarzschild Space-Time Penrose Diagram for Schwarzschild Space-Time Motion of Charged and Uncharged Particles in Reissner-Nordström Space-Time Penrose Diagram for Reissner-Nordström Space-Time Q < M Penrose Diagram for Reissner-Nordström Space-Time Q = M Penrose Diagram for Reissner-Nordström Space-Time Q > M Retarded Green s Function in Flat Space-time Advanced Green s Function in Flat Space-time Retarded Green s Function in Curved Space-time Advanced Green s Function in Curved Space-time Singular Green s Function in Curved Space-time Regular Green s Function in Curved Space-time vi

7 LIST OF FIGURES 3.1 Worldline Singular Field Components in D Singular Field Components in 3D Regularization of the Scalar Self-Force Radial Component in Schwarzschild Space-Time Regularization of the Electromagnetic Self-Force Radial Component in Schwarzschild Space-Time Regularization of the Gravitational Self-Force Radial Component in Schwarzschild Space-Time Regularisation of the Scalar Self-Force Radial Component in Kerr Space-Time Effective Source Components in D Effecitive Source Components in 3D m-mode Regularization of the Radial Component of the Self-Force Regularisation of the Second Radial Derivative of the Retarded Field vii

8 List of Tables 4.1 Regularization Parameters Notation Impact of Regularisation Parameters for Scalar Self-Force Regularization Parameters to Date viii

9 Abstract Gravitational waves are ripples in space-time and a prediction of Einsteins theory of relativity. The growing reality of gravitational wave astronomy is giving age-old problems a new lease of life; one such problem is that of the self-force. A charged or massive particle moving in a curved background space-time gives rise to a field that affects its motion, pushing it off its expected geodesic. This self-field gives rise to a so-called self-force acting on the particle. In modelling this motion, the self-force approach uses a perturbative expansion in the mass ratio. One of the most interesting sources of gravitational waves are extreme mass ratio inspirals. These systems have an extremely small mass ratio, making them perfectly suited to perturbative, gravitational self-force modelling. One of the key problems that immediately arises, within the self-force model, is the divergence of the field at the particle. To resolve this, the field is split into a singular component and a smooth regular field. This regular-singular split, introduced by Detweiler and Whiting, is used in most modern self-force calculations. In this thesis, we derive high order expansions of the Detweiler-Whiting singular field, and use these to push the boundaries on current precision limits of self-force calculations. Within the mode sum scheme, we give over 14 previously unknown regularisation parameters, almost doubling the current regularisation parameter database. We also produce smooth effective sources to high order, and propose an application of the higher terms to improve accuracy in the m-mode scheme. Finally, we investigate the status of the cosmic censorship conjecture and the role that the self-force plays. To this end, we give regularisation parameters for non-geodesic motion. Additionally, we show the necessity of our results in the exciting area of second order self-force calculations. Recently, second order self-force derivations have been developed, which benefit significantly from high-order coordinate expansions of the singular field, making them an immediate application of our current work. We calculate several parameters that these schemes require, and highlight the further advancements possible from the results of this thesis.

10 Statement of Original Authorship I hereby certify that the submitted work is my own work, was completed while registered as a candidate for the degree stated on the Title Page, and I have not obtained a degree elsewhere on the basis of the research presented in this submitted work.

11 Acknowledgements First and foremost, I would like to offer a massive thank you to my supervisor, Adrian Ottewill, the man of never-ending patience, for his constant support, supervision, good humour, and above all, enthusiasm in sharing both his brilliance and intuition in both the complex and mundane. There really was no question too big or small) and I count myself very lucky to have had such a supervisor. Throughout my PhD, I have had the good fortune of working in a fun and helpful group based in the Complex and Adaptive Systems Laboratory CASL) at the University College of Dublin UCD). In particular, my colleague, collaborator and good friend, Barry Wardell, has been the absolute best in all his roles. Switching the topic of my PhD would not have been possible without his constant support and friendship, be it in front of a desktop looking at Mathematica or in more social settings, and I am truly grateful for this. Marc Casals and Sam Dolan also deserve a special mention, for both their enthusiasm and support. Brien Nolan, Cliona Golden, David Brodigan, Patrick Nolan, Chris Kavanagh and all of my friends in CASL have also helped me on this journey. Ted Cox deserves a special mention, for his constant support throughout my education and research in Mathematical Physics, from suggesting my initial degree course in Theoretical Physics to initiating the process that became my PhD. There are few lecturers who share his enthusiasm for teaching and research and I know how lucky both myself and UCD are to have him in our lives. Nuria Garcia, Veronica Barker, and all of the support staff in both the School of Mathematical Sciences and CASL have been exemplary in their roles. I know I didn t make their lives easy, but regardless, they were tireless in both their support and good

12 humour throughout my PhD. A quick mention to my mathematics teachers through the years - Ms. Creaney, Mr Connelly and Mr. Gunning, your teachings and enthusiasm were never wasted, thank you. Through my research, I have had the opportunity and pleasure of meeting many helpful and fun researchers who have given me their time, in particular, Leor Barack and Niels Warburton have been especially helpful. Roland Haas, Sarp Akcay, Jonathan Thornburg and many others from the Capra meetings were also very generous with their data, assistance and interesting discussions. Finally I would like to thank all my friends and family, for their support as well as keeping me sane through the years. To all the girls, thank you. My sister, brothers, Mum and Dad have all been amazing throughout the years, from ensuring I don t dress like a boy, to driving me around, dragging me away from my mathematics when I needed it and supporting me both financially and emotionally, I got there in the end and it was thanks to you. This research was financially supported by the Irish Research Council for Science, Engineering and Technology funded by the National Development Plan) as well as the School of Mathematical Sciences, UCD.

13 To my Dad

14 Chapter 1 Introduction Every once in a while, science makes a ground-breaking discovery. This year we were lucky enough to witness such an event - it will be remembered as the year in which the Higgs Boson was finally detected. After decades of searching and non-stop research by both theorists and experimentalists of high energy particle physics, a detection was accomplished at CERN earlier this year. As is the nature with many scientific breakthroughs, the excitement of the Higgs Boson came in two waves. First is the theoretical wave, in this case, the production of the theory of electroweak unification 1,, 3) and with it, the prediction of the Higgs Boson 4). As with most revolutionary theories, it took several years for people to warm to the initial idea, but after much investigation, the theory spoke for itself and became recognised as part of the standard model, to be taught to physics students globally. As with all exciting theories, there then comes the search for physical evidence - a search which, in this case, would last almost half a decade, and result in the second wave of excitation - physical clarification that the theory is correct in the form of a direct detection of the Higgs Boson. It was a momentous occasion for every researcher who has given their time and patience to the area. And while this was all happening, those of us sitting in the gravitational research area, also thrilled by the result, couldn t help but think - it s our turn next. 1

15 1.1 Einstein s Theory of General Relativity 1.1 Einstein s Theory of General Relativity Testing the Theory Einstein s theory of general relativity GR) was a revolutionary step in fundamental physics 5). Like many of his era, Einstein was unsatisfied by the then accepted model of Newtonian physics, due to its inability to explain several observed effects in the world or universe around us and its unsatisfactory concept of absolute time and space. GR successfully united Newtonian Mechanics and Special Relativity and had an immediate success as it naturally explained the precession of the perihelion of Mercury - an observation for which Newtonian theory could not completely account. Depsite this initial success, there were many sceptics to the notion of curving space and time. However, since there were other predictions by GR that would differ from Newtonian mechanics, it would remain only a matter of time before the theory was fully accepted. A massive step in this direction was taken in 1919 by Sir Arthur Eddington. Having been one of the first to receive news of the theory of GR, he organised two expeditions to observe a solar eclipse. The reason was to measure the deflection of light by the sun, as Einstein s theory would predict a different value for this observation than that of Newtonian mechanics. The experiment was a success 6) and Einstein become world famous almost over night, while his theory started to overthrow its Newtonian counterpart. Since 1919, there have been many more experiments testing the various available observables that can be used to support GR. These have included verifying the gravitational redshift of light 7), gravitational lensing 8) and time delay 9), to name a few. One of the most exciting results to further fortify GR is the indirect detection of gravitational waves. Gravitational waves are ripples in space-time as predicted by GR; they can arise from various events - compact object binaries, black hole mergers and supernovae are just a few examples. In 1974, Hulse and Taylor discovered a new type of pulsar or radiating neutron star - one with another pulsar in its orbit 1). By observing the binary system, it was possible to calculate the orbit decay and show that the amount of energy being lost was consistent with the amount of energy that should be emitted as gravitational radiation as predicted by GR 11).

16 1.1 Einstein s Theory of General Relativity Figure 1.1: If a gravitational wave were to pass through this page at a 9 degree angle, with a plus polarisation, it would effect a circle of test masses as shown in the diagram. There also exists the cross polarization which would have the same effect but rotated 45 degrees. Regardless of polarisation, the magnitude of strain to be measured is 1 in 1 1 parts. The Nobel prize winning work of Hulse and Taylor has encouraged relativists to work on the possibility of a direct detection of gravitational waves. When a gravitational wave passes through space and time, it can be seen to stretch and squash the space it passes through, this is illustrated in Fig. 1.1, which shows a circle of test particles at rest being affected as a gravitational wave passes through this page. In order to detect the waves, it is therefore necessary to be able to measure this strain that is placed on the test particles. Due to the weakness of gravitational waves, however, this requires measuring a strain of 1 in 1 1 parts. Until the 199 s, this accuracy in measurement was believed to be impossible; however, advances in technology and research, have now made it a possibility. A direct detection of gravitational waves would mark a test of GR that would be the first of its kind - all previous tests of GR have measured the impact of GR on other observables in the weak regime while this would be a direct measurement of gravitational radiation predicted by GR in the strong field regime, i.e., when space and time are being strongly distorted. Such a detection would be analogous to the recent detection of the Higgs boson, and with it would come the same thrill of accomplishment that is currently being enjoyed by our particle physics counterparts, 3

17 1.1 Einstein s Theory of General Relativity albeit almost a century after Einstein revealed his theory Gravitational Wave Astronomy Amazing strides have been made in Astrophysics in the last 7 decades. We no longer rely solely on optical telescopes to inform us of the nature of our universe, instead there exists a network of satellites, antennas and telescopes that use optics, radio waves, infrared, X-ray and gamma rays to investigate the cosmos. With each new window, came surprises that dramatically changed our understanding of the universe, some were expected but the more exciting were the unexpected, like pulsars 1) or gamma ray bursts 13). We are now, once again, on the verge of opening a new window onto our universe - that of gravitational wave astronomy. The thrill of detecting gravitational waves is not solely in the success of the detection but also in the wealth of knowledge that we can extract from the waveforms. Gravitational waves can travel, relatively unaffected by any intervening matter, from their source to us, meaning they would carry first hand information about the violent processes that created them - processes that will often be invisible to all other types of detection available to us. This invisibility is often due to the amount of intervening matter that would affect all other types of radiation, but also, in some cases, such as those processes solely involving black holes, gravitational waves are the only type of classical radiation that will be emitted. Detection of gravitational waves is expected to occur in the next 5 years. A network of ground-based detectors LIGO 14), VIRGO 15), GEO6 16), TAMA 17)) have been operational for almost a decade - the first came online in. Although no detection has yet been made, hopes are high that the new advanced detectors will be successful. This optimism is not baseless - event rates for the gravitational wave detectors carry large error bars. It was known that the initial detectors may not be successful, whereas the advanced detectors are expecting greater event rates than their predecessors, by a factor of approximately 1. These, even with the more conservative estimates, predict that the advanced detectors should make positive detections 18). The aim was to get an array of detectors up and running and work on reducing the noise to obtain the highest signal-to-noise ratio SNR) possible. Considering these detectors are required to measure strain of one part in 4

18 1.1 Einstein s Theory of General Relativity 1 1, obtaining the optimal SNR was a learning curve - some noises, although unexpected were easily removed gunshots from hunters being such a source initially at the Louisiana LIGO site), others proved more difficult laser shot noise). In fact, during its final run, LIGO Laser Interferometer Gravitational-Wave Observatory) was able to obtain a strain sensitivity curve better than was anticipated 19). Gravitational wave detectors differ from their electromagnetic cousins in the sense that they have no ability to detect the direction from which the gravitational waves come. The detector will know when a gravitational wave passes through it, however it has no way of telling where it came from. For this reason, it has been crucial that there be a world network of detectors - by comparing what times each detector senses the incoming wave, we can figure out from what direction it came. The main detectors, LIGO and VIRGO are currently offline, as they undergo major upgrades which are expected to improve the sensitivity of the detectors in strain and hence distance, by more than a factor of 1 18). These advanced detectors are due to come online in 15, and are fully expected to make the first gravitational wave detection. One of the unavoidable noise sources for ground-based detectors is seismic activity. Together with other noise, this limits the range of the detectors, i.e., they can only see gravitational waves within a certain frequency range. For this reason, there has been a wealth of research into the area of space-based detectors. Such detectors, although free from seismic noise, are still susceptible to noise sources such as detector and acceleration noise shot noise in particular is responsible for the upward slope of all the sensitivity curves as they go towards higher frequencies as is seen in Fig. 1.). Their freedom from seismic noise opens these detectors to gravitational waves in a lower frequency range than their ground-based counterparts. NGO/eLISA New Gravitational-Wave Observatory/evolved Laser Interferometer Space Antenna)) is such a space-based detector. In Fig. 1., we can see the different noise curves attached to the detectors and what types of black hole binaries that they will be able to see. It should be noted that the figure attached is for LISA and not elisa/ngo which has a slightly higher noise curve. We can see from the curve that EMRIs are expected to be seen by LISA. NGO/eLISA is a modified version of the originally planned LISA which, due to cut backs in NASA, had to be redesigned on a smaller budget. It will be up 5

19 1. The Two-Body Problem Figure 1.: The space LISA) and ground based LIGO, VIRGO, GEO6) detectors have different ranges of sensitivity and will therefore see gravitational waves from different sources. Advanced LIGO is also shown as is ET - Einstein s Telescope - the third generation of gravitational wave detectors that will be underground. This figure was taken from 1) for selection as a L mission by the European Space Agency in 15. At the 1 L1 selection process, elisa did not get selected although it was ranked top by the scientific review committee. As the L decision will come after the launch of the LISA pathfinder ) as well as after the activation of advanced LIGO and VIRGO, the gravitational wave community are optimistic that the mission will be selected. 1. The Two-Body Problem The two-body problem in Newtonian theory is readily solvable. An isolated system of two point masses is governed by conserved integrals describing the energy and 6

20 1. The Two-Body Problem Figure 1.3: The orbits of six stars that were tracked with the Very Large Telescope of the European Southern Observatory, Chile. It is clearly seen that all 6 stars are orbiting a central mass that is invisible to the telescope. From the motion of S, the object s mass is estimated as 4.1x1 6 ll solar masses, recent observations have indicated that the radius of the object is no more than 6.5 light years - placing the object in the category of supermassive black hole. This figure was taken from 3) momentum resulting in periodic motion. The two-body problem in general relativity is somewhat different - it is a longstanding open problem going back to work by Einstein himself. With recent advances in gravitational wave detector technology, this age-old problem has been given a new lease of life. Some of the key sources expected to be seen by both space and ground based gravitational wave detectors are black hole binaries BHBs). These can be divided into 3 categories extreme mass ratio inspirals EMRIs), intermediate mass ratio inspirals IMRIs) and comparable mass BHBs. This development is today motivating numerical, analytical and experimental relativists to work together with the prospect of bringing about the reality of gravitational wave astronomy. 7

21 1. The Two-Body Problem 1..1 Black Hole Binary Sources Supermassive black holes SMBHs - black holes with mass approximately 1 6 times that of the sun) are believed to be located at the centre of galaxies; in fact it is known by indirect detection that one resides in the centre of our own galaxy 4). This is very clear in Fig. 1.3 where the orbits of several stars were tracked at the centre of our galaxy. It can be seen that they are all orbiting an invisible body that has dimensions that match that of a SMBH. Near central SMBHs, there are also a disproportionately large number of stellar-mass black holes, which have sunk there through random gravitational encounters. Every now and then, one of these stellar-mass black holes, through interactions with other bodies, will be bumped into the grasp of the SMBH, which will initiate the start of a slow inspiral into the SMBH. These inspirals are known as EMRIs. EMRIs are proving to be one of the most exciting prospects for space-based detectors. The smaller black hole can be expected to complete over 1 5 orbits in the relativistic regime of the Kerr rotating) black hole 5). The resulting emitted radiation will therefore carry information about both the inspiral parameters as well as the space-time geometry that in turn can be used to test General Relativity in the strong field regime. The existence of intermediate-mass black holes IMBHs) with masses ranging from 1 and 1 solar masses has not yet been confirmed but there is evidence that favours their existence 6, 7). These objects are of high astrophysical interest as their existence would impact current understandings of the formation and evolution of both SMBHs and galaxies. IMBHs are believed to reside in the centre of globular clusters GSs), which are difficult to resolve, making detection very difficult. Therefore, a key method of detecting an IMBH could be to detect an IMRI or comparable mass BHBs by use of gravitational wave detectors. IMRIs can be seen as falling into two categories - an IMBH falling into a SMBH that could be detected by space-based gravitational wave detectors or advanced ground-based detectors 8, 9), or a stellar-mass black hole falling into an IMBH, which is expected to be detectable by advanced ground-based detectors 3). IMRIs will also be interesting sources for gravitational wave detectors for similar reasons as EMRIs, they too will experience long inspirals and hence have the potential to reveal information about the space-time geometry of Kerr black holes 8, 9). 8

22 1. The Two-Body Problem Comparable mass BHBs as well as comparable mass compact body binaries are also expected to be key gravitational wave sources for both ground and space based detectors. Stellar mass BHBs are thought to form in GCs through 3 body interactions. Their attractiveness as a source for gravitational wave detection GWD) lies in the fact that they are not strongly bound to the cluster. This implies the possibility of the binary being expelled from the cluster due to interactions with other bodies, resulting in the system evolving in isolation away from the noise of the cluster, which in turn makes them an accessible source of gravitational waves for ground based detectors. SMBH binaries comparable mass BHBs where both black holes are supermassive), on the other hand, are expected to be seen by space-based detectors. SMBH binaries are of great interest to the gravitational wave detection community due to their expectantly large SNR, which should make them detectable with minimal use of data analysis. Accurate models of the inspiral and merger will still be required for using these signals to determine source parameters. 1.. Modelling Techniques Many data analysis techniques currently being used in the search of gravitational waves are based on matched filtering; this allows the extraction of signals buried deep in instrumental noise with significant SNR. For successful detection, matched filtering requires accurate waveform templates. In the case of BHBs, several methods are used to calculate the expected waveforms. Numerical relativity NR) has become an invaluable tool in these calculations; however, it does not come without its constraints. It is extremely computationally expensive and is not suited to BHBs with either a large separation or large mass ratios. In these instances, post-newtonian PN) and gravitational self-force GSF) techniques are required respectively - this sharing of the possible parameter space between the different techniques can be visualised in Fig GSF theory is closely related to black hole perturbation theory and uses a perturbation of Einstein s field equations in the mass ratio to describe the motion of a point particle in a given background space-time. At zeroth order in the small mass ratio, the point mass follows a geodesic of the background. At first order, it deviates from this geodesic due to its interaction with its own field. This deviation 9

23 1. The Two-Body Problem log 1 r 1 /m) 4 3 Post Newtonian Theory Post Newtonian Theory & Perturbation Theory 1 Numerical Relativity Perturbation Theory log 1 m /m 1 ) Figure 1.4: Application of BHB modelling methods for the parameter space of mass ratio m 1 /m and separation r 1. This schematic figure is taken from 31). is interpreted as a force acting on the mass, the so-called GSF. What makes these calculations difficult is that a point mass in curved space-time gives rise to a field that diverges at the particle. It is possible to isolate that part of the physical field that is responsible for its singular behaviour. By subtracting the singular component, the so-called Detweiler-Whiting singular field, from the retarded field, we are left with the regular part, which is by construction) wholly responsible for the selfforce. There are three main approaches to calculating the self-force in practice, and all involve this regular-singular split of the field. PN theory also uses a perturbation of Einsteins field equations by using two parameters - the typical velocity of the system divided by the speed of light) and a measure of the deviation of the curved space-time from a flat space-time i.e. the deviation from the flat metric). At lowest order, PN analysis gives a Newtonian description and general relativistic effects are described as higher order perturbations. 1

24 1. The Two-Body Problem As PN theory is a perturbation in the curvature of space-time and velocity of the system, it is effectively assuming that both are very small parameters, i.e. the theory is only applicable to slow systems in the weak field regime. PN approximations have proven their ability to model comparable mass BHBs as well as IMRIs and EMRIs in the inspiral stage, however PN breaks down at the merger stage where NR is required for solving the final orbits of the binaries. One can clearly see that PN and GSF by their nature are constrained to the modelling of certain systems. GSF requires an extreme mass ratio, while PN is only applicable to slow systems with weak fields, meaning it should not be expected to be very effective in the later stages of BHB inspirals. The word should is intentionally used in this description as PN theory has been applied to strong-field, fastmotion systems like BHBs with remarkable success. By going to higher orders, the PN community has shown impressive results that agree with computationally expensive NR simulations, proving the application of PN in strong fields with fast motions 3, 33, 34). PN theory does eventually become ineffective as the inspiral evolves in BHBs but at a much later point than previously expected. The reason for PNs ability to work outside its expected regime is largely unknown but welcomed by the PN community. GSF is also currently experiencing a similar inexplicable success outside its effective parameter space. Recent advances 35, 36) have shown how GSF can be applied to IMRIs and comparable mass binaries with encouraging results with comparisons to PN and NR. A great consequence of this work is extending the viability of the work from the GSF community to ground-based detector sources, which is also most welcomed by the community. Regardless of the method used, the endgame of BHB modelling is to have a complete waveform template bank available for use by both ground and spacebased detectors. To this end, researchers from each of the areas are beginning to come together to compare the different methods and use them to complement each other, making it a truly global effort to assist in the detection of gravitational waves. 11

25 1.3 The Self Force Problem - Thesis Outline 1.3 The Self Force Problem - Thesis Outline History of Theory This thesis will concentrate on the GSF technique, also known as the self-force problem. As described above, the main problem in this approach lies in the singularity of the field at the particle. Fortunately, producing expressions for such fields is nothing very new the self-acceleration of a charged point particle in flat space-time is given by the well known Abraham-Lorentz-Dirac formula 37). In this scenario, the charge produces a field that acts as radiation, which in turn, diverts the particle from its geodesic for this reason, it became known as the radiation reaction. It was almost three decades later when DeWitt and Brehme derived the formula for the self-force of a charged particle in curved space-time 38), generalizing the results of Dirac et al.. Their calculation did require a minor correction, which was provided by Hobbs several years later 39). It was not until the late 199s, however, that Mino, Sasaki and Tanaka produced the most physically relevant and interesting version of the result that of a point mass in curved space-time 4). This result, also obtained by Quinn and Wald 41) using a different approach, led to the famous MiSaTaQuWa equations, which identified the correct regularisation procedure to remove the problematic singularity. The method they formulated, however, was not practical for calculations, and so, was redesigned by Barack and Ori in 4). Quinn was also the first to produce results in the case of a point scalar charge 43) a simpler model, but one that has been used throughout the community as a test bed for new ideas and methods it is worth noting that Barack and Ori also considered this case initially for their mode-sum scheme 4)). There are several reviews that summarise very well all the work that has been done on this problem in particular those by Poisson 44), Detweiler 45) and Barack 5) Main Approaches The three main methods of calculating the self-force are known as matched expansions, mode sum and effective source. Like most complicated calculations, these GSF approaches are first attempted in toy-models. In the GSF context, the complexity 1

26 1.3 The Self Force Problem - Thesis Outline of the calculation increases with the spin of the field i.e., scalar is considered the simplest, followed by the electromagnetic and gravitational cases. The space-time can also increase in complexity, with the key increase arising from going from nonrotating black holes Schwarzschild, Reissner-Nordström) to rotating black holes Kerr, Kerr-Newman). The matched expansions concept was first suggested by Poisson and Wiseman 46). They suggested matching together two independent expansions for the Green s function one in the quasilocal regime and one in the distant past regime. The quasilocal approach was introduced by Anderson et al. 47, 48), this method uses the MiSaTaQuWa equations to compute the relevent Green s function via an analytic Hadamard expansion. This was built on by Ottewill and Wardell 49, 5), by obtaining a very high order of accuracy from the Hadamard expansion. Joining with Casals and Dolan, they successfully used their results to calculate the self-force on a charged particle, initially in Narai space-time a simple toy black hole space-time) 51), and more recently in Schwarzschild space-time 5). The effective source method was independently proposed by Barack and Goldburn 53, 54) and Detweiler and Vega 55). The methods they used were slightly different, but the concept was very much the same. That was to solve for the fully regularised field from the homogeneous wave equation in the near neighbourhood as well as that of the retarded field outside the near neighbourhood, and uniting the results at the boundary to give that part of the field responsible for the self force. In doing so, they were able to obtain an approximate regularised field that is fully derived from the singular field. The difference of their methods emerged in how they separated the two regions Detweiler and Vega developed the window function which effectively smeared the impact of the singular part of the field from full strength at the particle to zero outside the near neighbourhood; while Barack and Goldbourn introduced a world tube to separate the two regions and imposed boundary conditions to unite them. The most exciting result from the effective source method is the production of an outline to calculate the self-force to second order - a feat that has never before been accomplished, and so, is currently receiving much attention. This has led to another surge in excitement amongst the self force community, as second-order would no doubt lead to more accurate calculations of the self-force and resulting wave-forms 56, 57, 58, 59, 6). 13

27 1.3 The Self Force Problem - Thesis Outline To date, the mode sum method has been the most successful regularisation procedure for calculating the self-force, although the effective source is very clearly catching up. In the mode-sum method, one applies a spherical harmonic decomposition of the singular field; each of the multipole modes is then finite even at the particle, allowing to conveniently subtract the singular field mode by mode. One then numerically calculates the physical field multipoles for input into a mode-sum regularization formula one involving certain analytically given regularization parameters that characterizes the singular behaviour at large multipole numbers. The more regularization parameters one can derive, the faster the convergence of the mode sum becomes. Knowledge of high-order regularization parameters is crucial for assuring the efficiency and accuracy of the GSF calculation. The mode-sum was first introduced by Barack and Ori 4), and further developed by Barack, Ori, Nakamo and Sasaki 61, 6, 63, 64). The development of the Detweiler-Whiting singular field 65) furthered the approach even more, and was followed by a very clear decomposition of the scalar field into mode sums by Detweiler, Whiting and Messaritaki 66). Since its introduction, the mode-sum method has been successfully applied to the more complicated models - including a point electric charge and point mass in Schwarzschild space-time 67, 68), as well as a point scalar charge in Kerr space-time 69). The ultimate goal is to extend this to the astrophysically interesting case of a point mass in Kerr space-time Thesis Outline As self-force plays its part in BHB modelling, and BHB modelling plays its part in the search for gravitational waves, this thesis is also aimed to assist greater goals. We have mostly concerned ourselves with computing the singular field in the different scenarios, and using both the effective source and mode sum methods to obtain results that will assist our fellow researchers. By specialising solely on the singular field, we were able to bring it to an accuracy not conceived possible by even the founders of some of the methods used. To summarise, the results of this thesis enable more accurate and more efficient calculations of the self-force for all researchers in the field, thus making their lives a little bit easier. 14

28 1.3 The Self Force Problem - Thesis Outline Section contains the necessary background for calculating the singular field. This backgroud has been reviewed far more extensively in 44), however the scope of this section is on a need to know basis with respect to the rest of the thesis. Section 3 describes the methods used in calculating the singular field - this was done both covariantly and in coordinates, with both methods having advantages and disadvantages. In the different methods, we also expanded around different points to introduce as much independence as possible for the two methods. I played the the main role for the coordinate results, while my colleague, Barry Wardell, took the lead for the covariant results, which are also in this thesis for completeness. By working in this manner, we could independently check our results and, hence, have great confidence in the results produced. We found both methods produced the same singular field up to an order of ɛ 6, where ɛ is the order of distance in the calculations. A singular field to this accuracy has never before been calculated it assisted us in pushing the boundaries on both the matched expansion and mode-sum methods. Section 4 describes the mode sum method in detail and shows the regularisation parameters that we were able to produce in both Schwarzschild and Kerr spacetimes. These parameters have already been used by several groups and have resulted in self-force calculations to unprecedented accuracy. This work has resulted in over ten parameters, previously unknown, and greatly appreciated by our peers. Section 5 investigates the effective source method. As in the mode sum, we used our high-order singular field to push the boundaries on previous results producing a very smooth field in both Schwarzschild and Kerr space-times. We also extend on the m-mode method, which has evolved from the effective source model, and offer up parameters in both space-times for high-order calculations. The m-mode scheme is an alternative to the mode-sum scheme, introduced for Kerr black holes. It was found that mixing of the modes occurs when calculating the retarded field using the mode-sum method for the gravitational Kerr case, therefore, an alternative that avoids this mixing was introduced in the form of the m-mode method. Previously, researchers only used expansions of the singular field up to ɛ in the m-mode scheme, as the higher orders tend to slow the numerical calculations down. We introduce a method, whereby these higher orders can be used to further regularise the field, without slowing down the numerical calculations. 15

29 1.3 The Self Force Problem - Thesis Outline Section 6 describes further extensions of the high order expansions of the singular field. One of these is the investigation of the cosmic censorship conjecture, which involves the concept of overcharging or overspinning black holes. To assist in these investigations, we produce regularisation parameters for generic motion and radial infall in a spherically symmetric space-time, as well as the motion of a charged particle in Reissner-Nordström space-time. Another extension of this research is in the ongoing work towards calculating the second-order self-force. Such calculations require regularisation parameters of the second derivative of the singular field, which we provide. The final section summarises the results and accomplishments covered in this thesis. We discuss the impact and importance of our results and offer several avenues, down which, this work can be continued. Some parts of this thesis have been in collaboration with both Barry Wardell and my thesis supervisor, Adrian Ottewill. For clarity and completeness, that work has been included here in full. Sections that I was not the primary contributor are indicated by an asterisk *). While the primary focus of this thesis is on computing the singular field for specific space-times, many of the expressions we give are valid in more general spacetimes. In particular, where space allows, we do not make any assumptions about the spacetime being Ricci-flat. To make this distinction explicit, we use the Weyl tensor, C abcd, in expressions which are valid only in vacuum and the Riemann tensor, R abcd in expressions which are also valid for non-vacuum spacetimes. Note that this is done only for space reasons 1 ; our raw calculations include all non-vacuum terms in addition to those given in this thesis and we have made the full expressions available in electronic form 7). Throughout this thesis, we use units in which G = c = 1 and adopt the sign conventions of 71). We denote symmetrization of indices using parenthesis e.g. ab)), anti-symmetrization using square brackets e.g. [ab]), and exclude indices from anti- )symmetrization by surrounding them by vertical bars e.g. a b c), [a b c]). We denote pairwise anti-)symmetrization using an overbar, e.g. R ab cd) = 1R abcd+r cdab ), 1 The notable exception is the case of the gravitational singular field, as in that case the equations of motion have not yet been derived for non-ricci-flat spacetimes. 16

30 1.3 The Self Force Problem - Thesis Outline when multiple symmetries are required. Capital letters are used to denote the spinorial/tensorial indices appropriate to the field being considered. For convenience, we frequently make use of the shorthand notation of 7) by introducing definitions such as R uσuσ uσ = R abcd ef u a σ b u c σ d u e σ f. As is standard practice, commas denote partial differentiation whereas semi-colans represent covariant differentiation, however, these may sometimes be omitted when they are interchangeable, i.e., covariant derivative of a scalar σ ;a = σ,a σ a. 17

31 Chapter Background.1 Bitensors and Basics In this section, we will review the specific biscalars, bivectors and bitensors that are required to fully comprehend this thesis as well as concepts such as geodesics and Penrose diagrams. Throughout, we are primarily dealing with two points - x which is considered to be the source or base point and x, which is a field point, assumed to be in the normal convex neighbourhood of x - this concept will be explained in the next sections..1.1 Geodesics Before we look into the different categories of space-times, it is beneficial to understand how they are represented. Space-times are described by their metric, g ab or line-element ds which are related by ds = g ab x)dx a dx b,.1) where ds can be described as the infinitesimal space-time distance between two neighbouring points x a and x a + dx a. The line element can, therefore, be seen to specify a geometry, although it should be noted that many different line elements can describe the same geometry. The line element can be derived from the Lagrangian 18

32 .1 Bitensors and Basics 73), given by L = 1 g dx a dx b ab dλ dλ,.) where λ is some affine parameter along the geodesic - for time-like geodesics, λ may be proper time, τ. As the line element carries information about the infinitesimal space-time distance between two points, it can be used to determine whether the two points are time-like separated, null separated or space-like separated. If two point particles are time-like separated, it is possible for one particle that in the past of the other), to arrive at the same point in space and time as its partner. If they are null or light-like separated, one can only reach the position of the other in space and time if it can travel at the speed of light. While space-like separated means that unless one particle can travel faster than the speed of light, it can never occupy the same point in space and time as its partner. The line-element, by its nature, can tell us how two points are separated by, > space-like separated ds = null seperated..3) < time-like seprated This concept of separation in space and time can also be illustrated with the use of a light cone. Light cones are merely lines that represent the path of a particle travelling at the speed of light leaving and arriving at a point in space-time. As we take the speed of light c = 1, on a dimensional space-time diagram this represents lines of slope ±1, i.e., those that make a 45 degree angle with the axis. A example of their use to avoid confusion in the observation of events is illustrated in Fig..1. With light cones, when a particle is in the future or past light cone of another, they are said to be time-like separated, if they reside on each others light cones, they are null separated and if they are outside each others light cones, they are space-like separated. In space-time diagrams, the path a particle takes through space and time is known as a world line as it represents the points in space-time that the particle has occupied. Geodesics are world lines that extremise proper time, that is the curve for which an infinitesimal variation in space δx a, produces a vanishing variation in proper time. The flat space-time equivalent of this, is a straight line connecting two 19

33 .1 Bitensors and Basics t! Type%equa!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"! Type%equa!!" =!!!"!!!"!!!!!!!!!!!!! Type%equa!!!!!!!!!!!!!! Figure.1: A, C and D are particles at rest and so travel along straight world lines represented in blue by γ A, γ C and γ D. B starts at C and travels at a finite velocity to D and is represented by the green world line, γ B. When B passes the half way point, where A resides, at t = t O, light is emitted from C to D and from D to C. From above, we can see that A will observe the two light rays to be emitted at the same time t AC, while B will observe light being emitted from D first, then C. On our time line we can clearly see that B will first observe D occurring, then A will witness C and D occurring simultaneously, and then B will observe C happening. In this manner the concept of the light cone can be used to illustrate the relationship of events. points, however in four dimensional space-time, this concept, like many others, is slightly more complicated. If two points are time-like separated, the line element in Eq..1) can be used to describe the proper time between the two points in space-time from dτ = ds, that is τ = xλ1 ) xλ ) gab dx a dx b) 1/

34 .1 Bitensors and Basics = λ1 λ g ab dx a dλ ) dx b 1/ dλ..4) dλ World lines that extremise the proper time between two points must satisfy Lagrange s equations, d dl dλ dẋ dl =,.5) a dxa where the Lagrangian is given by Eq..) and ) refers to differentiation with respect to λ. Some straight forward algebra results in the geodesic equation, d x a dλ + dx a dx b Γa bc dλ dλ =,.6) where the Γ a bc s are called the Christoffel symbols and are given by, Γ a bc = 1 gad g db,c + g dc,b g bc,d ),.7) where A,b implies A x b. The line element can also be used to normalise the four-velocity, which is defined to be u a = dxa dτ..8) From the line element and dτ = ds, it is straight forward to show, g ab u a u b = g ab dx a dτ dx b dτ = 1..9) It is now possible to define the meaning of a normal convex neighbourhood: the normal convex neighbourhood of a point is the set of points that are connected to it by a unique geodesic. If we consider the geodesic which connects x and x, we can use zλ) to represent any point on this geodesic..1. Penrose Diagrams A Penrose diagram can be seen as a coordinate transformation that allows us to view our space-time geometry in a different light, giving us insights into the physical implications of the space-time. In emphasising the light cone structure of spacetime, it successfully maps all of space-time onto a finite space. This is a conformal mapping that compactifies the space-time whilst preserving the light-cone structure. 1

35 .1 Bitensors and Basics!!!!!!! Outgoing& radial& light&rays!!! Ingoing& radial& light&rays Figure.: The Minkowski flat) space-time in null coordinates u = t r and v = t + r. We note that outgoing radial light rays can be described by u = c 1 for any constant, c 1, while ingoing radial light rays are described by v = c 1. Minkowski space-time is limited to t, ) and r, ), this is equivalent to v > u in the above coordinate system. In Sec.., we will describe the various types of space-times that are essential for the understanding of this thesis. However for now we will consider a flat spacetime to introduce the concept of a Penrose diagrams. Flat space-time, known as the Minkowski space-time, is described in Cartesian coordinates, by the line element ds = dt + dx + dy + dz,.1) which can be rewritten in spherical polar coordinates as where we have used the transformation ds = dt + dr + r dθ + sin θ)dφ ),.11) x = r sin θ cos θ, y = r sin θ sin φ, z = r cos θ..1) Introducing null coordinates in the t-r section gives, u t r, v t + r, = t = 1v + u), r = 1 v u)..13)

36 3.1 Bitensors and Basics Figure.3: The function y = tan 1 x maps, ) on to π/, π/). -3 Substituting this transformation into Eq..11) gives a new line element, ds = dudv + 1 u 4 v) dθ + sin θ)dφ )..14) By considering the space time of constant θ and φ it is simple to illustrate the u, v) axes with respect to the x, y) axes, as is done in Fig... Radial light rays can be described as outgoing or ingoing, recalling that we have set c = 1 which implies that light rays are depicted by lines of 45 degrees to the x or r axis, we note that such rays are described by t = r + c 1, or t = r + c 1,.15) for any constant c 1. The slope then tells us if we are dealing with outgoing radial light rays slope = 1) or ingoing light rays slope = 1). Transforming these lines to our u, v) axes shows that outgoing light rays are described by u = c 1, while ingoing light rays are described by v = c 1, as is also depicted in Fig... Another way of finding the angle of radial light rays is to solve ds =, which integrates to give us u = c 1 or v = c 1. Another aspect, to consider, of the new coordinate system is its viable range and domain. In t, r) coordinates, we have t, ) and r, ), depicted as the shaded region in Fig... The equivalent of this in u, v) coordinates is the condition v > u as is easily seen from Fig... To illustrate the Penrose diagram for Minkowski space-time, we introduce another transformation, u tan 1 u 1 t r ), v tan 1 v 1 t + r )..16) 3

37 .1 Bitensors and Basics!!!!!!!!!!!!!!! Outgoing) radial) light)rays!!!!!!!!!!!!!!!!!!!!! Ingoing) radial) light)rays Figure.4: The Penrose diagram for Minkowski space-time is formed by u = tan 1 u and v = tan 1 v, for Minkowski null coordinates u and v. Outgoing and ingoing radial light rays are described by u = c 1 and v = c 1 respectively, for some constant, c 1. Minkowski space-time is limited to t, ), r, ), this is equivalent to v > u and represented by the shaded region. The darker region represents the past light cone for a particle at P, while γ P is its world line. The future and past null infinities are denoted J + and J respectively. I, I + and I are the past timelike infinity, future timelike infinity and spacelike infinity respectively. An immediate consequence of this transformation is that our coordinates now have a finite range due to the finite range of the function tan 1 x illustrated in Fig..3) - all values for u and v must lie in the range π/, π/). In fact, we can limit this further by recalling v > u for Minkowski space-time, from Fig..3, one can clearly see that the immediate implication of this is v > u, this is illustrated as the shaded area in Fig..4, which is also the Penrose diagram for flat space-time. In our null Minkowski coordinates, outgoing light rays were described by u = c 1 which transforms to u = c 1 in the Penrose diagram, where c 1 is also a constant, similarly ingoing light rays are described by v = c 1 v = c 1. This implies that light 4

38 .1 Bitensors and Basics rays are still described as lines parallel to the v, u ) axes or as lines of 45 degrees to the x, y ) axes. This preservation of angles defines the Pemrose diagram as a conformal transformation. The direction of the light rays also gives an immediate meaning to the boundaries. All outgoing light waves u = c 1 will end up on the boundary v = π/, which can now be described as the future null infinity and is denoted by J +. Similarly all ingoing light waves, v = c 1 will originate from the boundary u = π/, known as the past null infinity and denoted J. We can see from Fig..3 that as x, tan 1 x π/, similarly as x, tan 1 x π/. From this we can infer that as v, u), ), v, u ) π/, π/) and as u, v), ), u, v ) π/, π/). If we consider the world line of a particle, P, in Fig..4, denoted by γ P, and follow the particle s world line into its past light-cone, it will, therefore, tend to the point π/, π/) on v, u ). This means that all time-like) world lines originate at this point, which is known as the past time-like infinity, I. Similarly if we follow the particles world line into its future light cone, it will end up at π/, π/). We can, therefore, conclude that all time-like) world lines will end up at this point, known as the future time-like infinity, I +. If we consider space-like curves, we can see that their trajectory will be forced to the point, π/, π/) in v, u ), which is known as spacelike infinity, I. As Penrose diagrams describe a infinite space-time in a finite space, yet maintain the quality that light cones are 45 degrees with the axes, they are very useful in comprehending from which events an observer can receive information. This becomes extremely useful, in particular, for black-hole space-times, although these space-times are more complicated than the flat Minkowski space-time we considered here. We will take a close look at black-hole space-times and their Penrose diagrams in Sec Synge s World Function Synge s world function, σx, x ), is a biscalar defined as one half of the squared geodesic distance between x and x 74). As a biscalar, it holds the ability of a dual definition geometrically. If one was to calculate the derivative of σx, x ), they could do so at either x or x with the resulting vector being very different depending on where the derivative is taken. This is clearly illustrated in Fig..5. 5

39 .1 Bitensors and Basics!!!!!!!!!!!!!!!! Figure.5: The derivatives of sigma at x and x, σ a and σ a respectively. Once differentiated, σ a is a vector with respect to x but still a scalar with respect to x. Similarly, σ a is a vector with respect to x and a scalar with respect to x. This property leads to the ability, on taking further derivatives, of switching the order of primed and unprimed indices with respect to each other with no change to the bitensor, i.e., T ab cd efg h = T acefb d g h note that the indices must stay in order with respect to indices of the same variety - except in the case of the first two due to the scalar nature before derivatives are taken). Mathematically, Synge s world function is represented by, σ x, x ) = 1 λ 1 λ ) λ1 λ g ab z)ż a ż b dλ,.17) where λ affinely parameterises the geodesic connecting x and x, zλ ) = x and zλ 1 ) = x. The geodesic equation gives δ 1 = g ab ż a ż b, Eq..9) where 1 x and x are timelike related δ 1 = x and x are null or lightlike related..18) 1 x and x are spacelike related When x and x are timelike related, λ can be taken to be proper time τ, giving us σ = 1 δ 1 λ = 1 τ.19) To obtain an expression for σ a, we define δσ = σ x + δx, x ) σ x, x ). terms of z, the geodesic connecting x + δx to x is written as zλ) + δzλ), where In 6

40 .1 Bitensors and Basics zλ ) + δzλ ) = x δzλ ) = δx = and zλ 1 ) + δzλ 1 ) = x δzλ 1 ) = δx. δσ is now given by λ1 δσ = 1 [ λ1 λ g ab z + δz) ż a + δż a ) ż b + δż b) dλ λ λ [ g ab z)ż a δż b + 1 g ab,cz)ż a ż b δz c + O δz )] dλ λ1 = λ λ = λ [ g ab ż a δz b] λ 1 λ λ λ1 λ ] g ab z)ż a ż b dλ gab z b + Γ abc ż b ż c) δz a + O δz ) δλ = λg ab ż a δx b,.) where the second equality makes use of the Taylor expansion g ab z + δz) = g ab z) + g ab,c z)δz c + Oδz ) and the third equality involves integration by parts on the first term and a reshuffling of indices. The last equality makes use of the geodesic equation, Eq..6), to make the second term disappear and recalls δzλ ) = while δzλ 1 ) = dx and terms of order δz and higher have been neglected. This gives σ x a σ a x, x ) = λg ab ż b and σ a x, x ) = λż a,.1) where the second identity follows by simply multiplying the first by the inverse metric, g ac. By considering δσ = σ x, x + δx ) σ x, x ) and the geodesic connecting x + δx to x in the above calculation, so that δzλ ) = δx and δzλ 1 ) = δx =, it can also be shown that δσ = λg a b ża δx b σ x a σ a x, x ) = λg a b żb Multiplying Eqs..1) together gives and hence, and σ a x, x ) = λż a..) g ab σ a σ b = λ ż a ż a = λ δ 1 = σ,.3) where we have used Eq..19) in the final equality. Taking the limit of a biscalar, bivector and bitensor as x x is known as the coincident limit. Taking the coincidence limit of σ is easy enough, as we can see directly from Eq..19) that it would be zero, similarly Eq..1) in the coincident limit also gives zero. This is written as, [σ] =, [σ a ] = [σ a ] =..4) 7

41 .1 Bitensors and Basics Taking a derivative of Eq..3) and using Eq..1) gives, σ c = σ b σ bc g cb σ bc ) ż b =..5) As taking the coincident limit is independent of the geodesic, taking the coincident limit of σ ab will be independent of ż b, which leaves [σ ab ] = g a b to be a direct consequence of Eq..5). Similarly it can be found that, [σ ab ] = [σ a b ] = g a b and [σ ab ] = [σ a b] = g a b.6) Differentiating Eq..5) twice more gives σ abc = σ d abcσ d + σ d abσ dc + σ d acσ db + σ d aσ dbc..7) If we take the coincidence limit and use Eqs..6) and.4), we obtain [σ abc ] + [σ acb ] =.8) which can be rearranged to give [σ abc ] = 1 [ R d abc σ d ] =..9) Here, we have used Ricci s identity R d abcσ d = σ abc σ acb, where R d abc is the Riemann curvature tensor, a measure of the space-time curvature, defined by, R a bcd = Γa bd x c Γr bc x d + Γa ceγ e bd Γ a deγ e bc..3) Using Synge s rule 74), [σ...a ] = [σ... ] ;a [σ...a ],.31) it is now straight forward to also calculate [σ abc ] = [σ abc ] = [σ ab c ] = [σ a b c ] =..3) Differentiating Eq..7) again, gives, σ abcd = σ e abcdσ e + σ e abcσ ed + σ e abdσ ec + σ e abσ ecd + σ e acdσ eb + σ e acσ ebd + σ e adσ ebc + σ e aσ ebcd..33) 8

42 .1 Bitensors and Basics If we take the coincidence limit and use Eqs..6),.4) and.3), we arrive at [σ abcd ] + [σ acbd ] + [σ adbc ] =..34) Differentiating and taking the coincidence limit of the Ricci identity already used, R d abcσ d = σ abc σ acb, we arrive at, Substituting these into Eq..34) gives [σ acbd ] = [σ abcd ] R d a b c,.35) [σ adbc ] = [σ abdc ] R c a b d,.36) [σ abdc ] = [σ abcd ]..37) [σ abcd ] = 1 3 R a c b d + R a d b c ),.38) which shows how the Riemann curvature tensor naturally comes about with taking the coincidence limit of higher derivatives of σ. This will be useful when taking covariant Taylor series of σ to obtain our high order covariant expansions in Sec. 3. Using Synge s rule, Eq..31), it is straight forward to also obtain, [σ abcd ] = 1 3 R a c b d + R a d b c ), [σ abc d ] = 1 3 R a c b d + R a d b c ), [σ ab c d ] = 1 3 R a b c d + R a c b d ). [σ a b c d ] = 1 3 R a d b c + R a c b d )..39).1.4 Bivector of Parallel Transport The bivector of parallel transport, g a b, by definition takes a tensor at x and parallel transports it along the geodesic to x, it can also do the opposite, parallel transport from x to x. This is written as, Clearly as x x, g a b δa b, i.e., v a = g a b vb..4) [g a b ] = δa b..41) 9

43 .1 Bitensors and Basics As g a b parallel transports along the geodesic, we immediately have, If we differentiate this, we get g a a ;bσ b = g a a ;b σb =.4) g a b ;cdσ c + g a b ;cσ c d =,.43) which at coincidence gives [g a b ;c] = [g a b ;c ] =,.44) where we have once again made use of Synge s rule, Eq..31)..1.5 The Van Vleck Determinant The Van Vleck determinant is biscalar defined by, x, x ) = det [ a b x, x )] = det σ ab x, x ), g g a b = ga a x, x ) σ a b x, x ).45) where the second identity, when contracted with the bivector of parallel transport, g ca, can be rearranged to give σ ab = g aa a b..46) Taking the coincidence limit of Eq..45) and using Eqs..41) and.6) immediately gives us [ a b ] = δ a b it twice, gives and [ ] = 1. Recalling Eq..3), and differentiating σ ab = σ c aσ cb + σ cab σ c = σ c aσ cb + σ ab cσ c,.47) where the last equality follows as sigma is a biscalar. Combining Eqs..47) and.45) gives, a b = ga a σ ca σ cb + σ a b cσ c ) = g a ag c c σa c c b + a b ;cσ c,.48) 3

44 . Black Hole Space-Times where the second equality follows using Eq..46). Defining the inverse Van Vleck as a b 1 ) b c = δa c and multiplying Eq..48) by 1 ) b d gives δ a d = ga ag c d σa c + 1) b d a b ;cσ c..49) By taking the trace of Eq..49), i.e., setting d a gives, 4 = σ a a + 1) b a a b ;cσ c.5) which can be written as ln ),a σ a = 4 σ a a..51) From Eq..51), we can infer that increases or decreases along each geodesic from x according to whether the rate of divergence of the neighbouring geodesics from x measured by σ a a) is greater or lower than four. It therefore defines a transport equation for. If this divergence is largely negative, we can see that blows up.. Black Hole Space-Times When Einstein first published the full field equations of general relativity, often written as G ab = 8πT ab.5) where G ab and T ab are the Einstein and energy-momentum tensors respectively, they were so complex that he fathomed that it would be a very long time before anybody would be able to produce an exact solution, if at all. Therefore, you can imagine his astonishment when within months, Schwarzschild produced such a solution to the system of equations 75), which would become known as the Schwarzschild solution. Einstein had not known that an exact solution was possible as he, himself, was only able to produce an approximation in the weak field regime that would later become known as post-newtonian theory) to extract values for potential observables. What was more astonishing about this exact solution was that it contained a coordinate singularity. This singularity can be interpreted as a region of space from which nothing can escape, a region caused by an extremely compact object that would later become known as a black hole. 31

45 . Black Hole Space-Times Since the success of Schwarzschild solution, thousands more solutions have been found, however, very few are of actual physical relevance. The concept of a black hole has grown and is now considered a reality by most in the astrophysical world. When considering EMRIs, one can use the model of a test mass representing the stellar mass black hole), orbiting a supermassive black hole which is responsible for the background space-time. From astrophysical considerations, it is believed that this space-time will be either a static, spherically symmetric black hole space-time or the more interesting and more likely candidate of an axially symmetric, spinning black hole space-time. In this section we will look at several space-times that fall into these categories...1 Spherically Symmetric Space-times It can be shown that spherically symmetric vacuum solutions to the Einstein s equations are also static - this is known as Birkhoff s theorem 76). When considering space-times, stationary means that the geodesic is time independent, i.e., g ab x g ab, =.53) where we are taking x = t. In addition, static implies that the metric must remain invariant under time reversal. It, therefore, cannot have any dx α dx terms in the line element, where α {1,, 3}, which in turn implies g α = Schwarzschild Space-Time The Schwarzschild space-time represents the space-time outside a static, spherically symmetric black hole of mass, M. Its line element is given by, ds = 1 M r ) dt + 1 M r ) 1 dr + r dω.54) where dω = dθ + sin θdφ.55) and {t, r, θ, φ} are the standard Schwarzschild coordinates. 3

46 . Black Hole Space-Times The Lagrangian of a particle moving in Schwarzschild space-time is given by Eq..), this can be written as [ L = 1 1 M ) ṫ 1 M ) ] 1 ṙ r r r θ r sin θ φ,.56) where the ) represents differentiation with respect to τ. From Hamiltonian mechanics we have that the canonical momenta, p a = L q a given by p t = L ṫ = 1 M r p θ = L θ = r θ, The Hamiltonian itself is given by ) ṫ, p r = L ṙ = 1 M r where q a = {t, r, θ, φ}, are ) 1 ṙ, pφ = L φ = r sin θ φ..57) H = a p a q a L = L..58) As we are dealing with an isolated system, the Hamiltonian is constant, and therefore from Eq..58), so is the Lagrangian which we can set to be equal to 1, by rescaling the affine parameter, τ, for timelike geodesics. Hamiltonian mechanics also dictates that, L = ṗ a L q a t = p t τ, and L φ = p φ τ,.59) both of which we can set to zero due to the constancy of the Lagrangian. These imply p t = 1 M r ) ṫ = constant E ṫ = E ), 1 M r p φ = r sin θ φ = constant L φ = L r,.6) where we have set θ = π/ in the last equality, and the constants E and L correspond to the energy per unit mass and angular momentum per unit mass respectively. Taking the velocity to be in the plane of θ = π/ can be done without loss of 33

47 . Black Hole Space-Times generality due to the symmetry of the metric, this ensures θ = so the motion stays in that plane. Recalling that we ve set L = 1/, Eq..56) now gives, E ) + 1 M r which can be rearranged to give 1 M r ) 1 ṙ + L = 1.61) r ṙ = E 1 r 3 r M) L + r )..6) We now have expressions for the four velocity of a test particle in Schwarzschild space-time. The Schwarzschild solution describes the space-time of a non-rotating black hole, and as such, is very useful for describing the approach of particles or light rays towards that black hole. However, due to its singular nature at r = M, it is not ideal for understanding the nature of the event horizon or the singularity at r =. The Schwarzschild solution can be rewritten in different coordinates that avoid the r = M singularity and give us a clearer picture of the geometry associated with the region r < M. The Kruskal coordinates are one such transformation - they were introduced earlier in Sec..1., where we used their flat space-time equivalent to investigate the Minkowski space-time as well as to obtain a Penrose diagram for the space-time. As with Minkowski space-time, for the Schwarzschild solution, the Kruskal coordinates only change the r, t) components of the line element or metric. transformation for r > Mis given by, while the transformation for r < M is The r ) ) 1/ t U = M 1 e r/4m) cosh,.63) 4M r ) ) 1/ t V = M 1 e r/4m) sinh,.64) 4M U = V = 1 r M 1 r M ) 1/ t e r/4m) sinh 4M ) 1/ e r/4m) cosh t 4M ),.65) )..66) 34

48 . Black Hole Space-Times r"="m"! V" singularity+ r=)+ r"="1.75m"!! =!"!! =!!! =!"! U"! =!!! r"="m"! singularity+ r=)+! =!"!!! =!" Figure.6: The Kruskal diagram for Schwarzschild space-time is formed by the transformation in Eqs..63) and.65). Lines of constant t and r are highlighted in red and blue respectively. The r = M horizon green) with positive slope also represents t = while the horizon with negative slope also represents t =. t = is the positive U axis for r > M and the positive V axis for r < M. The area that represents only Schwarzschild coordinates is described by U > V and U > V. Regardless of the region, the line element transforms to, ds = 3M 3 e r/m) dv + du ) + r dθ + sin θdφ )..67) r To illustrate how r is represented on a U, V ) diagram, we square our U and V coordinates and subtract one from the other to obtain, r ) > r > M U V = M 1 e r/m) = r = M..68) < r < M This clearly illustrates that lines of constant r are hyperbolas on the U, V ) plane, 35

49 . Black Hole Space-Times in particular, r > M is represented by East-West opening hyperbola, while r < M will be shown as North-South opening hyperbola as is illustrated in Fig..6. We see by increasing r in the r > M region that the hyperbolas move out while decreasing r in the r < M regions has the same effect. It can also be seen from Eq..68) that r = M is not a singularity in Kruskal coordinates, which illustrates that r = M is only a coordinate singularity associated with the Schwarzschild coordinates. In fact, in Kruskal coordinates, it corresponds to the lines U = ±V. Eq..68) also tells us that r = corresponds to the hyperbola V = + U + 1, a hyperbola in the r < M region. The U, V ) plane can clearly be separated into quadrants, as shown in Fig..7, regions I and I represent that area where r > M while areas II and II are where r < M. To investigate the nature of t with respect to the U, V ) plane, we note that the above transformations give, V U = { tanh ) t r > M [ 4M tanh t )] 1 r < M..69) 4M t is, therefore, represented by a straight line in the U, V ) plane, the slope of which will have the same sign as t as is shown in Fig..6. We can also see that t = corresponds to V = for r > M and U = for r < M. By taking the limit of tanh x as x and getting 1, we note that t = corresponds to the line U = V. Similarly taking the limit of tanh x as x to get 1, gives t = corresponding to the line U = V. Schwarzschild coordinates cover the region t, ) and r M, ), immediately we can see that this corresponds to the areas labelled I and I Fig..7. However, further investigation into the transformation given by Eq..63) for r > M shows that U > V, this leaves us with region I, in Fig..7), representing the Schwarzschild coordinates. In Sec..1., we saw that we can investigate how the radial light cone is represented by looking at ds =. In our U, V ) line element Eq..67), this corresponds to, dv = du = V = ±U + c 1,.7) where c 1 is a constant. This means that radial light rays are depicted by straight lines at a 45 degree angle to the U, V ) axes. If we now look at timelike particles and in 36

50 . Black Hole Space-Times V" singularity+r=)+! r"="m" t"== )"! I " II " I" light" signals" observer" U" r"="m"" t"= )"! II " singularity+r=)+ World" line"of" infalling" particle" Figure.7: The Kruskal diagram showing an infalling particle for Schwarzschild space-time. emitting radial light signals at regular intervals of proper time. An observer receives the signals at increasing intervals until the particle crosses the horizon r = M whereupon the observer can no longer receive any light signals. The particle cannot return to region I as this would require it travelling faster than light, i.e., travelling outside its light cone. It will inevitably reach the singularity. their light cones in region I, we can see that their light cones allow them to cross the horizon into region II, however once in region II their future light cones never intersect with region I. This means that no information about the particle can ever be received in region I once the particle crosses the horizon, r = M. In fact once in region II, the particle s light cone ensures that the particle will eventually end up at the singularity, r =. This is clearly seen in Fig..7, where we see a particle releasing light signals to an observer, at a fixed r, at regular intervals of proper time. The signals are received at increasing intervals of proper time by the observer until the particle crosses the horizon whereupon no more information or light signals can 37

51 . Black Hole Space-Times!!!!!!!!!!!!!!!!! r"=""m"horizon)"! Outgoing) radial)!! light)rays!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Ingoing)! radial) light)rays r"=""m" horizon)" Figure.8: The Penrose diagram for Schwarzschild space-time. The singularity at r = is represented by V = ±π/4 yellow). I + and I are the future and past timelike infinities respectively, where all timelike world lines end up. J + and J are the future and past null infinities respectively as all light rays end and originate at these boundaries. Spacelike infinity I, is where all space-like world lines wind up. Once a particle crosses the horizon, the particle cannot return to the region described by Schwarzschild coordinates as this would require it travelling faster than light, i.e., travelling outside its light cone. It will inevitably reach the singularity. be received by the observer. To form the Penrose diagram of Schwarzschild space-time, we do two coordinate transformations from the Kruskal coordinates. We introduce null coordinates by rotating our axes by 45 degrees so light rays are parallel to our axes, U = 1 v u), V = 1 v + u),.71) and then as with the Minkowski space-time, we define our new coordinate system as 38

52 . Black Hole Space-Times u tan 1 u V U, v tan 1 v V + U,.7) U = 1 v u ), V = 1 v + u )..73) As was the case for Minkowski space-time, this transformation successfully maps the infinite ranges of u and v,, ), to the finite range π/, π/) for u and v and radial light rays are depicted by lines of constant u and v. Considering r = in Eq..68) and using Eq..71), we can see that this is equivalent to setting u = 1/v. Using Eq..7) and basic trigonometry gives V = π/4 for v >, while V = π/4 for v < and U = 1 tan 1 v/ 1/v) which, given the restrictions on the range of tan 1, implies U π/, π/). Similarly we can see that r = M, while giving us U = ±V in the Kruskal diagram, or u = and v =, we consequently have u = and v = in the Penrose diagram. This gives us enough to draw our Penrose diagram on the space where v and u both π/, π/), which can be seen in Fig..8. As in flat space we also have different types of infinity, I + and I are the future and past timelike infinities respectively, as they are where all world lines end up. J + and J are the future and past null infinities respectively as all light rays end and originate at these boundaries. We also have spacelike infinity I, where all space-like world lines wind up. There are two sets of all the infinities - one for each asymptotic region. One can clearly see from the Penrose diagram, Fig..8 that once a particle passes the horizon from the region covered by Schwarzschild coordinates into the black hole interior, there is no returning unless it can travel faster than the speed of light, i.e., travel outside its light cone. Similarly no light signals can leave the black hole interior so no information can ever be received from the particle to any observer remaining in the Schwarzschild region. The particle will eventually be forced into the singularity at r =. It should now be obvious that diagrams such as the Kruskal and Penrose diagrams give us further insight into the physical happenings of black holes, in particular, into the space-time geometry of the black hole interior. 39

53 . Black Hole Space-Times..1. Reissner-Nordström Space-Time Reissner-Nordström space-time is a static, asymptotically flat solution of the Einstein- Maxwell equations in general relativity. It describes charged, non-rotating spherical black holes or naked singularities. In general, the Reissner-Nordström space-time represents a gravitating source which is both electrically and magnetically charged. For the purpose of this thesis we need only concern ourselves with an electrically charged source, therefore we will take the magnetic charge to be zero. Both the Schwarzschild and Reissner-Nordström space-times have singularities at their origins. As we have seen in Sec...1.1, in the case of Schwarzschild spacetime, any particle that traverses the event horizon will inevitably be drawn into that singularity. However, in Reissner-Nordström space-time, this is not to be the case. Instead, due to the charge of the black hole, the test) particle can and will leave the vicinity of the singularity, passing through both the event and Cauchy horizons of the Reissner-Nordström black hole to arrive in another universe. In fact, the particle can continue to pass through further event and Cauchy horizons, and so, in a manner, can pass from universe to universe. This feature of the Reissner- Nordström solution is seen clearly in its Penrose diagram in Fig..1, and will be further explained later in this section. The Reissner-Nordström solution is the unique static, spherically symmetric solution of the Einstein-Maxwell equations. In units where the speed of light, Planck s constant, the Boltzmann constant and the Coulomb constant are set to unity, the metric is given by ds = r dt + r dr + r dω,.74) where = r Mr + Q r r )r r + ), dω is the metric on the two-sphere, Q is the electric charge of the black hole and M is the mass of the solution. The field strength, F ab = A b,a A a,b, of the electomagnetic field is produced by the only non-zero component of the vector potential, A t = Q r. There is a curvature singularity ar r = while for < Q < M, r) has two real roots given by r ± = M ± M Q..75) The Cauchy horizon and the event horizon are therefore defined to be at r and r + respectively. Between the horizons, the radial coordinate r becomes timelike and the 4

54 . Black Hole Space-Times time coordinate t spacelike. When Q = M, the horizons degenerate and the black hole is called extremal, while for Q > M the solution has a naked singularity. A test particle with a net electric charge will not describe a geodesic in the Reissner-Nordström due to the electromagnetic forces acting on it. The motion of the particle will be determined by the Lagrangian, L = ) dt r r dτ + qq r dr dτ ) ) dθ r r sin θ dτ ) dφ.76) dτ dt dτ,.77) where q denotes the electric charge per unit mass of the test particle. The equations of motion which follow from this Lagrangian are and dt r dτ + qq r r sin θ dφ dτ ) dr = E qq dτ r r 4 dθ dτ = E = constant.78) = L = constant.79) ) r + L ) r 4,.8) ) = L 1 sin θ L..81) As the motion of such a particle is spherically symmetric, we can say without loss of generality that the motion takes place in a plane which we may take to be θ = π/. This simplifies down Eqs..79) and.81) to dφ dτ = L r and dθ dτ =..8) If we know look at how these particles behave, i.e., a charged particle in Reissner- Nordström space-time, we can see, with the use of mathematical software 77), that particles with the same charge as the black hole are repelled from the horizon of the black hole as illustrated in Fig..9 where the particle is represented by the purple line. In a sense, the charge of the particle protects the particle from being swallowed by the black hole. If we look at the motion of an uncharged and charged with opposite charge to the black hole) particle in Fig..9, represented by the blue 41

55 . Black Hole Space-Times Figure.9: The motion of charged and uncharged particles in Reissner-Nordström space-time is represented by the purple, blue and green lines. The purple line represents a charged particle carrying the same charge as the black-hole which results in it being repelled from the event horizon red circle). The blue line and green lines represent the geodesic motion of an uncharged particle and the motion of a charged particle of opposite charge to the black hole respectively. These pass through both the event horizon and Cauchy horizon orange circle) and then escape out of the black hole by traversing both horizons again. However, the universe into which they re-emerge is not their original universe. As is also shown in Fig..1, this illustrates how a particle can move from universe to universe in Reissner-Nordström space-time. and green lines respectively, we can observe something far more interesting. The particles pass through the event horizon red circle) and the Cauchy horizon orange circle) and then traverses back through both horizons and essentially escape being absorbed by the black hole. Instead they emerge back into the universe, however the universe to which they escape is not the universe from which they came. As described before, this is an example of particles being able to move from universe to universe in the Reissner-Nordström space-time. This adventurous journey can be better explained with the use of a Penrose diagram. We, therefore, once again transform our line element into Kruskal coordinates 4

56 . Black Hole Space-Times to enable us to produce a Penrose diagram. To do this, it is necessary to define the coordinate, r = dr..83) r For our different values of Q, this gives us different definitions after integration, r + r + r + r log r r + ) r r + r log r r ) Q < M r = r + M log [ r M) ] r M ) Q = M..84) r + M log + tan 1 r M Q > M Q M Q M For now we will concentrate on the more interesting case of Q < M. We carry out the cordinate transformation, v = t + r, u = t r,.85) From Eq..84), it is straight forward to calculate dr = drr /, the above transformation therefore gives the line element, ds = r dudv + r dθ + sin θdφ ),.86) which the analogue of the Kruskal solution for Schwarzschild. To obtain a Penrose diagram, it is now necessary to carry out another coordinate transformation, namely, )] )] v = tan [exp 1 r+ r v, u = tan 1 r+ + r [ exp u..87) 4r+ 4r+ Taking the derivatives of either of these, i.e., A = {u, v}, gives, da = ± sec A tan A 4r + ) da r + r 4r = ± csc A ) + r + r ) da,.88) where ± is + for A = v and for A = u. The line element, from Eq..86), can now be given by ds = ) 64r + csc u ) csc v ) du dv + r dθ + sin θdφ ),.89) r r + r 43

57 . Black Hole Space-Times Figure.1: The Penrose diagram for Reissner-Nordström space-time when Q < M, has an infinite number of asymptotically flat regions where r > r +, labelled I, connected by regions II and III where r > r > r + and > r > r respectively. Region III contains a singularity, however it is timelike, which means it can be avoided by a timelike curve, i.e., the world line of a particle. In fact, this singularity acts repulsive in the sense that timelike curves cannot hit them. This Penrose diagram, therefore, hints at a very unusual scenario - it should be possible for a particle, P in Fig..1 to start in region I, cross the horizon r = r + into region II, continue from this region cross another horizon r = r into region III. Here, it will avoid the singularity and can pass through another r = r horizon in the next region II, and again through another r = r + horizon into a new region I. This physically translates to a particle crossing from one universe through a wormhole of sorts and arriving in a new universe. This diagram was taken from 78). 44

58 . Black Hole Space-Times where r is defined by ) tan v tan u r+ r = exp r 4r+ ) r+ r = exp r r r + ) 1/ r r ) r /r +..9) r + The Penrose diagram for this line element has an infinite number of asymptotically flat regions where r > r +, labelled I, connected by regions II and III where r > r > r + and > r > r respectively. Region III contains a singularity, however it is timelike, which means it can be avoided by a timelike curve, i.e., the world line of a particle. In fact, this singularity acts repulsive in the sense that timelike curves cannot hit them. This Penrose diagram, therefore, hints at a very unusual scenario - it should be possible for a particle, P in Fig..1 to start in region I, cross the horizon r = r + into region II, continue from this region cross another horizon r = r into region III. Here, it will avoid the singularity and can pass through another r = r horizon in the next region II, and again through another r = r + horizon into a new region I. This physically translates to a particle crossing from one universe through a wormhole of sorts and arriving in a new universe. However, once the particle has left its original region I, like in the Schwarzschild space-time, no information can be received from it by an observer in I, nor the possibility of return. The Q = M case of the Reissner-Nordström solution can be extended in the same manner as the Q < M case to produce the Penrose diagram of Fig..11. Here we see the we only have one horizon of the type r = m and as in the Q < M, we have a timelike singularity. For the Q > M, we have no horizons which suggests the concept of a naked singularity - a singularity that is not hidden behind an horizon. Such a singularity is thought not to exist in reality 79) - a belief known as the Cosmic Censorship Conjecture, however we will go into this in more detail in Sec If this conjecture is to be believed, this would give us an upper limit on the possible charge of a black hole, M Q..91) This naked singularity can be seen clearly with the Penrose diagam, Fig

59 . Black Hole Space-Times Figure.11: The Penrose diagram for Reissner-Nordström space-time when Q = M, has an infinite number of asymptotically flat regions where allowing the possibility of a particle crossing from one universe through a wormhole of sorts and arriving in a new universe. This diagram was taken from 78). 46

60 . Black Hole Space-Times Figure.1: The Penrose diagram for Reissner-Nordström space-time when the Q > M, we have no horizons which suggests the concept of a naked singularity - a singularity that is not hidden behind an horizon. Such a singularity is thought not to exist in reality Cosmic Censorship Conjecture). If this conjecture is to be believed, this would give us an upper limit on the possible charge of a black hole. This diagram was taken from 78). 47

61 . Black Hole Space-Times..1.3 A General Static, Spherically Symmetric Space-Time In Sec. 6.1, we investigate the singular field and the associated regularization parameters for non-geodesic motion. As our aim is to assist in calculations of the self-force, we tried to keep a very general space-time. To this end we introduce the fr) space-time with the line element, ds = fr)dt + fr) 1 dr + r dω..9) It is easily seen that this space-time is a more general version of both the Reissner- Nordström and Schwarzschild space-times... Axially Symmetric, Stationary Space-Time...1 Kerr Space-Time The Kerr solution describes rotating black holes. In Boyer-Lindquist coordinates, it has the line element, where ds = + 1 Mr ) dt 4aMr sin θ dtdφ + Σ Σ Σ dr + Σdθ [ + Mrr + a ) Σ ] sin θdφ,.93) Σ = r + a cos θ,.94) = r Mr + a,.95) M represents its mass, and a its angular momentum per unit mass. By setting a =, it is easily seen that the Kerr solution reduces to the Schwarzschild solution. As the metric coefficients are independent of t and φ, we know that it is an axially symmetric, stationary solution. As the metric is not time invariant, it is not static, however it is invariant under the simultaneous inversion of t and φ, that is it remains the same under the transformation, t t, φ φ..96) 48

62 . Black Hole Space-Times The consequence of this is that if we move forward in time with a positive spin direction, we will get the same field as if we moved backward in time with a negative spin direction - this, along with other factors, has lead to the understanding that the Kerr solution describes a rotating black hole. As in the Schwarzschild case, we can calculate the Lagrangian for a massive point particle moving in Kerr space-time using Eq..), looking at geodesics in the equatorial plane, i.e., setting θ =, L = 1 [ 1 M ) ṫ + 4aM ṫ r r φ r ṙ ) r + a + a M r φ ]..97) As before, the canonical momenta, p a = L q a, is given by p t = 1 M ) ṫ + am φ = E = constant,.98) r r p φ = am r p r = r.99) ṙ, ) ṫ r + a + a M φ = L = constant,.1) r where the constancy of p t and p φ comes from the independence of the Lagrangian from t and φ. Solving Eqs..98) and.1) for φ and ṫ gives, ṫ = 1 [ ) r + a + a M E am ] [ 1 L, φ = 1 M r r r ) L + am r ] E..11) Using the definition for the Hamiltonian in Eq..58) gives H = 1 1 M ) ṫ + am ṫ r r φ ṙ r 1 ) r + a + a M φ r = 1 {[ 1 M ) ṫ + am ] [ ) φ ṫ r + a + a M ṫ am ] } φ φ + r ṙ r r r r = 1 ) Eṫ L φ r ṙ..1) As H is independent of t we are allowed to set H = 1/ and using Eq..11), it is possible to solve for ṙ, to get, ṙ = 1 [r E + M ] r r ae L) + a E L..13) As in the case of Schwarzschild, by setting τ to be proper time, we have arrived at the four-velocity for geodesic Kerr space-time. 49

63 .3 The Detweiler-Whiting Singular Field... Kerr-Newman Space-Time As the Schwarzschild solution had its charged counterpart in the form of the Reissner-Nordström solution, the Kerr solution has a similar counterpart in the Kerr-Newman solution, which describes the space-time of a charged black hole. In Boyer-Lindquist coordinates, this has the line element, ds = a + z dt + r a )a z ) dtdφ + ρ ρ aρ dr + ρ dθ + a z a ρ a + r ) a z ) ) dφ,.14) where = r Mr + a + Q, ρ = r + a cos θ, z = a cos θ and Q is the charge per unit mass of the black hole. As in the Reissner-Nordström case, the Kerr- Newman solution has a vanishing Ricci scalar and a non-vanishing Ricci tensor, which becomes very useful when considering non-geodesic motion..3 The Detweiler-Whiting Singular Field When considering self-force and the potential fields involved, one is usually restricted to one of three cases, a point particle carrying a scalar or electric charge or a point mass. In each case, a field scalar, electromagnetic or gravitational depending on the case) is produced by the particle, which effects the motion of the particle. This is due to the particle interacting with its own field, causing the particle or mass to deviate from the geodesic of the background space-time. This deviation maybe interpreted as a force acting on the particle or mass, and is the so-called self-force. Calculating the self-force, therefore, requires knowledge of the field that produces it. The problem in producing an expression for such a field is that it is singular at the particle or mass. The traditional way to overcome this is to split the field into a direct and a tail part - such a decomposition originally proved very useful for describing the self-force, but as neither parts were solutions of the field equation, they did not give a meaningful explanation of the self-force. A novel solution to the same problem, first introduced by Detweiler and Whiting 65), is to similarly split the field into two parts - a regular part and a singular part. The singular field is 5

64 .3 The Detweiler-Whiting Singular Field designed so it will solve the same inhomogeneous wave equation as the actual field while the regular field will solve the homogeneous wave equation - in this manner the structure of the field is maintained. This singular-regular split, by design, then gives us two fields - the singular field which has no impact on the motion of the particle but completely contains the singular structure of the original field, and the regular field which is a smooth field, that is completely responsible for the self-force. Throughout this thesis, in the spirit of DeWitt 8), we shall use the notation that ϕ A refers to the field in all 3 cases with A specifying which case, that is Φ scalar field) ϕ A = A a electromagnetic field)..15) h ab gravitational field).3.1 The wave equation For each of the three cases, we have a slightly different version of the inhomogeneous wave equation with a distributional source. We will go through each of the cases here and then give a general expression which can incorporate all three scenarios. This is to allow the reader to become familiar and confident with our more general description of the similar structures associated with the three cases. A massless scalar field in curved space-time will satisfy the inhomogeneous wave equation ξr) Φx) = 4πµx),.16) with the distributional source, µx) = q δ 4 x, z) dτ, where δ 4 x, x ) = δ 4 x x ) = δ 4 x x ),.17) g g γ g ab a b, g ab is the contravariant) metric tensor, g is its determinant at x, g is its determinant at x, a is the covariant derivative defined by a connection A A Ba: a ϕ A = a ϕ A + A A Baϕ B, R is the Ricci scalar, ξ is an arbitrary coupling constant, zτ) describes the world line, q the scalar charge and δ 4 x x ) = δ ) ) ) ) x x δ x 1 x 1 δ x x δ x 3 x 3 is the ordinary coordinate) Dirac delta function. The solution of Eq..16) can be written in terms of a Green s function, G x, x ), 51

65 .3 The Detweiler-Whiting Singular Field Φx) = G x, x ) µx ) g d 4 x = q G x, z) dτ..18) Substituting this back into Eq..16), from this, it becomes clear that for Eq..16) to remain true, the following condition on the Green s function must be satisfied, ξr) Gx, x ) = 4πδ 4 x, x )..19) An electromagnetic field in curved space-time will satisfy, δ a b R a b) A b x) = 4πj a x),.11) in the Lorentz gauge a A a = ), with the distributional source, the current density, given by, j a x) = e g a c x, z) u c δ 4 x, z) dτ.111) γ where δ a b is the Kronecker delta function, R a b is the Ricci tensor, u b is the four velocity and e is the charge of the particle. As with the scalar case, we use a trial solution containing a Green s function of the type, A a x) = G a b x, x ) j b x ) g d 4 x = e G a c x, z) u c dτ,.11) in Eq..11) to produce the required equation for the Green s function, δ a b R a b) G b b x, x ) = 4πg a b x, x ) δ 4 x, x )..113) The propagation of gravitational perturbations in a vacuum space-time is described by, δ a c δ b d + C c a d b ) h cd x) = 16πT ab x),.114) in the Lorentz gauge h ab ;b = ), and with the energy-momentum tensor acting as the distributional source, T ab x) = m γ g a c x, z) g b d x, z) u c u d δ 4 x, z) dτ,.115) 5

66 .3 The Detweiler-Whiting Singular Field where C c a d b is the Weyl tensor and m the mass of the particle. As before we design a solution with the structure, h ab x) = 4 G ab c d x, x ) T c d x ) g d 4 x = 4m G ab cd x, z) u c u d dτ..116) When used in Eq..114), we find the following constraint on the Green s function, δ a c δ b d + C c a d b ) G cd c d x, x ) = 4πg a c x, x ) g b) d x, x ) δ 4 x, x ),.117) where the symmetry brackets in a, b) are introduced to maintain the symmetry on both sides of the equation. As one can see, all three cases are dealt with in the same manner, an attribute which applies to most calculations in this thesis. We therefore, come back to our previous notation introduced in Eq..15). That is, we now consider all of the field equations to be described by, D A Bϕ B [ δ A B P B] A ϕ B = 4πM A,.118) with the distributional source, M described by, M A x) = Q u A δ 4 x, z) dτ,.119) γ where the fields here are all for the scalar, electromagnetic and gravitational cases respectively, we have 1 ξr δ A B = δ a b, P A B = R a b,.1) δ a cδ b) a b) d C c d q 1 Q = e, u A = g a b x, z) u b..11) 4m g a c x, z) g b d x, z) u c u d The solution of these equations can be generally written in terms of a general Green s function, G A B x, x ), as 53

67 .3 The Detweiler-Whiting Singular Field ϕ A = G A B x, x ) M B x ) g d 4 x = Q G A B x, x ) u B dτ,.1) which gives us the general condition for our Green s function, D A BG B C x, x ) = 4πg A C x, x ) δ 4 x, x ),.13) where we have for the scalar, electromagnetic and gravitational waves respectively, G x, x ) 1 G A B x, x ) = G a b x, x ), g A C x, x ) = g G ab c d x, a b x, x )..14) x ) g a c x, x ) g b) d.3. Singular Field in Flat Space-time Before we start looking at the singular field and its associated Green s functions in curved space-time, it is beneficial to examine them in the simpler setting of flat space-time. When we look at the physical solutions to Eq..118), we find there are two key solutions in the form of Eq..1) - the retarded and advanced solutions given by, ϕ A ret)/adv) = G A B ret)/adv) x, x ) M B x ) g d 4 x..15) where G A B ret) x, x ) and G A B adv) x, x ) are the retarded and advanced Green s functions respectively. In this chapter, we will be taking x to be our source point on the world line, γ, and x to be a field point in the near neighbourhood of x. G A B ret) x, x ) is then only non-zero on γ when x is on the future light cone of x, we denote this point as x ret). Similarly G A B adv) x, x ) is only non-zero on γ when x is on the past light cone of x, which we label as x adv). From this nature of the Green s functions, it is quite clear that they have a reciprocity nature, that is G A B ret) x, x ) = G A B adv) x, x)..16) We can see from Figs..13 and.14, that x ret) is in the past of x while x adv) is in 54

68 .3 The Detweiler-Whiting Singular Field!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"#)!! Figure.13: The point at which the retarded Green s function is non-zero on γ is when x is on the future light cone of x, this point is highlighted in green and labelled x ret). its future. When considering the physical implications of our solutions Eq..15), we, therefore, take the retarded solution to be the physically meangingful one. From now on we refer to the field associated with the scalar or electric charge or point mass, through the retarded Green s function, as the retarded field. As was discussed at the beginning of Sec..3, the aim of this section is to obtain the desired regular-singular split of the retarded field of the form, ϕ A ret) = ϕ A S) + ϕ A R)..17) In flat space-time, this proves to be quite simple - we define the singular field to be half the sum of the retarded and advanced fields, while the regular field is half their difference. This give us ϕ A S) = 1 [ ] ϕ A ret) + ϕ A adv), ϕ A R) = 1 [ ] ϕ A ret) ϕ A adv)..18) We can see by applying D A B, that these definitions satisfy the criteria laid out 55

69 .3 The Detweiler-Whiting Singular Field!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"#)!!!!! Figure.14: The point at which the advanced Green s function is non-zero on γ is when x is on the past light cone of x, this point is highlighted in green and labelled x adv). earlier, i.e,, the singular field satisfies the inhomogeneous equation of Eq..118) while the regular field satisfies the homogeneous version of the same equation. When dealing with the singular and regular field, we often think of them in terms of their own Green s functions, therefore we make similar definitions to those in Eq..18) for the corresponding Green s functions. These are G A B S) = 1 [ ] G A B ret) + G A B adv). G A B R) = 1 [ ] G A B ret) G A B adv)..19).3.3 Singular Green s Function in Curved Space-time When we consider the regular-singular split of the retarded field in curved spacetime, our lives are not as easy as they were in flat space-time. In curved space time, energy waves, such as electromagnetic, don t just travel on the light cone Huygens 56

70 .3 The Detweiler-Whiting Singular Field!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"#)!!!!"#)!! Figure.15: The points at which the retarded Green s function is non-zero on γ is when x is in or on the future light cone of x, these points or that part of γ are highlighted in green. Principle) as they do in flat space-time. They scatter off the space-time curvature and so reach points inside the future light-cone. The ramifications of this for our Green s functions is that they now are not only non-zero on the light cone of x, but also within the light cone of x. For the retarded Green s function, this translates to it not being zero when x is on or within the future light cone of x, in other words, when x is in the chronological future of x or x I + x ). Similarly with the advanced Green s function, it is non-zero on γ when x is on or within the past light cone of x, or x I x ), where I x ) is the chronological past of x. These areas of non-zero values are shown in Figs..15 and.16. Now that we have our Green s functions defined we re-examine our flat spacetime definitions of our singular and regular fields. If we maintain these definitions, Eq..18), and hence the associated Green s functions definitions of Eq..19), we will find that both our singular and regular definitions are now non-zero for 57

71 .3 The Detweiler-Whiting Singular Field!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"#)!!!!!!!! Figure.16: The points at which the advanced Green s function is non-zero on γ is when x is in or on the past light cone of x, these points or that part of γ are highlighted in green. x ] [, x ret) xadv), ). In fact as we move our observation point, x, closer to the source represented by the world line, γ, we will notice that the distance between x ret) and x adv) tends to zero. The implication of this on our regular and singular functions are that they will now be dependent on the entire future and past of the x. Not only have our expressions for the singular and regular fields become mathematically impossible to calculate, but physically the system no longer makes sense causally. The conclusion is, therefore, that a different definition of the singular field and the resulting regular field are needed. The required revamp of the singular-regular field split of the retarded field was introduced by Detweiler and Whiting 65) and has since been called the Detweiler- Whiting singular field. The concept is to subtract a function, H A B x, x ), that is a biscalar in the case of the scalar field and a bitensor in the cases of the electromag- 58

72 .3 The Detweiler-Whiting Singular Field netic and gravitational fields, H x, x ) H A B x, x ) = H a b x, x ) H ab c d x, x ) scalar cases) electromagnetic case)..13) gravitational case) from our singular Green s function, i.e., G A B S) x, x ) = 1 [ G A B ret) x, x ) + G A B adv) x, x ) H A B x, x ) ] G A B R) x, x ) = 1 [ G A B ret) x, x ) G A B adv) x, x ) + H A B x, x ) ]..131) To avoid the dependence of our Green s functions on the entire past and future of x, but maintain our previous qualities of the singular field, we demand that H A B x, x ) has the following characteristics, 1. H A B x, x ) must satisfy the homogeneous version of Eq..118), the wave equation. this allows the singular, retarded and regular Green s functions to still satisfy their required wave equations.. H A B x, x ) must still have a reciprocity relation in respect to x and x so that the other Green s functions also maintain their required reciprocity, i.e., g B B H A B x, x ) = g A A HA B x, x) 3. H A B x, x ) = G A B adv) x, x ) when x I x ) 4. H A B x, x ) = G A B ret) x, x ) when x I + x ) The resulting structure can be seen in Figs..17 and.18 It will be shown in Sec..3.4 that within a normal neighbourhood, the retarded and advanced fields can be written in the form G A B ret)/adv) x, x ) = U A B x, x ) δ +/ σ) V A B x, x ) Θ +/ σ).13) where we remind the reader of the definitions and implications of the Dirac delta function, δ σ x, x )) and the Heaviside step function Θ σ x, x )) { δ σ x, x σ x, x ) = null geodesic connecting x and x )) = σ x and x are not connected by their light-cones,.133) 59

73 .3 The Detweiler-Whiting Singular Field!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"#)!!!!!!!"#)!!!"#)! Figure.17: The points at which the singular Green s function is non-zero on γ is when x and x are space and nulllike separated, that part of γ is highlighted in green. Θ σ x, x )) = { σ x, x ) > x and x are spacelike related 1 σ x, x ) x and x are timelike related,.134) where by, we mean that δ s support lies purely on the light-cone. introduce the definitions of δ +/ σ) and Θ +/ σ) as δ + σ x, x )) = { x is on the future light-cone of x elsewhere { x is on the past light-cone of x δ σ x, x )) = elsewhere { Θ + σ x, x )) = 1 σ x, x ) I + x ), Θ σ x, x )) = elsewhere We also,.135),.136) { 1 σ x, x ) I x ) elsewhere,.137) where we can see δ + σ x, x )) + δ σ x, x )) = δ σ x, x )) and Θ + σ x, x )) + Θ σ x, x )) = Θ σ x, x )). 6

74 .3 The Detweiler-Whiting Singular Field!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!"#)!!!!!!!!!"#)!!!"#)!!!! Figure.18: The points at which the regular Green s function is non-zero on γ is when x is not in the chronological past of x, that part of γ is highlighted in green. Recalling the required attributes of H A B x, x ), we examine those listed third and fourth earlier in this section, and using Eq..13), we have { H A B x, U A x B ) = x, x ) δ + σ) V A B x, x ) Θ + σ) x I + x ) U A B x, x ) δ σ) V A B x, x ) Θ σ) x I x ) { V A B = x, x ) x I + x ) V A B x, x ) x I x ),.138) where we note that H A B x, x ) is also defined when x and x are space-like related, however, there are only the first two constraints to impose in that domain. We can therefore, set H A B x, x ) = V A B x, x ) for all values of x, provided V A B x, x ) has the required reciprocity relation and satisfies the homogeneous wave equation of Eq..118), which we will prove is the case in the following section. Substituting this result for H A B x, x ) into our definitions for G A B S) x, x ) and G A B R) x, x ) of Eq..131), we arrive at expressions for both the regular and singlar Green s functions that are valid for all three cases, that is 61

75 .3 The Detweiler-Whiting Singular Field G A B S) x, x ) = 1 [ U A B x, x ) δσ)+v A B x, x ) Θ σ) ],.139) G A B R) x, x ) = 1 U A B x, x ) [δ + σ) δ σ)] V A B x, x ) [Θ + σ) + 1 ] Θ σ),.14) where we have used Θ σ) 1 = Θ σ)..3.4 Restrictions on U A B x, x ) and V A B x, x ) In the previous section, we stated that the retarded and advanced Green s function can be written in the form, G A B ret)/adv) x, x ) = U A B x, x ) δ +/ σ) V A B x, x ) Θ +/ σ)..141) By assuming this form, we can show how it is a suitable representation and deduce expressions and constraints on U A B x, x ) and V A B x, x ). This is done by applying the operator D A B to the Green s function form and proving its viability in Eq..13). An immediate problem arises as distributions δ +/ σ) and Θ +/ σ) are not defined at x = x and cannot be differentiated at this point. Their true meaning becomes clear at the boundary value of a function in the complex plane which we may encapsulate in the prescription where we shift σ σ + ɛ and take the limit from the right as ɛ +. This ensures that we are dealing with a positive ɛ, which in turn demands that x and x are time-like related, allowing us to take Θ +/ σ ɛ) = δ +/ σ + ɛ), where refers to differentiating with respect to σ. We label our shifted Green s functions as G ɛb B ret)/adv) and notice that they satisfy lim ɛ + G ɛb B ret)/adv) = G B B ret)/adv). We also note, by its definition, that σδ σ) = and by differentiating this with respect to σ we also get the identities σδ +/ σ) = δ +/ σ) and σδ +/ σ) = δ +/ σ). Adapting these to the shift that we have applied to σ gives σδ +/ σ + ɛ) = ɛδ +/ σ + ɛ), σδ +/ σ + ɛ) = δ +/ σ + ɛ) ɛδ +/ σ + ɛ), σδ +/ σ + ɛ) = δ +/ σ + ɛ) ɛδ +/ σ + ɛ)..14) 6

76 .3 The Detweiler-Whiting Singular Field Below for simplicity we drop the explicit x, x ) dependence notation and understand indices with imply a dependence on x and those without imply a dependence on x. Applying D A B, as defined in Eq..118), to the proposed shifted Green s function and recalling Eq..3), gives D A BG ɛb B ret)/adv) = G ɛa B ret)/adv) P A BG ɛb B ret)/adv) [ ) = g ab a b U A B δ+/ σ + ɛ) + U A B δ +/ σ + ɛ) σ b b V A B Θ +/ σ ɛ) +V A B δ ] +/ σ + ɛ) σ b [ P A B U A B δ +/ σ + ɛ) V A B Θ +/ σ ɛ) ] = ɛδ +/ σ + ɛ) U A B ɛδ +/ σ + ɛ) V A B + δ +/ σ + ɛ) [ ] U A B );a σa + σ a a 4) U A B + δ +/ σ + ɛ) [+ ] V A B );a σa σ a a) V A B + DA BU B B Θ +/ σ ɛ) D A BV B B..143) By taking the limit, we obtain D A BG ɛb B ret)/adv), D A BG B B ret)/adv) = lim DA BG ɛb B ɛ + ret)/adv) = 4πδ 4 x, x ) U A B + δ +/ σ) [ ] U A B );a σa + σ a a 4) U A B + δ +/ σ) [+ ] V A B );a σa σ a a) V A B + DA BU B B Θ +/ σ) D A BV B B = 4πg A B x, x ) δ 4 x, x ),.144) where the final equality is from equating the right hand side to that of Eq..13). We have also made use of the identities derived in 44), lim ɛδ ɛ + +/ σ + ɛ) =, lim ɛδ ɛ + +/ σ + ɛ) =, lim ɛδ ɛ + +/ σ + ɛ) = πδ 4 x, x )..145) Comparing the final two equalities of Eq..144) and equating coefficients of the distribution functions, we can immediately infer that δ 4 x, x ) U A B = ga B x, x ) δ 4 x, x )..146) 63

77 .3 The Detweiler-Whiting Singular Field As this only gives us information about U A B coincidence limit and denoted by [ U A B ], we have when x = x, also known as the [ ] [ ] U A B = g A B = δ A B..147) From Eqs..1) and.51), we use the identities dx a = σ a /λ)dλ and 1 x, x ) x, x ) σ a ) ;a = 4 to produce, σ a a + 1 σ a ;a = 4 σ a a 4 = 1 λ d dλ = λ d ln,.148) dλ where λ affinely parameterises the geodesic connecting x to x. The coefficient of δ +/ σ) in Eq..144) is to zero to give us ) U A 1 B U A B );a σa + σ a a 4) = λ d ln U A B ln ) =..149) dλ We, therefore, know that U A B ) / is constant on the geodesic connecting x to x. As it is constant we can use our values in the coincidence limit, namely Eq..147) and [ ] = 1 from Sec..1.5, to obtain an expression for U A B, U A B x, x ) = 1 x, x ) g A B x, x )..15) Due to the nature of the delta function, the coefficient of δ +/ σ) in Eq..144) only has to be zero when σ =. We therefore have for σ = [+ V A B );a σa σ a a) V A B + DA BU B B ] which can be rearranged to give σ= =,.151) V A B );a σa V A B 1 σ a,a V A B + DA BU B B =,.15) where we have used Eq..51) rearranged as σ a a = 1 σ a,a. This will be later used to obtain expressions for V A B B Similarly, due to the nature of Θ +/ σ), we also have as described in Sec. 3.1 and Appendices A and D A BV B B =,.153) 64

78 .3 The Detweiler-Whiting Singular Field when x I + x ) I x ), i.e., when x and x are time-like related. note that Eq..153) is one of the requirements that H A B We should satisfies the homogeneous wave equation), described in the previous section, that is necessary for setting H A B = V A B, which we have now shown. As we have established an expressions for U A B which has reciprocity property, as does the retarded and singular Green s functions, it also follows from the definition of Eq..13) that V A B will also have the required reciprocity, i.e., V A B x, x ) = V A B x, x). We have therefore fulfilled all requirements that were necessary for setting H A B with the singular Green s function, G A B S) x, x ) = 1 = V A B, which leaves us [ U A B x, x ) δσ)+v A B x, x ) Θ σ) ],.154) where U A B is defined by Eq..15) and V A B satisfies Eqs..151) and.153)..3.5 The Singular Field The singular field, by design, solves the inhomogeneous wave equation of Eq..118), this means like the retarded field, it can also be written in the form of Eq..1), ϕ A S) = G A B S) x, x ) M B x ) g d 4 x.155) γ Substituting M B x ) from Eq..119) and integrating over x gives ϕ A S) = Q G A B S) x, x ) u B x )dτ..156) γ Recalling Fig..17, we note that the singular Green s function is only supported in between the two points x ret) and x adv) on γ. Combining this with our definition of the singular Green s function, Eq..154), allows us to rewrite Eq..156) as 65

79 .3 The Detweiler-Whiting Singular Field xadv) ϕ A S) = Q G A B S) x, x ) u B x )dτ x ret) = Q xadv) U A B x, x ) δ σ x, x )) u B x )dτ x ret) xadv) + Q V A B x, x ) Θ σ x, x )) u B x )dτ x ret) [ ] U A B x, x ) u B x x ) =x adv) = Q + Q σ a x, x ) u a x ) xadv) x =x ret) V A B x, x ) u B x )dτ.157) x ret) where the integration of the first term follows from dσ = σ a u a dτ with the delta function forcing U A B x, x ) u B x ) onto the light cone, i.e., U A B x, x ) u B x ) is only non-zero for x = x ret) and x = x adv). The second term simply has the Θ σ x, x )) = 1 as its only integrating over that part of γ that is space-like..3.6 The Self-Force and the Singular Field We have seen that in an appropriate gauge, the retarded field, ϕ A x), of an arbitrary point particle satisfies the inhomogeneous wave equation with a distributional source, D A Bϕ B = 4πQ u A δ 4 x, zτ )) dτ,.158) where D A B = δ A B m ) P A B..159) The retarded solutions to this equation gives rise to a field, known as the retarded field which one might naïvely expect to exert a self-force F a = p a Aϕ A ret),.16) on the particle, where p a Ax) is a tensor at x, which depends on the type of charge. In the previous sections we illustrated how Detweiler and Whiting 65) showed how such a singular field can be constructed through a Green function decomposition. In four spacetime dimensions and within a normal neighbourhood, the 66

80 .3 The Detweiler-Whiting Singular Field Green function for the retarded/advanced solutions to Eq..158) may be given in Hadamard form, G ret)/adv) A B x, x ) = θ /+ x, x ) { U A B x, x ) δ [σ x, x )] V A B x, x ) θ [ σ x, x )] },.161) where U A B x, x ) and V A B x, x ) are symmetric bi-spinors/tensors which are regular for x x, defined by Eqs..15) and.15) respectively. The first term here, involving U A B x, x ), represents the direct part of the Green function while the second term, involving V A B x, x ), is known as the tail part of the Green function. Detweiler and Whiting proposed to define a singular Green function by taking A the symmetric Green function, G sym) B = 1G ret) A B + G adv) A B ) and adding V A B x, x ) a homogeneous solution to Eq..158)). This leads to the previously defined singular Green function, G S) A B x, x ) = 1 { U A B x, x ) δ [σ x, x )] + V A B x, x ) θ [σ x, x )] }..16) Note that this has support on and outside the past and future light-cone i.e. for points x and x spatially separated) and is only uniquely defined provided x and x are within a convex normal neighbourhood. Given this singular Green function, we may define the Detweiler-Whiting singular field, ϕ A S) = τret) A G S) B x, zτ )) u B dτ,.163) τ adv) which also satisfies Eq..158). Subtracting this singular field from the retarded field, we obtain the regularized field, ϕ A R) = ϕ A ret) ϕ A S),.164) which Detweiler and Whiting showed gives the correct finite physical self-force, F a = p a Aϕ A R)..165) Moreover, this regularized field is a solution of the homogeneous wave equation, D A Bϕ B R) =..166) This holds independently of whether one is considering a scalar or electromagnetically charged point particle or a point mass. To make this more explicit, in the following subsections we give the form these expressions take in each of scalar, electromagnetic and gravitational cases. 67

81 .3 The Detweiler-Whiting Singular Field Scalar Case In the scalar case the singular field, Φ S), is a solution of the inhomogeneous scalar wave equation, ξr m )Φ S) = 4πq δ 4 x, zτ))dτ,.167) where q is the scalar charge and ξ is the coupling to the background scalar curvature. An expression for Φ S) may be found by considering the scalar Green function obtained by taking U A B = Ux, x ) in Eq..16)), G S) = 1 {Ux, x )δ[σx, x )] + V x, x )θ[σx, x )]},.168) with Ux, x ) = 1/ x, x ) from Eqs..14) and.15), where 1/ x, x ) is the Van Vleck determinant as defined in Eq..45). This Green function is a solution of the equation ξr m )G S) = 4πδ 4 x, x )..169) Given this expression for the Green function, the scalar singular field is Φ S) x) = q G S) x, zτ))dτ = q [ γ Ux, x ) σ c u c ] x =x adv) x =x ret) + q τadv) τ ret) V x, zτ))dτ.17) and one computes the scalar self-force from the regularized scalar field Φ R) = Φ ret) Φ S) as F a = q g ab Φ R),b..171) Electromagnetic Case In Lorenz gauge, the electromagnetic singular field satisfies the equation A S) a R b a A S) b = 4πe g ab x, zτ)) u b δ 4 x, zτ))dτ,.17) where e is the electric charge. An expression for A S) a may be found by considering the electromagnetic Green function obtained by taking U A B = Ux, x ) aa in Eq..16)), G S) aa x, x ) = 1 {Ux, x ) aa δ σx, x )) + V x, x ) aa θ σx, x ))},.173) 68

82 .3 The Detweiler-Whiting Singular Field with U a a = 1/ g a a from Eqs..14) and.15), where ga a is the bi-vector of parallel transport defined in Eq..4). This Green function is a solution of the equation G S) aa R a b G S) ba = 4πg aa x, x )δ 4 x, x )..174) Given this expression for the Green function, the electromagnetic singular field is A S) a = e G S) ab x, zτ))ub zτ))dτ = e [ γ u a U aa x, x ) σ c u c ] x =x adv) x =x ret) + e τadv) τ ret) V ab x, zτ))u b zτ))dτ..175) One computes the electromagnetic self-force from the electromagnetic regular field, A R) a = A ret) a A S) a, as F a = e g ab u c A R) [c,b]..176) Gravitational Case In Lorenz gauge, the trace-reversed singular first order metric perturbation satisfies the equation h S) ab + C a c b d hs) cd = 16πm g a au a g b)b u b δ 4 x, x )dτ,.177) where µ is the mass of the particle and the trace-reversed singular field is related S) to the non-trace-reversed version by h ab = h S) ab 1 hs) g ab with h S) = g ab h S) ab. An expression for may be found by considering the gravitational Green function h S) ab obtained by taking U A B = Ux, x ) aba b in Eq..16)), G S) aba b x, x ) = 1 {Ux, x ) aba b δ [σx, x )] + V x, x ) aba b θ [σx, x )]},.178) with U ab a b = 1/ g a a gb) b solution of the equation from Eqs..14) and.15). This Green function is a G S) aba b + C a p b q G S) pqa b = 4πg a ag b)b δ 4 x, x )..179) 69

83 .3 The Detweiler-Whiting Singular Field Given this expression for the Green function, the trace-reversed singular first order metric perturbation is = 4µ h S) ab = µ [ γ G S) aba b x, zτ ))u a u b dτ u a u b U aba b x, x ) σ c u c ] x =x adv) x =x ret) + µ τadv) τ ret) V aba b x, zτ))ua u b dτ..18) One computes the gravitational self-force from the regularized trace-reversed singular first order metric perturbation, h R) ret) S) ab = h ab h ab, as where F a = µ k abcd hr) bc;d,.181) k abcd 1 gad u b u c g ab u c u d 1 ua u b u c u d ua g bc u d gad g bc..18) 7

84 Chapter 3 High-Order Expansions of the Singular Field 3.1 Coordinate and Covariant Expansions of Fundamental Bitensors In the previous chapter, we gave expressions for the singular field in terms of the bitensors U A B x, x ) and V A B x, x ). The first of these is given by U AB x, x ) = 1/ x, x ) g AB x, x ), 3.1) where x, x ) is the Van Vleck determinant from Eq..45), x, x ) = [ g x)] 1/ det [ σ ;αβ x, x )] [ g x )] 1/ [ ] = det g α α x, x ) σ ;α β x, x ), 3.) g AB is the bi-tensor of parallel transport appropriate to the tensorial nature of the field defined in Eq..14), e.g., 1 scalar) g AB = g ab electromagnetic), 3.3) g a a g b)b gravitational) 71

85 3.1 Coordinate and Covariant Expansions of Fundamental Bitensors and where the higher spin fields are taken in Lorentz gauge. Here, g α α x, x ) is the bivector of parallel transport defined by the transport equation σ α g ab ;α = = σ α g ab ;α. 3.4) The bitensor V AB x, x ) may be expressed in terms of a formal expansion in increasing powers of σ 38, 81): where the coefficients V AB n V AB x, x ) = AB V n x, x ) σ n x, x ), 3.5) n= x, x ) satisfy the recursion relations σ ;α 1/ Vn AB ) ;α + n + 1) 1/ Vn AB + 1 n 1/ D B AC C Vn 1 =, 3.6a) for n N, along with the initial condition σ ;α 1/ V AB ) ;α + 1/ V AB + 1 1/ D B C 1/ g AC ) =, 3.6b) derived from Eq..15). Looking at the above equations for U AB x, x ) and V n AB x, x ), we see that a key component of the present work involves the computation of several fundamental bitensors, in particular, the world function σx, x ), Van Vleck determinant 1/ x, x ), four-velocity u a x), and bivector of parallel transport g a b x, x ). This may be achieved by expressing them as expansions about some arbitrary point x which is close to x and x as shown in figure 3.1. We derive these here using both covariant and coordinate methods, each of which has its own advantages and disadvantages. The covariant expression is more elegant, allowing for compact formulae; however these formulae hide complex terms such as high order derivatives of the Weyl tensor that quickly become extremely time consuming to compute, even using computer tensor algebra packages such as GRTensorII 8) or xact 83). The coordinate approach is less elegant but more practical for explicit calculations and it avoids the need to use tensor algebra. Independently of the approach taken, these expansions may be used to compute expansions of U AB x, x ) and V n AB x, x ) by substituting into the above equations), and hence of the singular field. In the case of covariant 7

86 3.1 Coordinate and Covariant Expansions of Fundamental Bitensors x adv) x x x' x ret) γ Figure 3.1: We expand all bitensors which are functions of x and x ) about the arbitrary point x on the world line. expansions, for explicit calculations one must further expand the covariant expressions in coordinates, yielding an expression which may be directly compared with those obtained from the coordinate approach. The resulting expressions are long but are explicit functions of the coordinates, enabling them to be transformed directly into, for example, C functions, indeed we give them in such form online 7) Coordinate Approach In this section, we describe our method for obtaining coordinate expressions of the biscalars in Eqs..17),.175) and.18). We will start by considering two arbitrary points x and x near x as shown in Fig We will seek expansion where the coefficients are evaluated at x, so we introduce the notation x a = x a xā, δx a = x a x a = x a x a xā, 3.7) where we use the convention that the index carries the information about the point: x a = xā. In the calculations below x a and δx a are both assumed to be small, of 73

87 3.1 Coordinate and Covariant Expansions of Fundamental Bitensors order ɛ. The first item we require for our calculations is a coordinate expansion of the biscalar σ x, x ), the Synge world function. We start with a standard coordinate series expansion about x, see for example 5) note the difference in convention for x a in that paper), to get σx, x ) = 1 g abx)δx a δx b + A abc x)δx a δx b δx c + B abcd x)δx a δx b δx c δx d ) The coefficients are readily determined in terms of derivatives of the metric at x by use of the defining identity derived in Eq..3), σ = σ a σ a, see 5). To be explicit, the first few are given by A abc x) = 1 4 g ab,c)x), B abcd x) = 1 1 g ab,cd)x) 1 4 gpq x) [ g ab, p x)g cd),q x) 1g ab, p x)g q c,d) x) + 36g pa,b x)g q c,d) x) ]. We now go one step further by expanding the coefficients about x to give a double expansion in x a and δx a with coefficients at x. The first few terms are [ σx, x ) = 1g ā b x)δx a δx b 1 + g ā b, c x)δx a δx b x c + Aā b c x)δx a δx b δx c ] [ 1 + g 4 ā b, c d x)δx a δx b x c x d + Aā b c, d x)δx a δx b δx c x d where now we interpret δx a + Bā b c d x)δx a δx b δx c δx d ] + Oɛ 5 ), 3.9) as x a x a xā and we use square brackets to distinguish terms of different order in ɛ. Rather than disturb the flow here and throughout this section we just give the first few terms of each expansion for a general metric to make the structure clear and give explicit expressions in Schwarzschild space-time to much higher order in Appendix A. Now that as the coefficients are at the fixed point x, it is straightforward to take 74

88 3.1 Coordinate and Covariant Expansions of Fundamental Bitensors derivatives of σ at x and x, for example, σ a =gā bδx b + [gā b, c δx b x c + 3Aā b c δx b δx c ] [ 1 + g ā b, c dδx b x c x d + 3Aā b c, dδx b δx c x d + 4Bā b c dδx b δx c δx d ] + Oɛ 4 ), 3.1) [ σ a = σ a + 1g b c,āδx b δx c 1 + g b c, dāδx b δx c δx d + A b c dā δx b δx c δx d ] + Oɛ 4 ), 3.11) ] σ a b = gā b gā b, c x c + [3Aā c d, b 1Bā b c d)δx c δx d 1 g ā b, c d x c x d + Oɛ 3 ). 3.1) Likewise we can calculate the Van Vleck determinant directly from its definition Eq..45), giving 1 x, x ) =1 + 1 x, x ) = ) [ gx)] 1 σa bx, x ) [ gx )] 1 1, 3.13) { 1 g c dgā c, b d gā b, c d g c d,ā b) g c dgē f [ g cē,ā g d f, b + g cā, bgē f, d g dē, f) + gā c,ē g b d, f g b f, d) gā b, c gē f, d g dē, f) ]} δx a δx b + Oɛ 3 ). 3.14) To obtain an expressions at x adv) and x ret), we allow x a to be on the world line and again give it as an expansion around the point x, as shown in Fig Writing x a in terms of proper time τ gives 1 x a τ) = xā + uā τ + 1! uā τ + 1 3!üā τ 3 +, 3.15) where uā is the four velocity at the point xā, τ = τ τ, and an overdot denotes differentiation with respect to τ. 1 In principle this expression is valid for an arbitrary world line. However, later in this thesis, in our explicit calculations, we consider both geodesic and non-geodesic motion. Often we will make the assumption that is is geodesic and derive higher derivative terms from the equations of motion. 75

89 3.1 Coordinate and Covariant Expansions of Fundamental Bitensors We are interested in determining the points on the world line that are connected to x by a null geodesic, that is we want to solve σ x a, x a τ) ) = = 1 g ā buā τ x a )u b τ x b ) + [ 1 g ā buā τ x a ) u b τ + 1 g ā b, cuā τ x a )u b τ x b ) x c g ā b, c x)uā τ x a )u b τ x b )u c τ x c ) ] + Oɛ 4 ). 3.16) By writing τ = τ 1 ɛ + τ ɛ + τ 3 ɛ 3 + and explicitly inserting an ɛ in front of x a, we may equate coefficients of powers of ɛ to obtain which gives, τ 1 + gā buā x b τ 1 gā b x a x b =, 3.17) gā buāuāτ 1 x a )τ = [ 1 g ā buāτ 1 x a ) u bτ g ā b, cuāτ 1 x a )u bτ 1 x b ) x c g ā b, c x)uāτ 1 x a )u bτ 1 x b )u c τ 1 x c ) ]. 3.18) Equation 3.17) is a quadratic with two real roots of opposite sign for x spacelike separated from x) corresponding to the first approximation to our points x adv) and x ret), τ 1± = gā buā x b ± gā buā x b ) + gā b x a x b r 1) ± ρ, 3.19) where r 1) is the leading order term in the coordinate expansion of the quantity r that will be appearing in our covariant expansions in Section Equation 3.18) is typical of the higher order equations giving τ n in terms of lower order terms Covariant Approach * In this section, we briefly discuss our method for obtaining covariant expansions for the biscalars appearing in Eqs..17),.175) and.18). We eventually seek expansions about a point x on the worldline which we may treat as fixed in the majority of this thesis). In doing so, we follow the strategy of Haas and Poisson 44, 7): 76

90 3. Expansions of the Singular Field* For the generic biscalar Ax, zτ)), write it as Aτ) Ax, zτ)). Compute the expansion about τ = τ. This takes the form Aτ) = A τ) + Ȧ τ)τ τ) + 1 Ä τ)τ τ) +, 3.) where Ȧ τ) = A ;āuā, Ä τ) = A ;ā buāu b,. Compute the covariant expansions of the coefficients A τ), Ä τ), about τ. Evaluate the expansion at the desired point, eg. Ax ) = Ax, x ). The resulting expansion depends on τ through the powers of τ τ. Replace these by their expansion in ɛ about x), the distance between x and the worldline. A key ingredient of this calculation is the expansion of τ τ τ in ɛ. The leading orders in this expansion were developed by Haas and Poisson 7) for the particular choices τ + v τ and τ u τ, where u and v are the values of τ at x ret) and x adv) respectively. They found τ ± = r ± s) r ± s) r ± s) [ ] R uσuσ r ± s)ruσuσ u R uσuσ σ +Oɛ 5 ), 3.1) 6 s 4 s where r σāuā and s gā b+uāu b)σāσ b. In Appendix B we extend their calculation to the higher orders required in the present work. In the same Appendix, we also apply the above method to compute covariant expansions of all quantities appearing in the expression for the singular field. In order to obtain explicit expressions, we substitute in the coordinate expansion for σā as discussed in Sec ) along with the metric, Riemann tensor and 4- velocity all evaluated at x). In doing so, we only have to keep terms that contribute up to the required order and truncate any higher order terms. 3. Expansions of the Singular Field* In this section we list the covariant form of the singular field to order ɛ 4, where ɛ is the fundamental scale of separation, so, for example, r, s and σx, x)ā are all of 77

91 3. Expansions of the Singular Field* leading order ɛ. The coordinate forms of these expansions are too long to be useful in print form so instead they are available to download 7) with leading orders given in Appendix B. The structure of the singular field is found to be very similar for all three cases and so the scalar singular field is illustrated in Figs. 3. and 3.3 so the reader can get a feel for this form Scalar singular field To Oɛ 4 ), the scalar singular field is where { Φ S) = q 1 s + r s C 6 s 3 uσuσ + 1 [ ] r 3 s ) rc 4 s 3 uσuσ;u r s )C uσuσ;σ + 1 [ ] Φ S) 36 s + 1 [ ] } 5 3) Φ S) 43 s + 5 4) Oɛ5 ), 3.) [ Φ S) ] 3) = 15 [ r s ] Cuσuσ C uσuσ + s [ r s )3C uσuσ;σσ + 4C uσσā C uσσā) + r 4 and 6 r s 3 s 4 )4C uσuā C uσuā + 3C uσuσ;uu ) + r r 3 s )16C ] { uσuāc uσσā 3C uσuσ;uσ ) + s 4 ā C u u b[ r + s )C σāσ b + r r + 3 s )C uāσ b] ā + C u σ b[ rc σāσ b + r + s )C u bσā ] + r r s + s 4 ā )C uāu bc u u b + r } + s ā ā )C uāσ bc u σ b + C σāσ bc σ σ b 3.3) [ Φ S) ] 4) = 3C uσuσ[ r3 r 4 1 r s + 15 s 4 )C uσuσ;u 3 r s ) C uσuσ;σ ] + s { 3C uσuσ;uuu r r 4 1 r s 15 s 4) + 3 r r 3 s ) C uσuσ;uσσ 3 r s )C uσuσ;σσσ 3 r 4 6 r s 3 s 4 )C uσuσ;uuσ 9 r r 4 1 r s 15 s 4 )C uσuā;u C uσσā[18c uσσā;σ r s ) r r 3 s )1C uσσā;u 16C uσuā;σ 5C uσuσ;ā ) 3C uσuā;u r 4 6 r s 3 s 4 )] C uσuā[36 r r 3 s )C uσσā;σ + r 4 6 r s } 3 s 4 )13C uσuā;σ + 5C uσuσ;ā 5C uσσā;u )] 1 s 4 {C σ ā σ b[c σāσ b;σ rc σāσ b;u C uāσ b;σ ) r + s )C uāσ b;u C uāu b;σ ) r r + 3 s )C uāu b;u ] + C u ā σ b[c σāσ b;σ r + r + s )C uāσ b;σ + C u bσā;σ C σāσ b;u ) + r r + 3 s )C uāu b;σ C uāσ b;u C u bσā;u ) 78

92 3. Expansions of the Singular Field* C uāu bu r r s + s 4 ā )] + C u u b[ r + s )C σāσ b;σ r r + 3 s )C σāσ b;u C uāσ b;σ ) + r r s + s 4 )C uāu b;σ C uāσ b;u ) C uāu b;u r r r s } + 5 s 4 )]. 3.4) 3.. Electromagnetic singular field To Oɛ 4 ), the electromagnetic singular field is A S) a where ā uā = eg a s + 1 [ 6 s s 3 { s 5 [ ] 3 r s Cāuuσ + C uσuσ r s )uā ] } 4 s Cāuuσ;u rcāuuσ;σ ) + [ r r 3 s )C uσuσ;u r s )C uσuσσ uā ] A S) 1 [ ] ) ā + A S) 3) 59 s 5 ā + Oɛ 5 ), 3.5) 4) [A S) ā ] 3) = 1C uσuā;uσ s 4 r + s ) 1C uσuā;σσ r s 4 + 1C u c u dc σ cā d r s 6 4C uσuā C uσuσ r s r 3 s ) 36C uσā c C uσu c s 4 r + s ) and 1C uσuā;uu r s 4 r + 3 s ) + 1 s 6 C u c σ dc u dā c r + C σ dā c ) + 4C u cā dc u c u d s 6 3 r + s ) 1C uāσ c [C uσσ c r s 4 + C uσu c s 4 r + s )] 1C uāu c [C uσσ c s 4 r + s ) + C uσu c r s 4 r + 3 s )] + {8C σ cσ dc σ c σ d s 4 5C c dē fc c dē f s 8 4C uσuσ;uσ r s r 3 s ) + 18C uσu c C uσσ c r s r 3 s ) + 4C uσuσ;σσ s r s ) + 3C uσσ c C uσσ c s r s ) + 1C uσuσ C uσuσ r s ) + 16C u c u dc σ cσ d s 4 r + s ) + 4C uσuσ;uu s r 4 6 r s 3 s 4 ) + 3C uσu c C uσu c s r 4 6 r s 3 s 4 ) + 8C u cu dc u c u d s 4 r r s + s 4 ) + 16 s 4 C u c σ d[c σ cσ d r + C u dσ c r + s )] + 16 s 4 C u cσ d[c u c σ d r + s ) + C u c u d r r + 3 s )]}uā 3.6) [A S) ā ] 4) = 16C uσuā;σσσ r s C uσuσ;σ C uσuā r s r 3 s ) + 7C uσuā;σ C uσuσ r s r 3 s ) 16C uσuā;uσσ s 4 r + s ) + 151C uσā c ;σ C uσu c s 4 r + s ) + 16C uσuā;uuσ r s 4 r + 3 s ) 864C uσā c ;u C uσu c r s 4 r + 3 s ) 16C uσuā;uuu s 4 r r s + s 4 ) 54C uσuσ;u C uσuā s r 4 6 r s 3 s 4 ) + 7C uσuā;u C uσuσ s r r s 79

93 3. Expansions of the Singular Field* + 3 s 4 ) 43 s 6 d C σ cā ;σ C u cu d r + C u dσ c ) 648 s 6 d;σ C u c σ C u dā c r + C σ dā c ) + 16 s 8 C uā c d;ē C u c dē r C uē c d r + C σ c dē C σē c d) + 43 s 4 C uσσ c ;σ [C uāσ c r + C uāu c r + s )] + 54 s 4 C uāσ c ;σ [C uσσ c r + C uσu c r + s )] 34 s 4 C uσσ c ;u [C uāσ c r + s ) + C uāu c r r + 3 s )] + 18 s 4 C uσu c ;σ [5C uāσ c r + s ) + 6C uσā c r + s ) + 3C uāu c r r + 3 s )] 96 s 4 C uāσ c ;u [3C uσσ c r + s ) + 8C uσu c r r + 3 s )] + 96 s 4 C ; c uσuā [3C uσσ c r + s ) + 4C uσu c r r + 3 s )] + 4 s 4 C uāu c ;σ [9C uσσ c r + s ) + 17C uσu c r r + 3 s )] 16 s 4 C uσu c ;u [3C uāσ c r r + 3 s ) + 6C uσā c r r + 3 s ) + C uāu c s 4 r r s + s 4 )] 7 s 4 C uāu c ;u [8C uσσ c r r + 3 s ) + 7C uσu c r r s + s 4 )] 16 s 6 d;σ C u c u [3C σ cā d r + C u cā d3 r + s )] 144 s 6 d C u cā ;σ [3C u dσ c r + C u cu d3 r + s )] + 48 s 6 C d[3c uāσ c; u cσ d r + 3C u dσ c r + 3C σ cσ d + C u cu d3 r + s )] + 16 s 6 d C σ cā ;u [3C u dσ c r + C u cu d3 r + s )] s 6 d;u C u c σ [3C σ dā c r + C u dā c 3 r + s )] + 48 s 6 C d[3 rc uāu c; σ cσ d + 3 rc u cu d r + s ) + C u cσ d3 r + s ) + C u dσ c 3 r + s )] + 16 s 6 d C u cā ;u [3C u cu d r r + s ) + C u dσ c 3 r + s )] s 6 d;u C u c u [3C u cā d r r + s ) + C σ cā d3 r + s )] + {45C c dē f σ C c dē f s 8 45C c dē f u C c dē f r s C uσuσ;uσσ r s r 3 s ) 54C uσuσ;σ C uσuσ r s ) 36C uσuσ;σσσ s r s ) + 36C uσuσ;uuu r s r 4 1 r s 15 s 4 ) 36C uσuσ;uuσ s r 4 6 r s 3 s 4 ) + 18C uσuσ;u C uσuσ 3 r 5 1 r 3 s + 15 r s 4 ) 16 s C uσσ c ;σ [ rc uσu c r 3 s ) C uσσ c r + s )] 7 s 4 d;σ C σ c σ [C u cσ d r + C σ cσ d + C u cu d r + s )] 144 s 4 d;σ C u c σ [C σ cσ d r + C u cσ d r + s ) + C u dσ c r + s ) + C u cu d r r + 3 s )] + 7 s 4 d;u C σ c σ [C σ cσ d r + C u cσ d r + s ) + C u cu d r r + 3 s )] + 6 s C uσσ c ;u [C uσσ c r r 3 s ) + 5C uσu c r 4 6 r s 3 s 4 )] + 7 s C uσu c ;u [3C uσu c r r 4 1 r s 15 s 4 ) + 5C uσσ c r 4 6 r s 3 s 4 )] 7 s 4 d;σ C u c u [C σ cσ d r + s ) + C u cσ d r r + 3 s ) + C u cu d r r s + s 4 )] s 4 d;u C u c σ [C σ cσ d r + s ) + C u cσ d r r + 3 s ) + C u dσ c r r + 3 s ) + C u cu d r r s + s 4 )] 6 s C ; c uσuσ [C uσσ c r r 3 s ) + C uσu c r 4 6 r s 3 s 4 )] 1 s C uσu c ;σ [16C uσσ c r r 3 s ) + 13C uσu c r 4 6 r s 3 s 4 )] + 7 s 4 d;u C u c u [C σ cσ d r r + 3 s ) + C u cσ d r r s + s 4 ) + C u cu d r r r s + 5 s 4 )]}uā. 3.7) 8

94 3. Expansions of the Singular HSL DH-1L HSL DHL !!!!!!!!!!!!! Dr!!!!!!!!!! Dr!!!. HSL DH1L HSL DHL !!! Dr!!! Dr!! HSL DH3L !! Dr!! Dr!! HSL DH5L HSL HSL DH6L Dr! Figure 3.: Terms in the coordinate expansion of the singular field for Oɛ 1 ) top left) to Oɛ 6 ) bottom right). Shown is the scalar case of a circular geodesic of radius r = 1M in Schwarzschild space-time. 81

95 3. Expansions of the Singular Field*!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Figure 3.3: Terms in the coordinate expansion of the singular field, [Φ S) ] n), in the region of the particle for Oɛ 1 ) top left) to Oɛ 6 ) bottom right). Shown is the scalar case of a circular geodesic of radius r = 1M in Schwarzschild space-time. 8

96 3. Expansions of the Singular Field* 3..3 Gravitational singular field To Oɛ 4 ), the gravitational singular field is h S) ab where = µg uāu b aāg b b + 1 [ ] r s )C s 6 s 3 uσuσ uāu b 6 r s C uσuā u b) 6 s 4 Cāu bu + 1 { [ ] 1 s 4 4 s 3 Cāu bu;σ rcāu bu;u ) + 4 s u ā C b)uuσ;u r + s ) ru ā C b)uuσ;σ + uāu b [ r r 3 s )C uσuσ;u r s )C uσuσ;σ ]} + 1 [ hs) ] 648 s 5 ā b 4) + 1 [ hs) 144 s 5 ā b 3) ) + Oɛ 5 ), 3.8) [ hs) ] = 4Cāu bu;σσ s 6 + 4Cāu bu;uσ r s Cā c ā b uσ C buu c r s 6 1Cāuσ c C buσ c s 6 3) + 96Cāu b c C uσu c r s 6 + 4Cāu bu ; c c s 8 + Cāu c dc bu c d s 8 + 4Cā c b d C u cu d s Cāuuσ C buuσ s 4 r + s ) + 4Cāu bu C uσuσ s 4 r + s ) 8Cāu bu;uu s 6 3 r and + s ) 4Cāuu c s 6 [6C buσ c r + C buu c 3 r + s )] + {1C buuσ;σσ r s 4 + 4C buuσ C uσuσ r s r 3 s ) 1C buuσ;uσ s 4 r + s ) 36C b c uσ C uσu c s 4 r + s ) + 1C buuσ;uu r s 4 r + 3 s ) + 1C b c σ dc u cu d r + C u cσ d) s 6 + 1C buσ c [C uσσ c r + C uσu c r + s )] s 4 + 4C b c u d[3c u cσ d r + C u cu d3 r + s )] s 6 + 1C buu c [C uσσ c r + s ) + C uσu c r r + 3 s )] s 4 }uā + {4C σ cσ dc σ c σ d s 4 5C c dē fc c dē f s 8 + 1C uσuσ;σσ s r s ) + 16C uσσ c C uσσ c s r s ) 1C uσuσ;uσ r s r 3 s ) + 64C uσu c C uσσ c r s r 3 s ) + 6C uσuσ r s ) + 8C u c u dc σ cσ d s 4 r + s ) + 1C uσuσ;uu s r 4 6 r s 3 s 4 ) + 16C uσu c C uσu c s r 4 6 r s 3 s 4 ) + 4C u cu dc u c u d s 4 r r s + s 4 ) + 8C u c σ d[c σ cσ d r + C u dσ c r + s )] s 4 + 8C u cσ d[c u c σ d r + s ) + C u c u d r r + 3 s )] s 4 }uāu b 3.9) [ hs) ] = ā b 7Cāu bu;σσσ Cāu bu;uσσ r) s 6 18C uσu c ;σ Cāu b c + Cā c uσ;σ C buu c ) r s 6 4) 7Cāu b c ;σ rc uσu c s 6 9Cāu bu ; c cσ s 8 + 9Cāu bu ; c cu r s 8 d;σ 36C u c u Cā c b d s 8 + 7C u c u d;u rcā c b d s C u c σ d; dcāu b c s C u c u d; c rcāu b d s 8 1Cāu c d ;σ C bu c d s 8 6Cāuσ c; dc bu c d s 8 + 6Cāu c d ;u rc bu c d s 8 ] 83

97 3. Expansions of the Singular Field* d; 18Cā c σ dc buu c s 8 d; 18Cā c u d rc buu c s 8 d 7Cā c b ;σ C u cu d s 8 d ;u rc u cu d s 8 7C uσuσ;σ Cāu bu s 4 r + s ) 18Cāuuσ;σ C buuσ s 4 r + 36Cā c b + s ) 54Cāu bu;σ C uσuσ s 4 r + s ) 7Cāu bu;uuu r s 6 r + s ) + 7C uσuσ;u Cāu bu + 4Cāuuσ;u C buuσ + Cāu bu;u C uσuσ ) r s 4 r + 3 s ) + 54Cāuσ c ;σ C buu c r + C buσ c ) s 6 + 9Cāu bu;uuσ + 6C uσu c ;u Cāu b c + Cā c uσ;u C buu c + 6Cāu b c ;u C uσu c )3 r + s ) s 6 + 6Cāuu c; d3c b cu d r + 3C b cσ d C bu c d r) s Cāu b c; dc u cu d r + C u cσ d) s 8 18Cāuσ c ;u [3C buσ c r + C buu c 3 r + s )] s Cāuu c ;σ [ 3 rc b cuσ + 3 rc buσ c + C buu c 3 r + s )] s 6 + 9Cāu bu ; c [3C uσσ c r + C uσu c 3 r + s )] s Cāuu c ;u [3C b cuσ 3 r + s ) 3C buu c r r + s ) C buσ c 3 r + s )] s 6 + { 36C u c σ d;ē C b c dē s 8 18C buuσ;σσσ r s C u c u d;ē rc b c dē s 8 + 1C b c dē ;σ C u c dē s 8 + 4C b c σ d;ē C u c dē s 8 + 1C b c dē ;u C u c dē r s 8 + 4C b c u d;ē C u c dē r s 8 93C uσuσ;σ C buuσ + 4C buuσ;σ C uσuσ ) r s r 3 s ) + 18C buuσ;uσσ + 7C b c uσ;σ C uσu c ) s 4 r + s ) 18C buuσ;uuσ + 4C b c uσ;u C uσu c ) r s 4 r + 3 s ) + 9 s 3C uσuσ;u C buuσ + 4C buuσ;u C uσuσ ) r 4 6 r s 3 s 4 ) + 18C buuσ;uuu s 4 r r s + s 4 ) 16C u c σ d;σ C b cu d r + C b cσ d) s 6 34C b c σ d;σ C u cu d r + C u cσ d) s 6 +8C bu c d;ē C u c dē r C uē c d r + C σ c dē C σē c d) s 8 16C uσσ c ;σ [C buσ c r + C buu c r + s )] s 4 5C buσ c ;σ [ rc uσσ c + C uσu c r + s )] s C uσu c ;σ [6C b cuσ r + s ) 5C buσ c r + s ) 3 rc buu c r + 3 s )] s C uσσ c ;u [C buσ c r + s ) + C buu c r r + 3 s )] s 4 + 1C buu c ;σ [ 9C uσσ c r + s ) 17 rc uσu c r + 3 s )] s 4 48C buuσ ; c [3C uσσ c r + s ) + 4 rc uσu c r + 3 s )] s C buσ c ;u [3C uσσ c r + s ) + 8C uσu c r r + 3 s )] s C uσu c ;u [3 rc buσ c r + 3 s ) 6C b cuσ r r + 3 s ) + C buu c r r s + s 4 )] s C buu c ;u [8C uσσ c r r + 3 s ) + 7C uσu c r r s + s 4 )] s 4 7C u c u d;σ [3C b cσ d r + C b cu d3 r + s )] s C u c σ d;u [3C b cσ d r + C b cu d3 r + s )] s C u c u d;u [3C b cu d r r + s ) + C b cσ d3 r + s )] s 6 18C b c u d;σ [3C u cσ d r + C u cu d3 r + s )] s 6 + 4C buσ c; d[3c u cσ d r + 3C u dσ c r + 3C σ cσ d + C u cu d3 r + s )] s 6 + 7C b c σ d;u [3 rc u cσ d + C u cu d3 r + s )] s 6 + 7C b c u d;u [3C u cu d r r + s ) + C u cσ d3 r + s )] s 6 + 4C buu c; d[3c σ cσ d r + 3C u cu d r r + s ) + C u cσ d3 r + s ) + C u dσ c 3 r + s )] s 6} uā + { 18C σ c σ d;σ C σ c σ d;u r)c σ cσ d s 4 + 7C c dē f σ C c dē f s 8 18C c dē f u rc c dē f s 8 + 9C uσuσ;uσσ r s r 3 s ) 15C uσuσ ; c rc uσσ c s r 3 s ) 135C uσuσ;σ C uσuσ r s ) + 9C uσuσ;uuu r s r 4 1 r s 15 s 4 ) + 9C uσuσ;uuσ s 3 s r s r 4 ) + 15C uσuσ ; c C uσu c s r r s + 3 s 4 ) + C uσuσ;σσσ 9 r s + 9 s 4 ) + 45C uσuσ;u C uσuσ 3 r 5 1 r 3 s + 15 r s 4 ) 54C uσσ c ;σ [C uσu c r r 3 s ) 84

98 3. Expansions of the Singular Field* + C uσσ c r s )] s + 36C u cσ d s 4 d;u [ C σ c σ d;σ r + C σ c σ r + s )] 18C u cu d s 4 d;σ [C σ c σ r + s ) C σ c σ d;u r r + 3 s d;σ )] 36C u c σ [C σ cσ d r + C u cσ d r + s ) + C u dσ c r + s ) + C u cu d r r + 3 s )] s C uσσ c ;u [C uσσ c r r 3 s ) + 5C uσu c r 4 6 r s 3 s 4 )] s + 18C uσu c ;u [3 rc uσu c r 4 1 r s 15 s 4 ) + 5C uσσ c r 4 6 r s 3 s 4 )] s d;σ 18C u c u [C σ cσ d r + s ) + C u cσ d r r + 3 s ) + 18C u cu d r r s + s 4 )] s 4 d;u + 36C u c σ [C σ cσ d r + s ) + C u cσ d r r + 3 s ) + C u dσ c r r + 3 s ) + C u cu d r r s + s 4 )] s 4 + C uσu c ;σ [ 48C uσσ c r r 3 s ) + 39C uσu c r r s + 3 s 4 )] s d;u + 18C u c u [C σ cσ d r r + 3 s ) + C u cσ d r r s + s 4 ) + C u cu d r r r s + 5 s 4 )] s 4} uāu b. 3.3) 85

99 Chapter 4 Mode-Sum Decomposition of the Singular Field The singular field expansions derived in the previous sections have several applications in explicit self-force calculations. One of the most successful computational approaches to date is the mode-sum scheme of Barack and Ori 4, 61); the majority of existing calculations are based on it in one form or another 66, 67, 68, 69, 7, 84, 85, 86, 87, 88, 89, 9, 91, 9, 93, 94, 95, 96, 97, 98, 99, 1). The basic idea is to decompose the singular retarded field into spherical harmonic modes which are continuous and finite in general for the scalar case and in Lorenz gauge for the electromagnetic and gravitational cases. A key component of the calculation involves the subtraction of so-called regularization parameters - analytically derived expressions which render the formally divergent sum over spherical harmonic modes finite. In this section, we derive these parameters from our singular field expressions and show how they may be used to compute the self-force with unprecedented accuracy. 4.1 Mode Sum Concept The self-force, for each case can be represented by Eq..165), or alternatively as F a = p aa ϕ A R), where ϕ A R) = ϕ A ret) ϕ A S), 4.1) 86

100 4.1 Mode Sum Concept is the regularised field and p aa x) is a tensor at x and depends on the type of charge. We can therefore rewrite the self-force as F a = p aa ϕ A ret) p aa ϕ A S). 4.) Carrying out a spherical harmonic decomposition on the field terms, i.e., ϕ A ret)/s) = ϕ lma ret)/s)y lm, 4.3) lm allows the self-force to be rewritten as, ) F a = paa ϕ lma ret) p aa ϕ lma S). 4.4) lm By defining the l component of the retarded or singular self-force to be, l Faret)/S) l = p aa ϕ lma ret)/s), 4.5) the self-force can be expressed as F a = m= l l F l aret) F l as)). 4.6) It is the last term on the right that we calculate in this chapter for each of the 3 cases, in both Kerr and Schwarzschild space-times. As our singular field is an expansion, it is written in terms of order ɛ, and evaluated at x, that is F a l S) = F l a[-1] r, t ) + F l a[] r, t ) + F l a[] r, t ) + F l a[4] r, t ) + F l a[6] r, t ) +..., 4.7) where we are missing odd orders above 1, as these are zero - this will be shown to be the case later in this chapter. When summed over l, the contribution of F l a[] r, t ) and higher terms to the self-force is also zero. However, if we ignore these higher terms in the approximation of ϕ lma S), our resulting expression for ϕlma R) is only C 1, meaning we can only differentiate it once, which is not sufficient for ϕ lma R) to be a solution of the homogeneous wave equation, Eq..166). When it comes to numerically calculating the self-force using the mode-sum method, the inclusion of the higher order terms dramatically speeds up computation times. For this reason, every extra term or regularisation parameter that can be calculated is of great benefit to the self-force community. 87

101 4. Rotated Coordinates 4. Rotated Coordinates In order to obtain expressions which are readily written as mode-sums, previous calculations 4, 66, 7) found it useful to work in a rotated coordinate frame. We found it most efficient to carry out this rotation prior to doing any calculations. To this end, we introduce Riemann normal coordinates on the -sphere at x in the form α ) α ) w 1 = sin cos β, w = sin sin β, 4.8) where α and β are rotated angular coordinates given by sin θ cos φ = cos α, 4.9) sin θ sin φ = sin α cos β, 4.1) cos θ = sin α sin β. 4.11) In these coordinates, the Schwarzschild metric is given by the line element r M ds = r ) dt + r r M + dw 1 dw [ w1 w 8 w 1 w ) 4 4 w 1 w ) ) { [16 w dr + r 8 w1 w) 4 4 w1 w) ] [ 16 w w1 w) 4 4 w1 w) ] dw 1 ] } dw. 4.1) The algebraic form of the metric makes it very suitable for using with computer algebra programmes such as Mathematica. The apparent complexity of having a non-diagonal metric on S is in fact minimal since the determinant of that metric is simply 1. The Kerr metric in these coordinates is given by the line element [ ] [ ] ds 8Mr 4r = 4r + a w 4 w1 w) 1 dt + a w 4 w1 w + ) dr 4 r Mr + a ) { } amr [8 w 6 w1 w + dtdw )] 1 4 w 1 w [4r + a w 4 w1 w)] { } amrw 1 w 6 w1 w + dtdw ) 4 w 1 w [4r + a w 4 w1 w)] + g w1 w 1 dw 1 + g w1 w dw 1 dw + g w w dw, 4.13) 88

102 4. Rotated Coordinates where g w1 w 1 = g w1 w = g w w = 1 [ w 4 4 w1 w) [4 w 4 w1 w 1w )] 4r + a w 4 w 1 w)] + [ 8 w )] 6 w 1 w {r + a [ ] }) 4 w + Ma r 4 w1 w), 4r + a w 4 w1 w) 1 w 4 4 w1 w) [4 w 4 w1 w)] 1 w w 1 + w 4 ) [ 4r + a w 4 [ )] w1 w)] + w1 w 6 w 1 w) { 8 w 6 w 1 w r + a [ ] }) 4 w + Ma r 4 w1 w), 4r + a w 4 w1 w) w 1 w ) [4 w 4 w 1 w )] w )] + w 1 w ) 6 w 1 w {r + a 4 w 1 w ) [ 4r + a w 4 w 1 [ ] }) 4 w + Ma r 4 w1 w). 4.14) 4r + a w 4 w1 w) As in the Schwarzschild case, this algebraic form has an advantage over its trigonometric counterpart in computer algebraic programmes where trigonometric functions tend to slow calculations down. The complexity of the Kerr metric does slow down the calculation of the singular field. However, despite this, rotating the metric and then calculating the singular field and its resulting regularization parameters still remains faster than calculating the singular field in regular Kerr co-ordinates such as Boyer-Lindquist) and then rotating the resulting complicated expression to obtain the desired regularisation parameters. 89

103 4.3 Mode decomposition 4.3 Mode decomposition The method of regularization of the self force through l-mode decomposition is by now standard, see, for example, 4), 66) and 7). Having calculated the singular field, it is straightforward to calculate the component of the self-force that arises from the singular field 1, F a, for scalar, electromagnetic and gravitational cases using Eqs..171),.176) and.181) with the singular field substituted for the regular field. We study the multipole decomposition of F a by writing F a r, t, α, β) = lm F lm a r, t) Y lm α, β), 4.15) where Y lm θ, φ) are scalar spherical harmonics, and accordingly Fa lm r, t) = F a r, t, α, β) Y lm α, β) dω. 4.16) To calculate the l-mode contribution at x = t, r, α, β ), we have F l a r, t ) = lim r m F lm a r + r, t ) Y lm α, β ). 4.17) In previous calculations, Eq. 4.15) has naturally arisen in Schwarzschild coordinates with θ = π, and it was necessary to perform a rotation to move the coordinate location of the particle from the equatorial plane to a pole in the new coordinate system. However, by choosing to work in an S Riemann normal coordinate system from the start, our particle is already located on the pole. This saves us from further transformation and expansions at this stage. With the particle on the pole, Y lm α =, β ) = for all m. This also allows us, without loss of generality, to take β =. Taking α, β and m all to be equal to zero in Eq. 4.17) gives us F l a r, t ) = lim r = l + 1 4π l + 1 4π F a l,m= r + r, t ) F a r + r, t, α, β) P l cos α) dω. 4.18) lim r 1 In this section, for notational convenience we drop the implied S) superscript denoting singular as we are always referring to the singular component. 9

104 4.3 Mode decomposition For each spin field, the singular self-force, F a r, t, α, β), has the form B a 3n ) F a r, t, α, β) = ρ n+1 ɛn 3, 4.19) n=1 where B k) a = b a1 a...a k x) x a 1 x a... x a k. On identifying τ1 = r 1) ± ρ, this form can be easily seen to follow from the coordinate representation of the above expressions for the singular field. In using Eq. 4.19) to determine the regularization parameters, we only need to take the sum to the appropriate order: n = 1 for A a, n = for B a, etc. Explicitly, in our coordinates ρ = gā buā x b ) + gā b x a x b takes the form ρ r, t, α, β) = E r 3 L r M)) r + ) L + r r r M) w 1 ) Er ṙ r M r + EL w 1 t + Lr ṙ r M r w 1 + E + M ) 1 t + r r w, 4.) in Schwarzschild space-time, and ρ r, t, α, β) = r r [r a E L ) + ML ae) + E r] 3 a Mr + r) [ ] + t w 1 4aM ) EL rer ṙ r a Mr + r ) a + w1 M + a + L + r + r w 1Lrṙ r a Mr + r + t E + M ) 1 + w r r, 4.1) in Kerr space-time, where the α, β dependence is contained exclusively in w 1 and w, and E = u t and L = u φ are the energy per unit mass and angular momentum along the axis of symmetry, respectively. In particular, taking t = t t = ) allows us to write ρ r, t, α, β) = E r 4 L + r ) r M) r + L + r ) w 1 + Lr ṙ r M) L + r ) r ) + r w, 4.) 91

105 4.3 Mode decomposition in Schwarzschild space-time, and ρ r, t, α, β) = in Kerr space-time. r r [Er a + r ) + amae L)] a Mr + r ) [r a + L ) + a M + r 3 ] + a M r + a + L + r) [ rlr 3 w 1 + ṙ a Mr + r) a M + a r + L r + r) 3 + wr, 4.3) For mode-sum decomposition, it is favourable to get ρ α, β) ρ r, t, α, β) in the form ρ α, β) = 1 cos α) ζ 1 k sin β ) 4.4) This is can be done by rewriting Eqs. 4.) and 4.3) with r as where ρ α, β) = ζ w 1 + r w, 4.5) ζ = L + r and ζ = L + r + a M + a 4.6) r in Schwarzschild and Kerr space-times respectively, and rearranging to give ) ] ζ ρ α, β) = 1 cos α) ζ [1 r sin β. 4.7) which is equivalent to Eq. 4.4) with k = ζ r. Defining χβ) 1 k sin β, we ζ can rewrite our w s in the alternate forms ζ ] w1 = 1 cos α) cos β = ρ ζ χ cos β = [k 1 χ)], ζ r) χ 4.8) w = 1 cos α) sin β = ρ ρ ζ χ sin β = 1 χ), ζ r) χ 4.9) Suppose, for the moment, that we may take the limit in Eq. 4.18) through the integral sign, then using our alternate forms we have lim r B 3n ) a ρ n+1 ɛn 3 = b i 1 i...i 3n r ) w i1 w i... w i3n ɛ n 3 = ρ n 3 ρ n+1 ɛ n 3 c an) r, χ). 4.3) ρ 9

106 4.3 Mode decomposition In 6), it was shown that the integral and limit in Eq.4.18) are indeed interchangeable for all orders except the leading order, n = 1 term, where the limiting ρ 3 would not be integrable. Thus we find the singular self-force now has the form [ Fa l r, t ) = l + 1 1) B ɛ a r, t, α, β) lim 4π r ρ 3 r, t, α, β) P l cos α) dω ] + ɛ n 3 ρ n 3 c an) r, χ) P l cos α) dω n= F l a[-1] r, t ) ɛ + F l a[] r, t ) ɛ 1 + F l a[] r, t ) ɛ 1 + F l a[4] r, t ) ɛ 3 + F l a[6] r, t ) ɛ , 4.31) where the β dependence in the c n s are hidden in χ, and the α, β dependence of F a r, t, α, β) is hidden in both the ρ s and c n s. Note here that we use the convention that a subscript in square brackets denotes the term which will contribute at that order in 1/l. Furthermore the integrand in the summation is odd or even under w i w i according to whether n and so 3n ) is odd or even. This means only the even terms are non-vanishing, while F l a[1] r, t ) = F l a[3] r, t ) = F l a[5] r, t ) = etc. Some care is required in dealing with taking the limit in the first term. This has been addressed previously in Schwarzschild space-time 6, 64, 66, 7) which we will now extend to Kerr space-time. As is standard for the Schwarzschild first order, we shift our w 1 coordinate to enable us to remove the cross-terms r w 1, by setting w 1 w 1 + µ r, where µ can easily be read off from Eqs. 4.) and 4.3), to be µ Schwar = µ Kerr = Lr ṙ r M) L + r), 4.3) Lrṙ 3 a Mr + r) a M + a r + L r + r) ) for Schwarzschild and Kerr space-times respectively. This allows us to write ρ r, t, α, β) = ν r + ζ w 1 + r w = ν r + χζ 1 cos α) 4.34) 93

107 4.3 Mode decomposition where the expressions for ν for Schwarzschild and Kerr space-times can be read off from Eqs. 4.) and4.3) respectively. This can be easily rearranged to give ρ r, t, α, β) 3 = ζ 3 χ) 3/ δ + 1 cos α ) 3/ = ζ 3 χ) 3/ l= A 3/ l δ)p l cos α), 4.35) where δ = ν r ζ χ and A 3 l δ) = l + 1 δ 4.36) Here A 3 l δ) is derived from the generating function of the Legendre polynomials as shown in Eq. D1) of 66). ρ r, t, α, β) 3 can now be expressed as ρ r, t, α, β) 3 = 1 ζ νχ l + 1 r ) Pl cos α) 4.37) Bringing this result into our expression for F l a[-1] r, t ) from Eq. 4.31) and integrating over α gives Fa[-1] l r, t ) = 1 l + 1 π ) lim r = l + 1 ) lim r l= 1 ζ ν r bar r ζ ν 1 r π B1) a P l cos α) χ χ 1 dβ = l + ) 1 bar r lim r ζ ν χ 1 r = ) l + 1 bar sgn r) 4.38) ζνr where the first equality takes advantage of the orthogonal nature of the P l cos α) and last equality comes from taking the limit as r and noting from Appendix C of 66) that χ 1 is a special type of hypergeometric function given by χ 1 = F 1, 1 ; 1; k) = dω 1 1 k = ζ r 4.39) B a 1) and b ar also now carry a tilde to signify that they are not the same B a 1) b ar as Eq. 4.19), but rather the tilde represents that they have also undergone the coordinate shift w 1 w 1 + µ r. and 94

108 4.3 Mode decomposition In the higher order tems in Eq. 4.31), we may immediately work with ρ = χζ 1 cos α) so, ρ r, t, α, β) n = ζ n [χ 1 cos α)] n/ = ζ n χ) n/ l= A n/ l )P l cos α), 4.4) where A 1 l ) = from the generating function of the Legendre polynomials and, as derived in Appendix D of 66), for n + 1)/ N P n/ l + 1) l ) = l n) l n + )... l + n) l + n + ), 4.41) A n/ where P n/ = 1) n+1)/ 1+n/ n!!), 4.4) In this case the angular integrals involve 1 π dβ χβ) n/ = χ n/ β) = F 1 n, 1 ), 1, k 4.43) where n + 1)/ N {}. The resulting equations can then be tidied up using the following special cases of hypergeometric functions where Kk) π/ χ 1 χ 1 1 = F 1 k) = F 1, 1 ) ; 1; k = Kk), 4.44) π = F 1 k) = F 1 1, 1 ) ; 1; k = Ek), 4.45) π 1 k sin β) 1/ dβ, Ek) π/ 1 k sin β) 1/ dβ 4.46) are complete elliptic integrals of the first and second kinds, respectively. All other powers of χ can be integrated to hypergeometric functions that can then be manipulated to be one of the above by the use of the recurrence relation in Eq ) of 11); that is F p+1 k) = p 1 p k 1) F p 1k) + 1 p + p 1 ) k F p k). 4.47) p k 1) 95

109 4.3 Mode decomposition RP BO DMW F a[-1] A a A a F a[] B a B a F a[1] C a C a F a[] D a F a[4] Ea 1 F a[6] Ea Table 4.1: Relation between notational choices for the regularization parameters RPs). The most common choices are those of either Barack and Ori 4) or Detweiler, Messaritaki and Whiting 66). In the next sections, we give the results of applying this calculation to each of scalar, electromagnetic and gravitational cases in turn. In doing so, we omit the explicit dependence on l which in each case is F l a[-1] = l + 1)F a[-1], F l a[] = F a[], F l a[] = F l a[6] = F a[] l 1)l + 3), Fa[4] l F a[4] = l 3)l 1)l + 3)l + 5), F a[6] l 5)l 3)l 1)l + 3)l + 5)l + 7). 4.48) It is also worth pointing out that there exists in the literature several different notations for the regularization parameters. We have adopted a notation which is readily extensible to other orders and which makes the dependence on l explicit. To avoid confusion, in Table 4.1 we give the relation between our notation and other common notations. 96

110 4.4 Schwarzschild Space-time 4.4 Schwarzschild Space-time Scalar case In the Schwarzschild scalar case, the regularization parameters for the self-force, as described in Eq..171), are given by F t[-1] = ṙ sgn r) L + r), F Er sgn r) r[-1] = r M)L + r), F θ[-1] =, F φ[-1] =, 4.49) where Er ṙ F t[] = E K), 4.5) πl + r) 3/ 1 { } F r[] = F πr r M)L + r r[]e E + Fr[]K K, 4.51) ) 3/ F E r[] = [E r 3 r M)L + r )], F K r[] = [E r 3 + r M)L + r )], F θ[] =, r ṙ F φ[] = E K), 4.5) LπL + r) 1/ where F t[] = Eṙ E F πrl 4 + r) 7/ t[]e + Ft[]K), K 4.53) F E t[] = 8E L r )r 7 L + r )36L 6 M + 14L 4 Mr + 98L Mr 4 + L r Mr 6 7r 7 ), F K t[] = E r 7 3L 5r ) + r L + r )9L 4 M + 18L Mr + 13Mr 4 r 5 ), where F r[] = 1 πr 6 r M)L + r ) 7/ F E r[]e + F K r[]k), 4.54) 97

111 4.4 Schwarzschild Space-time F E r[] = 8E 4 r 1 L r ) + 4E r 3 L + r )9L 6 M + 6L 4 Mr + 3L Mr 4 + L r Mr 6 3r 7 ) r M)L + r ) 8L 6 M + 8L 4 Mr + 8L Mr 4 L r 5 + 3Mr 6 3r 7 ), F K r[] = E 4 r 1 3L 5r ) E r 5 L + r )18L 4 M + 34L Mr + L r 3 + 3Mr 4 7r 5 ) + r M)r L + r ) 14L 4 M + 8L Mr + 16Mr 4 r 5 ), F θ[] =, 4.55) where F φ[] = ṙ πlr 4 L + r ) 5/ F E φ[]e + F K φ[]k), 4.56) F E φ[] = E r 7 7L r ) + L + r )8L 6 M + 58L 4 Mr + 34L Mr 4 L r 5 + r 7 ), F K φ[] = E r 7 3L r ) r L + r )14L 4 M + 16L Mr + r 5 ), where F t[4] = 3Eṙ E F 4πr 11 L + r) 11/ t[4]e + F K t[4]k), 4.57) Ft[4] E = 3E 4 r 16 3L 4 8L r + 3r) 4 + E rl 5 + r)448l 1 M L 1 Mr + 485L 8 Mr L 6 Mr L 4 Mr L 4 r L Mr 1 13L r Mr r 13 ) L + r) 48L 14 M 978L 1 M r + 851L 1 Mr 3 783L 1 M r L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr 11 15L 4 r L M r L Mr 13 9L r M r Mr r 16 ), F K t[4] = 15E 4 r 16 15L 4 8L r + 31r 4 ) 4E r 7 L + r )11L 1 M L 8 Mr + 687L 6 Mr 4 98

112 4.4 Schwarzschild Space-time L 4 Mr L Mr 8 7L r Mr 1 + 1r 11 ) + r L + r ) 48L 1 M 115L 1 M r + 851L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr L M r L Mr 11 15L r 1 559M r Mr r 14 ), where F r[4] = 3 4πr 13 r M)L + r) F r[4]e E + Fr[4]K), K 4.58) 11/ Fr[4] E = 3E 6 r 19 3L 4 8L r + 3r) 4 E 4 rl 8 + r)896l 1 M L 1 Mr L 8 Mr L 6 Mr L 4 Mr L 4 r L Mr 1 369L r 11 9Mr r 13 ) + 8E rl 3 + r) 51L 14 M 35L 1 M r + 67L 1 Mr L 1 M r L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr L M r L Mr L r M r Mr r 16 ) r M)L + r) 3 496L 14 M 8616L 1 M r L 1 Mr 3 864L 1 M r L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr L 4 r L M r L Mr L r 14 31M r Mr r 16 ), F K r[4] = 15E 6 r 19 15L 4 8L r + 31r 4 ) + E 4 r 1 L + r )448L 1 M L 8 Mr + 798L 6 Mr L 4 Mr L 4 r L Mr 8 171L r Mr r 11 ) E r 5 L + r ) 48L 1 M 1587L 1 M r L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr L M r L Mr 11 45L r 1 558M r Mr r 14 ) + r M)r L + r ) 3 48L 1 M 698L 1 M r + 584L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr 7 99

113 4.4 Schwarzschild Space-time 35956L 4 M r L 4 Mr L M r L Mr L r 1 31M r Mr r 14 ), F θ[4] =, 4.59) where F φ[4] = 3ṙ 4πLr 11 L + r ) 9/ F E φ[4]e + F K φ[4]k), 4.6) F E φ[4] = 15E 4 r 16 43L 4 8L r + 3r 4 ) 1E r 5 L + r )435L 1 M L 1 Mr + 948L 8 Mr L 6 Mr L 4 Mr 8 9L 4 r 9 71L Mr L r 11 9r 13 ) + L + r ) 496L 14 M 9656L 1 M r L 1 Mr L 1 M r L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr L 4 r 1 676L M r L Mr L r M r 14 45r 16 ), F K φ[4] = 15E 4 r 16 L 3r )15L r ) + 1E r 7 L + r )176L 1 M + 635L 8 Mr + 618L 6 Mr L 4 Mr 6 3L Mr L r 9 9r 11 ) r L + r ) 48L 1 M 6648L 1 M r + 584L 1 Mr 3 938L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr L M r L Mr L r M r 1 45r 14 ), where F t[6] = 3Eṙ F E 56πr 16 L + r) 15/ t[6] E + Ft[6]K ) K, 4.61) 1

114 4.4 Schwarzschild Space-time t[6] = 8E 6 r 3 r L)L + r ) ) 11L 4 74L r + 11r 4 ) 5E 4 r 8 L + r 16563L 18 M L 16 Mr L 14 Mr L 1 Mr L 1 Mr L 8 Mr L 6 Mr L 6 r L 4 Mr 14 ) L 4 r L Mr L r Mr r 19 ) + E r 3 L + r 17616L M L 18 M r L 18 Mr L 16 M r L 16 Mr L 14 M r L 14 Mr L 1 M r L 1 Mr L 1 M r L 1 Mr 11 F E L 8 M r L 8 Mr L 6 M r L 6 Mr L 6 r L 4 M r L 4 Mr L 4 r L M r L Mr L r M r r ) + L + r ) L M L M r Mr L 18 M 3 r L 18 M r L 18 Mr L 16 M 3 r L 16 M r L 16 Mr L 14 M 3 r L 14 M r L 14 Mr L 1 M 3 r L 1 M r L 1 Mr L 1 M 3 r L 1 M r L 1 Mr L 8 M 3 r L 8 M r L 8 Mr L 6 M 3 r L 6 M r L 6 Mr L 6 r L 4 M 3 r L 4 M r L 4 Mr L 4 r L M 3 r L M r L Mr + 145L r 1 ), M r Mr 45875r M 3 r t[6] = 875E 6 r 3 15L L 4 r 1531L r r) 6 ) + 5E 4 r 1 L + r 4148L 16 M 3755L 14 Mr L 1 Mr L 1 Mr L 8 Mr L 6 Mr L 4 Mr L 4 r L Mr 14 F K 11

115 4.4 Schwarzschild Space-time ) 4956L r Mr r 17 ) E r 5 L + r 17616L 18 M 3668L 16 M r L 16 Mr L 14 M r L 14 Mr L 1 M r L 1 Mr L 1 M r L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr L 4 r L M r L Mr L r M r Mr 19 ) r ) + r L + r L 18 M L 18 M r L 16 M 3 r L 16 M r L 16 Mr L 14 M 3 r L 14 M r L 14 Mr L 1 M 3 r L 1 M r L 1 Mr L 1 M 3 r L 1 M r L 1 Mr L 8 M 3 r L 8 M r L 8 Mr L 6 M 3 r L 6 M r L 6 Mr L 4 M 3 r L 4 M r L 4 Mr L 4 r L M 3 r L M r L Mr L r M 3 r M r Mr + 77r 1 ), where F r[6] = 3 56πr 18 L + r ) 15/ r M) F E r[6] E + F K r[6]k ), 4.6) F E r[6] = 8E 8 r L)L + r ) 11L 4 74r L + 11r 4 ) r 6 ) 19565r Mr L r L Mr E 6 L + r L 4 r L 4 Mr L 6 r L 6 Mr L 8 Mr L 1 Mr L 1 Mr L 14 Mr L 16 Mr 8816L 18 M ) r 11 ) 19465r Mr 1 465L r E 4 L + r M r L Mr L 4 r L M r L 4 Mr L 6 r 16 1

116 4.4 Schwarzschild Space-time L 4 M r L 6 Mr L 6 M r L 8 Mr L 8 M r L 1 Mr L 1 M r L 1 Mr L 1 M r L 14 Mr L 14 M r L 16 Mr L 16 M r L 18 Mr L 18 M r 153L M ) r 6 ) r Mr 41915L r 1 + E L + r 1998M r M 3 r L Mr + 965L 4 r L M r L M 3 r L 4 Mr L 6 r L 4 M r L 4 M 3 r L 6 Mr L 6 M r L 6 M 3 r L 8 Mr L 8 M r L 8 M 3 r L 1 Mr L 1 M r L 1 M 3 r L 1 Mr L 1 M r L 1 M 3 r L 14 Mr L 14 M r L 14 M 3 r L 16 Mr L 16 M r L 16 M 3 r 4 437L 18 Mr L 18 M r L 18 M 3 r 97856L M r L M 3) r 3 ) r Mr + 975L r 1 + M r ) L + r 51944M r M 3 r L Mr L 4 r L M r L M 3 r L 4 Mr L 6 r L 4 M r L 4 M 3 r L 6 Mr L 6 M r L 6 M 3 r L 8 Mr L 8 M r L 8 M 3 r L 1 Mr L 1 M r L 1 M 3 r L 1 Mr L 1 M r L 1 M 3 r L 14 Mr L 14 M r L 14 M 3 r L 16 Mr L 16 M r L 16 M 3 r L 18 Mr L 18 M r L 18 M 3 r 5848L M r L M 3), F K r[6] = 875E 8 15L r L r 4 L + 47r 6 ) r 6 13

117 4.4 Schwarzschild Space-time 5E ) 6 L + r 755r Mr L r L Mr L 4 r L 4 Mr 1 65L 6 r L 6 Mr L 8 Mr L 1 Mr L 1 Mr L 14 Mr 8816L 16 M ) r 13 ) 3715r Mr L r 18 + E 4 L + r 6948M r L Mr L 4 r L M r L 4 Mr L 4 M r L 6 Mr L 6 M r L 8 Mr L 8 M r L 1 Mr L 1 M r L 1 Mr L 1 M r L 14 Mr L 14 M r L 16 Mr L 16 M r 5388L 18 M ) r 8 E ) L + r r Mr 35875L r M r M 3 r L Mr L 4 r L M r L M 3 r L 4 Mr L 4 M r L 4 M 3 r L 6 Mr L 6 M r L 6 M 3 r L 8 Mr L 8 M r L 8 M 3 r L 1 Mr L 1 M r L 1 M 3 r L 1 Mr L 1 M r L 1 M 3 r L 14 Mr L 14 M r L 14 M 3 r 4 437L 16 Mr L 16 M r L 16 M 3 r 97856L 18 M r L 18 M 3) r 5 + r M) ) L + r 4 715r Mr + 965L r M r M 3 r L Mr L 4 r L M r L M 3 r L 4 Mr L 4 M r L 4 M 3 r L 6 Mr L 6 M r L 6 M 3 r L 8 Mr L 8 M r L 8 M 3 r L 1 Mr L 1 M r L 1 M 3 r L 1 Mr L 1 M r L 1 M 3 r L 14 Mr L 14 M r L 14 M 3 r L 16 Mr L 16 M r L 16 M 3 r 194L 18 M r L 18 M 3) r, 14

118 4.4 Schwarzschild Space-time F θ[6] =, 4.63) where F φ[6] = 3ṙ F E 56πLr 16 L + r) 13/ φ[6] E + Fφ[6]K ) K, 4.64) Fφ[6] E = 875E 6 r L 4 r 337L 6 947L r r) 6 ) + 175E 4 r 8 L + r 9834L 18 M L 16 Mr L 14 Mr L 1 Mr L 1 Mr L 8 Mr L 6 Mr 1 + 5L 6 r L 4 Mr L 4 r 15 ) + 496L Mr L r 17 5r 19 ) + E r 3 L + r 153L M L 18 M r L 18 Mr L 16 M r L 16 Mr L 14 M r L 14 Mr L 1 M r L 1 Mr L 1 M r L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 6 r L 4 M r L 4 Mr L 4 r L M r 18 94L Mr L r 3148M r ) r ) L + r L M r L M L 18 M 3 r L 18 M r L 18 Mr L 16 M 3 r L 16 M r L 16 Mr L 14 M 3 r L 14 M r L 14 Mr L 1 M 3 r L 1 M r L 1 Mr L 1 M 3 r L 1 M r L 1 Mr L 8 M 3 r L 8 M r L 8 Mr L 6 M 3 r L 6 M r L 6 Mr 16 7L 6 r L 4 M 3 r L 4 M r L 4 Mr L 4 r L Mr L M 3 r L M r M 3 r 458M r r 3 ), 575L r 1 15

119 4.4 Schwarzschild Space-time Fφ[6] K = 875E 6 r 3 15L L 4 r 587L r r) 6 ) L + r 4915L 16 M 3948L 14 Mr 814L 1 Mr 4 758L 4 Mr E 4 r L 1 Mr L 8 Mr L 6 Mr 1 ) + 385L 4 r L Mr L r r 17 ) E r 5 L + r 17616L 18 M L 16 M r L 16 Mr L 14 M r L 14 Mr L 1 M r L 1 Mr L 1 M r L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr L 4 r L M r L Mr 17 ) 45675L r M r r ) + r L + r L 18 M L 18 M r L 16 M 3 r L 16 M r L 16 Mr L 14 M 3 r L 14 M r L 14 Mr L 1 M 3 r L 1 M r L 1 Mr L 1 M 3 r L 1 M r L 1 Mr L 8 M 3 r L 8 M r L 8 Mr L 6 M 3 r L 6 M r L 6 Mr L 4 M 3 r L 4 M r L 4 Mr 16 35L 4 r L M 3 r L M r L r M 3 r M r r L Mr 18 ) Electromagnetic case In the electromagnetic case, an ambiguity arises in the definition of u a in the angular directions away from the world-line. In.176) one is free to define u a x) as they wish provided lim x x u a x) = uā. A natural covariant choice would be to define this through parallel transport, u a x) = g a bu b. However, in reality it is more practical in numerical calculations to define u a such that its components in Schwarzschild coordinates are equal to the components of uā in Schwarzschild coordinates 68). 16

120 4.4 Schwarzschild Space-time Doing so, the regularization parameters are given by F t[-1] = ṙ sgn r) L + r ), F r[-1] = Er sgn r) r M)L + r ), F θ[-1] =, F φ[-1] =, 4.65) Eṙ F t[] = r πr r + L ) 3/ K + L E ), 4.66) where F r[] = 1 πr 3 r + L ) 3/ r M) F E r[] E + F K r[]k ), 4.67) Fr[] E = E L r 3 + L + r) L + r) M r ), ) Fr[] K = E r 5 + r L + r r M), F θ[] =, 4.68) F φ[] = ṙ πlr L + r [ ) ] E L + r Kr, 4.69) where Eṙ F t[] = F E πr 4 L + r) 7/ t[] E + Ft[]K ) K, 4.7) Ft[] E = E r 5 L 4 + 1L r + 3r) 4 + L + r) 6L 6 M + 168L 4 Mr + 18L Mr 4 13L r Mr 6 5r) 7, 11L + 3r) F K t[] = E r 7 ) ) + r L + r 1L 4 M 48L Mr + 3L r 3 3Mr 4 + r 5, where 1 F r[] = F E πr 6 L + r) 7/ r[] E + Fr[]K ) K, 4.71) r M) 17

121 4.4 Schwarzschild Space-time F E r[] = E 4 r 8 + E r 3 L 4 + 1L r + 3r 4 ) ) L + r 3L 6 M 86L 4 Mr + L 4 r 3 98L Mr 4 + 1L r 5 6Mr 6 + r 7 ) L + r) M r ) 44L 6 M + 94L 4 Mr + 54L Mr 4 + L r 5 + 3r) 7, 11L + 3r) Fr[] K = E 4 r 1 + E r 5 r L + r ) L + r 4L 4 M + 98L Mr 7L r 3 + 4Mr 4 + r) 5 ) r M) ) L 4 M + 4L Mr + L r 3 + 3r 5, F θ[] =, 4.7) where ṙ F φ[] = F E πlr 4 L + r) 5/ φ[] E + Fφ[]K ) K, 4.73) Fφ[] E = E r 5 L 4 7L r + 3r) 4 L + r) 44L 6 M + 94L 4 Mr + 54L Mr 4 + L r 5 + 3r) 7, Fφ[] K = E r 7 3r L ) ) ) + r L + r L 4 M + 4L Mr + L r 3 + 3r 5, where F t[4] = 3Eṙ F E 4πr 11 L + r) 11/ t[4] E + Ft[4]K ) K, 4.74) 18

122 4.4 Schwarzschild Space-time Ft[4] E = 3E 4 r 14 3L 6 1L 4 r + 43L r 4 + r) 6 ) + E r 5 L + r 3456L 1 M L 1 Mr L 8 Mr L 6 Mr L 4 Mr 8 145L 4 r L Mr L r 11 ) + 3Mr r 13 4 L + r) 115L 14 M 184L 1 M r + 336L 1 Mr 3 664L 1 M r L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr 9 393L 4 M r L 4 Mr 11 3L 4 r 1 45L M r L Mr L r M r 14 ) + 9Mr r 16, Ft[4] K = 15E 4 r 16 87L L r + 5r) 4 ) 4E r 7 L + r 864L 1 M L 8 Mr L 6 Mr 4 ) L 4 Mr 6 5L 4 r L Mr 8 + 4L r Mr r 11 ) + r L + r 34L 1 M 5664L 1 M r + 67L 1 Mr 3 468L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr L M r L Mr 11 ) + 615L r M r Mr r 14, where F r[4] = 3 4πr 13 L + r ) 11/ r M) F E r[4] E + F K r[4]k ), 4.75) Fr[4] E = 3E 6 r 17 3L 6 1L 4 r + 43L r 4 + r) 6 ) E 4 r 8 L + r 691L 1 M L 1 Mr L 8 Mr L 6 Mr 6 3L 6 r L 4 Mr 8 745L 4 r L Mr 1 ) + 48L r Mr r 13 ) + E r 3 L + r 34L 14 M 3648L 1 M r + 67L 1 Mr L 1 M r L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr L 4 r L M r L Mr L r M r 14 ) + 68Mr r 16 + L + r) 3 M r ) 468L 14 M L 1 M r L 1 Mr L 1 M r L 1 Mr L 8 M r L 8 Mr 7 19

123 4.4 Schwarzschild Space-time 3688L 6 M r L 6 Mr L 4 M r L 4 Mr L 4 r 1 544L M r L Mr L r M r 14 ) 75r 16, Fr[4] K = 15E 6 r 19 87L L r + 5r) 4 ) L + r 3456L 1 M L 8 Mr L 6 Mr 4 + E 4 r L 4 Mr 6 915L 4 r L Mr L r Mr 1 ) + 585r 11 ) E r 5 L + r 34L 1 M 5664L 1 M r + 67L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr L 4 r L M r 1 ) L Mr L r M r Mr r 14 ) r L + r 3 r M) 6376L 8 M r 4 34L 1 M 4768L 1 M r 688L 1 Mr 3 84L 8 Mr L 6 M r 6 879L 6 Mr L 4 M r L 4 Mr L M r 1 354L Mr L r 1 ) 16M r r 14, F θ[4] =, 4.76) where F φ[4] = 3ṙ F E 4πLr 11 L + r) 9/ φ[4] E + Fφ[4]K ) K, 4.77) Fφ[4] E = 15E 4 r 14 L L 4 r 18L r 4 + 5r) 6 ) + E r 7 L + r 496L 1 M L 8 Mr L 6 Mr 4 ) L 4 Mr 6 165L 4 r L Mr 8 81L r r 11 + L + r) 468L 14 M L 1 M r L 1 Mr L 1 M r L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr L 4 r 1 544L M r L Mr L r M r 14 ) 75r 16, Fφ[4] K = 15E 4 r 16 L 4 58L r + 5r) 4 ) E r 9 L + r 48L 8 M + 63L 6 Mr + 679L 4 Mr 4 9L 4 r 5 ) + 356L Mr 6 375L r r 9 11

124 4.4 Schwarzschild Space-time ) + r L + r 34L 1 M 4768L 1 M r 688L 1 Mr L 8 M r 4 84L 8 Mr L 6 M r 6 879L 6 Mr L 4 M r L 4 Mr L M r 1 354L Mr L r 1 ) 16M r r Gravitational case Self-force regularization The self force on a gravitational particle is given by F a = k abcd hr bc;d, 4.78) where k abcd 1 gad u b u c g ab u c u d 1 ua u b u c u d ua g bc u d gad g bc. 4.79) Note that, as in the electromagnetic case, an ambiguity arises here due to the presence of terms involving the four-velocity at x. One is free to arbitrarily choose how to define this provided lim x x u a = uā. Following Barack and Sago 68), we choose to take the Schwarzschild components of the four velocity at x to be exactly those at x. The regularisation parameters in the gravitational case are given by r ṙ [-1] = L + r)r M), F [-1] r E = L + r), F θ [-1] =, F φ [-1] =, F t 4.8) where ṙ E [] = πr M)L + r) 3/ L E + rk), 4.81) F[] r 1 ) = F πrl 4 + r E[]E r + FK[]K r, 4.8) ) 3/ F t F r E[] = E L r 3 + r M)L + r )L + r ), F r K[] = r [ E r 3 + r M)L + r ) ], 111

125 4.4 Schwarzschild Space-time F θ [] =, F φ [] = ṙ πlr 3 L + r ) 1/ [L + r )E r K], 4.83) where F t [] = Eṙ πr 3 r M)L + r ) 7/ F t E[]E + D t K[]K), 4.84) F t E[] = E r 5 11L L r + 15r 4 ) L + r )76L 6 M + 768L 4 Mr + 78L Mr 4 37L r Mr 6 9r 7 ), F t K[] = E r 5 1L L r + 15r 4 ) + r L + r )93L 4 M + 4L Mr 9L r Mr 4 7r 5 ), where F r [] = 1 πr 7 L + r ) 7/ F r E[]E + F r K[]K), 4.85) F r E[] = E 4 r 8 11L L r + 15r 4 ) E r 3 L + r )138L 6 M + 4L 4 Mr 19L 4 r 3 + 4L Mr 4 34L r 5 + 1Mr 6 7r 7 ), + r M)L + r ) 188L 6 M + 46L 4 Mr + L Mr L r r 7 ), F r K[] = E 4 r 8 1L L r + 15r 4 ) + E r 5 L + r )1L 4 M 1L 4 r + 41L Mr 19L r Mr 4 + r 5 )] r r M)L + r ) 94L 4 M + 96L Mr + 14L r r 5 ), F θ [] =, 4.86) where F φ [] = ṙ πl 3 rl 6 + r) F φ 5/ E[] E + F φ K[] K), 4.87) 11

126 4.4 Schwarzschild Space-time F φ E[] = E L r 5 38L L r 15r 4 ), L + r )188L 8 M + 46L 6 Mr 64L 6 r 3 + L 4 Mr 4 163L 4 r 5 145L r 7 48r 9 ) F φ K[] = E L r 5 1L 4 + L r 15r 4 ) + r L + r )94L 6 M + 96L 4 Mr 74L 4 r 3 11L r 5 48r 7 ), where F t [4] = 3Eṙ 4πr 1 r M)L + r ) F E[4]E t + FK[4]K), t 4.88) 11/ FE[4] t = 3E 4 r 1 64L L 8 r + 989L 6 r 4 + 1L 4 r L r 8 + 1r 1 ) + E rl 5 + r)916l 1 M L 1 Mr L 8 Mr 4 19L 8 r L 6 Mr 6 978L 6 r L 4 Mr L 4 r L Mr 1 987L r Mr 1 585r 13 ) + 4L + r) 468L 14 M + 43L 1 M r 48L 1 Mr L 1 M r L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 6 r L 4 M r L 4 Mr L 4 r L M r L Mr L r M r Mr r 16 ), F t K[4] = 15E 4 r 1 64L L 6 r + 495L 4 r 4 + L r 6 81r 8 ) 4E r 7 L + r )59L 1 M + 137L 8 Mr L 6 Mr 4 48L 6 r L 4 Mr L 4 r L Mr 8 1L r 9 455Mr r 11 ) r L + r ) 916L 1 M L 1 M r 1356L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr L 4 r L M r 1 618L Mr L r M r Mr r 14 ), where F r [4] = 3 4πr 14 L + r) F E[4]E r + FK[4]K), r 4.89) 11/ 113

127 4.4 Schwarzschild Space-time F r E[4] = 3E 6 r 13 64L L 8 r + 989L 6 r 4 + 1L 4 r L r 8 + 1r 1 ) + E 4 r 8 L + r )1843L 1 M L 1 Mr 19L 1 r L 8 Mr 4 184L 8 r L 6 Mr L 6 r L 4 Mr L 4 r L Mr 1 163L r Mr 1 645r 13 ) + E r 3 L + r ) 916L 14 M L 1 M r 11L 1 Mr L 1 M r L 1 Mr L 8 M r L 8 Mr L 8 r L 6 M r L 6 Mr L 6 r L 4 M r L 4 Mr L 4 r L M r L Mr L r M r 14 88Mr 15 15r 16 ) r M)L + r ) L 14 M L 1 M r 4576L 1 Mr L 1 M r L 1 Mr L 8 M r L 8 Mr L 6 M r 8 716L 6 Mr L 6 r L 4 M r L 4 Mr L 4 r L M r 1 81L Mr 13 15L r M r 14 51Mr 15 55r 16 ), F r K[4] = 15E 6 r 15 64L L 6 r + 495L 4 r 4 + L r 6 81r 8 ) E 4 r 1 L + r )1368L 1 M L 8 Mr 88L 8 r L 6 Mr 4 168L 6 r L 4 Mr L 4 r L Mr L r 9 87Mr r 11 ) E r 5 L + r ) 916L 1 M L 1 M r 578L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 6 r L 4 M r L 4 Mr 9 9L 4 r L M r 1 568L Mr L r M r Mr r 14 ) + r r M)L + r ) 3 916L 1 M L 1 M r 15744L 1 Mr L 8 M r L 8 Mr L 6 M r L 6 Mr L 4 M r L 4 Mr 9 96L 4 r L M r 1 388L Mr L r M r 1 51Mr 13 55r 14 ), F θ [4] =, 4.9) F φ [4] = 3ṙ 4πr 13 L 5 L + r ) F φ 9/ E[4] E + F φ K[4] K), 4.91) 114

128 4.4 Schwarzschild Space-time where F φ E[4] = 15E4 L 4 r 1 18L L 8 r + 66L 6 r L 4 r 6 8L r 8 35r 1 ) + E L rl 5 + r)115l 14 M L 1 Mr L 1 Mr 4 576L 1 r L 8 Mr 6 13L 8 r L 6 Mr L 6 r L 4 Mr L 4 r L r r 15 ) + L + r) 1843L 18 M L 16 M r 4576L 16 Mr L 14 M r L 14 Mr L 1 M r 6 197L 1 Mr L 1 M r L 1 Mr 9 576L 1 r L 8 M r L 8 Mr L 8 r L 6 M r L 6 Mr L 6 r L 4 M r 14 63L 4 Mr L 4 r 16 11L Mr L r 18 56r ), F φ K[4] = 15E4 L 4 r 1 19L L 6 r + 169L 4 r 4 3L r 6 35r 8 ) E L r 7 L + r )6336L 1 M L 1 Mr L 8 Mr 4 48L 8 r L 6 Mr 6 711L 6 r L 4 Mr L 4 r L r r 13 ) r L + r ) 916L 16 M L 14 M r 15744L 14 Mr L 1 M r 4 48L 1 Mr L 1 M r 6 391L 1 Mr L 8 M r L 8 Mr L 8 r L 6 M r L 6 Mr L 6 r L 4 M r 1 564L 4 Mr L 4 r 14 11L Mr L r 16 56r 18 ) huu regularization The quantity H R) = 1 hr) ab ua u b 4.9) was first proposed by Detweiler 45) as a tool for constructing gauge invariant measurements from self-force calculations. It has since proven invaluable in extracting gauge invariant results from gauge dependent self-force calculations 93, 1). Much the same as with self-force calculations, the calculation of H R) requires the subtraction of the appropriate singular piece, H S) = 1 hs) ab ua u b from the full retarded field. In this section, we give this subtraction in the form of mode-sum regularization parameters. In doing so, we keep with our convention that the term 115

129 4.4 Schwarzschild Space-time proportional to l + 1 is denoted by H [-1] = in this case), the constant term is denoted by H [], and so on. Note that, as in the self-force case, an ambiguity arises here due to the presence of terms involving the four-velocity at x. One is free to arbitrarily choose how to define this provided lim x x u a = uā. As before, we choose this in such a way that the Schwarzschild components of the four velocity at x are exactly those at x. The regularisation parameters are then given by where H [] = K π L + r, 4.93) H [] = HE [] E + HK [] K, 4.94) πrl 3 + r) 3/ H E [] = E r 5 + L + r )36L M 8L r + 38M rb 9r 3 ), H K [] = E r 3 16L + 17r ) L + r )16L M 4L r + 33Mr 1r 3 ), where H [4] = 3HE [4] E + HK [4] K), 4.95) πr 1 L + r ) 7/ H E [4] = 1E 4 r 1 8L L r + 7r 4 ) + E r 5 L + r )3584L 8 M + 171L 6 Mr L 4 Mr 4 + 1L 4 r L Mr L r Mr r 9 ) + L + r ) 1536L 1 M L 8 M r 16L 8 Mr L 6 M r 4 944L 6 Mr L 4 M r 6 141L 4 Mr 7 + 1L 4 r L M r 8 535L Mr L r M r Mr 11 9r 1 ), H[4] K = 15E 4 r 1 64L 6 + 4L 4 r + 59L r r) 6 4E rl 7 + r)1376l 6 M L 4 Mr + 4L 4 r L Mr L r 5 + 7Mr r) 7 r L + r ) 1536L 8 M L 6 M r 736L 6 Mr L 4 M r 4 193L 4 Mr L M r 6 114L Mr 7 7L r M r Mr 9 75r 1 ). 116

130 4.4 Schwarzschild Space-time Example - Scalar Self-Force As an example application of our high order regularization parameters, we consider the case of a scalar particle on a circular geodesic of the Schwarzschild space-time. In this case, the retarded field may be computed using the frequency domain method described in 66), along with improved asymptotics for the boundary conditions by expanding inside the exponential rather than outside) and with the use of the arbitrary precision differential equation solving support in Mathematica 77). These improvements allowed us to substantially increase the accuracy of the computed retarded field. We found this to be necessary to get the full benefit from the the higher order regularization parameters. If we consider the the scalar wave equation in Eq..16) with zero Ricci scalar as is the case for Schwarzschild space-time, Φx) = 4πµx), 4.96) where the distributional source, µx) = q γ δ 4 x x ) g dτ, 4.97) is representing a point charge, q moving along a world line γ described by z a τ), where τ is proper time. If we consider a circular orbit, i.e., x = {tτ), r, π/, Ωtτ)} with Ω = M/r 3 ) 1/, we can rewrite the distributional source accordingly, µ = q g) 1 δ t tτ)) δ r r ) δ θ π/) δ φ Ωtτ)) dτ γ ) 1 dt = qr δ r r ) δ θ π/) δ φ Ωt) 4.98) dτ where we have used dτ = dt dτ ) 1 dt. It is beneficial to now decompose δ θ π/) δ φ Ωt) into spherical harmonics, δ θ π/) δ φ Ωt) = lm D lm Y lm θ, φ), 4.99) where 117

131 4.4 Schwarzschild Space-time D lm = = l + 1) l m)! π π/ e imφ δ φ Ωt) δ θ π ) P lm cos θ) sin θdθdφ 4π l + m)! π π/ l + 1) l m)! 4π l + m)! e imωt P lm cos π/) = e imωt Y lm π/, ). 4.1) Incorporating this into Eq. 4.98) gives, µ = lm q lm 4πr δ r r ) e iωmt Y lm θ, φ), 4.11) where and dt dτ ω m = mω, q lm = 4πq r Y lm can be drived from g abẋ a ẋ b = 1 π/, ), dt dτ dt dτ = r r 3M, 4.1) Expanding out Φ = g ab a b Φ in Schwarzschild space-time gives r r M) Φ = Φ r M) Φ + r r M) r r Φ r M) t [ 1 Φ + r r M) θ + cos θ Φ sin θ θ + 1 ] Φ sin, 4.13) θ φ while decomposing the retarded field gives, Φ = lm Φ lm e iωmt Y lm θ, φ). 4.14) Using Eqs. 4.11), 4.13) and 4.14) in Eq. 4.96) allows us to write the lm part of the scalar wave equation as, [ ] d Φ lm r M) dφ lm ω r l l + 1) + + dr Φ lm = q lm r M) dr r M) r r M) r M δ r r ). 4.15) The tortoise coordinate was first introduced by Wheeler 13), it is designed to remove the single derivative, d/dr and is given by, d = 1 M ) d r ) dr r dr, or r = r + M log M ) 118

132 4.4 Schwarzschild Space-time To obtain appropriate boundary conditions, we consider ingoing waves at the horizon, Φ lm = eiωr, r for r M, 4.17) and outgoing waves at infinity, Φ lm = e iωr, for r. 4.18) r Our aim is now to solve the inhomogeneous equation of Eq. 4.15) using the above boundary conditions. To agree with these conditions, we assume that Φ lm admits an asymptotic expansion in 1/r at r and an asymptotic expansion in r M) as r M. We know what the leading order behaviour should look like from Eqs. 4.17) and 4.18), keeping this in mind, we assume expansions of Φ lm for Φ in) lm and Φout) lm to be, Φ out) lm r) = exp iωr r n= a n ) r n, and Φ in) e iωr lm r) = r b n r M) n n= 4.19) To determine the coefficients a n and b n, we use Eq. 4.15). For the outgoing waves, hence the a n s, we can use mathematical packages like Mathematica to substitute Φ out) lm into the homogeneous equation of Eq. 4.15) and solve for the coefficients using initial conditions, a = 1, a n< =. For the ingoing wave, it is possible to analytically obtain a recursion relation for the b n s by substituting in Φ in) lm into Eq. 4.15), that is, b n = 1iωM n 1) + n 3) n 1) l + l + 1) b M 4inωM + n n 1 ) 1iωM n ) + n ) n 3) l l + 1) b 4M 4inωM + n n ) iω n 3) + M 4inωM + n ) b n 3, 4.11) and use starting values, b = 1 and b n< = to determine the required b n s. Once values for the coefficients a n and b n have been obtained, we use our expressions for Φ in) lm and Φout) lm from Eqs. 4.19) in Eq. 4.15) to derive initial conitions away from the singular parts so we can numerically integrate to obtain homogeneous 119

133 4.4 Schwarzschild Space-time solutions, with the appropriate boundary conditions. The inhomogeneous solution to Eq. 4.15) is now of the type, Φ lm = A lm Φ in) lm r)θ r r) + B lm Φ out) lm r)θ r r ) = Φ lm r)θ r r) + Φ + lm r)θ r r ), 4.111) where Θ r r) is the Heaviside step function previously introduced in Eq..134) and we have introduced the notation Φ lm r) = A lmφ in) lm r) and Φ+ lm r) = B lmφ out) lm r). A lm and B lm can be determined by imposing suitable matching conditions. These are found by substituting Φ lm from Eq ) into the radial equation, Eq. 4.15), and gathering like terms; where = δ r r ) Φ + lm ) Φ lm [ + δ r r ) Φ + lm Φ ) ] r M) lm + Φ + lm r M ) q lm Φ lm + r M [ ] + Θ r r) Φ r M) lm + r M Φ lm + W Φ lm [ ] + Θ r r ) Φ + r M) lm + r M Φ+ lm + W Φ + lm, 4.11) W = and refers to differentiation with respect to r. ω r l l + 1) r M) r r M), 4.113) It can be clearly seen that the last two terms of Eq. 4.11) are zero as they are the homogeneous radial equation, which Φ lm r) and Φ+ lm r) solve by design. Our other two terms give, [ Φ + lm ] Φ lm r=r =, [ Φ + lm Φ ] lm r=r = 1 [ ] qlm ) r M r=r It is now possible to substitute Φ lm r) = A lmφ in) lm r) and Φ+ lm r) = B lmφ out) lm r) back into Eq ), to obtain a simple system of simultaneous equations in A lm and B lm. These are easily solved to give, 1

134 4.4 Schwarzschild Space-time fr Figure 4.1: Regularization of the radial component of the scalar self-force in Schwarzschild space-time for the case of a scalar particle on a circular geodesic of radius r = 1M in Schwarzschild space-time. In decreasing slope the above lines represent the unregularised self-force, self-force regularised by subtracting from it in turn the cumulative sum of Fr[ 1] l, F r[] l, F r[] l, F r[4] l, F r[6] l, F r[8] l, F r[1] l, F r[1] l. A lm = B lm = [ r M) [ r M) Φ out) lm Φ out) lm q lm Φ in) lm r ) ], r )Φ in) lm r ) Φ in) lm r )Φ out) lm r ) q lm Φ out) lm r ) ] ) r )Φ in) lm r ) Φ in) lm r )Φ out) lm r ) by The l component of the retarded self-force with the help of Eq. 4.5) is now given F l aret) = x a l m= l [ ] A lm Φ in) lm r)θ r r) + B lm Φ out) lm r)θ r r ) 11

135 4.4 Schwarzschild Space-time where we have numerically solved for Φ in lm r) and Φout lm r) and A lm and B lm are given by Eq ). We have also replaced p aa with its scalar operator. The l mode x a of the full self force can now be calculated from Eq ), F a = l F l aret) F l as)), 4.116) where we have shown the calculation of both F l aret) and F l as) in the last two sections Impact of Regularisation Parameters in Schwarzschild Space-Time The results and benefits) of the calculations in Sec. 4.4 are illustrated in Figs. 4.1, 4. and 4.3. There we show the effect of subtracting in turn the cumulative sums of the regularization parameters from the full retarded field. In Fig. 4.1, in order from top to bottom are Fr ret and the result of subtracting from it in turn the cumulative sum of the regularization terms Fr[-1] l, F r[] l, F r[] l, F l r[4], F l r[6], F l r[8], F l r[1] and F l r[1]. The parameters F r[-1], F r[], F r[], F r[4] and F r[6] are analytically derived in Sec , while F r[8], F r[1] and F r[1] were determined through a numerical fit to the data. The resulting rapid convergence with l enables the calculation of an extremely accurate value for the self-force. Summing over l, we find F r = ), where the uncertainty in the last digit is estimated by assuming that the only error comes from limiting the sum to a finite l max = 8. In addition to providing a highly accurate benchmark, the example in the previous section may be used to assess the benefits which can be obtained from the use of higher-order regularization parameters. The most obvious benefit is that with fixed computational resources i.e. fixed number of spherical harmonic modes) one can obtain a much more accurate value for the self-force. This is highlighted by comparison of our value for F r with that of the previous benchmark given in 66), F r = ). Both calculations consider the same case of a scalar charge in a circular orbit of radius 1M around a Schwarzschild black hole. Using 4 l-modes and regularization parameters up to Fr[] l, 66) obtained a value for the 1

136 4.4 Schwarzschild Space-time l max = 5, n = 1 l max = 8, n = 5 RPs used abs. rel. abs. rel. AB ABD ABDF ABDF H CPU time 155s 447s Table 4.: Table demonstrating the usefulness of the analytically derived higher order regularization parameters in practical self-force calculations. In this example we show the regularization of the radial component of the scalar self-force for a circular orbit at r = 1M about a Schwarzschild black hole. In the left most column we list the analytically derived regularization parameters employed in each calculation. For each calculation we numerically fit the higher order regularization parameters up to F r[15]. The next two wide columns show the result of computing the scalar self-force for 5 and 8 l-modes respectively l max = number of modes) and numerically fitting the unknown regularization parameters. We show the resulting absolute value of the selfforce and its relative difference verses the highly accurate value provided in the main text. The CPU time taken to compute the l-modes of the retarded field using a code running on 1 cores of a machine with a 3GHz clock speed is also given. This shows the improvement in run-time is over a factor of 36 self-force with a fractional accuracy of 1 9. The inclusion of the next two regularization parameters improves this to a fractional error of 1 1 which increases to a fractional error of 1 17 when 8 modes are used. To illustrate this further we have included Table 4.4.5, where the physical time taken to run numerical code is recorded as is the accuracy obtained. The results in this table show that, by using the regularization parameters derived in this thesis, it is possible to calculate a scalar self-force using only 5 l-modes which is as accurate as a calculation made using 8 l-modes and regularizing with only the A and B parameters. As the high l-modes are computationally expensive to calculate, 13

137 4.4 Schwarzschild Space-time 1 Electromagnetic self force F l r Figure 4.: Regularization of the radial component of the self-force for the case of a electromagnetic particle on a circular geodesic of radius r = 1M in Schwarzschild space-time. In decreasing slope the above lines represent the unregularised self-force and the self-force regularised by subtracting from it in turn the cumulative sum of Fr[ 1] l, F r[] l, F r[] l, F r[4] l, F r[6] l, F r[8] l, F r[1] l. using the higher order regularization parameters offers a substantial improvement in code run-time for a fixed level of accuracy in the final result. In this example the improvement in run-time is over a factor of 36. This example represents a somewhat extreme case: it uses highly accurate frequency domain methods combined with high-precision numerical integration and a relatively large number of spherical harmonic modes. In more typical time-domain calculations, numerical data up to l 15 is used and it is common that the dominant source of error comes from the tail fit. While it may seem that one merely needs to compute more modes to reduce this error, this is not a realistic solution. In a mode-sum calculation, the number of spherical harmonic modes required for each 14

138 4.4 Schwarzschild Space-time Regularization of the gravitational self force.1.1 F r Figure 4.3: Regularization of the radial component of the self-force for the case of a gravitational particle on a elliptic geodesic of radius r = 1M in Schwarzschild space-time. The graph is plotting Fr l against l. In decreasing slope the above lines represent the unregularised self-force and the self-force regularised by subtracting from it in turn the cumulative sum of Fr[ 1] l, F r[] l, F r[] l, F r[4] l. l scales as l, meaning that simply running simulations for larger and larger l rapidly becomes prohibitively expensive in terms of computational cost. Additionally, the improvement with each additional mode falls off as an inverse power in l, meaning that many more l modes are required for an increasingly small benefit. In this case, the inclusion of higher order regularization parameters essentially eliminates this problem: without them the tail fit is the dominant source of error, with them sufficiently accurate results may be obtained without even fitting for a tail. In the electromagnetic and gravitational cases, data for the retarded field, the ansatz of what was calculated in Sec , is increasingly more complicated with higher spins. As outlined in Sec , the aim of this thesis is to concentrate on 15

139 4.5 Kerr Space-time the singular component of the self-force. We therefore did not derive the retarded field for the higher spin cases. However, with the assistance of Roland Haas, Patrick Nolan and the Southampton group Sarp Ackay, Niels Warburton, Leor Barack), who gave us access to their data for the retarded field in these cases, we are able to show the success of our regularisation parameters. Figs. 4. and 4.3 depict the retarded field in the electromagnetic and gravitational cases respectively; unregularised and then regularised by the sum of the parameters for different n, where n can be considered as the order of ɛ to which we calculated the singular field - this is currently being increased with ongoing work from the Dublin Self-Force group 14). In all three cases, we can clearly see that with the regularisation parameters comes a jump in fractional accuracy with the same number of l-modes - meaning it is possible to calculate more accurate data with the same number of l modes. The importance of this lies with the exponential increase in computation time with the higher l modes, making the higher regularisation parameters invaluable to those numerically calculating the self-force. It should be pointed out that there is one caveat to our conclusions. The use of high order regularization parameters requires the subtraction of increasingly relatively) large numbers to obtain a small regularized remainder. It is therefore essential that any numerically provided data for the retarded field must be of sufficient accuracy for the subtraction to yield meaningful results. As a result, calculations which were previously deemed sufficient would not necessarily gain an immediate benefit from higher order regularization parameters. 4.5 Kerr Space-time The Kerr cases follow the same necessary calculations, for the self-force, as their Schwarzschild counterparts. I therefore won t go through the calculations again but remind the reader that the necessary calculations for the scalar, electromagnetic, gravitational and huu regularisation parameters is explained in detail in Secs , 4.4., and respectively. The regularisation parameters for Kerr at the 16

140 4.5 Kerr Space-time higher orders prove to be too large for paper format, in these instances they have been made available in electronic form 7) Scalar Case The regularisation parameters in the Kerr scalar case are given by, r ṙ sgn r F t[-1] =, r a + L ) + a M + r 3 F r[-1] = sgn r Er a + r) + amae L)) a Mr + r) r a + L ) + a M + r), 3 F θ[-1] =, F φ[-1] =, 4.117) where F t[] = πr ṙ r + L + a M r + a ) 3/ F E t[] E + F K t[]k ), 4.118) Ft[] E = 4aLM ) 4a 4 M + a 4 Mr + a L Mr a Mr 3 a r 4 L r 4 + E 1a 6 M 3 16a 6 M r 7a 6 Mr a 6 r 3 4a 4 L M r 6a 4 L Mr a 4 L r 3 6a 4 M r 3 5a 4 Mr 4 a 4 r 5 + a L 4 Mr a L 4 r 3 5a L Mr 4 3a L r 5 L 4 r) 5, F K t[] = alm a 4 M a 4 Mr a 4 r a L Mr a L r a Mr 3 a r 4 ) L 4 r L r 4 + E 4a 6 M 3 + 4a 6 M r + a 6 Mr a 4 L M r a 4 L Mr + a 4 M r 3 + a 4 Mr 4 a L 4 Mr + a L Mr 4 + a L r 5 + L 4 r 5 ), F r[] = where Fr[] E E + F r[] K K ) 3/, πr 3 a M + a r + L r ) r + L + a M r + a r Mr + a ) 4.119) 17

141 4.5 Kerr Space-time F E r[] = 4a 8 M 3 r 3a 8 M r 14a 8 Mr 3 a 8 r 4 + 4a 6 L M 4 + 1a 6 L M 3 r 3a 6 L M r 7a 6 L Mr 3 6a 6 L r a 6 M 4 r + 4a 6 M 3 r 3 16a 6 M r 4 a 6 Mr 5 4a 6 r 6 + 8a 4 L 4 M 3 r 1a 4 L 4 Mr 3 6a 4 L 4 r a 4 L M 3 r 3 + 1a 4 L M r 4 1a 4 L Mr 5 9a 4 L r 6 + 4a 4 M 3 r 5 + 8a 4 M r 6 6a 4 Mr 7 a 4 r 8 a L 6 M r + a L 6 Mr 3 a L 6 r 4 + 6a L 4 M r 4 + a L 4 Mr 5 6a L 4 r 6 + 1a L M r 6 3a L r 8 ) + L 6 Mr 5 L 6 r 6 + L 4 Mr 7 L 4 r 8 aelm 4a 6 M 3 + 8a 6 M r + 1a 6 Mr + a 6 r 3 + 8a 4 L M r + 8a 4 L Mr + a 4 L r 3 4a 4 Mr 4 a 4 r 5 a L 4 Mr + a L 4 r 3 ) a L Mr 4 a L r 5 6a Mr 6 3a r 7 + L 4 r 5 3L r 7 + E a M + a r + r 3 ) 1a 6 M a 6 M r + 7a 6 Mr + a 6 r 3 + 4a 4 L M r + 6a 4 L Mr + a 4 L r 3 + 6a 4 M r 3 + 5a 4 Mr 4 + a 4 r 5 a L 4 Mr + a L 4 r 3 + 5a L Mr 4 + 3a L r 5 + L 4 r 5 ), F K r[] = 8a 8 M 3 r + 8a 8 M r + a 8 Mr 3 8a 6 L M 4 + 4a 6 L M 3 r + 1a 6 L M r + 4a 6 L Mr 3 16a 6 M 4 r 8a 6 M 3 r 3 + 8a 6 M r 4 + 4a 6 Mr 5 + 4a 4 L 4 M 3 r + 8a 4 L 4 M r + a 4 L 4 Mr 3 + 8a 4 L M 3 r 3 + 1a 4 L M r 4 + a 4 L Mr 5 a 4 L r 6 8a 4 M 3 r 5 + a 4 Mr 7 + 4a L 6 M r + 16a L 4 M r 4 a L 4 r 6 ) + 4a L M r 6 a L r 8 + L 6 Mr 5 L 6 r 6 + L 4 Mr 7 L 4 r 8 + aelm 8a 6 M 3 + 4a 6 M r a 6 Mr a 6 r 3 4a 4 L M r 6a 4 L Mr a 4 L r 3 8a 4 M r 3 1a 4 Mr 4 4a 4 r 5 4a L 4 Mr a L 4 r 3 ) 1a L Mr 4 5a L r 5 6a Mr 6 3a r 7 L 4 r 5 3L r 7 E a M + a r + r 3 ) 4a 6 M 3 + 4a 6 M r + a 6 Mr a 4 L M r a 4 L Mr + a 4 M r 3 + a 4 Mr 4 a L 4 Mr + a L Mr 4 + a L r 5 + L 4 r 5 ), F θ[] =, 4.1) F φ[] = where Lṙ πr a M + a r + L r ) r + L + a M r + a ) 1/ F E φ[] E + F K φ[]k ), 4.11) 18

142 4.5 Kerr Space-time Fφ[] E = a 4 M a 4 Mr a L Mr 4a Mr 3 a r 4 L r, 4 ) Fφ[] K = r 3 4a M + a r + L r. The regularisation parameters for F a[] and F a[4] are too large for paper format and have instead been made available electronically 7). For the reader to get an understanding of the form and size of these expressions, we have included F t[] for eccentric orbits, and F r[] for a circular orbit. Expressions online are for eccentric orbits all expressions in Kerr space-time are for the case equatorial plane). condense F t[], the resulting expression we have used the notation, F t[] is described by, L = L + a + r + Ma r. 4.1) To F t[] = where ṙ 6πr 7 a M + a r + L r ) 5 r + L + a M r + a ) 7/ F E t[] E + F K t[]k ), 4.13) Ft[] E = [ 1aLMM r )rl aLMr 4 4r Mr 3 + 5a r + a M ) L 16 almr 3 399r 7 984Mr 6 1a r 5 8a Mr 4 + 4a M r 3 138a 4 Mr + 17a 4 M r + 168a 4 M 3) L 14 + almr 73r 1 87Mr 9 181a r 8 73a Mr a M r 6 36a 4 Mr a 4 M r a 4 M 3 r 3 15a 6 M r + 84a 6 M 3 r + 676a 6 M 4) L 1 almr 18r 13 31Mr 1 15a r 11 9a Mr a M r 9 1a 4 Mr a 4 M r a 4 M 3 r a 6 M r a 6 M 3 r a 6 M 4 r 3 6a 8 M 3 r + 1a 8 M 4 r + 46a 8 M 5) L 1 + alm 138r Mr a r 14 + a Mr a M r 1 + 9a 4 Mr a 4 M r a 4 M 3 r a 6 M r a 6 M 3 r 7 978a 6 M 4 r a 8 M 3 r 5 81a 8 M 4 r a 8 M 5 r 3 + 4a 1 M 5 r + 84a 1 M 6) L 8 + almr 165r 16 31Mr a r 14 3a Mr a M r 1 339a 4 Mr a 4 M r a 4 M 3 r a 6 M r a 6 M 3 r a 6 M 4 r 6 a 8 M 3 r 5 19

143 4.5 Kerr Space-time + 98a 8 M 4 r a 8 M 5 r a 1 M 5 r a 1 M 6) L 6 almr 4 99r Mr a r 14 53a Mr a M r 1 59a 4 Mr a 4 M r a 4 M 3 r a 6 M r a 6 M 3 r 7 1a 6 M 4 r 6 117a 8 M 3 r 5 588a 8 M 4 r 4 44a 8 M 5 r a 1 M 5 r + 396a 1 M 6) L aLMr 6 a M r) 3 3r Mr 1 3a r a Mr a M r 9 + 6a 4 Mr a 4 M r 7 18a 4 M 3 r 6 3a 6 M r 5 16a 6 M 3 r 4 18a 6 M 4 r 3 + 1a 8 M 4 r + 4a 8 M 5) L 384a 3 LM 3 r 8 a M r) 3 3 r 3 + a r + a M ) ] + E [ 1Mr 5 9r 3 + 4a M ) L Mr 4 9r a r 4 + 4a Mr 3 + 4a 4 Mr + 8a 4 M ) L 16 + r 3 3r Mr 9 4a Mr 7 51a M r 6 6a 4 Mr a 4 M r a 4 M 3 r 3 9a 6 Mr a 6 M r + 34a 6 M 3 r a 6 M 4) L 14 r 39r 13 11Mr 1 4a r a Mr 1 9a 4 r 9 538a M r 9 19a 4 Mr a 4 M r a 4 M 3 r 6 198a 6 Mr a 6 M r a 6 M 3 r a 6 M 4 r 3 16a 8 M r a 8 M 3 r + 114a 8 M 4 r a 8 M 5) L 1 + r 15r 16 36Mr 15 16a r 14 43a Mr 13 66a 4 r a M r 1 13a 4 Mr a 4 M r a 4 M 3 r 9 63a 6 Mr a 6 M r 8 + 4a 6 M 3 r a 6 M 4 r a 8 M r a 8 M 3 r a 8 M 4 r a 8 M 5 r 3 6a 1 M 3 r a 1 M 4 r + 178a 1 M 5 r + 145a 1 M 6) L 1 + 7r 19 37Mr a r 17 19a Mr a 4 r a M r 15 5a 4 Mr 14 48a 4 M r 13 71a 4 M 3 r 1 43a 6 Mr a 6 M r a 6 M 3 r a 6 M 4 r 9 159a 8 M r 9 399a 8 M 3 r 8 615a 8 M 4 r a 8 M 5 r 6 453a 1 M 3 r a 1 M 4 r a 1 M 5 r a 1 M 6 r 3 16a 1 M 5 r 54a 1 M 6 r 54a 1 M 7) L 8 55r Mr a r a Mr a 4 r 15 r 714a M r 15 58a 4 Mr a 4 M r a 4 M 3 r 1 45a 6 Mr a 6 M r a 6 M 3 r a 6 M 4 r a 8 M r a 8 M 3 r a 8 M 4 r a 8 M 5 r 6 a 1 M 3 r a 1 M 4 r a 1 M 5 r a 1 M 6 r 3 13r a 1 M 5 r a 1 M 6 r a 1 M 7) L 6 + r 4 46Mr a r a Mr a 4 r a M r 15 37a 4 Mr a 4 M r a 4 M 3 r 1 79a 6 Mr 1 13

144 4.5 Kerr Space-time 33a 6 M r a 6 M 3 r a 6 M 4 r a 8 M r a 8 M 3 r a 8 M 4 r 7 1a 8 M 5 r 6 117a 1 M 3 r a 1 M 4 r a 1 M 5 r a 1 M 6 r a 1 M 5 r + 376a 1 M 6 r + 376a 1 M 7) ) L 4 8r 6 a M r 3 r 3 + a r + a M ) 3r Mr 1 3a r a Mr a M r 9 + 6a 4 Mr a 4 M r 7 1a 4 M 3 r 6 3a 6 M r 5 5a 6 M 3 r 4 55a 6 M 4 r a 8 M 4 r + 7a 8 M 5) L + 576a M r 8 a M r) 3 3 r 3 + a r + a M ) ] + E [ 4aLMr 5 9r 3 + 5a r + 1a M ) L aLMr 4 77r 6 3a r 4 49a Mr 3 + 5a 4 r + a 4 Mr + a 4 M ) L 14 almr 3 73r 9 17a r 7 76a Mr a 4 r a 4 Mr a 4 M r a 6 Mr + 183a 6 M r + 17a 6 M 3) L 1 + almr 348r a r 1 54a Mr 9 + a 4 r a 4 Mr a 4 M r a 6 Mr a 6 M r a 6 M 3 r a 8 M r + 7a 8 M 3 r + 94a 8 M 4) L 1 almr 61r a r a Mr a 4 r a 4 Mr 1 + 5a 4 M r a 6 Mr a 6 M r a 6 M 3 r 6 179a 8 M r 5 179a 8 M 3 r 4 364a 8 M 4 r a 1 M 3 r + 16a 1 M 4 r + 16a 1 M 5) L 8 almr 3 141r 15 7a r a Mr 1 674a 4 r a 4 Mr a 4 M r a 6 Mr a 6 M r a 6 M 3 r 6 9a 8 M r a 8 M 3 r a 8 M 4 r a 1 M 3 r a 1 M 4 r a 1 M 5) L 6 + almr 5 138r a r a Mr 1 963a 4 r a 4 Mr a 4 M r a 6 Mr a 6 M r 7 1a 6 M 3 r 6 351a 8 M r 5 613a 8 M 3 r 4 6a 8 M 4 r a 1 M 3 r + 594a 1 M 4 r + 594a 1 M 5) L 4 16aLMr 7 a M r) 3 r 3 + a r + a M ) 87r a r 7 13a 4 Mr 4 159a 4 M r 3 + 9a 6 M r + 18a 6 M 3) ) L + 88aLMr 9 a M r 3 3 r 3 + a r + a M ) ] + E 3[ 6a MrL a Mr 7 44r 3 + 9a r + a M ) L r 3 8r 1 + 1a r 1 131a Mr a 4 Mr 7 a 6 r 6 + 8a 4 M r 6 + 7a 6 Mr a 6 M r a 6 M 3 r 3 + 3a 8 Mr 3 + 8a 8 M r + 56a 8 M 3 r + 448a 8 M 4) L 1 3r 56r a r 13 + a Mr a 4 r a 4 Mr 1 a 6 r 9 + 6a 4 M r a 6 Mr a 6 M r 7 131

145 4.5 Kerr Space-time + 911a 6 M 3 r a 8 Mr a 8 M r a 8 M 3 r a 8 M 4 r 3 + 6a 1 M r 3 + 4a 1 M 3 r + 88a 1 M 4 r + 64a 1 M 5) L 1 + 3r 16r a r a Mr 15 + a 4 r a 4 Mr a 6 r a 4 M r a 6 Mr a 6 M r 1 163a 6 M 3 r a 8 Mr a 8 M r a 8 M 3 r 7 a 8 M 4 r 6 + 4a 1 M r 6 173a 1 M 3 r 5 74a 1 M 4 r 4 756a 1 M 5 r a 1 M 3 r 3 + 1a 1 M 4 r + 4a 1 M 5 r + 8a 1 M 6) L 8 + 3r 3 r 3 + a r + a M ) 4r 15 38a r 13 17a Mr 1 98a 4 r 11 + a 4 Mr a 4 M r a 6 Mr a 6 M r 7 + 3a 6 M 3 r 6 76a 8 M r 5 + 3a 8 M 3 r a 8 M 4 r a 1 M 3 r + 776a 1 M 4 r + 776a 1 M 5) L 6 3r 5 r 3 + a r + a M ) r 15 3a r a Mr 1 14a 4 r a 4 Mr a 4 M r a 6 Mr a 6 M r 7 19a 6 M 3 r 6 84a 8 M r 5 78a 8 M 3 r 4 3a 8 M 4 r a 1 M 3 r + 66a 1 M 4 r + 66a 1 M 5) ) L 4 + 4r 7 a M r 3 r 3 + a r + a M ) 11r 9 + 9a r 7 a Mr 6 15a 4 Mr 4 18a 4 M r 3 + 1a 6 M r + a 6 M 3) L ) 48r 9 a M r 3 3 r 3 + a r + a M ) 3 ], Ft[] K = [ 6aLMM r )rl aLMr 6 1r 4 + Mr 3 + a r + a Mr + a M ) L 14 + almr 5 18r 7 456Mr 6 99a r 5 116a Mr 4 138a M r 3 34a 4 Mr 9a 4 M r + 894a 4 M 3) L 1 almr 4 13r 1 48Mr 9 14a r 8 9a Mr a M r 6 535a 4 Mr 5 36a 4 M r a 4 M 3 r 3 495a 6 M r + 855a 6 M 3 r + 598a 6 M 4) L 1 + almr 3 78r 13 5Mr 1 54a r 11 33a Mr a M r 9 41a 4 Mr a 4 M r a 4 M 3 r 6 353a 6 M r a 6 M 3 r a 6 M 4 r 3 345a 8 M 3 r + 8a 8 M 4 r + 816a 8 M 5) L 8 almr 7r Mr a r 14 3a Mr a M r 1 9a 4 Mr a 4 M r a 4 M 3 r a 6 M r a 6 M 3 r a 6 M 4 r 6 6a 8 M 3 r a 8 M 4 r a 8 M 5 r 3 + 1a 1 M 5 r + 4a 1 M 6) L 6 + almr 4 48r 16 96Mr a r 14 1a Mr a M r 1 15a 4 Mr 11 89a 4 M r a 4 M 3 r a 6 M r a 6 M 3 r 7 114a 6 M 4 r 6 57a 8 M 3 r 5 117a 8 M 4 r 4 13a 8 M 5 r a 1 M 5 r + 187a 1 M 6) L 4 8aLMr 6 a M ) r 3 3r Mr 1 3a r a Mr a M r 9 + 6a 4 Mr 8 13

146 4.5 Kerr Space-time + 8a 4 M r 7 18a 4 M 3 r 6 3a 6 M r 5 154a 6 M 3 r 4 171a 6 M 4 r a 8 M 4 r + 34a 8 M 5) ) L + 19a 3 LM 3 r 8 a M r 3 3 r 3 + a r + a M )] + E [ 6Mr 7 9r 3 + 4a M ) L 16 6Mr 6 45r a r a Mr 3 + 3a 4 r + 76a 4 Mr + a 4 M ) L 14 r 5 1r 1 564Mr a r 8 73a Mr 7 151a M r 6 1a 4 Mr 5 613a 4 M r a 4 M 3 r 3 16a 6 Mr 3 678a 6 M r + 848a 6 M 3 r a 6 M 4) L 1 + r 4 6r 13 66Mr a r a Mr a 4 r a M r 9 88a 4 Mr a 4 M r a 4 M 3 r 6 333a 6 Mr 6 873a 6 M r a 6 M 3 r a 6 M 4 r 3 378a 8 M r a 8 M 3 r + 186a 8 M 4 r a 8 M 5) L 1 r 3 1r 16 51Mr a r 14 46a Mr a 4 r 1 34a M r 1 37a 4 Mr a 4 M r a 4 M 3 r 9 351a 6 Mr a 6 M r a 6 M 3 r a 6 M 4 r 6 58a 8 M r a 8 M 3 r a 8 M 4 r a 8 M 5 r 3 345a 1 M 3 r a 1 M 4 r a 1 M 5 r a 1 M 6) L 8 + r 1r 19 94Mr a r a Mr a 4 r a M r a 4 Mr a 4 M r a 4 M 3 r 1 43a 6 Mr 1 351a 6 M r a 6 M 3 r a 6 M 4 r 9 + a 8 M r a 8 M 3 r a 8 M 4 r a 8 M 5 r 6 6a 1 M 3 r a 1 M 4 r a 1 M 5 r a 1 M 6 r a 1 M 5 r + 5a 1 M 6 r + 5a 1 M 7) L 6 r 4 6r 19 1Mr a r a Mr a 4 r 15 95a M r 15 18a 4 Mr a 4 M r a 4 M 3 r 1 135a 6 Mr a 6 M r a 6 M 3 r a 6 M 4 r a 8 M r a 8 M 3 r a 8 M 4 r a 8 M 5 r 6 57a 1 M 3 r 6 43a 1 M 4 r 5 837a 1 M 5 r 4 18a 1 M 6 r a 1 M 5 r + 113a 1 M 6 r + 113a 1 M 7) ) L 4 + 4r 6 a M r 3 r 3 + a r + a M ) 3r Mr 1 3a r a Mr a M r 9 + 6a 4 Mr 8 + 7a 4 M r 7 1a 4 M 3 r 6 3a 6 M r 5 48a 6 M 3 r 4 55a 6 M 4 r a 8 M 4 r + 7a 8 M 5) L 88a M r 8 a M r) 3 3 r 3 + a r + a M ) ] + E [ 1aLMr 7 9r 3 + 5a r + 1a M ) L aLMr 6 77r 6 + 3a r 4 36a Mr 3 + 4a 4 r + 8a 4 Mr + 16a 4 M ) L 1 + almr 5 435r 9 86a r a Mr a 4 r a 4 Mr a 4 M r 3 133

147 4.5 Kerr Space-time + 384a 6 Mr + 65a 6 M r + 346a 6 M 3) L 1 almr 4 58r 1 346a r 1 + 4a Mr 9 + 8a 4 r a 4 Mr a 4 M r a 6 Mr a 6 M r a 6 M 3 r a 8 M r a 8 M 3 r + 93a 8 M 4) L 8 + almr 3 78r 15 96a r 13 43a Mr 1 4a 4 r a 4 Mr a 4 M r a 6 Mr a 6 M r a 6 M 3 r a 8 M r a 8 M 3 r a 8 M 4 r a 1 M 3 r + 63a 1 M 4 r + 63a 1 M 5) L 6 almr 5 651r a r a Mr 1 448a 4 r a 4 Mr a 4 M r a 6 Mr a 6 M r 7 585a 6 M 3 r 6 179a 8 M r 5 77a 8 M 3 r 4 198a 8 M 4 r 3 + 7a 1 M 3 r + 88a 1 M 4 r + 88a 1 M 5) ) L 4 + almr 7 a M r 3 r 3 + a r + a M ) 339r a r 7 5a 4 Mr 4 69a 4 M r a 6 M r + 7a 6 M 3) L 144aLMr 9 a M r) 3 3 r 3 + a r + a M ) ] + E 3[ 3a Mr 1 L 14 3r 6 3r 9 + 3a r a Mr 6 3a 4 Mr 4 11a 4 M r a 6 Mr + 96a 6 M r + 1a 6 M 3) L 1 3r 5 3r 1 37a r 1 11a Mr 9 15a 4 r a 4 Mr 7 a 6 r a 4 M r 6 + 4a 6 Mr a 6 M r a 6 M 3 r 3 + 3a 8 Mr a 8 M r + 76a 8 M 3 r + 74a 8 M 4) L 1 + 3r 4 7r 15 14a r 13 8a Mr 1 86a 4 r 11 68a 4 Mr 1 14a 6 r 9 + 7a 4 M r a 6 Mr a 6 M r a 6 M 3 r a 8 Mr a 8 M r a 8 M 3 r a 8 M 4 r a 1 M r a 1 M 3 r + 54a 1 M 4 r + a 1 M 5) L 8 11r 18 58a r 16 31a Mr 15 19a 4 r 14 58a 4 Mr 13 3r 3 44a 6 r a 4 M r a 6 Mr a 6 M r a 6 M 3 r a 8 Mr a 8 M r a 8 M 3 r 7 + 6a 8 M 4 r a 1 M r a 1 M 3 r a 1 M 4 r a 1 M 5 r a 1 M 3 r a 1 M 4 r + 1a 1 M 5 r + 14a 1 M 6) L 6 + 3r 5 r 3 + a r + a M ) 95r 15 15a r 13 78a Mr 1 58a 4 r a 4 Mr 1 + 3a 4 M r a 6 Mr a 6 M r 7 17a 6 M 3 r 6 131a 8 M r 5 33a 8 M 3 r 4 136a 8 M 4 r a 1 M 3 r + 31a 1 M 4 r + 31a 1 M 5) L 4 ) 3r 7 a M r 3 r 3 + a r + a M ) 43r a r 7 4a Mr 6 58a 4 Mr 4 69a 4 M r a 6 M r + 78a 6 M 3) L ) + 4r 9 a M r 3 3 r 3 + a r + a M ) 3 ], 134

148 4.5 Kerr Space-time F r[] = 1 6πr 6 where F E r[] = a Mr 5 3Mr 3 + r 4 ) [a + r r M)] [ a M + ] ) r r 3M) F Er[] E + F Kr[] K a M + r 3/ [ a 4 M + a 3 Mr 3 + a r Mr + r M ) 4aM 3/ r 5/ + Mr 4 r 3/ + a ) M 3M 3 r 33/ 36M 4 r 31/ + 9a M r 31/ 78aM 7/ r ) + 93M 5 r 9/ + 7a M 3 r 9/ + 9a 4 Mr 9/ + 744aM 9/ r a 3 M 5/ r a 6 r 7/ 6M 6 r 7/ 168a M 4 r 7/ 1a 4 M r 7/ 189aM 11/ r a 3 M 7/ r a 5 M 3/ r a M 5 r 5/ + 913a 4 M 3 r 5/ + 15a 6 Mr 5/ + 133aM 13/ r 1 74a 3 M 9/ r 1 44a 5 M 5/ r 1 + 4a 7 Mr 1 + 6a 8 r 3/ 451a M 6 r 3/ 3846a 4 M 4 r 3/ 357a 6 M r 3/ + 47a 3 M 11/ r a 5 M 7/ r a 7 M 3/ r a M 7 r 1/ a 4 M 5 r 1/ + 73a 6 M 3 r 1/ + 81a 8 Mr 1/ 3336a 3 M 13/ r 1 75a 5 M 9/ r 1 138a 7 M 5/ r a 9 Mr 1 786a 4 M 6 r 19/ + 141a 6 M 4 r 19/ 366a 8 M r 19/ 17a 3 M 15/ r a 5 M 11/ r a 7 M 7/ r a 9 M 3/ r a 4 M 7 r 17/ 5631a 6 M 5 r 17/ 1651a 8 M 3 r 17/ + 156a 1 Mr 17/ + 94a 5 M 13/ r 8 1a 7 M 9/ r 8 136a 9 M 5/ r a 4 M 8 r 15/ 5611a 6 M 6 r 15/ + 49a 8 M 4 r 15/ + 88a 1 M r 15/ a 5 M 15/ r 7 186a 7 M 11/ r 7 151a 9 M 7/ r a 11 M 3/ r a 6 M 7 r 13/ + 59a 8 M 5 r 13/ 1686a 1 M 3 r 13/ 88a 5 M 17/ r a 7 M 13/ r 6 + 9a 9 M 9/ r a 11 M 5/ r a 6 M 8 r 11/ 1a 8 M 6 r 11/ 49a 1 M 4 r 11/ + 46a 1 M r 11/ + 358a 7 M 15/ r a 9 M 11/ r 5 79a 11 M 7/ r 5 348a 6 M 9 r 9/ 5398a 8 M 7 r 9/ + 495a 1 M 5 r 9/ + 414a 1 M 3 r 9/ + 36a 7 M 17/ r 4 174a 9 M 13/ r 4 598a 11 M 9/ r 4 + 1a 13 M 5/ r a 8 M 8 r 7/ a 1 M 6 r 7/ 16a 1 M 4 r 7/ 78a 9 M 15/ r a 11 M 11/ r a 13 M 7/ r 3 543a 1 M 7 r 5/ 116a 1 M 5 r 5/ + 4a 14 M 3 r 5/ + 554a 11 M 13/ r 5a 13 M 9/ r + 186a 1 M 6 r 3/ + 144a 14 M 4 r 3/ 188a 13 M 11/ r 4a 14 M 5 r ] 3, 135

149 4.5 Kerr Space-time F K r[] = r 7/ + 4a 15 M 9/), r 3/ 3M r + a M ) 3M 3 r 9/ 4M 4 r 7/ + 9a M r 7/ 78aM 7/ r M 5 r 5/ + 57a M 3 r 5/ + 9a 4 Mr 5/ + 486aM 9/ r a 3 M 5/ r 1 + 3a 6 r 3/ 885a M 4 r 3/ 1a 4 M r 3/ 61aM 11/ r a 3 M 7/ r a 5 M 3/ r a M 5 r 1/ + 688a 4 M 3 r 1/ + 9a 6 Mr 1/ 186a 3 M 9/ r 1 354a 5 M 5/ r 1 + 4a 7 Mr 1 16a M 6 r 19/ 185a 4 M 4 r 19/ 5a 6 M r 19/ + 58a 3 M 11/ r a 5 M 7/ r a 7 M 3/ r 9 58a 4 M 5 r 17/ 158a 6 M 3 r 17/ + 78a 8 Mr 17/ + 1a 3 M 13/ r a 5 M 9/ r 8 998a 7 M 5/ r 8 189a 4 M 6 r 15/ + 358a 6 M 4 r 15/ + 171a 8 M r 15/ 4a 5 M 11/ r 7 + a 7 M 7/ r a 9 M 3/ r a 4 M 7 r 13/ 368a 6 M 5 r 13/ 163a 8 M 3 r 13/ 196a 5 M 13/ r a 7 M 9/ r a 9 M 5/ r 6 639a 6 M 6 r 11/ 347a 8 M 4 r 11/ + 13a 1 M r 11/ + 684a 5 M 15/ r a 7 M 11/ r 5 59a 9 M 7/ r 5 879a 6 M 7 r 9/ + 3a 8 M 5 r 9/ + 186a 1 M 3 r 9/ 684a 7 M 13/ r 4 83a 9 M 9/ r 4 + 6a 11 M 5/ r a 6 M 8 r 7/ + 148a 8 M 6 r 7/ 9a 1 M 4 r 7/ 18a 7 M 15/ r 3 + 5a 9 M 11/ r a 11 M 7/ r 3 147a 8 M 7 r 5/ 45a 1 M 5 r 5/ + 1a 1 M 3 r 5/ + 146a 9 M 13/ r 38a 11 M 9/ r + 69a 1 M 6 r 3/ + 7a 1 M 4 r 3/ 7a 11 M 11/ r 1a 1 M 5 r + 1a 13 M 9/) Electromagnetic Case The regularisation parameters for Kerr space-time in the electromagnetic case are given by r ṙ sgn r F t[-1] =, r a + L ) + a M + r 3 sgn r Er a + r F r[-1] = ) + amae L)) a Mr + r) r a + L ) + a M + r), 3 F θ[-1] =, F φ[-1] =, 4.15) 136

150 4.5 Kerr Space-time where F t[] = πr ṙ r + L + a M r + a ) 3/ a M + a r + L r ) F E t[] E + F K t[]k ), 4.16) Ft[] E = 4aLM ) 4a 4 M + a 4 Mr + a L Mr a Mr 3 a r 4 L r 4 + E 1a 6 M 3 16a 6 M r 7a 6 Mr a 6 r 3 8a 4 L M r a 4 L Mr 4a 4 L r 3 6a 4 M r 3 5a 4 Mr 4 a 4 r 5 15a L 4 Mr 5a L 4 r 3 5a L Mr 4 a L r 5 L 6 r) 3, F K t[] = alm a 4 M a 4 Mr a 4 r a L Mr a L r a Mr 3 a r 4 ) L 4 r L r 4 + E 4a 6 M 3 + 4a 6 M r + a 6 Mr + 1a 4 L M r + 5a 4 L Mr + a 4 M r 3 + a 4 Mr 4 + 4a L 4 Mr + a L Mr 4 a L r 5 L 4 r 5 ), F Er[] E + F Kr[] K ) F r[] = where πr 3 ) 3/, r + L + a M r + a a M + a r + L r ) a Mr + r) 4.17) F E r[] = L 4a 6 M 4 + 8a 6 M 3 r 6a 6 M r 11a 6 Mr 3 a 6 r a 4 L M 3 r + 4a 4 L M r 18a 4 L Mr 3 6a 4 L r 4 + 5a 4 M 3 r 3 + a 4 M r 4 11a 4 Mr 5 3a 4 r 6 + 3a L 4 M r 3a L 4 Mr 3 6a L 4 r 4 + 4a L M r 4 5a L Mr 5 6a L r 6 + 8a M r 6 a Mr 7 a r 8 + 4L 6 Mr 3 L 6 r 4 ) + 6L 4 Mr 5 3L 4 r 6 + L Mr 7 L r 8 aelm 4a 6 M a 6 M r + 18a 6 Mr + 3a 6 r a 4 L M r + 48a 4 L Mr + 1a 4 L r 3 + 4a 4 M r 3 + 4a 4 Mr 4 + 6a 4 r 5 + 3a L 4 Mr + 11a L 4 r 3 + a L Mr 4 + 9a L r 5 + 6a Mr 6 + 3a r 7 + 4L 6 r 3 + 3L 4 r 5 ) + 3L r 7 + E a M + a r + r 3 ) 1a 6 M a 6 M r + 7a 6 Mr + a 6 r 3 + 8a 4 L M r + a 4 L Mr + 4a 4 L r 3 + 6a 4 M r 3 + 5a 4 Mr 4 + a 4 r a L 4 Mr + 5a L 4 r 3 + 5a L Mr 4 + a L r 5 + L 6 r 3 ), 137

151 4.5 Kerr Space-time F K r[] = L 8a 6 M 4 + 1a 6 M 3 r a 6 Mr 3 + a 4 L M 3 r + 8a 4 L M r 4a 4 L Mr 3 + 3a 4 M 3 r 3 + 1a 4 M r 4 6a 4 Mr 5 a 4 r 6 + 8a L 4 M r a L 4 Mr 3 + 4a L M r 4 4a L Mr 5 a L r 6 + 8a M r 6 a Mr 7 ) a r 8 + L 4 Mr 5 L 4 r 6 + L Mr 7 L r 8 + aelm 8a 6 M 3 + 1a 6 M r + 6a 6 Mr + a 6 r 3 + a 4 L M r + 14a 4 L Mr + a 4 L r a 4 M r a 4 Mr 4 + 4a 4 r 5 + 8a L 4 Mr ) + a L 4 r a L Mr 4 + 5a L r 5 + 6a Mr 6 + 3a r 7 + L 4 r 5 + 3L r 7 E a M + a r + r 3 ) 4a 6 M 3 + 4a 6 M r + a 6 Mr + 1a 4 L M r + 5a 4 L Mr + a 4 M r 3 + a 4 Mr 4 + 4a L 4 Mr + a L Mr 4 a L r 5 L 4 r 5 ), F θ[] = 4.18) F φ[] = Lṙ F E φ[] E + F K φ[] K ) πr r + L + a M r + a ) 1/ a M + a r + L r ), 4.19) where F E φ[] = 14a 4 M + 11a 4 Mr + a 4 r + 11a L Mr + 4a L r + 4a Mr 3 + a r 4 + L 4 r + L r 4, F K φ[] = 4a 4 M a 4 Mr a L Mr 4a Mr 3 a r 4 L r 4. As with the Scalar case, F a[] proves too large to include in paper format and so is available electronically 7), although we provide F r[] for circular orbits below. F r[] = where 1 6πMr 11/ am 5/ r 3/ + M rr 3M) a ) M + r 3/ [a + r r M)] a M + ) ) r r 3M) F Er[] E + F Kr[] K [ a 4 M + a 3 Mr 3 + a r M + Mr + r ) 4aM 3/ r 5/ + Mr ) ] 3, 138

152 4.5 Kerr Space-time Fr[] E = r 3/ + a ) M 3M 3 r 33/ 96M 4 r 31/ 9a M r 31/ + 438aM 7/ r 15 F E r[] = r M 5 r 9/ 399a M 3 r 9/ 9a 4 Mr 9/ 88aM 9/ r a 3 M 5/ r 14 3a 6 r 7/ 588M 6 r 7/ + 18a M 4 r 7/ 8a 4 M r 7/ + 58aM 11/ r 13 37a 3 M 7/ r a 5 M 3/ r a M 5 r 5/ + 599a 4 M 3 r 5/ 13a 6 Mr 5/ + 88aM 13/ r 1 754a 3 M 9/ r a 5 M 5/ r 1 4a 7 Mr 1 6a 8 r 3/ 1165a M 6 r 3/ a 4 M 4 r 3/ + 69a 6 M r 3/ a 3 M 11/ r a 5 M 7/ r 11 58a 7 M 3/ r a M 7 r 1/ a 4 M 5 r 1/ a 6 M 3 r 1/ 189a 8 Mr 1/ + 768a 3 M 13/ r 1 65a 5 M 9/ r a 7 M 5/ r 1 48a 9 Mr a 4 M 6 r 19/ 169a 6 M 4 r 19/ 33a 8 M r 19/ 184a 3 M 15/ r a 5 M 11/ r 9 353a 7 M 7/ r 9 75a 9 M 3/ r a 4 M 7 r 17/ a 6 M 5 r 17/ a 8 M 3 r 17/ 156a 1 Mr 17/ 998a 5 M 13/ r a 7 M 9/ r a 9 M 5/ r 8 576a 4 M 8 r 15/ 16697a 6 M 6 r 15/ 16348a 8 M 4 r 15/ 1368a 1 M r 15/ + 138a 5 M 15/ r a 7 M 11/ r a 9 M 7/ r 7 64a 11 M 3/ r a 6 M 7 r 13/ a 8 M 5 r 13/ + 346a 1 M 3 r 13/ 396a 5 M 17/ r 6 : 3546a 7 M 13/ r a 9 M 9/ r 6 13a 11 M 5/ r a 6 M 8 r 11/ a 8 M 6 r 11/ + 654a 1 M 4 r 11/ 46a 1 M r 11/ + 588a 7 M 15/ r a 9 M 11/ r a 11 M 7/ r 5 516a 6 M 9 r 9/ 9146a 8 M 7 r 9/ 3435a 1 M 5 r 9/ + 54a 1 M 3 r 9/ + 6a 7 M 17/ r 4 418a 9 M 13/ r 4 453a 11 M 9/ r 4 1a 13 M 5/ r a 8 M 8 r 7/ + 699a 1 M 6 r 7/ + 7a 1 M 4 r 7/ 188a 9 M 15/ r a 11 M 11/ r a 13 M 7/ r 3 969a 1 M 7 r 5/ 336a 1 M 5 r 5/ 4a 14 M 3 r 5/ + 13a 11 M 13/ r 116a 13 M 9/ r + 354a 1 M 6 r 3/ + 88a 14 M 4 r 3/ 364a 13 M 11/ r 48a 14 M 5 r + 48a 15 M 9/), r 3/ 3M r + a M ) 3M 3 r M 4 r a M r aM 7/ r 7/ 55M 5 r a M 3 r a 4 Mr aM 9/ r 5/ 388a 3 M 5/ r 5/ + 3a 6 r 1 69a M 4 r 1 + 1a 4 M r 1 + 5aM 11/ r 3/ 46a 3 M 7/ r 3/ 54a 5 M 3/ r 3/ 763a M 5 r 11 3a 4 M 3 r a 6 Mr a 3 M 9/ r 1/ 738a 5 M 5/ r 1/ + 954a M 6 r a 4 M 4 r 1 4a 6 M r 1 94a 3 M 11/ r 19/ 139

153 4.5 Kerr Space-time + 161a 5 M 7/ r 19/ + 474a 7 M 3/ r 19/ 3648a 4 M 5 r 9 314a 6 M 3 r 9 + a 8 Mr 9 564a 3 M 13/ r 17/ + 656a 5 M 9/ r 17/ 366a 7 M 5/ r 17/ a 4 M 6 r a 6 M 4 r a 8 M r 8 196a 5 M 11/ r 15/ 338a 7 M 7/ r 15/ + 996a 9 M 3/ r 15/ 1881a 4 M 7 r 7 764a 6 M 5 r 7 933a 8 M 3 r a 5 M 13/ r 13/ a 7 M 9/ r 13/ + 156a 9 M 5/ r 13/ 187a 6 M 6 r a 8 M 4 r a 1 M r 6 144a 5 M 15/ r 11/ 16a 9 M 7/ r 11/ a 6 M 7 r a 8 M 5 r a 1 M 3 r 5 868a 7 M 13/ r 9/ + 494a 9 M 9/ r 9/ + 94a 11 M 5/ r 9/ 174a 6 M 8 r 4 171a 8 M 6 r a 1 M 4 r a 7 M 15/ r 7/ + 149a 9 M 11/ r 7/ 36a 11 M 7/ r 7/ 111a 8 M 7 r a 1 M 5 r a 1 M 3 r 3 6a 9 M 13/ r 5/ 5774a 11 M 9/ r 5/ + 175a 1 M 6 r 48a 1 M 4 r + 986a 11 M 11/ r 3/ + 864a 13 M 7/ r 3/ 876a 1 M 5 r 17a 7 Mr ) 11/ + 144a 14 M 4 156a 13 ) M 9 r + 4a Mr Gravitational Case Self-Force Regularisation r ṙ sgn r F t[-1] =, r a + L ) + a M + r ) F r[-1] = sgn r Er a + r) + amae L)) a Mr + r) r a + L ) + a M + r), ) F θ[-1] =, F φ[-1] =, 4.133) where F t[] = πr ) ṙ Ft[] E E + F t[] K K ) 3/, 4.134) r + L + a M r + a a M + a r + L r ) 14

154 4.5 Kerr Space-time Ft[] E = 4aLM ) 4a 4 M + a 4 Mr + a L Mr a Mr 3 a r 4 L r 4 + E 1a 6 M a 6 M r + 7a 6 Mr + a 6 r 3 + 8a 4 L M r + a 4 L Mr + 4a 4 L r 3 + 6a 4 M r 3 + 5a 4 Mr 4 + a 4 r a L 4 Mr + 5a L 4 r 3 + 5a L Mr 4 + a L r 5 + L 6 r) 3, F K t[] = alm a 4 M a 4 Mr a 4 r a L Mr a L r a Mr 3 a r 4 ) L 4 r L r 4 + E 4a 6 M 3 4a 6 M r a 6 Mr 1a 4 L M r 5a 4 L Mr a 4 M r 3 a 4 Mr 4 4a L 4 Mr a L Mr 4 + a L r 5 + L 4 r 5 ), F Er[] E + F Kr[] K ) F r[] = where πr 6 ) 3/, r + L + a M r + a a M + a r + L r ) a Mr + r) 4.135) F E r[] = +L r 3 4a 6 M 4 8a 6 M 3 r + 6a 6 M r + 11a 6 Mr 3 + a 6 r 4 56a 4 L M 3 r 4a 4 L M r + 18a 4 L Mr 3 + 6a 4 L r 4 5a 4 M 3 r 3 a 4 M r a 4 Mr 5 + 3a 4 r 6 3a L 4 M r + 3a L 4 Mr 3 + 6a L 4 r 4 4a L M r 4 + 5a L Mr 5 + 6a L r 6 8a M r 6 + a Mr 7 + a r 8 ) 4L 6 Mr 3 + L 6 r 4 6L 4 Mr 5 + 3L 4 r 6 L Mr 7 + L r 8 + aelmr 3 4a 6 M a 6 M r + 18a 6 Mr + 3a 6 r a 4 L M r + 48a 4 L Mr + 1a 4 L r 3 + 4a 4 M r 3 + 4a 4 Mr 4 + 6a 4 r 5 + 3a L 4 Mr + 11a L 4 r 3 + a L Mr 4 + 9a L r 5 + 6a Mr 6 + 3a r 7 + 4L 6 r 3 + 3L 4 r 5 ) + 3L r 7 ) E r 3 a M + a r + r 3 1a 6 M a 6 M r + 7a 6 Mr + a 6 r 3 + 8a 4 L M r + a 4 L Mr + 4a 4 L r 3 + 6a 4 M r 3 + 5a 4 Mr 4 + a 4 r a L 4 Mr + 5a L 4 r 3 + 5a L Mr 4 + a L r 5 + L 6 r 3 ), F K r[] = L r 3 8a 6 M 4 1a 6 M 3 r + a 6 Mr 3 a 4 L M 3 r 8a 4 L M r + 4a 4 L Mr 3 3a 4 M 3 r 3 1a 4 M r 4 + 6a 4 Mr 5 + a 4 r 6 8a L 4 M r + a L 4 Mr 3 4a L M r 4 + 4a L Mr 5 + a L r 6 8a M r 6 + a Mr 7 ) + a r 8 L 4 Mr 5 + L 4 r 6 L Mr 7 + L r 8 aelmr 3 8a 6 M 3 + 1a 6 M r + 6a 6 Mr + a 6 r 3 + a 4 L M r 141

155 4.5 Kerr Space-time + 14a 4 L Mr + a 4 L r a 4 M r a 4 Mr 4 + 4a 4 r 5 + 8a L 4 Mr ) + a L 4 r a L Mr 4 + 5a L r 5 + 6a Mr 6 + 3a r 7 + L 4 r 5 + 3L r 7 ) E r 3 a M + a r + r 3 4a 6 M 3 4a 6 M r a 6 Mr 1a 4 L M r 5a 4 L Mr a 4 M r 3 a 4 Mr 4 4a L 4 Mr a L Mr 4 + a L r 5 + L 4 r 5 ), F θ[] =, 4.136) F Eφ[] E + F Kφ[] K ) F φ[] = where πr 6 ) 3/, r + L + a M r + a a M + a r + L r ) a Mr + r) 4.137) F E φ[] = 14a 4 M 11a 4 Mr a 4 r 11a L Mr 4a L r 4a Mr 3 a r 4 L 4 r L r 4, F K φ[] = 4a 4 M + a 4 Mr + a L Mr + 4a Mr 3 + a r 4 + L r 4. As with the scalar and electromagnetic cases, F a[] is too large for paper format and so is available electronically 7). As outlined above, we do give F r[] for circular orbits, F r[] = 1 where πr 11/ am 5/ r 3/ + M rr 3M) a ) M + r 3/ [a + r r M)] [ a M + ] r r 3M) ) [a + r r M)] 1 Fr[] E E + F r[] K K [ a 4 M + a 3 ], Mr 3 + a r M + Mr + r) 4aM 3/ r 5/ + Mr 4 14

156 4.5 Kerr Space-time Fr[] E = a ) M + r 3/ a Mr + r) 3a 11 M 5/ 48a 1 M 3 r F K r[] = r a 1 M r 3/ + 478a 9 M 3/ r 3 14a 9 M 5/ r 19a 9 M 7/ r + 88a 8 M 4 r 3/ 99a 8 M 3 r 5/ 473a 8 M r 7/ + 779a 8 Mr 9/ 9a 7 M 3/ r 5 496a 7 M 5/ r 4 + 8a 7 M 7/ r a 7 M 9/ r + 54a 7 Mr 6 576a 6 M 5 r 5/ + 4a 6 M 4 r 7/ + 489a 6 M 3 r 9/ 463a 6 M r 11/ + 7a 6 Mr 13/ + 17a 6 r 15/ 1846a 5 M 3/ r a 5 M 5/ r a 5 M 7/ r 5 195a 5 M 9/ r 4 56a 5 M 11/ r a 5 Mr a 4 M 6 r 7/ 1664a 4 M 5 r 9/ 395a 4 M 4 r 11/ + 638a 4 M 3 r 13/ 75a 4 M r 15/ 46a 4 Mr 17/ + 15a 4 r 19/ 14a 3 M 3/ r a 3 M 5/ r a 3 M 7/ r a 3 M 9/ r a 3 M 11/ r 5 144a 3 Mr a M 5 r 13/ + 674a M 4 r 15/ + 588a M 3 r 17/ 131a M r 19/ + 348a Mr 1/ a r 3/ 148aM 3/ r aM 5/ r 1 19M 4 r 19/ + 87M 3 r 1/ a M 3M r + r 3/ 137aM 7/ r aM 9/ r 8 131M r 3/ + 18Mr 5/ ), ) 144a 1 M 96a 11 M 5 r + 576a 11 Mr ) 3/ 91a 1 M 3 r + 81a 1 M r + 864a 1 Mr a 9 M 3/ r 7/ + 58a 9 M 7/ r 3/ + 576a 9 Mr 9 41a 9 Mr ) 5/ + 5a 8 M 4 r + 313a 8 M 3 r 3 534a 8 M r a 8 Mr a 8 r 6 38a 7 M 3/ r 11/ 3816a 7 M 5/ r 9/ 1176a 7 M 9/ r 5/ + 9a 7 Mr a 7 Mr ) 7/ 1764a 6 M 5 r 3 6a 6 M 4 r a 6 M 3 r a 6 M r 6 + 6a 6 Mr a 6 r 8 58a 5 M 3/ r 15/ + 44a 5 M 5/ r 13/ + 354a 5 M 7/ r 11/ + 79a 5 M 11/ r 7/ + 4a 5 Mr a 5 Mr ) 9/ + 384a 4 M 6 r 4 + 6a 4 M 5 r a 4 M 4 r a 4 M 3 r a 4 M r 8 111a 4 Mr a 4 r 1 4a 3 M 3/ r 19/ + 98a 3 M 5/ r 17/ 4866a 3 M 7/ r 15/ + 137a 3 M 9/ r 13/ 144a 3 Mr a 3 Mr ) 11/ + 178a M 5 r a M 4 r a M 3 r 9 175a M r a Mr 11 a r 1 148aM 3/ r 3/ + 718aM 5/ r 1/ 14aM 7/ r 19/ + 384aM 9/ r 17/ ) 19M 4 r 1 + 8M 3 r 11 18M r Mr

157 4.5 Kerr Space-time huu Regularisation H [] = K ) 1/, 4.138) π r + L + a M r + a H [1] =, 4.139) H E[] E + HK[] K ) H [] = 3πr 7 r + L + a M r + a ) 3/ a M + a r + L r ) 3, 4.14) where H E [] = 1Mr 5 a 1 + 9M r 4 a M 3 r 3 a M 4 r a M 5 r a 1 4r 8 a 1 4Mr 7 a 1 114L M 6 a 1 114M r 6 a 1 88M 3 r 5 a L Mr 5 a 1 47M 4 r 4 a 1 + 3L M r 4 a 1 364M 5 r 3 a L M 3 r 3 a 1 96M 6 r a L M 4 r a 1 96L M 5 r a 1 48r 1 a 8 4Mr 9 a 8 1L r 8 a M r 8 a 8 87M 3 r 7 a 8 88L Mr 7 a 8 448M 4 r 6 a 8 781L M r 6 a 8 67M 5 r 5 a 8 477L M 3 r 5 a 8 + 7L 4 Mr 5 a 8 47L M 4 r 4 a 8 + 9L 4 M r 4 a 8 144L M 5 r 3 a 8 144L 4 M 3 r 3 a 8 98L 4 M 4 r a 8 11L 4 M 5 r a 8 4r 1 a 6 168Mr 11 a 6 195L r 1 a 6 456M r 1 a 6 48M 3 r 9 a 6 186L Mr 9 a 6 4L 4 r 8 a 6 96M 4 r 8 a 6 119L M r 8 a 6 158L M 3 r 7 a 6 198L 4 Mr 7 a 6 84L M 4 r 6 a L 4 M r 6 a 6 696L 4 M 3 r 5 a L 6 Mr 5 a L 4 M 4 r 4 a 6 19L 6 M r 4 a 6 13L 6 M 3 ra L 6 M 4 ra 6 75L r 1 a 4 46L Mr 11 a 4 97L 4 r 1 a 4 84L M r 1 a L M 3 ra L 4 Mra 9 4 4L 6 ra L 4 M ra L 4 M 3 ra 7 4 4L 6 Mra L 6 M ra L 6 M 3 r 5 a 4 + 1L 8 Mr 5 a 4 19L 8 M r 4 a 4 444L 8 M 3 r 3 a 4 78L 4 r 1 a + 36L 4 Mr 11 a 1L 6 r 1 a + 384L 4 M r 1 a + 198L 6 Mr 9 a 1L 8 r 8 a + 19L 6 M r 8 a + 16L 8 Mr 7 a + 67L 8 M r 6 a 48L 1 M ra 4 7L 6 r L 6 Mr 11 51L 8 r 1 + L 8 Mr 9 ) 4L 1 r L 1 Mr 7 + 4aELM 71a 4 r 11 75L 4 r a L r 11 6a 4 Mr 1 7a L Mr 1 147a 6 r 9 159L 6 r 9 465a L 4 r 9 453a 4 L r 9 4a 4 M r 9 811a 6 Mr 8 144

158 4.5 Kerr Space-time 873a L 4 Mr a 4 L Mr 8 76a 8 r 7 78L 8 r 7 31a L 6 r 7 46a 4 L 4 r 7 36a 6 L r 7 149a 6 M r 7 154a 4 L M r 7 91a 6 M 3 r 6 5a 8 Mr 6 549a L 6 Mr 6 16a 4 L 4 Mr a 6 L Mr 6 133a 8 M r a 4 L 4 M r 5 691a 6 L M r 5 146a 8 M 3 r a 6 L M 3 r a 1 Mr 4 + 4a L 8 Mr a 4 L 6 Mr a 6 L 4 Mr a 8 L Mr 4 588a 8 M 4 r a 1 M r 3 + a 4 L 6 M r a 6 L 4 M r a 8 L M r a 1 M 3 r + 738a 6 L 4 M 3 r a 8 L M 3 r + 114a 1 M 4 r + 156a 8 L M 4 r + 55a 1 M 5) 3E 368M 6 a 1 + 4Mr 5 a M r 4 a 1 + 8M 3 r 3 a M 4 r a 1 + 8M 5 r a 1 9r 8 a 1 14Mr 7 a 1 64M r 6 a M 3 r 5 a 1 + 1L Mr 5 a 1 147M 4 r 4 a L M r 4 a 1 68M 5 r 3 a L M 3 r 3 a L M 4 r a L M 5 r a 1 18r 1 a 8 4Mr 9 a 8 35L r 8 a 8 91M r 8 a 8 149M 3 r 7 a 8 438L Mr 7 a 8 884M 4 r 6 a L M r 6 a 8 64L M 3 r 5 a 8 + 1L 4 Mr 5 a 8 141L M 4 r 4 a L 4 M r 4 a L 4 M 3 r 3 a L 4 M 4 r a 8 9r 1 a 6 96Mr 11 a 6 55L r 1 a 6 78M r 1 a 6 44M 3 r 9 a 6 596L Mr 9 a 6 51L 4 r 8 a L M r 8 a L M 3 r 7 a 6 57L 4 Mr 7 a L 4 M r 6 a 6 154L 4 M 3 r 5 a 6 + 4L 6 Mr 5 a 6 + 8L 6 M ra L 6 M 3 ra 3 6 L r 1 a 4 176L Mr 11 a 4 56L 4 r 1 a 4 66L M r 1 a 4 51L 4 Mra L 6 ra L 4 M ra L 6 Mra L 6 M ra L 8 M ra L 4 r 1 a 8L 4 Mr 11 ) 19L 6 r 1 a 14L 6 Mra 9 8L 8 ra 8 68L 8 Mra 7 L 6 r 1, a H[] K = r 3 4r 5 a Mra M ra M 3 ra 1 + 3M 4 r a ra Mra L M 5 a 8 + 1L ra M ra M 3 r 4 a L Mr 4 a 8 91M 4 r 3 a L M r 3 a 8 115M 5 r a 8 434L M 3 r a 8 17L M 4 r a r 9 a Mr 8 a L r 7 a M r 7 a 6 168M 3 r 6 a L Mr 6 a 6 + 4L 4 r 5 a 6 576M 4 r 5 a L M r 5 a 6 7L M 3 r 4 a L 4 Mr 4 a 6 337L M 4 r 3 a 6 45L 4 M r 3 a 6 896L 4 M 3 r a 6 5L 4 M 4 r a 6 + 1L r 9 a L Mr 8 a L 4 r 7 a 4 384L M r 7 a 4 14L M 3 r 6 a L 4 Mr 6 a 4 + 4L 6 r 5 a 4 17L 4 M r 5 a 4 376L 4 M 3 r 4 a L 6 Mr 4 a 4 155L 6 M r 3 a 4 186L 6 M 3 r a L 4 r 9 a 18L 4 Mr 8 a + 381L 6 r 7 a 888L 4 M r 7 a 16L 6 Mr 6 a + 1L 8 r 5 a 145

159 4.5 Kerr Space-time 1776L 6 M ra 5 19L 8 Mra 4 696L 8 M ra 3 + 7L 6 r 9 198L 6 Mr 8 ) + 96L 8 r 7 94L 8 Mr 6 + 4L 1 r 5 96L 1 Mr 4 + 4arELM 3 456M 4 a ra Mra M ra M 3 r a ra Mra L ra M ra M 3 ra L Mra L M ra L M 3 r a ra Mra L ra M ra L Mra L 4 ra L M ra L 4 Mra L 4 M ra L ra 8 + 7L Mra L 4 ra 6 ) + 69L 4 Mra L 6 ra L 6 Mra L 4 r L 6 r L 8 r 4 3rE 3 34M 5 a 1 + 8ra Mra M ra M 3 ra M 4 r a 1 + 3ra Mra L ra M ra M 3 ra L Mra M 4 ra L M ra L M 3 ra L M 4 r a 8 + 4ra Mra L ra M ra M 3 ra L Mra L 4 ra L M ra L M 3 ra L 4 Mra L 4 M ra L 4 M 3 ra L ra L Mra L 4 ra L M ra L 4 Mra L 6 ra L 4 M ra L 6 Mra L 6 M ra L 4 ra L 4 Mra L 6 ra 7 + 3L 6 Mra 6 + 8L 8 ra 5 + 3L 8 Mra L 6 r L 8 r) Results Using the mode sum, to date, researchers have only been able to produce values for the Kerr retarded field and resulting self-force in the scalar case. As there is no known decomposition of the metric perturbation into tensor harmonic modes that separates the field equations in Kerr Lorentz gauge metric perturbation, higher spins are yet to be calculated. Therefore it is not possible to demonstrate the validity of our electromagnetic and gravitational parameters. However, as with our Schwarzschild parameters, deriving the expressions by independent methods gives us confidence in our results. The only other check we can do is to set the spin of the black hole to be zero, i.e., a =, and compare with our Schwarzschild results which we know to be correct. All of the parameters given have passed these tests. For the scalar case we were able to use the numerical data from Barack and Warburton 69), for calculating the retarded field for Kerr scalar eccentric orbits 146

160 4.6 Regularisation Parameters to Date 1.1 modes of the scalar self force M q F r F r F r F r 1 F r F r 1 F r F r 1 F r F r F r F r F r 1 F r F r F r F r 1 F r F r F r 4 F r Figure 4.4: Regularization of the radial component of the scalar self-force in Kerr space-time for the case of a scalar particle as F l r against l. This data is for an eccentric geodesic with energy E =.95549, spin a = 1/M and angular momentum L = In decreasing slope the above lines represent the unregularised self-force black), self-force regularised by subtracting from it in turn the cumulative sum of F l r[ 1] blue), F l r[] red), F l r[] yellow), F l r[4] green). in the equatorial plane. regularised their data as illustrated in Fig We were able to show that our parameters successfully 4.6 Regularisation Parameters to Date In Table 4.3, we summarise all the regularisation parameters that are known to date and the authors of each parameter. We can see that the work of this thesis has greatly added to the already existing data base of parameters, which in turn 147

161 4.6 Regularisation Parameters to Date Cases F a[-1] F a[] F a[] F a[4] F a[6] Schwarzschild scalar BO BO DMW / HP HOW HOW Schwarzschild electromagnetic BO BO HP HOW HOW Schwarzschild gravity BO BO HOW HOW Schwarzschild huu BO HOW HOW Kerr scalar BO BO HOW HOW Kerr electromagnetic BO BO HOW Kerr gravity BO BO HOW Kerr huu HOW HOW Table 4.3: This table represents the current regularisation parameters that are known for each case in Schwarzschild and Kerr space-times. We have indicated which authors first derived the regularisation parameters. BO is Barack and Ori, DMW is Detweiler, Messaritaki and Whiting, HP is Haas and Poisson and HOW is Heffernan, Ottewill and Wardell. DMW and HP share the authorship for Schwarzschild scalar F a[-] as DMW produced the first expressions for circular orbits while HP extended this to the elliptic orbits. The HOW results were produced as part of this thesis. is dramatically reducing the computation time necessary for accurate predictions of the self-force. By deriving our expressions using two independent methods we were able to be confident in our results. In particular cases, this was boosted by the success of the parameters in regularising numerical data for the retarded field. The alternative approaches used both came with their advantages and disadvantages. For the lower orders, computation time as well as the personal time of the researcher were very much on par, however, the covariant method produces more elegant expressions. However, when we increase the order of our expansions, the covariant technique becomes more time consuming, both computationally and personally for the researcher. In comparison, working in coordinates has a sense of automation attached to it for the higher orders. Some finesse is required for the very high orders 148

162 4.6 Regularisation Parameters to Date when the calculation may get stuck, but this usually only requires minutes of the researchers time. Such blips usually occur due to the high number of unknown parameters, memory issues and certain formalisms that sufficiently slow our software Mathematica) down, most of which can be circumnavigated with experience and knowledge of the software. 149

163 Chapter 5 Effective Source As another application of our high-order expansions of the Detweiler-Whiting singular field, we consider its use in the effective source approach to calculating the self-force. The effective source approach independently proposed by Barack and Golbourn 54) and by Vega and Detweiler 55) relies on knowledge of the singular field to derive an equation for a regularized field that gives the self-force without any need for post-processed regularization. Where the mode-sum method solves for the retarded field, which is in turn regularised by using the singular field after it is numerically calculated, the effective source method numerically calculates the regularised field directly from the wave equation. There are numerous advantages for this method. By design, there are no delta functions or singularities - a desirable attribute for numerical calculations. There is no post-processed construction of the regularized field - that is we don t have to cancel two large quantities ϕ A ret) and ϕa S) which when carried out numerically can lead to large round-off errors. It does not rely on the separability of the perturbation equations, an advantage for time-domain calculations in Kerr space-time where the perturbation equations are not fully separable. 15

164 5.1 Methods 5.1 Methods Following directly from Eq..164), we have D A Bϕ B R) = D A Bϕ B ret) D A Bϕ B S), 5.1) where the definition of the singular field and Eq..158) gives us D A Bϕ B R) = 4πQ u A δ 4 x, zτ )) dτ Dϕ A S) = 5.) However, as the singular field is only well-defined in the neighbourhood of the particle, it should be noted that Eq. 5.) is also only valid in this neighbourhood. Therefore a similar equation is required for when we are not in this region. The two methods introduced to do this are the window function as first described by Vega and Detweiler 55) and the world tube, which was proposed by Barack and Golbourn 54) World Tube Method The world tube method resolves the problem of lack of global definition of the singular field by introducing a world tube - inside the world tube one solves for the regular field, ϕ B R), and outside the world tube one solves for the retarded field, ϕb R). If we consider the boundary of this world tube to be at x and world line at x, we then have D A Bϕ B R) = 4πQ u A δ 4 x, zτ )) dτ Dϕ A S) =, x x < x D A Bϕ B ret) =, x x > x 5.3) where the matching condition ϕ A R) = ϕa ret) ϕa S) is imposed on the boundary. As the use of the singular field is now restricted to the region near the particle, its non-global definition is no longer a problem Window Function Method The window function method involves the use of a globally defined approximate singular field, ϕ A S). This is obtained by the use of a smooth window function, W, 151

165 5.1 Methods that is chosen so that at the particle, ϕ A S) is given by the exact singular field, while away from the particle ϕ A S) quickly diminishes and only the retarded field remains on and past the boundary, x, i.e., ϕ A S) = W ϕa S). In this manner, the world function method unites the two regions of the world tube method avoiding the use of two seperate computational domains. The definition of an approximate singular field, gives a new definition of an approximate regular field, ϕ A R) = ϕa ret) ϕa S), and with it, a slightly altered wave equation, D A B ϕ B R) = 4πQ ) u A δ 4 x, zτ )) dτ D A B W ϕ B S) 5.4) Due to the nature of the window function, Eq. 5.4) is not restricted to a certain region of space unlike Eqs. 5.) and 5.3). By design, the window function will have certain restrictions - the resulting approximate regular field must give the correct self-force, the approximate singular field must be equal to the exact singular field at the particle and equal to zero far from the particle. These give W = 1 + f, where f = O ɛ n ) and p a Af = O ɛ n 1 ) where p a A is that of Eq..16). W is smooth W = for x > x where n is an integer >. This lower bound comes from the equation defining the self-force as F a = p a A ϕ A R) = p a A ϕ A ret) ϕs)) A p a A f) ϕ A S) f ) p a Aϕ A S) ) p a A ϕ A ret) ϕ A S) as x x 5.5) where the last equality holds at the particle which requires the last two terms of the second equality as ɛ. As stated before, one of the key advantages of the effective source is its avoidance of δ functions - this can easily be shown when considering the window function approach. The effective source is defined by, 15

166 5. Effective Source from the Singular Field Seff A = D A B ϕ B R) = D A B ϕ B S) 4πQ ) = D A B fϕ B S) u A δ 4 x, zτ )) dτ = f) ϕ A S) a f a ϕ A S) f ϕ A S) = f) ϕ A S) a f a ϕ A S) 5.6) where the last two equalities are for a vacuum space-time, i.e., D A B = δ A B and we take advantage of f being zero at the particle and ϕ A S) the particle. begin zero away from 5. Effective Source from the Singular Field The effective source, defined by Eq. 5.6), for sufficiently good approximations gives { Seff A at the particle) = D A B ϕ B S) away from the particle), 5.7) where at the particle follows from f = in Eq. 5.6) and the fact that ϕ B S) solves the same inhomogeneous wave equation as ϕ B ret)). Away from the particle holds as D A Bϕ B ret) = in the first line of Eq. 5.6). If the singular field is known exactly, then the regularized field is totally regular and is a solution of the homogeneous wave equation. In reality, exact expressions for the singular field can only be obtained for very simple space-times. More generally, the best one can do is an approximation such as that given in Sec. 3.1, which we now define as ϕ A S), so we have ϕ A S) = [1 + O ɛ) n ] ϕ A S) 5.8) which satisfies the above conditions for W. Our calculated high-order expansion of the singular field can, therefore, be placed into Eq. 5.7) as the approximate singular field. In Figs. 3.,3.3, 5.1 and 5. we show the result of applying our expansions to the case of a scalar particle on a circular geodesic of radius r = 1M in Schwarzschild 153

167 5. Effective Source from the Singular Field 6. S eff H-1L !!!!!!!!!!!!!! Dr!!!!!!!!!! Dr!. S eff H1L S eff HL !! Dr! Dr! S eff H3L Dr!!! Dr! Dr Dr!! Figure 5.1: Effective source for the approximate singular field of order Oɛ 1 ) top S eff H5L. 1-9 S eff H4L S eff H6L S eff HL left) to Oɛ 6 ) bottom right). Shown is the scalar case of a circular geodesic of radius r = 1M in Schwarzschild space-time. 154

168 5. Effective Source from the Singular Field!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Figure 5.: Effective source for the approximate singular field, S eff n), in the region of the particle of order Oɛ 1 ) top left) to Oɛ 6 ) bottom right). Shown is the scalar case of a circular geodesic of radius r = 1M in Schwarzschild space-time. 155