Propagation of Chaos in Rényi Divergences

• The paper establishes a sharp O(d q²/N²) convergence rate for propagation of chaos using log-Sobolev inequalities and a Donsker–Varadhan splitting approach.
• It extends the entropic framework beyond KL divergence by rigorously characterizing fluctuations and concentration phenomena in weakly interacting diffusions.
• The methodology relies on strong isoperimetry and weak interaction assumptions, with illustrative Gaussian examples confirming the optimality of the derived scaling laws.

Propagation of chaos in Rényi divergences quantifies the extent to which finite-particle approximations to mean-field interacting diffusion systems decorrelate as the system size $N$ grows, using the $q$-Rényi divergence as a measure of marginal independence. Recent work establishes sharp $O(d q^2/N^2)$ rates for the stationary measures of weakly interacting diffusions, extending the entropic framework previously developed for Kullback–Leibler (KL) divergence. This approach rigorously characterizes fluctuations and concentration phenomena in high-dimensional and high-particle-number regimes, under conditions of strong isoperimetry and weak particle interaction (Zhang, 15 Jan 2026).

1. $q$-Rényi Divergence and Its Role

For probability measures $\mu,\nu$ on $\mathbb{R}^d$ with $\mu\ll\nu$ and real $q > 1$, the $q$-Rényi divergence is defined as

$$\mathsf R_q(\mu\Vert\nu) = \frac{1}{q-1} \log\Bigl( \mathbb E_\nu\Big[(\frac{d\mu}{d\nu})^q \Big] \Bigr).$$

Setting $q$0, this simplifies to

$q$1

Rényi divergences interpolate between various risk-sensitive divergences and play a central role in non-asymptotic information-theoretic analysis. In the context of interacting diffusions, propagation of chaos in Rényi divergence precisely captures the rate at which empirical marginals approach statistical independence as $q$2.

2. Main Theorem: Sharp Rényi Propagation of Chaos Rate

For an $q$3-particle interacting diffusion system at stationarity, let $q$4 denote the $q$5-marginal of the stationary law and $q$6 the mean-field Gibbs minimizer. The main result asserts that, under specified regularity and weak interaction assumptions, there exist constants $q$7—independent of $q$8—such that

$q$9

where hidden constants depend only on the log-Sobolev constant $O(d q^2/N^2)$0 and interaction smoothness parameter $O(d q^2/N^2)$1, and $O(d q^2/N^2)$2 notation suppresses polylogarithmic terms. The result establishes an optimal $O(d q^2/N^2)$3 rate for first marginals; the arguments also yield bounds for general $O(d q^2/N^2)$4-marginals of the form $O(d q^2/N^2)$5 under analogous conditions.

3. Technical Conditions and Setting

The result holds for stationary solutions of interacting diffusion systems characterized by:

• Potentials and Interactions: Confinement potential $O(d q^2/N^2)$6 and interaction kernel $O(d q^2/N^2)$7, with mean-field Gibbs law $O(d q^2/N^2)$8.
• Assumption 2.1 (Smoothness): $O(d q^2/N^2)$9 for all $q$0 (interaction gradient Lipschitz with parameter $q$1).
• Assumption 2.2 (Uniform log-Sobolev/Isoperimetry): There exists $q$2 such that for all $q$3,

$q$4

and this holds uniformly for all conditional measures of $q$5.

• Assumption 2.3 (Very Weak Interaction): $q$6; i.e., the interaction is a small perturbation on the log-Sobolev scale.

These assumptions guarantee uniqueness of the relevant Gibbs measures and well-posedness of both McKean–Vlasov and $q$7-particle stationary dynamics.

4. Proof Structure and Key Components

The analysis proceeds through several layers:

1. LSI-to-Rényi Inequality: For any $q$8 satisfying a log-Sobolev inequality, Lemma 3.1 shows

$q$9

where the Rényi–Fisher information is

$\mu,\nu$0

2. Donsker–Varadhan Splitting: The Rényi–Fisher information can be decomposed as

$\mu,\nu$1

with

$\mu,\nu$2

where $\mu,\nu$3 is a tilted law and $\mu,\nu$4 captures variance-type fluctuation.

3. Propagation-of-Chaos in Fisher Information: Lemma 3.4 yields

$\mu,\nu$5

under stationarity.

4. Control of Tilted-KL Term: Applying a second LSI-to-Rényi type argument:

$\mu,\nu$6

absorbing half of the $\mu,\nu$7 term to the left.

5. Exponential Moment Bound: Via concentration of Lipschitz functionals under the stationary law and LSI,

$\mu,\nu$8

with a hierarchical coupling controlling the relevant Lipschitz constant as $\mu,\nu$9.

Aggregating these, one concludes the sharp

$\mathbb{R}^d$0

for general $\mathbb{R}^d$1-marginals, so for $\mathbb{R}^d$2,

$\mathbb{R}^d$3

5. Dependence on System Parameters and Sharpness

The established bound reveals several sharp regimes:

• Scaling Laws: The bound scales as $\mathbb{R}^d$4 up to polylog factors. In the Gaussian case, the $\mathbb{R}^d$5 and $\mathbb{R}^d$6 dependences are sharp, with an explicit expansion yielding $\mathbb{R}^d$7.
• Optimality: The only apparent suboptimality is the $\mathbb{R}^d$8 dependence; the true scaling is conjectured to be linear in $\mathbb{R}^d$9.
• Reduction to KL: In the limit $\mu\ll\nu$0, the bound recovers the $\mu\ll\nu$1 KL divergence propagation-of-chaos law.
• Concentration Transfer: The corollary

$\mu\ll\nu$2

for $\mu\ll\nu$3 demonstrates high-probability transfer of concentration phenomena.

The stationary setting is essential for (i) closed-form score functions, (ii) uniform application of LSI to all conditional marginals, and (iii) sub-Gaussian concentration under the stationary law. Extensions to time-dependent (non-stationary) dynamics remain open.

6. Illustrative Gaussian Example

Consider $\mu\ll\nu$4, $\mu\ll\nu$5. The stationary law $\mu\ll\nu$6 is Gaussian with explicit block-covariance. For $\mu\ll\nu$7,

$\mu\ll\nu$8

For $\mu\ll\nu$9, this results in $q > 1$0, certifying the optimality of the $q > 1$1 and $q > 1$2 dependence, with the only discrepancy in the $q > 1$3 exponent relative to the general bound.

7. Broader Context and Implications

Propagation of chaos in Rényi divergence generalizes classical mean-field limit results beyond KL or Wasserstein metrics, providing a fine-grained, risk-sensitive quantification of independence in high-dimensional particle systems. The entropic and functional-analytic approach employed relies crucially on log-Sobolev inequalities, Donsker–Varadhan large deviations, and hierarchical couplings, synthesizing ideas from kinetic theory, concentration of measure, and information theory. These results further clarify the rates and mechanisms underlying empirical measure convergence in statistical physics and related domains (Zhang, 15 Jan 2026). High-dimensional sharpness and explicit scaling laws lay the foundation for subsequent investigations into temporal dynamics and non-stationary extensions.