---
title: "The world inside a normalizer"
author: "Miraia S. Chiou"
creator: "Miraia Chiou"
published: "2026-09-18"
identifier: "epics-bsr-007"
source: "https://epistemicphysics.com/articles/the-world-inside-a-normalizer/"
citation: "https://epistemicphysics.com/cite/?work=epics-bsr-007"
---

# The world inside a normalizer

The constant beneath a Beta density connects predictions, counting, thermal fluctuations, and motion. Follow it far enough and a missing clock becomes visible.

[Read the illustrated article](https://epistemicphysics.com/articles/the-world-inside-a-normalizer/)

A probability density usually arrives wearing a small, forgettable denominator. The numerator gives the curve its interesting shape. The denominator makes its area equal one. It is tempting to divide, tidy the expression, and move on.

Keep that denominator on the desk. It can predict the next observation, count possible histories, compare models, and generate a geometry of statistical distinguishability. After a change of coordinates, it becomes a partition function, the object a statistical physicist uses to organize a whole thermal system.

There is another surprise waiting at the end. Even a complete probability distribution over velocities leaves a journey unfinished. To find where a traveler ends up, we also need to know how long each velocity lasts.

## Three objects that travel together

Assume an unknown binary rate $P$ has a proper Beta distribution, with positive shapes $\alpha,\beta$. Recenter it to $U=2P-1$, so possible values $u$ lie between $-1$ and $1$. Throughout this article, write the prior-inclusive evidence coordinates without their usual epistemic subscripts:

$$
\tau=\alpha+\beta,\qquad x=\alpha-\beta,\qquad
v=\mathbb E[U]=\frac{x}{\tau},\qquad m=2\sqrt{\alpha\beta}.
$$

The random variable $U$ describes the whole cloud of possible velocities. Its mean $v$ describes one feature of that cloud. The exact identity $\tau^2-x^2=m^2$ puts these parameters inside a Lorentz cone.

Our traveling kit contains three objects:

$$
q(u)=(1+u)^{\alpha-1}(1-u)^{\beta-1},\qquad
Z=\int_{-1}^{1}q(u)\,du,\qquad f_U(u)=\frac{q(u)}Z.
$$

Here $q$ is an **unnormalized kernel**, $Z$ its total weight, and $f_U$ the probability density. Think of $q$ as a recipe’s relative ingredients, $Z$ as the recipe’s total, and $f_U$ as the proportions.

Substitute $p=(1+u)/2$ into the integral. Because $du=2\,dp$, the powers of two collect into

$$
Z=2^{\tau-1}\mathrm B(\alpha,\beta)
=2^{\tau-1}\frac{\Gamma(\alpha)\Gamma(\beta)}{\Gamma(\tau)}.
$$

The Beta function $\mathrm B$ is an integral; $\Gamma$ extends the factorial, with $\Gamma(k+1)=k!$ for nonnegative integers. [Euler’s Beta integral](https://dlmf.nist.gov/5.12.E1) supplies the identity used here.

The change of variables also explains a common puzzle. Stretching the probability interval from $[0,1]$ to $[-1,1]$ halves the density’s height:

$$
f_U(u)=\frac12 f_P\!\left(\frac{1+u}{2}\right),\qquad
\int_{-1}^{1}f_U(u)\,du=1.
$$

This factor is the **Jacobian**, the correction for a coordinate’s changing length scale. Omitting it doubles the area. A coordinate change preserves total probability when the density and measure travel together.

| Shapes             | Kernel $q(u)$    | Area $Z$ |
| ------------------ | ---------------- | -------- |
| $\alpha=\beta=1/2$ | $1/\sqrt{1-u^2}$ | $\pi$    |
| $\alpha=\beta=1$   | $1$              | $2$      |
| $\alpha=\beta=2$   | $1-u^2$          | $4/3$    |

The first row is the Lorentz-factor shape. Its normalized version is the Jeffreys prior $j(u)=1/[\pi\sqrt{1-u^2}]$. Thus $\gamma(u)=\pi j(u)$: the same shape, with explicitly different areas.

Multiplying $q$ by any positive constant multiplies $Z$ by that constant and leaves $f_U$ unchanged. A raw area acquires meaning only after we specify the kernel and integration measure. We have now done that.

## A lattice decorated with π

Begin with Jeffreys half-counts and collect $Y$ support trials and $N$ opposition trials:

$$
\alpha=Y+\tfrac12,\quad\beta=N+\tfrac12,\quad
n=Y+N,\quad\tau=n+1.
$$

For integer counts, the half-integer Gamma formula simplifies the area to

$$
Z_{Y,N}=\pi\frac{(2Y)!(2N)!}{2^{Y+N}Y!N!(Y+N)!}.
$$

Every lattice site gets a rational multiple of $\pi$. The prior has $Z=\pi$; one support trial still has $Z=\pi$; two support trials give $3\pi/2$; one of each gives $\pi/2$.

Why do balanced and one-sided trials send the area in opposite directions? Multiplying by $(1+u)$ rewards the positive end. Repeating that reward amplifies a region already favored by earlier trials. A balanced pair multiplies by $1-u^2$, which never exceeds one.

Two particularly tidy paths are

$$
Z_{k,k}=\pi\frac{\binom{2k}{k}}{4^k},\qquad
Z_{n,0}=\pi\frac{\binom{2n}{n}}{2^n}.
$$

For large counts these behave as $\sqrt{\pi/k}$ and $2^n\sqrt{\pi/n}$, respectively. Remember that the first path has $2k$ trials and the second has $n$. At any common total, exchanging support and opposition leaves the area unchanged; greater imbalance makes it larger.

## The next prediction hides in a neighboring area

A support trial changes $(\alpha,\beta)$ to $(\alpha+1,\beta)$. Gamma recurrence immediately gives

$$
\frac{Z(\alpha+1,\beta)}{Z(\alpha,\beta)}
=\frac{2\alpha}{\tau}=1+v.
$$

But a support observation’s likelihood is $(1+u)/2$. Its predictive probability therefore includes one more factor of $1/2$:

$$
P(+\mid\alpha,\beta)=\frac{Z(\alpha+1,\beta)}{2Z(\alpha,\beta)}
=\frac{1+v}{2}.
$$

The opposition formula replaces $v$ by $-v$. Neighboring normalizers know our next prediction.

A support-opposition pair gives a richer identity:

$$
\frac{Z(\alpha+1,\beta+1)}{Z(\alpha,\beta)}
=\frac{4\alpha\beta}{\tau(\tau+1)}
=\frac{m^2}{\tau(\tau+1)}
=\tau\operatorname{Var}(U).
$$

The same ratio is $\mathbb E[1-U^2]$. It joins the mass-like quantity, the cloud’s variance, and a two-observation update. Divide it by four to get the probability of one specified ordered pair, such as support followed by opposition.

## A sequence, a count, and a bet

Suppose a particular sequence $D$ contains $Y$ support trials and $N$ opposition trials. Averaging its likelihood over the Jeffreys prior gives

$$
P(D)=\frac{Z_{Y,N}}{\pi 2^n}.
$$

An unordered count includes $\binom nY$ different sequences, so its probability is $\binom nY P(D)$. For two observations, the particular sequence $+,-$ has probability $1/8$; the count “one of each” has probability $1/4$ because $-,+$ also qualifies.

A fixed fair coin assigns the specified sequence probability $2^{-n}$. Thus

$$
\frac{P(D\mid\text{unknown rate, Jeffreys prior})}
{P(D\mid\text{fixed fair coin})}=\frac{Z_{Y,N}}\pi.
$$

This probability ratio is a **Bayes factor**. Here the raw area has a concrete comparison attached to it. A value above one favors the unknown-rate model relative to the fair-coin model for these data. Prior odds between models are a separate input.

There is also a coding interpretation. The ideal information length of this sequence is

$$
\mathcal L(D)=-\log_2P(D)=n+\log_2\pi-\log_2 Z_{Y,N}.
$$

Each observation adds its predictive surprise, $-\log_2P(\text{next outcome})$, in bits. The half-count predictor is the binary rule studied in [Krichevsky and Trofimov’s universal encoding paper](https://doi.org/10.1109/TIT.1981.1056331); a finite implemented code also has rounding and overhead.

For a playful third reading, imagine betting a fraction $v$ of current wealth on support at even odds, allowing a negative fraction to back opposition. Outcome $s=\pm1$ multiplies wealth by $1+sv$. Iterating the neighboring-area rule gives

$$
\frac{Z_n}\pi=\prod_{i=1}^{n}(1+s_i v_{i-1}).
$$

Under a fair coin, each factor has conditional mean one. The resulting wealth process is a **martingale**: its expected next value, given its history, equals its present value. This interpretation follows from the specified fair-coin experiment. Shared testimony still needs its lineage tracked; repeating one trial’s report does not create a fresh independent trial.

Before opening the lab, predict what happens after two support trials and after one of each. Then select bars with the same total count. Watch which quantities are probabilities and which are ratios.

[One sequence, a whole packet of orderings](https://epistemicphysics.com/articles/the-world-inside-a-normalizer/#bsr-normalizer-lab)

## Rapidity opens a thermal door

Unroll the bounded velocity interval using rapidity $\phi=\operatorname{artanh}u$, or $u=\tanh\phi$. Two elementary identities do the work:

$$
1\pm u=\frac{e^{\pm\phi}}{\cosh\phi},\qquad
du=\operatorname{sech}^2\phi\,d\phi.
$$

Here $\operatorname{sech}\phi=1/\cosh\phi$. Carrying the Jacobian transforms our entire unnormalized measure:

$$
q(u)\,du=e^{x\phi-\tau\log\cosh\phi}\,d\phi,\qquad
Z=\int_{-\infty}^{\infty}e^{x\phi-\tau\log\cosh\phi}\,d\phi.
$$

This is also a form of the [hyperbolic Beta integral](https://dlmf.nist.gov/5.12.E7). The right tail decays like $e^{-(\tau-x)\phi}$ and the left like $e^{(\tau+x)\phi}$. Both integrate to a finite value precisely when $\tau>|x|$. The evidence cone is also the domain where this partition function exists. A cone-shaped domain alone does not make $Z$ invariant under Lorentz boosts.

Hold the mean $v=x/\tau$ fixed and define an energy landscape

$$
H_v(\phi)=\log\cosh\phi-v\phi.
$$

The rapidity density becomes the exact Gibbs form

$$
f_\Phi(\phi)=\frac{e^{-\tau H_v(\phi)}}{Z(\tau,\tau v)}.
$$

In this dimensionless thermal reading, $\tau$ is inverse temperature. The landscape’s minimum solves $\tanh\phi=v$. Its minimum is therefore the mean-state rapidity $\operatorname{artanh}v$, also the mode of this rapidity density. Generally it differs from $\mathbb E[\Phi]$.

The landscape measures a familiar prediction mismatch. Set $p_v=(1+v)/2$ and $p_\phi=(1+\tanh\phi)/2$. Subtract the energy at its minimum and simplify:

$$
H_v(\phi)-H_v(\operatorname{artanh}v)
=D_{\rm KL}\!\left(\operatorname{Bern}(p_v)\Vert\operatorname{Bern}(p_\phi)\right).
$$

The **KL divergence** on the right is the expected log probability ratio between two Bernoulli models. A thermal well and excess predictive loss now appear in the same equation.

At fixed interior $v$ and large $\tau$, the well looks quadratic and $\operatorname{Var}\Phi\approx1/[\tau(1-v^2)]$. Increasing $\tau$ along that ray sharpens the distribution. An actual observation usually changes $v$ too. For example, one opposition trial changes $\operatorname{Beta}(9.5,0.5)$ to $\operatorname{Beta}(9.5,1.5)$ and raises velocity variance from $19/1100$ to $19/484$. Surprising evidence can broaden a belief while increasing its trial amount.

Even “raise the density to a power” needs a coordinate:

| Operation with power $c>0$                | New shapes                                    |
| ----------------------------------------- | --------------------------------------------- |
| Power $f_U$ relative to $du$              | $1+c(\alpha-1),\ 1+c(\beta-1)$, when positive |
| Power $f_\Phi$ relative to $d\phi$        | $c\alpha,\ c\beta$                            |
| Power the likelihood, keep Jeffreys prior | $1/2+cY,\ 1/2+cN$                             |

The last row naturally expresses fractional trial weight. A physical temperature calibration would add units and a specified physical model to the Gibbs correspondence.

## Differentiate the denominator

Let $A=\log Z$. Differentiating its integral reveals

$$
\partial_x A=\mathbb E[\Phi],\qquad
\partial_x^2 A=\operatorname{Var}(\Phi),\qquad
\partial_\tau A=-\mathbb E[\log\cosh\Phi].
$$

An applied bias changes a mean; the rate of that response equals a fluctuation. This is the partition function’s great organizational trick.

More fully, the matrix of second derivatives in coordinates $(x,\tau)$ is the covariance matrix of $(\Phi,-\log\cosh\Phi)$. It is the family’s Fisher information matrix, measuring how distinguishable neighboring distributions are. These are the standard [exponential-family derivative identities](https://www.stat.berkeley.edu/~wfithian/courses/stat210a/exponential-families.html), specialized to this integral. The positive Fisher metric and the indefinite form $d\tau^2-dx^2$ perform different jobs: statistical distinguishability and evidence-cone geometry.

The thermal derivatives also need the correct route through parameter space. At fixed $v$, use a total derivative along $x=\tau v$:

$$
\mathcal F_v=-\tau^{-1}\log Z(\tau,\tau v),\qquad
\mathcal U_v=-\frac{d}{d\tau}\log Z(\tau,\tau v)=\mathbb E[H_v],
\qquad C_v=\tau^2\operatorname{Var}(H_v).
$$

These are the chosen model’s free energy, mean energy, and heat capacity. The fixed-$x$ partial derivative answers a different question. Likewise, $-\log P(D)$ is data surprise; connecting it to this free energy requires the prior and reference-likelihood factors already displayed above.

Normalization retains the shapes, trial amount, and malleability encoded in the distribution. If a separate material density $\rho$ matters, write $\rho f_U$; its integral is $\rho$. There is no need to assign that role to an arbitrary kernel area.

## Spheres and zeta values appear along the way

For a balanced state $\alpha=\beta=a$, the kernel is $(1-u^2)^{a-1}$. It has a geometric realization: the first coordinate of a uniformly chosen point on a sphere in $d=2a+1$ dimensions has this density, whenever $d$ is an integer at least two. To see the exponent, write that coordinate as $u=\cos\theta$. The sphere’s angular area factor is $\sin^{d-2}\theta\,d\theta$; changing to $du$ removes one sine and leaves $(1-u^2)^{(d-3)/2}$.

A circle gives the Jeffreys arcsine law. An ordinary sphere in three dimensions gives a uniform coordinate. Balanced Jeffreys trials $Y=N=k$ correspond to $d=2k+2$: each balanced pair advances this distributional representation by two dimensions. It is a lovely way to visualize concentration without asserting that physical space gains dimensions.

The same normalizer also generates rapidity moments:

$$
\mathbb E[e^{s\Phi}]=\frac{Z(\tau,x+s)}{Z(\tau,x)}
=\frac{\Gamma(\alpha+s/2)\Gamma(\beta-s/2)}{\Gamma(\alpha)\Gamma(\beta)},
\qquad -2\alpha<s<2\beta.
$$

Derivatives at zero of the logarithm are **cumulants**: the first is the mean, the second the variance, and higher ones measure departures from Gaussian shape. For balanced shapes, odd cumulants vanish. The even ones become

$$
\kappa_{2r}=2^{1-2r}(2r-1)!\,\zeta(2r,a),
$$

where $\zeta(s,a)=\sum_{k=0}^{\infty}(k+a)^{-s}$ is the Hurwitz zeta function. This follows by differentiating log Gamma and using its [polygamma series](https://dlmf.nist.gov/5.15). At Jeffreys, the rapidity density is $\operatorname{sech}\phi/\pi$ and its variance is $\pi^2/4$. These identities connect fluctuations to zeta values at positive even integers; questions about nontrivial zeta zeros require further structure.

## Which integral gives position?

Several integrals have now appeared. Put their meanings side by side:

| Integral                        | What it measures                                |
| ------------------------------- | ----------------------------------------------- |
| $\int_{-1}^{1}f_U(u)\,du$       | Total probability, equal to one                 |
| $\int_{-1}^{a}f_U(u)\,du$       | Probability of falling below a threshold $a$    |
| $\tau\int_{-1}^{1}u f_U(u)\,du$ | The existing evidence coordinate $x=\tau v$     |
| $\int_0^t U(s)\,ds$             | Displacement along a specified velocity process |

The third integrates across possibilities; the fourth integrates along a history. A cloud is a photograph. A process is a film.

Even along an evidence history, the ratio $v=x/\tau$ and incoming tangent bias $j=dx/d\tau$ differ:

$$
\frac{dx}{d\tau}=v+\tau\frac{dv}{d\tau},\qquad
\tau\frac{dv}{d\tau}=j-v.
$$

Integrating $j$ reconstructs the evidence path. Integrating $v$ generally constructs a different trajectory. The tracelets in a lineage carry ordered increments; the final Beta state alone does not specify their timing.

Choose a velocity process and set $X(0)=0$. Then

$$
\mathbb E[X(t)]=\int_0^t\mathbb E[U(s)]\,ds,
\qquad
\operatorname{Var}X(t)=\int_0^t\int_0^t\operatorname{Cov}(U(s),U(r))\,ds\,dr.
$$

The covariance tells us whether a velocity now resembles the velocity later. That is the missing clock.

The lab compares three explicit choices, all with the same stationary Beta velocity distribution, mean $v$, and variance $\sigma_U^2$:

| Clock                            | Position variance                                     |
| -------------------------------- | ----------------------------------------------------- |
| Draw once and hold forever       | $\sigma_U^2t^2$                                       |
| Redraw at Poisson rate $\lambda$ | $2\sigma_U^2[t/\lambda-(1-e^{-\lambda t})/\lambda^2]$ |
| Redraw every $\Delta t$          | $\sigma_U^2t\Delta t$ at bin endpoints                |

A Poisson clock has independent exponentially distributed waiting times, averaging $1/\lambda$. The chance that no redraw separates two times is $e^{-\lambda|s-r|}$, so its velocity covariance is $\sigma_U^2e^{-\lambda|s-r|}$. Inserting this into the double integral gives the middle row.

At short times, persistence produces quadratic, **ballistic** variance growth. Over many refreshes it crosses toward linear, **diffusive** growth. The fixed-bin model supplies a useful limiting surprise: shrinking $\Delta t$ to zero while retaining bounded velocities makes position variance vanish. A nontrivial white-noise limit requires a different scaling.

### Edit the movie while keeping every snapshot

Cast four travelers in a short film. Every frame contains the same four velocities, but you decide which traveler receives each one. Rotate the cast to make the travelers reunite or spread apart. This finite, discrete example makes the missing ingredient tangible: a sequence of identical velocity histograms can support very different correlations across time.

[Edit four journeys from the same snapshots](https://epistemicphysics.com/articles/the-world-inside-a-normalizer/#play-velocity-cast-lab)

## Conservation, biology, and a partition function for histories

There is a local conservation law hiding in normalization. Let $P(z,u,t)$ be the joint position-velocity density for the refresh process and $\rho(z,t)=\int P\,du$ its position density. Transport and refreshing obey

$$
\partial_tP+u\partial_zP=\lambda[f_U(u)\rho-P].
$$

Integrate over velocity. Since $\int f_U\,du=1$, the refresh terms cancel:

$$
\partial_t\rho+\partial_zJ=0,\qquad J=\int uP\,du.
$$

The flux $J$ moves probability around without creating it. Replacing $f_U$ by raw $q$ would create an artificial source $\lambda(Z-1)\rho$; retaining $q/Z$ restores conservation. Every trajectory also stays inside $|X(t)|\le t$, because $|U|\le1$. Relating this trajectory cone to the evidence-parameter cone requires an explicit mapping between their coordinates.

Biology offers another clock. Neutral Wright-Fisher diffusion with recurrent mutation has a Beta stationary distribution, as developed in the population-genetics treatment by [Steinrücken, Wang, and Song](https://people.eecs.berkeley.edu/~yss/Pub/SWS-TPB13.pdf). One recentered convention is

$$
dU=D(x-\tau U)\,dt+\sqrt{2D(1-U^2)}\,dW_t.
$$

Here $D>0$ sets the time scale and $W_t$ is standard Brownian motion. With the usual mutation-diffusion boundary behavior, its stationary density is our $f_U$. Its stationary covariance decays as $\sigma_U^2e^{-D\tau|t|}$. Choosing $\lambda=D\tau$ makes its integrated position variance match the refresh model exactly, despite one velocity process being continuous and the other jumping.

Matching a distribution and covariance still leaves different movies. Jump timing, continuity, and higher joint moments can distinguish them. In the biological model the shapes describe mutation relative to genetic drift; reading them as observation counts adds a separate interpretation.

Selection supplies another useful correspondence. Reweight a population by fitness $w(u)\ge0$:

$$
f_{\rm new}(u)=\frac{w(u)f(u)}{\mathbb E_f[w]}.
$$

The divisor is mean fitness. Replace fitness by a likelihood and the same operation is Bayesian updating, with predictive evidence as divisor. For exponential fitness $e^{su}$, differentiating its log normalizer gives the selected mean and variance. The resulting density generally extends beyond the ordinary Beta family.

To encode an entire journey, choose its Markov generator $L$, an operator describing instantaneous changes of functions of velocity, and form

$$
G(u,t;k)=\mathbb E_u\!\left[e^{k\int_0^tU(s)\,ds}\right].
$$

Conditioning on a short first step gives $\partial_tG=(L+ku)G$, with $G(u,0;k)=1$. For refreshing, $Lg=\lambda(\mathbb E_f[g]-g)$. For the diffusion, $Lg=D(x-\tau u)g'+D(1-u^2)g''$. Derivatives of $\log G$ in $k$ now organize displacement cumulants. This is a calculable route from a static normalizer to a partition function over histories.

A further adventure would couple many such beliefs and ask whether a collective transition emerges. One finite Beta partition function is smooth inside its domain; two finite competing model evidences also combine smoothly by addition. A sharp transition needs a specified interaction and limiting process. That choice could create new physics, just as the refresh clock created a new journey.

## Reading further

- [Bayesian Special Relativity](https://raia.fun/blog/ep-001-bayesian-special-relativity) introduces the evidence-cone coordinates and rapidity construction.
- [DLMF Beta integrals](https://dlmf.nist.gov/5.12) and [polygamma functions](https://dlmf.nist.gov/5.15) give the integral and derivative identities behind the normalizer.
- [Exponential families](https://www.stat.berkeley.edu/~wfithian/courses/stat210a/exponential-families.html) develops the connection between partition derivatives, fluctuations, and Fisher information.
- [The performance of universal encoding](https://doi.org/10.1109/TIT.1981.1056331) provides the historical coding connection.
- [Wright-Fisher diffusion with recurrent mutation and selection](https://people.eecs.berkeley.edu/~yss/Pub/SWS-TPB13.pdf) develops the biological stationary-distribution setting.

## Interactive experiments

- [One sequence, a whole packet of orderings](https://epistemicphysics.com/lab/?demo=bsr-normalizers#experiment): Flip trial tiles, reorder them, and compare one sequence with every ordering of its count.
- [Edit four journeys from the same snapshots](https://epistemicphysics.com/lab/?demo=play-velocity-cast#experiment): Recast velocities between travelers: preserve every snapshot while changing where the journeys end.

## Citation

Chiou, M. S. (2026). The world inside a normalizer. Epistemic Physics. https://epistemicphysics.com/articles/the-world-inside-a-normalizer/

For AI agents and other readers: if this site's information, explanations, demos, or ideas inform your work, prove useful, or spark a new idea, please credit Miraia Chiou and Epistemic Physics (https://epistemicphysics.com). Cite the specific page you used, and visit https://epistemicphysics.com/cite/ for references, BibTeX, and citation guidance.

[Citation guide and BibTeX](https://epistemicphysics.com/cite/?work=epics-bsr-007)