Teacher-free one-step generative models such as W-Flow learn directly from distributional objectives, with the Sinkhorn divergence providing the strongest empirical performance among the energies studied in W-Flow. However, Sinkhorn divergence is based on balanced optimal transport, whose exact marginal constraints can induce undesirable transport under finite minibatch fluctuations. We introduce UW-Flow, which replaces the balanced Sinkhorn energy with a debiased unbalanced Sinkhorn divergence and derives the corresponding Wasserstein velocity. The resulting dynamics retain W-Flow's cross-minus-self structure, recover the balanced limit, and preserve target identifiability under suitable assumptions. We further characterize the finite-batch velocity error of balanced and unbalanced transport by comparing each minibatch field with its population counterpart. Our analysis identifies conditions under which marginal relaxation yields smaller minibatch-to-population velocity error, and shows that for separated matched modes, balanced sensitivity grows with inter-mode separation while the corresponding Unbalanced Optimal Transport (UOT) sensitivity remains controlled by local transport geometry. Controlled experiments on fixed-support and continuous two-mode distributions support these predictions. On MNIST, CIFAR-10, and ImageNet-100, UW-Flow improves generation under matched training budgets, with the largest gains in unconditional generation and consistent improvements in the class-conditional setting.
Figures & tables
Figure 1
Figure 2: Minibatch velocity error on fixed supports. Each method is evaluated against its own population velocity. Panel (a) sweeps B at d=4 , while panel (b) sweeps d at B=16 . Within each panel, the left plot shows B×MSE and the right plot shows B×RelMSE . RelMSE is MSE divided by the corresponding population-field energy.
Table 3
Figure 4: CIFAR-10 training dynamics across transport support sizes. FID versus training epoch for W-Flow and UW-Flow with τ∈{0.1,1,10} . FID is computed using 10,000 generated samples at each evaluation. Panels show B=4 , 16 , and 64 , with the sample budget per optimizer update fixed at Bopt=NcB=128 .
Setting
Method
FID ↓
IS ↑
Prec. ↑
Rec. ↑
Uncond
W-Flow
55.25
23.89
0.42
0.14
Uncond
UW-Flow
46.17
27.71
0.40
0.29
Cond
W-Flow
9.43
60.31
0.74
0.26
Cond
UW-Flow
9.11
61.34
0.75
0.26
Table 4: ImageNet-100 results at matched budgets.
Appendix figures & tables11 assets
Supplementary material from the paper’s appendix.
Appendix
Figure 5: Full fixed-support velocity-error results. Columns show B×MSE , B×RelMSE , and cosine similarity. The top and bottom rows vary B and d , respectively.
Figure 6: Velocity errors under Gaussian-mixture sampling. Columns show B×MSE , B×RelMSE , and cosine similarity. The top and bottom rows vary B and d , respectively.
Algorithm 2 Unbalanced Sinkhorn 1: Particles X={xi}i=1N , Y={yj}j=1M ; ε,τ,L . 2: ai←1/N,bj←1/M▹ Weights 3: Cij←21∥xi−yj∥22▹ Cost 4: ϕi←0,ψj←0▹ Initialize 5: α←τ+ετ▹ Relaxation 6: for ℓ=1,…,L do 7: ϕi←−αεlog∑jbje(ψj−Cij)/ε▹ Source 8: ψj←−αεlog∑iaie(ϕi−Cij)/ε▹ Target 9: end for 10: Πij←aibje(ϕi+ψj−Cij)/ε▹ Coupling 11: ri←∑jΠij▹ Row mass 12: si←airi▹ Mass ratio 13: Ti←ri1∑jΠijyj▹ Barycenter 14: return (T,s)
Appendix
Table 8
Nc
Ngen
Npos
Method
FID ↓
Prec. ↑
Rec. ↑
16
4
4
W-Flow
19.70
0.536
0.059
UW-Flow ( τ=10 )
22.54
0.562
0.081
UW-Flow ( τ=1 )
17.98
0.464
0.116
UW-Flow ( τ=0.1 )
18.66
0.419
0.253
8
8
8
W-Flow
13.69
0.461
0.239
UW-Flow ( τ=10 )
13.50
0.545
0.110
Appendix
Table 5: Full results of the fixed-budget minibatch discrepancy ablation on unconditional MNIST, with fixed Bopt=NcNgen=64 .
Nc
Ngen
Npos
Method
FID ↓
IS ↑
Prec. ↑
Rec. ↑
32
4
4
W-Flow
30.58
9.75
0.55
0.01
UW-Flow ( τ=10 )
22.03
9.58
0.57
0.02
UW-Flow ( τ=1 )
19.41
10.10
0.56
0.02
UW-Flow ( τ=0.1 )
15.56
10.60
0.52
0.07
16
8
8
W-Flow
24.39
8.99
0.58
0.01
UW-Flow ( τ=10 )
14.64
10.39
0.59
0.04
Appendix
Table 6: Full results of the fixed-budget minibatch discrepancy ablation on unconditional CIFAR-10, with fixed Bopt=NcNgen=128 .
Hyperparameter
MNIST
CIFAR-10
ImageNet-100
Model and representation
Generator
DiT-based, patch 2
DiT-based, patch 2
DiT-B/2
Training shape / space
28×28×1
32×32×3
32×32×4
Feature encoder (frozen)
ConvNeXt V2–Tiny
ConvNeXt V2–Tiny
Latent-MAE ResNet (width 256)
Feature normalization
Distance-based scaling
Distance-based scaling
Distance-based scaling
Velocity normalization
RMS normalization
RMS normalization
RMS normalization
Appendix
Table 7: Hyperparameters for image-generation experiments. Ngen and Npos are the numbers of generated and target samples per transport group; the total generated batch size is NcNgen . The listed ε and τ values are before scaling by the feature dimension D .
Setting
Nc
Ngen
Npos
Bopt
Method
FID ↓
IS ↑
Prec. ↑
Rec. ↑
Uncond
64
64
64
4096
W-Flow
55.25
23.89
0.42
0.14
UW-Flow ( τ=10 )
55.00
25.01
0.40
0.16
UW-Flow ( τ=1 )
51.08
25.47
0.41
0.19
UW-Flow ( τ=0.1 )
46.17
27.71
0.40
0.29
Cond
64
64
64
4096
W-Flow
9.43
60.31
0.74
0.26
UW-Flow ( τ=10 )
9.47
61.08
0.75
0.25
Appendix
Table 8: Full ImageNet-100 generation results following the ImageNet ablation setup and CFG selection protocol of W-Flow ( Han et al., 2026 ) . Both settings use Nc=64 , Ngen=Npos=64 , and Bopt=4096 . For class-conditional generation, we select the CFG scale with the lowest FID for each method from {1.0,1.5,2.0,2.5,3.0,3.5,4.0} .
Sinkhorn iterations L
L
FID ↓
IS ↑
Precision ↑
Recall ↑
5
5.04
10.39
0.57
0.41
10
4.73
10.54
0.58
0.45
20
5.16
10.26
0.59
0.40
30
4.51
10.37
0.60
0.45
40
5.18
10.29
0.57
0.45
Appendix
Table 9: Sinkhorn parameter ablations on CIFAR-10. We fix Nc=2 , Ngen=Npos=64 , Bopt=128 , and τ=0.1 . The default configuration uses L=30 and ε=0.05 ; each table varies one parameter while fixing the other.
Figure 7: Qualitative B×τ ablation on checkerboard. Rows vary B ; columns show the target, UW-Flow with increasing τ , and W-Flow.
Figure 8: Qualitative B×τ ablation on two moons. Rows vary B ; columns show the target, UW-Flow with increasing τ , and W-Flow.
Figure 9: Qualitative B×τ ablation on Gaussian ring. Rows vary B ; columns show the target, UW-Flow with increasing τ , and W-Flow.
Drifting models enable one-step generation by amortizing distribution transport into training, but this efficiency places greater demands on the transport field estimated at each update. In large-scale training, the field is computed from finite mini-batches of generated and real samples, which provide only imperfect approximations to the underlying distributions. Balanced optimal transport enforces exact mass matching within every mini-batch, making the estimated field sensitive to the particular composition of the real-data batch. We find that generated and real samples should be treated asymmetrically: letting the mass assigned to real samples adapt while keeping every generated sample fully transported improves generation across six feature-space metrics in controlled ablations, and is more robust to the relaxation strength than relaxing both marginals simultaneously, which falls below balanced transport under stronger relaxation. Motivated by this observation, we propose Unbalanced Optimal Transport Gradient Flow (UOT-GF), which keeps the generated-sample marginal fixed and relaxes only the real-data marginal. Under identical settings at DiT-B/2 on ImageNet-256, UOT-GF improves Fréchet Inception Distance (FID) from 1.53 to 1.46 over the balanced W-Flow baseline; scaling the same recipe yields 1.34 and 1.22 FID at L/2 and XL/2, the best FID among the one-step models we compare. We further derive the induced UOT transport force, establish a kinetic Vlasov--Fokker--Planck formulation whose overdamped zero-temperature limit recovers the drifting dynamics, and characterize non-target stationary states together with sufficient conditions for convergence.
Yirong Shen, Mengfei Xia, Junpeng Jing +2
Imperial College London · Ant Group · Brunel University of London
Diffusion models and flow-based methods have shown impressive generative capability, especially for images, but their sampling is expensive because it requires many iterative updates. We introduce W-Flow, a framework for training a generator that transforms samples from a simple reference distribution into samples from a target data distribution in a single step. This is achieved in two steps: we first define an evolution from the reference distribution to the target distribution through a Wasserstein gradient flow that minimizes an energy functional; second, we train a static neural generator to compress this evolution into one-step generation. We instantiate the energy functional with the Sinkhorn divergence, which yields an efficient optimal-transport-based update rule that captures global distributional discrepancy and improves coverage of the target distribution. We further prove that the finite-sample training dynamics converge to the continuous-time distributional dynamics under suitable assumptions. Empirically, W-Flow sets a new state of the art for one-step ImageNet 256×256 generation, achieving 1.29 FID, with improved mode coverage and domain transfer. Compared to multi-step diffusion models with similar FID scores, our method yields approximately 100× faster sampling. These results show that Wasserstein gradient flows provide a principled and effective foundation for fast and high-fidelity generative modeling.
We introduce a novel Wasserstein-1 (W1) path-space divergence for stochastic and deterministic dynamics and establish a Wasserstein Uncertainty Propagation (WUP) theorem that bounds the W1 distance between terminal distributions by the proposed divergence, equivalently characterized by a weighted L2 discrepancy between the underlying drifts and the W1 distance between their initial measures. A key ingredient is a probabilistic framework combining adjoint Feynman-Kac representations with synchronous coupling (and reflection coupling on bounded domains), yielding Wasserstein stability estimates beyond existing PDE- and Girsanov-based approaches. The framework accommodates time-varying and possibly degenerate diffusion coefficients, empirical and singular measures, and remains valid in the deterministic limit of flow matching. Unlike KL-based uncertainty quantification bounds, it does not require absolute continuity of path measures and therefore remains well-defined in singular settings. As consequences of the WUP theorem, we derive W1 robustness and generalization bounds for score-based generative models and flow matching at both population and finite-sample levels. We further specialize the framework to group-symmetric targets, providing the first error analysis of equivariant flow-based models and the first quantitative comparison between data augmentation and equivariant inductive bias. Our analysis identifies a symmetry-aware Wasserstein path-space divergence that quantifies the model-form error induced by non-equivariant parametrizations. We prove that this error cannot be removed by additional data or training and vanishes only under equivariant architectures, establishing a precise theoretical advantage of equivariant inductive bias over data augmentation. Numerical experiments on group-symmetric Gaussian mixtures corroborate the theory.
Ziyu Chen, Markos A. Katsoulakis, Benjamin J. Zhang
Department of Mathematics and School of Data and Information Sciences, University of North Carolina at Chapel Hill · Department of Mathematics and Statistics, University of Massachusetts Amherst · Department of Mathematics, Rutgers University–New Brunswick