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.

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]

The complementary error function describes the Gaussian integral area remaining beyond the endpoint. For a real $x$,It is an entire function for complex arguments when defined by $1-\mathrm{erf}(z)$. On the real line it decreases from two to zero, with $\mathrm{erfc}(0)=1$ and $\mathrm{erfc}(-x)=2-\mathrm{erfc}(x)$.^[ref]

For example, $\mathrm{erfc}(1)\approx 0.157299$ is a sizeable tail, whereas $\mathrm{erfc}(8)\approx 1.12243\times 10^{-29}$ is extremely small. In the tested floating-point arithmetic, a direct evaluation of $1-\mathrm{erf}(8)$ gives zero because the computed error function has already rounded to one. Evaluating the complementary error function directly preserves this tail. It is the useful quantity in an upper-tail normal distribution probability and in the decaying side of an error-function diffusion front.

For large positive real $x$, the leading asymptotic expansion isThe first corrections are $1-1/(2x^2)+3/(4x^4)-\cdots$ multiplying the leading term. This is an asymptotic expansion, not an assertion that summing arbitrarily many terms converges at a fixed $x$.^[ref] The complementary error function tail bound provides finite-$x$ inequalities instead.

For a diffusion-controlled mushy layer of constant solid fraction $\varphi$, a sharp moving edge $h=2\lambda \sqrt{\kappa t}$ reduces the thermal heat equation to ${\theta }^{\prime \prime }+2\eta {\theta }^{\prime }=0$ and the exterior solute diffusion equation to ${\mathcal{C}}^{\prime \prime }+2Le\eta {\mathcal{C}}^{\prime }=0$, where $Le=\kappa /D$ is the Lewis number. Wall values select an error function in the mush; far-field decay selects complementary error functions in the liquid. Continuity fixes their amplitudes in terms of the interfacial temperature and $\lambda$. The Stefan condition, solute conservation and marginal equilibrium at a mush-liquid boundary give three algebraic/transcendental equations for these two constants and $\varphi$. Derivation follows by ${\eta }_t=-\eta /(2t)$ and ${\eta }_z=1/(2\sqrt{\kappa t})$, followed by one integration $Y^{\prime }=K{e}^{-a{\eta }^2}$. Dividing interface equations can discard degenerate branches, so the undivided equations must be retained when a denominator vanishes. Physical admissibility additionally requires $\lambda >0$ and $0\le \varphi \le 1$.

The error function records the signed area under a Gaussian function between zero and a variable endpoint:For a negative real number $x$, reversing the integration limits makes this integral negative. The normalization uses the Gaussian integral: the full area of $e^{-t^2}$ is $\sqrt{\pi }$, so its positive half contributes one after multiplication by $2/\sqrt{\pi }$. The same definition extends to a complex number $z$ along any path from zero to $z$, giving an entire function.^[ref]

For example, $\mathrm{erf}(0)=0$, $\mathrm{erf}(1)\approx 0.842701$, and $\mathrm{erf}(-1)\approx -0.842701$. A finite positive endpoint captures only part of the positive Gaussian integral; $\mathrm{erf}(x)$ approaches one as $x\to +\infty$ but never reaches it at a finite real $x$.

The useful idea is to name a recurring integral once instead of carrying it through every calculation. The error function and the normal distribution express probabilities and quantiles through this area, while an error-function diffusion front describes how diffusion smooths a sharp jump. The complementary error function retains the small area left in the tail.

The fundamental theorem of calculus givesThus the real-axis error function is strictly increasing, with values in $(-1,1)$. Its integrand is an even function, so the error function is an odd function, $\mathrm{erf}(-x)=-\mathrm{erf}(x)$. The derivative also gives ${\mathrm{erf}}^{\prime \prime }(x)=-4x{e}^{-x^2}/\sqrt{\pi }$: the curve is concave for $x>0$ and convex for $x<0$. These real-axis properties do not imply boundedness throughout the complex plane; the existing quadrant asymptotics of the error function describe a different behavior in other directions.

For a standard normal random variable $Z$, its standard normal distribution function isThe change of variables $t=u/\sqrt{2}$ and the symmetric distribution of the standard normal random variable give the second expression. For $X\sim N(\mu ,{\sigma }^2)$ with $\sigma >0$, standardization of a normal random variable then givesIn particular, the intervals within one, two, and three standard deviations have calculated probabilities about $0.682689$, $0.954500$, and $0.997300$. These values depend on a normal distribution; they are not universal coverage guarantees for an arbitrary random variable. This is the direct connection between a Gaussian function area and the familiar central intervals.

For a small upper-tail probability, use the complementary error function directly:Computing this as $1-\Phi (z)$ instead can lose the tail to catastrophic cancellation. The Gaussian tail bound gives another useful estimate without requiring an exact error function value.

Consider the one-dimensional diffusion equation on the whole real line,with initial value one for $x<0$ and zero for $x>0$. This is an idealized Heaviside step function profile, not a measured material or an apparatus specification. The heat kernel gives a convolution with the initial profile:The change of variables $s=(x-y)/(2\sqrt{Dt})$ turns the convolution into a Gaussian integral tail. The resulting similarity variable $\eta =x/(2\sqrt{Dt})$ shows that the transition width grows like $\sqrt{Dt}$: multiplying elapsed time by four doubles every fixed-profile distance from the original interface.

At every positive time, $u(0,t)=1/2$, while $u(x,t)\to 1$ as $x\to -\infty$ and $u(x,t)\to 0$ as $x\to +\infty$. For every fixed $x\ne 0$, the limit as $t\to 0^{+}$ returns the specified initial value. The value assigned at the single initial discontinuity does not change this convolution. The error function therefore explains smoothing without pretending that a discontinuous initial profile remains sharp.

Direct differentiation checks the partial differential equation. Writing $\eta =x/(2\sqrt{Dt})$ givesThe model assumes constant diffusion coefficient, no advection, and no finite-domain boundary effects. For an interpretation as heat conduction, the heat equation applies to a normalized temperature difference under its own material assumptions.

For conductivity $\lambda ={\lambda }_{0}T$, no cooling in the hot phase and a zero-temperature boundary at $m=0$, the mass-coordinate equation is linear diffusion. Its similarity profile solves $f^{\prime \prime }+\xi f^{\prime }/2=0$, $f(0)=0$, $f(\infty )=1$, giving the displayed error function. The thermal layer broadens as the square root of time. Conducted energy enters the cold radiating reservoir at rate per area $L={\lambda }_{0}A{T}_m(0)=A{T}_0\sqrt{{\lambda }_0}/\sqrt{\pi \tau }$, hence $L\propto t^{-1/2}$.

The Faddeeva function is an entire function built from the complementary error function. Its integral representation for $\mathrm{Im}z>0$ is $w(z)=(i\pi {)}^{-1}{\int }_{-\infty }^{\infty }{e}^{-t^2}/(t-z)dt$, as normalized in NIST's integral representation. Its reflection identity is $w(-z)=2e^{-z^2}-w(z)$, directly from the odd function property of the error function.

The Gaussian integral gives the exact area under a Gaussian function. For a positive real number $a$,This is a finite improper integral over the whole real line. For example, $I(1)=\sqrt{\pi }\approx 1.772454$ and $I(4)=\sqrt{\pi }/2\approx 0.886227$. Increasing $a$ from one to four narrows the curve and halves its area, while its height at zero stays one. The change of variables $u=\sqrt{a}x$ explains the factor $a^{-1/2}$.

The useful consequence is an exact normalization. For $a>0$ and any real center $\mu$,This is a probability density function for a normal distribution with expected value $\mu$ and variance $1/(2a)$. A narrower probability density function must be taller to keep total probability one. In the usual notation $a=1/(2{\sigma }^2)$, the prefactor is $1/(\sigma \sqrt{2\pi })$, where $\sigma$ is the standard deviation.

The same area calculation makes the one-dimensional heat kernelintegrate to one in $x$, with $D$ the diffusion coefficient. A spreading Gaussian function can therefore preserve total mass while its height falls. This links the Gaussian integral with both the normal distribution and the diffusion equation.

Finite intervals lead to the error function, rather than this single constant:The complementary tail is described by the complementary error function. These definitions and their real-axis limits are given in the NIST reference. ^[ref]

The polar-coordinate proof of the Gaussian integral explains why $\pi$ appears in an integral of a function defined on a line: squaring its area creates a rotationally symmetric integral on a plane. Substitution $t=x^2$ also gives $\Gamma (1/2)=2{\int }_0^{\infty }{e}^{-x^2}dx=\sqrt{\pi }$, connecting the result with the gamma function.

The normalized Gaussian mass of an interval isIt is the convolution of an interval indicator with the heat kernel, and has an error function representation. As the variance tends to zero it converges to the interval indicator away from its endpoints, with value one half at an endpoint. This is a concrete approximate identity.

For the upper-pole integral $I={\int }_{{\mathcal{C}}_{+}}{t}^{\gamma -1}{e}^{\lambda (1-t)}(1-t{)}^{-1}dt$, let $\gamma \to \infty$ and $\lambda =\gamma +i\mu \sqrt{\gamma }+O(1)$. Scaling $t=1+s/\sqrt{\gamma }$ gives $-{\int }_{{\mathcal{C}}_{+}}{e}^{-s^2/2-i\mu s}ds/s$ at leading order. At $\mu =0$ its indentation is $i\pi$. Differentiation with respect to $\mu$ gives $i\sqrt{2\pi }{e}^{-{\mu }^2/2}$, so integration gives the displayed error function. This is an optimal truncation transition; the underlying analytic integral has no literal discontinuity. The general mechanism is discussed in the DLMF account of Stokes smoothing.

For $z=2n+i\mu \sqrt{n}+O(1)$, the factor $t^{2n}{e}^{z(1-t)}$ has saddle scale $t=1+u/\sqrt{n}$ and becomes $e^{-u^2-i\mu u}$. If the amplitude has a simple pole there, its leading local integral contains $e^{-u^2-i\mu u}/u$. Differentiating its Cauchy principal value with respect to $\mu$ reduces it to a Gaussian integral, producing an error function. The nonsingular amplitude $(1+t{)}^{-1}$ gives $\frac{1}{2}\sqrt{\pi /n}{e}^{-{\mu }^2/4}$.

For initial data $u_0$ on the real line, Fourier transformation or convolution with the fundamental solution gives

For a bounded function $u_0$, the convolution is well defined at every positive time because the heat kernel has integral one. The error-function diffusion front uses a Heaviside step function initial condition and gives one smoothed interface. The heat equation with interval-indicator initial data instead has two interfaces; subtracting their two error functions produces the finite-interval profile.

In ${\epsilon }^2{y}^{\prime \prime }+a(x)y^{\prime }+\cdots =f$, a simple zero $a(x_{\ast })=0$, $a^{\prime }(x_{\ast })\ne 0$, balances diffusion and advection on width $\epsilon$. With $a(x)\simeq 2(x-x_{\ast })$, the leading homogeneous inner equation is $Y^{\prime \prime }+2z{Y}^{\prime }=0$, whose bounded solutions are constants and the error function. Thus distinct order-one outer solutions can match across an interior transition. Reversing the diffusion sign makes the nonconstant homogeneous solution grow at both ends, changing the matching structure.

Because the real error function is strictly increasing from $-1$ to $1$, its real inverse function maps $(-1,1)$ onto the real line. For $0<q<1$, the standard normal quantile is thereforeMore generally, the $q$-quantile of a normal distribution with expected value $\mu$ and positive standard deviation $\sigma$ is $\mu +\sigma \sqrt{2}{\mathrm{erf}}^{-1}(2q-1)$. This reverses the area calculation: instead of asking how much probability lies below an endpoint, a quantile asks which endpoint encloses a specified probability. The domain restriction matters; $q=0$ and $q=1$ correspond to infinite endpoints, rather than finite real values of the inverse error function.

The Neumann solution places a planar phase boundary at $h(t)=2\lambda \sqrt{\kappa t}$. Similarity reduction of the heat equation gives temperature profiles in terms of the error function, while the Stefan condition determines $\lambda$ implicitly.

The normalized quantum harmonic oscillator ground state has $|A{|}^2=\sqrt{m\omega /(\pi \hslash )}$ and energy $E_0=\hslash \omega /2$. Its classical turning points are $x=\pm \sqrt{\hslash /(m\omega )}$. With $y=x\sqrt{m\omega /\hslash }$, the ground-state oscillator probability beyond the turning points isThe stated rounded error function gives $0.16$; the more precise Gaussian integral is $0.15730$.

On bounded intervals the exponential series converges uniformly, allowing termwise integration. Since $e^{-t^2}=1-t^2+t^4/2+O(t^6)$,The requested first two nonzero terms are thereforeThe function is the error function, with $I(0)=0$ and derivative $2e^{-x^2}/\sqrt{\pi }$.

The fundamental theorem of calculus differentiates the variable upper endpoint:At zero, $I(0)=0$, $I^{\prime }(0)=2/\sqrt{\pi }$, $I^{\prime \prime }(0)=0$ and $I^{\prime \prime \prime }(0)=-4/\sqrt{\pi }$. The integrand is even, so its integral from zero is an odd function; all even coefficients in the Maclaurin series vanish. HenceThis is the local expansion of the error function.