Stochastic Interpolants: A Unifying Framework for Flows and Diffusions

Michael S. Albergo, Nicholas M. Boffi, Eric Vanden-Eijnden

2023 · arXiv

Stochastic Interpolants: A Unifying Framework for Flows and Diffusions

Problem

Framing

Flow matching and score-based diffusion gave separate recipes for ODE and SDE generative transport, with different objectives and likelihood guarantees. This paper unifies them with stochastic interpolants: design any finite-time bridge ρ(t)\rho(t) between ρ0\rho_0 and ρ1\rho_1, then sample it with either a probability-flow ODE or an SDE using quadratic regression losses.

Currently Used Methods

Foundational

Proposed Method

Architecture

The core object is the interpolant

xt=I(t,x0,x1)+γ(t)z,t[0,1],  zN(0,Id).\mathbf{x}_t = I(t, \mathbf{x}_0, \mathbf{x}_1) + \gamma(t)\mathbf{z}, \qquad t \in [0,1],\; \mathbf{z} \sim \mathcal{N}(\mathbf{0}, \mathbf{I}_d).

It separates bridge design from sampler choice. The same learned velocity and score define either deterministic probability flow or stochastic dynamics.

Verified figure: the framework bridges \rho_0 to \rho_1 through a shared density path \rho(t), then samples it with either smooth ODE trajectories or noisy SDE trajectories.

Loss / Objective

Velocity and score are learned by least-squares regression on interpolant samples.

Lb[b^]=01E[12b^(t,xt)2(tI(t,x0,x1)+γ˙(t)z)b^(t,xt)]dtL_b[\hat{\mathbf{b}}] = \int_0^1 \mathbb{E}\left[ \frac{1}{2}\|\hat{\mathbf{b}}(t, \mathbf{x}_t)\|^2 - \left(\partial_t I(t, \mathbf{x}_0, \mathbf{x}_1) + \dot{\gamma}(t)\mathbf{z}\right) \cdot \hat{\mathbf{b}}(t, \mathbf{x}_t) \right] dt Ls[s^]=01E[12s^(t,xt)2+γ1(t)zs^(t,xt)]dtL_s[\hat{\mathbf{s}}] = \int_0^1 \mathbb{E}\left[ \frac{1}{2}\|\hat{\mathbf{s}}(t, \mathbf{x}_t)\|^2 + \gamma^{-1}(t)\mathbf{z} \cdot \hat{\mathbf{s}}(t, \mathbf{x}_t) \right] dt

Sampling Rule / Algorithm

The same bridge density admits both a probability-flow ODE and a forward SDE with tunable diffusion ϵ(t)\epsilon(t).

dXtdt=b(t,Xt)\frac{d\mathbf{X}_t}{dt} = \mathbf{b}(t, \mathbf{X}_t) dXt=(b(t,Xt)+ϵ(t)s(t,Xt))dt+2ϵ(t)dWtd\mathbf{X}_t = \left(\mathbf{b}(t, \mathbf{X}_t) + \epsilon(t)\mathbf{s}(t, \mathbf{X}_t)\right)dt + \sqrt{2\epsilon(t)}\, d\mathbf{W}_t

Training Procedure

Evaluation

Datasets

Metrics

Headline results

Ablations

Verified figure: grid of checkerboard density estimates comparing exact samples against models across several \gamma(t) schedules and diffusion levels \epsilon; SDE columns look closer to the exact target than the \epsilon=0 probability-flow cases.

Method Strengths and Weaknesses

Strengths

Weaknesses

Suggestions from the authors

Links

Prior Papers

Further Papers