---
address: 'Massachusetts Institute of Technology, Cambridge, MA 02139'
author:
- 'Maurice H.P.M. van Putten'
title: 'Nonlinear wave equations for numerical relativity: towards the computation of gravitational wave forms of black hole binaries'
---

\#1\#2\#3\#4[[\#1]{} [**\#2**]{}, \#3 (\#4)]{}

Introduction
============

The advanced stages of construction of LIGO and VIRGO has stimulated some new developments in hyperbolic formulations for numerical relativity[@mvp1; @Abr]. These developments describe nonlinear gravitational wave motion in the late stages of spiral in of compact binaries of black holes and neutron stars. Binary black hole and neutron star systems remain today the prime candidates for astrophysical sources of gravitational radiation, and statistical calculations indicate their event rates to be more than a few per year at the level of the sensitivity of the initial LIGO [@ph; @na]. Other known burst sources with potentially detectable levels of gravitational radiation are the collapse of white dwarfs into neutron stars[@usov], and secularly unstable single neutron stars[@lai].

$\dot{f}(f)-$diagram for burst sources
======================================

The time evolution of the above mentioned three classes of burst sources of gravitational radiation can be illustrated (at the present level of our understanding) in an $\dot{f}(f)-$diagram (Figure 1). Recall that compact binaries show an initial spiral in set by the Newtonian relationship $$\begin{aligned}
\dot{f}\sim f^{\frac{11}{3}},\end{aligned}$$ where $\dot{f}=\mbox{d}f/\mbox{d}t$. The frequency $f$ refers to the base frequency of the gravitational wave-form, initially of the “chirp," and finally of the quasi-normal mode ring down. During the intermediate merger phase, $f$ and $\dot{f}$ either follow by interpolation, or from $f=\frac{1}{\pi}\dot{\phi}$ with $\phi$ the orbital phase. Numerical relativity on black hole binaries specifically targets the merger phase, where post-Newtonian calculations break down. The present challenges are found to be the computation of stable, multiple orbits and the treatment of horizon boundary conditions.

[**FIGURE 1.**]{} [ Three cartoons of burst sources of gravitational waves in an $\dot{f}(f)$-diagram. The [**top**]{} diagram shows the three phases of binary coalescence: the chirp (with the Newtonian result $\dot{f}(f)\sim
f^{11/3}$), the merger phase, and the ring down phase with quasi-normal mode frequency $f_{QNR}\approx865 \frac{37M}{M}$ Hz (see, $e.g.$, Flanagan & Hughes[@fla]). The merger phase of two black holes and two neutron stars is expected to be distinct, but in as yet unknown ways. Sketched is the expected trend with increasing spin $S$ of the binary objects. Order of magnitude estimates[@fla; @tho; @sch] of the radiated energy in the merger phase is also indicated. (A possible break up and disk formation in neutron star mergers has been by Zhuge[@zhu], and could also be included here.) The [**middle**]{} diagram sketches the expected gravitational wave form resulting from a dynamical instability in an (accretion induced) collapse scenario, such as the collapse of a white dwarf to a neutron star. The arrow indicates the direction of time in this trajectory. The first stage shows a rapid sweep of increasing frequency, followed by a phase of rapid damping due to gravitational radiation. This scenario assumes the ratio $\beta$ of kinetic energy to gravitational energy to exceed the value which triggers a dynamical instability ($\beta>0.27$). It has been suggested as a model for GRBs by Usov[@usov]. The [**lower**]{} diagram sketches a translation of results on the secular instability of single neutron stars, $\beta_{\mbox{sec}}(=0.1375)
<\beta<\beta_{\mbox{dyn}}(=0.27)$, of Lai & Shapiro[@lai] into the $\dot{f}(f)-$diagram. The secular instability describes gravitational radiation driven transitions from unstable Jacobi ($J_u$) ellipsoids to either stable Mclaurin spheroids ($M_s$) or stable Dedekind ellipsoids ($D_s$). The Jacobi ellipsoids are rotation $(\Omega)$ dominated (rel. vorticity), and the Dedekind ellipsoids are vorticity ($\omega$) dominated (rel. rotation) stars. The trajectories are extended to the left, where the shape is in counter rotation relative to the fluid’s vorticity. Note that the end point of the trajectory in the middle diagram connects to $J_u$ of the lower diagram.]{}

Nonlinear wave equations for relativity
=======================================

Fully covariant, nonlinear wave equations have recently been obtained general relativity in the tetrad approach. They are Yang-Mills type equations for the Utiyama-connections $\omega_{a\mu\nu}=(e_\mu)_c\nabla_a(e_\nu)^c$ (for a tetrad field $\{(e_\mu)_a\}_{\mu=1}^4$) in the Lorentz gauge $\nabla^a\omega_{a\mu\nu}$. In vacuo, they reduce to $$\begin{aligned}
%\hat{\nabla}_c\hat{\nabla}^c
\hat{\Box}
\omega_{a\mu\nu}-
[\omega^c,\nabla_a\omega_{c}]_{\mu\nu}=0,\end{aligned}$$ were $\hat{\Box}=
\hat{\nabla}_c\hat{\nabla}^c$ is the $SO(3,1)$ gauged wave-operator. The full set of equations further involves the equations of structure, which can be integrated in the form of $$\begin{aligned}
(e_\mu)_a(\tau)=\Lambda_\mu^{\hskip.1in\nu}(\tau,\tau_0)(e_\mu)_c(\tau_0)
+\int_{\tau_0}^\tau \Lambda_\mu^{\hskip.1in\nu}(\tau,s)\hat{\partial}_b
N_\nu(s)\mbox{d}s,\end{aligned}$$ where the spatial coordinates have been suppressed. Here, the $$\begin{aligned}
N_\mu=(e_\mu)_\tau\end{aligned}$$ are the [*tetrad lapse functions*]{}, which are algebraically equivalent to the Hamiltonian lapse and shift functions. Strict hyperbolicity is maintained with arbitrary $N_\mu$, provided they define space-like foliations. There is full generality in this regard. By comparison, the otherwise closely related hyperbolic formulation of Abrahams $et$ $al.$[@Abr] shows some restrictions on the Hamiltonian lapse function for maintaining strict hyperbolicity. (This difference was left unnoticed in the comparison by Finn[@Finn]).

A numerical implementation is given in van Putten,[@mvp2] and applied to a nonlinear Gowdy-wave in a Gowdy $T^3$-cosmology. Stable and accurate results are obtained with proper second-order convergence by comparison with analytical solutions. No indication of any “coordinate shocks" as reported by Alcubierre & Masso [@Alc] have been observed in these numerical experiments.

References {#references .unnumbered}
==========

[99]{}

E.S. Phinney,. R. Narayan, T. Piran & A. Shemi, . K.S. Thorne, gr-qc/9506086 (1995). B.F. Shutz,. X. Zhuge, J.M. Centrella & L.W. MacMillan, gr-qc/9610039 (1996). D. Lai & S.L.Shapiro,. E.E. Flanagan & S.A. Hughes, gr-qc/9701039 (1997). V.V. Usov,. M.H.P.M. van Putten & D.M. Eardley, . M.H.P.M. van Putten, . L.S. Finn, gr-qc/9603004, (1996). A. Abrahams, A. Anderson, Y. Choquet-Bruhat & J.W. York, ; in these proceedings. M. Alcubierre, J. Masso, gr-qc/9709024, (1997).
