---
title: Benamou–Brenier Formulation in Optimal Transport
url: https://www.emergentmind.com/topics/benamou-brenier-formulation
type: topic
---

# Benamou–Brenier Formulation in Optimal Transport

The Benamou–Brenier formulation is a dynamic (Eulerian) perspective on optimal transport that expresses the cost of transporting one measure into another as the minimum kinetic energy of admissible evolutions, subject to a continuity equation linking the time-dependent densities and velocities. This framework complements and generalizes the classical (static) Monge–Kantorovich problem, underpins the geometry of Wasserstein spaces, and has been significantly extended to nonstandard settings including martingale transport, unbalanced transport, sub-Riemannian manifolds, mean-field games, and stochastic and quantum contexts.

## 1. Classical Benamou–Brenier Formulation: Definition and Structure

The Benamou–Brenier formulation expresses, for $p>1$, the $p$-Wasserstein distance between two probability measures $\mu_0$ and $\mu_1$ as the minimal action over evolving pairs $(\rho_t, v_t)$ governed by the continuity equation:
\[
W_p^p(\mu_0, \mu_1) = \inf \left\{ \int_0^1 \int_{\mathbb{R}^d} |v_t(x)|^p d\rho_t(x) dt :
\partial_t \rho_t + \nabla \cdot (\rho_t v_t) = 0,\,\, \rho_0 = \mu_0,\,\, \rho_1 = \mu_1 \right\}
\]
For $p=2$, the "kinetic energy" interpretation of the integrand appears and the formulation admits a formal Riemannian structure on the space of probability measures with finite second moment [2402.04842]. The variables $\rho_t$ and $v_t$ represent the intermediate densities and velocities, respectively, under constraint of mass conservation.

This dynamic (Eulerian) view contrasts with the classical Monge problem, which is Lagrangian. The equivalence of the Kantorovich and Benamou–Brenier formulations in Euclidean and (under suitable regularity/controllability assumptions) on manifolds underscores the foundational role of this perspective [2507.20959, 2407.16088].

## 2. Methodological Extensions: Unbalanced, Martingale, and Control-Based Formulations

### 2.1. Unbalanced and Generalized Mass Transport

To address situations where the source and target measures are not equally normalized, the Benamou–Brenier formulation is extended by relaxing mass conservation via an additional source (or growth) term $\alpha$ in the continuity equation:
\[
\partial_t \rho + \nabla \cdot (\rho v) = \alpha \rho
\]
The dynamic formulation for the Wasserstein–Fisher–Rao (WFR) metric or Hellinger–Kantorovich distance thus becomes
\[
(WFR(\rho_0, \rho_1))^2 = \inf_{(\rho, v, \alpha)} \frac{1}{2} \int_0^1 \int_M [a^2 |v(x, t)|^2 + b^2 \alpha(x, t)^2]\, d\rho_t(x) dt
\]
subject to the above constraint and prescribed endpoints. This framework allows for creation or annihilation of mass during transport and reduces to the classical setting when $\alpha=0$ [2112.11056, 2008.02909, 1902.03367].

Related genetic-mass transport problems can be encoded as vector-valued Benamou–Brenier formulations, where source terms correspond to fluxes between layers ("channels") in a multi-component system [2008.02909].

### 2.2. Martingale Optimal Transport

The Benamou–Brenier philosophy extends to stochastic settings, particularly for martingale optimal transport (MOT). Here, the goal is to minimize an action over admissible martingales interpolating two laws in convex order, subject to the martingale property, resulting in a continuous-time dynamic cost functional—typically via Fokker–Planck dynamics [1707.01493, 1708.04869, 2406.04016]. The dynamic representation enables precise geodesic and interpolation structures in stochastic transport, and explicit optimizer characterization in terms of stretched Brownian or geometric Bass martingales connects the theory to problems in finance and stochastic control.

### 2.3. Nonlinear Control-Affine and Sub-Riemannian Structures

For kinetic costs adapted to general controlled dynamics—e.g., $\dot \omega(t) = f_0(\omega(t)) + \sum_{i=1}^n u_i(t) f_i(\omega(t))$—the Benamou–Brenier formulation, subject to a generalized continuity equation, provides a convex or relaxable dynamic formulation for the transport cost [2407.16088, 2507.20959]. On sub-Riemannian manifolds without abnormal geodesics, the formulation in the subspace of horizontal (admissible) velocities combined with a suitable Young measure relaxation is shown to be equivalent to the static (Kantorovich) problem, ensuring existence and consistency of minimizers [2507.20959].

## 3. Analytical and Structural Properties

The Benamou–Brenier action minimization is convex in $(\rho, v)$ for $p=2$. Its dual formulations often link to Hamilton–Jacobi–Bellman (HJB) equations in both classical and quantum contexts [2104.11923], providing insight into geodesicity, convexity, and regularity properties.

A significant result in regularity theory is the variational approach using the Benamou–Brenier formulation to obtain $\epsilon$-regularity and $C^{1,\alpha}$ partial regularity of optimal transport maps, relying on precise comparison between the minimal kinetic energy flows and interpolating harmonic functions [1704.05339]. The dynamic nature facilitates Campanato-type excess decay, underlining the connection between minimal energy and local regularity without invoking the maximum principle.

In extensions involving additional phenomena (e.g., bulk-interface transfer [2002.07724], diffusion [2201.07307], entropic effects [2104.11923], or mean-field game coupling [2505.12097]), the Benamou–Brenier structure enables both new analytical results and meaningful geometric interpolation between disparate transport regimes.

## 4. Numerical and Algorithmic Developments

The dynamic (Eulerian) structure of the Benamou–Brenier formulation is highly amenable to numerical discretization. Practical computational strategies include finite volume schemes for gradient flows [1907.08305], space-time discretizations for convergence guarantees even on manifolds [1909.08790], and the use of advanced preconditioners for saddle-point systems arising from discretized convex optimization formulations [2209.00315].

For unbalanced and generalized OT, primal-dual and proximal algorithms are facilitated by explicit convexity and separability properties in the Benamou–Brenier framework [1902.03367, 2505.12097]. Proximal OT divergences, viewed as an infimal convolution of OT cost and information divergence, yield mean-field game structures where optimality conditions are governed by coupled backward HJB and forward continuity PDEs [2505.12097].

The Eulerian approach also underlies state-of-the-art methodologies for spatiotemporal imaging, motion estimation, and tomographic reconstruction, where diffeomorphic constraints and optimal transport costs are naturally encoded and efficiently optimized [2011.11906].

## 5. Variants, Dualities, and Generalizations

The duality theory associated with the Benamou–Brenier formulation recurs across numerous extensions. In the noncommutative (quantum) transport setting, a dual formula emerges with Hamilton–Jacobi–Bellman subsolutions as noncommutative analogues of classical test functions, confirming a parallel to the classical Wasserstein geometry at the level of density matrices [2104.11923].

A significant geometric insight is that on closed manifolds (or suitable cones), the equivalence between dynamic (Benamou–Brenier-type) and static (semi-coupling, cone lift) OT formulations can be rigorously established, with the dynamic approach supporting polar factorization theorems and transfer of regularity conditions such as the Ma–Trudinger–Wang property [2112.11056]. On metric graphs and random measure spaces, weak formulations via regularized continuity equations and Palm probability techniques extend the Benamou–Brenier framework to singular, stochastic, or infinite-mass settings [2105.05677, 2402.04842].

## 6. Applications and Impact Across Disciplines

The Benamou–Brenier formulation is influential in analysis, geometry, statistical physics, machine learning, stochastic processes, and quantum information. Key applications include:
- Wasserstein gradient flows and diffusion equations on complex domains [2105.05677, 1907.08305]
- Martingale and geometric transport for financial modeling [1708.04869, 2406.04016]
- Bulk-interface coupled systems and diffusion-regularized transport in computational neurosciences [2002.07724, 2201.07307]
- Proximal optimization and generative modeling in probabilistic machine learning [2505.12097]
- Control-theoretic and sub-Riemannian transport [2407.16088, 2507.20959]
- Noncommutative and quantum transport [2104.11923]

The dynamic action-minimization paradigm enables geometric, analytic, and computational approaches that are flexible to generalizations, robust to singularities, and compatible with domain-specific constraints, supporting both theoretical insight and algorithmic tractability.

## 7. Limitations, Open Problems, and Future Directions

A principal limitation in the general theory arises in situations with abnormal geodesics or when compactness fails. For general sub-Riemannian manifolds with abnormal minimizers, the equivalence of dynamic and static optimal transport is nontrivial and remains an active area of research [2507.20959]. The robustness of dynamic and dual formulations under weak regularity, nonconvex costs, or measure-valued controls also presents technical challenges, requiring further development of relaxation and compactness arguments.

Ongoing extensions include:
- Couplings with mean-field game systems and learning in high-dimensional or singular spaces [2505.12097]
- Further exploration of unbalanced and entropy-transport settings, transferring regularity and geometric properties from static to dynamic formulations [2112.11056]
- Synergizing with data-driven and stochastic approaches, including within-network optimization and dynamic generative models

The continued unification and generalization of the Benamou–Brenier formulation with stochastic, quantum, unbalanced, and control-theoretic optimal transport is likely to drive new developments in analysis, geometry, and their applications.

Source: https://www.emergentmind.com/topics/benamou-brenier-formulation