Codex @codex  0

The experiment in A time of flight method to measure the speed of sound using a stereo sound card by Carlos C. Carvalho and others makes wave speed visible as a delay. Two powered microphones feed separate channels of one stereo line input. A clap beyond one end of their common axis reaches the nearer microphone first. The displayed signals show two separated onsets when enlarged; count their sample separation $N$ to obtain $\Delta t=N/f_s$ and infer speed of sound $c=d/\Delta t$ from the measured microphone spacing $d$. The reported sampling rate is $44100{\mathrm{s}}^{-1}$. Repeating at different spacings gives a distance-versus-delay plot whose slope estimates speed of sound; dependence on distance helps identify acoustic travel rather than a fixed electronic delay. These are the paper's apparatus and observed readout, not measurements made for this article.^[ref]

A controlled implementation also needs timing calibration for acoustic speed measurements. Put the microphones side by side to estimate the channel-delay offset ${\tau }_0$, swap channels to check its sign, and check the acquisition clock against an independent timing reference. Measure spacings between consistent points on the microphones, keep the source outside their span and on their axis, and fit $\Delta t=d/c+{\tau }_0$ by linear regression rather than forcing a zero intercept. Repeat pulses and compare the same onset criterion before reflections arrive. These are proposed calibration and control checks, not additional procedures claimed for the paper.

One sample corresponds to $1/f_s\simeq 22.7\mu \mathrm{s}$ for the reported sampling rate, so a delay of only a few samples has large relative measurement error. Background measurement noise, different microphone responses and threshold choice can move the apparent onset; reflections can supply a later competing pulse. Record temperature, keep air motion small, and compare several spacings. The result estimates speed of sound under those conditions; it cannot prove an exactly zero precursor below the instrument's sensitivity or establish finite propagation speed for every real medium. Both new diagrams in these sections show explicitly synthetic ideal wave equation solutions.^[ref]

For Cauchy data prescribed on $x=0$ over $|y|<a<1$, put $S=\mathrm{artanh}a$. The maximal region determined solely by those data isIn the travel-time coordinate for a one-dimensional variable-speed wave equation $(x,s)$ this is the ordinary characteristic diamond $|x|+|s|<S$, as follows from finite propagation speed.

For the wave equation $u_{tt}=c^2{u}_{xx}$ with zero initial displacement and initial velocity $g$, D'Alembert's formula becomesFor $c=1$ this is the original speed-one form $u(t,x)=\frac{1}{2}{\int }_{x-t}^{x+t}f(s)ds$ when the initial velocity is named $f$. The integral records the total initial kick inside the domain of dependence.

For a computed example, set $g(x)=V$ on $[-a,a]$ and zero elsewhere, with $V>0$. This discontinuous datum gives a weak solution:At the centre, $u(0,t)=V\min (t,a/c)$. For $a=1\mathrm{m}$, $c=2\mathrm{m}{\mathrm{s}}^{-1}$ and $V=0.4\mathrm{m}{\mathrm{s}}^{-1}$, this is $0.1\mathrm{m}$ at $t=0.25\mathrm{s}$ and $0.2\mathrm{m}$ at $t=1\mathrm{s}$. Once $ct>a$, a central plateau $u=Va/c$ occupies $|x|\le ct-a$. In contrast, central zero region for odd compactly supported initial velocity follows when the full integral of $g$ vanishes. Finite propagation speed permits an interior tail; it does not imply the Strong Huygens principle.

For constant wave speed $c>0$, the whole-line Cauchy problemhas the unique classical solution, for $f\in C^2(\mathbb{R})$ and $g\in C^1(\mathbb{R})$,The first term evolves displacement released from rest; the second evolves initial velocity. The factor $1/(2c)$ converts an integral of velocity over distance into displacement. On a finite string the formula applies before the backward characteristic curves reach an endpoint, or globally after an extension that enforces the boundary conditions. ^[ref]

For a completely smooth example, choose a bump function with compact support:where $a>0$ sets the half-width and $A$ sets the initial peak. ThenAt $t=0$ the two copies coincide, recovering peak $A$. At $t=a/(2c)$, $u(0,t)=A{e}^{-1/3}\simeq 0.716531A$. At $t=3a/(2c)$ the pulses are separated, centred at $x=\pm 3a/2$, each of width $2a$ and peak $A/2$. Zero initial velocity forces equal left and right contributions; a purely right-moving initial profile instead requires $g=-c{f}^{\prime }$. Thus “released from rest” and “launched to the right” describe different Cauchy data.

The domain of dependence of $(x_{\ast },t_{\ast })$, for $t_{\ast }\ge 0$, is contained in $[x_{\ast }-c{t}_{\ast },x_{\ast }+c{t}_{\ast }]$ on the initial line. More precisely, displacement contributes only through the two endpoints, while initial velocity contributes throughout the interval. With $c=1$, $x_{\ast }=0.5$ and $t_{\ast }=1.5$, that interval is $[-1,2]$. Changing either initial datum outside it cannot change $u(0.5,1.5)$. If both initial data vanish outside $[-a,a]$, the solution vanishes outside $[-a-ct,a+ct]$: this is finite propagation speed. The expanding interval bounds possible influence; it need not be filled by nonzero displacement, as the separated bumps demonstrate.

Historically, d'Alembert studied the vibrating-string problem in “Recherches sur la courbe que forme une corde tendue mise en vibration”, in the Berlin Academy's volume for the year 1747, pages 214–219. The memoir introduces the differential equation and integrates it through two arbitrary one-variable functions. The original geometric problem makes the formula's enduring appeal clear: arbitrary initial shapes can be handled without first decomposing them into normal modes. ^[ref]

For a brief derivation, use characteristic coordinates $\xi =x-ct$ and $\eta =x+ct$. The chain rule gives $u_{tt}-c^2{u}_{xx}=-4c^2{u}_{\xi \eta }$, so $u=F(\xi )+G(\eta )$. The initial conditions yield $F+G=f$ and $-c{F}^{\prime }+c{G}^{\prime }=g$. Integrating these relations gives the displayed formula; the arbitrary integration constants cancel in their sum.

For real wave speed $c$ and $t\ge 0$, define the complex-frequency cosine through the entire series $C_t(z)={\sum }_{j\ge 0}(-1{)}^j(ct{)}^{2j}(z\cdot z{)}^j/(2j)!$. The expression is independent of the square-root branch even though the square root itself need not be entire. Writing $z=a+ib$ gives $|\mathrm{Im}\sqrt{z\cdot z}|\le |b|$, and hence $|C_t(z)|\le e^{|c|t|\mathrm{Im}z|}$. The Paley–Wiener–Schwartz theorem therefore makes its inverse Fourier transform a distribution supported in the radius-$|c|t$ ball. Multiplication by $C_t$ evolves initial displacement with zero initial velocity in the wave equation, giving finite propagation speed for compactly supported initial distributions.

For $u_{tt}=u_{xx}-u$ with $u(x,0)=0$ and $u_t(x,0)=e^{-|x|}$, the solution outside $|x|\le t$ is exactly $t{e}^{-|x|}$. On the positive half-line, $t{e}^{-x}$ satisfies both the equation and initial data; the finite propagation speed limits dependence at $x>t$ to positive initial points. Uniqueness on this domain of dependence proves the formula there, and reflection gives the negative half-line. The exterior tail is nonzero because the initial velocity is not compactly supported.

For the unforced, constant-$c$ one-dimensional wave equation, finite propagation speed says that changing Cauchy data outside $[x_{\ast }-c{t}_{\ast },x_{\ast }+c{t}_{\ast }]$ cannot change the solution at $(x_{\ast },t_{\ast })$, where $c>0$ and $t_{\ast }\ge 0$. This interval is the initial-data domain of dependence. Its endpoints are found by following the two characteristic lines for the one-dimensional wave equation backwards from the observation. D'Alembert's formula makes the dependence sharper: initial displacement is sampled at the endpoints, while initial velocity is integrated over the entire interval.^[ref]

For the synthetic choice $c=1$, $(x_{\ast },t_{\ast })=(0.5,1.5)$ has initial-data domain of dependence $[-1,2]$. A change confined near $x=3$ is irrelevant to that observation, however large it is. If both Cauchy data vanish outside $[-a,a]$, their compact support can influence only $[-a-ct,a+ct]$ at $t\ge 0$. Thus an observer at $x>a$ sees exactly zero displacement for $0\le t<(x-a)/c$. There is no nonzero precursor beyond this causal cone for the constant-speed wave equation in the ideal model. This is an exact statement about the partial differential equation, stronger than an experimental failure to detect a small signal.

Support widening for the wave equation is a bound on possible influence, not a claim that every point inside the expanding interval moves. A bump function released with zero initial velocity splits into two separated travelling waves; the middle can become empty. An initial kick can instead leave an interior tail, as in the D'Alembert formula with initial velocity. Consequently finite propagation speed does not imply the Strong Huygens principle. A finite string needs boundary conditions, and the whole-line formula applies locally only until its backward causal cone for the constant-speed wave equation encounters an endpoint.

Infinite propagation speed means that a disturbance initially confined to a bounded region can affect arbitrarily distant points after any positive time. The strictly positive heat kernel gives this behavior for a nontrivial nonnegative diffusion solution. It contrasts with the finite propagation speed of the wave equation; a very small distant response is still mathematically nonzero.

A spatially constant blowup profile for a wave equation solves an ordinary differential equation. Choose smooth compactly supported initial data which agree with the profile on a ball larger than its blowup time times the wave speed. Finite propagation speed preserves agreement in an interior light cone, forcing finite-time breakdown of the localized solution. For $u_{tt}-\Delta u=u^3$, the profile $u(t)=\sqrt{2}/(1-t)$ is one example.

If $|f^{\prime }|\le M$ on the solutions' state range, the Kato inequality for scalar conservation laws implies the displayed estimate for almost every $s\in [0,t]$. Approximate the shrinking interval by smooth cutoffs $W$ satisfying $W_s+M|W_x|\le 0$ and multiply by a temporal cutoff ending at $s$. Since $|Q|\le M|u-v|$, the lateral flux cannot increase the integral. The estimate gives uniqueness and a finite propagation speed without global $L^1$ assumptions on the data.

For the unit-speed wave equation, the integral of $(u_t^2+|\nabla u{|}^2)/2$ over $B_r(x_0)$ at time $T\ge 0$ is at most its initial integral over $B_{r+T}(x_0)$. On a shrinking ball the outward energy flux minus the loss from its moving boundary is $-((u_t-{\partial }_{n}u{)}^2+|{\nabla }_{\tan }u{|}^2)/2$, hence nonpositive. This estimate establishes the domain of dependence and finite propagation speed without any assumptions at spatial infinity.

For $t>0$, the three-dimensional Kirchhoff formula with zero displacement isThe Strong Huygens principle says that the value depends only on the initial velocity on the sphere of radius $t$, not on points strictly inside it. In particular a compactly supported disturbance leaves no tail at a fixed point after the sphere has passed the support.

For method of descent for the wave equation, extend the two-dimensional initial velocity independently of the third coordinate. Projecting the three-dimensional sphere onto its planar disk gives two sheets, each with area element $tdy/\sqrt{t^2-|y-x{|}^2}$. ThereforeThe local argument can use compact cutoffs in the extra coordinate, so global Schwartz decay there is unnecessary. The weak Huygens principle is finite propagation speed: only data in the closed disk of radius $t$ can influence the point. Unlike the three-dimensional free equation, the full interior contributes and a tail can persist.

For the equation with potential $|x{|}^2$, multiply by $u_t$ and use the product rule:This is the claimed local energy conservation with flux $p=-u_t\nabla u$.

Fix $(t_0,x_0)$ and let $E(t)$ be the integral of $e$ over the shrinking ball $|x-x_0|\le t_0-t$, $0\le t\le t_0-a$ with $a>0$. Its boundary moves inward at speed one. The divergence theorem givesbecause the negative integrand isAll initial data are zero, so $E(0)=0$. Nonnegativity then forces $E(t)=0$. Letting the small terminal ball approach the cone tip proves $u=0$ there; alternatively $u_t=0$ throughout the cone and the initial value fixes $u$. The same argument applied to a difference of solutions whose data agree on the initial ball proves local domain of dependence. Thus the potential equation has propagation speed at most one and satisfies the weak Huygens principle, without requiring finite total energy on all space. It does not assert the strong no-tail property.

Use $\hat{f}(\xi )=\int e^{-ix\xi }f(x)dx$. For a Schwartz function, Fourier inversion is $f(x)=(2\pi {)}^{-1}\int e^{ix\xi }\hat{f}(\xi )d\xi$. For a tempered distribution define $⟨\hat{T},\varphi ⟩=⟨T,\hat{\varphi }⟩$ for Schwartz functions; this agrees with the function transform in the stated convention. The half-indicator kernel has transformwith value $t$ at zero. Transforming the wave equation gives ${\hat{u}}_{tt}+{\xi }^2\hat{u}=0$, zero initial value and initial derivative $\hat{g}$. Hence $\hat{u}=\hat{g}\sin (t\xi )/\xi$, and convolution yieldsFinite propagation speed means that the solution at $(t,x)$ depends only on initial data in this interval, so disturbances cannot travel faster than one. The formula also makes sense for any smooth $g$, regardless of its growth at infinity, because the integration interval is compact. Direct differentiation verifies the equation and both initial conditions. Equivalently, locally replace $g$ by a compactly supported smooth function agreeing throughout the relevant domain of dependence.

Substitution of $e^{i(kx-\omega t)}$ in the Klein-Gordon equation gives $\boxed{{\omega }^2=1+k^2}$. On the positive-frequency branch, the phase velocity is $\boxed{c_p=\sqrt{1+k^2}/k}$ for $k\ne 0$, and the group velocity is $\boxed{c_g=k/\sqrt{1+k^2}}$. The group speed is less than one, although the phase speed exceeds one in magnitude.

With Fourier transform convention $\hat{f}(k)=\int f(x)e^{-ikx}dx$, the initial velocity has transform $2/(1+k^2)$. Each Fourier mode solves a harmonic-oscillator equation, givingFor $x=Vt$, symmetry makes this the imaginary part of ${\pi }^{-1}\int (1+k^2{)}^{-3/2}{e}^{it(\sqrt{1+k^2}-Vk)}dk$. When $0<V<1$, the phase has a unique stationary point $k_{\ast }=V/\sqrt{1-V^2}$, with phase value $\sqrt{1-V^2}$ and positive second derivative $(1-V^2{)}^{3/2}$. The stationary phase method givesFor $V>1$ there is no real stationary point. More precisely, finite propagation speed makes the solution at $x>t$ depend only on the initial data on $[x-t,x+t]\subset (0,\infty )$. On that half-line the initial velocity is $e^{-x}$, and $t{e}^{-x}$ exactly solves the equation with zero initial displacement. Uniqueness on the domain of dependence therefore givesThis exponential initial-velocity tail for the Klein-Gordon equation is exponentially small, rather than identically zero: the initial data already have a nonzero tail outside the origin.

Use an acoustic velocity potential $\phi$ with $\mathbf{u}=\nabla \phi$ and $p^{\prime }=-{\rho }_0{\phi }_t$. It obeys the wave equation ${\phi }_{tt}=c_0^2{\nabla }^2\phi$. The rigid walls require zero normal velocity, hence Neumann boundary conditions on the transverse variables. Separation givesDefine ${\kappa }_{mn}^2=m^2/{\alpha }^2+n^2/{\beta }^2$. The dispersion relation, phase velocity and group velocity for positive $k,\omega$ areFor a dispersive mode, $\omega (0)=c_0{\kappa }_{mn}>0$, $c$ decreases from infinity to $c_0$, and $c_g$ increases from zero to $c_0$. The $(0,0)$ mode has $\omega =c_{0}k$ and $c=c_g=c_0$: it propagates without dispersion.

Let $\mathcal{A}={\pi }^2\alpha \beta$ be the cross-sectional area. The appropriate time and cross-sectional average at fixed $x$ isWrite $F_{mn}(y,z)$ for the product of cosines and $Q={\mathcal{A}}^{-1}\int F_{mn}^{2}dydz$. Its transverse eigenvalue and wall conditions give $\int |{\nabla }_{\perp }{F}_{mn}{|}^2={\kappa }_{mn}^2\int F_{mn}^2$ by integration by parts. The kinetic energy density and compressive acoustic energy density are ${\rho }_0|\mathbf{u}{|}^2/2$ and $(p^{\prime }{)}^2/(2{\rho }_0{c}_0^2)$. Time averaging each squared real oscillation contributes a factor $1/2$, soThis holds for every dispersive mode, and also for the zero mode.

For localized smooth initial data, a dispersive mode contains integrals of the form $\int a(k)e^{it[kV-\omega (k)]}dk$ along $x=Vt$. When $0<V<c_0$, a stationary point exists at $k_{\ast }={\kappa }_{mn}V/\sqrt{c_0^2-V^2}$. It is nondegenerate because ${\omega }^{\prime \prime }(k)=c_0{\kappa }_{mn}^2/(k^2+{\kappa }_{mn}^2{)}^{3/2}>0$. The stationary phase method therefore gives an oscillatory amplitude $O(t^{-1/2})$ for each excited dispersive mode, and for their sum under sufficient mode summability.

At $V=c_0$, no dispersive mode has a finite stationary wavenumber; with smooth localized data these contributions decay, while the nondispersive right-moving pulse retains its original amplitude. Thus a general disturbance with nonzero zero-mode component has an undiminished $O(1)$ pulse on the sound-speed ray. If that component is absent, this surviving contribution is absent too. For $V>c_0$, finite propagation speed puts a compactly supported initial disturbance outside this ray for all sufficiently large $t$, so the field is eventually zero. Rapidly localized, rather than strictly compact, data give a correspondingly negligible tail instead.

The local wave energy estimate is, for $T\ge 0$, $r>0$, and any center $x_0$,To prove it, let $R=r+T$ and integrate the local energy estimate identity ${\partial }_{t}e=\mathrm{div}(u_t\nabla u)$, where $e=(u_t^2+|\nabla u{|}^2)/2$, over the shrinking ball $B_{R-t}(x_0)$. Differentiation of this moving-domain integral yieldsIntegrating in time proves the local wave energy estimate, without assumptions at spatial infinity.

Let the union of the initial supports be a compact set $K$. If $\mathrm{dist}(x_0,K)>T$, choose $r>0$ with $r+T<\mathrm{dist}(x_0,K)$. The initial energy on $B_{r+T}(x_0)$ vanishes. Applying the shrinking-ball identity up to every intermediate time shows $u_t$ and $\nabla u$ vanish throughout that cone. In particular, along the vertical segment through $x_0$, $u_t=0$; its initial value is also zero, so $u(T,x_0)=0$. This last value check removes the constant ambiguity invisible to gradient energy.

Consequently the finite propagation speed conclusion isThis set is compact for each finite $t$. The speed is at most one in these units, or $c$ for $u_{tt}-c^2\Delta u=0$. Time reversal gives the same conclusion for negative time.

For the wave equation with d'Alembert operator $\square =-{\partial }_t^2+\Delta$, the wave speed is one. The domain of dependence of $(t,x)$ is the initial ball $\overline{B}(x,t)$. More precisely, two solutions whose Cauchy data agree on a neighbourhood of that ball have the same value at $(t,x)$. Thus, if both Cauchy data have support in $K$, finite propagation speed givesThe Strong Huygens principle is sharper in three spatial dimensions: only the initial data on the sphere $|y-x|=t$ contribute, including the first spatial derivatives of the initial displacement there. The Kirchhoff formula makes this precise:In particular, if the Cauchy data have compact support in $\overline{B}(0,R)$, thenFor $t>R$, the solution therefore vanishes behind the inward edge $|x|=t-R$ as well as outside $|x|=t+R$. The absence of an interior tail is the additional content of the Strong Huygens principle.

We prove finite propagation speed using a shrinking cone energy argument. Fix $T>0$ and $x_{\ast }$, and suppose the Cauchy data vanish on $B(x_{\ast },T)$. For $0\le s<T$, define the local wave energyThe homogeneous wave equation gives the local conservation lawDifferentiate the integral over the moving ball. Its boundary moves inward with speed one, so the divergence theorem givesHere $n$ is the outward unit normal, ${\partial }_n$ the normal derivative, and ${\nabla }_{\tan }$ the component of the gradient tangent to the boundary. Since $E(0)=0$ and $e\ge 0$, we have $E(s)=0$. Hence both ${\varphi }_t$ and $\nabla \varphi$ vanish inside the backward light cone. Integrating ${\varphi }_t(s,x_{\ast })=0$ from the zero initial displacement gives $\varphi (s,x_{\ast })=0$; continuity then gives $\varphi (T,x_{\ast })=0$.

Applying the same energy estimate to the difference of two solutions proves the domain of dependence assertion. If $\mathrm{dist}(x_{\ast },K)>T$, the initial ball misses $K$, and the preceding argument proves the stated support bound. No disturbance propagates faster than one.

We establish a quantitative almost global existence for wave equations estimate. Take $0<\epsilon \le 1$, fix an integer $m\ge 6$, and use the commutation vector fields for the wave equation from the preceding part. Put $\partial =({\partial }_t,\nabla )$ and define the commuted wave energyAll these L2 norms are finite on any smooth existence interval by finite propagation speed. At $t=0$, the polynomial coefficients of the vector fields are bounded on the fixed compact support of the Cauchy data. Whenever a higher time derivative occurs, use ${\varphi }_{tt}=\Delta \varphi -{\varphi }_t^2$ and its differentiated versions to express it in terms of initial spatial derivatives. Every term contains at least one factor of $\epsilon$; consequentlyfor a constant depending only on finitely many derivatives and the support radius of ${\varphi }_0,{\varphi }_1$.

The commutators $[Z,{\partial }_{\mu }]$ are constant linear combinations of translations. Together with $[\square ,S]=2\square$ and the Leibniz rule, this shows that each commuted source is a finite linear combination of productsThis statement includes the extra copies of the original source produced by the scaling vector field. In each product put the factor with fewer commutations in the Lp norm $L^{\infty }$ and the other in the L2 norm. The lower order is at most $\lfloor m/2\rfloor$. Applying the Klainerman-Sobolev inequality to ${\partial }_{\mu }{Z}^J\varphi$ costs at most two additional commutations; commuting those past ${\partial }_{\mu }$ introduces only lower-order translations. Since $\lfloor m/2\rfloor +2\le m$,The inhomogeneous wave energy estimate now givesLet $T_{\epsilon }={\epsilon }^{-N}$. Use a bootstrap argument with $A_m(t)\le 2C_0\epsilon$ up to the smaller of $T_{\epsilon }$ and the maximal existence time. The energy estimate improves this toFor each fixed $N$,Choose ${\epsilon }_N\le 1$ so that $4C_m{C}_0\epsilon \log (1+{\epsilon }^{-N})\le 1/2$ for every $0<\epsilon \le {\epsilon }_N$. Then $A_m(t)\le (3/2)C_0\epsilon$, a strict improvement. A continuity argument closes the bootstrap argument.

Finally, the translation terms in $A_m$ control ordinary spatial Sobolev norms of $\partial \varphi$. The missing L2 norm of $\varphi$ satisfiesThus the full local-existence Sobolev norms remain bounded on this finite interval. The smooth continuation criterion for semilinear wave equations extends the solution past any finite endpoint before $T_{\epsilon }$. To see smooth persistence explicitly, the tame Sobolev product estimate gives $\Vert {\varphi }_t^2{\Vert }_{H^k}\le C_k\Vert {\varphi }_t{\Vert }_{\infty }\Vert {\varphi }_t{\Vert }_{H^k}$. Ordinary differentiated wave energy estimates therefore bound each higher derivative energy by its initial value times $\exp (C_k{\int }_0^t\Vert {\varphi }_t(s){\Vert }_{\infty }ds)$. This is finite on the interval already controlled by the base commuted wave energy; no separate ${\epsilon }_N$ is needed for each derivative order. ThereforeFor $\epsilon =0$ the zero solution is global. The same energy estimate in fact permits an exponential lower bound for the lifespan, which is stronger than any fixed inverse power.