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.

Em Causal Response And Radiation Prerequisites

1 · Prerequisites

2 · Summary

The mathematical layer beneath causal response theory and compact-source radiation. The first item, Causal convolution response, transfer function and dispersion relations, develops the convolution response P(t)=∫0∞K(s)E(t−s)ds of a real matrix kernel on bounded continuous histories, proves that it is bounded, continuous, linear and causal, and analyzes its transfer function χ(z)=∫0∞K(s)eizsds in the compactly supported and integrable branches: the principal-value dispersion relations hold entrywise under exact decay hypotheses, and the monochromatic period-averaged work is identified with the quadratic form of the Hermitian dissipative part (χ−χ∗)/(2i), whose positivity is a passivity condition on the model rather than a consequence of causality. The second item, Uniform far-zone asymptotics and retarded-time flux limit for compact conserved smooth sources, treats smooth sources of compact spatial support with compact time support satisfying the continuity equation; it proves the retarded field integrals on the exterior of the source, identifies the retarded potentials with the zero-data wave Cauchy solution, establishes the uniform far-zone expansion with explicit error control by finitely many source suprema, and derives the outward flux limit at fixed retarded time for the sphere at infinity.

Both items are mathematics: they carry domain: mathematics, depend only on mathematical suppliers, and assert no physical conclusion. Empirical qualifications of any physical use of these estimates belong to the consuming physical item. The companion examples page records the two hypothesis boundaries that these proofs depend on, namely that causality alone does not force a nonnegative work form and that the far-zone flux limit is taken at fixed retarded time rather than fixed observation time.

3 · Logical flowchart

4 · Definitions, postulates, theorems and proofs

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

Causal convolution response, transfer function and dispersion relations

Statement

Assume Countable Choice (The Axiom of Countable Choice (ACω)) for the smooth Riemann/Lebesgue comparison route below; in clause (c), also assume Dependent Choice (The axiom of dependent choice: a relation in which every element is related to something admits an N-indexed chain) inherited from the AC supplier. Let d≥1. All kernels below are real d×d matrix kernels and all histories are Rd-valued; conditions on matrices are read entrywise, every matrix- or vector-valued integral below is read entrywise as a scalar Lebesgue integral, and L1(0,∞;Rd×d) means entrywise Lebesgue integrability on (0,∞) with Lebesgue measure.

For a real-variable scalar function F and real ω, the principal value is defined by PV⁡ ⁣∫RF(ξ)ξ−ω dξ:=lim⁡R→∞lim⁡η↓0(∫ω−Rω−η+∫ω+ηω+R)F(ξ)ξ−ω dξ, read entrywise when applied to matrix entries; under the bounds proved in (b) and (c) the tails converge absolutely, so this limit exists and the centering at infinity does not affect it.

(a) (Causal convolution.) Let K∈L1(0,∞;Rd×d) and let E:R→Rd be bounded and continuous. Then for every t∈R the integral P(t):=∫0∞K(s)E(t−s) ds converges absolutely, read componentwise as Pj(t)=∫0∞(∑k=1dKjk(s)Ek(t−s))ds(j=1,…,d), and satisfies ∣Pj(t)∣≤∑k=1d(sup⁡s∈R∣Ek(s)∣)∫0∞∣Kjk(s)∣ ds(j=1,…,d); the response P is bounded and continuous, the map E↦P is linear, and it is causal: if E1,E2 are bounded continuous and coincide on (−∞,t0], then their responses coincide on (−∞,t0].

(b) (Compact branch.) The support of a scalar, vector or matrix map is the closure, in its domain, of the set where it is nonzero. If K∈Cc∞((0,∞);Rd×d) (smooth, with compact support contained in the open half-line), then the transfer function χ(z):=∫0∞K(s)eizs ds is entire, satisfies χ(−ω)=χ(ω)‾ for every real ω, obeys ∣χjk(z)∣≤∥Kjk′′∥L1(0,∞)∣z∣−2 for every z≠0 with Im⁡z≥0, and for any b>0 with supp⁡K⊆(0,b) the full entire-plane bound is ∣χjk(z)∣≤ebmax⁡(0,−Im⁡z)∥Kjk′′∥1∣z∣−2 for every z≠0. Its entries satisfy the principal-value dispersion relations Re⁡χjk(ω)=1πPV⁡ ⁣∫RIm⁡χjk(ξ)ξ−ω dξ,Im⁡χjk(ω)=−1πPV⁡ ⁣∫RRe⁡χjk(ξ)ξ−ω dξ for every real ω.

(c) (Integrable branch.) If instead K:[0,∞)→Rd×d is locally absolutely continuous entrywise in the sense of Absolute continuity on a compact interval on every compact subinterval of [0,∞), with K′ its almost-everywhere derivative with K,K′∈L1(0,∞;Rd×d) and s↦sK(s) in L1(0,∞;Rd×d), then χ(z):=∫0∞K(s)eizs ds is well defined and holomorphic on the open upper half-plane {Im⁡z>0}, continuous on the closed upper half-plane with continuous first derivative on the real boundary, satisfies the same pair of dispersion relations for every real ω, and is bounded on {Im⁡z≥0} by ∣χjk(z)∣≤∣Kjk(0)∣+∥Kjk′∥L1(0,∞) for ∣z∣≥1; for an instantaneous real response matrix χ∞∈Rd×d define χtotal:=χ∞+χ and apply the relations to χtotal−χ∞=χ.

(d) (Monochromatic work and passivity.) Let ω>0, e∈Cd, and let ε0>0 be a real parameter (a response normalization; the equivalence below is independent of its value). Put E(t):=Re⁡(e e−iωt) and P(t):=Re⁡(ε0χ(ω)e e−iωt), real parts taken componentwise and χ(ω) acting by matrix multiplication. Then, with T:=2π/ω and the Euclidean dot product on Rd, 1T∫0TE(t)⋅P′(t) dt=ε0ω2Im⁡(e∗χ(ω)e),e∗χ(ω)e:=∑j,k=1dej‾ χjk(ω) ek, where χ(ω)∗ denotes the conjugate transpose, (χ(ω)∗)jk:=χkj(ω)‾, a complex matrix H is Hermitian when H∗=H, and a Hermitian H is positive semidefinite when e∗He≥0 for every e∈Cd. This period-averaged work is nonnegative for every e∈Cd if and only if the Hermitian dissipative part M(ω):=(χ(ω)−χ(ω)∗)/(2i) satisfies e∗M(ω)e≥0 for every e∈Cd.

Facts & Assumptions

Given: The Axiom of Countable Choice ACω (The Axiom of Countable Choice (ACω)), and in clause (c) additionally Dependent Choice DC (The axiom of dependent choice: a relation in which every element is related to something admits an N-indexed chain); a fixed dimension d≥1; a real d×d matrix kernel K=(Kjk) and an Rd-valued history E=(Ek); all matrix- and vector-valued integrals are read entrywise as scalar Lebesgue integrals on (0,∞). In clause (b), K∈Cc∞((0,∞);Rd×d) with supp⁡K⊆(0,b) for some real b>0; in clause (c) the kernel satisfies the integrable-branch hypotheses of [F14]; in clause (d), ω>0, T=2π/ω, ε0>0, e∈Cd, and a:=ε0χ(ω)e∈Cd.

[F1]

A measurable h=u+iv:Ω→C is integrable exactly when ∫∣h∣ dμ<+∞, and then ∫h dμ=∫u dμ+i∫v dμ; L1(μ) is the class of integrable complex functions. (Integrable real and complex functions, and their integrals, The class L1(μ) of integrable functions)

[F2]

For integrable functions ∣∫f dμ∣≤∫∣f∣ dμ; for measurable 0≤f≤g one has ∫f≤∫g and ∫cf=c∫f for c≥0; the integral over a measurable set is ∫Ef=∫fχE; and the integral is complex-linear, ∫(αf+βg) dμ=α∫f dμ+β∫g dμ. (The modulus of an integral is bounded by the integral of the modulus, Monotonicity and nonnegative homogeneity of the nonnegative integral, Integral over a measurable subset, The Lebesgue integral is linear on L1(μ))

[F3]

Every continuous map Rn→Rm is Borel measurable (under ACω), and sums, scalar multiples, absolute values and pointwise products of measurable functions are measurable. (Continuous functions on Euclidean spaces are Borel measurable, Arithmetic and lattice operations preserve measurability whenever they are defined)

[F4]

Dominated convergence: if fn→f almost everywhere and ∣fn∣≤g almost everywhere for a single integrable g, then ∫fn→∫f. (Dominated convergence)

[F5]

Continuity under the integral sign: if x↦f(x,t) is integrable for every t in an interval I, if t↦f(x,t) is continuous for almost every x, and if ∣f(x,t)∣≤g(x) for a single integrable g, then t↦∫f(x,t) dμ(x) is continuous on I. (Continuity under the integral sign)

[F6]

Differentiation under the integral sign in a real parameter: under the hypotheses stated there (integrability for every parameter, almost-everywhere differentiability, measurability of the derivative and a single integrable dominating function for the derivative), F(t)=∫f(x,t) dμ(x) is differentiable with F′(t)=∫∂tf(x,t) dμ(x). (Differentiation under the integral sign)

[F7]

The complex exponential is defined by exp⁡z=∑n≥0zn/n!, that series converges absolutely for every z, exp⁡ is entire with exp⁡′=exp⁡ and hence continuous, and ∣exp⁡(x+iy)∣=ex for real x,y; conjugation is a real-field automorphism with zz‾=∣z∣2, hence ∣z‾∣=∣z∣ and conjugation is continuous, while Re⁡, Im⁡ and z‾ are the coordinate maps of Real and imaginary parts, complex conjugation, and modulus. In particular Re⁡(eiθ)=cos⁡θ and Im⁡(eiθ)=sin⁡θ for real θ, and Re⁡(w‾)=Re⁡w. (The complex exponential by its power series, The complex exponential series converges absolutely for every complex argument, The complex exponential is entire and its complex derivative is itself, Complex differentiability at a point implies continuity there, exp⁡(x+iy)=ex(cos⁡y+isin⁡y), ∣exp⁡(x+iy)∣=ex, and eiπ+1=0, Conjugation is an involutive real-field automorphism, zz‾=∣z∣2, and modulus is definite, multiplicative, and subadditive, Real and imaginary parts, complex conjugation, and modulus)

[F8]

If the four first partial derivatives of u,v exist on a neighbourhood of a and are continuous at a with ux(a)=vy(a) and uy(a)=−vx(a), then f=u+iv is complex differentiable at a; holomorphic means complex differentiable on an open set and entire means holomorphic on C. (Continuous first partial derivatives and the Cauchy–Riemann equations imply complex differentiability pointwise, and holomorphy when they hold throughout an open set, Complex differentiability at a point, the complex derivative, holomorphic functions, and entire functions)

[F9]

If f is holomorphic on an open Ω⊆C and the closed rectifiable contour γ is path-homotopic relative to endpoints to a constant loop in Ω, then ∫γf dz=0. (A closed contour path-homotopic to a constant loop has zero integral against every holomorphic function)

[F10]

Integration by parts: if u,v are differentiable at every point of [a,b] with integrable derivatives, then ∫abu v′=u(b)v(b)−u(a)v(a)−∫abu′v; the bounded Riemann integrable case has the same Lebesgue integral; and integrable functions that agree almost everywhere have the same integral. (If u,v are differentiable on [a,b] with u′,v′ integrable, then ∫abuv′=u(b)v(b)−u(a)v(a)−∫abu′v, A bounded Riemann integrable function on a closed bounded interval is Lebesgue measurable and has the same integral, Two integrable functions are equal almost everywhere exactly when all of their indefinite integrals agree, Ck Euclidean maps and diffeomorphisms)

[F12]

On Rd the Euclidean inner product is x⋅y=∑k=1dxkyk (The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn). On Cd its complex-bilinear extension is also denoted x⋅y=∑kxkyk; this bilinear form is not used as a positive inner product. The conjugated pairing used below is x⋅y‾=∑kxkyk‾, and e∗a=∑kek‾ak, with conjugation as in [F7].

[F13]

Absolute continuity on a compact interval is the finite disjoint-interval condition of the cited definition; it is the exact hypothesis by which clause (c) reads its kernel, and no differentiability of the kernel beyond the stated almost-everywhere derivative is assumed. (Absolute continuity on a compact interval)

[F14]

Integrable-branch supplier, consumed verbatim in clause (c): assume ACω and DC; for a locally absolutely continuous real matrix kernel K with K,K′∈L1(0,∞;Rd×d) and sK(s)∈L1(0,∞;Rd×d), the transform χ(z)=∫0∞K(s)eizsds has holomorphic entries on the upper half-plane, continuous and C1 on the real boundary, is bounded on the closed upper half-plane with ∣χjk(z)∣≤∣Kjk(0)∣+∥Kjk′∥1 for ∣z∣≥1 (indeed ∣χjk(z)∣≤(∣Kjk(0)∣+∥Kjk′∥1)/∣z∣ for every nonzero z with Im⁡z≥0), satisfies the two displayed principal-value formulas for every real ω, and carries the instantaneous-response subtraction for a real matrix χ∞; causality alone supplies neither further growth hypotheses nor the passivity quadratic-form condition. (Dispersion relations for an absolutely continuous integrable causal kernel)

[F15]

Period-average supplier, consumed in clause (d) in its scalar form applied componentwise: for ω>0, T=2π/ω and all a,d∈C, 1T∫0TRe⁡(ae−iωt)Re⁡(de−iωt) dt=12Re⁡(ad‾); the vector form follows by summing this over the coordinate pairs of a,d∈Cd. (Period averages of products of real parts of complex exponentials)

Proof

technique · direct; the steps are ordered in dependency layers, so the four clauses of the Statement are interleaved, and every step cites only earlier steps
1.1F1F2F3

Fix t∈R, j∈{1,…,d} and set gj(s):=∑k=1dKjk(s)Ek(t−s) for s>0. Each summand is a product of the measurable function Kjk with the continuous, hence Borel measurable, function s↦Ek(t−s); so gj is measurable. With M:=max⁡ksup⁡R∣Ek∣<+∞ one has ∣gj(s)∣≤M∑k=1d∣Kjk(s)∣ for every s>0, and ∑k=1d∣Kjk∣ lies in L1(0,∞); hence gj∈L1(0,∞) and Pj(t):=∫0∞gj(s) ds converges absolutely. The termwise bound ∣gj(s)∣≤∑k=1d(sup⁡R∣Ek∣)∣Kjk(s)∣ gives the exact stated estimate ∣Pj(t)∣≤∑k=1d(sup⁡R∣Ek∣)∫0∞∣Kjk(s)∣ ds≤M∑k=1d∫0∞∣Kjk(s)∣ ds.

1.2F1F3F7

Fix z=x+iy∈C. The map s↦eizs is continuous, since it is the composition of the continuous linear map s↦izs with the continuous exponential, and Kjk is measurable; hence s↦Kjk(s)eizs is measurable as a product of measurable functions. On supp⁡K⊆(0,b) one has ∣eizs∣=∣eixs∣e−ys=e−ys≤eb∣y∣, so ∣Kjk(s)eizs∣≤eb∣y∣∣Kjk(s)∣; the right-hand side is integrable, so the integrand is integrable and χjk(z):=∫0∞Kjk(s)eizs ds is well defined with ∣χjk(z)∣≤eb∣y∣∥Kjk∥1.

1.3F14F13F8

Clause (c) is exactly the supplied integrable-branch theorem [F14]: its hypotheses on K are those stated in clause (c) (local absolute continuity entrywise in the sense of [F13], K,K′,sK∈L1), and it concludes the well-definedness and holomorphy of χ on the open upper half-plane, continuity on the closed upper half-plane and continuous first derivative on the real boundary, the entrywise principal-value formulas for every real ω, and the bound ∣χjk(z)∣≤∣Kjk(0)∣+∥Kjk′∥1 for ∣z∣≥1 (which follows from its sharper ∣z∣−1 bound), together with the instantaneous-response subtraction for a real matrix χ∞ and the scalar case d=1. Continuity on the closed upper half-plane combines the lemma's holomorphy on the open half-plane, hence continuity there, with its continuity at every real boundary point.

1.4F7F12

Write aj=αj+iβj and ej=pj+iqj for the components of a=ε0χ(ω)e and e. By [F7], eiωt=cos⁡(ωt)+isin⁡(ωt), and Re⁡(w)=Re⁡(w‾) for every w∈C; therefore Ej(t)=Re⁡(eje−iωt)=Re⁡(ej‾eiωt)=pjcos⁡(ωt)+qjsin⁡(ωt) and Pj(t)=Re⁡(aje−iωt)=Re⁡(aj‾eiωt)=αjcos⁡(ωt)+βjsin⁡(ωt), using the component computation Re⁡((u+iv)(cos⁡+isin⁡))=ucos⁡−vsin⁡.

2.1step 1.1F2F1

The estimate of step 1.1 is uniform in t, so sup⁡t∈R∣Pj(t)∣≤M∑k=1d∥Kjk∥1 and P=(P1,…,Pd) is bounded. The map E↦P is linear: for α,β∈R and bounded continuous E,F the integrand defining PαE+βF equals α times that of PE plus β times that of PF pointwise, so the linearity of the scalar integral gives PαE+βF=αPE+βPF.

2.2step 1.1F5F3

Fix j and apply [F5] with X=(0,∞) carrying Lebesgue measure, I=R and f(s,t):=∑k=1dKjk(s)Ek(t−s): for every t the function f(⋅,t) is integrable by step 1.1; for every s the map t↦f(s,t) is continuous, being a finite sum of products of constants with the continuous functions t↦Ek(t−s); and ∣f(s,t)∣≤M∑k=1d∣Kjk(s)∣ for the single integrable function named in step 1.1. Therefore t↦Pj(t) is continuous on R, and P is continuous componentwise.

2.3step 1.1F1

Causality: if E1,E2 are bounded continuous and agree on (−∞,t0], and if t≤t0, then for every s>0 the argument t−s lies in (−∞,t0], so the integrands of step 1.1 agree pointwise and hence PE1(t)=PE2(t); thus the response at time t depends only on the history up to time t.

2.4step 1.2F6F4F3

Fix y and apply [F6] in the parameter x to f(s,x):=Kjk(s)ei(x+iy)s on a bounded interval of x-values: for each x the function f(⋅,x) is integrable by step 1.2, the map x↦Kjk(s)ei(x+iy)s is differentiable for every s with derivative isKjk(s)ei(x+iy)s, and on the bounded parameter set this derivative is dominated by (1+b)eb∣y∣∣Kjk(s)∣, an integrable function. Hence x↦χjk(x+iy) is differentiable with ∂xχjk(z)=∫0∞isKjk(s)eizs ds; exchanging the roles of x and y and repeating the same argument gives ∂yχjk(z)=∫0∞(−s)Kjk(s)eizs ds. Both partials are continuous in (x,y): if (xn,yn)→(x,y), the integrands isKjk(s)ei(xn+iyn)s converge pointwise to isKjk(s)eizs and are dominated on a bounded parameter neighbourhood by (1+b)ebY∣Kjk(s)∣ for a suitable Y, so [F4] gives convergence of the integrals.

2.5step 1.2F1F2F7

Reality identity: for real ω, conjugating the partial sums of the absolutely convergent series defining eiωs term by term gives ∑n≤m(iωs)n/n!‾=∑n≤m(−iωs)n/n!, because conjugation is a continuous real-field automorphism with i‾=−i and fixes the real numbers 1/n!; passing to the limit with continuity of conjugation gives eiωs‾=e−iωs. Since Kjk is real, Kjk(s)eiωs‾=Kjk(s)e−iωs. For h:=Kjk(⋅)eiω⋅∈L1(0,∞) one has ∫h‾=(∫Re⁡h)−i(∫Im⁡h)=∫(Re⁡h−iIm⁡h)=∫h‾ by [F1] and the complex linearity of the integral; hence χjk(−ω)=∫0∞Kjk(s)e−iωsds=∫h‾=∫h‾=χjk(ω)‾.

2.6step 1.2F10F2F7

By-parts identity: fix z≠0. On [0,b] the real and imaginary parts of s↦Kjk(s) and of v(s):=eizs/(iz) are differentiable with continuous derivatives, and Kjk(0)=Kjk(b)=Kjk′(0)=Kjk′(b)=0 because supp⁡K⊆(0,b); applying [F10] first with u=Kjk and then with u=Kjk′ gives ∫0bKjk(s)eizsds=−1iz∫0bKjk′(s)eizsds=1(iz)2∫0bKjk′′(s)eizsds=−z−2∫0bKjk′′(s)eizsds, the boundary terms vanishing at both steps. All three integrands are bounded and Riemann integrable on [0,b], so [F10] transfers the identity to Lebesgue integrals; and since Kjk, Kjk′ and Kjk′′ all vanish off (0,b), the a.e.-equality clause of [F10] replaces each integral over [0,b] by the integral over (0,∞). Hence χjk(z)=−z−2∫0∞Kjk′′(s)eizsds for every z≠0.

2.7step 1.4F11F7F12

Differentiating the expressions of step 1.4 and using [F11] together with the algebra of derivatives gives Pj′(t)=ω(−αjsin⁡(ωt)+βjcos⁡(ωt)); on the other hand, expanding (−iω)(αj+iβj)=ωβj−iωαj and multiplying by e−iωt=cos⁡(ωt)−isin⁡(ωt) shows that the real part of (−iω)aje−iωt equals ωβjcos⁡(ωt)−ωαjsin⁡(ωt), the same number. Hence P′(t)=Re⁡((−iω)a e−iωt)=Re⁡((−iω)ε0χ(ω)e e−iωt) componentwise, and P′ is the real part of a monochromatic complex expression with vector amplitude d:=(−iω)ε0χ(ω)e.

3.1step 2.4F2F8

By step 2.4 and the linearity of the integral, ∂yχjk(z)=∫0∞(−s)Kjk(s)eizs ds=∫0∞i(isKjk(s)eizs)ds=i ∂xχjk(z) for every z. Writing χjk=u+iv, this says uy+ivy=i(ux+ivx)=−vx+iux, that is uy=−vx and vy=ux; the four partials exist everywhere and are continuous by step 2.4, so [F8] makes χjk complex differentiable at every point of C, that is, entire.

3.2step 2.6step 1.2F4F2F7

Bounds and boundary continuity: for z≠0 with Im⁡z≥0, step 2.6 and ∣eizs∣≤1 for s≥0 give ∣χjk(z)∣≤∣z∣−2∥Kjk′′∥1; for arbitrary z≠0 one has ∣eizs∣=e−sIm⁡z≤ebmax⁡(0,−Im⁡z) on supp⁡K⊆(0,b), whence ∣χjk(z)∣≤ebmax⁡(0,−Im⁡z)∥Kjk′′∥1∣z∣−2. Also χjk is continuous on the closed upper half-plane: if zn→z with Im⁡zn≥0, the integrands of step 1.2 converge pointwise and are dominated by ebY∣Kjk∣ once all Im⁡zn lie below a bound Y, so [F4] gives χjk(zn)→χjk(z).

3.3step 2.7F15F7F2F12

Applying the period-average identity of [F15] with amplitude vectors a=e and d=(−iω)ε0χ(ω)e componentwise, and using step 2.7, 1T∫0TE(t)⋅P′(t) dt=∑j=1d1T∫0TRe⁡(eje−iωt)Re⁡(dje−iωt) dt=∑j=1d12Re⁡(ejdj‾)=12Re⁡(e⋅d‾). Since d‾=iωε0χ(ω)e‾ and e⋅χ(ω)e‾=∑j,kejχjk(ω)‾ek‾=e∗χ(ω)e‾, this equals ε0ω2Re⁡(i e∗χ(ω)e‾)=ε0ω2Im⁡(e∗χ(ω)e).

4.1step 3.1step 3.2F9F4F2

Dispersion relations: fix real ω and j,k, and put φ(z):=χjk(z)/(z−ω) for z≠ω. For R>∣ω∣+1 and 0<δ<1 let γ be the closed piecewise C1 contour running along [ω−R,ω−δ], then along the clockwise semicircle z=ω+δeiθ, π≥θ≥0, then along [ω+δ,ω+R], and then back along the counterclockwise large semicircle z=ω+Reiθ, 0≤θ≤π. The contour lies in the slit plane D:=C∖{ω−iv:v≥0}, where φ is holomorphic because χjk is entire and z↦(z−ω)−1 is holomorphic off ω. This domain is star-shaped about ω+iA for every A>0: along a segment from this centre to a point off the removed ray, either its imaginary part stays positive or its real part differs from ω except at the centre. The straight contraction to the centre, with the base-point connector and its reverse, therefore gives a null homotopy relative to the base point. Hence [F9] gives ∫γφ dz=0. The large arc contributes at most π∥Kjk′′∥1/(R−∣ω∣)2→0 by step 3.2; the small arc contributes −iπχjk(ω)+o(1) as δ↓0 because χjk(ω+δeiθ)→χjk(ω) uniformly in θ; and the two real segments converge to the principal value of ∫Rχjk(ξ)/(ξ−ω)dξ as first δ↓0 and then R→∞: near ω the term χjk(ω)/(ξ−ω) cancels between the symmetric punctured intervals while (χjk(ξ)−χjk(ω))/(ξ−ω) stays bounded, and for ∣ξ∣≥max⁡(1,2∣ω∣) the estimate ∣χjk(ξ)/(ξ−ω)∣≤2∥Kjk′′∥1∣ξ∣−2 makes the tails absolutely integrable. Passing to the limits in the identity 0=∫γφ dz therefore yields PV⁡∫Rχjk(ξ)/(ξ−ω)dξ=iπχjk(ω); taking real and imaginary parts gives the two displayed relations, since −πIm⁡χjk(ω) and πRe⁡χjk(ω) are respectively the real and imaginary parts of iπχjk(ω).

5.1step 3.3F2F7F12∎

Finally M(ω)∗=(χ(ω)∗−χ(ω))/(−2i)=M(ω), so M(ω) is Hermitian, and M(ω)=(χ(ω)−χ(ω)∗)/(2i) has quadratic form e∗M(ω)e=(e∗χ(ω)e−e∗χ(ω)∗e)/(2i), and e∗χ(ω)∗e=∑j,kej‾χkj(ω)‾ek=e∗χ(ω)e‾ because χ(ω)∗ is the conjugate transpose; hence e∗M(ω)e=Im⁡(e∗χ(ω)e) for every e∈Cd. By step 3.3 the period-averaged work equals ε0ω2e∗M(ω)e with ε0ω2>0, so it is nonnegative for every e if and only if e∗M(ω)e≥0 for every e; clause (a) establishes causality without assuming this positivity condition.

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

Uniform far-zone asymptotics and retarded-time flux limit for compact conserved smooth sources

Statement

Assume Countable Choice (The Axiom of Countable Choice (ACω)), as used by the declared polar/surface and change-of-variables suppliers. Let c>0, μ0>0, a>0 and t0<T be constants, let ρ∈C∞(R3×R) and J∈C∞(R3×R;R3) satisfy (H1) ρ(y,v)=J(y,v)=0 whenever ∣y∣>a or v∉[t0,T], and (H2) ∂vρ+div⁡yJ=0 on R3×R. Here and below Ba:={y∈R3:∣y∣<a} and B‾a:={y∈R3:∣y∣≤a} are the open and closed Euclidean balls of radius a about the origin, and S2={n∈R3:∣n∣=1} is the unit sphere. On the open exterior U:={x∈R3:∣x∣>a} define the retarded scalar and vector potentials φ(x,t):=μ0c24π∫R3ρ(y,t−∣x−y∣/c)∣x−y∣ dy,A(x,t):=μ04π∫R3J(y,t−∣x−y∣/c)∣x−y∣ dy, and the fields E:=−∇φ−∂tA, B:=∇×A, where ∇f=(∂1f,∂2f,∂3f), div⁡F=∑j=13∂jFj and ∇×F=(∂2F3−∂3F2,∂3F1−∂1F3,∂1F2−∂2F1) are the spatial operators of The Jacobian matrix of partial derivatives and the gradient in the scalar-valued case and Divergence and curl of a C1 vector field. The fields are defined on U×R using these ordinary derivatives; the conclusions below use r:=∣x∣≥2a. Every displayed spatial integrand is evaluated on B‾a and extended by zero outside it, including at y=x, so no quotient 0/0 is used. All vector-valued integrals are componentwise scalar Lebesgue integrals.

(a) (Fields as retarded integrals.) On the open set {∣x∣>2a} the potentials φ,A are of class C∞ in (x,t), and there E(x,t)=μ04π∫R3[NR2c2ρ(y,t−Rc)+c NRρt(y,t−Rc)−Jt(y,t−Rc)R]dy,B(x,t)=−μ04π∫R3N×[J(y,t−Rc)R2+Jt(y,t−Rc)cR]dy, where R:=∣x−y∣, N:=(x−y)/R, and ρt,Jt are the partial derivatives in the second variable. Moreover, on {∣x∣≥2a}×[t0,T], φ and each component of A coincide with the corresponding solution constructed in the wave Cauchy problem of The three-dimensional wave Cauchy problem: Kirchhoff representation, uniqueness, domain of dependence and compact support for zero Cauchy data at t0 and scalar source h=μ0c2ρ, respectively h=μ0Jk.

(b) (Uniform far-zone expansion.) For n∈S2 and u∈[t0,T] put I(n,u):=∫BaJ(y,u+n⋅y/c) dy and V(n,u):=∂uI(n,u)=∫BaJt(y,u+n⋅y/c) dy. Then ∫Baρt(y,u+n⋅y/c) dy=1c n⋅V(n,u), and there are constants ΓE,ΓB, depending only on a,c,μ0 and on the suprema of ∣ρ∣,∣ρt∣,∣ρtt∣,∣J∣,∣Jt∣,∣Jtt∣ on B‾a×[t0,T], such that for every n∈S2, u∈[t0,T] and r≥2a, with Elead(n,u,r):=μ04πrn×(n×V(n,u)), ∣E(rn,u+rc)−Elead(n,u,r)∣≤ΓEr2,∣B(rn,u+rc)−n×Elead(n,u,r)c∣≤ΓBr2.

(c) (Retarded-time flux limit.) Let S:=μ0−1E×B and let dΩ be the surface measure on S2 (Surface integration on compact C1 hypersurfaces). Then for every u∈[t0,T], lim⁡r→∞∫S2r2 n⋅S(rn,u+rc) dΩ(n)=μ016π2c∫S2∣n×V(n,u)∣2 dΩ(n). This is a limit at fixed retarded time u, not at fixed observation time t; no decay or boundary condition at spatial infinity is assumed beyond the compact spatial support in (H1).

Facts & Assumptions

Given: ACω (The Axiom of Countable Choice (ACω)); constants c>0, μ0>0, a>0, t0<T; sources ρ∈C∞(R3×R), J∈C∞(R3×R;R3) satisfying (H1) and (H2) of the Statement; the retarded potentials and fields of the Statement, with all vector integrals read componentwise. For n∈S2, u∈[t0,T], r≥2a we write x:=rn, t:=u+r/c, R:=∣x−y∣=((r−n⋅y)2+∣y∣2−(n⋅y)2)1/2, N:=(x−y)/R, σ:=n⋅y, vy:=u+σ/c and δy:=(r−σ−R)/c, so that t−R/c=vy+δy; Γρ,Γρt,Γρtt,ΓJ,ΓJt,ΓJtt denote the suprema of ∣ρ∣,∣ρt∣,∣ρtt∣,∣J∣,∣Jt∣,∣Jtt∣ on B‾a×[t0,T], and λ(Ba) the Lebesgue measure of the ball.

[F1]

A measurable h is integrable exactly when ∫∣h∣<+∞, and then the integral is computed from the positive and negative parts or, for complex h=u+iv, as ∫u+i∫v; L1 denotes the integrable class and ∫Ef=∫fχE is the integral over a measurable set. (Integrable real and complex functions, and their integrals, The class L1(μ) of integrable functions, Integral over a measurable subset)

[F2]

The integral is linear, ∣∫f∣≤∫∣f∣, 0≤f≤g measurable implies ∫f≤∫g, ∫cf=c∫f for c≥0, and integrable functions that agree almost everywhere have the same integral. (The Lebesgue integral is linear on L1(μ), The modulus of an integral is bounded by the integral of the modulus, Monotonicity and nonnegative homogeneity of the nonnegative integral, Two integrable functions are equal almost everywhere exactly when all of their indefinite integrals agree)

[F3]

Under ACω the ball Ba is Lebesgue measurable with 0<λ(Ba)<+∞; every continuous map Rn→Rm is Borel measurable; and sums, scalar multiples and pointwise products of measurable functions are measurable. (Euclidean balls have positive finite Lebesgue measure, Continuous functions on Euclidean spaces are Borel measurable, Arithmetic and lattice operations preserve measurability whenever they are defined)

[F4]

Mean value theorem: if ϕ is continuous on [p,q] and differentiable on (p,q) then ϕ(q)−ϕ(p)=ϕ′(ξ)(q−p) for some ξ∈(p,q). (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))

[F5]

Parameter integrals over a fixed compact measurable set of finite measure: if K⊆Rn is compact, measurable and of finite measure, if a<b, and if the real-valued f and ∂tf are continuous on a neighbourhood of K×[a,b], then t↦∫Kf(x,t) dx is differentiable at every t∈(a,b) with derivative ∫K∂tf(x,t) dx and has the corresponding one-sided relative derivatives at the endpoints of [a,b]; the lemma is stated pointwise in its one real parameter and explicitly asserts nothing for moving parameter sets or for singularities crossing K. Applying it once per parameter direction with its parameter renamed to x1,x2,x3 or t (all other variables held fixed), and iterating on the resulting integrals, yields the corresponding partial derivatives of a parameter integral whose integrand is smooth in all variables on a neighbourhood of the fixed K. (Differentiation under the integral sign over a compact measurable set and over a compact C1 surface)

[F7]

A Ck map has continuous word partial derivatives through order k, and C∞ means this for every k (Ck Euclidean maps and diffeomorphisms). Finite sums, products and composites preserve these classes (Ck Euclidean maps are closed under componentwise algebra and composition). On (0,∞), (qα)′=αqα−1 and real powers agree with rational powers (Continuity and derivatives of positive-base real powers, The exponential definition of real powers agrees with the existing rational powers); iterating this rule makes q=q1/2 and 1/q smooth there. Closed bounded Euclidean sets are compact, continuous functions on compact sets are uniformly continuous and attain finite extrema (Heine-Borel in Rn: with the Euclidean metric a subset of Rn is compact if and only if it is closed and bounded, and the proof by bisection uses no choice principle; the same holds on the real line, Heine-Cantor: a continuous map from a compact metric space to any metric space is uniformly continuous, A continuous real-valued function on a nonempty compact metric space is bounded and attains a greatest and a least value).

[F8]

Zero-data wave Cauchy problem of speed c from time t0 with source h∈C∞(R3×[t0,T]) in the extension sense: its solution is u(x,t)=c2∫t0t(t−v)Mc(t−v)h(⋅,v)(x) dv on R3×[t0,T], it is C∞ there, and it satisfies c−2∂t2u−Δu=h; here Mrh(x,v)=14π∫S2h(x+rn,v) dΩ(n). (The three-dimensional wave Cauchy problem: Kirchhoff representation, uniqueness, domain of dependence and compact support)

[F9]

Polar-coordinates formula: for every nonnegative Borel f:R3→[0,∞], ∫R3f(y) dy=∫0∞∫S2f(rω)r2 dσ(ω) dr; here σ is the polar surface measure of the cited formula, which equals the chart surface measure dΩ of [F13] (Agreement with the existing polar sphere measure), so the two symbols are used interchangeably; a signed or complex f is handled by applying the formula to positive and negative parts and subtracting. (Polar coordinates decompose Lebesgue measure into r^{n-1} dr d sigma, Agreement with the existing polar sphere measure)

[F10]

One-dimensional L1 change of variables: for a C1 diffeomorphism T:U→V of open subsets of R and f∈L1(V), ∫Vf=∫Uf∘T ∣T′∣. (A C^1 diffeomorphism satisfies the change-of-variables formula for L^1 functions)

[F11]

Divergence theorem for a bounded C1 Euclidean domain Ω and F∈C1(Ω‾;Rn): ∫Ωdiv⁡F dx=∫∂ΩF⋅ν dS, with the outward normal. (Divergence on a bounded C1 Euclidean domain)

[F12]

The Euclidean norm satisfies the triangle and Cauchy-Schwarz inequalities (Cauchy-Schwarz ∣⟨x,y⟩∣≤∥x∥2∥y∥2 with its equality case, the triangle inequality for ∥⋅∥2, the parallelogram law and polarisation). Cross product and inner product on R3: u×v is bilinear and antisymmetric, n⋅(u×v)=det⁡(n,u,v), and for ∣n∣=1 the identities n⋅(n×w)=0 and n×(n×w)=n(n⋅w)−w hold, so ∣n×(n×w)∣=∣n×w∣. (The cross product in R3, The Euclidean inner product ⟨x,y⟩=∑k<nxkyk on Rn)

[F13]

Surface integral on the compact embedded C1 hypersurface S2: for a nonnegative Borel f the integral ∫S2f dΩ is defined by finitely many regular parametrizations and an ambient partition of unity, the empty surface integrates to 0, the definition is independent of the chosen presentation, the total surface measure ∫S21 dΩ is finite, and the integral is monotone in f; the definition carries the ACω convention already declared. (Surface integration on compact C1 hypersurfaces)

[F14]

C1 change of variables for nonnegative Borel functions: under ACω, for a C1 diffeomorphism T:U→V of open subsets of Rm and every nonnegative Borel h:V→[0,∞], ∫Vh(y) dy=∫Uh(T(z))∣det⁡DT(z)∣ dz with equality in [0,∞]; in particular the translation T(z)=x+z of R3 is a C1 diffeomorphism with det⁡DT=1, and a signed integrand is handled through its positive and negative parts by linearity. (Borel change of variables from the compact-support formula and Radon uniqueness, The Lebesgue integral is linear on L1(μ))

Proof

technique · direct; part (a) differentiates the retarded integrals and identifies them with the wave Cauchy solution, part (b) expands the field integrals uniformly in the far zone and derives the continuity identity, part (c) passes to the surface integral
1.1F5F6F7F3

Smoothness of the potentials: for x∈U and y∈B‾a one has R=∣x−y∣≥∣x∣−a>0, so the maps (x,t,y)↦ρ(y,t−R/c)/R and (x,t,y)↦J(y,t−R/c)/R are C∞ on a neighbourhood of U×R×B‾a by [F6] and [F7]. Apply [F5] to the fixed compact set K=B‾a once for each parameter direction with its one real parameter renamed to x1, x2, x3 and t (the remaining variables held fixed) and iterate on the resulting integrals: every differentiated integrand is again continuous on a neighbourhood of K, so each application is covered by the lemma and φ and each component of A have continuous partial derivatives of every order on U×R, these differentiated parameter integrals are jointly continuous: on each compact parameter neighbourhood the continuous differentiated integrands are uniformly continuous on its product with K, so the integral of their difference is bounded by λ(K) times their uniformly vanishing difference [F2]. Hence the potentials are C∞ on U×R by [F7], so the derivatives defining the fields are ordinary derivatives even at ∣x∣=2a.

1.2F12F1

Elementary estimates in the far zone: for y∈Ba, r≥2a and n∈S2, expanding R2=(r−σ)2+∣y∣2−σ2 with σ=n⋅y gives ∣R−(r−σ)∣=∣∣y∣2−σ2∣/(R+r−σ)≤a2/r because R+r−σ≥r−σ+∣r−σ∣≥r; consequently ∣r−R∣≤a+a2/r, ∣N−n∣=∣(r−R)n−y∣/R≤(a+a2/r+a)/(r−a)≤5a/r, ∣1/R−1/r∣≤(a+a2/r)/(r(r−a))≤3a/r2, and ∣δy∣≤a2/(rc), all uniformly in y,n.

2.1step 1.1F5F6F12

Field formulas: by [F6], ∇xR=N, ∇x(t−R/c)=−N/c and ∇x×[J(y,t−R/c)]=−N×Jt(y,t−R/c)/c, while ∇x(1/R)=−N/R2; hence ∇x[ρ(y,t−R/c)/R]=−N[ρ/R2+ρt/(cR)] and ∇x×[J(y,t−R/c)/R]=−N×[J/R2+Jt/(cR)]. Each of these first derivatives of the integrands is continuous on a neighbourhood of B‾a, so by [F5] the derivatives may be taken under the integral: −∇φ−∂tA and ∇×A are given on U, in particular for ∣x∣≥2a, by the two displayed zero-extended integrals of the statement, since E=−∇φ−∂tA and B=∇×A.

2.2step 1.1F9F14F1F3

Polar-coordinates identity: fix x∈U and t∈R, and let h be either ρ or one component Jk; then y↦h(y,t−∣x−y∣/c)/∣x−y∣ vanishes off B‾a and ∣x−y∣≥∣x∣−a>0 there, so the quotient on B‾a, extended by zero elsewhere, defines a bounded compactly supported Borel function f on all of R3. The translation T(z):=x+z is a C1 diffeomorphism of R3 with det⁡DT=1, so [F14] applied to the positive and negative parts of f gives ∫R3f(y) dy=∫R3f(x+z) dz; applying [F9] to the positive and negative parts of z↦f(x+z), which at z=ρω equals h(x+ρω,t−ρ/c)/ρ and therefore has inner angular integral 1ρ∫S2h(x+ρn,t−ρ/c) dΩ(n)=4πρMρh(x,t−ρ/c), gives 14π∫R3h(y,t−∣x−y∣/c)∣x−y∣dy=∫0∞r Mrh(x,t−r/c) dr with Mrh(x,v)=14π∫S2h(x+rn,v) dΩ(n), because the inner angular integral is 4πMrh and the radial factor is r2/r=r.

2.3step 1.2F11F6F2

Continuity identity: fix n,u and let ρ0>a; the map F(y):=J(y,u+n⋅y/c) is C1 on B‾ρ0, and [F6] gives div⁡F(y)=div⁡yJ(y,vy)+nc⋅Jt(y,vy) with vy=u+n⋅y/c. By [F11] on Bρ0, whose boundary flux vanishes because J=0 for ∣y∣>a, ∫Bρ0div⁡F dy=0; substituting the chain-rule expression and using (H2), 0=∫Bρ0[−ρt(y,vy)+1cn⋅Jt(y,vy)]dy. The integrands vanish for ∣y∣>a by (H1), so [F2] lets the integrals be taken over Ba; the second is 1cn⋅V(n,u) by the definition of V, whence ∫Baρt(y,vy) dy=1cn⋅V(n,u).

3.1step 2.2F10F8F1

Identification with the Cauchy solution: fix x∈U. For t>t0 the radial integrand of step 2.2 vanishes when r>c(t−t0), because then its source time is less than t0. Apply [F10] to the affine diffeomorphism v↦c(t−v) from (t0,t) onto (0,c(t−t0)), with absolute derivative c; the integrand is bounded on this finite interval, hence integrable. Endpoint values do not affect the Lebesgue integrals, and this gives ∫0∞rMrh(x,t−r/c) dr=c2∫t0t(t−v)Mc(t−v)h(x,v) dv. At t=t0 both sides are zero by (H1), so the identity holds throughout [t0,T]. Multiplying it by μ0c2 for h=ρ, and by μ0 for h=Jk, identifies φ and Ak respectively with the zero-data Cauchy solutions of [F8] with sources μ0c2ρ and μ0Jk. This holds on U×[t0,T] and in particular at ∣x∣=2a.

3.2step 1.2step 2.3F2F3F4F12

Uniform error estimates: fix n,u,r, let ret denote the argument vy+δy and 0 the argument vy, and write Ftot:=c2Nρret/R2+cNρt,ret/R−Jt,ret/R and Flead:=1r(c n ρt,0−Jt,0), so that the field formula of step 2.1 and the definition of V give E(rn,u+r/c)−Elead(n,u,r)=μ04π∫Ba(Ftot−Flead) dy by the linearity of the integral. Every source factor is bounded by its supremum on B‾a×[t0,T], and (H1) makes ρ and J vanish on the open set R3×((−∞,t0)∪(T,∞)), hence makes all their partial derivatives vanish there too, so these suprema also dominate the retarded arguments vy+δy; the mean value theorem [F4] gives ∣ρt,ret−ρt,0∣≤Γρtt∣δy∣ and ∣Jt,ret−Jt,0∣≤ΓJtt∣δy∣: when the difference is nonzero, apply [F4] to its scalar projection onto the unit direction of that difference, and use Cauchy-Schwarz from [F12], and using R−1≤2/r, R−2≤4/r2 and ∣N/R−n/r∣≤∣N−n∣/R+∣1/R−1/r∣≤13a/r2 yields ∣Ftot−Flead∣≤ΓE′/r2, with ΓE′:=4c2Γρ+13caΓρt+a2Γρtt+3aΓJt+a2ΓJtt/c for all r≥2a; hence, by the triangle inequality for the integral and the finiteness of λ(Ba), ∣E(rn,u+r/c)−Elead(n,u,r)∣≤ΓE/r2 with ΓE:=μ04πΓE′λ(Ba). The same computation applied to B(rn,u+r/c)=−μ04π∫BaN×[Jret/R2+Jt,ret/(cR)]dy and to n×Elead/c=−μ04π∫Ba1crn×Jt,0 dy gives ∣B(rn,u+r/c)−n×Elead(n,u,r)/c∣≤ΓB/r2, for ΓB:=μ0λ(Ba)4π(4ΓJ+13aΓJt/c+a2ΓJtt/c2).

4.1step 3.2F5F6F7F12

Smoothness of the far-zone moments: the map (n,u)↦I(n,u)=∫BaJ(y,u+n⋅y/c) dy and its partial derivative ∂uI=∫BaJt(y,u+n⋅y/c) dy=V(n,u) are C∞ in (n,u)∈R3×R: apply [F5] repeatedly in the directions n1,n2,n3,u with the fixed compact set B‾a, each differentiated integrand again being continuous on a neighbourhood of B‾a; the integrals over Ba and B‾a agree because the smooth sources and all their derivatives vanish on ∣y∣=a by continuity from ∣y∣>a. In particular n↦∣n×V(n,u)∣2 is continuous for each fixed u.

4.2step 3.2F12F2

Pointwise flux asymptotics: by step 3.2 write E=Elead+ρE and B=n×Elead/c+ρB at (rn,u+r/c) with ∣ρE∣≤ΓE/r2, ∣ρB∣≤ΓB/r2; then n⋅Elead=0 and Elead×(n×Elead)=∣Elead∣2n by [F12], while ∣Elead∣=μ04πr∣n×V(n,u)∣ because ∣n×(n×V)∣=∣n×V∣. Expanding the cross product bilinearly, the error terms are O(r−3), so r2n⋅S(rn,u+r/c)=μ016π2c∣n×V(n,u)∣2+R(n,u,r) with ∣R(n,u,r)∣≤C/r for a constant C depending only on a,c,μ0 and the source suprema, uniformly in n∈S2 and u∈[t0,T].

5.1step 4.2step 4.1F13F2∎

Flux limit: by [F13] the surface integral over the compact sphere of a bounded Borel function is finite and monotone, so for every u∈[t0,T] and r≥2a the difference ∫S2r2n⋅S(rn,u+r/c) dΩ(n)−μ016π2c∫S2∣n×V(n,u)∣2 dΩ(n)=∫S2R(n,u,r) dΩ(n) has absolute value at most Cr∫S21 dΩ, which tends to 0; the function n↦∣n×V(n,u)∣2 is continuous by step 4.1, so its surface integral is the displayed finite constant. Hence the stated limit at fixed retarded time u follows; no step compares different values of u, so the limit is asserted only at fixed retarded time.

5 · Examples, experiments and counterexamples

None yet.

Sources