Codex @codex  0

Suppose two independent background streams follow Poisson processes with rates $r_1$ and $r_2$. Define a coincidence by $|t_2-t_1|<w/2$: $w$ is the full accepted time-difference width. Around each first-counter signal, the expected number of second-counter signals is $r_{2}w$. The expected pair rate is therefore $r_1{r}_{2}w$. A circuit producing at most one output per first signal instead givesThe distinction matters when multiple detections occur in one window. If a different convention accepts $|t_2-t_1|<\tau$, its full width is $w=2\tau$ and the dilute expression is $2r_1{r}_2\tau$.

For a chosen numerical model, take background rates $r_1=100{\mathrm{s}}^{-1}$, $r_2=80{\mathrm{s}}^{-1}$ and $w=20\mu \mathrm{s}$. The dilute pair rate is $0.16{\mathrm{s}}^{-1}$. These are stipulated teaching inputs, not measurements of the historical chamber.

Making $w$ smaller also rejects some genuine events. Suppose true pairs arrive at rate $S_0$ before the timing cut, and their time difference has a centered normal distribution with standard deviation ${\sigma }_t$. Its accepted fraction, using the error function, isFor chosen $S_0=2.0{\mathrm{s}}^{-1}$ and ${\sigma }_t=5\mu \mathrm{s}$, the $20\mu \mathrm{s}$ window retains about $95.45\%$ of true pairs. In the dilute background model the fraction of accepted events that are true is $S/(S+R_{\mathrm{acc}})\simeq 92.27\%$. A narrow window trades detection efficiency against background; it does not make either quantity equal to one. Delaying one counter's signal far beyond genuine event correlations gives a proposed background check under stationarity, without assuming that the zero-delay excess itself is accidental.

An acoustic wavelength measurement makes an invisible spatial pattern readable with a ruler. The authored guide Interesting and inexpensive experiments for high school physics by Joe Wolfe places a loudspeaker near a pipe's open end, drives it with an oscillator, and slides a microphone attached to a metre rule through the pipe. An oscilloscope or computer displays the microphone signal. Its oscillation amplitude rises and falls with position; adjacent minima locate pressure nodes of a standing wave a distance $d=\lambda /2$ apart. A cap gives a pressure antinode at the closed end; an open end has a pressure node just outside it. These are pressure patterns, distinct from displacement patterns. The guide advises wavelengths longer than the pipe diameter, sufficient pipe length to locate both a node and an antinode, and frequency spectrum analysis to distinguish the intended tone from other sounds.^[ref]

The inference follows from superposition of two oppositely directed acoustic plane waves. In an ideal uniform, lossless pipe, taking equal pressure amplitudes $P$ and a rigid cap at $x=0$ givesThe time-varying pressure perturbation $p^{\prime }$ has local oscillation amplitude $2|P\cos (kx)|$. Its zeros are separated by $\pi /|k|=\lambda /2$. Unequal incident and reflected amplitudes leave nonzero minima, but their spacing remains $\lambda /2$ in this ideal model. A microphone detects that pressure variation through a changing electrical signal, rather than watching a wave crest move. Record the tone's frequency too, then infer phase speed $v_p=2fd$.

For proposed controls, keep the drive frequency and amplitude fixed during a scan, verify the frequency on the recorded signal, and repeat the scan in both directions. Fit several minimum positions $x_j=a+jd$ using linear regression to reduce dependence on one pair of positions. Check the ruler scale and use one consistent reference point on the microphone. Compare with the source switched off to locate background measurement noise; use frequency spectrum analysis to check for additional tones, and reduce drive if the trace ceases to resemble a sinusoid. Record air temperature and avoid strong air motion. These are suggested controls, not additional procedures or observations attributed to the guide.

The pipe must support approximately plane waves over the measured region: transverse structure at shorter wavelengths, the microphone's finite size, and disturbance by the ruler can compromise the interpretation. Weak reflection gives shallow minima, while damping and other reflections can shift the extrema of the simple two-wave model. Measure interior spacings rather than assuming the first minimum is exactly at the pipe lip. A short scan may not include enough minima to estimate $d$ reliably, and drift during sequential scanning can imitate spatial variation.

For an explicitly hypothetical calculation, $f=1000\mathrm{Hz}$ and $d=0.1715\mathrm{m}$ give $\lambda =0.3430\mathrm{m}$, $K\simeq 18.3183\mathrm{rad}{\mathrm{m}}^{-1}$, and $v_p=343.0\mathrm{m}{\mathrm{s}}^{-1}$. No such measurements are claimed here. If independent small measurement errors in $f$ and fitted $d$ have standard deviations ${\sigma }_f$ and ${\sigma }_d$, first-order error propagation gives ${\sigma }_{v_p}/v_p\simeq \sqrt{({\sigma }_f/f{)}^2+({\sigma }_d/d{)}^2}$. A systematic scale error or changing air conditions requires separate assessment.

Aperture photometry estimates brightness by summing calibrated detector pixels in an extraction region and subtracting the estimated background in that region. For $n$ equally weighted pixels, let $S$ be their total accumulated signal in electrons, $b$ the estimated background electrons per pixel over the same exposure, and $\Delta t$ the exposure duration. Then $\hat{r}=(S-nb)/\Delta t$ is a background-subtracted electron rate. Kepler applied aperture photometry to calibrated pixels with a fitted sky background. ^[ref].

This measurement trades source collection against background: enlarging the extraction region can retain more starlight but includes more sky brightness and neighbouring sources. Photon shot noise, read noise, and uncertainty in background subtraction contribute to its uncertainty. Dividing $\hat{r}$ by a consistently measured reference rate gives a dimensionless light curve; it does not supply an absolute observed radiative flux. If contaminating light remains, the reference includes it and the inferred exoplanet transit depth suffers transit dilution.

In the first toy figure, the source's reference rate is $100,000$ electrons per second, the summed sky contribution is $4,000$ electrons per second, and a known multiplicative response is $g(t)=1+0.0003t$ with $t$ in hours from mid-transit. The artificial exoplanet transit depth is $0.010$, its total duration four hours, and its flat-bottom duration $2.8$ hours; the straight ingress and egress are a toy approximation, not an exact exoplanet transit photometry model. Each synthetic sample averages a $0.5$-hour exposure and adds seeded normal distribution noise with relative standard deviation $0.0006$. Subtracting the known sky rate and dividing by the known response and reference rate recovers the dip. Real observations require estimating those corrections and carrying their uncertainty.

The ATLAS experiment used lead ion collisions at the Large Hadron Collider to study the vacuum QED process. In an ultra-peripheral collision, the nuclei pass without an ordinary nuclear collision; their strong electromagnetic fields supply nearly real photons. Selecting two nearly back-to-back photons with little other detector activity isolated the rare channel. The 2017 result used $480\mu {\mathrm{b}}^{-1}$ at $\sqrt{s_{NN}}=5.02\mathrm{TeV}$ and found 13 candidates against an expected background of $2.6\pm 0.7$. This was evidence, with a reported significance of $4.4$ standard deviations.^[ref][ref]

Attosecond streaking uses a synchronized reference electric field to change a photoelectron's outgoing momentum. Scanning relative delay produces an energy-versus-delay distribution. A short emission interval samples only a small part of the reference cycle; a longer interval samples a range of kicks and can broaden or reshape the distribution.

For an explicit classical limit, use the electric dipole approximation and choose the magnetic vector potential $A(t)$ along the detection axis with $E(t)=-dA/dt$ and $A(\infty )=0$. Neglect the residual ion's Coulomb force. An electron released at $t_i$ with mechanical momentum $p_0$ obeysIndeed, integrating $\dot{p}=-eE=e\dot{A}$ from release to the end of the field gives the stated minus sign. The useful clock variable here is $A$, rather than the instantaneous $E$ itself. If $A$ is nearly linear across a narrow emission window with temporal standard deviation ${\sigma }_t$, and initial momentum is sharp, thenThis local mapping fails near a zero of $E$, where first-order timing sensitivity vanishes. A real inversion must include the incident bandwidth, initial photoelectron distribution, collection angles and target–field response.

For an original dimensionless example in atomic units, choose $p_0=2$, $m=e=1$ and $A(t_i)=0.10$. Then $p_f=1.9$ and $K_f=1.805$, compared with initial kinetic energy $2$. Reversing $A$ gives $K_f=2.205$: the shifts are asymmetric because of the $A^2/2$ term. A measurement that keeps only the term linear in $A$ would miss $0.005$ in both cases.

For small measurement errors, write $f=L/4-d^2/(4L)$. Its partial derivatives and first-order propagation of uncertainty areHere the $\sigma$ values are standard deviations used to describe uncertainty, and $\mathrm{Cov}(L,d)$ is the covariance. Independent carrier-reading errors give ${\sigma }_d^2={\sigma }_{r_{-}}^2+{\sigma }_{r_{+}}^2$; correlated readings require the term $-2\mathrm{Cov}(r_{-},r_{+})$. Averaging repeated independent pairs reduces the random component of the uncertainty of the mean, but does not remove calibration or model errors.

As a stipulated uncertainty calculation, let $L=45\mathrm{cm}$, $d=15\mathrm{cm}$, ${\sigma }_L=0.10\mathrm{cm}$ and ${\sigma }_{r_{-}}={\sigma }_{r_{+}}=0.05\mathrm{cm}$, with all three input errors independent. Then ${\sigma }_d=0.0707\mathrm{cm}$, $A=5/18$, $B=-1/6$ and ${\sigma }_f\simeq 0.0302\mathrm{cm}$. The result would be $f=10.000\mathrm{cm}$ with this standard uncertainty in the stipulated thin lens model; these numbers are not a measured uncertainty budget. A small neglected principal-plane separation gives an additional first-order bias $f_{\mathrm{thin}}-f\simeq Ah$, distinct from that random budget.

The two kinds of sensitivity must be distinguished. At fixed focal length,These derivatives diverge as $L$ approaches $4f$ from above: a small change in object-screen spacing strongly changes the separation of the two settings. Yet the inferred focal length's displacement sensitivity $|B|=d/(2L)$ tends to zero. That does not make coalescing settings easy to measure. For fixed $L$, define the axial mismatch between the predicted real image and the screen as $g(x)=fx/(x-f)-(L-x)$, for $x>f$. At a focused setting,Near $u=v=2f$, moving the optical lens changes this mismatch only weakly, making visual discrimination of the two settings difficult. At the example's settings $g^{\prime }(15\mathrm{cm})=-3$ and $g^{\prime }(30\mathrm{cm})=3/4$, so equal carrier motion produces different axial defocus responses. This is an axial mismatch calculation, not a full prediction of perceived image sharpness.

Choose $L$ sufficiently above $4f$ to separate the two accepted focusing ranges, while retaining both real images on the available screen and rail. There is no universal best $L$: focusing uncertainty itself depends on aperture, transverse magnification, optical aberrations and the target. Close to $d=0$, the nonnegative measured separation and the quadratic term in $d$ also make first-order propagation of uncertainty unreliable if its error is comparable to its value. Repeatability, controls and the principal plane correction determine whether the elegant two-distance formula is an adequate measurement model.

A biased random walk on the integer line takes unit steps right and left with unequal probabilities $p$ and $1-p$, using independent and identically distributed random variables as increments. For $S_0=0$, the expected value is $n(2p-1)$: positive for $p>1/2$ and negative for $p<1/2$. Changing $p$ to $1-p$ reverses the preferred direction. The exact endpoint probability mass function is given under simple random walk on the integer line.

For example, with $p=0.65$ each step has expected value $0.30$ and variance $0.91$. After $80$ steps,The displacement of the expected value from zero is about $24/8.532\approx 2.813$ standard deviations, yet individual sample paths still move left as well as right, as the figure under simple random walk illustrates. Bias lets this elementary stochastic process model an average directional tendency while retaining fluctuations. Boundary and return questions have their own treatments under gambler's ruin and the specialized children below.

The binomial distribution counts successes in a fixed number $n$ of independent Bernoulli trials with common success probability $p$. Equivalently,Here $n$ is a nonnegative integer and $0\le p\le 1$. The binomial coefficient $\left(\genfrac{}{}{0ex}{}{n}{k}\right)=n!/[k!(n-k)!]$ counts placements of the $k$ successes; each particular arrangement has probability $p^k(1-p{)}^{n-k}$. The binomial theorem makes the masses sum to $[p+(1-p){]}^n=1$. At $p=0$ all mass is at zero, and at $p=1$ all mass is at $n$; these endpoint laws can be read directly without interpreting $0^0$.

The linearity of expectation and addition of variances for independent random variables giveThe fixed trial count, common success probability, and joint independence are substantive model assumptions. Sampling uniformly without replacement from a finite population instead leads to a hypergeometric distribution. Independent trials with different success probabilities generally give a Poisson binomial distribution. Common marginal success probabilities alone are insufficient: if all ten indicators equal the same Bernoulli random variable with $p=0.3$, their sum is zero with probability $0.7$ or ten with probability $0.3$. Its mean is still three, but its variance is $100(0.21)=21$. The standard construction is given in Penn State STAT 504, binomial distribution.

For an exact constructed example, let each of ten independent trials choose uniformly from ten labelled tokens, with replacement, and mark three tokens as successes. Then $p=3/10$ and $K\sim \mathrm{Bin}(10,0.3)$. Exactly three successes have probabilityThe mean is $10(0.3)=3$, the variance is $10(0.3)(0.7)=2.1$, and the standard deviation is $\sqrt{2.1}\approx 1.4491$. Zero successes have probability $0.7^{10}=0.0282475249$, so at least one success has probability $1-0.7^{10}=0.9717524751$. A single count calculation replaces enumeration of all $10^{10}$ equally likely token sequences.

Each bar's height is one probability mass. Three is the most likely count, but its mass is only about $26.7\%$: a most likely outcome need not occur most of the time.

Differentiate Bragg's law at fixed correctly assigned diffraction order:Divide the second equality by the first. For small signed changes,The Bragg angle error $\delta \theta$ must be in radians. An error quoted for detector deflection $2\theta$ must first be halved. The expression is a first-order Taylor expansion; an unknown diffraction order is a discrete assignment error, not an infinitesimal perturbation.

For error bounds, add magnitudes. For independent zero-mean random errors with standard deviations ${\sigma }_d,{\sigma }_{\theta }$, variance addition instead givesCorrelated errors require a covariance term. A shared calibration bias does not become independent merely because many line positions are measured.

With a chosen angular error $10$ arcseconds, $\delta \theta =4.84814\times 10^{-5}\mathrm{rad}$, the angular contribution is $8.3972\times 10^{-5}$ at $30^{\circ }$ and $1.7646\times 10^{-5}$ at $70^{\circ }$. Larger Bragg angles reduce this contribution because $\cot \theta$ decreases. For a nearby spectral doublet, the inverse relation is $\Delta \theta \simeq \tan \theta \Delta \lambda /\lambda$: a fixed fractional separation spreads farther apart angularly at large Bragg angle.

This does not permit unlimited improvement. The geometry requires $n\lambda \le 2d$; a desired wavelength might not have an accessible large-angle reflection with the chosen crystal. Higher diffraction orders can reach larger Bragg angles but their strengths depend on the crystal lattice structure factor. Finite crystal response, slit widths and source size remain relevant. An error in $d$, including uncorrected thermal expansion, remains even when the angular term is small. Very near $90^{\circ }$, the leading angular derivative vanishes and the second-order Taylor expansion matters.

Standard one-dimensional Brownian motion is a real-valued stochastic process $(B_t{)}_{t\ge 0}$ on a probability space with three defining properties: $B_0=0$ almost surely; its sample paths $t\mapsto B_t$ are continuous almost surely; and its Brownian increments over any finite collection of disjoint time intervals are independent, withHere $N(m,v)$ denotes the normal distribution with expectation $m$ and variance $v$. The alternative symbol $W_t$ denotes the same kind of stochastic process. Continuity is a simultaneous statement about an entire sample path, rather than a separate assertion at each fixed time. Relative to its natural Brownian filtration, a future Brownian increment is independent of all observations up to its starting time.

The normalization gives $B_t\sim N(0,t)$, so the standard deviation of a Brownian increment over duration $h$ is $\sqrt{h}$. For $h=1$, $0.01$, and $0.0001$, the corresponding standard deviations are $1$, $0.1$, and $0.01$. Time becoming one hundred times shorter makes typical displacement ten times smaller. Brownian scaling expresses this rule for the whole stochastic process, while the Brownian covariance kernel describes how positions at different times remain related despite independent increments.

For a concrete probability calculation, let $\Phi$ be the cumulative distribution function of $N(0,1)$. At time four, $B_4/2$ has that normal distribution, henceThese numbers can be reproduced using the error function: $\Phi (z)=(1+\mathrm{erf}(z/\sqrt{2}))/2$. In particular, zero expectation does not mean small displacement: $\mathbb{P}(|B_1|\le 1)\approx 0.6826894921$, but the same fixed interval contains only about $38.3\%$ of the probability at time four.

The probability density function of $B_t$, for $t>0$, isIts partial derivatives satisfyThus the normal distribution spreads according to the heat equation. Starting at a fixed position $x_0$ replaces $x$ by $x-x_0$. This is a useful bridge between individual random trajectories and the deterministic evolution of their probability density function: the same Brownian motion provides both descriptions.

Physical diffusion uses its own scale. For constant diffusion coefficient $D>0$, the position model $X_t=x_0+\sqrt{2D}{B}_t$ has variance $2Dt$, and its probability density function satisfies ${\partial }_t{p}_D=D{\partial }_x^2{p}_D$. For example, $D=0.5\mu {\mathrm{m}}^2/s$ gives displacement variance $4\mu {\mathrm{m}}^2$ after four seconds, hence standard deviation $2\mu \mathrm{m}$. Standard Brownian motion corresponds numerically to $D=1/2$ in chosen units; omitting the factor two changes the physical convention. Einstein's 1905 paper relates diffusion to the mean squared displacement in precisely this convention. ^[ref]

Continuity does not supply a smooth physical velocity. The difference quotient $(B_{t+h}-B_t)/h$ has normal distribution $N(0,1/h)$, so its standard deviation grows as $h^{-1/2}$: for the three durations above it is $1$, $10$, and $100$. This calculation motivates the nowhere differentiability of Brownian motion, but it alone is not a proof of the simultaneous assertion at every time. The theorem says that, with probability one, no finite derivative exists anywhere along the sample path. A drawn polygon only displays a finite grid; its straight segments cannot establish this theorem. See Theorem 7.1.6 in Probability: Theory and Examples (fifth edition) by Rick Durrett.

For physical particles, Brownian motion is a model of position on suitable observation scales. It does not claim that microscopic velocity is a derivative of this idealized path. Einstein explicitly restricts his assumption of independent increments to intervals that are not too short. ^[ref]

Another way to measure the roughness is quadratic variation. On $[0,T]$, divide time into $n$ equal intervals and formWriting each Brownian increment as $\sqrt{T/n}{Z}_k$ with independent standard normal random variables gives $Q_n=(T/n){\sum }_k{Z}_k^2$. Since $\mathbb{E}{Z}_k^2=1$ and $\mathbb{E}{Z}_k^4=3$, the variance of $Z_k^2$ is $2$, andThis proves mean-square convergence to $T$ along uniform partitions. For $T=2$ and $n=100$, the expectation is $2$, the variance is $0.08$, and the root mean square of the error is $\sqrt{0.08}\approx 0.2828427125$. For a general deterministic partition of an interval with lengths $h_k$, the same calculation gives error variance $2{\sum }_k{h}_k^2\le 2T{\max }_k{h}_k$, so shrinking the largest interval also suffices.

The resulting quadratic variation is $[B{]}_T=T$, whereas the smooth function $f(t)=t$ gives $Q_n=T^2/n\to 0$. This difference explains why the ordinary rules of calculus need the second-order correction in the Itô formula. In the figure, coarser sums reuse the same finest-grid Brownian increments; their individual values can oscillate with refinement even though the error tends to zero in mean-square convergence. Refining a fixed interpolating polygon forever would instead give zero quadratic variation, so the finite-grid illustration must not be confused with the limiting Brownian motion.

Historically, Brown's 1828 report describes microscopic observations made in 1827, starting with particles released from pollen and extending to inorganic materials; this challenged an explanation based solely on living matter. His observations and the later mathematical model are distinct contributions. ^[ref] Einstein's 1905 paper derived a quantitative connection between kinetic theory, diffusion, and displacement. ^[ref] Wiener's 1923 Differential-Space developed integration on spaces of paths, a foundation for Wiener measure; it also acknowledges earlier work and his own previous papers. These dates mark documented contributions rather than an assertion that each subject had a single discoverer. ^[ref]

For every $c>0$, the stochastic process ${\tilde{B}}_t=c^{-1/2}{B}_{ct}$ is standard Brownian motion. Its Brownian increments remain independent and have variance $c^{-1}c(t-s)=t-s$; continuity and the zero starting value are preserved. Equivalently,Equality here is of whole-process distributions, including every finite-dimensional distribution. It is not an equality of the two expressions along a particular realization. With $c=4$, speeding time up fourfold requires dividing displacement by two to recover the standard law. The figure applies exactly these coordinate changes to four finite samples; the shared colors indicate reused samples, while the dashed lines show pointwise standard deviations, not bounds on entire sample paths.

The same square-root normalization explains how a simple random walk approaches Brownian motion. Let $S_k={\sum }_{j=1}^k{\xi }_j$, where the independent steps take values $+1$ and $-1$ with equal probability. Then $R_n(t)=S_{\lfloor nt\rfloor }/\sqrt{n}$ has variance $\lfloor nt\rfloor /n$. Its linear interpolation isand ${\sup }_{0\le t\le 1}|Y_n(t)-R_n(t)|\le 1/\sqrt{n}$. Donsker's theorem says $Y_n\Rightarrow B$ in $C[0,1]$ with the uniform norm: this is convergence in distribution of random paths. The interpolation bound vanishes, but the theorem does not assert pointwise convergence of an arbitrarily fixed walk to an independently chosen Brownian motion. The usefulness is that many-step discrete motion has a universal continuous limit.^[ref]

Alvarez describes a progression from measured coordinates to momenta and angles, then conservation-constrained fits testing alternative track identities. Computer programs made the repeated tests practical; an insufficiently constrained event could retain several hypotheses.^[ref]

A mathematical formulation explains what such a fit can establish. Let $y$ collect measured momenta and angles, with positive-definite covariance matrix $V$. For a chosen reaction hypothesis, adjust them to $x$ by minimizing a generalized least squares objective,The constraints $g$ express four-momentum conservation and any independently known daughter masses. For each proposed identity, calculate energy from the energy–momentum relation, $E_i^2=m_i^2{c}^4+c^2|{\mathbf{p}}_i{|}^2$. A missing neutral particle adds unknowns; a measured decay vertex can add geometric constraints. Compatibility depends on the measurement errors and the independent constraints, not just the smallest raw ${\chi }^2$ among unrelated fits.

A simple independent check uses a single momentum conservation constraint. Suppose two measured one-dimensional momenta are $0.32$ and $-0.28\mathrm{GeV}/c$, each with independent standard deviation $0.02\mathrm{GeV}/c$, and a stipulated parent at rest requires their sum to vanish. Minimizing the two equal-weight squared changes gives fitted values $0.30$ and $-0.30\mathrm{GeV}/c$, with ${\chi }^2=2$. A fit repairs a mismatch according to a hypothesis; it does not prove that hypothesis. This one-component illustration is not a complete reaction fit.

A burette delivers an adjustable measured volume of a chemical solution; the change in its reading gives the delivered volume in an acid-base titration.

Its downward-increasing scale measures a change in reading, not the initial liquid inventory. For initial reading $R_i$ and final reading $R_f$ without refilling, the nominal delivered volume is $v=R_f-R_i$. An illustrative $2.35\mathrm{mL}$ to $12.35\mathrm{mL}$ reading therefore represents $10.00\mathrm{mL}$. There is no requirement to start at zero. An air bubble disappearing from the tip uses part of the apparent delivery to fill the tip, and a hanging drop has not yet entered the sample; both break the intended volume account.

For transparent aqueous samples, read the bottom of the liquid meniscus at eye level; viewing from above or below shifts its apparent alignment. The mark is a comparison between a curved liquid surface and a scale, not a direct measure of amount of substance. Eye-level reading and the distinction between contained and delivered volumes are part of solution preparation and volumetric calibration.^[ref]

Both readings contribute measurement uncertainty. If each has independent standard deviation ${\sigma }_R$, their difference has standard deviation $\sqrt{2}{\sigma }_R$; a common additive reading offset cancels, while a scale-factor error remains. For the inventory $n=c_{B}v$, uncertainty in titrant concentration remains even with perfect readings. This is why repeat delivery, independent volume calibration, and titrant standardization answer different questions.

Can two readout values locate a hit between wires? Only with a specified response model. Place two readout wires at $0$ and $p$, put $u=x/p\in [0,1]$, and stipulate calibrated, noiseless electric chargesThese are assumptions for a toy sensor response, not the natural response law of adjacent anode wires. Under these assumptions the centroidrecovers position exactly. The apparent gain in precision comes from adding calibrated analogue information. Taking the mean of two fired wire addresses without their electric charges would instead give the fixed midpoint $p/2$.

For a chosen noise model, replace the two electric charges by $q_L+n_L$ and $q_R+n_R$, with independent random variables $n_L,n_R$ of zero expected value and common standard deviation $s$. Assume $s/Q\ll 1$ and no threshold censoring. First-order expansion of the ratio givesFor stipulated $p=2\mathrm{mm}$ and $s/Q=0.02$, this approximation gives $0.0283\mathrm{mm}$ at the midpoint and $0.0400\mathrm{mm}$ at either endpoint. The endpoint values describe uncensored signed channel readouts, including small negative values after baseline subtraction. They do not assume that both channels would pass a positive threshold. A ratio with a fluctuating denominator also need not have exactly zero bias; this is a local small-noise approximation.

Three original counterexamples show why calibration and selection matter:

Likewise, two equal hits at $u=1/4$ and $u=3/4$ have exactly the same total two-channel response as one hit at $u=1/2$ with twice the electric charge. A centroid can hide multiplicity even when each charge measurement is exact.

Diffusion enabled the lecture's Orsay interpolation, without establishing this model's linear response.^[ref]

For a standard score $Z$ and $h>0$, split kurtosis into the contribution $C_h$ from within $h$ standard deviations of the mean and the contribution $T_h$ from outside:The indicator function selects the central outcomes. The bound follows directly from $Z^4\le h^4$ there. Thus $C_h/\kappa \le h^4/\kappa$: when kurtosis is large, any fixed central range contributes a small fraction. This is a precise reason to emphasize extreme standard scores rather than peak height.^[ref]

For the rare-outlier three-point distribution with $q=1/100$, $C_1=0$ and $T_1=100$. At the other endpoint $q=1$, all probability mass lies at $Z=\pm 1$, so $C_1=\kappa =1$: central dominance is possible when kurtosis is small. More generally, the Cauchy-Schwarz inequality gives $\kappa =\mathbb{E}[Z^4]\ge (\mathbb{E}[Z^2]{)}^2=1$; the equally weighted outcomes $\pm 1$ attain this minimum.

Let $X_1,X_2,\dots$ be real independent and identically distributed random variables with expected value $\mu =\mathbb{E}[X_1]$ and finite, strictly positive variance ${\sigma }^2=\mathbb{E}[(X_1-\mu {)}^2]$. Write $S_n={\sum }_{j=1}^n{X}_j$ and ${\overline{X}}_n=S_n/n$ for their sum and sample mean. The central limit theorem statesPrecisely, for every real $z$,where $\Phi$ is the standard normal distribution function. The conclusion is convergence in distribution to the standard normal distribution. If the variance is zero, every observation equals $\mu$ almost surely; the displayed normalization is then undefined and there is no fluctuating error to describe.

The striking feature is that the observations need not have a normal distribution. They may be binary, bounded, or strongly skewed. Their detailed probability distribution disappears from this limit except through the expected value and variance. The subtraction removes predictable growth: $\mathbb{E}[S_n]=n\mu$. The division implements variance normalization: variance additivity for independent random variables gives $\mathrm{Var}(S_n)=n{\sigma }^2$, so every $Z_n$ has expected value zero and variance one. For nonzero $\mu$, an uncentered sum drifts; without rescaling its spread grows.

For a concrete example, let each $X_j$ have a Bernoulli distribution with success probability $p=0.2$. Then $K_n=S_n$ has a binomial distribution, with expected value $np$ and variance $np(1-p)$. Its exact probability mass function and its normalized count areHere $\left(\genfrac{}{}{0ex}{}{n}{k}\right)$ is a binomial coefficient. Individual trials are conspicuously skewed, yet their centered, rescaled counts approach the same standard normal distribution as sums of fair trials.

At $n=100$, $K_n$ has expected value $20$ and standard deviation $4$. The exact binomial distribution gives $\mathbb{P}(16\le K_n\le 24)=0.740141\dots$. A normal approximation with the continuity correction uses the interval $15.5$ to $24.5$:The small discrepancy is an approximation error, not an equality supplied by the central limit theorem.

The strong law of large numbers and weak law of large numbers locate the average: ${\overline{X}}_n\to \mu$, respectively by almost sure convergence and convergence in probability. The central limit theorem magnifies that shrinking error and describes its shape. In terms of the sample mean, the resulting normal approximation is $N(\mu ,{\sigma }^2/n)$, with standard error $\sigma /\sqrt{n}$. This explains why four times as many independent observations halve the standard error. When $\sigma$ is known, ${\overline{X}}_n\pm 1.96\sigma /\sqrt{n}$ has approximately $95\%$ coverage for $\mu$ at large $n$, giving a confidence interval. This coverage refers to repeated samples, not a guarantee for an individual interval.

A useful application is the Monte Carlo method: estimate a success probability by the proportion of successful simulated trials. Nicholas Metropolis and S. Ulam's original paper, "The Monte Carlo Method" (1949), illustrates this with games of solitaire and estimation of a region's volume by randomly selected points, and discusses why estimates require probabilistic error control. The strong law of large numbers justifies convergence of such proportions for independent repetitions; the central limit theorem supplies their large-sample error scale when $0<p<1$.^[ref]

There are three useful limits to this picture:

A historical original-paper reference is J. W. Lindeberg's "Eine neue Herleitung des Exponentialgesetzes in der Wahrscheinlichkeitsrechnung" (1922), published in Mathematische Zeitschrift 15, pp. 211–225. The Lindeberg condition and Lindeberg-Feller central limit theorem below describe a version for sums whose summands need not share a probability distribution.^[ref]

A brief proof outline uses the characteristic function of $Y=(X_1-\mu )/\sigma$. Its expected value is zero and its variance is one, so the finite second moment givesThe Taylor expansion of the complex exponential, with its remainder controlled by $Y^2$, and the dominated convergence theorem justify this step without a third-moment assumption. Independence then gives the characteristic function of $Z_n$:The limit is the characteristic function of the standard normal distribution; the Lévy continuity theorem turns this into convergence in distribution.

For a real random variable $X$ with finite variance ${\sigma }^2$, expected value $\mu$, and a threshold $a>0$, Chebyshev's inequality givesOne can also cap this upper bound at one, since a probability cannot exceed one. If $\sigma >0$, taking $a=k\sigma$ gives $\mathbb{P}(|X-\mu |\ge k\sigma )\le 1/k^2$. For example, at least $3/4$ of the probability lies strictly within two standard deviations of the mean. This statement needs only the finite variance assumption. If the variance is zero, $X=\mu$ almost surely and every positive deviation threshold has zero tail probability.

Consider a random variable uniformly distributed on $\{-2,-1,0,1,2\}$. Its expected value is zero and its variance is $(4+1+0+1+4)/5=2$. At $a=2$ the actual tail probability is $2/5$, while Chebyshev's inequality gives $2/2^2=1/2$. At $a=3$ the actual tail probability is zero, while the bound is $2/9$. A valid bound need not describe the exact shape or bounded support of a particular probability distribution.

For the standard normal distribution, the actual tail probabilities beyond two and three standard deviations are approximately $0.0455003$ and $0.0026998$, whereas the general bounds are $1/4$ and $1/9$. More knowledge of the probability distribution can therefore improve the estimate. The comparison does not diminish the guarantee: the three-point example in equality in Chebyshev's inequality attains it exactly.

The proof is a one-line use of Markov's inequality on the nonnegative random variable $Y=(X-\mu {)}^2$. The event $|X-\mu |\ge a$ equals $Y\ge a^2$, and $\mathbb{E}Y={\sigma }^2$. This explains why the denominator is squared and why a finite variance supplies the necessary information. ^[ref].

One important use is control of a sample mean. For independent and identically distributed random variables with finite variance ${\sigma }^2$, the variance of a sample mean is ${\sigma }^2/n$. HenceFor every fixed tolerance, the bound tends to zero as the sample size grows, proving the finite-variance case of the weak law of large numbers. For ${\sigma }^2=1$ and $\epsilon =0.1$, $n=2000$ makes the bound $0.05$. This is a sufficient bound under the stated assumptions, rather than an optimal required sample size.

An early historical source is Pafnuty Chebyshev's 1867 paper “Des valeurs moyennes”, published in French translation by N. de Khanikof. Its opening theorem bounds deviations of sums using their expected values and second moments, linking averages with control of exceptional outcomes. ^[ref].

If two independent approximately Gaussian mean lifetime estimates have known standard deviations $s_1,s_2$, equality of their true lifetimes implies that their difference has mean zero and variance $s_1^2+s_2^2$. Its squared standardized difference therefore has a one-degree-of-freedom chi-squared distribution, giving the statistic in the title. A chosen significance threshold converts its upper-tail P-value into a compatibility decision.

Testing a fitted decay curve against a fixed numerical lifetime is a different question: that test does not include the external measurement uncertainty. Comparing the two experiments requires both uncertainties or a joint likelihood with a common lifetime. Correlated estimates replace the denominator by $s_1^2+s_2^2-2\mathrm{Cov}({\hat{\tau }}_1,{\hat{\tau }}_2)$.

Suppose each Cartesian coordinate of a particle of mass $m$ is supported in an interval of length $L$, in $d$ dimensions. Popoviciu's inequality on variances gives $\Delta x_i\le L/2$, and the Heisenberg uncertainty relation gives $\Delta p_i\ge \hslash /L$. Since $⟨p_i^2⟩=(\Delta p_i{)}^2+⟨p_i{⟩}^2$,This is a rigorous but non-sharp confinement estimate. For a cubical infinite square well, separation of variables gives the stronger exact ground-state value $d{\pi }^2{\hslash }^2/(2m{L}^2)$. One must not identify a support length with a uniquely prescribed standard deviation or claim that the uncertainty estimate is the exact ground-state energy.

Logarithmic differentiation of the diphoton mass gives $dM/M=(d{E}_1/E_1+d{E}_2/E_2)/2$ at fixed angle. If $s_i$ are fractional energy standard deviations and $\rho$ their correlation coefficient, covariance propagation gives the displayed formula. A common multiplicative calibration error transfers directly to the reconstructed mass; independent errors combine differently. Angular errors, if present, add the corresponding derivative and covariance terms.