Source-linked AI summary

The information geometry of product-reference discrete diffusion: Interaction growth complexity and optimal scheduling

Martin J. Wainwright

arXiv:2608.28949v1stat.MLcs.AIcs.LGmath.ST

TL;DR

The paper addresses how to quantify and optimize sampling accuracy for product-reference diffusion on discrete distributions. It introduces interaction growth complexity, analyzes its bivariate and univariate forms, and uses them to characterize KL error and stepsize-dependent iteration complexity. The resulting framework covers arbitrary product references and connects aggregate complexity to total correlation and dual total correlation.

  • Problem

    The paper asks whether discrete diffusion sampling performance can be explained by data geometry and whether that geometry can guide, optimize, and certify sampling schemes.

  • Method

    The paper analyzes product-reference diffusion samplers through interaction growth complexity along a path from a factorized reference distribution to the target.

  • Results

    In the fine-grid limit, the log-SR IGC density determines exact KL discretization error up to a factor of two, while aggregate IGC mass and refined stepsizes characterize iteration complexity.

  • Takeaways & Limitations

    IGC provides a geometric basis for adapting stepsizes to the target distribution and relating sampling complexity to total correlation and dual total correlation.

  • Takeaways & Limitations

    Refined stepsizes based on IGC increments require reliably estimating those increments from samples, unlike the data-independent stepsize choice of the simpler guarantee.

Abstract

from arXiv · show

We study a class of product-reference diffusion algorithms for sampling from a discrete distribution. We show that their sampling performance can be characterized using a path-based measure of data geometry that we call the interaction growth complexity (IGC). We show that a bivariate IGC kernel gives an exact representation of both the KL discretization error and a simple one-step upper bound. The simpler univariate IGC density can be used to study the effect of stepsize choices on the iteration complexity required to obtain $ε$-accurate samples in KL divergence. Samplers that traverse the path with equi-spaced steps in log-squared-reliability-odds have performance that depends on the aggregate IGC mass, whereas refined choices of stepsizes have a lower complexity depending on a square-root functional. In the fine-grid limit, both of these characterizations become sharp. We also allow general product reference distributions and show that the reference law can substantially reshape the IGC profile and the resulting sampling complexity; in particular, references far from both the uniform and the data marginals can yield dimension-dependent improvements. Finally, the aggregate IGC mass admits bounds in terms of total correlation and dual total correlation, thereby connecting the pathwise geometry to classical measures of multivariate dependence.

1 Introduction

The paper develops interaction growth complexity (IGC) to characterize and optimize product-reference discrete diffusion sampling. It extends geometric analyses of diffusion to discrete distributions, general product references, and product-kernel approximations.

  • Background: Discrete diffusion replaces Gaussian perturbations with Markovian corruption kernels and approximately reverses them on the underlying alphabet.The paper analyzes one product-kernel class within the broader design space of discrete diffusion and relates it to coordinatewise reverse updates.
  • Motivation: The paper asks whether discrete diffusion performance can be quantified by data geometry and used to design and certify sampling schemes.These questions continue the geometric perspective developed for Gaussian and masking diffusion.
  • Contributions: IGC extends information-geometric complexity measures to product-reference diffusion samplers for discrete distributions.The samplers trace a noising path from an arbitrary factorized reference distribution to the target using product approximations.
  • Contributions: The true IGC structure is bivariate, unlike the essentially univariate complexity measures used for Gaussian and masking diffusion.A univariate IGC density still provides intuition and guarantees, while the bivariate kernel connects to total correlation and dual total correlation.
  • Related work: Product-reference diffusion permits arbitrary coordinate-dependent reference components, extending beyond uniform and absorbing-state references.Related work reports performance differences between these standard choices and matched marginal references.
  • Related work: The analysis connects product-kernel sampling errors to classical multivariate-dependence measures, including total correlation and dual total correlation.This places the paper's pathwise complexity perspective alongside existing non-asymptotic and geometric analyses of discrete diffusion.

2 Overview: Interaction growth complexity and sampling

The paper studies product-reference diffusion samplers through interaction growth complexity (IGC), which tracks information along a path from a product reference to the target. IGC characterizes discretization error and guides stepsize choices, while examples show how its profile reflects dependence geometry.

  • Product-reference diffusion: PRD samplers move from a coordinate-factorized reference distribution to the target using product-based approximations to exact transition kernels.The reference may have coordinate-dependent marginals, with the uniform law as a special case.
  • Product-reference diffusion: The squared-reliability path starts at the product reference and ends at the target, with coordinates independently retained or replaced by fresh reference samples at each transition.The semigroup identity makes the construction a consistent Markov process.
  • Interaction growth complexity: The univariate IGC density q summarizes pathwise information growth, while the underlying exact complexity is bivariate and determines KL discretization error and one-step bounds.In the fine-stepsize limit, q determines the exact KL discretization error up to a factor of two.
  • Sampling complexity: Constant-spacing steps in log-SR-odds yield complexity controlled by aggregate IGC mass, whereas optimized stepsizes exploit uneven q through a lower square-root-based complexity.The fine-grid iteration estimate is sharp up to a constant prefactor, and the optimization gain grows when q is unevenly distributed.
  • Reference choice: The reference distribution is a design parameter because changing it reshapes q and can substantially alter the coarse and fine-partition complexities.The paper therefore analyzes complexity separately for each product reference ν.
  • Geometric examples: For the noisy repeated-bit model, q peaks under strong shared latent dependence and vanishes when η = 1/2 makes the target product; Curie–Weiss profiles develop modes near criticality and can become bimodal.Hierarchical prototypes similarly produce multiple IGC peaks associated with distinct tree-level Hamming separations.

3 Product-reference diffusion and KL guarantees

The section introduces PRD samplers and explains that KL error bounds depend on IGC increments and squared-reliability odds, making log-SR-odds spacing a key design choice.

  • PRD samplers use product approximations to transition kernels and can be parameterized through single-site denoisers or leave-one-out versions.

3.1 Sampling using product-reference diffusion (PRD)

PRD sampling approximates exact reference-to-output transitions with coordinatewise product kernels over a squared-reliability grid, whose spacing is the main design parameter.

  • The exact kernel K_a,b is approximated by the product of its one-coordinate transition marginals.
  • Single-site transitions can use ordinary posteriors or leave-one-out posteriors, with the ordinary posterior admitting a one-site bridge representation.
  • For a fixed iteration budget N, the stepsize grid is the key design parameter, with Theorem 1 and Corollary 1 providing guidance for arbitrary and equi-spaced grids.

3.2 Controlling the Kullback–Leibler error

The section develops KL-error guarantees for PRD samplers, including arbitrary-grid bounds, log-SR-odds schedules, learned-posterior penalties, and reference-dependent complexity effects.

  • Theorem 1 gives a KL error upper bound for any grid, involving IGC increments and odds factors, together with initialization and termination costs.
  • Initialization and termination costs can be controlled by symmetric endpoints, while sampler complexity depends only logarithmically on r_0.
  • Replacing exact posteriors with learned estimates adds an additive block denoiser penalty to the exact-posterior discretization bound.
  • An equi-spaced log-SR-odds grid yields an ε-accurate sample in KL divergence under an iteration guarantee governed by aggregate IGC complexity.
  • Aggregate IGC admits bounds through total correlation and dual total correlation, but these bounds can be conservatively loose by an Ω(d) gap.
  • Reference choice can change iteration complexity with dimension: in Ensemble B, uniform-reference IGC stays constant-order while absorbing-state IGC grows as Θ(log d).
  • For Ensemble A with ε_d = d^-1/2, the marginal-matched reference becomes increasingly unfavorable relative to the uniform and absorbing-state references.

3.3 The bivariate IGC kernel and its structure

The bivariate IGC kernel provides the paper’s central geometric representation: it exactly characterizes one-step KL error and IGC increments, while connecting aggregate IGC to dependence measures.

  • Unlike prior diffusion complexity measures, PRD complexity is fundamentally bivariate, although a univariate IGC density remains useful for intuition and guarantees.
  • The kernel Q_IGC is formed by mixed partial derivatives of log-SR-odds cut mutual information aggregated across coordinates.
  • On the diagonal, the bivariate kernel equals the univariate density: q(η) = Q_IGC(η, η).
  • Proposition 1 gives integral representations of the KL product-reference error, IGC increments, total correlation, and dual total correlation.
  • An exponential off-diagonal sandwich for Q_IGC yields upper bounds on aggregate IGC in terms of total correlation and dual total correlation.

4 Multi-block analysis and certification

The multi-block analysis adapts PRD stepsizes to blockwise IGC increments, yielding near-optimal guarantees and data-certified schedules. Bregman-divergence identities support estimating both IGC increments and KL error from samples.

  • Adaptive schedules: Adaptive multi-block schedules tailor log-odds spacing to the underlying IGC density and estimate its increments from data.The analysis introduces a sandwich estimator based on a Bregman divergence that also supports practical KL-error estimation.
  • Adaptive schedules: A K-block sampler with explicitly chosen within-block odds multipliers uses at most N rounds and returns a KL guarantee.The allocation is within a factor of four of the best bound for the given partition.
  • Partition refinement: Partition refinement can only improve the coarse complexity, with strict improvement when blockwise IGC mass is not proportional to log-odds length.For the full symmetric interval, the unpartitioned value equals the single-block complexity.
  • Data certification: Under mild moment conditions, the confidence bound for block k decreases as O(m^-1/2) with the number m of clean samples.The bound depends on both the block and the number of clean samples.
  • Data certification: The data-certified guarantee replaces unknown IGC increments with estimated increments and confidence tolerances under a simultaneous confidence event.The resulting certified odds multipliers determine the within-block log-odds spacing.

P CIGC(P), (34)

Refining the partition leads to a fine-partition complexity identified with a square-root integral, and this quantity governs the sharp leading-order error of optimal discretizations.

  • Fine-partition limit: The fine-partition limit is the infimum over finite partitions and equals the square-root integral defining the partition complexity.The same quantity governs the sharp leading-order Euler error.
  • Fine-partition limit: Theorem 3 establishes this equivalence under continuity of the relevant mixed partial derivatives and strict positivity of q on the log-odds region.These are the stated regularity conditions for the fine-partition result.
  • Euler optimality: For the log-odds-uniform grid, the sum of exact KL one-step errors has the fine-partition characterization.The theorem also gives the corresponding optimal N-step discretization result.

5 Proofs

The proofs connect PRD’s global KL error to one-step transition discrepancies and represent those discrepancies exactly through the bivariate IGC kernel. Bregman and information-theoretic identities then establish the required bounds and pathwise controls.

  • Global KL control: Theorem 1 follows by decomposing the sampler’s error across grid intervals and bounding each one-step product-transition KL error.The proof uses data processing, the KL chain rule, and Lemma 2.
  • IGC representation: The one-step KL divergence has an exact representation through the bivariate IGC kernel Q_IGC, with a pointwise small-step expansion.The representation also yields an equivalence used in the proof.
  • Anchored IGC: The anchored IGC increment is related to both the exact KL error and the interval IGC increment through its Bregman-based properties.The proof derives these relations by integrating the relevant derivatives over the log-odds interval.
  • Pathwise control: The pathwise control follows from a sandwich relation for divergences indexed by log-SR-odds, together with nonnegativity obtained from data processing.The resulting bounds control the evolution of the path-dependent divergence across reliability levels.
  • Bregman structure: The proof identifies the KL derivative along the path with a Bregman divergence because the path map is linear and KL itself is a Bregman divergence.This correspondence supplies the explicit representation needed for the anchored IGC analysis.

6 Discussion

The discussion presents IGC as a pathwise information-geometric measure that links PRD sampling complexity, optimal scheduling, KL discretization, and multivariate dependence. It also emphasizes that the reference distribution changes the induced complexity profile and sampling behavior.

  • Main implications: IGC measures the structure of PRD paths and provides iteration-complexity bounds for constant and square-root-optimized stepsizes.The measure is defined through the way dependence is removed along the noising path.
  • Main implications: The theory supplies estimators for IGC increments and KL error using a Bregman divergence, unifying local factorization error, discretization, scheduling, and dependence.These are presented as connected consequences of the pathwise framework.
  • Reference distributions: Allowing arbitrary product references makes each reference induce a different IGC complexity and provides a way to study how reference choice affects sampling complexity.The discussion identifies reference choice as an important design parameter.

A Parameterization in terms of the LOO posterior

The ordinary and leave-one-out posteriors are connected through Bayes’ rule, enabling an equivalent leave-one-out representation after substitution.

  • Bayes’ rule relates the ordinary posterior to the leave-one-out posterior.
  • Substituting the Bayes-rule relation into (8a) produces an equivalent leave-one-out representation.

B Proof of Proposition 1

The proof expresses total correlation and dual total correlation through the bivariate IGC density by differentiating their pathwise definitions and integrating the resulting identities.

  • The proof targets representations of total correlation and dual total correlation in terms of the bivariate IGC density (20b).
  • Defining T(ξ) and D(η) as total correlation and dual total correlation along the path, mutual information gives an expression for Ci(ξ, η).
  • Differentiating the definitions of total correlation and dual total correlation yields identities involving the bivariate IGC density.
  • Boundary relations for Ci allow the differentiated identities to be combined into integral identities.
  • Because T(−∞)=D(−∞)=0 and the endpoint values equal TC(Z) and DTC(Z), the claims follow from the fundamental theorem of calculus.
Loading 2608.28949v1…