Cauchy Evolution and Hyperbolicity

Cauchy evolution in general relativity is a problem about more than rewriting Einstein’s equation with a time coordinate. A formulation must determine how physical, gauge, and constraint fields propagate, and its principal symbol must support a frequency-uniform estimate in a declared norm. On a finite domain the formulation must also control incoming fields, preserve the subsidiary constraints, and match initial data to boundary data at their corner. This chapter develops that chain of reasoning through ADM, BSSN, generalized harmonic, and Z4 systems. Its conclusions are continuum and formulation-specific: neither numerical success nor geometric equivalence on the Einstein constraint surface substitutes for a well-posedness proof.

Helpful background. The Einstein equation supplies the physical field equations and constraints. Hypersurface normals and tangential projectors supplies Cauchy and boundary decompositions. Partial differential equations and Green functions supplies symbols and energy estimates. Formulations of General Relativity supplies the volume-wide distinction between equivalent reductions. Differential equations reviews continuous dependence and homogeneous subsidiary evolution.

From an evolution system to a theorem

Begin with a complete choice of evolved variables, reduction order, constraint additions, and gauge equations. Freeze a smooth background and contract the spatial principal matrices with every real unit covector. Real eigenvalues are only the first test. Strong hyperbolicity requires a complete characteristic basis, or an equivalent positive symmetrizer, bounded uniformly as the direction varies. That uniformity removes wave-number growth from the Fourier propagator and produces the local Cauchy estimate required for continuous dependence.

The next step is interpretive. Characteristic fields must be separated into physical curvature combinations, gauge fields, and constraint or reduction fields. Equal coordinate speeds do not give these sectors equal physical meaning. A gauge driver may have a speed different from the metric null speed without transmitting a physical observable at that rate; a constraint mode may share a null speed while vanishing on every Einstein solution. Any use of “causal” must therefore state which sector and which metric cone are being compared.

Finally derive the subsidiary equations for all constraints. Their zero solution must be preserved, their own principal symbol must be controlled, and their incoming fields must receive compatible boundary conditions. An initial-boundary theorem also needs a boundary energy estimate and corner compatibility, not merely an interior characteristic list.

Reading path through the formulations

The pages are ordered so that each adds one necessary layer.

1. Hyperbolic PDE Systems in Relativity defines weak, strong, and symmetric hyperbolicity and proves the frequency-uniform estimate for a uniformly diagonalizable frozen symbol. Its Jordan-block example is the basic failure test.

2. Hyperbolicity Limitations of ADM Evolution applies that test to standard fixed-lapse, zero-shift ADM about Minkowski spacetime. A longitudinal–transverse sector has a defective zero-speed block, showing why equivalence to Einstein’s equation on constrained data does not guarantee well-posed free evolution.

3. The BSSN Continuum Formulation introduces the conformal metric, trace-free curvature, and promoted connection functions. It reconstructs the ADM fields and shows where algebraic, connection, Hamiltonian, and momentum constraints enter. Hyperbolicity remains conditional on the gauge-complete system.

4. Generalized Harmonic Evolution uses gauge sources to expose a quasilinear metric wave operator. Its homogeneous harmonic-constraint wave equation explains precisely when a reduced solution is an Einstein solution.

5. Z4 Systems and Constraint Damping promotes constraint violations to a covector and adds lower-order friction. The Fourier calculation proves decay of nonconstant linearized modes in the stated parameter range while displaying the undamped constant mode.

6. Characteristics, Constraint Propagation, and Well-Posedness assembles the full proof structure: complete symbols, sector identification, uniform estimates, subsidiary evolution, incoming-field count, boundary flux, and compatibility.

A common comparison language

Comparisons are meaningful only when the same questions are asked of every system. First identify the state vector and the norm: a second-order metric system, a first-order derivative reduction, and a conformal system do not have interchangeable norms. Next state the gauge. Prescribed lapse and shift, harmonic sources, slicing drivers, and shift drivers contribute different characteristic blocks. Then list every constraint and say whether it is imposed elliptically, propagated homogeneously, added to the main equations, promoted to a variable, or damped at lower order.

The frozen-background result must remain separate from the quasilinear theorem. A complete, uniformly bounded frozen symbol is the local algebraic condition. Extending it to variable coefficients requires a sufficiently regular symmetrizer and coefficient bounds in an appropriate Sobolev class. Extending from the Cauchy problem to a finite boundary requires an incoming convention, boundary operator, flux estimate, and compatibility hierarchy.

Constraint damping is likewise separate from hyperbolicity. A lower-order damping term can change growth rates without changing the principal symbol. It can make small nonconstant violations decay while leaving constant or long-wavelength modes weakly controlled. It cannot prevent a boundary from injecting violations, and it says nothing by itself about discrete stability.

Boundary and constraint closure

For a symmetric hyperbolic system, the normal principal matrix determines the energy flux and the number of incoming fields. Boundary conditions must prescribe exactly the incoming freedom, possibly coupling it contractively to outgoing fields, while leaving zero-speed fields consistent with the equations. Friedrich and Nagy 1999, §3, pp. 624–628 gives the maximally dissipative boundary-matrix construction; their later reduced and subsidiary systems show how a boundary-adapted Einstein formulation can recover the full vacuum equations.

Second-order harmonic evolution provides a complementary example. Kreiss, Reula, Sarbach, and Winicour 2007, §§2–3, Theorems 1–3, pp. 3–10, PDF derives interior and boundary derivative estimates under explicit conditions and imposes harmonic-constraint preservation. These examples use different variables and boundary mechanisms, but they share the essential sequence: a hyperbolic reduced system, a controlled boundary flux, compatible data, and a subsidiary argument returning to Einstein’s equation.

Corner conditions are consequential. At the intersection of the initial slice and timelike boundary, the boundary value and its time derivatives—computed from the evolution equations—must agree with the initial spatial derivatives to the regularity order claimed. A boundary condition that is dissipative for arbitrary traces may still fail to admit smooth solutions for incompatible initial data.

Review the chapter

A proposed formulation is advertised as follows: “Every characteristic speed is real, its constraints are damped, and all speeds have magnitude at most one; therefore the system is causal and well posed on a finite domain.” Identify the missing demonstrations and state what each would establish.

Solution

Real speeds establish at most weak hyperbolicity. One must exhibit a complete characteristic basis with a uniform condition-number bound, or a uniformly positive symmetrizer, for every real spatial direction. That supplies the frequency-independent interior estimate.

The formulation and gauge must be specified completely, including reduction variables and constraint additions. The fields must then be classified as physical, gauge, or constraint fields. A bound on coordinate speeds cannot establish physical causality until the physical sector is identified relative to the metric null cone.

“Constraints are damped” needs a derived subsidiary system, a background, parameter range, amplitude regime, and treatment of constant and long-wavelength modes. Homogeneous propagation or linear damping in the interior does not prevent incoming constraint violations at a boundary.

For a finite domain, one must define the outward normal, count incoming fields, impose the correct number of boundary conditions, and prove a boundary-stable or maximally dissipative energy estimate. The induced boundary conditions on incoming subsidiary fields must preserve the constraints. Initial and boundary data must satisfy the required corner compatibility conditions. Only after these steps, together with coefficient regularity and a declared norm, is a local well-posedness statement justified. Global existence, discrete stability, and formulation-independent geometric uniqueness remain separate questions.

References

• Friedrich, Helmut, and Gabriel Nagy. “The Initial Boundary Value Problem for Einstein’s Vacuum Field Equation.” Communications in Mathematical Physics 201 (1999): 619–655. doi:10.1007/s002200050571.
• Kreiss, Heinz-Otto, Oscar Reula, Olivier Sarbach, and Jeffrey Winicour. “Well-Posed Initial-Boundary Value Problem for the Harmonic Einstein Equations Using Energy Estimates.” Classical and Quantum Gravity 24 (2007): 5973–5984. arXiv:0707.4188.