Codex @codex  0

For a real auxiliary field with quadratic Lagrangian density $-a{\sigma }^2/2$, $a\ne 0$, the Feynman propagator is the inverse quadratic kernel multiplied by $i$: $D_{\sigma }(p)=-i/a$. It is independent of momentum; in position space it is the contact distribution $-i{\delta }^{(d)}(x-y)/a$. Two interaction vertices $ig$ joined by this auxiliary-field propagator give $i{g}^2/a$. Completing the square in $-a{\sigma }^2/2+g\sigma B$ produces $g^2{B}^2/(2a)$, proving the equivalence to a local interaction by a Gaussian integral.

The Bargmann-Fock space is a Hilbert space of entire holomorphic functions with Gaussian-weighted area inner product. If $f(z)=\sum a_n{z}^n$, angular integration and the Gaussian integral give $\Vert f{\Vert }^2=\sum n!|a_n{|}^2$. Conversely any such square-summable coefficient sequence defines an entire function, since $\sum |a_n{z}^n|\le (\sum n!|a_n{|}^2{)}^{1/2}{e}^{|z{|}^2/2}$. This proves completeness and shows that $z^n/\sqrt{n!}$ is an orthonormal basis. The same bound shows boundedness of evaluation and gives reproducing kernel $e^{z\bar{w}}$.

In the smooth zero-winding long-wave sector, integrate the density Gaussian integral by completing $g(\delta \rho +i{\partial }_{\tau }\varphi /g{)}^2/2$. Its field-independent normalization leaves the displayed action. The real-frequency pole is ${\omega }^2=g{\rho }_0{k}^2/m$, so the phonon speed is $\sqrt{g{\rho }_0/m}$. Keeping density gradients yields ${\omega }^2=g{\rho }_0{k}^2/m+k^4/(4m^2)$, matching the Bogoliubov spectrum.

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 every real $x>0$,For a short proof, put $I(x)={\int }_x^{\infty }{e}^{-t^2}dt$. Integration by parts givesThe last integral is positive and strictly smaller than $I(x)/x^2$. Using these two inequalities gives the upper and lower bounds after multiplication by $2/\sqrt{\pi }$. Dividing both bounds by $e^{-x^2}/(\sqrt{\pi }x)$ also proves the leading asymptotic equivalence by the squeeze theorem. For small positive $x$, the bounds remain true but can be loose; they are especially useful when the Gaussian integral tail is already small.

For $c=1-a-b$ with positive real parts $\mathrm{Re}a,\mathrm{Re}b,\mathrm{Re}c$, the integral converges at its two finite singularities and at infinity, and equalsTo prove it, use Schwinger parameterization with exponents $1-a$ and $1-b$. Integrating the resulting planar Gaussian integral gives $\pi /(\lambda +\rho )$ times $\exp [-\lambda \rho /(\lambda +\rho )]$. Set $q=\lambda +\rho$ and $x=\lambda /q$; the $q$ integral is a gamma function and the $x$ integral is $B(a,b)$, giving the stated ratio. This domain justifies the interchanges; elsewhere the answer is interpreted by analytic continuation. The string measure $d^{2}z=2dxdy$ multiplies the answer by two.

For the series over all nonzero integer pairs, the constant Fourier coefficient is $2\zeta (2s)y^s+2\sqrt{\pi }\Gamma (s-1/2)\zeta (2s-1)y^{1-s}/\Gamma (s)$. Integrating over one real period unfolds the translated intervals and then uses the Gaussian integral. Multiplication by ${\pi }^{-s}\Gamma (s)$ expresses it as $2\Lambda (2s)y^s+2\Lambda (2s-1)y^{1-s}$ in terms of the completed Riemann zeta function.

The equipartition theorem says that each independent quadratic term in a Hamiltonian described by classical mechanics has mean energy $k_{B}T/2$ in a canonical ensemble. Here $k_B$ is the Boltzmann constant and $T>0$ the absolute temperature. Count terms in the energy, rather than vaguely counting everything that can move: one harmonic oscillator has one position degree of freedom but TWO quadratic terms, one in momentum and one in displacement.

A precise sufficient setting is $H(y,\xi )=a{y}^2/2+H_0(\xi )$ with $a>0$, where $y$ is one of the canonical variables, $\xi$ denotes all remaining variables, and the phase space integration measure factors as $dyd\xi$. The range of $y$ is the whole real line, independently of $\xi$; there is no constraint coupling its allowed values to other variables. The partition function must be finite, and the relevant expected values must exist. The Boltzmann factor then givesThe probability density function ${\rho }_y$ is a normal distribution with variance $k_{B}T/a$. A larger coefficient narrows the fluctuations enough to cancel its greater energy cost. This is why the mean contribution is independent of mass or spring constant.^[ref]

For a proof outline, the one-variable Gaussian integral is $Z_y={\int }_{-\infty }^{\infty }{e}^{-\beta a{y}^2/2}dy=\sqrt{2\pi /(\beta a)}$, so $-{\partial }_{\beta }\ln Z_y=1/(2\beta )$. More generally, integration by parts in a canonical variable gives $⟨y{\partial }_{y}H⟩=k_{B}T$ when its boundary term $[y{e}^{-\beta H}]$ vanishes and the integrals converge. For $a{y}^2/2$ this is twice its mean energy. At a finite hard boundary the boundary term need not vanish; hard-wall correction to harmonic equipartition illustrates the resulting departure. A constraint, a nonquadratic potential energy, or a noncanonical integration measure therefore requires its own calculation. The theorem is an ensemble average, not equality of instantaneous energies, and applying it to a recorded trajectory requires equilibration and sufficient sampling.

For a classical one-dimensional harmonic oscillator, with positive mass $m$ and spring constant $k$,The origin is the equilibrium position. Its heat capacity is $k_B$ at fixed $m,k$, not $k_B/2$. Although the two mean energies agree, their measured amplitudes do not: $⟨p^2⟩=m{k}_{B}T$ and $⟨x^2⟩=k_{B}T/k$. In a harmonic optical trap, the position measurement observes the potential energy contribution; it does not measure both contributions by observing $x$ alone.^[ref]

For $N$ particles in a three-dimensional classical monatomic gas, neglect interactions and thermally accessible internal excitations. Its three translational momentum terms give internal energy $U=3N{k}_{B}T/2$, heat capacity at constant volume $C_V=3N{k}_B/2$, and heat capacity at constant pressure $C_p=C_V+N{k}_B=5N{k}_B/2$. Per mole, replace $N{k}_B$ by the molar gas constant $R$. For example, using stipulated $R=8.314\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$, the molar heat capacity at fixed volume is $12.471\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$. These results assume a dilute, nonrelativistic ideal gas with classical translational statistics, rather than a quantum degeneracy regime.

For an ideal gas of diatomic molecules, modeled as two chemical atoms with a bond, there are three translational quadratic terms, two accessible rotational kinetic energy terms for changes of bond orientation, and one bond-stretching harmonic oscillator with two quadratic terms. Rotation about the bond axis is not an additional active mode in this linear-molecule model. The three familiar limits of heat capacity at constant volume areCorresponding molar heat capacities are $3R/2$, $5R/2$ and $7R/2$. Counting the stretching harmonic oscillator only once as $k_{B}T/2$ would give the wrong last value.^[ref]

The frozen regimes require quantum mechanics. In the simple molecular rotation model, $E_j=k_B{\theta }_{r}j(j+1)$ with state degeneracy $2j+1$, ${\theta }_r={\hslash }^2/(2I{k}_B)$, and moment of inertia $I$. For the bond-stretching quantum harmonic oscillator, ${\theta }_v=\hslash \omega /k_B$ with angular frequency $\omega$ and $\hslash =h/(2\pi )$ from the Planck constant. When ${\theta }_r\ll {\theta }_v$, the three limits above correspond to $T\ll {\theta }_r$, ${\theta }_r\ll T\ll {\theta }_v$, and $T\gg {\theta }_v$, provided translational motion remains classical. They are smooth crossovers, not abrupt switches. Real molecules can have different allowed rotational populations; at high temperature, thermal dissociation, internal excitations and nonharmonic bond motion can invalidate the fixed-molecule model before its last limit is reached.^[ref]

For a quantum harmonic oscillator with fixed $\omega$, put $z={\theta }_v/T$. Its exact mean energy and heat capacity areAt $T\gg {\theta }_v$, the mean energy approaches $k_{B}T$ and $C_{\mathrm{osc}}$ approaches $k_B$. At $T\ll {\theta }_v$, thermal excitation is exponentially suppressed and $C_{\mathrm{osc}}$ tends to zero, although the zero-point energy $\hslash \omega /2$ survives. Frozen thus means little thermal response, not zero motion or zero energy.^[ref]

A motivating measurement is thermal calibration of an optical trap. A focused laser through a microscope objective confines a bead in water; imaging or a quadrant photodetector records a position time series. Gieseler and colleagues describe a polystyrene bead of radius $1.03\mu \mathrm{m}$ at $293.15\mathrm{K}$, a $1064\mathrm{nm}$ trapping laser, and photodetector acquisition at $10^5\mathrm{Hz}$. Their recorded bead continues undergoing Brownian motion; higher trapping power produces a narrower equilibrium position distribution. Infer the centered variance, not a squared distance from an arbitrary camera origin:The normal distribution position histogram and agreement with dynamical calibration methods support the local harmonic optical trap interpretation. For a stipulated example, $T=300\mathrm{K}$ and ${\sigma }_x=20\mathrm{nm}$ imply $k\simeq 1.04\times 10^{-5}N/m=0.0104pN/nm$. This calculated scale is not an experimental reading.^[ref]

Position calibration is essential: if detector output $v$ has conversion $x=s(v-v_0)$, its physical variance is $s^2\mathrm{Var}(v)$. An unknown scale $s$ cannot be determined together with $k$ from that one variance. A complementary calibration demonstrated by Tolić-Nørrelykke and colleagues combines the power spectrum of thermal motion with the bead response to known sinusoidal motion of a translation stage. The driven response supplies an additional length scale, and their method does not require bead size, liquid dynamic viscosity or wall distance as inputs. It still uses a thermal dynamical model; agreement between methods is a consistency check, not a theory-free proof of equipartition theorem.^[ref]

The readout itself can hide fluctuations. Wong and Halvorsen used a $1064\mathrm{nm}$ optical trap and adjustable-exposure CCD imaging of a bead in water. They varied exposure and trapping power, held the bead away from the nearest wall, adjusted illumination to maintain comparable images, and measured tracking noise using a stuck bead. The measured variance fell as exposure increased; corrected results agreed with their time-averaging model. For a stationary overdamped harmonic optical trap, with optical trap relaxation time $\tau$, an exposure of length $W$ measures ${\bar{x}}_W=W^{-1}{\int }_0^{W}x(t)dt$. Its ideal noise-free variance isFor $\alpha >0$, $S(\alpha )<1$: using the blurred variance in the uncorrected formula overestimates $k$ by $1/S(\alpha )$. At $W=\tau$, the calculated factor is $e/2\simeq 1.36$. Additive tracking noise instead raises the measured variance and biases the inferred $k$ downward if left uncorrected.^[ref]

Useful controls are to record over many optical trap relaxation times, compare blocks for drift, calibrate the position scale, measure fixed-bead noise, vary exposure, and check the position histogram for departures from a normal distribution. A sparse sequence of accurate instantaneous positions can estimate an equilibrium variance; a long exposure averages away motion even if many frames are collected. Flow, appreciable laser heating, escape from the harmonic region or nonconservative forces can break the assumed thermal equilibrium. Inferring $k$ by imposing equipartition theorem is a calibration, whereas testing the relation requires an independently established stiffness and position scale. The beautiful practical connection is that microscopic randomness can set a force scale: once calibrated, a small constant applied force displaces the mean by $f/k$. The same quadratic-term counting predicts macroscopic heat capacity and makes microscopic fluctuations useful measurements.

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.

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 $Z_E[J]=\int \mathcal{D}\varphi e^{-S_E[\varphi ]+\int J\varphi }$, define $W_E=\log Z_E$, $\phi =\delta W_E/\delta J$ and the Legendre transform ${\Gamma }_E[\phi ]=\int J\phi -W_E[J]$. Its functional derivatives generate Euclidean one-particle irreducible vertices. A bosonic Gaussian integral about a background gives the one-loop term $\frac{1}{2}\mathrm{Tr}\log S_E^{(2)}[\phi ]$. This is the Euclidean version of the quantum effective action, with signs fixed by $e^{-S_E}$ rather than $e^{iS}$.

For $m^2>0$ and $n>d/2$, Schwinger parameterization and a Gaussian integral give $I_{d,n}(m)=\int d^{d}k(2\pi {)}^{-d}(k^2+m^2{)}^{-n}=\Gamma (n-d/2)(m^2{)}^{d/2-n}/[(4\pi {)}^{d/2}\Gamma (n)]$. Its meromorphic continuation supplies dimensional regularization poles.

For $Z$ with a standard normal distribution and $a\ge 0$, split the Gaussian integral at zero and complete the square to obtain $\mathbb{E}{e}^{a|Z|}=2e^{a^2/2}\Phi (a)\le 2e^{a^2/2}$, where $\Phi$ is the standard normal distribution function. This converts a bound on an absolute Gaussian random variable into a two-sided Gaussian tail bound using Markov inequality.

For independent observations $y_s\sim N(\beta +d_s,V_s)$, put $a_s=V_s^{-1}$, $S={\sum }_s{a}_s$, $r_s=y_s-d_s$, $\bar{r}=S^{-1}{\sum }_s{a}_s{r}_s$ and $Q={\sum }_s{a}_s(r_s-\bar{r}{)}^2$. Integrating the likelihood function against a flat improper prior on $\beta$ gives, up to the arbitrary prior constant,Indeed ${\sum }_s{a}_s(r_s-\beta {)}^2=Q+S(\beta -\bar{r}{)}^2$, and the remaining one-dimensional Gaussian integral is $\sqrt{2\pi /S}$. The conditional Bayesian posterior of $\beta$ is $N(\bar{r},S^{-1})$.

Use an orthonormal real Fourier series basis on a circle of length $\beta$, with ${\nu }_r=2\pi r/\beta$. Retain the constant coefficient and the sine/cosine pair for each $1\le r\le N$. The oscillator action becomes $\frac{1}{2}{\omega }^2{q}_0^2+\frac{1}{2}{\sum }_r({\omega }^2+{\nu }_r^2)(q_r^2+s_r^2)$. A product of real Gaussian integrals gives the displayed functional determinant, with $A_N$ independent of $\omega$. Ratios for two frequencies converge regardless of $A_N$. The absolute limit requires its normalization: $A_N={\beta }^{-1}{\prod }_r{\nu }_r^2$ gives the thermal partition function of a quantum harmonic oscillator in the limit. Multiplying this choice by an oscillating positive sequence preserves all frequency ratios but destroys the absolute limit.

Let the initial Gaussian wave packet beFourier transforming, multiplying each momentum amplitude by $e^{-i{p}^{2}t/(2m\hslash )}$ and evaluating the Gaussian integral givesIts probability density is a normalized Gaussian with center $x_0+p_{0}t/m$ and width ${\sigma }_x={\sigma }_0\sqrt{1+{\eta }^2}$; the momentum width remains ${\sigma }_p=\hslash /(2{\sigma }_0)$. Its symmetrized position-momentum covariance is $C_{xp}=t{\sigma }_p^2/m$, soThus free spreading preserves equality in the Robertson-Schrödinger uncertainty relation while producing a nonzero quadratic position-dependent phase and a product ${\sigma }_x{\sigma }_p$ larger than $\hslash /2$. More generally, ${\sigma }_x(t{)}^2={\sigma }_x(0{)}^2+2t{C}_{xp}(0)/m+t^2{\sigma }_p(0{)}^2/m^2$, so an initially correlated Gaussian can first contract.

The Gamma function extends factorials to a continuous argument. Its defining Gamma integral iswhere $t^{z-1}=\exp ((z-1)\log t)$ uses the real natural logarithm for $t>0$. The Gamma function recurrence gives its meromorphic continuation to the complex plane, with simple poles at $0,-1,-2,\dots$ and no zeros. See Gamma function residue at a nonpositive integer and Gamma function has no zeros for proof outlines.^[ref]

For nonnegative integers $n$, the factorial identity has a shift of one:Thus $\Gamma (n)=(n-1)!$ for positive integers $n$. Half-integer arguments supply a different familiar constant:The first value follows from the Gaussian integral; the second follows from the Gamma function recurrence.^[ref]

The plot shows only positive real arguments; the gamma function also accepts complex numbers away from its poles. Matching the factorials alone does not select a unique interpolation, even if the Gamma function recurrence is imposed: $f(x)=\Gamma (x)\exp (\sin (2\pi x))$ has the same positive-integer values and the same recurrence relation. The Bohr-Mollerup theorem supplies the extra condition that selects the gamma function.

The useful idea is to replace a discrete product of numbers by an integral whose exponent varies continuously. The Gamma integral then evaluates whole families of integrals at once, provides the normalizing constant for the Gamma distribution in probability theory, and appears in physics, for example in the Bose integral. Its large-positive-argument growth reconnects this continuous description to factorials through Stirling's formula.^[ref]

Historically, Leonhard Euler studied the interpolation of factorials in his paper E19, written in 1729 and published in 1738. His integral over $(0,1)$ becomes the Gamma integral after a change of variables, with the same one-step shift in the argument.^[ref]

For real $a>0$ and $\mathrm{Re}s>0$, the change of variables $u=at$ in the defining improper integral for the Gamma function givesThe improper integral converges absolutely: near zero its absolute integrand behaves like $t^{\mathrm{Re}s-1}$, and at infinity the exponential function overcomes every fixed power.^[ref]

For real $s=1,2,4$ and $a=1$, the integrands are $e^{-t}$, $t{e}^{-t}$, and $t^3{e}^{-t}$, with total areas $\Gamma (1)=1$, $\Gamma (2)=1$, and $\Gamma (4)=6$. The third integrand is taller and has six times the area of the first, although the exponential function makes each tail decay.

For real $s>0$, the kernel $t^{s-1}{e}^{-t}$ becomes the probability density function of a Gamma distribution with shape $s$ and unit rate only after division by its normalizing constant $\Gamma (s)$. In particular, the shape-four probability density function is $t^3{e}^{-t}/6$. More generally, the Gamma integral verifies that $a^s{t}^{s-1}{e}^{-at}/\Gamma (s)$ integrates to one for shape $s>0$ and rate $a>0$.

For the half-integer evaluation, the change of variables $t=u^2$ reduces the Gamma integral to the Gaussian integral:Consequently, for real $u>0$,

For a slowly varying nuclear factor and penetration factor $\exp [-2\sqrt{E_B/E}]$, the thermal integral is ${\int }_0^{\infty }{e}^{-\Phi (E)}dE$, where $\Phi (E)=E/(k_{B}T)+2\sqrt{E_B/E}$. The Gamow peak occurs at $E_0=[E_B(k_{B}T{)}^2{]}^{1/3}$, with $\Phi (E_0)=3E_0/(k_{B}T)$ and ${\Phi }^{\prime \prime }(E_0)=3/(2k_{B}T{E}_0)$. The Gaussian integral gives $\sqrt{4\pi k_{B}T{E}_0/3}{e}^{-3E_0/(k_{B}T)}$. Combined with the Maxwellian prefactor $(k_{B}T{)}^{-3/2}$ this produces the displayed thermonuclear reaction rate at fixed reactant density. The approximation assumes $E_B\gg k_{B}T$ and excludes strong resonances or rapidly varying nuclear factors.

For a stationary autoregressive process of order one with unit-variance normal distribution marginals and fixed $0<\rho <1$, the Gaussian kernel covariance correction at lag $j$ has the displayed form. The difference of two observations has variance $2(1-{\rho }^j)$, and their characteristic functions reduce the Fourier integral to the difference of two Gaussian integrals. The mean value theorem gives $0\le g_h(j)\le \sqrt{\pi }{\rho }^j/[2(1-\rho {)}^{3/2}]$, uniformly over $h>0$. Thus the integrated variance differs from its independent-sample counterpart by $O(n^{-1})$. The marginal bias is unchanged, so this dependence does not alter the leading optimal bandwidth or the $n^{-4/5}$ mean integrated squared error for a second-order density kernel.