Schwarzschild Tests of the Wahlquist-Estabrook-Buchman-Bardeen Tetrad Formulation for Numerical Relativity

L. T. Buchman,J. M. Bardeen

DOI: https://doi.org/10.1103/PhysRevD.72.124014

2005-12-07

Abstract:A first order symmetric hyperbolic tetrad formulation of the Einstein equations developed by Estabrook and Wahlquist and put into a form suitable for numerical relativity by Buchman and Bardeen (the WEBB formulation) is adapted to explicit spherical symmetry and tested for accuracy and stability in the evolution of spherically symmetric black holes (the Schwarzschild geometry). The lapse and shift which specify the evolution of the coordinates relative to the tetrad congruence are reset at frequent time intervals to keep the constant-time hypersurfaces nearly orthogonal to the tetrad congruence and the spatial coordinate satisfying a kind of minimal rate of strain condition. By arranging through initial conditions that the constant-time hypersurfaces are asymptotically hyperbolic, we simplify the boundary value problem and improve stability of the evolution. Results are obtained for both tetrad gauges (``Nester'' and ``Lorentz'') of the WEBB formalism using finite difference numerical methods. We are able to obtain stable unconstrained evolution with the Nester gauge for certain initial conditions, but not with the Lorentz gauge.

General Relativity and Quantum Cosmology

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is, in numerical relativity, to use the Wahlquist - Estabrook - Buchman - Bardeen (WEBB) four - vector formula to simulate and test the evolutionary stability of spherically symmetric black holes (i.e., Schwarzschild geometry). Specifically, the author hopes to verify the feasibility and accuracy of this four - vector formula in numerical relativity and explore the influence of different gauge conditions (such as Nester gauge and Lorentz gauge) on the evolution results.

### Main research contents:
1. **Application of the WEBB four - vector formula**:
- The author applies the WEBB four - vector formula to the spherically symmetric vacuum black - hole spacetime to test its accuracy and stability.
- The WEBB formula is a first - order symmetric hyperbolic system, which is suitable for numerical relativity.

2. **Coordinate selection and resetting**:
- By frequently resetting the lapse function and shift vector, keep the constant - time hypersurface close to being orthogonal to the four - vector homology lines and make the spatial coordinates satisfy the minimum - strain condition.
- The initial conditions are arranged so that the constant - time hypersurface is asymptotically hyperbolic, simplifying the boundary - value problem and improving the stability of the evolution.

3. **Comparison of different gauge conditions**:
- Using the finite - difference numerical method, tests are carried out under two different four - vector gauges (Nester gauge and Lorentz gauge).
- It is found that for some initial conditions, a stable unconstrained evolution can be obtained using the Nester gauge, while the Lorentz gauge cannot do this.

4. **Boundary - condition handling**:
- The inner boundary is located inside the event horizon, and a constant R value is maintained through the shift condition.
- The outer boundary is set at approximately \(R = 20M\), ensuring that variables such as \(\hat{n}_r, K_T, K_R\) reach approximately constant values.

### Summary of mathematical formulas:
- **Spacetime metric**:
\[
ds^2=\left[-\alpha^2+\beta_r^2 e^{2\lambda}(1 - \hat{A}_r^2)+2e^\lambda\alpha\beta_r\hat{A}_r\right]dt^2
+ 2e^\lambda\left[\alpha\hat{A}_r+\beta_r e^\lambda(1 - \hat{A}_r^2)\right]dr dt
+ e^{2\lambda}(1 - \hat{A}_r^2)dr^2+R^2d\theta^2+R^2\sin^2\theta d\phi^2
\]

- **Spatial metric**:
\[
dl^2 = e^{2\lambda}(1 - \hat{A}_r^2)dr^2+R^2d\theta^2+R^2\sin^2\theta d\phi^2
\]

- **Characteristic velocities**:
\[
s_1(r,t)=\frac{e^{-\lambda}\alpha}{1 + \hat{A}_r}-\beta_r
\]
\[
s_2(r,t)=-\frac{e^{-\lambda}\alpha}{1 - \hat{A}_r}-\beta_r
\]

- **Constraint equations**:
\[
B_R\partial_r K_T=(K_R - K_T)\left(\frac{1}{R}-\hat{n}_r\right)-\frac{\hat{A}_r}{2}\left[\frac{2}{R}(\hat{a}_r-\hat{n}_r)-3K_T^2 - 2\hat{a}_r\hat{n}_r+\hat{n}_r^2\right]
\]
\[
B_R\partial