Damped Harmonic Gauge in Numerical Relativity

• Damped harmonic gauge is a coordinate condition that dynamically damps gauge-induced growth and coordinate distortions, ensuring well-posed evolution of Einstein’s equations.
• It employs dynamic gauge-source functions and a first-order gauge-driver system to maintain symmetric hyperbolicity and robust numerical performance.
• The method enables construction of accurate binary black hole initial data using spectral elliptic solvers and dual-frame evolution techniques.

The damped harmonic gauge is a class of coordinate conditions for the generalized harmonic (GH) formulation of Einstein’s equations, central to stable and accurate numerical relativity simulations involving highly dynamical spacetimes such as binary black hole (BBH) mergers. By promoting the gauge-source functions to dynamical fields and prescribing their evolution toward suitable damping targets, the damped harmonic gauge effectively suppresses pathological gauge dynamics and ensures well-posedness (symmetric hyperbolicity) of the combined Einstein–gauge-driver system. The approach has been foundational in enabling robust evolutions, especially in multiframe (dual-frame) computational strategies and in the construction of BBH initial data directly in damped harmonic coordinates (0904.4873, Varma et al., 2018).

1. Formulation within the Generalized Harmonic Einstein System

In the generalized harmonic formalism, the four spacetime coordinates $x^a$ are cast as inhomogeneous wave equations with gauge-source functions $H^a$: $$H^a = \nabla^c \nabla_c x^a, \qquad H_a = \psi_{ab} H^b,$$ where $\psi_{ab}$ is the spacetime metric. Einstein’s equations are then written as a quasilinear system for the metric components: $$\psi^{cd} \partial_c \partial_d \psi_{ab} + \partial_a H_b + \partial_b H_a = Q_{ab}(H, \psi, \partial \psi),$$ with $Q_{ab}$ collecting terms quadratic in the metric, its derivatives, and $H^a$. As long as the gauge-source functions $H_a$ are specified as functions of the coordinates and metric, the evolution system is manifestly hyperbolic.

2. Damped Harmonic Gauge Condition

The damped harmonic (or damped-wave) gauge introduces a class of target gauge-source functions designed to dynamically damp gauge-induced growth in the lapse function and coordinate distortions in the shift vector. The target is

$$F_a = \mu_L \ln\left( \frac{g^p}{N} \right) t_a - \mu_S N^{-1} g_{ai} N^i,$$

where $N$ is the lapse, $H^a$0 the shift, $H^a$1 is the 3-metric with determinant $H^a$2, $H^a$3 is the unit normal to slices of constant time, and $H^a$4 are positive damping parameters (typically with $H^a$5). This yields a damped-wave evolution for the lapse and controls spatial coordinate drift.

In practice, for BBH simulations, these damping parameters may be chosen as functions of spatial radius and time to smoothly activate the gauge (e.g., $H^a$6, with $H^a$7 itself modulated in space and time) (Varma et al., 2018).

3. First-Order Gauge-Driver System and Hyperbolicity

To enforce $H^a$8 even when $H^a$9 depends on derivatives of the metric, $$H^a = \nabla^c \nabla_c x^a, \qquad H_a = \psi_{ab} H^b,$$0 is promoted to a dynamical field satisfying a first-order evolution system: $$H^a = \nabla^c \nabla_c x^a, \qquad H_a = \psi_{ab} H^b,$$1

$$H^a = \nabla^c \nabla_c x^a, \qquad H_a = \psi_{ab} H^b,$$2

where $$H^a = \nabla^c \nabla_c x^a, \qquad H_a = \psi_{ab} H^b,$$3 is an auxiliary field, and $$H^a = \nabla^c \nabla_c x^a, \qquad H_a = \psi_{ab} H^b,$$4 are positive driver parameters. In stationary regimes, $$H^a = \nabla^c \nabla_c x^a, \qquad H_a = \psi_{ab} H^b,$$5. The full Einstein–gauge-driver system, when reduced to first-order, is provably symmetric hyperbolic and thus well-posed (0904.4873).

4. Boundary Conditions

Well-posedness and constraint preservation require carefully constructed boundary conditions. The characteristic decomposition of the coupled system yields:

• Gauge sector: On $$H^a = \nabla^c \nabla_c x^a, \qquad H_a = \psi_{ab} H^b,$$6, a “Bjørhus driver” form is imposed at the boundary, $$H^a = \nabla^c \nabla_c x^a, \qquad H_a = \psi_{ab} H^b,$$7, ensuring the gauge target is enforced at boundaries.
• Metric sector: Constraint-preserving conditions project out incoming constraint-violating modes, while physical (radiation) conditions enforce that the incoming Newman–Penrose Weyl field vanishes, corresponding to no incoming gravitational radiation.

5. Construction of Damped Harmonic Black Hole Solutions

For initial data construction, especially for BBH systems, analytic single BH solutions are needed in damped harmonic coordinates. The method proceeds as follows (Varma et al., 2018):

• A coordinate transformation from Kerr–Schild coordinates to damped harmonic coordinates is constructed.
• Correction functions $$H^a = \nabla^c \nabla_c x^a, \qquad H_a = \psi_{ab} H^b,$$8 are introduced so that the transformed coordinates are stationary, accommodate specified boosts and spins, and result in the metric satisfying $$H^a = \nabla^c \nabla_c x^a, \qquad H_a = \psi_{ab} H^b,$$9.
• This yields a set of four coupled, nonlinear elliptic PDEs for $\psi_{ab}$0, solved with boundary conditions: $\psi_{ab}$1 (asymptotic flatness) and Robin-type regularity at the horizon.
• Spectral elliptic solvers with domain decomposition provide efficient, exponentially convergent solutions for these corrections, ensuring the constraint norm

$\psi_{ab}$2

is minimized below a specified tolerance.

These analytic single-damped-harmonic-BH solutions serve as building blocks for BBH initial data by superposing two such solutions (each with their own spin and boost), extracting their induced 3-metrics and extrinsic curvatures, and using the sum as conformal free data in the extended conformal thin sandwich (XCTS) constraint equations.

6. Implementation in Dual-Frame Evolution and Binary Black Hole Initial Data

Dual-frame evolution is essential in BBH simulations, with both an inertial frame (asymptotically non-rotating) and a co-moving frame (tracking BHs). The gauge-driver system accommodates either:

• Enforcing the gauge in inertial or co-moving frame,
• Or using a hybrid driver (with a radial weighting function) interpolating between the two.

With the construction of single BHs in damped harmonic gauge, one can now solve directly for quasi-equilibrium BBH initial data in which both black holes are in damped harmonic coordinates ab initio, eliminating gauge transients from early-time evolution. The approach uses spectral AMR and the analytic solution of the elliptic transformation equations, providing high accuracy and regularity throughout the computational domain (Varma et al., 2018).

7. Numerical Performance and Stability

Numerical tests in single BH evolutions demonstrate that the damped harmonic gauge:

• Rapidly suppresses gauge dynamics on timescales $\psi_{ab}$3,
• Maintains GH constraint violations $\psi_{ab}$4 in $\psi_{ab}$5 norm, with exponential convergence under increased resolution,
• Controls the gauge-mismatch norm $\psi_{ab}$6 to $\psi_{ab}$7 at late times,
• Remains stable under a wide range of driver parameters and perturbations, both in single- and dual-frame scenarios (0904.4873).

A plausible implication is that the robust control of gauge dynamics by the damped harmonic gauge is critical for long, stable evolutions of inspiraling and merging black hole binaries—with well-posed initial-boundary value structure, robust performance under grid refinement, and the elimination of artificial gauge transition transients in initial data setups.