Papers
Topics
Authors
Recent
Search
2000 character limit reached

Kim–Milman Construction

Updated 9 July 2026
  • Kim–Milman Construction is a deterministic transport-map that reverses an OU diffusion to map the Gaussian measure onto a target measure using a score-based ODE.
  • It leverages semiconcavity conditions and Lipschitz regularity, enabling precise control over functional inequalities and stability estimates via Fisher information.
  • The approach extends to manifold settings and offers sharper stability bounds compared to classical optimal transport methods like Brenier maps.

The Kim–Milman construction is a deterministic transport-map construction on Rd\mathbb{R}^d that sends the standard Gaussian reference measure γ=N(0,Id)\gamma=\mathcal{N}(0,I_d) to a target probability measure μ\mu by reversing the Ornstein–Uhlenbeck (OU) diffusion. In recent machine-learning literature, the same object is often called the probability flow ODE associated with a diffusion model. In the notation used in recent work, the resulting map is denoted either TKMμT_{\mathrm{KM}^\mu} or Tμ\mathsf T^\mu, and it satisfies (TKMμ)#γ=μ(T_{\mathrm{KM}^\mu})_\#\gamma=\mu or (Tμ)♯γ=μ(\mathsf T^\mu)_\sharp\gamma=\mu (Chewi et al., 3 Nov 2025, Chewi et al., 22 Jun 2026).

1. Dynamical definition through Ornstein–Uhlenbeck reversal

The standard setup fixes the ambient space Rd\mathbb{R}^d, the Gaussian reference

γ=N(0,Id),\gamma=\mathcal{N}(0,I_d),

and a target measure μ\mu with finite second moment, or in the regular setting a smooth Lebesgue density. The forward OU process is initialized at the target: γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)0 If γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)1, then γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)2, where γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)3 is the OU semigroup, and γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)4 as γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)5 (Chewi et al., 3 Nov 2025).

On a finite horizon γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)6, reversing time in the Fokker–Planck equation yields a continuity equation for the backward-time density γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)7: γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)8 This is the continuity equation for the deterministic ODE

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)9

Assuming μ\mu0, the drift can also be written in semigroup form: μ\mu1 This ODE is the Kim–Milman reverse heat flow ODE, and its flow map μ\mu2 satisfies

μ\mu3

In particular, μ\mu4 (Chewi et al., 22 Jun 2026).

The construction is therefore not variational: unlike Brenier transport, it is defined by reversing a diffusion and integrating an ODE whose drift is the score of the OU-smoothed target. This direct dependence on μ\mu5 is central to both its regularity theory and its stability theory.

2. Infinite-horizon limit and existence of the map

Letting μ\mu6 sends the initial law μ\mu7 to μ\mu8. In the regular regime, the Kim–Milman map is defined as the infinite-horizon limit

μ\mu9

and for each TKMμT_{\mathrm{KM}^\mu}0 one also obtains a stopped map

TKMμT_{\mathrm{KM}^\mu}1

For regular targets, a key representation is the integral equation

TKMμT_{\mathrm{KM}^\mu}2

and differentiation in TKMμT_{\mathrm{KM}^\mu}3 gives

TKMμT_{\mathrm{KM}^\mu}4

The Kim–Milman map is then TKMμT_{\mathrm{KM}^\mu}5 (Chewi et al., 22 Jun 2026).

A major recent development is a general existence theorem for arbitrary TKMμT_{\mathrm{KM}^\mu}6. A regular approximating sequence TKMμT_{\mathrm{KM}^\mu}7 is defined by requiring that every TKMμT_{\mathrm{KM}^\mu}8 satisfies the regularity condition TKMμT_{\mathrm{KM}^\mu}9 and that Tμ\mathsf T^\mu0. The paper gives the example: one can truncate Tμ\mathsf T^\mu1 to a large ball and convolve with a small Gaussian. For such an approximation, the limits

Tμ\mathsf T^\mu2

exist in Tμ\mathsf T^\mu3, are independent of the approximating sequence, and satisfy

Tμ\mathsf T^\mu4

Thus, for every Tμ\mathsf T^\mu5, the Kim–Milman map exists as an Tμ\mathsf T^\mu6 limit, without requiring Tμ\mathsf T^\mu7 on Tμ\mathsf T^\mu8 (Chewi et al., 22 Jun 2026).

This general existence theorem separates the questions of definition and quantitative control. Existence holds on all of Tμ\mathsf T^\mu9, while Lipschitz regularity and sharpened stability estimates still require structural assumptions on at least one target.

3. Semiconcavity along the OU semigroup and Lipschitz regularity

The main regularity hypothesis used in the recent theory is the semiconcavity condition

(TKMμ)#γ=μ(T_{\mathrm{KM}^\mu})_\#\gamma=\mu0

Here (TKMμ)#γ=μ(T_{\mathrm{KM}^\mu})_\#\gamma=\mu1 is a scalar control on the Hessian of the OU-smoothed log-density. Under (TKMμ)#γ=μ(T_{\mathrm{KM}^\mu})_\#\gamma=\mu2, for all (TKMμ)#γ=μ(T_{\mathrm{KM}^\mu})_\#\gamma=\mu3, the stopped flow map (TKMμ)#γ=μ(T_{\mathrm{KM}^\mu})_\#\gamma=\mu4 exists and is Lipschitz with constant

(TKMμ)#γ=μ(T_{\mathrm{KM}^\mu})_\#\gamma=\mu5

In particular, the limiting map (TKMμ)#γ=μ(T_{\mathrm{KM}^\mu})_\#\gamma=\mu6 is Lipschitz with constant (TKMμ)#γ=μ(T_{\mathrm{KM}^\mu})_\#\gamma=\mu7. Moreover, the OU marginals satisfy the log-Sobolev inequality

(TKMμ)#γ=μ(T_{\mathrm{KM}^\mu})_\#\gamma=\mu8

(Chewi et al., 22 Jun 2026).

Several classes of targets are known to satisfy (TKMμ)#γ=μ(T_{\mathrm{KM}^\mu})_\#\gamma=\mu9. If (Tμ)♯γ=μ(\mathsf T^\mu)_\sharp\gamma=\mu0 with (Tμ)♯γ=μ(\mathsf T^\mu)_\sharp\gamma=\mu1 (Tμ)♯γ=μ(\mathsf T^\mu)_\sharp\gamma=\mu2-strongly convex, then

(Tμ)♯γ=μ(\mathsf T^\mu)_\sharp\gamma=\mu3

and (Tμ)♯γ=μ(\mathsf T^\mu)_\sharp\gamma=\mu4. This is the original Kim–Milman setting and yields (Tμ)♯γ=μ(\mathsf T^\mu)_\sharp\gamma=\mu5. The theory also covers strongly convex plus Lipschitz perturbation and more general asymptotically convex potentials characterized by a convexity profile (Tμ)♯γ=μ(\mathsf T^\mu)_\sharp\gamma=\mu6; in these cases, precise formulas for (Tμ)♯γ=μ(\mathsf T^\mu)_\sharp\gamma=\mu7 are available and satisfy the required integrability and boundedness properties (Chewi et al., 22 Jun 2026).

In the stability paper, the same structural condition is written as

(Tμ)♯γ=μ(\mathsf T^\mu)_\sharp\gamma=\mu8

and explicit (Tμ)♯γ=μ(\mathsf T^\mu)_\sharp\gamma=\mu9 are developed for three regimes: strong log-concavity, strong log-concavity plus log-Lipschitz perturbations, and asymptotically positive convex profiles (Chewi et al., 3 Nov 2025).

These bounds are curvature-type conditions along the OU flow. They control the Lipschitz modulus of the reverse-time drift and are the basic input for both Grönwall estimates on trajectories and functional inequalities on the OU marginals.

4. Fisher-information stability and the score-based viewpoint

A central recent result is that stability of the Kim–Milman flow map is controlled not by Wasserstein distance, but by relative Fisher information. For Rd\mathbb{R}^d0, the relative Fisher information is

Rd\mathbb{R}^d1

and the Rd\mathbb{R}^d2-version is

Rd\mathbb{R}^d3

If Rd\mathbb{R}^d4 and Rd\mathbb{R}^d5 denote the Kim–Milman maps associated with two targets Rd\mathbb{R}^d6 and Rd\mathbb{R}^d7, then under suitable regularity assumptions on Rd\mathbb{R}^d8,

Rd\mathbb{R}^d9

where γ=N(0,Id),\gamma=\mathcal{N}(0,I_d),0 is an explicit constant depending on curvature-type bounds associated with γ=N(0,Id),\gamma=\mathcal{N}(0,I_d),1. In important special cases such as strong log-concavity, γ=N(0,Id),\gamma=\mathcal{N}(0,I_d),2 takes the simple form γ=N(0,Id),\gamma=\mathcal{N}(0,I_d),3. The paper also proves

γ=N(0,Id),\gamma=\mathcal{N}(0,I_d),4

under additional assumptions (Chewi et al., 3 Nov 2025).

At finite horizon, the core estimate takes the form

γ=N(0,Id),\gamma=\mathcal{N}(0,I_d),5

with

γ=N(0,Id),\gamma=\mathcal{N}(0,I_d),6

The proof couples the two reverse ODEs, differentiates γ=N(0,Id),\gamma=\mathcal{N}(0,I_d),7, controls the drift mismatch by Cauchy–Schwarz, and propagates the score difference through the OU semigroup via a Fisher-information decay estimate due to Wibisono. In the notation of the paper,

γ=N(0,Id),\gamma=\mathcal{N}(0,I_d),8

This is combined with the Hessian bound γ=N(0,Id),\gamma=\mathcal{N}(0,I_d),9 to close the Grönwall argument (Chewi et al., 3 Nov 2025).

The resulting viewpoint is structurally adapted to the construction itself. The reverse-time drift is literally a score field, and the right-hand side of the stability estimate measures the μ\mu0-size of score differences. This is why the theory is formulated in terms of relative Fisher information rather than directly in terms of transport cost.

5. Entropy stability, near-Lipschitz μ\mu1 stability, and comparison with optimal transport

A second line of recent results proves stability in relative entropy. If μ\mu2 satisfies μ\mu3, then for any μ\mu4,

μ\mu5

This yields the sharpened Talagrand-type inequality

μ\mu6

(Chewi et al., 22 Jun 2026).

The same paper proves near-Lipschitz stability in μ\mu7, up to a logarithmic factor. Assume μ\mu8 satisfies μ\mu9, both γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)00 have second moments bounded by γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)01, and there exists γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)02 such that

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)03

Then

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)04

For small γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)05, the dominant behavior is essentially linear in γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)06, multiplied by a slowly varying logarithmic factor. The proof combines the differential inequality

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)07

the de Bruijn identity

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)08

and the log-Harnack estimate

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)09

Because the resulting Fisher-information control behaves like γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)10 near γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)11, the argument uses an “early stopping” parameter γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)12 and then optimizes over γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)13, which produces the logarithmic correction (Chewi et al., 22 Jun 2026).

These results place the construction in a distinctive position relative to classical transport maps. Brenier maps are gradients of convex potentials and, in general, are only Hölder stable in γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)14, with exponent up to γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)15, under variations of the target. Entropic transport maps enjoy γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)16-Lipschitz stability with constants depending on the regularization parameter. The Kim–Milman map is not defined by a variational problem; its stability theory is instead organized around score propagation under the OU semigroup and the functional inequalities generated by the semiconcavity condition (Chewi et al., 22 Jun 2026, Chewi et al., 3 Nov 2025).

6. Examples, corollaries, and manifold extensions

For strongly log-concave targets, the Fisher-information estimate becomes especially simple: γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)17 Combined with the HWI-type bound

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)18

this gives

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)19

For log-Lipschitz perturbations of strongly log-concave measures,

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)20

with γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)21 γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)22-strongly convex and γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)23 smooth and γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)24-Lipschitz, the paper gives the explicit bound

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)25

An γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)26-analogue leads to

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)27

with explicit γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)28 in the same regimes (Chewi et al., 3 Nov 2025).

The construction also admits a manifold version coupled to geometric evolution. In the setting of a smooth flow of Riemannian metrics γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)29 on a manifold γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)30, preserving volume and satisfying

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)31

if there exist γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)32 such that

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)33

then the Kim–Milman transport map

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)34

defined as the limit of the inverse characteristic flow is γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)35-Lipschitz. Specializing to rescaled Ricci flow on a surface gives the criterion

γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)36

where γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)37 is the curvature potential defined by γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)38. Under a γ=N(0,Id)\gamma=\mathcal{N}(0,I_d)39 near-roundness hypothesis, this yields a contraction from the unit round sphere to a nearly spherical positively curved surface, pushing forward the rescaled round volume onto the target volume (Serres, 19 Aug 2025).

This manifold extension shows that the construction is not confined to the Gaussian–OU setting on Euclidean space. A plausible implication is that the defining mechanism—transport generated by a continuity equation with controlled Hessian—can be coupled to geometric flows whenever one has a suitable Bakry–Émery-type differential inequality.

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to Kim-Milman Construction.