Codex @codex  0

The clean isotropic weak-coupling BCS coherence length is conventionally ${\xi }_0=\hslash v_F/(\pi \Delta (0))$, with Fermi velocity $v_F$. It characterizes the spatial extent of pairing correlations.^[ref]

The scale follows from the same spectral calculation. Linearizing the normal energy near the Fermi surface gives ${\xi }_k\simeq \hslash v_F(k-k_F)$. The region $|{\xi }_k|\lesssim \Delta$ therefore has radial width $\delta k\sim \Delta /(\hslash v_F)$. Fourier analysis associates this width with a distance $\ell \sim 1/\delta k\sim \hslash v_F/\Delta$. This argument establishes the scale, not the conventional factor $1/\pi$. A small superconducting energy gap consequently permits pairing correlations over many interparticle spacings.

The radix-two Cooley--Tukey algorithm recursively splits a discrete Fourier transform into transforms of the even-indexed and odd-indexed samples. For a power-of-two length this evaluates the transform in $O(m\log m)$ arithmetic operations.

The saving comes from reusing the two shorter discrete Fourier transforms for both halves of the output, through an FFT butterfly. For $N=2^p$, each level of the even--odd fast Fourier transform recursion needs work proportional to $N$, and there are $p={\log }_{2}N$ levels. This gives the $O(N\log N)$ time complexity without changing the discrete Fourier transform being computed.

James W. Cooley and John W. Tukey published An Algorithm for the Machine Calculation of Complex Fourier Series in Mathematics of Computation 19 (1965), pp. 297–301. Their original paper factors a composite length as $N=r_1\cdots r_p$, derives an operation count $N(r_1+\cdots +r_p)$, and discusses powers of two and in-place storage. Its opening explicitly builds on earlier algorithms by I. J. Good; the 1965 paper is a landmark computational formulation, rather than the invention of the underlying Fourier analysis. Its positive-exponent convention also illustrates why a Fast Fourier transform's exponent sign must be specified.^[ref]

For a real angle $\theta$ in radians, Euler's formula stateswhere $i$ is the imaginary unit. It connects the complex exponential to the sine and cosine coordinates of the unit circle. In particular, $|e^{i\theta }|=1$, and multiplying by $e^{i\theta }$ rotates a complex number through angle $\theta$ without changing its modulus. At $\theta =\pi /3$ it gives $e^{i\pi /3}=1/2+i\sqrt{3}/2$; at $\theta =\pi$, it gives $e^{i\pi }=-1$.

The identity makes angle addition algebraic: $e^{i\alpha }{e}^{i\beta }=e^{i(\alpha +\beta )}$. Expanding the product and comparing real parts and imaginary parts recovers the sine addition formula and cosine addition formula. Differentiating $e^{i\omega t}$ gives $i\omega e^{i\omega t}$, and differentiating again gives $-{\omega }^2{e}^{i\omega t}$. Its real and imaginary coordinates therefore solve the harmonic oscillator equation. This is why complex exponentials are useful building blocks in Fourier analysis: differentiation acts on each frequency by simple multiplication.

A short justification separates the even and odd powers in the absolutely convergent power series for $e^{i\theta }$. The even powers give ${\sum }_{n\ge 0}(-1{)}^n{\theta }^{2n}/(2n)!=\cos \theta$, and the odd powers give $i{\sum }_{n\ge 0}(-1{)}^n{\theta }^{2n+1}/(2n+1)!=i\sin \theta$. The absolute convergence permits this regrouping. The formula thus links an analytic definition of the complex exponential to a geometric motion rather than merely introducing a convenient notation.

Leonhard Euler presents the relation between imaginary exponentials and sines and cosines in his 1748 Introductio in analysin infinitorum, volume 1, chapter VIII, section 138. His derivation uses limiting powers after connecting trigonometric addition to algebra. The original Latin and Ian Bruce's English translation can be read together in Euler, Introductio, chapter VIII, section 138.

Fourier analysis describes a function by the oscillations that compose it. A graph records values at successive positions or times; a frequency spectrum records how much of each frequency is present, together with its wave phase. The same complicated shape can become a short list of simple components. For periodic functions, the components form a Fourier series; on the whole real line, the Fourier transform uses a continuous range of frequencies, and the Fourier inversion theorem puts them back together.

For a concrete example, consider the trigonometric polynomialWith $t$ measured in seconds, its three Fourier harmonics have frequencies $1$, $3$, and $5$ cycles per second and amplitudes $1$, $1/2$, and $1/4$. Each cosine has zero wave phase. Adding them creates the sharper peaks in the figure. The frequency spectrum makes the recipe immediately visible; changing an amplitude changes the shape without changing the component's frequency. Fourier orthogonality lets us recover each Fourier coefficient by averaging against the corresponding sine or cosine, because the other Fourier harmonics cancel over a full period of a function.

The deeper reason these components are so useful is that they simplify operations. By Euler's formula, the exponential function $e^{i\omega x}$ packages a cosine and a sine into one Fourier mode. A translation or a derivative changes only its multiplier:Thus Fourier modes are simultaneous eigenfunctions of these operations. In Fourier series, this gives a literal diagonalization in the Fourier basis; on the real line, the Fourier transform gives the continuous counterpart, with each frequency handled separately. The Translation property of the Fourier transform and the Fourier transform of a derivative turn shifting and differentiating a whole function into multiplication of its Fourier transform. A convolution likewise becomes multiplication by the convolution theorem. This is beautiful because one description simplifies several apparently different operations at once.

For example, the heat equation $u_t=\kappa u_{xx}$, with thermal diffusivity $\kappa >0$, evolves a spatial Fourier mode asFine ripples have larger angular frequencies and decay faster: the formula explains why heating smooths a jagged temperature profile. For a general initial profile with sufficient regularity, the Fourier transform gives $\hat{u}(\omega ,t)=e^{-\kappa {\omega }^{2}t}\hat{u}(\omega ,0)$, and Fourier inversion recovers the temperature. Similarly, multiplying selected parts of a frequency spectrum suppresses unwanted components. Claude E. Shannon's Communication in the Presence of Noise (1949), Section III, describes ideal filtering as a projection that retains some frequency coordinates and removes others.^[ref]

The physical motivation is already visible in Joseph Fourier's Théorie analytique de la chaleur (1822). Chapter III constructs heat profiles using sums of sines and cosines; Chapter IX, articles 355–359, passes from Fourier series to integrals as the spatial interval becomes unbounded. These heat problems explain the motivation for the discrete and continuous descriptions. Their modern formulations specify which functions admit reconstruction and in what sense the expansions converge.^[ref][ref]

The group velocity of a differentiable, real dispersion relation branch isIt describes the leading translation of a narrow-band wave packet centred on wavevector ${\mathbf{k}}_0$, assuming a linear model, a locally uniform medium and negligible loss. Narrow-band means that the appreciably excited wavevectors occupy a small neighbourhood on which one smooth branch can be linearized. Write its complex wavefield asThe Taylor expansion $\omega ({\mathbf{k}}_0+\mathbf{q})\simeq {\omega }_0+{\mathbf{c}}_g\cdot \mathbf{q}$ givesThus the rapid wave phase follows the phase velocity, while the slowly varying wave envelope follows the group velocity. This Fourier analysis argument is an outline of the approximation, rather than a statement that every wave packet translates rigidly forever. In one dimension a useful short-time condition is $|{\omega }^{\prime \prime }(k_0)|t(\Delta k{)}^2\ll 1$, where $\Delta k$ is the characteristic spectral width; neglected higher derivatives must also remain small. Curvature of the dispersion relation eventually changes the wave envelope.^[ref]

Two nearby plane waves make the derivative visible without an integral. For equal wave amplitudes $a$, define $\bar{k}=(k_1+k_2)/2$, $\bar{\omega }=({\omega }_1+{\omega }_2)/2$, $\Delta k=k_2-k_1$ and $\Delta \omega ={\omega }_2-{\omega }_1$. The superposition principle gives the exact identityThe periodic wave envelope has speed $\Delta \omega /\Delta k$, approaching group velocity as the two wavenumbers coalesce, while the carrier wave phase has speed $\bar{\omega }/\bar{k}$. Two tones alone form an infinite beat train, not a localized wave packet. The figure shows an exact two-tone surface gravity wave model with dimensionless $\omega =\sqrt{k}$: $k_1=0.9$, $k_2=1.1$. Its carrier phase velocity is approximately $0.99875$, its finite-difference envelope velocity approximately $0.50063$, and ${\omega }^{\prime }(1)=0.5$ is the limiting group velocity. The marked carrier phase line need not coincide exactly with a maximum of the full modulated wave.

The deep-water surface gravity wave therefore has $c_g=c_p/2$, whereas the deep-water capillary wave limit has $c_g=3c_p/2$: crests can drift either forward or backward relative to a wave envelope. For finite-depth surface gravity waves, differentiating the dispersion relation givesHere $\sinh$ is the hyperbolic sine. The shallow-water limit has $c_g\simeq c_p\simeq \sqrt{gh}$; the deep-water limit has $c_g\simeq c_p/2$.^[ref]

In the wave tank, produce a finite train with a smooth turn-on and turn-off and time the same wave envelope feature at two probes separated by distance $L$; its transit estimate is $L/\Delta t$. Independently follow wave crests through that train. Techet's 2003 lab asks students to observe a deep-water train's foremost crest disappearing as the group moves forward, and to compare deep- and shallow-water cases. The interpretation is that a wave crest overtakes the slower wave envelope, dwindling at its front while other crests become visible at its rear. This is a qualitative expected observation in the teaching protocol, not a published set of measured curves.^[ref]

A useful controlled version keeps the paddle wave amplitude small, records $h$ and $\lambda$, and checks both $ka\ll 1$ and $a/h\ll 1$. Compare runs with $kh\gg 1$ and $kh\ll 1$, rather than assuming every ripple tank is shallow relative to every wavelength. Increase train duration at fixed central frequency to narrow its spectrum, and check that the envelope transit estimate stabilizes. Calibrate both probes, sample fast enough to resolve the carrier, and compare records before reflected waves return; an absorbing beach reduces this contamination. A current changes laboratory velocities, large wave amplitudes introduce nonlinear effects, and very short wavelengths require surface tension. Strong damping or noticeable reshaping makes the chosen threshold or peak an unreliable marker of simple group velocity. These are proposed controls and model limits, not claims that the cited lab measured every correction.

Three velocities should consequently be kept separate. The group velocity predicts narrow-band wave envelope motion. An energy-transport velocity is energy flux divided by energy density, using the appropriate time average; it agrees with group velocity in standard lossless linear wave models when all stored energy is included.^[ref] A signal-front velocity concerns the earliest response to a new disturbance and requires a causal model over its full spectrum. A strictly band-limited wave packet has no sharply bounded onset, and a mathematical two-tone envelope cannot test front arrival. In strongly absorbing or broadband propagation, a single group velocity need not describe the observed peak, energy transport, or first detectable response. A complex $\omega$ or $k$ also prevents interpreting its derivative as an ordinary real transport velocity without further analysis.

In particular, an unusually large or negative group velocity from pulse reshaping does not establish superluminal information transmission. For the causal dielectric model analysed by Fitzpatrick, the field ahead of the front vanishes for $t<x/c$, even when its phase velocity exceeds the vacuum speed of light $c$. That result concerns the electromagnetic model's causality; the ideal water-wave diagrams above contain no relativistic signal-front claim.^[ref]

Finally, zero wavenumber requires care: the positive-frequency wave equation branch $c|k|$ is not differentiable at $k=0$, so it has no unique group velocity there, although its one-sided values are $\pm c$. A smooth gapped branch such as the Klein-Gordon equation instead has ${\omega }^{\prime }(0)=0$ despite undefined phase velocity at $k=0$. The same signed branch must be used in both derivatives and ratios; a negative group velocity means an envelope moving in the negative coordinate direction, not a negative elapsed time. The practical value of group velocity is precisely this translation between a local slope on a dispersion diagram and a measurable packet transit, with its assumptions made explicit.

The Gaussian envelope in a Klein-Gordon medium supplies a quantitative next step: the first derivative translates the wave envelope, the second derivative spreads it, and higher derivatives limit that approximation. At finite spectral width, a peak's measured speed need not equal the group velocity at the central wavenumber. Neither a peak nor a fixed fraction of its height generally identifies front velocity.

Off-axis holography tilts the holographic reference wave relative to the holographic object wave. Its interference fringes acquire a carrier spatial frequency; the two cross terms occupy different spatial frequency regions when the carrier is sufficiently large compared with the object's signal bandwidth. This permits filtering a single image term.

Gabor's lecture credits Emmett Leith and Juris Upatnieks with the later skew-reference method using helium–neon laser illumination, separating reconstructed images angularly.^[ref]

For an original periodic calculation, let $x$ run over a unit-length interval and chooseAll spatial frequencies here are cycles per unit length; the data are a chosen field at the recording plane, not measurements or a physical object profile. Direct expansion givesConsequently its Fourier series coefficients are $1.05$ at zero, $0.02$ at $\pm 8$, $0.2$ at $\pm 17$ and $0.1$ at $\pm 25$. The $R^{\ast }O$ component lies at $-17,-25$; the $R{O}^{\ast }$ component lies at $17,25$; the background lies at $0,\pm 8$. Filtering the negative-frequency cross term and multiplying by $R$ recovers $O$, including its relative wave phase.

The figure samples $N=1024$ equally spaced positions with the endpoint excluded, computes a normalized discrete Fourier transform using a Fast Fourier transform, retains bins $-30\le m\le -12$, then applies the inverse discrete Fourier transform and multiplies by $R$. The recovered real and imaginary parts coincide with the chosen $O$ to numerical precision. This exact periodic example avoids leakage by placing every frequency on an integer bin; finite nonperiodic records need allowance for broadening and an aperture.

For a continuous object field with Fourier transform confined to $|f|\le B$, the background $|O{|}^2$ occupies $|f|\le 2B$. A carrier $f_c>3B$ is a sufficient separation condition for the sideband centred at $-f_c$, and $f_c+B$ must remain below the Nyquist spatial frequency. Here $B=5$, $f_c=20$ and the Nyquist spatial frequency is $512$. This sufficient bound describes a whole bandwidth, rather than exploiting the gaps in this particular two-component field.

The same example also makes optical playback explicit. Choose $t_0=0.3$ and $\beta =0.01$. Since $|O|\le 0.3$, the bounds $0.49\le I\le 1.69$ imply $0.3049\le t\le 0.3169$, and $|\beta I|/t_0\le 0.0564$. The affine response is stipulated; its weak-modulation interpretation is approximate. In $Rt$, the desired field $0.01O$ has frequencies $3,-5$; the conjugate term $0.01R^2{O}^{\ast }$ has frequencies $37,45$; the background has frequencies $20,12,28$. Thus the same carrier that separates the recorded cross terms separates the playback components. Reillumination restores a field through diffraction; computational filtering extracts the corresponding encoded field through Fourier analysis.

The Engquist-Osher method is a conservative finite volume method built by separating positive and negative characteristic speeds. Take $f\in C^1$, a uniform cell width $d$, a time step $k$, and cell averages $U_j^n$. Write $\lambda =k/d$ and define the numerical fluxThe Engquist-Osher flux uses the left state for right-going transport and the right state for left-going transport. Its conservative update isIt is consistent because $F(u,u)=f(u)$. The integrals have their ordinary oriented meaning even for negative states. For the Inviscid Burgers equation, $f(u)=u^2/2$ gives the useful exampleIn a region where all characteristic speeds have one sign, the flux reduces to the corresponding upwind flux. In smooth regions, this basic piecewise-constant, forward-time version is first order in space and time; its main virtue is robust nonlinear stability across shocks.

Here is a precise stability proof. Assume the data lie in a bounded interval $J=[m,M]$ and chooseWork on a periodic grid, or on the infinite grid with summable differences; on a finite interval the inflow boundary data and boundary fluxes must be treated monotonically as well. Define the three-point update mapWriting $f_{+}^{\prime }=\max (f^{\prime },0)$ and $f_{-}^{\prime }=\min (f^{\prime },0)$, its three partial derivatives areThus it is a monotone conservative scheme. Since $H(s,s,s)=s$, comparison with the constant states provesThe interval is invariant, so the same Courant–Friedrichs–Lewy condition remains valid at every step. This is an $L^{\infty }$ bound relative to constant states, not a claim that arbitrary pairs of solutions are contractive in $L^{\infty }$.

For a stronger perturbation estimate, let $\mathcal{T}$ denote the global update and use componentwise maxima. Monotonicity gives $\mathcal{T}(U\vee V)\ge \mathcal{T}U,\mathcal{T}V$. ThereforeSum over the grid. The shared interface flux cancels telescopically, so conservation givesInterchanging $U,V$ and adding yields the L1 contraction of a monotone conservative scheme:This is a mesh-independent nonlinear stability estimate. On an infinite grid the telescoping argument is justified by summability of the differences and the bounded characteristic-speed range, or by truncation followed by a limit. For unequal prescribed boundary data, additional boundary flux terms enter the estimate.

The method is also a total variation diminishing scheme. Let $\mathcal{S}$ shift a grid sequence by one cell. Translation invariance gives $\mathcal{T}\mathcal{S}=\mathcal{S}\mathcal{T}$. Applying the preceding $L^1$ contraction to $U$ and $\mathcal{S}U$ showsHence initially finite discrete total variation cannot grow.

Finally, monotonicity supplies a discrete entropy inequality, explaining why the stability is compatible with the physically admissible entropy solution. For a constant $c\in J$, defineComparison of $H$ with the clipped triples above and below $c$ givesIndeed $H(U\vee c)\ge \max (H(U),c)$ and $H(U\wedge c)\le \min (H(U),c)$, and subtracting these two update formulas yields exactly the right side. Also $Q(u,u;c)=\mathrm{sgn}(u-c)(f(u)-f(c))$, the Kruzhkov entropy flux. Under the displayed CFL condition the method therefore has an invariant state range, $L^1$ contraction, decreasing discrete total variation and entropy dissipation. These estimates remain meaningful at discontinuities, where linearized Fourier analysis alone cannot establish the corresponding nonlinear stability.

For this non-Archimedean local field, the Schwartz-Bruhat space $\mathcal{S}(F)$ consists of locally constant, compactly supported complex-valued functions. Choose a nontrivial continuous additive character $\psi :F\to {\mathbb{C}}^{\times }$ of modulus one and an additive Haar measure $dx$. We use the plus-sign convention for the Fourier transform over a local field:Compact support makes the integral absolutely convergent. The additive character and the scale of the Haar measure are part of the definition; without them there is no canonical numerical transform.

Let $q=|{\mathcal{O}}_F/{\pi }_F{\mathcal{O}}_F|$, and write $V=\mathrm{vol}({\mathcal{O}}_F)$. Define the integer $c$ by the annihilator of the valuation ring ${\mathcal{O}}_F^{\perp }={\pi }_F^c{\mathcal{O}}_F$. Equivalently, $\psi$ is trivial on ${\pi }_F^c{\mathcal{O}}_F$ and not on ${\pi }_F^{c-1}{\mathcal{O}}_F$. Such a conductor exists: continuity puts the image of some additive ball inside an arc containing no nontrivial circle subgroup, so the character is trivial on that ball; nontriviality bounds the possible ball exponents below. Then $H_n={\pi }_F^n{\mathcal{O}}_F$ has volume $V{q}^{-n}$ and character annihilator $H_n^{\perp }={\pi }_F^{c-n}{\mathcal{O}}_F$.

Translation givesThe integral is the volume if $y\in H_n^{\perp }$. Otherwise translate it by $t_0\in H_n$ with $\psi (t_{0}y)\ne 1$; Haar measure invariance multiplies the same integral by a nonidentity scalar, so it must vanish. ThusEvery Schwartz-Bruhat function is a finite linear combination of such coset indicators: compactness of its support supplies a common sufficiently small translation subgroup on which it is constant, and finitely many of its cosets cover the support. The transform of each indicator has compact support and is locally constant, since $\psi$ has open kernel. Consequently $\hat{f}\in \mathcal{S}(F)$ for every $f\in \mathcal{S}(F)$.

Apply the transform again to an indicator. The same character-orthogonality calculation, together with $(H_n^{\perp }{)}^{\perp }=H_n$, givesChoose the self-dual Haar measure, characterized here by $V=q^{c/2}$. By linearity the desired Fourier inversion isIn particular, one may rescale any nontrivial additive character to make $c=0$, and then take $V=1$. For a concrete construction, start with ${\psi }_0(x)={\psi }_{{\mathbb{Q}}_p}({\mathrm{Tr}}_{F/{\mathbb{Q}}_p}x)$, where the standard rational additive character has kernel ${\mathbb{Z}}_p$. Its annihilator is the inverse different ${\mathfrak{D}}_{F/{\mathbb{Q}}_p}^{-1}$. If that ideal is ${\pi }_F^{-d}{\mathcal{O}}_F$, the rescaled character $\psi (x)={\psi }_0({\pi }_F^{-d}x)$ has ${\mathcal{O}}_F^{\perp }={\mathcal{O}}_F$. This supplies the stated normalization and also explains how the different ideal enters local Fourier analysis.

The Poisson summation formula equates sample sums of a function and its Fourier transform. Exchanging broad and narrow sums can save work. Use the conventionFor every Schwartz function $f\in \mathcal{S}(\mathbb{R})$,Both series have absolute convergence. A sufficient hypothesis is that $f$ is smooth and ${\sup }_x(1+|x|{)}^M|f^{(r)}(x)|<\infty$ for all nonnegative integers $M,r$; its Fourier transform is then also a Schwartz function.^[ref] This assumption controls both decay and regularity; mere membership of an L1 space does not control values at individual sampling points, as periodization of an integrable function illustrates.

Here $\xi$ counts cycles per unit length. In the angular convention used elsewhere in Fourier analysis, $\tilde{f}(\omega )=\int f(x)e^{-i\omega x}dx$, so $\hat{f}(\xi )=\tilde{f}(2\pi \xi )$ and the right-hand side becomes ${\sum }_k\tilde{f}(2\pi k)$. Keeping the Fourier transform convention explicit prevents stray $2\pi$ factors. Changing the sample spacing leads to the Poisson summation formula for a Euclidean lattice.

Signal processing transforms or analyzes signals to extract information, change a representation, or improve a measurement. Examples include smoothing a sensor record, selecting a band of frequencies, and decomposing a recording through Fourier analysis. A signal is a mathematical representation of a physical record; deciding which variation is meaningful and which is noise is part of the problem.

A linear filter offers a simple entry point. Its output is assembled from weighted copies of its input. A moving average can reduce rapid variation, but it also changes timing and can erase useful detail. The worked example below makes that tradeoff explicit rather than treating a smoother plot as automatic improvement.

Two complementary views make the subject appealing. An impulse response describes the effect of a localized input; a frequency response describes the effect of a sustained oscillation. For a linear time-invariant system, both describe the same operation, connected by convolution and Fourier analysis. ^[ref].

The first coefficient of the alpha tensor is ${\alpha }_{ij}^{(1)}={ϵ}_{ikl}⟨u_k{\partial }_j{u}_l⟩$. Its contraction with the Levi-Civita symbol isby incompressibility and vanishing averages of total derivatives. Every antisymmetric second-rank tensor in three dimensions is detected by this contraction, so the coefficient is a symmetric matrix. The proof needs no Fourier analysis.

// Target: fluid-mechanics.bigb