PDEs and Green Functions for Gravitation

Relativistic field equations are meaningful only together with their principal part, admissible data, boundary conditions, and support properties. This bridge checks the limited PDE language needed to distinguish evolution from constraint problems and to use Green functions without hiding a boundary prescription. It does not replace a course in PDE theory or spectral analysis.

Required background. Calculus and Linear Algebra for Relativity and Differential Equations for Relativity supply derivatives, Fourier modes, initial and boundary data, and linear ODE methods. Helpful background. The Wave Equation and Characteristics develops wave propagation and geometric optics.

You should leave able to classify representative elliptic and hyperbolic problems, identify characteristics, and construct a retarded solution with its data and causal support explicit.

PDE and Green-function diagnostic

1. Classify the principal parts of $$u_{xx}+u_{yy}=f, \qquad u_{tt}-u_{xx}=f, \qquad u_t-u_{xx}=f.$$ State the natural kind of data for each.
2. Factor the $1+1$ wave operator $\partial_t^2-\partial_x^2$ and identify its characteristic curves.
3. Find the decaying Green function on the real line satisfying $$\left(\frac{d^2}{dx^2}-\kappa^2\right)G(x)=\delta(x), \qquad \kappa\gt0.$$ Derive its derivative jump rather than quoting it.
4. For $u_{tt}-u_{xx}=0$ with smooth compactly supported data, show that $$E(t)=\frac12\int_{\mathbb R}(u_t^2+u_x^2)\,dx$$ is conserved.

Hints and solutions

Hints for 1

Look only at the highest derivatives. Elliptic problems are boundary-value problems; hyperbolic equations evolve Cauchy data; the heat equation is parabolic and evolves one datum forward in time.

Solution to 1

$u_{xx}+u_{yy}=f$ is elliptic and normally receives boundary data on a spatial domain. $u_{tt}-u_{xx}=f$ is hyperbolic and receives $u$ and its normal time derivative on a noncharacteristic initial curve. $u_t-u_{xx}=f$ is parabolic and receives one initial profile plus spatial boundary or decay conditions. The classification concerns the principal symbol; lower-order terms do not usually change it, while degeneracy can.

Hints for 2

Use $(\partial_t-\partial_x)(\partial_t+\partial_x)$. Curves of constant $t-x$ and $t+x$ carry the two traveling families.

Solution to 2

$$\partial_t^2-\partial_x^2 =(\partial_t-\partial_x)(\partial_t+\partial_x).$$

The characteristics are $x-t=\text{constant}$ and $x+t=\text{constant}$, or $dx/dt=\pm1$. Data propagate along the null lines of the flat metric; the domain of dependence is bounded by them.

Hints for 3

Away from $x=0$, solve the homogeneous equation and impose decay at both infinities. Integrate the equation from $-\epsilon$ to $+\epsilon$ to determine the derivative jump.

Solution to 3

Decay and continuity give $G=C e^{-\kappa|x|}$. Integrating across the source yields

$$G'(0^+)-G'(0^-)=1.$$

Since the jump is $-2\kappa C$, $C=-1/(2\kappa)$. Thus

$$G(x)=-\frac{e^{-\kappa|x|}}{2\kappa}.$$

The exponential decay checks the chosen boundary condition; changing the domain or boundary data changes the Green function.

Hints for 4

Differentiate under the integral, replace $u_{tt}$ by $u_{xx}$, and integrate one spatial derivative by parts. Compact support removes the boundary flux.

Solution to 4

$$\frac{dE}{dt}=\int_{\mathbb R}(u_tu_{tt}+u_xu_{xt})\,dx =\int_{\mathbb R}(u_tu_{xx}+u_xu_{xt})\,dx.$$

Integrating the first term by parts gives

$$\frac{dE}{dt} =\left[u_tu_x\right]_{-\infty}^{+\infty} -\int u_{xt}u_x\,dx+\int u_xu_{xt}\,dx=0.$$

Without compact support or decay, the boundary flux must be retained. The identity supplies uniqueness and stability information, not merely a conserved number.

Interpret the result

• Secure: you distinguish PDE type and data, derive characteristics, normalize a Green function from its jump, and retain boundary flux in an energy estimate.
• Uncertain: the physical picture is right but one normalization or data statement needs a hint. Change the source or boundary condition and repeat that part.
• Not yet: elliptic and hyperbolic data are interchanged, a Green function is quoted without its operator and boundary condition, or causal support is inferred from a Fourier denominator alone.

Questions 1–2 primarily sample PDE and wave structure; questions 3–4 sample Green functions, boundary prescriptions, and energy control.

Focused review

Principal symbols and characteristics

For a second-order scalar operator $P=A^{\mu\nu}\partial_\mu\partial_\nu+\cdots$, the principal symbol is $A^{\mu\nu}\xi_\mu\xi_\nu$. Characteristic covectors make it vanish. Hyperbolic equations admit causal evolution from suitable noncharacteristic data; elliptic equations communicate boundary information throughout a spatial domain. Gauge systems may need constraints or gauge fixing before this classification becomes useful.

Distributions and Green operators

A Green function $G$ is an inverse only after the operator, domain, and boundary or support prescription are fixed. Distributional normalization is checked by integrating across the source or applying the operator to a test function. Retarded and advanced inverses solve the same differential equation but have different support. Review Retarded and Advanced Green Functions for the flat-space construction.

Energy estimates and causal support

Multiplying a hyperbolic equation by a suitable time derivative and integrating by parts produces an energy identity or inequality. Boundary flux, lower-order terms, and source work must be accounted for. On curved spacetime, global retarded propagation additionally needs geometric hypotheses; Advanced, Retarded, and Causal Propagators states the stronger setting.

Transfer task: retarded response of a field mode

Spatial Fourier transformation of a massive wave equation gives one mode

$$\ddot q_k(t)+\omega_k^2q_k(t)=f_k(t), \qquad \omega_k^2=k^2+m^2,$$

with $q_k(t_0)=\dot q_k(t_0)=0$. Construct the retarded solution and evaluate it for an impulse $f_k(t)=J\delta(t-t_*)$ with $t_*\gt t_0$. State the support and the derivative jump.

Hint 1

The retarded Green function of $d^2/dt^2+\omega_k^2$ vanishes for negative time difference and is proportional to $\sin(\omega_k\Delta t)$ for positive difference.

Hint 2

Integrate the mode equation across $t=t_*$. $q_k$ remains continuous while $\dot q_k$ jumps by $J$.

Solution

The retarded kernel is

$$G_k^{\rm ret}(t,t') =\Theta(t-t')\frac{\sin[\omega_k(t-t')]}{\omega_k}.$$

Therefore

$$q_k(t)=\int_{t_0}^{t}dt'\, \frac{\sin[\omega_k(t-t')]}{\omega_k}f_k(t').$$

For the impulse,

$$q_k(t)=J\,\Theta(t-t_*) \frac{\sin[\omega_k(t-t_*)]}{\omega_k}.$$

It vanishes before the source acts, is continuous at $t_*$, and satisfies

$$\dot q_k(t_*^+)-\dot q_k(t_*^-)=J.$$

The jump fixes normalization independently. A single Fourier mode is spatially delocalized, so temporal retardedness of $G_k^{\rm ret}$ is not by itself the full spacetime support statement; that emerges only after inverse transforming all modes with the correct contour and geometry.

Common errors

• Classifying an equation by all its terms instead of its principal symbol.
• Writing “the Green function” without naming the operator, domain, and boundary or causal prescription.
• Losing a sign or normalization because the delta-function jump was not checked.
• Dropping a boundary flux and then claiming conservation on a domain with an open boundary.

Return to your route

After the transfer task, return to the route that sent you here: Black-Hole Perturbation Theory, Computational Relativity Self-Study, Relativistic Cosmology, Quantum Fields in Curved Spacetime, Gravitational Radiation and Compact Binaries, or the Computational Relativity course.