Numerical relativity is one of the branches of general relativity that uses numerical methods and algorithms to solve and analyze problems. To this end, supercomputers are often employed to study black holes, gravitational waves, neutron stars and many other phenomena governed by Einstein's theory of general relativity. A currently active field of research in numerical relativity is the simulation of relativistic binaries and their associated gravitational waves. Other branches are also active.
A primary goal of numerical relativity is to study spacetimes whose exact form is not known. The spacetimes so found computationally can either be fully dynamical, stationary or static and may contain matter fields or vacuum. In the case of stationary and static solutions, numerical methods may also be used to study the stability of the equilibrium spacetimes. In the case of dynamical spacetimes, the problem may be divided into the initial value problem and the evolution, each requiring different methods.
Numerical relativity is applied to many areas, such as cosmological models, critical phenomena, perturbed black holes and neutron stars, and the coalescence of black holes and neutron stars, for example. In any of these cases, Einstein's equations can be formulated in several ways that allow us to evolve the dynamics. While Cauchy methods have received a majority of the attention, characteristic and Regge calculus based methods have also been used. All of these methods begin with a snapshot of the gravitational fields on some hypersurface, the initial data, and evolve these data to neighboring hypersurfaces.
Like all problems in numerical analysis, careful attention is paid to the stability and convergence of the numerical solutions. In this line, much attention is paid to the gauge conditions, coordinates, and various formulations of the Einstein equations and the effect they have on the ability to produce accurate numerical solutions.
Numerical relativity research is distinct from work on classical field theories as many techniques implemented in these areas are inapplicable in relativity. Many facets are however shared with large scale problems in other computational sciences like computational fluid dynamics, electromagnetics, and solid mechanics. Numerical relativists often work with applied mathematicians and draw insight from numerical analysis, scientific computation, partial differential equations, and geometry among other mathematical areas of specialization.
Foundations in theory
Albert Einstein published his theory of general relativity in 1915. It, like his earlier theory of special relativity, described space and time as a unified spacetime subject to what are now known as the Einstein field equations. These form a set of coupled nonlinear partial differential equations (PDEs). After more than 100 years since the first publication of the theory, relatively few closed-form solutions are known for the field equations, and, of those, most are cosmological solutions that assume special symmetry to reduce the complexity of the equations.
The field of numerical relativity emerged from the desire to construct and study more general solutions to the field equations by approximately solving the Einstein equations numerically. A necessary precursor to such attempts was a decomposition of spacetime back into separated space and time. This was first published by Richard Arnowitt, Stanley Deser, and Charles W. Misner in the late 1950s in what has become known as the ADM formalism. Although for technical reasons the precise equations formulated in the original ADM paper are rarely used in numerical simulations, most practical approaches to numerical relativity use a "3+1 decomposition" of spacetime into three-dimensional space and one-dimensional time that is closely related to the ADM formulation, because the ADM procedure reformulates the Einstein field equations into a constrained initial value problem that can be addressed using computational methodologies.
At the time that ADM published their original paper, computer technology would not have supported numerical solution to their equations on any problem of any substantial size. The first documented attempt to solve the Einstein field equations numerically appears to be Hahn and Lindquist in 1964, followed soon thereafter by Smarr and by Eppley. These early attempts were focused on evolving Misner data in axisymmetry (also known as "2+1 dimensions"). At around the same time Tsvi Piran wrote the first code that evolved a system with gravitational radiation using a cylindrical symmetry. In this calculation Piran has set the foundation for many of the concepts used today in evolving ADM equations, like "free evolution" versus "constrained evolution",[clarification needed] which deal with the fundamental problem of treating the constraint equations that arise in the ADM formalism. Applying symmetry reduced the computational and memory requirements associated with the problem, allowing the researchers to obtain results on the supercomputers available at the time.
The first realistic calculations of rotating collapse were carried out in the early eighties by Richard Stark and Tsvi Piran in which the gravitational wave forms resulting from formation of a rotating black hole were calculated for the first time. For nearly 20 years following the initial results, there were fairly few other published results in numerical relativity, probably due to the lack of sufficiently powerful computers to address the problem. In the late 1990s, the Binary Black Hole Grand Challenge Alliance successfully simulated a head-on binary black hole collision. As a post-processing step the group computed the event horizon for the spacetime. This result still required imposing and exploiting axisymmetry in the calculations.
Some of the first documented attempts to solve the Einstein equations in three dimensions were focused on a single Schwarzschild black hole, which is described by a static and spherically symmetric solution to the Einstein field equations. This provides an excellent test case in numerical relativity because it does have a closed-form solution so that numerical results can be compared to an exact solution, because it is static, and because it contains one of the most numerically challenging features of relativity theory, a physical singularity. One of the earliest groups to attempt to simulate this solution was Anninos et al. in 1995. In their paper they point out that
- "Progress in three dimensional numerical relativity has been impeded in part by lack of computers with sufficient memory and computational power to perform well resolved calculations of 3D spacetimes."
Maturation of the field
In the years that followed, not only did computers become more powerful, but also various research groups developed alternate techniques to improve the efficiency of the calculations. With respect to black hole simulations specifically, two techniques were devised to avoid problems associated with the existence of physical singularities in the solutions to the equations: (1) Excision, and (2) the "puncture" method. In addition the Lazarus group developed techniques for using early results from a short-lived simulation solving the nonlinear ADM equations, in order to provide initial data for a more stable code based on linearized equations derived from perturbation theory. More generally, adaptive mesh refinement techniques, already used in computational fluid dynamics were introduced to the field of numerical relativity.
In the excision technique, which was first proposed in the late 1990s, a portion of a spacetime inside of the event horizon surrounding the singularity of a black hole is simply not evolved. In theory this should not affect the solution to the equations outside of the event horizon because of the principle of causality and properties of the event horizon (i.e. nothing physical inside the black hole can influence any of the physics outside the horizon). Thus if one simply does not solve the equations inside the horizon one should still be able to obtain valid solutions outside. One "excises" the interior by imposing ingoing boundary conditions on a boundary surrounding the singularity but inside the horizon. While the implementation of excision has been very successful, the technique has two minor problems. The first is that one has to be careful about the coordinate conditions. While physical effects cannot propagate from inside to outside, coordinate effects could. For example, if the coordinate conditions were elliptical, coordinate changes inside could instantly propagate out through the horizon. This then means that one needs hyperbolic type coordinate conditions with characteristic velocities less than that of light for the propagation of coordinate effects (e.g., using harmonic coordinates coordinate conditions). The second problem is that as the black holes move, one must continually adjust the location of the excision region to move with the black hole.
The excision technique was developed over several years including the development of new gauge conditions that increased stability and work that demonstrated the ability of the excision regions to move through the computational grid. The first stable, long-term evolution of the orbit and merger of two black holes using this technique was published in 2005.
In the puncture method the solution is factored into an analytical part, which contains the singularity of the black hole, and a numerically constructed part, which is then singularity free. This is a generalization of the Brill-Lindquist  prescription for initial data of black holes at rest and can be generalized to the Bowen-York prescription for spinning and moving black hole initial data. Until 2005, all published usage of the puncture method required that the coordinate position of all punctures remain fixed during the course of the simulation. Of course black holes in proximity to each other will tend to move under the force of gravity, so the fact that the coordinate position of the puncture remained fixed meant that the coordinate systems themselves became "stretched" or "twisted," and this typically led to numerical instabilities at some stage of the simulation.
In 2005 researchers demonstrated for the first time the ability to allow punctures to move through the coordinate system, thus eliminating some of the earlier problems with the method. This allowed accurate long-term evolutions of black holes. By choosing appropriate coordinate conditions and making crude analytic assumption about the fields near the singularity (since no physical effects can propagate out of the black hole, the crudeness of the approximations does not matter), numerical solutions could be obtained to the problem of two black holes orbiting each other, as well as accurate computation of gravitational radiation (ripples in spacetime) emitted by them.
The Lazarus project (1998–2005) was developed as a post-Grand Challenge technique to extract astrophysical results from short lived full numerical simulations of binary black holes. It combined approximation techniques before (post-Newtonian trajectories) and after (perturbations of single black holes) with full numerical simulations attempting to solve General Relativity field equations. All previous attempts to numerically integrate in supercomputers the Hilbert-Einstein equations describing the gravitational field around binary black holes led to software failure before a single orbit was completed.
The Lazarus approach, in the meantime, gave the best insight into the binary black hole problem and produced numerous and relatively accurate results, such as the radiated energy and angular momentum emitted in the latest merging state, the linear momentum radiated by unequal mass holes, and the final mass and spin of the remnant black hole. The method also computed detailed gravitational waves emitted by the merger process and predicted that the collision of black holes is the most energetic single event in the Universe, releasing more energy in a fraction of a second in the form of gravitational radiation than an entire galaxy in its lifetime.
Adaptive mesh refinement
Adaptive mesh refinement (AMR) as a numerical method has roots that go well beyond its first application in the field of numerical relativity. Mesh refinement first appears in the numerical relativity literature in the 1980s, through the work of Choptuik in his studies of critical collapse of scalar fields. The original work was in one dimension, but it was subsequently extended to two dimensions. In two dimensions, AMR has also been applied to the study of inhomogeneous cosmologies, and to the study of Schwarzschild black holes. The technique has now become a standard tool in numerical relativity and has been used to study the merger of black holes and other compact objects in addition to the propagation of gravitational radiation generated by such astronomical events.
In the past few years, hundreds of research papers have been published leading to a wide spectrum of mathematical relativity, gravitational wave, and astrophysical results for the orbiting black hole problem. This technique extended to astrophysical binary systems involving neutron stars and black holes, and multiple black holes. One of the most surprising predictions is that the merger of two black holes can give the remnant hole a speed of up to 4000 km/s that can allow it to escape from any known galaxy. The simulations also predict an enormous release of gravitational energy in this merger process, amounting up to 8% of its total rest mass.
- "Archived copy". Archived from the original on 2006-07-12. Retrieved 2005-12-01.CS1 maint: archived copy as title (link)
- Einstein, Albert. Der Feldgleichungen der Gravitation. Sitzungsberiche der Deutschen Akademie der Wissenschaften zu Berlin, Klasse fur Mathematik, Physik, und Technik.
- Arnowitt, R.; Deser, S.; Misner, C. W. (1962). "The dynamics of general relativity". In Witten, L. (ed.). Gravitation: An Introduction to Current Research. New York: Wiley. pp. 227–265.
- Hahn, S. G.; Lindquist, R. W. (1964). "The two-body problem in geometrodynamics". Ann. Phys. 29 (2): 304–331. Bibcode:1964AnPhy..29..304H. doi:10.1016/0003-4916(64)90223-4.
- Smarr, Larry (1975). The Structure of General Relativity with a Numerical Example. Ph.D. Dissertation, University of Texas, Austin. Austin, Texas.
- Smarr, Larry (1977). "Spacetimes generated by computers: Black holes with gravitational radiation". Ann. N.Y. Acad. Sci. 302: 569–. Bibcode:1977NYASA.302..569S. doi:10.1111/j.1749-6632.1977.tb37076.x. S2CID 84665358.
- Eppley, K. (1975). The numerical evolution of the collision of two black holes. Ph.D. Dissertation, Princeton University. Princeton, New Jersey.
- Piran, T. (1978). "Cylindrical general relativistic collapse". Phys. Rev. Lett. 41 (16): 1085–1088. Bibcode:1978PhRvL..41.1085P. doi:10.1103/PhysRevLett.41.1085.
- Stark, R. F.; Piran, T. (1985). "Gravitational-wave emission from rotating gravitational collapse". Phys. Rev. Lett. 55 (8): 891–894. Bibcode:1985PhRvL..55..891S. doi:10.1103/PhysRevLett.55.891. PMID 10032474.
- Matzner, Richard A.; Seidel, H. E.; Shapiro, Stuart L.; Smarr, L.; Suen, W.-M.; Teukolsky, Saul A.; Winicour, J. (1995). "Geometry of a black hole collision" (PDF). Science. 270 (5238): 941–947. Bibcode:1995Sci...270..941M. doi:10.1126/science.270.5238.941. S2CID 121172545.
- Anninos, Peter; Camarda, Karen; Masso, Joan; Seidel, Edward; Suen, Wai-Mo; Towns, John (1995). "Three dimensional numerical relativity: the evoluation of black holes". Phys. Rev. D. 52 (4): 2059–2082. arXiv:gr-qc/9503025. Bibcode:1995PhRvD..52.2059A. doi:10.1103/PhysRevD.52.2059. PMID 10019426. S2CID 15501717.
- Alcubierre, Miguel; Brugmann, Bernd (2001). "Simple excision of a black hole in 3+1 numerical relativity". Phys. Rev. D. 63 (10): 104006. arXiv:gr-qc/0008067. Bibcode:2001PhRvD..63j4006A. doi:10.1103/PhysRevD.63.104006. S2CID 35591865.
- Bona, C.; Masso, J.; Seidel, E.; Stela, J. (1995). "New formalism for numerical relativity". Phys. Rev. Lett. 75 (4): 600–603. arXiv:gr-qc/9412071. Bibcode:1995PhRvL..75..600B. doi:10.1103/PhysRevLett.75.600. PMID 10060068. S2CID 19846364.
- Cook, G. B.; et al. (1998). "Boosted three-dimensional black-hole evolutions with singularity excision". Phys. Rev. Lett. 80 (12): 2512–2516. arXiv:gr-qc/9711078. Bibcode:1998PhRvL..80.2512C. doi:10.1103/PhysRevLett.80.2512. S2CID 14432705.
- Alcubierre, Miguel (2003). "Hyperbolic slicings of spacetime: singularity avoidance and gauge shocks". Classical and Quantum Gravity. 20 (4): 607–623. arXiv:gr-qc/0210050. Bibcode:2003CQGra..20..607A. doi:10.1088/0264-9381/20/4/304. S2CID 119349361.
- Alcubierre, Miguel; Brugmann, Bernd; Diener, Peter; Koppitz, Michael; Pollney, Denis; Seidel, Edward; Takahashi, Ryoji (2003). "Gauge conditions for long-term numerical black hole evolutions without excision". Phys. Rev. D. 67 (8): 084023. arXiv:gr-qc/0206072. Bibcode:2003PhRvD..67h4023A. doi:10.1103/PhysRevD.67.084023. S2CID 29026273.
- Brugmann, Bernd; Tichy, Wolfgang; Jansen, Nina (2004). "Numerical simulation of orbiting black holes". Phys. Rev. Lett. 92 (21): 211101. arXiv:gr-qc/0312112. Bibcode:2004PhRvL..92u1101B. doi:10.1103/PhysRevLett.92.211101. PMID 15245270. S2CID 17256720.
- Shoemaker, Deirdre; Smith, Kenneth; Sperhake, Ulrich; Laguna, Pablo; Schnetter, Erik; Fiske, David (2003). "Moving black holes via singularity excision". Class. Quantum Grav. 20 (16): 3729–3744. arXiv:gr-qc/0301111. Bibcode:2003CQGra..20.3729S. doi:10.1088/0264-9381/20/16/313. S2CID 118897417.
- Pretorius, F. (2005). "Evolution of Binary Black-Hole Spacetimes". Phys. Rev. Lett. 95 (12): 121101. arXiv:gr-qc/0507014. Bibcode:2005PhRvL..95l1101P. doi:10.1103/PhysRevLett.95.121101. PMID 16197061. S2CID 24225193.
- Brandt, Steven; Bruegmann, Bernd (1997). "A simple construction of initial data for multiple black holes". Physical Review Letters. 78 (19): 3606–3609. arXiv:gr-qc/9703066. Bibcode:1997PhRvL..78.3606B. doi:10.1103/PhysRevLett.78.3606. S2CID 12024926.
- Brill, D.; Lindquist, R. (1963). "Interaction energy in geometrostatics". Phys. Rev. 131 (1): 471–476. Bibcode:1963PhRv..131..471B. doi:10.1103/PhysRev.131.471.
- Bowen, J.; York, J. W. (1980). "Time-asymmetric initial data for black holes and black-hole collisions". Phys. Rev. D. 21 (8): 2047–2056. Bibcode:1980PhRvD..21.2047B. doi:10.1103/PhysRevD.21.2047.
- Campanelli, M.; Lousto, C. O.; Marronetti, P.; Zlochower, Y. (2006). "Accurate Evolutions of Orbiting Black-Hole Binaries without Excision". Phys. Rev. Lett. 96 (11): 111101. arXiv:gr-qc/0511048. Bibcode:2006PhRvL..96k1101C. doi:10.1103/PhysRevLett.96.111101. PMID 16605808. S2CID 5954627.
- Baker, John G.; Centrella, Joan; Choi, Dae-Il; Koppitz, Michael; van Meter, James (2006). "Gravitational-Wave Extraction from an Inspiraling Configuration of Merging Black Holes". Phys. Rev. Lett. 96 (11): 111102. arXiv:gr-qc/0511103. Bibcode:2006PhRvL..96k1102B. doi:10.1103/PhysRevLett.96.111102. PMID 16605809. S2CID 23409406.
- Baker, J.; Campanelli, M.; Lousto, C. O. (2002). "The Lazarus project: A pragmatic approach to binary black hole evolutions". Phys. Rev. D. 65 (4): 044001. arXiv:gr-qc/0104063. Bibcode:2002PhRvD..65d4001B. doi:10.1103/PhysRevD.65.044001. S2CID 11080736.
- Baker, J.; Brügmann, B.; Campanelli, M.; Lousto, C. O.; Takahashi, R. (2001). "Plunge wave forms from inspiralling binary black holes". Phys. Rev. Lett. 87 (12): 121103. arXiv:gr-qc/0102037. Bibcode:2001PhRvL..87l1103B. doi:10.1103/PhysRevLett.87.121103. PMID 11580497. S2CID 39434471.
- Baker, J.; Campanelli, M.; Lousto, C. O.; Takahashi, R. (2002). "Modeling gravitational radiation from coalescing binary black holes". Phys. Rev. D. 65 (12): 124012. arXiv:astro-ph/0202469. Bibcode:2002PhRvD..65l4012B. doi:10.1103/PhysRevD.65.124012. S2CID 39834308.
- Campanelli, Manuela (2005). "Understanding the fate of merging supermassive black holes". Class. Quantum Grav. 22 (10): S387–S393. arXiv:astro-ph/0411744. Bibcode:2005CQGra..22S.387C. doi:10.1088/0264-9381/22/10/034. S2CID 119011566.
- Baker, J.; Campanelli, M.; Lousto, C. O.; Takahashi, R. (2004). "Coalescence remnant of spinning binary black holes". Phys. Rev. D. 69 (2): 027505. arXiv:astro-ph/0305287. Bibcode:2004PhRvD..69b7505B. doi:10.1103/PhysRevD.69.027505. S2CID 119371535.
- Choptuik, M. W. (1989). "Experiences with an adaptive mesh refinement algorithm in numerical relativity". In Evans, C.; Finn, L.; Hobill, D. (eds.). Frontiers in numerical relativity. Cambridge: Cambridge University Press. ISBN 0521366666.
- Choptuik, M. W. (1993). "Universality and scaling in gravitational collapse of massless scalar field". Phys. Rev. Lett. 70 (1): 9–12. Bibcode:1993PhRvL..70....9C. doi:10.1103/PhysRevLett.70.9. PMID 10053245.
- Choptuik, Matthew W.; Hirschmann, Eric W.; Liebling, Steven L.; Pretorius, Frans (2003). "Critical collapse of the massless scalar field in axisymmetry". Phys. Rev. D. 68 (4): 044007. arXiv:gr-qc/0305003. Bibcode:2003PhRvD..68d4007C. doi:10.1103/PhysRevD.68.044007. S2CID 14053692.
- Hern, Simon David (1999). Numerical relativity and inhomogeneous cosmologies. Ph.D. Dissertation, Cambridge University.
- Belanger, Z. B. (2001). Adaptive mesh refinement in the T2 symmetric spacetime. Master's Thesis, Oakland University.
- Schnetter, Erik; Hawley, Scott H.; Hawke, Ian (2004). "Evolutions in 3D numerical relativity using fixed mesh refinement". Class. Quantum Grav. 21 (6): 1465–1488. arXiv:gr-qc/0310042. Bibcode:2004CQGra..21.1465S. doi:10.1088/0264-9381/21/6/014. S2CID 52322605.
- Imbiriba, Breno; Baker, John; Choi, Dae-Il; Centrella, Joan; Fiske, David R.; Brown, J. David; van Meter, James R.; Olson, Kevin (2004). "Evolving a puncture black hole with fixed mesh refinement". Phys. Rev. D. 70 (12): 124025. arXiv:gr-qc/0403048. Bibcode:2004PhRvD..70l4025I. doi:10.1103/PhysRevD.70.124025. S2CID 119376660.
- Fiske, David R.; Baker, John G.; van Meter, James R.; Choi, Dae-Il; Centrella, Joan M. (2005). "Wave zone extraction of gravitational radiation in three-dimensional numerical relativity". Phys. Rev. D. 71 (10): 104036. arXiv:gr-qc/0503100. Bibcode:2005PhRvD..71j4036F. doi:10.1103/PhysRevD.71.104036. S2CID 119402841.
- Etienne, Zachariah B.; Liu, Yuk Tung; Shapiro, Stuart L.; Baumgarte, Thomas W. (2009). "Relativistic Simulations of Black Hole-Neutron Star Mergers: Effects of black-hole spin". Phys. Rev. D. 76 (4): 104021. arXiv:0812.2245. Bibcode:2009PhRvD..79d4024E. doi:10.1103/PhysRevD.79.044024. S2CID 119110932.
- Lousto, Carlos O.; Zlochower, Yosef (2008). "Foundations of multiple-black-hole evolutions". Phys. Rev. D. 77 (2): 024034. arXiv:0711.1165. Bibcode:2008PhRvD..77b4034L. doi:10.1103/PhysRevD.77.024034. S2CID 96426196.
- Campanelli, Manuela; Lousto, Carlos O.; Zlochower, Yosef; Merritt, David (2007). "Maximum Gravitational Recoil". Phys. Rev. Lett. 98 (23): 231102. arXiv:gr-qc/0702133. Bibcode:2007PhRvL..98w1102C. doi:10.1103/PhysRevLett.98.231102. PMID 17677894. S2CID 29246347.
- Healy, James; Herrmann, Frank; Hinder, Ian; Shoemaker, Deirdre M.; Laguna, Pablo; Matzner, Richard A. (2009). "Superkicks in Hyperbolic Encounters of Binary Black Holes". Phys. Rev. Lett. 102 (4): 041101. arXiv:0807.3292. Bibcode:2009PhRvL.102d1101H. doi:10.1103/PhysRevLett.102.041101. PMID 19257409. S2CID 9897187.
- Campanelli, Manuela; Lousto, Carlos O.; Zlochower, Yosef; Krishnan, Badri; Merritt, David (2007). "Spin flips and precession in black-hole-binary mergers". Phys. Rev. D. 75 (6): 064030. arXiv:gr-qc/0612076. Bibcode:2007PhRvD..75f4030C. doi:10.1103/PhysRevD.75.064030. S2CID 119334687.
- Initial Data for Numerical Relativity — A review article which includes a technical discussion of numerical relativity.
- Rotating Stars in Relativity — A technical review article about rotating stars, with a section on numerical relativity applications.
- A Relativity Tutorial at Caltech — A basic introduction to concepts of Numerical Relativity.