Alphabeta Math
Pipeline-generated
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

Wave Equation Representation Formulas

1 · Prerequisites

2 · Summary

This page develops the classical representation formulas for the Cauchy problem of the wave equation □cu=∂t2u−c2Δu=0 in Rn. After fixing the vocabulary of wave operator, Cauchy data and wave speed, together with the unit-speed rescaling convention that transfers the unit-speed formulas of the sources to general speed, the page builds the one-dimensional theory: factorisation of the operator, the general solution as a sum of a right- and a left-travelling wave, d'Alembert's formula with uniqueness and data attainment, the domain of dependence, the forced formula over the characteristic triangle. The companion examples page contains a counterexample showing that data on one characteristic line do not determine the solution. In parallel it develops spherical means: their smooth signed-radius extension and parity, the vanishing first moment of the sphere measure, the ball–sphere radial identity, the Euler–Poisson–Darboux equation, the projection of sphere integrals of cylindrical functions onto weighted ball integrals, and the two iterated radial-derivative lemmas behind the odd-dimensional reduction. These tools yield Kirchhoff's formula in three dimensions, Poisson's formula in two dimensions by descent, the odd-dimensional formula by iterated spherical means, the even-dimensional formula by descent, the attainment of the Cauchy data, the support dichotomy between sphere-supported odd-dimensional kernels and interior even-dimensional kernels, time-reversal invariance, Duhamel's principle with the forced three-dimensional retarded potential, local determination by the Cauchy data, and a remark separating the wave Poisson formula from the harmonic Poisson kernel.

All statements of the page are read under the Axiom of Countable Choice, which supplies the polar surface measure and the sphere integrals; the one-dimensional construction is choice-free. Uniqueness of classical solutions in dimensions n≥2 is deliberately deferred to the energy page that consumes this one.

For odd dimensions, sphere dependence means dependence on the data and finitely many transverse derivatives along the sphere; agreement on an open neighbourhood suffices, while bare equality of values on the sphere generally does not. Duhamel's principle assumes joint continuity of the source's spatial derivatives through order ⌊n/2⌋+1, without requiring time derivatives.

3 · Logical flowchart

4 · Definitions, theorems and proofs

DefinitionDefinition: Literature-sourcedProof: Not applicablejudge pass (gpt-6.1-sol)Open item page →

Wave equation, Cauchy data and wave speed

Definition

Let n≥1, T>0 and c>0, let ∂t,∂j be the partial derivatives of Directional derivatives and partial derivatives of a map U⊆Rm→Rn and let Δ=∑j=1n∂j2 be the Laplacian of The Laplacian of a C2 function and of a C2 vector field. The wave operator of speed c is □c:=∂t2−c2Δ, and the wave equation with source f is the equation □cu=f on a slab Rn×(0,T); the equation is homogeneous when f=0. A classical solution on the time slab Rn×[0,T) is a function u of class C2 on Rn×(0,T) in the sense of Ck maps and multi-index derivative notation in Euclidean space which satisfies the equation at every point of Rn×(0,T). This is the classical-solution and Cauchy-data vocabulary of Scalar partial differential equations, order, and classical solutions, and □c is a linear second-order operator in the sense of Linear, semilinear, quasilinear, and fully nonlinear partial differential equations.

The Cauchy problem for □cu=f prescribes a displacement u0:Rn→R and a velocity u1:Rn→R and asks for a classical solution attaining them as t↓0 in the pointwise sense: u(x,t)→u0(x) and ∂tu(x,t)→u1(x) for every x∈Rn as t↓0. Only the limits are part of the data, so u need not be defined at t=0 a priori.

Unit-speed rescaling convention. A function u solves □cu=0 on Rn×[0,T) with data (u0,u1) if and only if v(x,s):=u(x,s/c) solves ∂s2v=Δv on Rn×[0,cT) with data (u0,u1/c); equivalently u(x,t)=v(x,ct). Indeed continuous coordinate partial derivatives imply total differentiability (If all partial derivatives exist on a neighbourhood and are continuous at a point, then the map is totally differentiable there with Jacobian derivative), so the chain rule (The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)) gives ∂tu(x,t)=c ∂sv(x,ct) and ∂t2u(x,t)=c2∂s2v(x,ct) while spatial derivatives are unchanged, so □cu=c2(∂s2v−Δv) evaluated at s=ct, and the initial displacement limits agree, while the initial velocity limit of u is c times that of v. The cited treatments state their formulas at unit speed, and this convention is the only place where the general speed enters those formulas.

Sphere normalisation. Let σ be the polar surface measure on the unit sphere (The polar surface set function on the unit sphere). Write Vm:=∣B1m∣ for the unit ball and ωm:=σ(Sm)=(m+1)Vm+1>0 for the total polar measure of the unit sphere Sm⊆Rm+1 (Sphere and ball measures scale in Rn); thus the boundary sphere of a ball in Rn has total measure ωn−1=nVn. For an integer m≥0 put (2m+1)!!:=(2m+1)(2m−1)⋯3⋅1 and (2m)!!:=(2m)(2m−2)⋯2, so that 0!!=1!!=1 and (n−2)!! is defined for every odd n≥3. All sphere and weighted-ball integrals on this page are read under the Axiom of Countable Choice of The Axiom of Countable Choice (ACω), under which the polar measure and its integrals are supplied.

LemmaStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

The iterated radial-derivative identity behind the odd-dimensional reduction

Statement

Let k≥1 and let Dr:=r−1∂r be the radial derivative on (0,∞). For every f∈Ck+1((0,∞)), ∂2∂r2Drk−1(r2k−1f(r))=Drk−1[r2k−11r2k∂∂r(r2k∂rf(r))] on (0,∞). Both sides are finite combinations of derivatives of f computed by the product and quotient rules; no integral and no differential equation for f is used.

Facts & Assumptions

Given: an integer k≥1, a function f∈Ck+1((0,∞)), and the operator Dr=r−1∂r acting on functions on (0,∞).

[F1]

Sums, products, quotients of differentiable functions are differentiable, with the usual sum, product, quotient rules; nonnegative integer powers are differentiated by repeated product rules, and reciprocals by the quotient rule on nonzero domains (Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0).

Proof

1.1F1given

Put g:=r2k−1f, so that g∈Ck+1((0,∞)) and f=r1−2kg. The left-hand side is ∂r2Drk−1g, and for the right-hand side the product rule gives r2k∂rf=(1−2k)g+rg′, hence r−1∂r(r2k∂rf)=g′′+(2−2k)r−1g′; since r2k−1r−2k=r−1, the right-hand side is Drk−1(g′′+(2−2k)r−1g′). It therefore suffices to prove ∂r2Drk−1g=Drk−1(g′′+(2−2k)r−1g′) for every g∈Ck+1((0,∞)), which is what the following steps do.

1.2F1algebra

On (0,∞) the operators satisfy ∂r=rDr, hence ∂r2=∂r(rDr)=Dr+r2Dr2; and for every j≥1 and every v∈Cj((0,∞)) one has Drj(r2v)=r2Drjv+2j Drj−1v. The last identity is proved by induction on j: for j=1, Dr(r2v)=r−1(2rv+r2v′)=2v+r2Drv; and if it holds for j, then Drj+1(r2v)=Dr(r2Drjv+2j Drj−1v)=r2Drj+1v+2Drjv+2j Drjv=r2Drj+1v+2(j+1)Drjv.

1.3F1algebra

Fix g∈Ck+1((0,∞)). If k=1, then Drk−1(r2Dr2g)=r2Dr2g=r2Drk+1g+2(k−1)Drkg directly; if k≥2, the second identity of the previous step with j=k−1 and v=Dr2g gives the same equality. Also g′′=∂r2g=Drg+r2Dr2g by the first identity of the previous step, so Drk−1g′′=Drkg+Drk−1(r2Dr2g)=r2Drk+1g+(2k−1)Drkg. Moreover r−1g′=Drg, so Drk−1((2−2k)r−1g′)=(2−2k)Drkg.

2.1F1algebra∎

Adding the two pieces of the previous step gives Drk−1(g′′+(2−2k)r−1g′)=r2Drk+1g+(2k−1+2−2k)Drkg=r2Drk+1g+Drkg, and by the first identity of the second step this last quantity is ∂r2Drk−1g. This proves the equivalent identity for g and hence, by the substitution of the first step, the identity of the statement.

LemmaStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

Radial-derivative expansion of the Euler–Poisson–Darboux transform and its zero-radius limit

Statement

Let k≥1 and let Dr=r−1∂r on (0,∞). For every f∈Ck−1((0,∞)) there are real constants αk,j, 0≤j≤k−1, such that Drk−1(r2k−1f(r))=∑j=0k−1αk,j rj+1f(j)(r),αk,0=(2k−1)!!. Consequently, if the finite limit f(0):=lim⁡r↓0f(r) exists and each derivative f(j) with 1≤j≤k−1 is bounded on (0,1] (an empty condition when k=1) — in particular if f extends to a Ck−1 function on [0,∞) — then r−1Drk−1(r2k−1f)(r)⟶(2k−1)!! f(0)(r↓0). Every coefficient αk,0=(2k−1)!! is nonzero, so the transformed average recovers the value at r=0 with the exact dimensional constant. The limit assumption is needed even for k=1, when f(r)=sin⁡(1/r) is bounded but has no limit. Boundedness of the higher derivatives is also substantive: for k=2 the function f(r)=rsin⁡(r−2), which is continuous at 0, has rf′(r)=rsin⁡(r−2)−2r−1cos⁡(r−2) unbounded, and then r−1Dr(r3f) does not tend to 3f(0).

Facts & Assumptions

Given: an integer k≥1, a function f∈Ck−1((0,∞)), and the operator Dr=r−1∂r on (0,∞).

Proof

1.1F1given

Base case. For k=1 and f∈C0((0,∞)) one has Dr0(r2k−1f)=rf=α1,0 r1f(0) with α1,0=1=(2⋅1−1)!!; under the limit hypothesis, r−1Dr0(rf)=f(r)→f(0) as r↓0.

1.2given

Induction hypothesis. Fix k≥1 and assume that for every h∈Ck−1((0,∞)) there are real constants αk,0,…,αk,k−1 such that Drk−1(r2k−1h)=∑j=0k−1αk,jrj+1h(j) and αk,0=(2k−1)!!.

1.3F1algebra

Shape of the (k+1)-st transform. Let f∈Ck((0,∞)) and put h:=r2f, so that h∈Ck((0,∞))⊆Ck−1((0,∞)) and r2k+1f=r2k−1h. The induction hypothesis gives Drk−1(r2k+1f)=∑j=0k−1αk,jrj+1h(j), and the product rule gives h(j)=r2f(j)+2jrf(j−1)+j(j−1)f(j−2) for 0≤j≤k−1, where the last two terms are read as 0 for j=0 and j=1.

1.4F1algebra

Applying Dr once to the finitely many resulting terms, and using Dr(rmv)=mrm−2v+rm−1v′, produces a finite sum ∑i=0kβiri+1f(i) with constants βi independent of f: the term αk,jrj+1⋅r2f(j) contributes (j+3)αk,jrj+1f(j) and αk,jrj+2f(j+1); the term 2jαk,jrj+2f(j−1) contributes 2j(j+2)αk,jrjf(j−1) and 2jαk,jrj+1f(j); and the term j(j−1)αk,jrj+1f(j−2) contributes j(j−1)(j+1)αk,jrj−1f(j−2) and j(j−1)αk,jrjf(j−1). Every displayed monomial rmf(i) has m=i+1 with 0≤i≤k, and negative powers do not occur because the terms with j=0,1 have one or both of the last two summands read as zero.

1.5F1algebra

Leading coefficient. Evaluating the identity of the previous step at the constant function f≡1, which lies in Ck, gives β0r=Drk(r2k+1)=[(2k+1)(2k−1)⋯3]r=(2k+1)!! r: indeed Dr(r2k+1)=(2k+1)r2k−1, each further application of Dr lowers the exponent by 2 and multiplies the coefficient by the previous exponent, and k applications leave the exponent 1. Hence β0=(2k+1)!!, and setting αk+1,i:=βi for 0≤i≤k extends the conclusion of the induction hypothesis from k to k+1.

2.1givenalgebra∎

Conclusion. By the base case and the induction step, the expansion Drk−1(r2k−1f)=∑j=0k−1αk,jrj+1f(j) with αk,0=(2k−1)!! holds for every k≥1 and every f∈Ck−1((0,∞)). If f(0):=lim⁡r↓0f(r) is finite and each f(j) with 1≤j≤k−1 is bounded on (0,1], then dividing by r gives r−1Drk−1(r2k−1f)(r)=∑j=0k−1αk,jrjf(j)(r)→αk,0f(0)=(2k−1)!!f(0) as r↓0: the j=0 term converges by the assumed limit, and each term with j≥1 is bounded by ∣αk,j∣rjsup⁡(0,1]∣f(j)∣→0.

LemmaStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

Reflection invariance and vanishing first moment of the sphere measure

Statement

Assume the Axiom of Countable Choice. Let n≥1 and let σ be the polar surface measure on the unit sphere Sn−1⊆Rn (The polar surface set function on the unit sphere). Then the reflection Rω=−ω preserves σ, and the first moment vanishes: ∫Sn−1ω dσ(ω)=0∈Rn,so∫Sn−1a⋅ω dσ(ω)=0(a∈Rn). The integrals are finite because σ is a finite Borel measure.

Facts & Assumptions

Given: the Axiom of Countable Choice, an integer n≥1 and the polar surface measure σ on Sn−1.

[F1]

σ(E)=nλn{rω:ω∈E, 0<r≤1} for Borel E⊆Sn−1 (The polar surface set function on the unit sphere).

[F2]

For an invertible linear T:Rn→Rn with matrix A and every Lebesgue measurable E, λn(T[E])=∣det⁡A∣λn(E); in particular λn(−E)=λn(E) (A linear map T of Rn sends Lebesgue measurable sets to Lebesgue measurable sets, with λn(T[E])=∣det⁡T∣ λn(E) when T is invertible and T[E] Lebesgue null when it is not).

[F3]

Under Countable Choice, σ is a finite Borel measure on Sn−1 (Polar coordinates decompose Lebesgue measure into r^{n-1} dr d sigma).

Proof

1.1F1F2algebra

Reflection invariance. For Borel E⊆Sn−1 the set −E is Borel and {rθ:θ∈−E, 0<r≤1}=−{rω:ω∈E, 0<r≤1}, so [F1] and [F2] give σ(−E)=nλn(−{rω:ω∈E, 0<r≤1})=nλn{rω:ω∈E, 0<r≤1}=σ(E).

1.2F3algebra

Vanishing of the moment. Each coordinate function gi(ω)=ωi is Borel and bounded by 1 on Sn−1, so ∫gi dσ is a finite real number by [F3]. Because σ is invariant under the bijection ω↦−ω, the substitution formula for a measure-preserving bijection — valid for indicators by definition, for simple functions by linearity, and for bounded Borel functions by the supremum definition of the integral — gives ∫Sn−1gi dσ=∫Sn−1gi(−ω) dσ(ω)=−∫Sn−1gi dσ; hence ∫gi dσ=0 for every i, that is ∫Sn−1ω dσ(ω)=0.

2.1algebra∎

For a∈Rn, linearity of the integral in the integrand gives ∫Sn−1a⋅ω dσ(ω)=a⋅∫Sn−1ω dσ(ω)=0.

LemmaStatement: Literature-sourcedProof: AI-adaptedOpen item page →

Differentiating an integral with moving endpoints

Statement

Let I⊆R be an open interval, let α,β∈C1(I) with α(t)<β(t) for t∈I, let J⊆R be an open interval containing the closure of the union of the intervals [α(t),β(t)] over t∈I, and let F:I×J→R be continuous with continuous partial derivative ∂tF. Then G(t):=∫α(t)β(t)F(t,y) dy is C1 on I and G′(t)=F(t,β(t))β′(t)−F(t,α(t))α′(t)+∫α(t)β(t)∂tF(t,y) dy.

Facts & Assumptions

Given: open intervals I,J, functions α,β∈C1(I) with α<β, and a continuous F:I×J→R whose partial derivative ∂tF exists and is continuous on I×J, with J containing the closure of the union of the intervals [α(t),β(t)].

[F1]

If I0⊆R is order-convex with at least two elements and f:I0→R is continuous, then for c0∈I0, F(x)=∫c0xf is a primitive of f and ∫abf=G(b)−G(a) for every primitive G of f and a<b in I0 (Every continuous function on an interval has a primitive; two primitives differ by a constant; and ∫abf=G(b)−G(a) for any primitive G).

[F2]

Let a<b, c<d and let g,h:[a,b]×[c,d]→R be continuous with x↦g(x,t) differentiable on (a,b) and derivative h(x,t) for every fixed t∈[c,d]. Then G(x)=∫cdg(x,t) dt is differentiable on [a,b] with G′(x)=∫cdh(x,t) dt (Leibniz's rule on a compact rectangle: an interior parameter derivative with a continuous extension may be passed through a Riemann integral).

[F3]

If f:U→V⊆Rn is totally differentiable at a and g:V→Rp is totally differentiable at f(a), then g∘f is totally differentiable at a with D(g∘f)(a)=Dg(f(a))∘Df(a) (The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)). The required total differentiability follows from continuous coordinate partial derivatives (If all partial derivatives exist on a neighbourhood and are continuous at a point, then the map is totally differentiable there with Jacobian derivative).

Proof

1.1F1given

Localisation. Fix t0∈I. Choose a compact interval K⊂I with t0 in its interior I0. The endpoint functions are bounded on K by A continuous real-valued function on a nonempty compact metric space is bounded and attains a greatest and a least value, so choose a<b in J with a<α(t)<β(t)<b for t∈K, and choose r∗∈J with r∗<a. Define Φ(u,t):=∫r∗uF(t,y) dy on (a,b)×I0. By [F1], G(t)=Φ(β(t),t)−Φ(α(t),t).

1.2F1F2algebra

The primitive is C1 on an open neighbourhood of the endpoint curves. By [F1], ∂uΦ(u,t)=F(t,u); by [F2] on compact rectangles inside I×J, ∂tΦ(u,t)=∫r∗u∂tF(t,y) dy. This last expression is jointly continuous: on a fixed compact rectangle, uniform continuity of ∂tF (Heine-Cantor: a continuous map from a compact metric space to any metric space is uniformly continuous) bounds the change in t by the interval length times a uniform error, and boundedness bounds the change in u by a constant times ∣u−u0∣. Thus both partial derivatives are continuous, and If all partial derivatives exist on a neighbourhood and are continuous at a point, then the map is totally differentiable there with Jacobian derivative makes Φ totally differentiable.

2.1F1F3F4step 1.2algebra∎

Apply [F3] to the curves t↦(β(t),t) and t↦(α(t),t) in the open domain of Φ, and subtract using [F4]. Evaluation of the primitive gives G′(t)=F(t,β(t))β′(t)−F(t,α(t))α′(t)+∫α(t)β(t)∂tF(t,y) dy. The same uniform estimate as in step 1.2 shows this derivative is continuous. Since t0 was arbitrary, the formula holds throughout I.

DefinitionDefinition: Literature-sourcedProof: Not applicablejudge pass (gpt-6.1-sol)Open item page →

Spherical means and the weighted ball integral of space-dependent data

Definition

Assume the Axiom of Countable Choice. Let n≥1, let f:Rn→R be continuous, let σ be the polar surface measure on Sn−1 (The polar surface set function on the unit sphere) and let ωn−1=σ(Sn−1)=nVn>0, where Vn=∣B1n∣ (Sphere and ball measures scale in Rn). The spherical mean of f is Mf(x,r):=1ωn−1∫Sn−1f(x+rω) dσ(ω)(x∈Rn, r>0), the average of f over the sphere of centre x and radius r with respect to the polar measure; this is the mean of Spherical averages and local ball means in Rn with u=f, restricted to continuous data. One sets Mf(x,0):=f(x); that value is a convention whose consistency as a limit is proved later on this page, not assumed here. The unnormalised sphere integral is Sf(x,r):=∫Sn−1f(x+rω) dσ(ω)=ωn−1Mf(x,r)(x∈Rn, r>0). For even n and r>0 put Wf(x,r):=1n!!Vn∫Br(x)f(y)r2−∣y−x∣2 dy, where n!!=n(n−2)⋯2 and Vn=∣B1n∣ is the volume of the unit ball (Sphere and ball measures scale in Rn); for n=2 this is the weighted disk integral appearing in the two-dimensional Poisson formula below. The integral defining Wf(x,r) is absolutely convergent and hence well defined: the weight y↦(r2−∣y−x∣2)−1/2 is integrable over Br(x) — by translation invariance (Lebesgue outer measure, Lebesgue measurability and Lebesgue measure are unchanged by translation) and the polar formula its integral is ωn−1∫0rsn−1(r2−s2)−1/2ds<∞ (Polar coordinates decompose Lebesgue measure into r^{n-1} dr d sigma) — while f is bounded on the closed ball Br(x)‾ because that ball is compact (For n≥1, every Euclidean closed ball and every Euclidean sphere of positive radius is compact) and a continuous function is bounded on a compact set (A continuous real-valued function on a nonempty compact metric space is bounded and attains a greatest and a least value); the product of an integrable function and a bounded function is Lebesgue integrable. The ball average used on this page is the normalised mean Ag(x,r) of The average of a locally integrable function over a Euclidean ball, and every sphere or ball integral below is read under the Axiom of Countable Choice of The Axiom of Countable Choice (ACω), which supplies the polar measure and its integrals.

LemmaStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

Factorisation of the one-dimensional wave operator

Statement

Let c>0 and let u be C2 on an open subset of R2 (Wave equation, Cauchy data and wave speed). Then ∂t2u−c2∂x2u=(∂t−c∂x)(∂t+c∂x)u=(∂t+c∂x)(∂t−c∂x)u. In the characteristic coordinates ξ=x−ct, η=x+ct one has ∂t2u−c2∂x2u=−4c2 ∂ξ∂ηu, so u solves the homogeneous one-dimensional wave equation exactly on the open set where uξη=0.

Facts & Assumptions

Given: a speed c>0 and a C2 function u on an open subset of R2, with coordinates (x,t) and characteristic coordinates ξ=x−ct, η=x+ct.

[F1]

If f is C2 on an open subset of Rm, then ∂i∂jf=∂j∂if for every pair of coordinate indices (Clairaut--Schwarz theorem for continuous second partial derivatives).

[F2]

If f:U→V⊆Rn is totally differentiable at a and g:V→Rp is totally differentiable at f(a), then g∘f is totally differentiable at a with D(g∘f)(a)=Dg(f(a))∘Df(a) (The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)). The required total differentiability follows from continuous coordinate partial derivatives (If all partial derivatives exist on a neighbourhood and are continuous at a point, then the map is totally differentiable there with Jacobian derivative).

Proof

1.1F1F3algebra

Expanding the two compositions and using [F1] for the mixed terms, (∂t−c∂x)(∂t+c∂x)u=∂t2u+c∂t∂xu−c∂x∂tu−c2∂x2u=∂t2u−c2∂x2u, and (∂t+c∂x)(∂t−c∂x)u=∂t2u−c∂t∂xu+c∂x∂tu−c2∂x2u=∂t2u−c2∂x2u; hence both factorisations equal the operator applied to u.

1.2F2F3algebra

Write U(ξ,η):=u(x,t) with ξ=x−ct, η=x+ct. By [F2] the chain rule for the substitution (x,t)↦(ξ,η)=(x−ct,x+ct) gives ∂xu=Uξ+Uη and ∂tu=−cUξ+cUη, with the right-hand sides evaluated at (x−ct,x+ct), hence ∂t−c∂x=−2c∂ξ and ∂t+c∂x=2c∂η as operators on U; composing, ∂t2u−c2∂x2u=(∂t−c∂x)(∂t+c∂x)u=(−2c∂ξ)(2c∂η)U=−4c2Uξη.

2.1algebra∎

Since c>0 the factor −4c2 is nonzero, so ∂t2u=c2∂x2u at a point if and only if Uξη=0 there; this proves the claimed equivalence and completes the factorisation identities.

LemmaStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

The radial recursion between dimensions n and n+2

Statement

Let c>0, n≥1 and let w be a C3 function of (r,t) on a domain with r>0. Write Lmw:=wtt−c2(wrr+m−1rwr) for the radial m-dimensional wave operator and Δmg:=grr+m−1rgr for the radial Laplacian, in the sign convention of the wave operator of Wave equation, Cauchy data and wave speed. Then Ln+2[1r∂rw]=1r∂r[Lnw](r>0). Consequently w↦r−1∂rw maps radial classical solutions of the m-dimensional homogeneous wave equation to radial classical solutions of the (m+2)-dimensional one; at m=1 this is the correspondence between one-dimensional waves and radial three-dimensional waves. It is the mechanism by which the kernels in dimension n+2 are radial derivatives of the kernels in dimension n.

Facts & Assumptions

Given: a speed c>0, an integer n≥1, and a C3 function w of (r,t) on a domain with r>0.

[F2]

If f is C2 on an open subset of Rm, then ∂i∂jf=∂j∂if for every pair of coordinate indices (Clairaut--Schwarz theorem for continuous second partial derivatives).

Proof

1.1F1algebra

Put h:=r−1wr. Differentiating the product and quotient, hr=r−1wrr−r−2wr and hrr=r−1wrrr−2r−2wrr+2r−3wr, so the radial (n+2)-dimensional Laplacian of h is Δn+2h=hrr+n+1rhr=r−1wrrr+n−1r2wrr−n−1r3wr.

1.2F1algebra

On the other hand Δnw=wrr+n−1rwr, whence ∂r(Δnw)=wrrr+n−1rwrr−n−1r2wr and r−1∂r(Δnw)=r−1wrrr+n−1r2wrr−n−1r3wr, the same expression as the radial Laplacian of the previous step.

2.1F2algebra∎

Since w is C3, [F2] gives wrtt=wttr, so time differentiation commutes with r−1∂r: htt=(r−1wr)tt=r−1∂r(wtt). Subtracting c2 times the identity Δn+2h=r−1∂r(Δnw) from this equality gives Ln+2h=htt−c2Δn+2h=r−1∂r(wtt)−c2r−1∂r(Δnw)=r−1∂r(wtt−c2Δnw)=r−1∂r(Lnw); if in particular Lnw=0 then Ln+2(r−1wr)=0, which is the stated mapping of radial solutions.

LemmaStatement: Literature-sourcedProof: AI-adaptedOpen item page →

Ball means and sphere means are related by a radial derivative

Statement

Assume the Axiom of Countable Choice, let n≥1, let g∈C0(Rn), and for x∈Rn, r>0 let Ag(x,r):=∣Br(x)∣−1∫Br(x)g be the ball average of The average of a locally integrable function over a Euclidean ball and Mg the spherical mean of Spherical means and the weighted ball integral of space-dependent data. Then for every r>0 Ag(x,r)=nrn∫0rsn−1Mg(x,s) ds,Mg(x,r)=Ag(x,r)+rn ∂rAg(x,r).

Facts & Assumptions

Given: Countable Choice, n≥1, g∈C0(Rn), and the means Ag and Mg of the statement.

[F1]

Under Countable Choice, ∫Rnf dλn=∫0∞∫Sn−1f(sθ)sn−1 dσ(θ) ds for every Borel measurable f:Rn→[0,∞] (Polar coordinates decompose Lebesgue measure into r^{n-1} dr d sigma).

[F2]

Mg(x,s)=ωn−1−1∫Sn−1g(x+sω) dσ(ω) and ωn−1=σ(Sn−1) for x∈Rn, s>0 (Spherical means and the weighted ball integral of space-dependent data).

[F3]

∣Br(x)∣=ωn−1rn/n for r>0 (Sphere and ball measures scale in Rn).

[F5]

If I⊆R is order-convex with at least two elements and f:I→R is continuous on I, then for r0∈I, r↦∫r0rf is a primitive of f on I (Every continuous function on an interval has a primitive; two primitives differ by a constant; and ∫abf=G(b)−G(a) for any primitive G).

Proof

1.1F1F2F3algebra

Polar coordinates applied to the positive and negative parts of the function u↦g(x+u)1Br(0)(u) give, using translation invariance (Lebesgue outer measure, Lebesgue measurability and Lebesgue measure are unchanged by translation) to set y=x+u, ∣Br(x)∣Ag(x,r)=∫Br(x)g=∫0r∫Sn−1g(x+sω)sn−1 dσ(ω) ds=ωn−1∫0rsn−1Mg(x,s) ds, where the inner integral is ωn−1Mg(x,s) by the definition of the spherical mean; dividing by ∣Br(x)∣=ωn−1rn/n gives the first identity.

1.2F4algebra

The function s↦Mg(x,s) is continuous on (0,∞): for s,s0>0 with ∣s−s0∣<1 one has ∣Mg(x,s)−Mg(x,s0)∣≤sup⁡ω∈Sn−1∣g(x+sω)−g(x+s0ω)∣, and the two points x+sω and x+s0ω lie in the compact ball Bs0+1(x)‾ at distance ∣s−s0∣; since g is uniformly continuous on that ball by [F4], the supremum tends to 0 as s→s0.

1.3F5algebra

Hence s↦sn−1Mg(x,s) is continuous on (0,∞), so locally at any r>0 split the integral at a fixed r0∈(0,r); the first part is constant and [F5] differentiates the second part. Thus r↦∫0rsn−1Mg(x,s) ds is a primitive of sn−1Mg(x,s), and the first identity gives ∂r(rnAg(x,r))=n rn−1Mg(x,r) for r>0.

2.1algebra∎

Differentiating the product gives ∂r(rnAg(x,r))=nrn−1Ag(x,r)+rn∂rAg(x,r), so equating with the previous display and dividing by nrn−1>0 yields Mg(x,r)=Ag(x,r)+rn∂rAg(x,r), the second identity.

LemmaStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

Sphere integrals of a cylindrical function project to weighted ball integrals

Statement

Assume the Axiom of Countable Choice. Let n≥1 and g∈C0(Rn), and let G(ξ,z):=g(ξ) be the extension of g to Rn+1 independent of the last coordinate. For r>0 and x∈Rn, with Srn(x)={(ξ,z)∈Rn+1:∣ξ−x∣2+z2=r2} the sphere of radius r in Rn+1 and ωn=∣Sn∣ its total polar measure (The polar surface set function on the unit sphere), ∫Srn(x)G dS=2r∫Brn(x)g(y)r2−∣y−x∣2 dy, and the spherical mean of G over Srn(x) equals 2ωnrn−1∫Brn(x)g(y)r2−∣y−x∣2dy. In particular, for even n=2k, every c>0 and t>0, with Wg the weighted ball integral of Spherical means and the weighted ball integral of space-dependent data, tn−1MG((x,0),ct)=(n−1)!!cn−1 Wg(x,ct), where MG is the (n+1)-dimensional spherical mean of G with centre (x,0). The connecting constant identity is 2 n!!Vn/ωn=(n−1)!!.

Facts & Assumptions

Given: Countable Choice, n≥1, g∈C0(Rn), the cylindrical extension G(ξ,z)=g(ξ), and r>0, x∈Rn.

[F1]

Surface measure is chart-independent; in graph coordinates (y,h(y)) its density is 1+∣Dh(y)∣2 (Surface integration on compact C1 hypersurfaces, Chart and partition independence of surface measure). On spheres it agrees with polar measure and scales by the appropriate radius power (Agreement with the existing polar sphere measure).

[F2]

Under Countable Choice, ∫Rnf dλn=∫0∞∫Sn−1f(sθ)sn−1 dσ(θ) ds for every Borel f:Rn→[0,∞] (Polar coordinates decompose Lebesgue measure into r^{n-1} dr d sigma).

[F3]

∣∂Brm∣=ωm−1rm−1 and ∣Brm∣=ωm−1rm/m, so in particular ∣Srn(x)∣=ωnrn for the sphere in Rn+1 and ωm=(m+1)Vm+1 for the unit sphere Sm (Sphere and ball measures scale in Rn).

[F4]

Vm=πm/2/Γ(m/2+1) for every m≥1 (The closed form for the volume of the unit n-ball).

[F5]

Γ(s+1)=sΓ(s) for s>0 with Γ(1)=1 (The real Gamma functional equation Γ(s+1)=sΓ(s)); Γ(1/2)=π (Γ(1/2)=π from the Gaussian integral); Γ(m+1)=m! for every integer m≥0 (Gamma at the positive integers).

Proof

1.1F1algebra

Graph sheets. The upper and lower open hemispheres of Srn(x) are the graphs y↦(y,±r2−∣y−x∣2) over Brn(x). Their graph density is 1+∣y−x∣2/(r2−∣y−x∣2)=r/r2−∣y−x∣2 by [F1]. The equator has zero surface measure: near each of its points choose a sphere graph omitting a nonzero one of the first n coordinates. Its parameter set for the equator lies in the coordinate hyperplane z=0, which has Lebesgue measure zero (for n=1, it is a singleton, null because it lies in intervals of arbitrarily small length); Fubini's theorem for L^1 functions on a sigma-finite product applied to its indicator proves nullity, and the continuous graph density preserves it. A finite chart cover suffices by compactness of the sphere (For n≥1, every Euclidean closed ball and every Euclidean sphere of positive radius is compact).

1.2F1F2algebra

Integrating the two sheets gives ∫Srn(x)G dS=2r∫Brn(x)g(y)(r2−∣y−x∣2)−1/2 dy. These integrals are absolutely convergent: g is bounded on the closed ball by A continuous real-valued function on a nonempty compact metric space is bounded and attains a greatest and a least value, and [F2] reduces the weight integral to ωn−1∫0rsn−1(r2−s2)−1/2ds, bounded by ωn−1rn−1∫0r(r2−s2)−1/2ds=ωn−1rn−1π/2. Thus the graph computation applies separately to positive and negative parts.

1.3F3algebra

The mean. By [F3], ∣Srn(x)∣=ωnrn, so the spherical mean of G over Srn(x) is ωn−1r−n∫Srn(x)G dS=2ωnrn−1∫Br(x)g(y)r2−∣y−x∣2dy.

1.4F3F4F5algebra

The even-dimensional form. Let n=2k be even, c>0, t>0 and r=ct. Substituting r=ct in the mean identity, tn−1MG((x,0),ct)=2ωncn−1∫Bct(x)g(y)c2t2−∣y−x∣2dy=2n!!Vnωncn−1Wg(x,ct) by the definition of Wg. For the constant: Vn=πn/2/Γ(k+1)=πk/k! by [F4]; the functional equation and Γ(1/2)=π give by induction Γ(k+1/2)=(2k−1)(2k−3)⋯12kπ=(2k)!π4kk!, so with ωn=∣Sn∣=(n+1)Vn+1=2πk+1/2Γ(k+1/2)=22k+1πkk!(2k)! one gets 2n!!Vn=2⋅2kk!⋅πk/k!=2k+1πk and 2n!!Vnωn=2k+1πk(2k)!22k+1πkk!=(2k)!2kk!=(2k−1)!!=(n−1)!!.

2.1algebra∎

Substituting the constant gives tn−1MG((x,0),ct)=(n−1)!!cn−1Wg(x,ct), which together with the two integral identities proves all assertions.

LemmaStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

General solution of the one-dimensional wave equation

Statement

Let c>0 and let R⊆R2 be a nonempty open rectangle. If u∈C2(R) satisfies utt=c2uxx on R, then there are intervals I,J⊆R and functions F∈C2(I), G∈C2(J) with u(x,t)=F(x−ct)+G(x+ct)((x,t)∈R), where I and J are the projections of R onto the ξ- and the η-axis under ξ=x−ct, η=x+ct. Conversely, every such sum is a C2 solution of utt=c2uxx on R. The pair (F,G) is unique up to the replacement F↦F+k, G↦G−k with k∈R, and no other freedom remains.

Facts & Assumptions

Given: a speed c>0, a nonempty open rectangle R⊆R2, and a C2 function u on R.

[F1]

On the domain of u, ∂t2u−c2∂x2u=−4c2∂ξ∂ηU where U(ξ,η):=u(x,t) and ξ=x−ct, η=x+ct (Factorisation of the one-dimensional wave operator).

[F2]

If f:U→V⊆Rn is totally differentiable at a and g:V→Rp at f(a), then g∘f is totally differentiable at a with D(g∘f)(a)=Dg(f(a))∘Df(a) (The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)). The required total differentiability follows from continuous coordinate partial derivatives (If all partial derivatives exist on a neighbourhood and are continuous at a point, then the map is totally differentiable there with Jacobian derivative).

[F3]

Let I⊆R be order-convex and f:I→R continuous and differentiable at every interior point with f′=0 there. Then f is constant; if moreover f,g are continuous with f′=g′ at every interior point then f−g is constant (A function continuous on an interval I whose derivative vanishes at every interior point of I is constant on I; consequently two such functions with the same derivative differ by a constant).

[F4]

Let I be order-convex with at least two elements and f:I→R continuous. Fix c0∈I; then F(x)=∫c0xf is a primitive of f on I, and primitives differ by constants (Every continuous function on an interval has a primitive; two primitives differ by a constant; and ∫abf=G(b)−G(a) for any primitive G).

Proof

1.1F1F2algebra

The affine change of variables (x,t)↦(ξ,η)=(x−ct,x+ct) is a bijection of R2 with inverse x=(ξ+η)/2, t=(η−ξ)/(2c), and the image R′ of the rectangle R is a nonempty open convex affine image of R (a parallelogram when R is bounded); the projections I of R′ onto the ξ-axis and J onto the η-axis are nonempty open intervals, and every section {η:(ξ,η)∈R′} and {ξ:(ξ,η)∈R′} is a nonempty open interval. By [F2] the function U(ξ,η):=u((ξ+η)/2,(η−ξ)/(2c)) is C2 on R′, and [F1] gives utt−c2uxx=−4c2Uξη on R′.

1.2F3F4algebra

Suppose u solves utt=c2uxx on R; since c>0, Uξη=0 on R′. For fixed ξ∈I the section Jξ:={η:(ξ,η)∈R′} is a nonempty open interval and ∂η(Uξ(ξ,⋅))=Uξη(ξ,⋅)=0 there, so by [F3] the value Uξ(ξ,η) is independent of η in that section; call it p(ξ). To see p is C1, fix ξ0∈I and choose η0 with (ξ0,η0)∈R′; openness gives an interval I0 about ξ0 on which (ξ,η0)∈R′, so p(ξ)=Uξ(ξ,η0) on I0. Since U is C2, this local representative is C1, and hence p∈C1(I). Choose a primitive P of p on I by [F4]; since p∈C1, P∈C2(I). Then ∂ξ(U−P)=0 on R′, and the same section argument in the ξ-direction gives that U−P is independent of ξ: there is G:J→R with U(ξ,η)−P(ξ)=G(η) for all (ξ,η)∈R′. To see G is C2, fix η0∈J and choose ξ0 with (ξ0,η0)∈R′; openness gives a neighbourhood J0 on which G(η)=U(ξ0,η)−P(ξ0), a C2 function. Thus G∈C2(J) and u(x,t)=P(x−ct)+G(x+ct) on R.

1.3F2F5algebra

Conversely, if F∈C2(I) and G∈C2(J), then (x,t)↦F(x−ct)+G(x+ct) is C2 on R by [F2] and [F5], and two applications of the chain rule give ∂t2[F(x−ct)+G(x+ct)]=c2F′′(x−ct)+c2G′′(x+ct)=c2∂x2[F(x−ct)+G(x+ct)]; hence every such sum solves the wave equation.

1.4F3F5algebra

Uniqueness of the pair. If F1(ξ)+G1(η)=F2(ξ)+G2(η) for all (ξ,η)∈R′, then Φ(ξ,η):=F1(ξ)−F2(ξ)−G2(η)+G1(η) vanishes on R′. For fixed ξ∈I the section in η is a nonempty interval, so ∂ηΦ=G1′−G2′=0 on J, and for fixed η∈J similarly ∂ξΦ=F1′−F2′=0 on I; hence F1−F2 and G2−G1 are the same constant k by [F3], that is F1=F2+k and G1=G2−k.

2.1given∎

Therefore every C2 solution of the one-dimensional homogeneous wave equation on a nonempty open rectangle has the form F(x−ct)+G(x+ct) with F,G of class C2 on the projections, every such sum is a solution, and the decomposition is unique up to the additive shift F↦F+k, G↦G−k.

LemmaStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

Smoothness, parity and zero-radius limits of spherical means

Statement

Assume the Axiom of Countable Choice of Spherical means and the weighted ball integral of space-dependent data. Let n≥1, k≥1 and f∈Ck(Rn), and let Mf be the spherical mean with the convention Mf(x,0):=f(x). Then: (i) (x,r)↦Mf(x,r) is Ck on Rn×(0,∞), and every derivative is obtained by differentiating f under the sphere integral: for every multi-index α and every m≥0 with ∣α∣+m≤k, Dxα∂rmMf(x,r)=1ωn−1∫Sn−1Dxα∂rm[f(x+rω)] dσ(ω)(x∈Rn, r>0), where ∂rm[f(x+rω)] is the m-th r-derivative of the composed function r↦f(x+rω). (ii) Mf extends to a continuous function on Rn×[0,∞) with Mf(x,0)=f(x), and ∂rMf(x,r)→0 as r↓0, uniformly for x in compact subsets of Rn. (iii) The signed-radius integral ωn−1−1∫Sn−1f(x+rω) dσ(ω) for r∈R is an even Ck extension of Mf to Rn×R. Its differentiated-integral formula holds also at r=0; in particular every available odd-order radial derivative vanishes there. (iv) ∣∂rMf(x,r)∣≤sup⁡Br(x)∣Df∣ for every x∈Rn and r>0.

Facts & Assumptions

Given: Countable Choice, n≥1, k≥1, f∈Ck(Rn), and the spherical mean Mf with Mf(x,0):=f(x).

[F1]

The spherical mean is integration against the finite measure σ/ωn−1 of total mass one (Spherical means and the weighted ball integral of space-dependent data).

[F4]

Under Countable Choice reflection ω↦−ω preserves the polar measure, and its first moment vanishes: ∫Sn−1ω dσ(ω)=0 (Reflection invariance and vanishing first moment of the sphere measure).

Proof

1.1F1F2F3algebra

Differentiation, including signed radii. Put H(x,r):=ωn−1−1∫Sn−1f(x+rω) dσ(ω) for all real r. On a bounded parameter neighbourhood, all points x+rω lie in a fixed compact ball. For a coordinate parameter p and a continuous derivative integrand q(p,ω) with continuous ∂pq, [F2] gives ∣(q(p+h,ω)−q(p,ω))/h−∂pq(p,ω)∣≤sup⁡∣s−p∣≤∣h∣,ω∣∂pq(s,ω)−∂pq(p,ω)∣. Uniform continuity on the compact ball makes this bound tend to zero uniformly in ω; integration against the probability measure [F1] preserves the bound. Iterating through total order k therefore gives every stated derivative under the integral, with continuity again following from the same uniform estimate. This works for n=1 as well, since it requires only a finite measure, not positive-dimensional surface charts.

1.2F3algebra

Part (ii), continuity at r=0. Let K⊆Rn be compact and choose R>0 with K⊆BR(0); the ball BR+1(0)‾ is compact, so by [F3] f is bounded there and uniformly continuous on it. Given ε>0 choose δ∈(0,1) such that ∣f(u)−f(v)∣<ε whenever u,v∈BR+1(0)‾ and ∣u−v∣<δ. For x∈K, 0≤r<δ and ω∈Sn−1 one has x+rω∈BR+1(0)‾ and ∣x+rω−x∣=r<δ, so ∣Mf(x,r)−f(x)∣≤sup⁡ω∣f(x+rω)−f(x)∣<ε. Hence Mf(x,r)→f(x) as r↓0 uniformly on compact subsets, and the extension by Mf(⋅,0)=f is continuous because f is.

1.3F3F4algebra

Part (ii), the derivative limit. By part (i) with ∣α∣=0, m=1, ∂rMf(x,r)=ωn−1−1∫Sn−1Df(x+rω)⋅ω dσ(ω) for r>0, so adding and subtracting Df(x) inside the integral gives ∂rMf(x,r)=ωn−1−1∫Sn−1[Df(x+rω)−Df(x)]⋅ω dσ(ω)+ωn−1−1Df(x)⋅∫Sn−1ω dσ(ω), where the second term vanishes by [F4]; the first is bounded by sup⁡ω∈Sn−1∣Df(x+rω)−Df(x)∣, which tends to 0 as r↓0 uniformly for x in a fixed compact set by uniform continuity of the continuous function Df on a large compact ball, again by [F3]. Hence ∂rMf(x,r)→0 uniformly on compact subsets of Rn.

2.1F1F3F4step 1.1algebra

Parity and the bound. Reflection preserves σ by [F4], hence the substitution ω↦−ω gives H(x,−r)=H(x,r). Since H is Ck by step 1.1, its odd-order radial derivatives at zero vanish whenever their orders are at most k. For (iv), ∣∂rMf(x,r)∣≤ωn−1−1∫∣Df(x+rω)∣ dσ(ω). By continuity, each boundary value of ∣Df∣ is at most sup⁡Br(x)∣Df∣, and integration gives the stated bound.

3.1step 1.1step 1.2step 1.3step 2.1∎

Thus Mf has the differentiated-integral formula, the stated uniform zero-radius limits and gradient bound, and an even Ck signed-radius extension.

LemmaStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

The Euler–Poisson–Darboux equation for spherical means

Statement

Assume the Axiom of Countable Choice, let n≥1, let f∈C2(Rn) and let M=Mf be the spherical mean of Spherical means and the weighted ball integral of space-dependent data. Then for every x∈Rn and r>0 ΔxM(x,r)=Mrr(x,r)+n−1rMr(x,r). The right-hand side extends continuously to r=0 with value Δf(x): with M(x,0):=f(x) one has Mr(x,r)→0 and (Mrr+n−1rMr)(x,r)→Δf(x) as r↓0. If u∈C2(Rn×I) is a classical solution of □cu=0 on an interval I, then its space-time mean (x,r,t)↦Mu(x,r,t):=ωn−1−1∫Sn−1u(x+rω,t) dσ(ω) satisfies ∂t2Mu=c2(Mrr+n−1rMr). No equation for f is needed for the identity itself.

Facts & Assumptions

Given: Countable Choice, n≥1, f∈C2(Rn), the spherical mean M=Mf with M(x,0)=f(x), and, when stated, a C2 solution u of □cu=0 on Rn×I.

[F1]

For k≥1 and h∈Ck(Rn), (x,r)↦Mh(x,r) is Ck on Rn×(0,∞) with all derivatives obtained by differentiating h under the sphere integral; it extends continuously to r=0 with Mh(x,0)=h(x), and its radial derivative tends to 0 uniformly on compacta (Smoothness, parity and zero-radius limits of spherical means).

[F2]

If g∈C2(Rn) and r>0, then ∂rMg(x,r)=rnAΔg(x,r), where AΔg(x,r)=∣Br(x)∣−1∫Br(x)Δg is the ball average (Radial derivative of a spherical average, The average of a locally integrable function over a Euclidean ball).

[F3]

For h∈C0(Rn) one has Mh(x,r)=Ah(x,r)+rn∂rAh(x,r) for all r>0 (Ball means and sphere means are related by a radial derivative).

Proof

1.1F1algebra

The spatial Laplacian passes under the sphere integral: by [F1] the map (x,r)↦M(x,r) is C2, and differentiating the defining integral twice in x, ΔxM(x,r)=ωn−1−1∫Sn−1Δf(x+rω) dσ(ω)=MΔf(x,r) for every r>0.

1.2F2F3algebra

Radial identities. Put A(x,r):=AΔf(x,r). Applying [F2] only to f gives Mr(x,r)=rnA(x,r). Since Δf is continuous, [F3] applied to Δf gives MΔf(x,r)=A(x,r)+rn∂rA(x,r). Differentiating the first identity in r gives Mrr(x,r)=1nA(x,r)+rn∂rA(x,r); therefore Mrr+n−1rMr=1nA+rn∂rA+n−1nA=A+rn∂rA=MΔf(x,r)=ΔxM(x,r). The radial-derivative formula is used only with the C2 function f; no differentiability of Δf is assumed.

1.3F4algebra

The limits at r↓0. For 0<r≤1 the bound ∣AΔf(x,r)∣≤sup⁡B1(x)∣Δf∣<∞ holds, the supremum being finite by [F4]; hence Mr(x,r)=rnAΔf(x,r)→0. Moreover the identity just proved gives Mrr(x,r)+n−1rMr(x,r)=MΔf(x,r)=ΔxM(x,r). Since Δf is continuous at x, ∣MΔf(x,r)−Δf(x)∣≤sup⁡ω∈Sn−1∣Δf(x+rω)−Δf(x)∣→0 as r↓0.

1.4F1algebra

The space-time version. If u∈C2 solves □cu=0, applying the differentiation-under-the-sphere-integral computation of [F1] to the C2 function (x,t)↦u(x,t) yields ΔxMu(x,r,t)=MΔxu(x,r,t) and ∂t2Mu(x,r,t)=Mutt(x,r,t) for all r>0 and t∈I. Since utt=c2Δu, this gives ∂t2Mu=c2MΔu=c2ΔxMu; and the identity of the first two steps applied to the C2 spatial function u(⋅,t) gives ΔxMu(x,r,t)=Mrr(x,r,t)+n−1rMr(x,r,t). Substituting, ∂t2Mu=c2(Mrr+n−1rMr).

2.1given∎

Collecting: the Euler–Poisson–Darboux identity holds for every C2 datum, its right-hand side has the stated continuous extension at r=0 with value Δf(x), and the space-time means of solutions satisfy the same radial equation with two time derivatives on the left.

TheoremStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

d'Alembert's formula and uniqueness in one dimension

Statement

Let c>0, u0∈C2(R) and u1∈C1(R). Then u(x,t):=12(u0(x−ct)+u0(x+ct))+12c∫x−ctx+ctu1(y) dy is a C2 function on R×[0,∞) and is the unique classical solution of the homogeneous Cauchy problem utt=c2uxx,u(⋅,0)=u0,ut(⋅,0)=u1. Its value depends on u0 only through the two endpoints x∓ct and on u1 only through its integral over [x−ct,x+ct].

Facts & Assumptions

Given: a speed c>0, data u0∈C2(R), u1∈C1(R), and the displayed function u.

[F1]

A continuous function on an interval has the primitive P(z)=∫0zu1(y) dy when the interval is R; since u1∈C1, this primitive is C2, and ∫abu1=P(b)−P(a) (Every continuous function on an interval has a primitive; two primitives differ by a constant; and ∫abf=G(b)−G(a) for any primitive G).

[F2]

If f is totally differentiable at a and g at f(a), then g∘f is totally differentiable at a with D(g∘f)(a)=Dg(f(a))∘Df(a) (The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)). The required total differentiability follows from continuous coordinate partial derivatives (If all partial derivatives exist on a neighbourhood and are continuous at a point, then the map is totally differentiable there with Jacobian derivative).

[F4]

Every C2 solution of utt=c2uxx on a nonempty open rectangle has the form F(x−ct)+G(x+ct) with F,G of class C2 on the projections, and conversely; the pair is unique up to F↦F+k, G↦G−k (General solution of the one-dimensional wave equation).

Proof

1.1F1F2F3algebra

Regularity and the equation. The first summand is C2 because u0∈C2, and the integral term is (P(x+ct)−P(x−ct))/(2c) with P∈C2 by [F1], so it is C2 even across t=0 by [F2]: ∂tu(x,t)=12(−cu0′(x−ct)+cu0′(x+ct))+12(u1(x+ct)+u1(x−ct)) and ∂xu(x,t)=12(u0′(x−ct)+u0′(x+ct))+12c(u1(x+ct)−u1(x−ct)). Differentiating once more by [F2] and [F3], ∂t2u=c22(u0′′(x−ct)+u0′′(x+ct))+c2(u1′(x+ct)−u1′(x−ct)) and c2∂x2u=c22(u0′′(x−ct)+u0′′(x+ct))+c2(u1′(x+ct)−u1′(x−ct)), so utt=c2uxx on R×[0,∞).

1.2F1algebra

Both data are attained. At t=0 the displacement terms give u(x,0)=12(u0(x)+u0(x))+0=u0(x), and the velocity formula of the previous step gives ∂tu(x,0)=12(−cu0′(x)+cu0′(x))+12(u1(x)+u1(x))=u1(x).

1.3F2F3F4algebra

Uniqueness. Let v be any classical solution with the same pointwise displacement and velocity limits at zero and put w=u−v. Apply [F4] to the open rectangle R×(0,T), where T>0 is arbitrary. Both characteristic projections are all of R, so w(x,t)=F(x−ct)+G(x+ct) with F,G∈C2(R). At each fixed x, letting t↓0 gives F(x)+G(x)=0 and −cF′(x)+cG′(x)=0. Differentiating the first identity and combining with the second yields F′=G′=0 everywhere; their sum is zero, so w=0 throughout this rectangle. As T is arbitrary, uniqueness holds on the whole time slab.

2.1given∎

The display of u therefore defines the unique classical solution, and its value at (x,t) involves u0 only through the endpoint values u0(x∓ct) and u1 only through the integral over [x−ct,x+ct].

LemmaStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

The d'Alembert expression attains both initial data

Statement

Let c>0, u0∈C2(R), u1∈C1(R) and let A be the d'Alembert expression of d'Alembert's formula and uniqueness in one dimension, A(x,t):=12(u0(x−ct)+u0(x+ct))+12c∫x−ctx+ctu1(y) dy. Then A∈C2(R×[0,∞)), A(x,0)=u0(x) for every x, and ∂tA(x,t)=12(−cu0′(x−ct)+cu0′(x+ct))+12(u1(x+ct)+u1(x−ct)), so ∂tA(x,0)=u1(x). The orientation of the velocity integral is the + sign in the second bracket: both ends of the characteristic base are traversed with speed c, and the two endpoint contributions add.

Facts & Assumptions

Given: a speed c>0, data u0∈C2(R), u1∈C1(R), and the expression A of the statement.

[F1]

Let α,β∈C1(I) with α<β and let F be continuous on I×J with continuous ∂tF, where J contains the closure of the union of the intervals [α(t),β(t)]. Then G(t)=∫α(t)β(t)F(t,y) dy is C1 with G′(t)=F(t,β(t))β′(t)−F(t,α(t))α′(t)+∫α(t)β(t)∂tF(t,y) dy (Differentiating an integral with moving endpoints).

[F3]

If f is totally differentiable at a and g at f(a), then g∘f is totally differentiable at a with D(g∘f)(a)=Dg(f(a))∘Df(a) (The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)). The required total differentiability follows from continuous coordinate partial derivatives (If all partial derivatives exist on a neighbourhood and are continuous at a point, then the map is totally differentiable there with Jacobian derivative).

Proof

1.1algebra

Displacement at t=0. At t=0 the two displacement terms are both u0(x) and the integral has equal endpoints, so A(x,0)=12(u0(x)+u0(x))+0=u0(x) for every x∈R.

1.2F1F2F3algebra

Velocity. For t>0, by [F1] applied to the integral term with α(t)=x−ct, β(t)=x+ct, α′(t)=−c, β′(t)=c and inner integrand u1, ∂tA(x,t)=12(−cu0′(x−ct)+cu0′(x+ct))+12c(c u1(x+ct)+c u1(x−ct))=12(−cu0′(x−ct)+cu0′(x+ct))+12(u1(x+ct)+u1(x−ct)); the displacement term is differentiated by [F2] and [F3]. Continuity of the data and the C2 regularity supplied by d'Alembert's formula and uniqueness in one dimension extend this derivative formula to t=0.

2.1F1algebra∎

Setting t=0 in the velocity formula gives ∂tA(x,0)=12(−cu0′(x)+cu0′(x))+12(u1(x)+u1(x))=u1(x); the regularity A∈C2(R×[0,∞)) is established in d'Alembert's formula and uniqueness in one dimension. Hence A attains both initial data, and the two endpoint contributions of the velocity integral add with a + sign as displayed.

CorollaryStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

The one-dimensional value depends on the characteristic interval

Statement

Let c>0 and let u be the solution of d'Alembert's formula and uniqueness in one dimension for data u0∈C2(R), u1∈C1(R). For every (x,t)∈R×[0,∞), the value u(x,t) is determined by the restrictions of u0 and u1 to the closed interval [x−ct,x+ct]: if (u~0,u~1) are admissible data agreeing with (u0,u1) there, then the corresponding solution satisfies u~(x,t)=u(x,t). In particular, changing the data outside [x−ct,x+ct] does not change the value at (x,t).

Facts & Assumptions

Given: a speed c>0, data u0∈C2(R), u1∈C1(R), a point (x,t)∈R×[0,∞), and admissible data (u~0,u~1) with u~0=u0 and u~1=u1 on [x−ct,x+ct].

[F1]

The d'Alembert formula of d'Alembert's formula and uniqueness in one dimension reads u(x,t)=12(u0(x−ct)+u0(x+ct))+12c∫x−ctx+ctu1(y) dy.

[F2]

The d'Alembert expression attains both initial data and defines the unique classical solution of the corresponding Cauchy problem (The d'Alembert expression attains both initial data, d'Alembert's formula and uniqueness in one dimension).

Proof

1.1F1algebra

The formula of [F1] evaluates u0 only at the two endpoints x−ct and x+ct of the interval and integrates u1 only over that interval; hence replacing (u0,u1) by any admissible pair with the same restrictions to [x−ct,x+ct] leaves the right-hand side unchanged, so the d'Alembert expression of the new data equals u(x,t) at the given point.

2.1F2algebra∎

By [F2] both expressions are the solutions of their respective Cauchy problems, so the solution u~ of the data (u~0,u~1) satisfies u~(x,t)=u(x,t); in particular the value at (x,t) is unchanged by altering the data off [x−ct,x+ct].

TheoremStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

The forced one-dimensional wave formula over the characteristic triangle

Statement

Let c>0, u0∈C2(R), u1∈C1(R) and let f be of class C1 on R×[0,∞), so that f, ∂xf and ∂tf are continuous. Then u(x,t)=12(u0(x−ct)+u0(x+ct))+12c∫x−ctx+ctu1(y) dy+12c∫0t∫x−c(t−s)x+c(t−s)f(y,s) dy ds is the unique classical solution of utt=c2uxx+f with u(⋅,0)=u0, ut(⋅,0)=u1. The source integral is over the backward characteristic triangle with vertex (x,t): 0≤s≤t, ∣y−x∣≤c(t−s), and its coefficient is 1/(2c).

Facts & Assumptions

Given: a speed c>0, data u0∈C2(R), u1∈C1(R), a source f∈C1(R×[0,∞)), and the displayed function u.

[F1]

Let α,β∈C1(I) with α<β and let F be continuous on I×J with continuous ∂tF, where J contains the closure of the union of the intervals [α(t),β(t)]. Then G(t)=∫α(t)β(t)F(t,y) dy is C1 with G′(t)=F(t,β(t))β′(t)−F(t,α(t))α′(t)+∫α(t)β(t)∂tF(t,y) dy (Differentiating an integral with moving endpoints).

[F2]

For continuous φ, the function x↦∫a(x)b(x)φ has derivative φ(b(x))b′(x)−φ(a(x))a′(x); more generally the primitive of a continuous function is recovered by evaluation at the endpoints (Every continuous function on an interval has a primitive; two primitives differ by a constant; and ∫abf=G(b)−G(a) for any primitive G).

[F3]

With data u0,u1 the homogeneous d'Alembert expression is the unique C2 solution of utt=c2uxx with those data (d'Alembert's formula and uniqueness in one dimension).

Proof

1.1F1F2F4algebra

The source term. Extend f(y,s) to s<0 by f(y,0); its value and first spatial derivative remain continuous. Define Φ(x,t,s):=∫x−c(t−s)x+c(t−s)f(y,s) dy as an oriented integral, also when t<s. Primitives [F2] give Φt=c[f(x+c(t−s),s)+f(x−c(t−s),s)] and Φx=f(x+c(t−s),s)−f(x−c(t−s),s) on an open rectangle in (t,s), including s=t. Since Φ(x,t,t)=0, [F1] gives Dt=12∫0t[f(x+c(t−s),s)+f(x−c(t−s),s)]ds for D=(2c)−1∫0tΦ ds. The compact-rectangle theorem Leibniz's rule on a compact rectangle: an interior parameter derivative with a continuous extension may be passed through a Riemann integral permits spatial differentiation under this fixed s-integral. Thus Dx=(2c)−1∫0t[f(x+c(t−s),s)−f(x−c(t−s),s)]ds and Dxx=(2c)−1∫0t[fx(x+c(t−s),s)−fx(x−c(t−s),s)]ds.

2.1F1F4step 1.1algebra

Equation and regularity. Differentiating Dt by [F1] gives Dtt=f(x,t)+c2∫0t[fx(x+c(t−s),s)−fx(x−c(t−s),s)]ds=f+c2Dxx. The chain rule The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a) gives the displayed integrand derivatives; differentiation of Dt in x gives Dtx=12∫0t[fx(x+c(t−s),s)+fx(x−c(t−s),s)]ds, also equal to Dxt by differentiating Dx. All these derivatives are continuous because their integrands are continuous on local compact rectangles. At zero, D,Dx,Dt,Dxx,Dxt tend to zero and Dtt→f(x,0) locally uniformly, proving the asserted C2 regularity up to the initial time.

3.1F3step 1.1step 2.1algebra

Data and uniqueness. At t=0 the source integral vanishes, so u(⋅,0)=u0 and ut(⋅,0)=u1 are exactly the statements of [F3] for the homogeneous part. If v is any classical solution of the forced problem with the same data, then w:=u−v is C2 with wtt=c2wxx and zero data, so w=0 by the uniqueness clause of [F3]; hence u is the unique classical solution.

4.1given∎

The double integral runs over 0≤s≤t and ∣y−x∣≤c(t−s), the backward characteristic triangle with vertex (x,t), and its coefficient is 1/(2c); this completes the identification of the displayed solution.

TheoremStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

Kirchhoff's formula in three dimensions

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let c>0, u0∈C3(R3), u1∈C2(R3) and let M be the spherical mean of Spherical means and the weighted ball integral of space-dependent data. Then u(x,t):=∂∂t[t Mu0(x,ct)]+t Mu1(x,ct) defines a C2 function on R3×(0,∞) solving utt=c2Δu. Equivalently, for t>0, u(x,t)=14πc2t2∫∂Bct(x)(u0(y)+∇u0(y)⋅(y−x)+t u1(y)) dS(y), the sphere integral being the unnormalised form of the mean because ω2=4π (Sphere and ball measures scale in Rn). The formula uses the values of u0,u1 and the normal derivative of u0 on ∂Bct(x); pointwise agreement of the two data only on that sphere need not give the same solution value.

Facts & Assumptions

Given: Countable Choice, c>0, u0∈C3(R3), u1∈C2(R3), and the means A(x,r):=Mu0(x,r), B(x,r):=Mu1(x,r).

[F1]

For k≥1 and h∈Ck(R3), the spherical mean Mh is Ck on R3×(0,∞), every derivative being obtained by differentiating h under the sphere integral (Smoothness, parity and zero-radius limits of spherical means).

[F2]

For h∈C2(R3) one has ΔxMh(x,r)=Mrr(x,r)+2rMr(x,r) for r>0 (The Euler–Poisson–Darboux equation for spherical means).

[F3]

If f is totally differentiable at a and g at f(a), then g∘f is totally differentiable at a with D(g∘f)(a)=Dg(f(a))∘Df(a) (The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)). The required total differentiability follows from continuous coordinate partial derivatives (If all partial derivatives exist on a neighbourhood and are continuous at a point, then the map is totally differentiable there with Jacobian derivative).

[F5]

If h is C2 on an open set of Rm, then ∂i∂jh=∂j∂ih (Clairaut--Schwarz theorem for continuous second partial derivatives).

[F6]

∣∂Br3∣=ω2r2 with ω2=3V3, and the map ω↦a+Rω multiplies surface measure by R2: ∫S2h(ω) dσ(ω)=R−2∫∂BR(a)h((y−a)/R) dS(y) (Sphere and ball measures scale in Rn, Agreement with the existing polar sphere measure).

Proof

1.1F1F3F4algebra

Time derivatives of the candidate. Write r=ct. By [F1] the means A,B are C3 respectively C2 on R3×(0,∞), so [F3] and [F4] give u=A+ctAr+tB, ∂tu=2cAr+c2tArr+B+ctBr and ∂t2u=3c2Arr+c3tArrr+2cBr+c2tBrr on R3×(0,∞), all derivatives being evaluated at (x,ct).

1.2F2F5F4algebra

Spatial derivatives. By [F2] applied to u0 and to u1, ΔxA=Arr+2rAr and ΔxB=Brr+2rBr; since A is C3 and ∂r commutes with the x-derivatives by [F5], ΔxAr=∂r(ΔxA)=Arrr+2rArr−2r2Ar. Hence Δu=ΔxA+ctΔxAr+tΔxB=(Arr+2rAr)+ct(Arrr+2rArr−2r2Ar)+t(Brr+2rBr) at (x,ct).

1.3F1F3F6F7algebra

The unnormalised form. Since ∣S2∣=ω2=4π by [F6] and [F7], Mu0(x,r)=14π∫S2u0(x+rω) dσ(ω)=14πr2∫∂Br(x)u0 dS, and likewise for u1; also ∂rMu0(x,r)=14π∫S2∇u0(x+rω)⋅ω dσ(ω)=14πr2∫∂Br(x)∇u0(y)⋅y−xr dS(y) by [F1], [F3] and the scaling in [F6]. Substituting r=ct into u=A+ctAr+tB and collecting the common factor 14πc2t2 gives u(x,t)=14πc2t2∫∂Bct(x)(u0(y)+∇u0(y)⋅(y−x)+tu1(y))dS(y).

2.1F4step 1.1step 1.2algebra

Comparison. Multiplying step 1.2 by c2 and substituting r=ct gives c2Δu=c2Arr+(2c/t)Ar+c3tArrr+2c2Arr−(2c/t)Ar+c2tBrr+2cBr=3c2Arr+c3tArrr+c2tBrr+2cBr. This equals utt from step 1.1, proving the equation. The C3 and C2 regularity of A and B gives u∈C2.

3.1step 2.1step 1.3∎

Both displays define the same C2 solution. The unnormalised display uses the data values and the normal derivative of u0 on the sphere, or equivalently the data on an open neighbourhood of that sphere.

TheoremStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

Poisson's formula in two dimensions by descent

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let c>0, u0∈C3(R2), u1∈C2(R2) and let W be the weighted ball integral of Spherical means and the weighted ball integral of space-dependent data. Then u(x,t):=12πc∂∂t∫Bct(x)u0(y)c2t2−∣y−x∣2 dy+12πc∫Bct(x)u1(y)c2t2−∣y−x∣2 dy(t>0) defines a C2 function on R2×(0,∞) solving utt=c2Δu; here the first ∂t differentiates the C1 function t↦∫Bct(x)u0(y)(c2t2−∣y−x∣2)−1/2dy, which is legitimate by the projection identity, not by termwise differentiation of a singular integrand. Equivalently, for t>0, u(x,t)=12πct∫Bct(x)u0(y)+∇u0(y)⋅(y−x)+t u1(y)c2t2−∣y−x∣2 dy. The extension to the initial time and the attainment of the data are treated later on this page.

Facts & Assumptions

Given: Countable Choice, c>0, u0∈C3(R2), u1∈C2(R2), the weighted ball integral W of Spherical means and the weighted ball integral of space-dependent data, and the extension Uj(ξ,z):=uj(ξ) of uj to R3.

[F1]

In three dimensions the Kirchhoff expression ∂t[tMU0(3)((ξ,z),ct)]+tMU1(3)((ξ,z),ct) is a C2 solution of Vtt=c2Δ3V (Kirchhoff's formula in three dimensions).

[F2]

For even n, tn−1MG(n+1)(x,0,ct)=(n−1)!!cn−1Wg(x,ct) for the cylindrical extension G(ξ,z)=g(ξ), where M(n+1) is the (n+1)-dimensional spherical mean (Sphere integrals of a cylindrical function project to weighted ball integrals).

[F3]

Wf(x,ct)=(n!!Vn)−1∫Bct(x)f(y)(c2t2−∣y−x∣2)−1/2dy for even n, with n!!=n(n−2)⋯2 and Vn the unit-ball volume (Spherical means and the weighted ball integral of space-dependent data); in particular, the weight w(z)=(1−∣z∣2)−1/2 is integrable on B1⊂R2 by the polar-coordinate integrability statement in that definition.

[F4]

If measurable functions converge pointwise almost everywhere and are dominated by one nonnegative integrable function, their integrals converge (Dominated convergence).

[F7]

A real function continuous on a closed interval and differentiable on its interior has a difference quotient equal to a derivative at an interior point (The mean value theorem, as the case g(x)=x of Cauchy's: for f continuous on [a,b] with a<b and differentiable on (a,b) there is c∈(a,b) with f(b)−f(a)=f′(c)(b−a)).

Proof

1.1F1algebra

Descent. Extend u0,u1 to U0,U1 on R3 by Uj(ξ,z):=uj(ξ); these are C3 respectively C2. By [F1] the function V(ξ,z,t):=∂t[tMU0(3)((ξ,z),ct)]+tMU1(3)((ξ,z),ct) is C2 on R3×(0,∞) with Vtt=c2Δ3V. Since Uj(ξ+ctω′,z+ctω3)=uj(ξ+ctω′) does not depend on z, the sphere means of U0,U1 at centre (ξ,z) are independent of z; hence V is independent of z, so ∂z2V=0, Δ3V=Δ2V, and the restriction u(x,t):=V(x,0,t) is a C2 function on R2×(0,∞) with utt=c2Δ2u.

1.2F2F3algebra

The weighted-ball form. By [F2] with n=2 and G=Uj, t MUj(3)((x,0),ct)=1!!cWuj(x,ct)=1cWuj(x,ct); therefore u(x,t)=1c[∂tWu0(x,ct)+Wu1(x,ct)]. With n=2, n!!Vn=2π, [F3] reads Wf(x,ct)=12π∫Bct(x)f(y)(c2t2−∣y−x∣2)−1/2dy, so u is the displayed Poisson expression; the derivative ∂t acts on the C1 function t↦Wu0(x,ct), since Wu0(x,ct)=ctMU0(3)((x,0),ct) by [F2], and the spherical mean is C3, and not by differentiating a singular integrand.

1.3F3F4F5F6F7algebra

The integrated equivalent form. For fixed x put w(z):=(1−∣z∣2)−1/2 on B1⊂R2 and A(x,t):=∫B1u0(x+ctz)w(z) dz. By [F3], w∈L1(B1). To differentiate in t, fix compact sets K⊂R2 and J⊂(0,∞) for x and t. Choose a closed ball Q containing x+c(t+s)z for x∈K, t∈J, ∣s∣ sufficiently small and ∣z∣≤1. By [F6], C:=sup⁡Q∣∇u0∣<∞. The difference quotients of u0(x+ctz) in t, for these parameters and sufficiently small increments h, are bounded in absolute value by cC∣z∣ by [F7]. After multiplication by w(z) they are dominated by the locally uniform integrable function cC∣z∣w(z). They converge pointwise to c∇u0(x+ctz)⋅z w(z), so [F4] gives ∂tA(x,t)=c∫B1∇u0(x+ctz)⋅z w(z) dz. The same argument for spatial difference quotients, and dominated convergence applied to convergent parameter sequences with the same compact-set majorant, shows that these first derivatives are continuous locally; in particular A is C1 in (x,t). Now [F5] gives I(x,t):=∫Bct(x)u0(y)(c2t2−∣y−x∣2)−1/2dy=ctA(x,t), hence ∂tI=cA+ct∂tA. The same change of variables gives ∫Bct(x)∇u0(y)⋅(y−x)(c2t2−∣y−x∣2)−1/2dy=(ct)2∫B1∇u0(x+ctz)⋅z w(z) dz=(ct)2∂tA/c. Therefore 12πc[∂tI+∫Bct(x)u1(y)(c2t2−∣y−x∣2)−1/2dy]=12πct∫Bct(x)u0(y)+∇u0(y)⋅(y−x)+tu1(y)c2t2−∣y−x∣2dy, which is the equivalent form.

2.1given∎

Both displays therefore define the same C2 solution of the two-dimensional homogeneous wave equation on R2×(0,∞); the limit at t↓0 and the attainment of the data are the subject of the data-attainment lemma below.

TheoremStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

The odd-dimensional wave formula by iterated spherical means

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let n=2k+1≥3 be odd, c>0, u0∈Ck+2(Rn), u1∈Ck+1(Rn), put Dt:=t−1∂t, and let M be the spherical mean of Spherical means and the weighted ball integral of space-dependent data. Then u(x,t):=1(n−2)!![∂∂tDtk−1(tn−2Mu0(x,ct))+Dtk−1(tn−2Mu1(x,ct))] defines a C2 function on Rn×(0,∞) solving utt=c2Δu; each Dtj is applied to the t-dependent function t↦tn−2Mf(x,ct). For n=3 (k=1) this is exactly Kirchhoff's formula Kirchhoff's formula in three dimensions, and the displayed constant (n−2)!! is the one required by the leading coefficient of Radial-derivative expansion of the Euler–Poisson–Darboux transform and its zero-radius limit.

Facts & Assumptions

Given: Countable Choice, n=2k+1≥3, c>0, u0∈Ck+2(Rn), u1∈Ck+1(Rn), and the spherical means Hf(x,r):=Mf(x,r).

[F1]

For j≥1 and every φ∈Cj+1((0,∞)), ∂r2Drj−1(r2j−1φ(r))=Drj−1[r2j−1r−2j∂r(r2j∂rφ(r))] on (0,∞), where Dr=r−1∂r (The iterated radial-derivative identity behind the odd-dimensional reduction).

[F2]

For m≥1 and h∈Cm(Rn), (x,r)↦Mh(x,r) is Cm on Rn×(0,∞) with all derivatives obtained by differentiating h under the sphere integral (Smoothness, parity and zero-radius limits of spherical means).

[F3]

For h∈C2(Rn) the Euler–Poisson–Darboux identity ΔxMh(x,r)=Mrr(x,r)+n−1rMr(x,r) holds for every r>0 (The Euler–Poisson–Darboux equation for spherical means).

[F4]
[F5]

If h is C2 on an open subset of Rm, then ∂i∂jh=∂j∂ih (Clairaut--Schwarz theorem for continuous second partial derivatives).

Proof

1.1F1F2

Fix f∈Ck+1(Rn) and put W(x,t):=Dtk−1(tn−2Hf(x,ct)), where n−2=2k−1. Applying [F1] with j=k and φ(t):=Hf(x,ct) — admissible for each fixed x because Hf(x,⋅)∈Ck+1 by [F2] — gives ∂t2Dtk−1(t2k−1Hf(x,ct))=Dtk−1[t2k−1t−2k∂t(t2k∂tHf(x,ct))].

1.2F3F4algebra

The inner expression. By [F4], ∂tHf(x,ct)=c (Hf)r(x,ct) and ∂t2Hf(x,ct)=c2(Hf)rr(x,ct), so with r=ct and 2k=n−1, t−1∂t(t2k∂tHf(x,ct))=t−1∂t(tn−1c(Hf)r(x,ct))=c tn−3[(n−1)(Hf)r(x,ct)+ct(Hf)rr(x,ct)]=c tn−3[(n−1)(Hf)r(x,r)+r(Hf)rr(x,r)]=c2tn−2ΔxHf(x,ct), the last equality because r((Hf)rr+n−1r(Hf)r)=rΔxHf by [F3].

1.3F5algebra

Hence ∂t2W(x,t)=c2Dtk−1(tn−2ΔxHf(x,ct)). The operator Δx acts only on the x-variables while Dtk−1 and multiplication by tn−2 act only on t, and the mixed partials involved commute by [F5] since Hf is Ck+1; therefore Dtk−1(tn−2ΔxHf(x,ct))=ΔxDtk−1(tn−2Hf(x,ct))=ΔxW(x,t), that is ∂t2W=c2ΔxW on Rn×(0,∞).

1.4F2algebra

Regularity and superposition. For f=u0∈Ck+2(Rn) the function t↦tn−2Hu0(x,ct) is Ck+2 by [F2], so Wu0=Dtk−1(⋅) is C3; for f=u1∈Ck+1(Rn) similarly Wu1 is C2. Hence u=(n−2)!!−1[∂tWu0+Wu1] is C2 on Rn×(0,∞) and utt=(n−2)!!−1[∂t∂t2Wu0+∂t2Wu1]=(n−2)!!−1[∂t(c2ΔxWu0)+c2ΔxWu1]=c2Δxu.

2.1given∎

For k=1, that is n=3, the formula reads ∂t[tMu0(x,ct)]+tMu1(x,ct) with (n−2)!!=1, which is Kirchhoff's expression of Kirchhoff's formula in three dimensions. The prefactor (n−2)!! is the leading coefficient of Drk−1(r2k−1φ) by Radial-derivative expansion of the Euler–Poisson–Darboux transform and its zero-radius limit: with r=ct and with the evenness of the means, which kills the first-order term, that expansion makes the normalised combination the data-carrying normalisation; the precise attainment of u0 and u1 is proved by the data-attainment lemma below.

TheoremStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

The even-dimensional wave formula by descent

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let n=2k≥2 be even, c>0, u0∈Ck+2(Rn), u1∈Ck+1(Rn), let W be the weighted ball integral of Spherical means and the weighted ball integral of space-dependent data and put Dt=t−1∂t. Then u(x,t):=c1−n[∂∂tDtk−1Wu0(x,ct)+Dtk−1Wu1(x,ct)] defines a C2 function on Rn×(0,∞) solving utt=c2Δu. For n=2 (k=1) this is exactly Poisson's formula Poisson's formula in two dimensions by descent, and the factor c1−n is the exact rescaling of the unit-speed formula obtained by substituting s=ct in the (n+1)-dimensional odd-dimensional formula.

Facts & Assumptions

Given: Countable Choice, n=2k≥2, c>0, u0∈Ck+2(Rn), u1∈Ck+1(Rn), and the cylindrical extensions Uj(ξ,z):=uj(ξ) to Rn+1.

[F1]

For odd m=2k+1 and data in Ck+2 respectively Ck+1, the formula of The odd-dimensional wave formula by iterated spherical means with Dt=t−1∂t defines a C2 solution of vtt=c2Δv on Rm×(0,∞).

[F2]

For even n, tn−1MU(n+1)((x,0),ct)=(n−1)!!cn−1Wu(x,ct) for the cylindrical extension of u, where M(n+1) is the spherical mean in n+1 variables (Sphere integrals of a cylindrical function project to weighted ball integrals).

[F3]

For n=2 and k=1 the formula reduces to Poisson's formula Poisson's formula in two dimensions by descent with Wf(x,ct)=12π∫Bct(x)f(y)(c2t2−∣y−x∣2)−1/2dy (Spherical means and the weighted ball integral of space-dependent data).

Proof

1.1F1algebra

Descent. The number n+1=2k+1 is odd, and U0∈Ck+2(Rn+1), U1∈Ck+1(Rn+1), so [F1] applies in dimension n+1 with the same k: V(ξ,z,t):=1(n−1)!![∂tDtk−1(tn−1MU0(n+1)((ξ,z),ct))+Dtk−1(tn−1MU1(n+1)((ξ,z),ct))] is C2 on Rn+1×(0,∞) with Vtt=c2Δn+1V. Since Uj(ξ+ctω,z+ctωn+1)=uj(ξ+ctω) is independent of z, the means of U0,U1 at centre (ξ,z) do not depend on z; hence V does not depend on z, so ∂z2V=0, Δn+1V=ΔnV, and the restriction u(x,t):=V(x,0,t) is a C2 function on Rn×(0,∞) with utt=c2Δnu.

1.2F2F4algebra

Rewriting with the weighted ball integral. By [F2] with the cylindrical extensions, tn−1MUj(n+1)((x,0),ct)=(n−1)!!cn−1Wuj(x,ct) for j=0,1; substituting into the definition of V and using [F4] to move the constant (n−1)!!cn−1 and the factor 1(n−1)!! through the t-derivatives gives u(x,t)=c1−n[∂tDtk−1Wu0(x,ct)+Dtk−1Wu1(x,ct)].

1.3F3algebra

The two-dimensional case. For n=2, k=1 and c1−n=c−1, so the formula reads c−1[∂tWu0(x,ct)+Wu1(x,ct)], which by [F3] is exactly the Poisson expression of Poisson's formula in two dimensions by descent.

2.1given∎

Therefore the displayed even-dimensional formula defines a C2 solution of the homogeneous wave equation, it reduces to Poisson's formula when n=2, and its prefactor c1−n is the one produced by substituting s=ct in the (n+1)-dimensional odd formula.

LemmaStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

The dimension formulas attain the Cauchy data

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let c>0, n≥2 and let u be one of the functions constructed from data u0,u1 of the regularity required by the corresponding formula: Kirchhoff (n=3, Kirchhoff's formula in three dimensions), Poisson (n=2, Poisson's formula in two dimensions by descent), the odd-dimensional formula (The odd-dimensional wave formula by iterated spherical means) or the even-dimensional formula (The even-dimensional wave formula by descent). Then, as t↓0, for every x u(x,t)⟶u0(x),∂tu(x,t)⟶u1(x), and the extension u(x,0):=u0(x) is continuous on Rn×[0,∞).

Facts & Assumptions

Given: Countable Choice, c>0, n≥2, and one of the four representation formulas with its data classes.

[F1]

For j≥1 and every φ∈Cj+1((0,∞)), Drj−1(r2j−1φ(r))=∑i=0j−1αj,iri+1φ(i)(r) with αj,0=(2j−1)!! (Radial-derivative expansion of the Euler–Poisson–Darboux transform and its zero-radius limit).

[F2]

For m≥1 and h∈Cm(Rn), the spherical mean Mh is Cm on Rn×(0,∞) with ∂riMh obtained by differentiating h under the sphere integral, and r↦Mh(x,r) is even (Smoothness, parity and zero-radius limits of spherical means).

[F3]

The odd- and even-dimensional formulas define C2 solutions of utt=c2Δu on Rn×(0,∞) (The odd-dimensional wave formula by iterated spherical means, The even-dimensional wave formula by descent); for n=3 the odd formula is Kirchhoff's expression and for n=2 the even formula is Poisson's expression (Kirchhoff's formula in three dimensions, Poisson's formula in two dimensions by descent).

Proof

1.1F1F2F4algebra

Differentiated finite expansion. Suppose n=2k+1 and put hf(x,t)=Mf(x,ct), using the signed-radius extension of [F2]. This is Ck+1 and even when f∈Ck+1, so hf(x,0)=f(x) and ∂thf(x,0)=0. By [F1], Tf(x,t):=Dtk−1(t2k−1hf(x,t))=∑j=0k−1αk,jtj+1∂tjhf(x,t), where a:=αk,0=(2k−1)!!. Differentiate this finite sum itself: Tf′=∑jαk,j((j+1)tjhf(j)+tj+1hf(j+1)) and Tf′′=∑jαk,j(j(j+1)tj−1hf(j)+2(j+1)tjhf(j+1)+tj+1hf(j+2)), with the first summand omitted for j=0. Every derivative used has order at most k+1. Continuity and hf′(x,0)=0 give Tf→0, Tf′→af(x) and Tf′′→0, uniformly for x in compact sets. For Tf′′, the only terms without a positive power of t are constant multiples of hf′, which vanish at zero. These formulas never differentiate an unspecified error term.

2.1F3step 1.1algebra

Odd-dimensional data. The odd formula is u=a−1(Tu0′+Tu1) and ut=a−1(Tu0′′+Tu1′). Step 1.1 applies to both data, since their classes are at least Ck+1. It follows that u→u0 and ut→u1, locally uniformly in x. The even smooth signed-radius means in step 1.1 also show that Tf and its first two derivatives extend continuously through zero.

3.1F3step 2.1

Even-dimensional data by descent. For n=2k, extend the data cylindrically to Rn+1. Their differentiability classes are exactly those of the odd formula in dimension n+1=2k+1. The construction in The even-dimensional wave formula by descent identifies the even solution with the restriction of that odd solution to the last coordinate zero. The limits of step 2.1 therefore apply without differentiating a singular ball weight.

4.1F3step 2.1step 3.1∎

The locally uniform displacement limit and continuity of u0 give joint continuity of the extension u(x,0)=u0(x). The velocity limit holds as stated. The cases n=3 and n=2 are Kirchhoff and Poisson by [F3].

TheoremStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

Duhamel's principle for the wave equation

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let c>0, n≥2, T>0 and let f:Rn×[0,T]→R be continuous, compactly supported in x for each fixed s, with all spatial derivatives Dxαf through order q:=⌊n/2⌋+1 existing and jointly continuous on Rn×[0,T]. Thus every slice is in the velocity-data class, and the additional derivatives needed to differentiate the launched solution twice are jointly continuous. No time derivative of f is required. For an admissible velocity datum g let W[g](x,τ) be the homogeneous solution constructed from the formulas above with zero displacement and velocity datum g, so that, by The dimension formulas attain the Cauchy data, W[g](x,0)=0 and ∂τW[g](x,0)=g(x). Then u(x,t):=∫0tW[f(⋅,s)](x,t−s) ds(0≤t≤T) is a C2 function with u(⋅,0)=ut(⋅,0)=0 and utt=c2Δu+f on Rn×(0,T).

Facts & Assumptions

Given: Countable Choice, c>0, n≥2, a source f of the stated class, and for each admissible g the launched solution W[g] with W[g](x,0)=0, ∂τW[g](x,0)=g(x).

[F1]

The formulas of the page define C2 solutions of the homogeneous equation on Rn×(0,∞) for admissible data, and the data are attained in the limit sense (The dimension formulas attain the Cauchy data).

[F2]

Let α,β∈C1(I) with α<β and let F be continuous on I×J with continuous ∂tF, where J contains the closure of the union of the intervals [α(t),β(t)]. Then G(t)=∫α(t)β(t)F(t,y) dy is C1 with G′(t)=F(t,β(t))β′(t)−F(t,α(t))α′(t)+∫α(t)β(t)∂tF(t,y) dy (Differentiating an integral with moving endpoints).

[F3]

In odd dimension n=2k+1, a zero-displacement launch is W[g](x,τ)=a−1∑j=0k−1αk,jτj+1∂τjMg(x,cτ), a=(2k−1)!!, by Radial-derivative expansion of the Euler–Poisson–Darboux transform and its zero-radius limit. The signed-radius mean is Ck+1 (Smoothness, parity and zero-radius limits of spherical means). Differentiating this finite sum through total order two involves at most k+1 spatial derivatives of g and nonnegative powers of τ. For g=f(⋅,s), the uniform integral estimate in the smoothness lemma applies also with the continuous parameter s: the assumed joint continuity on compact spatial-time sets makes W and its first two (x,τ) derivatives jointly continuous, including at τ=0. In even dimension use the cylindrical launch in dimension n+1 from The even-dimensional wave formula by descent, with the same derivative order q=k+1.

Proof

1.1F1F2F3

First derivative. The launched solutions vanish at τ=0 and have velocity g there by [F1]: W[g](x,0)=0 and ∂τW[g](x,0)=g(x). Applying [F2] to the moving-endpoint integral u(x,t)=∫0tW[f(⋅,s)](x,t−s) ds in the form G(t)=∫0tF(t,s) ds with F(t,s)=W[f(⋅,s)](x,t−s) — defined also for negative t−s by the signed-radius finite sum in [F3]; extend the source slices constantly for s<0 and s>T. Then F and its first two t-derivatives are continuous on a rectangular neighbourhood of the integration region by [F3] — gives ut(x,t)=∫0t∂tW[f(⋅,s)](x,t−s) ds+W[f(⋅,t)](x,0)=∫0t∂tW[f(⋅,s)](x,t−s) ds, since W[⋅](x,0)=0; in particular u(⋅,0)=0 and ut(⋅,0)=0.

1.2F1F2F3algebra

Second derivative. Differentiating once more with [F2], utt(x,t)=∂tW[f(⋅,t)](x,0)+∫0t∂t2W[f(⋅,s)](x,t−s) ds=f(x,t)+∫0tc2ΔxW[f(⋅,s)](x,t−s) ds=f(x,t)+c2Δxu(x,t), where the last equality uses ∂τW[f(⋅,t)](x,0)=f(x,t) and the homogeneous equation for every launched solution from [F1]. To justify moving Δx through the integral, fix any compact set of x-values and a compact time interval [0,T0]⊆[0,T]. By [F3], the integrand and its first two x-derivatives are jointly continuous on the resulting compact (x,s,t−s) parameter set; applying [F2] twice with the fixed s-interval endpoints therefore permits differentiating under the s-integral locally in x. The mixed derivative utx is obtained similarly from the integral formula for ut. These derivative integrals and their boundary terms are continuous; their bounds on local compact sets give continuous one-sided derivatives also at t=0,T. No common compact support of all source slices is needed.

2.1given∎

Hence u is C2 with zero Cauchy data and utt−c2Δu=f on Rn×(0,T); the constructed u is a classical solution of the forced problem.

TheoremStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

The forced three-dimensional version as a retarded potential

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let c>0 and let f be continuous on R3×[0,∞) with f(⋅,t)∈C2(R3) and compactly supported in x for each t. Assume that ∇xf and every second spatial partial derivative ∂xi∂xjf are jointly continuous in (x,t); no time derivatives of f are required. Then the Duhamel construction gives a classical solution of utt=c2Δu+f with zero Cauchy data, namely u(x,t)=14πc2∫Bct(x)f(y, t−∣y−x∣/c)∣y−x∣ dy(t>0), the retarded potential over the backward light cone of (x,t). In particular the value uses f only on {(y,s):0≤s≤t, ∣y−x∣=c(t−s)}, and the radius factor is 1/(4πc2).

Facts & Assumptions

Given: Countable Choice, c>0, a source f of the stated class, and the launched Kirchhoff solutions with zero displacement.

[F1]

For admissible sources the Duhamel principle gives the forced solution as u(x,t)=∫0tW[f(⋅,s)](x,t−s) ds, where W[g] is the homogeneous solution with zero displacement and velocity datum g (Duhamel's principle for the wave equation).

[F2]

In three dimensions W[g](x,τ)=τMg(x,cτ), and the mean is the normalised sphere integral with ∣S2∣=ω2=4π: Mg(x,ρ)=1ω2∫S2g(x+ρz) dσ(z)=14πc2τ2∫∂Bcτ(x)g dS, where ω2=3V3=4π follows from V3=π3/2/Γ(5/2) and the Gamma values (Kirchhoff's formula in three dimensions, The dimension formulas attain the Cauchy data, Sphere and ball measures scale in Rn, The closed form for the volume of the unit n-ball, The real Gamma functional equation Γ(s+1)=sΓ(s), Γ(1/2)=π from the Gaussian integral).

[F3]

Under Countable Choice, ∫R3F dλ3=∫0∞∫S2F(ρz)ρ2 dσ(z) dρ for every Borel F≥0 (Polar coordinates decompose Lebesgue measure into r^{n-1} dr d sigma).

Proof

1.1F1F2

Duhamel form. By [F1] the solution is u(x,t)=∫0tW[f(⋅,s)](x,t−s) ds, and by [F2] the launched solution is W[f(⋅,s)](x,τ)=τMf(⋅,s)(x,cτ); substituting τ=t−s gives u(x,t)=∫0t(t−s) Mf(⋅,s)(x,c(t−s)) ds.

1.2F2F4algebra

Sphere-integral form. Writing the mean over S2 as the normalised integral with ω2=4π from [F2], u(x,t)=14π∫0t(t−s)∫S2f(x+c(t−s)z, s) dσ(z) ds; substituting ρ=c(t−s), so s=t−ρ/c, (t−s)=ρ/c and ds=−dρ/c, gives u(x,t)=14πc2∫0ctρ∫S2f(x+ρz, t−ρ/c) dσ(z) dρ.

2.1F3step 1.2algebra

Ball form. The source is bounded on the compact backward cone by A continuous real-valued function on a nonempty compact metric space is bounded and attains a greatest and a least value and For n≥1, every Euclidean closed ball and every Euclidean sphere of positive radius is compact. The weight ∣y−x∣−1 is integrable on Bct(x), since [F3] gives its integral as 4π∫0ctρ dρ<∞. Give the integrand any value at y=x, a null singleton. Apply [F3] separately to the positive and negative parts after translating by x (Lebesgue outer measure, Lebesgue measurability and Lebesgue measure are unchanged by translation). This gives ∫Bct(x)f(y,t−∣y−x∣/c)∣y−x∣−1dy=∫0ctρ∫S2f(x+ρz,t−ρ/c) dσ(z) dρ. Multiplication by 1/(4πc2) identifies this with step 1.2.

3.1F1algebra∎

The integrand is evaluated at ∣y−x∣=c(t−s) with s=t−∣y−x∣/c∈[0,t], that is on the backward light cone of (x,t), and the coefficient is 1/(4πc2); this is the retarded potential.

TheoremStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

Sphere-supported versus interior-supported free wave kernels

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let c>0 and let u be the free solution assigned to admissible data (u0,u1) by the formulas of the page (Kirchhoff, Poisson, odd- and even-dimensional formulas). (i) If n≥3 is odd and t>0, the value u(x,t) depends on the displacement and velocity data only through their restrictions to a neighbourhood of the sphere ∂Bct(x): if two admissible data pairs agree on such a neighbourhood, their solutions agree at (x,t). Pointwise agreement only on the sphere is not asserted. (ii) If n≥2 is even, t>0 and u1 is supported in a compact subset of the open ball Bct(x), then the velocity contribution has the kernel form c1−nDtk−1Wu1(x,ct)=∫Bct(x)u1(y) K(∣y−x∣,t) dy,K(ρ,t):=c1−nDtk−1[1n!!Vn(c2t2−ρ2)−1/2] for 0≤ρ<ct, where K is smooth on that region and K(0,t)=a t−(n−1) with a=c−n(−1)k−1(2k−3)!!/(n!!Vn)≠0 (read (2k−3)!!=1 for k=1); consequently some data supported strictly inside Bct(x) give a nonzero velocity contribution, so the kernel fills the interior of the ball rather than sitting on the sphere. For displacement data also supported in a compact subset of the open ball, the displacement term has the corresponding kernel ∂tK(ρ,t).

Facts & Assumptions

Given: Countable Choice, c>0, even n=2k≥2, and velocity data u1 supported in a compact subset of Bct(x).

[F1]

The odd-dimensional formula expresses u(x,t) as a finite combination of t-derivatives of t↦tn−2Mu0(x,ct) and t↦tn−2Mu1(x,ct) (The odd-dimensional wave formula by iterated spherical means).

[F2]

For m≥1 and h∈Cm(Rn), (x,r)↦Mh(x,r) is Cm with all derivatives obtained by differentiating h under the sphere integral (Smoothness, parity and zero-radius limits of spherical means).

[F3]

The even-dimensional formula reads u=c1−n[∂tDtk−1Wu0(x,ct)+Dtk−1Wu1(x,ct)], with Wf(x,ct)=(n!!Vn)−1∫Bct(x)f(y)(c2t2−∣y−x∣2)−1/2dy (The even-dimensional wave formula by descent, Spherical means and the weighted ball integral of space-dependent data).

[F4]

If G(y,s) is continuous with continuous ∂sG on a compact rectangle, then s↦∫G(y,s) dy is C1 with derivative ∫∂sG; iterating gives the higher derivatives when they are continuous (Leibniz's rule on a compact rectangle: an interior parameter derivative with a continuous extension may be passed through a Riemann integral).

Proof

1.1F1F2algebra

Part (i). By [F1] the value u(x,t) is computed from the functions r↦Mu0(x,r) and r↦Mu1(x,r) and their r-derivatives at r=ct, and by [F2] these are obtained by differentiating the defining sphere integrals. For r in a neighbourhood of ct, Mf(x,r) is the average of f over ∂Br(x), a sphere contained in the chosen neighbourhood of ∂Bct(x); hence all these quantities depend on f only through its restriction to that neighbourhood, and so does u(x,t).

1.2F3F4algebra

Part (ii), kernel form. Fix t>0 and let u1 vanish on a neighbourhood of ∂Bct(x); the assumed compact support lies inside the open ball, so there is ε∈(0,t) with supp⁡u1⊆Bc(t−ε)(x). For ∣s−t∣<ε/2 the integrand u1(y)(c2s2−∣y−x∣2)−1/2 and all its s-derivatives are continuous on a fixed box containing the support, with the product integrand extended by zero where u1=0, for s∈[t−ε/2,t+ε/2], so by [F4] applied successively in each coordinate of the box, with Fubini (Fubini's theorem for L^1 functions on a sigma-finite product), the derivative Dtk−1 passes under the integral sign: c1−nDtk−1Wu1(x,ct)=∫Bct(x)u1(y)K(∣y−x∣,t) dy with K(ρ,t)=c1−nDtk−1[(n!!Vn)−1(c2t2−ρ2)−1/2]. For displacement data u0 with the same compact-interior support condition, one additional time differentiation under the fixed-box integral gives the kernel ∂tK(∣y−x∣,t).

1.3F5algebra

The interior kernel. The identity Dt(c2t2−ρ2)−p=−2pc2(c2t2−ρ2)−p−1 gives by induction K(ρ,t)=b(c2t2−ρ2)−(2k−1)/2, where b=c1−n(−1)k−1(2k−3)!!c2k−2/(n!!Vn)≠0 and (−1)!!=1 for k=1. Thus K has a fixed nonzero sign throughout 0≤ρ<ct. In particular K(0,t)=c−n(−1)k−1(2k−3)!!t−(n−1)/(n!!Vn), the stated value.

2.1F3step 1.2step 1.3algebra

Admissible interior data. Choose 0<δ<ct/2 and, by A smooth bump between concentric Euclidean balls translated to centre x, choose u1∈Cc∞(B2δ(x)) with 0≤u1≤1 and u1=1 on B‾δ(x). This is admissible in every dimension here. The actual contribution is ∫B2δ(x)u1(y)K(∣y−x∣,t) dy, including the transition annulus. By step 1.3 its integrand has one sign and is strictly of that sign on the inner ball, of positive volume, so the contribution is nonzero. The same construction around any point strictly inside the ball shows that the interior kernel is nonzero throughout, rather than just at the centre.

3.1given∎

Collecting: in odd dimensions the value depends only on data near the sphere ∂Bct(x) (a statement about open neighbourhoods, not about pointwise traces), while in even dimensions the velocity kernel is the explicitly displayed smooth function of ρ<ct, nonzero at the centre, so interior data contribute.

CorollaryStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

The constructed classical solutions are locally determined by the Cauchy data

Statement

Assume the Axiom of Countable Choice (The Axiom of Countable Choice (ACω)). Let c>0 and let u be a free solution constructed by the formulas of Kirchhoff's formula in three dimensions, Poisson's formula in two dimensions by descent, The odd-dimensional wave formula by iterated spherical means or The even-dimensional wave formula by descent from compactly supported admissible data. Then for every (x,t) with t>0 the value u(x,t) is determined by the data restricted to B‾ct(x): two admissible data pairs agreeing there produce the same value at (x,t). For the forced solution of The forced three-dimensional version as a retarded potential, the value is determined by the source on the backward cone {(y,s):0≤s≤t, ∣y−x∣≤c(t−s)}; two sources agreeing there produce the same value at (x,t).

Facts & Assumptions

Given: Countable Choice, c>0, a point (x,t) with t>0, and two admissible configurations agreeing on the stated set.

[F1]

The four free formulas express u(x,t) through the spherical means Muj(x,ct) and their r-derivatives (odd case) or through Wuj(x,ct) and its t-derivatives with the substitution y=x+ctz (even case), and the forced solution is the retarded potential (The odd-dimensional wave formula by iterated spherical means, The even-dimensional wave formula by descent, Kirchhoff's formula in three dimensions, Poisson's formula in two dimensions by descent, The forced three-dimensional version as a retarded potential).

[F2]

For m≥1 and h∈Cm(Rn), the spherical mean Mh and its derivatives are obtained by differentiating h under the sphere integral (Smoothness, parity and zero-radius limits of spherical means).

[F4]

Difference quotients converging pointwise almost everywhere under one integrable majorant have convergent integrals (Dominated convergence).

[F5]

A scalar function continuous on a closed interval and differentiable on its interior has a difference quotient equal to one of its derivatives on the interior (The mean value theorem, as the case g(x)=x of Cauchy's: for f continuous on [a,b] with a<b and differentiable on (a,b) there is c∈(a,b) with f(b)−f(a)=f′(c)(b−a)).

Proof

1.1F1F2F3F4F5F6algebra

Free case. Let (u0,u1) and (u~0,u~1) be admissible data agreeing on B‾ct(x) and put δj:=uj−u~j. Each δj vanishes on the open ball, so all its derivatives through the orders in the formulas vanish on the closed ball by continuity. In the odd-dimensional formula, every sphere-average ingredient is an average of a derivative of δj evaluated on ∂Bct(x), hence is zero by [F1, F2]. For even n, put Gδj(t′):=∫B1δj(x+ct′z)w(z) dz. Fix a compact interval J about t and a closed ball containing all x+ct′z for t′∈J, ∣z∣≤1. For each derivative order 0≤m<k, the difference quotients in t′ of cmDmδj(x+ct′z)[z,…,z]w(z) are bounded by cm+1Cm+1w(z) on J, where Cm+1 bounds the next derivative on that ball by [F5, F6]. This is an integrable majorant independent of t′; [F4] therefore justifies differentiating under the integral successively through order k. At t′=t, all integrand derivatives vanish because x+ctz∈B‾ct(x), so Gδj(m)(t)=0 for m≤k. By [F3] the even-formula terms are finite combinations of these derivatives and hence vanish. Thus replacing the data by (u~0,u~1) changes no term of the formula and leaves u(x,t) unchanged.

2.1F1given∎

Forced case and conclusion. The retarded potential of the forced three-dimensional formula is an integral of the source over the backward cone ∣y−x∣≤c(t−s), 0≤s≤t; sources agreeing there give equal integrals, hence equal values at (x,t). This proves the local determination of the constructed solutions by the stated data or source; no uniqueness claim for arbitrary C2 solutions is made, that being the energy statement of the wave-energy page.

CorollaryStatement: Literature-sourcedProof: AI-adaptedjudge pass (gpt-6.1-sol)Open item page →

Time reversal of the homogeneous wave equation

Statement

Let c>0, let I⊆R be an open interval and let u∈C2(Rn×I) satisfy □cu=0 on Rn×I. For τ∈I define v(x,t):=u(x,2τ−t) for t∈2τ−I. Then □cv=0 on Rn×(2τ−I), and v(⋅,τ)=u(⋅,τ),∂tv(⋅,τ)=−∂tu(⋅,τ). Thus the homogeneous wave flow is reversible: the same equation propagates the time-reversed state, and the velocity is negated.

Facts & Assumptions

Given: an open interval I, a parameter τ∈I, a C2 function u on Rn×I with □cu=0, and v(x,t)=u(x,2τ−t) on Rn×(2τ−I).

[F1]

If g is totally differentiable at a and f at g(a), then f∘g is totally differentiable at a with D(f∘g)(a)=Df(g(a))∘Dg(a) (The chain rule for total derivatives: D(g∘f)(a)=Dg(f(a))∘Df(a)). The required total differentiability follows from continuous coordinate partial derivatives (If all partial derivatives exist on a neighbourhood and are continuous at a point, then the map is totally differentiable there with Jacobian derivative).

Proof

1.1F1F2algebra

Apply [F1] to the composition (x,t)↦(x,2τ−t)↦u(x,2τ−t): the inner map is affine with differential (h,s)↦(h,−s), so ∂tv(x,t)=−∂tu(x,2τ−t); applying the same rule once more, ∂t2v(x,t)=(−1)2∂t2u(x,2τ−t)=∂t2u(x,2τ−t). In the spatial directions the inner map is the identity, so ∂xjv(x,t)=∂xju(x,2τ−t) for each j and hence Δv(x,t)=Δu(x,2τ−t).

2.1F2algebra∎

Therefore for every (x,t) with t∈2τ−I, □cv(x,t)=∂t2u(x,2τ−t)−c2Δu(x,2τ−t)=□cu(x,2τ−t)=0, while at t=τ the substitutions 2τ−τ=τ give v(x,τ)=u(x,τ) and ∂tv(x,τ)=−∂tu(x,τ). This is the reversibility of the homogeneous wave flow.

RemarkRemark: Literature-sourcedProof: Not applicablejudge pass (gpt-6.1-sol)Open item page →

Two different objects are called Poisson's formula

Remark

The name "Poisson's formula" denotes two different objects that should not be conflated. The two-dimensional wave formula of Poisson's formula in two dimensions by descent is the weighted disk integral u(x,t)=12πc[∂∂t∫Bct(x)u0(y) dyc2t2−∣y−x∣2+∫Bct(x)u1(y) dyc2t2−∣y−x∣2], whose kernel (2πc)−1(c2t2−∣y−x∣2)−1/2 is an interior singular weight on the expanding disk, while the harmonic Poisson kernel of Poisson kernel of a Euclidean ball integrates boundary data against a positive density on the sphere. The two solve different problems — an initial-value (Cauchy) problem in space-time versus the Dirichlet boundary-value problem — use respectively a time parameter t and a fixed radius R. The shared name does not identify their kernels or transfer estimates between these problems.

The wave formula is stated for the speed-c convention of Wave equation, Cauchy data and wave speed; the harmonic kernel is the ball boundary-value kernel of the Poisson-problem page. The remark asserts no new mathematics: it isolates the naming collision so that no consumer imports a boundary-value estimate into the wave representation.

5 · Examples, counterexamples and false statements

None yet.

Sources