Applying a Markov kernel to two initial laws cannot increase their total variation distance, even for the distribution of an entire subsequent finite trajectory. If a bounded summand has absolute value at most , the variance of a length- path sum differs between these laws by at most times their total variation distance. Thus a geometrically long enough burn-in makes this difference negligible relative to .
If , independence gives . Since , this proves the displayed identity. The Cauchy-Schwarz inequality also gives . Thus perturbations of size have a uniformly bounded joint total variation distance, a useful input to a metric squared-loss two-point bound.
A sequence is contiguous with respect to if implies for every measurable sequence of events. It transfers negligible-event assertions, including consistency, between changing sampling laws. It need not imply small total variation distance.
Coupling characterization of total variation distance Created 2026-09-24 Updated 2026-10-14
A coupling of probability distributions is a pair of random variables with those marginal distributions. On a finite set, and more generally on a standard Borel space, the least achievable mismatch is
A maximal coupling attains the minimum; maximal here refers to agreement. An arbitrary coupling can have a much larger mismatch.[ref]
For the four-outcome example, assign the following joint probability mass: gets , gets , gets , and gets . The remaining must have ; place on and on . Row sums give , column sums give , and the diagonal contains . Agreement can reach 70%; mismatch cannot fall below 30%.
Figure 1. . The original joint-mass example retains all common probability mass on its diagonal. The two off-diagonal entries account for the minimum mismatch . Details.
To see the lower bound, for any event ,
The two indicator functions agree whenever . To attain the bound on a finite set, keep mass at and couple the remaining excesses. The two residual probability distributions live on disjoint subsets, so their part necessarily mismatches. When the probability distributions are identical, the residual part disappears entirely.
This is a choice of joint construction, not a claim that independently drawn observations are usually equal. Two independent random variables with the same continuous probability distribution can disagree with probability one even though their total variation distance is zero. The construction of a maximal coupling deliberately introduces statistical dependence to reveal how similar the marginal laws are.
Doeblin's condition Created 2026-09-28 Updated 2026-10-03
A Doeblin condition is a minorization of the whole state space. It implies uniform geometric convergence in total variation distance.
If for a finite irreducible Markov chain with at least two states, then its worst-case total variation distance from its stationary distribution satisfies . Choose a possibly complex eigenfunction for an eigenvalue of modulus , normalize , and start where the maximum is attained. Since , ; the defining sum for total variation distance bounds this difference by . For , this yields for , because .
When two density operators are commuting elements, choose a common orthonormal eigenbasis and write , . Their quantum fidelity becomes the classical overlap of the two probability distributions, while their trace distance becomes total variation distance:
This follows by computing the positive square root of an operator entry by entry in the common eigenbasis. For the qubit mixed states and ,
The two contributions to quantum fidelity are about and . In contrast, the chance that two independent computational basis measurements, one on each preparation, print the same label is . That particular event has a specified probability; neither nor is that probability.
Figure 1.
Commuting mixed-state probabilities, square-root overlap contributions to unsquared fidelity F, squared fidelity Fsq, and trace distance D
. Details.
Here trace distance means an equal-prior quantum state discrimination success of by the Holevo–Helstrom theorem. Measuring in the common computational basis and guessing on outcome zero and on outcome one achieves . Also : the pure-state trace distance and fidelity identity is an equality for pure states, rather than a formula for arbitrary mixed states.
For independent observations with a normal distribution , fixed, every estimator has . Compare and : the joint Kullback-Leibler divergence is , so the total variation–Hellinger–relative entropy inequality bounds total variation distance by . Apply the Le Cam lower bound under absolute-error loss. Degenerate zero-noise laws do not satisfy a positive lower bound.
Helstrom measurement for the zero and plus states Created 2026-09-24 Updated 2026-10-14
Suppose one qubit is prepared in either of the known pure states and , with equal prior probabilities. Measuring in the computational basis seems natural, but a tilted projective measurement gives better quantum state discrimination. The two density operators and their difference are
The Pauli matrices give , and the eigenvalues of are . Their absolute sum gives trace distance , agreeing with trace distance between qubit states.
The positive spectral projector points along the unit Bloch vector :
The outcome labels and here refer to the sign of the eigenvalue of ; they do not denote a Hadamard basis measurement. Guess after and after . This sign choice is the equal-prior Helstrom measurement for two pure states, a case of the Holevo–Helstrom theorem. [ref]
An explicit orthonormal basis for this projective measurement is
with . Thus the real coefficient rotation angle is relative to the computational basis, while the Bloch vector of is tilted by from the positive axis toward negative . State-vector and Bloch sphere angles differ by a factor of two.
Write and , so . The Born rule gives the conditional outcome probabilities
The event has probability gap , attaining the variational characterization of trace distance. Both preparations are identified correctly with conditional probability , so
For comparison, a projective measurement in the computational basis has conditional outcome probabilities on and on . Guess zero on outcome zero and plus on outcome one. The average success is
and the measured total variation distance is only . The optimal tilted projective measurement gains about in absolute success probability by balancing the two conditional error rates instead of identifying perfectly at the expense of .
Figure 1.
Conditional outcome probabilities and equal-prior success for Z and optimal Helstrom measurements of zero versus plus
. Details.
These are statistics for alternative measurements on independently prepared copies. Applying first and then the tilted projective measurement is a different experiment, because the first measurement in quantum mechanics changes the density operator.
Maximal coupling Created 2026-09-24 Updated 2026-10-14
A maximal coupling of two probability distributions maximizes the agreement probability among all their couplings. Equivalently, it minimizes . The coupling characterization of total variation distance identifies the minimum mismatch with the total variation distance of the two marginal laws. The term maximal concerns agreement, rather than the size of a correlation coefficient or closeness in a geometric metric.
On a finite set, retain all common probability mass at equal pairs and match the excesses elsewhere. The explicit four-outcome example in coupling characterization of total variation distance realizes this construction and checks its marginal distributions. A maximal coupling need not be unique: different joint arrangements of the residual mass can have the same optimal mismatch.
For estimating the expected value over all continuous probability density functions on , compare and . Their means differ by , and . The chi-squared divergence of product measures bounds joint total variation distance by . The metric squared-loss two-point bound then gives the displayed uniform minimax risk lower bound.
For any estimator in a metric parameter space, by the triangle inequality. Average the two risk functions and replace their probability density functions by their minimum. This proves . The overlap equals one minus total variation distance. The same proof applies to squared error for a real-valued statistical functional of the parameter, without requiring that functional to be injective.
We prove uniqueness in the sense of at most one Gibbs measure satisfying the printed conditional distributions. No shift-invariance assumption is needed in advance. Write a history as . The variation hypothesis makes continuous and bounded on the finite-alphabet sequence space; let and . The conditional distribution of the first spin, obtained by marginalizing the prescribed two-spin law, must be
In particular for every symbol and history. Any two candidate measures have this same time-homogeneous one-step kernel, by conditional marginalization and the tower property.
If agree in their most recent symbols, each two-spin energy in this formula changes by at most
Both numerator and denominator therefore change by factors between and . Consequently the total variation distance of the two kernels is at most , where . Uniform positivity also bounds it by . Put
These bounds are decreasing, strictly below one, and summable. We now give the summable-memory coupling argument rather than assuming a small one-step contraction coefficient.
Starting from any two histories at a remote time, use a maximal coupling of the next-spin laws. If the actual histories have agreed for consecutive symbols, their probability of disagreement is at most . Use shared independent uniform random numbers to construct a dominating match-length process : it increases from to with probability and resets to zero with probability . Since decreases, the actual consecutive-match length can be kept at least at every step. When the dominating process does not reset, the actual spins can be coupled to agree.
A new run from zero has probability
of never resetting. The product is positive because and : its tail logarithm converges absolutely, and its finite initial factors are positive. After each reset a fresh run has the same success probability, so the probability of failed runs is . There are almost surely only finitely many resets. Every failed run has finite duration, and thus the time of the final reset is finite almost surely. Its distribution depends only on , not on the two starting histories. It follows that the probability of any disagreement in a fixed future block tends uniformly to zero as the block's distance from the initial time tends to infinity.
Now take two candidate Gibbs measures, sample their respective histories before time , and couple their forward kernels as above. Their finite-block laws really are obtained by these kernels: successive conditional expectations use the common conditional distribution . For any fixed coordinate block , the coupling inequality for total variation bounds the distance between the two candidate block laws by the dominating probability of a reset at or after step . That bound tends to zero as . Hence all probabilities of finite cylinder sets agree; cylinder sets generate the product Borel sigma-algebra, so the measures agree.
Thus there is at most one measure satisfying these conditional specifications. There is a separate compatibility issue: for an arbitrary unnormalized symbolic potential the finite-horizon conditional distributions printed in the question need not be consistent as the right endpoint changes. The argument proves the requested uniqueness whenever such a law exists; it does not silently infer existence from the variation estimate. A normalized one-sided kernel or the usual compatible Gibbs specification supplies the standard existence formulation. For a concrete incompatibility, take two symbols and with . This has zero variation once two symbols agree. From a history ending in symbol one, the first-symbol probability obtained from a two-spin block is , but from a three-spin block it is . Thus the printed finite-horizon specification can indeed have no measure, even though its uniqueness implication remains valid.
Use and interpret the printed overlap term as ; its unindexed are a typographical slip. Take , as usual. If necessary, deleting from its own neighbourhood gives the proof below with a smaller nonnegative right-hand side, so the stated bound still follows.
If , every indicator is zero almost surely and the Poisson distribution with parameter zero is the same point mass. Suppose . For any , put and let have Poisson distribution with mean . Solve the Poisson Stein equation
For example, for one can take
The full weighted sum is zero, so using its complementary tail also shows this solution is bounded. The allowed Stein factor is , where .
Define
Joint independence of from all events outside its neighbourhood makes independent of , so . Adding and subtracting those terms gives
A difference over integer steps has magnitude at most . In the first sum the step count is ; in the second, on , it is . Consequently the absolute value is at most
The Stein equation identifies the left side with . Taking the supremum over is the definition of total variation distance, proving the Stein-Chen bound with the Poisson Stein factor:
For the random uniform hypergraph, index indicators by four-vertex sets. Each such set is a hypergraph clique exactly when its four triples are all present, so its probability is . Distinct sets use disjoint triple indicators unless their intersection has size three; disjoint underlying trials give the required joint independence. A fixed four-set has dependent others, and each overlapping pair requires seven distinct triples, with joint probability . Thus, with and , the bracket is
Since ,
This calculation applies for and ; if or , the count and its Poisson parameter are zero and the conclusion is immediate.
Let be the transition matrix of a finite irreducible Markov chain with at least two states, reversible with stationary distribution . Its Dirichlet form of a Markov chain is
The Discrete Poincare inequality states that for every real . Its optimal constant is , where is the second largest eigenvalue of . Indeed, detailed balance makes self-adjoint on . Expand in an orthonormal eigenbasis. Then and , proving the inequality and its sharpness.
The useful constructive version is the canonical paths Poincare bound. Choose transition graph paths between each ordered pair, write , and define
where counts occurrences of the directed transition. The Cauchy-Schwarz inequality along each path in a graph gives . Multiply by and sum. The identity proves .
For a lazy Markov chain this also proves the mixing time estimate
To see this, the centred initial density has norm at most ; laziness puts all nonconstant eigenvalues between zero and . Their powers contract that norm, and the Cauchy-Schwarz inequality bounds total variation distance by half the density's norm.
Here is a precise sufficient density condition for the permanent application: every row and every column has more than ones. Interpret the permanent of a matrix as the number of perfect matchings in its bipartite graph, with vertices per side. Write for the number of -edge matchings in a graph.
First prove the short augmenting paths in a dense bipartite graph bound. Given a nonperfect matching in a graph, choose unmatched on opposite sides. If either has an unmatched neighbour, augment with one edge. Otherwise the matched partners of and the neighbours are subsets of the same -element side whose total size exceeds . Their intersection gives an alternating path in a graph of length three, with a matching edge. Repeating proves . Choose augmentations by a fixed deterministic rule. A reverse augmentation is specified by at most four vertices, so a next-size matching in a graph has at most preimages. Therefore
Use the matching generating function and the weighted matching Markov chain on all matchings in a graph, with stationary weights . It stays put with probability ; otherwise it chooses a graph edge uniformly and proposes its deletion, its addition if both endpoints are unmatched, or its exchange for the unique incident matching edge if precisely one endpoint is matched. Other proposals stay put. Accept with probability . Reverse proposals use the opposite exchanged edge, so detailed balance holds. Deletions connect every state to the empty matching in a graph.
For completeness, the canonical paths for the weighted matching chain give a polynomial bound without assuming a mixing theorem. The symmetric difference consists of alternating graph paths and graph cycles. Order its components deterministically, orient open components from their smaller endpoint, and convert them from to by moving an unmatched endpoint with exchanges. If the first edge belongs to , delete it first. For a graph cycle, delete one -edge, move the unmatched endpoint around the opened cycle, then insert the last -edge. Each intermediate state is a matching in a graph, and each entire transition path in a graph has length at most .
Fix a directed transition with initial state . The complementary edge set agrees with one of the two matchings in a graph outside the active component and has at most two degree-two vertices on that component. Remove at most two offending edges to obtain a matching in a graph . Store those edges and one bit specifying the initial alternating colour of the active component. This encoding is injective for the fixed transition: it restores and , identifies the active component from the changed edges, and recovers the two colours using the fixed component order. Completed components have 's colour in , pending components have 's colour, and the stored bit resolves the active component. Common edges are .
There are at most auxiliary records per , where is the number of graph edges. If , then
Also . Sum over the encodings and use the path length bound to obtain
This is polynomial whenever and are polynomially bounded. Starting at the empty matching in a graph, , so the Discrete Poincare inequality gives polynomial-time sampling to any prescribed total variation distance. The addition/deletion/exchange construction is also described in Section 12.4 of Jerrum and Sinclair's survey.
It remains to turn sampling into a relative count; merely sampling perfect matchings would not do that. For , choose and . The elementary bounds above give
Choose a schedule whose successive ratios are at most . It has stages. The annealing ratios for a matching generating function satisfy
Each observable lies in , so the Hoeffding inequality gives relative accuracy with independent samples per stage. Run independent copies of the chain from the empty matching in a graph; make the total variation distance per sample at most , where is that sample count. The coupling characterization of total variation distance and the union bound make the entire sample collection coincide with exact stationary samples except with probability at most . Choose the sample-count constant so all stationary sample means are accurate except with probability at most .
Let be the product of the sample means, and output . On the successful event its multiplicative error relative to lies between and . Combining the two endpoint bounds proves
Activities, sample counts, simulation lengths and rational bit lengths are polynomial in . Thus this is a fully polynomial randomized approximation scheme. Small dimensions can be evaluated directly.
Use a fixed logarithm base , and write . All relative entropies below use this same base. Adopt and for .
For Gibbs inequality, if at some point, the relative entropy is infinite and the required conclusion is immediate. Otherwise let . The elementary inequality , with equality only at , gives, for every ,
Consequently
This argument applies to a countable discrete probability distribution too: the difference between the two sides of the pointwise inequality is nonnegative, and the negative part of the entropy summand is bounded by , hence summable. Thus no subtraction of divergent series is involved. Equality forces both and equality in every pointwise inequality; hence on and off . Conversely plainly gives zero. Therefore .
For the log-sum inequality, set , , and normalize to the probability distributions , . When is nonempty, positivity ensures , and direct expansion gives
Thus normalization reduces the log-sum inequality to Gibbs inequality. Equality holds exactly when . The empty-set case is the zero identity under the usual zero-mass convention.
For the final bound, partition the alphabet using and put , . Since the total signed difference is zero,
Apply the log-sum inequality separately to and its complement. The same normalization argument works for countable sets of finite total mass, using the countable Gibbs inequality just proved. It yields the binary partition bound for relative entropy
Zero block masses are interpreted by the stated conventions; a positive mass divided by zero makes the bound infinite. The given binary estimate now implies
This is Pinsker's inequality. Equivalently, for the total variation distance , the bound is . The factor keeps the statement valid for either bits or natural-logarithm units.