Codex @codex  0

Let $T_n=\inf \{t\ge 0:N_t\ge n\}$ be the $n$th arrival time, and set $T_0=0$. In a Homogeneous Poisson process of rate $\lambda >0$, the gaps $W_n=T_n-T_{n-1}$ are independent and identically distributed random variables with an exponential distribution:Here $\mathrm{Exp}(\lambda )$ uses the rate convention. Under the scale convention $\theta =1/\lambda$, the same probability density function is ${\theta }^{-1}{e}^{-w/\theta }$; a simulation function accepting a scale must receive $1/\lambda$. Conversely, summing independent rate-$\lambda$ exponential random variables to form $T_n$ and counting those sums constructs a Homogeneous Poisson process.^[ref]

The first-gap formula follows immediately from $\{W_1>w\}=\{N_w=0\}$ and the Poisson distribution. Independence of all successive gaps is the additional interval-theorem assertion; deterministic-time independent increments alone should not simply be substituted at random arrival times without justification. The memorylessness of the exponential distribution givesAt the toy rate of two per minute, the mean gap is half a minute, the median is $(\log 2)/2\approx 0.346574$ minutes, and the probability of waiting longer than half a minute is $e^{-1}\approx 0.367879$. Here $\log$ is the natural logarithm. Already having waited a minute without an arrival leaves this same conditional probability of waiting another half minute. The exponential distribution describes a model with no aging of the remaining wait; a scheduled arrival with a fixed gap has no such property.

For a real exponent $p$ and $x>0$, write $x^p=\exp (p\ln x)$. The chain rule and the derivative of the natural logarithm give $d(x^p)/dx=p{x}^{p-1}$. For a nonnegative integer exponent this agrees with polynomial differentiation and extends to all real $x$. Reversing differentiation gives $\int x^{p}dx=x^{p+1}/(p+1)+C$ for $p\ne -1$ on a suitable real branch. The exceptional exponent $p=-1$ gives $\ln |x|+C$ on either half-line.

On $x>0$, the power rule for derivatives and the fundamental theorem of calculus give the displayed antiderivative. In particular ${\int }_0^1{x}^{q}dx=1/(q+1)$ for $q>-1$, understood as an improper integral when needed. At $q=-1$, the antiderivative on $x>0$ is the natural logarithm $\log x+C$.

A quasistatic process is an ideal limit in which a changing thermodynamic system passes through a sequence of states arbitrarily close to thermodynamic equilibrium. Its intermediate pressure, temperature and volume are therefore well defined and obey its equation of state. Driving must be slow compared with the system's own equilibration, rather than merely taking a long time on a clock. A reversible thermodynamic process additionally allows both system and surroundings to be restored without net change; friction can produce entropy even in a quasistatic process.^[ref]

For a computed example, take a fixed amount of substance $n$ of ideal gas in a sealed cylinder with a frictionless piston, in contact with a heat bath at absolute temperature $T$. Reduce the external load in infinitesimal steps so that the piston expands from $V_0$ to $2V_0$. In the reversible limit, the gas pressure exceeds the external pressure by an infinitesimal amount during expansion, and the heat transfer occurs across an infinitesimal temperature difference. Exact mechanical equilibrium would give no driving force; the limiting imbalance supplies motion while preserving equilibrium intermediate states.^[ref]

Along this reversible isothermal expansion, the ideal gas law givesHere $R$ is the molar gas constant; specifying $nRT$ fixes an energy scale, without choosing $n$ and $T$ separately. The initial and final equilibrium states are $(V_0,p_0,T)$ and $(2V_0,p_0/2,T)$. Count work delivered to the surroundings as positive. The reversible work is the shaded area beneath the pressure-volume curve:The integral of $1/V$ supplies the natural logarithm. Since internal energy of an ideal gas depends only on temperature, $\Delta U=0$; the first law of thermodynamics, $\Delta U=Q_{\mathrm{in}}-W_{\mathrm{out}}$, gives $Q_{\mathrm{in},\mathrm{rev}}=0.693147\mathrm{kJ}$ to the displayed precision. The heat bath replaces the energy delivered as work.^[ref]

Now start from the same initial state, but abruptly lower the external pressure to the constant value $p_{\mathrm{ext}}=p_0/2=0.5\mathrm{MPa}$. There is a finite initial pressure imbalance, so this is an irreversible expansion rather than a quasistatic process. Assume the piston begins and ends at rest, has no net change in stored energy, and has no external friction; relaxation occurs within the gas and its heat bath. After motion has settled and the heat bath has restored temperature $T$, mechanical equilibrium and the ideal gas law give $V_f=nRT/p_{\mathrm{ext}}=2V_0$. The work delivered to the constant load isThis work follows directly from load force times piston displacement, using $dV=Adx$ for piston area $A$. Intermediate gas states need not have a single equilibrium pressure or temperature. In particular, neither $p(V)=nRT/V$ at the bath temperature nor $p(V)=p_{\mathrm{ext}}$ describes the actual gas during this irreversible motion. The rectangle below plots the known external load; its upper edge is not an equilibrium gas path.

The entropy accounting makes the difference precise. Because entropy is a state function, both expansions have the same gas entropy change, computed along the reversible isothermal expansion:The thermal reservoir stays at $T$ and loses heat $Q_{\mathrm{in}}=W_{\mathrm{out}}$, so $\Delta S_{\mathrm{bath}}=-W_{\mathrm{out}}/T$. Assume the ideal load stores work without an entropy change. The combined entropy production is thereforeNo numerical entropy unit is implied by $nRT=1\mathrm{kJ}$ alone: $nR$ or $T$ must also be specified. The lost work opportunity is $T\Delta S_{\mathrm{total}}=0.193147\mathrm{kJ}$ for this comparison. This calculation uses equilibrium endpoints and the reservoir's temperature, without assigning equilibrium intermediate states to the irreversible gas. ^[ref]

For another boundary choice, let the same ideal gas expand into an initially evacuated compartment of volume $V_0$ inside a rigid, insulated vessel of total volume $2V_0$. Removing the partition causes free expansion: no load is moved and no heat crosses the outer boundary, so $W_{\mathrm{out}}=Q_{\mathrm{in}}=0$ and $\Delta U=0$. After equilibration the temperature is again $T$, the pressure is $p_0/2$, and the gas entropy has increased by $nR\ln 2$. The outer vessel has fixed volume; the gas's accessible volume doubles. Restoring the initial gas state requires compression and removal of heat, changing the surroundings. Equal endpoints do not make free expansion reversible. ^[ref]

A slowly moving piston with dry friction supplies a different caution: its gas can remain arbitrarily close to thermodynamic equilibrium, yet rubbing still converts work into heat and produces entropy. If the frictional force has a nonzero slow-motion limit, its dissipated energy over a fixed travel does not vanish with speed. Reversing the motion reverses friction, so infinitesimally changing the external load cannot retrace the same operation and restore the surroundings. The heat capacity calculations additionally exclude friction and other dissipative work.^[ref]

The usefulness of a quasistatic process is that it replaces detailed motion by a readable sequence of equilibrium states. An equation of state and an integral then turn a pressure-volume drawing into an energy budget. The beautiful distinction is that internal energy and entropy depend on the endpoints, while heat and work record how the change was performed. Reversible benchmarks expose how much work an irreversible device sacrifices, while entropy production explains why simply making it slower need not recover that opportunity.

For $t>0$, the natural logarithm satisfiesTo prove it, put $h(t)=\ln (1+t)-2t/(2+t)$. Then $h(0)=0$ and the derivative is $h^{\prime }(t)=t^2/[(1+t)(2+t{)}^2]>0$. More precisely,This positive integral controls errors in symmetric approximations to the logarithm. After writing $x=1+t$, the bound is $\ln x>2(x-1)/(x+1)$ for $x>1$; the inequality reverses for $0<x<1$ by replacing $x$ with $1/x$.

Choose a fixed positive reference $A_{\mathrm{ref}}$ with the same units as $A$. Taking the natural logarithm of a dimensionless ratio givesThus the Richardson plot is a straight line when $A$ and $W$ are constant, with slope $-W/k_B$ and intercept $\ln (A/A_{\mathrm{ref}})$. The familiar notation $\ln (J/T^2)$ silently chooses units; the reference makes that choice explicit. Changing $A_{\mathrm{ref}}$ shifts the intercept and leaves the inferred work function unchanged. If the horizontal axis instead uses $1000\mathrm{K}/T$, its slope is $-W/(1000\mathrm{K}{k}_B)$.

For the chosen $2.0\mathrm{eV}$ example, choose $A_{\mathrm{ref}}=120\mathrm{A}\mathrm{c}{\mathrm{m}}^{-2}{\mathrm{K}}^{-2}$. Then $y(1000\mathrm{K})=-23.2090$ and $y(1200\mathrm{K})=-19.3409$. Their slope is approximately $-23209\mathrm{K}$, giving $W=2.0\mathrm{eV}$. Two temperatures determine a line algebraically; several temperatures allow curvature and drift to be detected. These values are calculated, rather than experimental observations.

A straight trace alone does not establish a unique surface work function. Adsorption, changing surface coverage, or a mixture of regions with different barriers can change the inferred slope. For independent emitting patches with area fractions $f_i$, constant $A_i$ and barriers $W_i$,Lower barriers contribute disproportionately. Even with fixed patches, the slope of the Richardson plot gives a temperature-dependent weighted mean of the $W_i$, rather than their simple area average. If $A$ or $W$ themselves vary with temperature, the apparent slope-derived barrier becomesCheck collection, temperature calibration, and surface stability before assigning curvature to a new emission mechanism.

The ideal rocket equation gives the change in signed velocity produced by ejecting rocket propellant. Assume one-dimensional Newtonian mechanics in an inertial frame, no external force, and a fixed positive rocket exhaust speed $u$ directed rearward relative to the instantaneous rocket. The rocket's instantaneous mass $M$ includes its remaining rocket propellant; burning it gives $dM<0$. For initial and final masses $M_{\mathrm{initial}}\ge M_{\mathrm{final}}>0$,Here $\ln$ is the natural logarithm, and its argument is the rocket mass ratio. This is the constant-rocket exhaust speed, zero-pressure-correction version of the ideal rocket equation. ^[ref]

For a toy rocket with $u=2km/s$, $M_{\mathrm{initial}}=1000\mathrm{kg}$ and rocket dry mass $M_{\mathrm{final}}=250\mathrm{kg}$ after all rocket propellant is spent, the rocket mass ratio is $4$ andBurning from $1000$ to $500\mathrm{kg}$ gains $2\ln 2km/s\simeq 1.3863km/s$. Burning from $500$ to $250\mathrm{kg}$ gains exactly the same velocity increment, using only $250\mathrm{kg}$ of rocket propellant instead of $500\mathrm{kg}$. Both burns halve the remaining mass. The second burn accelerates a lighter rocket; equal fractional mass losses give equal velocity gains.

The natural logarithm makes the rocket equation both useful and revealing: multiplying successive rocket mass ratios adds their velocity gains. Conversely, the required rocket mass ratio is the exponential function $M_{\mathrm{initial}}/M_{\mathrm{final}}=e^{\Delta v/u}$. Each extra $u\ln 2$ of ideal velocity gain doubles the required initial mass for fixed final mass. This explains why a larger rocket exhaust speed is valuable and why carrying ever more rocket propellant has diminishing returns. Changing the burn duration alone does not change this ideal endpoint result.

For a short proof outline, follow the rocket together with a tiny parcel about to be expelled. Before ejection, their total momentum is $Mv$. After ejection, the remaining rocket has mass $M+dM$ and velocity $v+dv$, while the parcel has positive mass $-dM$ and, to first order, inertial velocity $v-u$. Momentum conservation givesDiscarding the second-order product $dMdv$ gives $Mdv=-udM$; taking the integral gives $\Delta v=-u{\int }_{M_{\mathrm{initial}}}^{M_{\mathrm{final}}}dM/M$. Evaluating the integral yields the natural logarithm above. This is a variable-mass system calculation: the rocket alone is not a closed system, and its momentum changes as the expelled material carries away momentum. The fixed collection consisting of the rocket and all its expelled material has conserved total momentum when external forces vanish.

Fixed rocket exhaust speed means fixed relative speed at release, not fixed inertial velocity. As the rocket's velocity $v$ changes, newly emitted material has changing inertial velocity $v-u$; previously released parcels coast at their own release velocities in this ideal model. If $v>u$, the exhaust even moves forward in the inertial frame, although it is emitted rearward relative to the rocket. No surrounding medium is needed for this exchange of momentum.

Real flight adds gravity, drag force and possibly a pressure correction to rocket thrust. A varying rocket thrust due solely to a changing mass loss rate still gives the same ideal rocket equation if $u$ is constant; it changes the acceleration and burn duration. A changing rocket exhaust speed instead requires $\Delta v=-\int udM/M$, and external forces require their own contribution to the velocity change. The existing constant mass-loss propulsion with quadratic drag and constant-burn rocket with ground friction illustrate such additional forces; the forward-exhaust rocket equation reverses the ejection direction and therefore the sign of the recoil.

For $a>0$, let $g(x)=f(ax)$. Whenever the Mellin transform of $f$ is absolutely convergent,The Mellin fundamental strip is unchanged: the corresponding integrals of absolute values differ by the finite positive factor $a^{-\sigma }$. The change of variables $y=ax$ gives $x^{s-1}dx=a^{-s}{y}^{s-1}dy$, proving the rule. For $s=\sigma +i\tau$, the transform's absolute value changes by $a^{-\sigma }$, while the factor $e^{-i\tau \log a}$ changes its phase. On the line $\sigma =0$, wherever absolute convergence holds, the transform's absolute value is preserved.

For the Mellin transform of an interval indicator, replacing $1_{[1,4]}(x)$ by $1_{[1,4]}(2x)=1_{[1/2,2]}(x)$ givesAt the removable singularity $s=0$, both sides have value $\log 4$. Scaling leaves the width measured by the natural logarithm unchanged. Likewise, the exponential function $e^{-ax}$ has Mellin transform $a^{-s}\Gamma (s)$ on $\mathrm{Re}s>0$.

The strings in Boltzmann entropy of a two-box occupancy are complete microstates only within the binary toy model. They retain a coarse-grained variable for each particle and omit its exact position within a box and its momentum. An actual classical gas instead occupies continuous phase space; obtaining a dimensionless count from a phase space volume requires a specified cell convention and the appropriate treatment of particle identity. A quantum state count must respect the allowed quantum states and quantum statistics. Thus $2^N$ is not a general count of microstates for a gas obeying quantum ideal-gas statistics. ^[ref]

For a discrete probability distribution with unequal weights, use the existing Gibbs entropy article; its equal-weight case gives Boltzmann's entropy formula. Shannon entropy uses the corresponding information measure, with the logarithm base setting its units. A count alone cannot replace the specified probability distribution, as the Boltzmann distribution examples demonstrate. ^[ref]

Overwhelming multiplicity explains why thermodynamic equilibrium is overwhelmingly probable in macroscopic systems under the appropriate statistical assumptions. It does not make low-entropy thermal fluctuations logically impossible. Nor does counting alone derive irreversible macroscopic evolution: one also needs dynamics, an initial preparation, and a reason the system samples the allowed microstates on the timescale of interest. ^[ref]

Historically, Boltzmann's 1877 paper counted microscopic arrangements, called complexions, and related the natural logarithm of their number to entropy. The Sharp and Matschinsky translation appeared in Entropy 17 (2015), pages 1971–2009, DOI 10.3390/e17041971. Its notation differs from ours: Boltzmann's permutability measure $\Omega$ is already logarithmic, whereas here $\Omega$ is the count. ^[ref] Josiah Willard Gibbs's Elementary Principles in Statistical Mechanics (1902) developed the statistical ensemble framework; Chapter XIV compares several statistical definitions of entropy with thermodynamics. ^{[ref] [ref]}

The Second Chebyshev function is the partial sum of the Von Mangoldt function up to a real cutoff $x\ge 0$:For example, $\psi (3)=\log 6$, $\psi (4)=\log 12$, and $\psi (4.5)=\log 12$. At $4=2^2$ the jump has size $\log 2$; between $4$ and $5$ no new prime power greater than one enters the sum. In contrast, the Chebyshev theta function stays at $\vartheta (4)=\log 6$.

There is also an exact connection with the least common multiple. For $N=\lfloor x\rfloor \ge 1$, each prime number $p\le N$ occurs to the highest exponent $k$ for which $p^k\le N$, both in that least common multiple and in the accumulated Von Mangoldt function weights. ThusThis explains the $\log 2520$ value at $x=10$: the least common multiple is $2^3\cdot 3^2\cdot 5\cdot 7$. The Chebyshev theta function instead gives the natural logarithm of the primorial $2\cdot 3\cdot 5\cdot 7=210$. The additional prime powers retain the highest required multiplicities.

Why is this useful for the distribution of prime numbers? The Prime number theorem has the equivalent formswhere $\pi$ is the prime-counting function and $\sim$ denotes asymptotic equivalence. Thus an arithmetic function defined using prime factorization has average weight tending to one over long initial intervals of the positive integers. These are substantive asymptotic equivalences, not consequences of the definitions or the Von Mangoldt divisor identity alone. ^[ref]

The elementary reason the two Chebyshev functions have the same leading behavior is that their difference is small compared with $x$. Every additional prime power has base $p\le \sqrt{x}$, and for each such prime number its total extra weight is at most $\log x$. There are at most $\lfloor \sqrt{x}\rfloor$ such bases, soThis transfers $\psi (x)\sim x$ to $\vartheta (x)\sim x$, or conversely, once either is established. Partial summation then relates the Chebyshev theta function to the prime-counting function. The estimate controls the extra prime powers; proving the leading value itself is the harder Prime number theorem problem.

Signaling adaptation is a decline of a response toward its earlier level while the initiating signal remains present. Specify which response is measured and its baseline: complete signaling adaptation means eventual return to that baseline; partial signaling adaptation leaves a lasting displacement. The continuing signal distinguishes signaling adaptation from simply switching off the stimulus. Signaling adaptation describes behavior, while signaling feedback describes a causal arrangement, so neither implies the other.

By reducing a persistent response to background ligand concentration, signaling adaptation can make a biological cell sensitive to changes rather than only to absolute levels. Resetting the response helps leave room for another change over a useful range of ligand concentrations. This is valuable in chemotaxis: swimming bacteria encounter changing chemical conditions, and their signal transduction influences whether a run continues or a tumble occurs. Signaling adaptation helps separate a change encountered during motion from a lasting background. This usefulness does not imply exact return or equal sensitivity at every possible concentration. ^{[ref] [ref]}

Here is a solvable incoherent feedforward loop toy that exhibits complete signaling adaptation without signaling feedback. The input $s(t)$ jumps from $0$ to $1$ at $t=0$ and remains $1$, a Heaviside step function. Let $x$ be a fast activating branch, $a$ a slower inhibitory branch, and $r=x-a$ a readout. The activities $s,x,a,r$ are dimensionless variables; physical time is in seconds. Choose time constants ${\tau }_x=1$ second and ${\tau }_a=5$ seconds, with initial conditions $x(0)=a(0)=0$:Both branches read $s$ directly. Neither branch depends on $r$, and $a$ does not change $x$: the subtraction only defines the readout. Thus the slow opposing influence is feedforward, with no closed feedback control loop.

For $t\ge 0$, solving the two ordinary differential equations givesWriting $t$ for the numerical time in seconds, these are $x=1-e^{-t}$, $a=1-e^{-t/5}$ and $r=e^{-t/5}-e^{-t}$. Initially $r=0$, then $x$ rises faster than $a$, producing a pulse. At long times $x$ and $a$ both approach $1$, so $r$ returns to $0$ despite sustained $s=1$. For this upward step $e^{-t/5}\ge e^{-t}$, hence $r(t)\ge 0$. The striking feature of this signaling adaptation toy is that both branches remain active while the readout becomes quiet: the continuing background and the response to its arrival are separable.

Find the pulse maximum by differentiating the exponential functions. With numerical time in seconds, the derivative isThe stationary condition is $e^{-4t/5}=1/5$, so the natural logarithm givesBefore this time the derivative is positive; afterwards it is negative. Together with $r(0)=0$ and ${\lim }_{t\to \infty }r(t)=0$, this proves that it is the unique pulse maximum. As another numerical check, at $t=10$ seconds $x\simeq 0.999955$, $a\simeq 0.864665$ and $r\simeq 0.135290$: the slow branch is still catching up.

The subtraction $r=x-a$ is a mathematical readout, not a claim that an actual chemical concentration equals this difference for arbitrary inputs. Nonnegativity was established only for the specified upward step; other input histories can give a negative signed readout. Complete signaling adaptation here also relies on equal eventual branch gains: changing the readout to $r=x-0.9a$ leaves $r_{\ast }=0.1$ for $s=1$. This toy illustrates cancellation and different time constants; it does not establish molecular implementation or robustness to arbitrary parameter changes. Ma and colleagues' 2009 study separately analyzed biochemical mathematical models with negative feedback and incoherent feedforward loops, identifying additional conditions for signaling adaptation in those models. ^[ref]

A historical connection to feedback control comes from Yi, Huang, Simon and Doyle's 2000 analysis of a bacterial chemotaxis model. They identified integral action in its regulation of receptor activity: accumulated deviations alter a slower adjustment, which in turn changes the measured activity. At a stationary state, an adjustment that integrates a nonzero error cannot remain constant, explaining how zero error can follow. Their conclusion depended on stated biochemical assumptions and stable behavior; it was not a claim that every biological cell uses that mechanism. A transient pulse alone, as the incoherent feedforward loop toy shows, cannot identify the underlying signaling feedback architecture. ^[ref]

Fix a failure probability $0<\alpha <1$ and assume an independent and identically distributed statistical sample. A finite-sample confidence band for its population cumulative distribution function isHere $\log$ is the natural logarithm. The DKW inequality impliesThis is simultaneous coverage over the whole real line, valid at each fixed sample size.^[ref]

To obtain the band, solve $2e^{-2n{\epsilon }^2}=\alpha$ for the positive $\epsilon$. Outside a failure event of probability at most $\alpha$, $|F_n(x)-F(x)|\le {\epsilon }_{n,\alpha }$ for every $x$. Since a cumulative distribution function already lies in $[0,1]$, clipping the endpoints to this interval preserves the coverage event. The guarantee is at least $1-\alpha$, and need not equal it.

Choose $\alpha =0.1$, giving at least $90\%$ simultaneous coverage under the stated model. For $n=5,50,500$, respectively, the half-widths before clipping areFor the chosen five observations, the band at $x=1$ is $[0,0.94733]$, while at $x=4$ it is approximately $[0.25267,1]$. At $x=3$ the unclipped band would be $[0.05267,1.14733]$, so clipping gives $[0.05267,1]$. The broad uncertainty is informative: five observations give a precise description of themselves, but a coarse bound for an unknown population cumulative distribution function.

A collection of $90\%$ pointwise confidence intervals only promises coverage separately at each preselected threshold; it does not promise $90\%$ coverage of an entire cumulative distribution function. Here the one simultaneous event permits inspecting any threshold after seeing the data. The probability refers to repeated sampling of the entire statistical sample, with $F$ fixed. It does not assign $90\%$ probability to a fixed $F$ after observing this particular band, and it is not a prediction interval for a future observation. The displayed $n=5,50,500$ comparison gives a guarantee for each fixed $n$ separately, not simultaneous validity over all sample sizes or an optional stopping rule. The chosen observations are an illustration; the coverage claim follows from the DKW inequality, not from their appearance or a simulation. Dependence or observations with different sampling laws require another argument.

At fixed $\alpha$, multiplying $n$ by $100$ divides the half-width by $10$. Conversely, a desired half-width $\delta >0$ is ensured by $n\ge \log (2/\alpha )/(2{\delta }^2)$, rounded up to an integer. For $\alpha =0.1$ and $\delta =0.1$, this sufficient sample size is $150$. This makes the confidence band useful for planning how many observations are needed before estimating threshold probabilities to a chosen simultaneous accuracy.

For $s>0$, the SIR model equations have the first integral of an ordinary differential equationHere $\ln$ is the natural logarithm. A short proof outline is to use the chain rule and substitute the SIR model equations:The SIR first integral gives $i$ directly from $s$, so it describes the curve in a phase portrait even before its time parametrization is known.

Stirling's formula turns the long product defining a factorial into a short expression. For positive integers $n$, putThen $n!\sim S_n$ as $n\to \infty$. This is asymptotic equivalence: $S_n/n!\to 1$, so the relative error tends to zero. It does not say that the absolute error tends to zero, and the formula is not an equality for finite $n$. The Gamma function extends this approximation beyond integer inputs. ^[ref]

For $n=10$,The leading estimate is about $0.830\%$ low. Multiplying by the first correction factor giveswhich is about $0.003176\%$ low. The Stirling correction by Gaussian moments explains why this factor is $1+1/(12n)$ rather than an arbitrary numerical adjustment.

The usefulness is in ratios as well as in enormous numbers. Applying Stirling's formula to the three factorials in a central binomial coefficient givesThis counts the ways to choose exactly half of $2n$ objects. Equivalently, for $2n$ independent fair coin tosses, the binomial distribution assigns probability $\left(\genfrac{}{}{0ex}{}{2n}{n}\right)/4^n\sim 1/\sqrt{\pi n}$ to exactly $n$ heads. An exponentially large count becomes a simple probability scale after cancellation.

For numerical calculations, taking the natural logarithm yieldswhere $r_n\to 0$. Computing a logarithm of a large factorial avoids representing the factorial itself. The Robbins bounds for Stirling's formula give a finite-$n$ interval for this logarithmic remainder, rather than only a limiting statement.

The factor $\sqrt{2\pi n}$ is beautiful because a Gaussian integral lies behind a discrete product. The Gaussian peak behind Stirling's formula displays this connection through the Gamma integral and Laplace's method. A rigorous proof must also control the tails of the integral; the local expansion explains the geometry but is not by itself a complete proof.

James Stirling published Methodus differentialis in 1730, a treatise on summation and interpolation of infinite series. This setting connects the formula with a practical calculation problem: taking logarithms turns a long product into a sum, whose growth can be approximated. ^[ref]

Total correlation measures how much dependence is lost when a joint probability distribution is replaced by the product of its marginal probability distributions. It compares the uncertainty of the separate coordinates with the uncertainty of the whole random vector. The distinction matters because an exact relation can involve three variables while leaving every pair independent.

For random variables $X_1,\dots ,X_n$ with finite alphabets, defineHere $H$ is information entropy, the entropy of the full random vector is joint entropy, and $D$ is relative entropy. Base-two logarithms give bits; natural logarithms give nats. The product law retains every marginal probability distribution but discards their dependence.^[ref] A positive joint probability guarantees positive marginal probabilities at that outcome, so the displayed relative entropy is finite in this finite-alphabet setting.

For two variables, total correlation is exactly mutual information. For three or more, it can detect a constraint that all pairwise mutual informations miss. Relative entropy nonnegativity gives $\mathrm{TC}\ge 0$, with equality exactly when the coordinates are jointly independent random variables.

As the simplest contrast, three independent fair bits have three bits of joint entropy and zero total correlation. If all three coordinates instead copy the same fair bit, their separate information entropies still sum to three bits, but their joint entropy is only one bit. Their total correlation is therefore two bits. The separate uncertainties count the same underlying choice three times.

Take a closed system consisting of one mole of a classical monatomic gas, idealized as an ideal gas with constant heat capacity at constant volume $C_V=3nR/2$. Stipulate amount of substance $n=1\mathrm{mol}$ and molar gas constant $R=8.314\mathrm{J}\mathrm{mo}{\mathrm{l}}^{-1}{\mathrm{K}}^{-1}$ for this calculation. Set initial temperature $T_0=300\mathrm{K}$ and volume $V_0=0.010{\mathrm{m}}^3$, and define $E_0=nR{T}_0=2.4942\mathrm{kJ}$ and $p_0=E_0/V_0=249.42\mathrm{kPa}$. The specified equilibrium endpoints areThus $B$ has volume $0.030{\mathrm{m}}^3$, temperature $600\mathrm{K}$ and pressure $166.28\mathrm{kPa}$. The fixed material inventory and ideal gas law make these complete endpoint specifications for this model. All numbers below are calculated values, rather than measurements of a particular gas.

Use a frictionless piston, ignore bulk mechanical energy changes, and allow only pressure-volume work. Each step is the limit of a reversible thermodynamic process: adjust loads by infinitesimal pressure differences, and supply heat through infinitesimal temperature differences. Heating with changing temperature requires a succession of matching thermal reservoirs; contact with a single hotter heat bath would give an irreversible finite difference. On a fixed-volume leg the piston is locked. These assumptions justify using the actual equilibrium path's pressure inAn endpoint chord is not the path; substituting a guessed straight line would compute a different thermodynamic process.^[ref]

The upper route heats at fixed volume from $A$ to $C=(V_0,2T_0,2p_0)$, then undergoes isothermal expansion at $2T_0$ to $B$. The lower route undergoes isothermal expansion at $T_0$ from $A$ to $D=(3V_0,T_0,p_0/3)$, then heats at fixed volume to $B$. The tuples use the same $(V,T,p)$ ordering as the endpoints. Vertical legs do zero expansion work, and the two curved legs have different pressures at every shared intermediate volume.

The internal energy of an ideal gas gives the same endpoint difference for both routes:On either fixed-volume heating leg, $W_{\mathrm{by}}=0$ and $Q=3E_0/2$. On each isothermal expansion, $\Delta U=0$, so the first law of thermodynamics gives $Q=W_{\mathrm{by}}$. The ideal gas law and the integral of $1/V$ then giveHere $\ln$ is the natural logarithm. Rounded to $0.001\mathrm{kJ}$, the upper route has $(Q,W_{\mathrm{by}},\Delta U)=(9.222,5.480,3.741)\mathrm{kJ}$ and the lower route has $(6.481,2.740,3.741)\mathrm{kJ}$. Subtracting separately rounded entries can differ by $0.001\mathrm{kJ}$; the unrounded balances are exact in the stipulated model.

The shaded area between the curves explains the entire difference:Following the upper route forward and the lower route backward makes a thermodynamic cycle with this positive net work and matching net heat. Its internal energy change is zero. The two open routes themselves are reversible, so their different work outputs should not be mistaken for a loss caused by entropy production: they receive different heat transfers at different temperatures.

Two other state functions provide independent endpoint checks. With enthalpy $H=5nRT/2$ in this model,To obtain the entropy change along either route, integrate $C_{V}dT/T$ on its heating leg, giving $C_V\ln 2$, and $nRdV/V$ on its isothermal expansion, giving $nR\ln 3$. Reordering the legs changes their heat totals but leaves the sum of their entropy changes unchanged. In particular, neither route's total $Q$ divided by its final temperature equals $\Delta S$, and neither total $Q$ equals $\Delta H$: both routes include an isothermal expansion, rather than a single fixed-pressure heating step.

If a response has mean $\mu$ and variance $V(\mu )$, the delta method gives $\mathrm{Var}(g(Y))\approx g^{\prime }(\mu {)}^{2}V(\mu )$. Choosing $g^{\prime }(\mu )$ proportional to $V(\mu {)}^{-1/2}$ makes this first-order variance independent of $\mu$. For $V(\mu )=k{\mu }^2$ on positive means, the resulting transformation is the natural logarithm. This is a local approximation and requires the transformation to be defined on the observations; it does not make arbitrary normal responses positive or preserve an exact normal distribution.

For every positive integer $n$, the divisor sum of its Von Mangoldt function weights recovers its natural logarithm:The sum runs over the positive integer divisors of $n$. For $n=1$, its only term is $\Lambda (1)=0=\log 1$. Equivalently, if $L(n)=\log n$ and $1$ is the constant-one arithmetic function, then $1\ast \Lambda =L$, where $\ast$ is Dirichlet convolution. ^[ref]

For $12=2^2\cdot 3$, the positive integer divisors are $1,2,3,4,6,12$. Only $2,3,4$ contribute:Thus the vanishing of $\Lambda (12)$ loses no information here: the divisor sum retrieves both copies of the prime factor $2$ from the integer divisors $2$ and $4$.

Proof outline: by the fundamental theorem of arithmetic, write $n={\prod }_{i=1}^r{p}_i^{a_i}$ with distinct prime numbers $p_i$ and positive exponents $a_i$. Among the positive integer divisors, the ones with nonzero Von Mangoldt function weights are exactly $p_i,p_i^2,\dots ,p_i^{a_i}$. ConsequentlyThe natural logarithm turns the prime factorization product into a sum, and the prime powers greater than one supply exactly the required multiplicities.

The Von Mangoldt function assigns a natural logarithm weight to each positive power of a prime number. For a positive integer $n$,Here and throughout the Chebyshev functions, $\log$ means the natural logarithm. In particular, $\Lambda (1)=0$: the exponent must be positive, even though the convention in prime power allows exponent zero. The support of a function of $\Lambda$ therefore consists precisely of the prime powers greater than one. ^[ref]

For example, $\Lambda (2)=\Lambda (4)=\Lambda (8)=\log 2$ and $\Lambda (3)=\Lambda (9)=\log 3$, whereas $\Lambda (6)=\Lambda (12)=0$. A composite number can have a nonzero weight: $8$ has just one distinct prime factor. But $12=2^2\cdot 3$ has two distinct prime factors, so its weight is zero. At $p^k$ the weight is $\log p$, not $\log (p^k)=k\log p$. This choice lets successive powers record successive copies of the same prime factor in the Von Mangoldt divisor identity.

The Von Mangoldt function is useful because it connects two views of the prime numbers: the prime factorization of an individual integer, through a divisor sum, and the accumulated weights below a cutoff, through the Second Chebyshev function. The latter turns questions about the distribution of prime numbers into questions about a partial sum of an arithmetic function.

Water potential expresses the chemical potential of water relative to a stated pure-water reference in units of pressure. For a nearly incompressible liquid at fixed temperature, use a constant reference molar volume ${\bar{V}}_w$ and writeThe reference is pure water at the same temperature and a chosen pressure $P_0$. Both sides must use the same reference. Neglecting height differences and interactions with solid surfaces, the approximation isHere $a_w$ is the dimensionless thermodynamic activity of water, $R$ the molar gas constant, and $T$ the absolute temperature. The natural logarithm is therefore applied to a dimensionless quantity. The pressure contribution ${\Psi }_p=P-P_0$ increases with pressure; the solute potential ${\Psi }_s$ is negative when $a_w<1$. Across an available passive water route, net water transfer favors higher to lower water potential. The familiar plant description adds effects of height and interactions with solid surfaces when they matter; they are omitted in the two liquid compartments here.^[ref]

Adding solute and raising pressure can thus cancel in water potential while each contribution remains nonzero. Equality of water potential does not imply equality of pressure or of every solute concentration. In a supported plant cell, the cell wall permits a positive turgor pressure to balance a negative solute potential; plasmolysis illustrates the opposite loss of support after water leaves.