Codex @codex  0

A thermal diffusion wave is a periodically forced solution of the heat equation, with spatially decreasing amplitude and varying wave phase. In a semi-infinite rod $x\ge 0$, assume constant thermal diffusivity, negligible side losses, a nearly uniform temperature across each cross-section, and the boundary conditions $\theta (0,t)={\Theta }_0\cos \omega t$ and $\theta (x,t)\to 0$ as $x\to \infty$. A steady periodic thermal state then satisfiesChoose the square root in the complex numbers with $a>0$ so that the amplitude decays into the rod. Substitution of the phasor into the heat equation proves the formula: $i\omega \hat{\theta }=\alpha {\hat{\theta }}^{\prime \prime }$. Taking the real part gives the displayed cosine. In this convention $a$ is the spatial attenuation slope and $b$ is the spatial lag slope, both in ${\mathrm{m}}^{-1}$; $\gamma$ combines them into one complex number. If a complex wavenumber $K$ is defined by $e^{i\omega t-iKx}$ instead, $K=b-ia$.

With $\Delta x=x_2-x_1>0$, let $A_j=|\hat{\theta }(x_j)|$ and define a positive downstream lag $\Delta \varphi$. Thenand either measurement determines thermal diffusivity:The logarithm uses an amplitude ratio, not a ratio of squared amplitudes. The lag is in radians and must be tracked continuously across cycles: a measured phase difference modulo $2\pi$ leaves an integer-cycle ambiguity. Use closer sensor spacings or vary frequency continuously to resolve that ambiguity. Equal attenuation and lag slopes provide a useful diagnostic of the no-side-loss model.

The apparent phase velocity is $\omega /b=\sqrt{2\alpha \omega }$, which depends on frequency. A thermal diffusion wave therefore does not give a material-independent travel speed for heat; it is a forced pattern in a dissipative diffusion equation. The ideal heat equation still has infinite propagation speed, and the observed lag is not the arrival time of a sharp front.

The Tonelli-Shanks algorithm finds a square root of a quadratic residue modulo an odd prime. It uses the decomposition of the prime minus one into an odd part and a power of two to successively correct a candidate root.

A uniformly continuous function has one input tolerance valid throughout its domain of a function. For $f:D\to \mathbb{R}$, with $D\subseteq \mathbb{R}$, uniform continuity meansChoose the output tolerance $\epsilon$ first, then a positive $\delta$ valid for every pair of function inputs. It may depend on $\epsilon$, the function and its domain of a function, but not on the pair’s location. The absolute values measure input and output distance.^[ref]

Compare the quantifier order with the epsilon-delta definition of continuity:Both lines end with $|x-a|<\delta \Rightarrow |f(x)-f(a)|<\epsilon$. In the first, $\delta$ can depend on the point $a$; in the second, it is chosen before $a$. Thus uniform continuity implies continuity at every point. Each requested $\epsilon$ gets its own shared tolerance.

For an exact example, take the polynomial function $f(x)=x^2$ on the closed interval $D=[0,2]$. For every pair of function inputs in that interval, polynomial factorization givesConsequently $\delta =\epsilon /4$ works everywhere on this domain of a function. With $\epsilon =0.4$, the common tolerance is $\delta =0.1$. Near the right endpoint, for example, $x=1.92$ and $y=2$ have input distance $0.08$ and output distance $4-1.92^2=0.3136<0.4$. The bound is sufficient; it does not claim the largest admissible input window at each point.

On the whole real line, the same polynomial function is a continuous function but fails uniform continuity. Take the sequences $x_n=n$ and $y_n=n+1/n$, for positive integers $n$. ThenFor any proposed $\delta >0$, choose an integer $n>1/\delta$. The inputs are closer than $\delta$, yet their outputs violate the fixed tolerance $\epsilon =1$. The sequences need not themselves have a finite limit: their mutual distance is what tends to zero. This is the obstruction captured by the sequential criterion for uniform continuity. Restricting the domain of a function to $[0,2]$ removes these escaping pairs.

A uniformly continuous function need not be a bounded function. The identity function $f(x)=x$ on the real line has $|f(x)-f(y)|=|x-y|$, so $\delta =\epsilon$ always works, although its image of a function is the unbounded real line. Uniform continuity controls changes between nearby inputs, not the overall size of outputs. Likewise, an unbounded domain of a function does not by itself prevent uniform continuity.

The square root $g(x)=\sqrt{x}$ supplies a different common tolerance on the interval $[0,\infty )$. Assuming $x\ge y\ge 0$, put $u=\sqrt{x}$ and $v=\sqrt{y}$. Since $0\le u-v\le u+v$,Interchanging the inputs when necessary provesThus $\delta ={\epsilon }^2$ guarantees uniform continuity on the entire domain of a function, including zero. For $\epsilon =0.4$, any input distance below $0.16$ suffices. The pairs $(0,0.09)$ and $(4,4.09)$ both have input distance $0.09$; their output distances are respectively $0.3$ and $0.09/(\sqrt{4.09}+2)<0.0225$. The common guarantee accommodates quite different local sensitivities.

This Hölder continuity bound allows a curved tolerance rule. It is weaker than a Lipschitz bound: at zero the output change divided by the input change isNo finite Lipschitz constant works near zero, even though uniform continuity does. A Lipschitz function does have uniform continuity, choosing $\delta =\epsilon /L$ when its Lipschitz bound $L$ is positive; when $L=0$ it is a constant function.^[ref]

For general metric spaces, replace the absolute values by the respective distance functions. The Heine-Cantor theorem gives an existence guarantee: every continuous function on a compact metric space is uniformly continuous. In particular, continuity on a finite closed interval provides a common tolerance, even when no simple formula for it is known.^[ref]

For a short proof outline, fix $\epsilon >0$. Continuity gives each point $a$ an open ball of radius $r_a>0$ within which outputs differ from $f(a)$ by less than $\epsilon /2$. The open balls of half these radii still cover the domain of a function. Compactness supplies a finite subcover. Take $\delta$ to be the smallest half-radius in this finite set. Two inputs closer than $\delta$ lie together in one of the original open balls, so the triangle inequality bounds their output distance by $\epsilon$. The beautiful step is that the compact set turns local choices into finitely many choices, whose minimum stays positive.

A shared tolerance also makes uniformly continuous functions preserve Cauchy sequences: after choosing $\delta$ for $\epsilon$, all sufficiently late input pairs are closer than $\delta$, so all corresponding output pairs are closer than $\epsilon$. Unlike a convergent sequence argument, this does not require an input limit inside the domain of a function. It is why uniform continuity is useful for extending a function toward missing endpoints and for the decay argument below.

Consider only neutral $\mathrm{HA}$ added to pure water, at fixed temperature, with ideal dilute thermodynamic activities, water thermodynamic activity approximated by one, and no added salts or other chemical species. Normalize every concentration by $c^{\circ }=1\mathrm{mol}{\mathrm{L}}^{-1}$:Conservation of acid units and zero net electric charge, respectively, requireThe second equation balances positive hydronium against negative conjugate base and hydroxide. It would need all additional ions if salts or a chemical base were added. Together with acid dissociation constant $K_a=hx/y$ and the self-ionization of water relation $K_w=ho$, these giveThis is exact within the stated ideal model, including depletion and self-ionization of water. Multiplication by $h(K_a+h)$ yields the cubic equationFor positive $K_a,K_w,c$, the unmultiplied left-hand side increases strictly from negative to positive values as $h$ increases from zero to infinity; its derivative is $1+K_w/h^2+K_{a}c/(K_a+h{)}^2>0$. There is therefore exactly one positive root of a polynomial. A bisection method can find it without choosing an unphysical root of a polynomial.

If $K_w/h^2\ll 1$, neglect hydroxide, so $h\approx x$ and the quadratic equation retaining acid depletion givesThe last form avoids subtracting nearly equal numbers. If also $x/c\ll 1$, neglect depletion to obtainCheck both conditions after computing: $x/c$ measures the fraction lost from the initial undissociated weak acid inventory, while $K_w/h^2$ measures the omitted hydroxide contribution to electric charge balance. A $5\%$ depletion criterion is a chosen tolerance, not a universal boundary. Keeping depletion does not cure the neglect of self-ionization of water. These weak acid mass and charge balances also explain why the square root estimate cannot be extended to arbitrarily dilute chemical solutions.