---
title: Multi-Marginal Monge Problem
url: https://www.emergentmind.com/topics/multi-marginal-monge-problem
type: topic
---

# Multi-Marginal Monge Problem

The multi-marginal Monge problem is the deterministic branch of multi-marginal optimal transport. Given marginals \(\mu_1,\dots,\mu_m\) on spaces \(X_1,\dots,X_m\) and a cost \(c:X_1\times\cdots\times X_m\to\mathbb{R}\), it asks for measurable maps \(T_i:X_1\to X_i\), \(i=2,\dots,m\), such that \((T_i)_\#\mu_1=\mu_i\) and the graph-induced plan \((\mathrm{Id},T_2,\dots,T_m)_\#\mu_1\) minimizes the total cost. Its relaxed counterpart is the multi-marginal Kantorovich problem, which minimizes over all couplings \(\gamma\in\Pi(\mu_1,\dots,\mu_m)\). The central structural question is when an optimal Kantorovich plan is concentrated on a graph, when it is unique, and when the best one can hope for is concentration on several graphs or on a higher-dimensional set [1406.0026][1307.6293].

## 1. Variational formulation and basic geometry

For Polish spaces or smooth manifolds \(X_i\) with prescribed marginals \(\mu_i\), the Kantorovich problem is
\[
\inf_{\gamma\in \Pi(\mu_1,\dots,\mu_m)} \int c(x_1,\dots,x_m)\,d\gamma,
\]
while the Monge formulation restricts admissible couplings to plans of the form
\[
\gamma=(\mathrm{Id},T_2,\dots,T_m)_\#\mu_1,
\qquad (T_i)_\#\mu_1=\mu_i.
\]
Thus the Monge problem is a constrained version of the Kantorovich problem: it minimizes over deterministic couplings, whereas Kantorovich allows arbitrary plans with the prescribed marginals [1601.05608][1007.0424].

A Monge solution exists precisely when an optimal Kantorovich plan is supported on a graph over the first marginal. This is straightforward in formal terms but highly nontrivial in structure, because for \(m\ge 3\) the geometry of \(\Pi(\mu_1,\dots,\mu_m)\) differs sharply from the two-marginal case. The literature repeatedly emphasizes that the general multi-marginal Monge problem is largely open for \(m\ge 3\), even though important positive classes are known, most notably the quadratic cost of Gangbo–Świȩch [1212.1680].

The same problem is often written in maximization form when one works with a surplus \(b\) instead of a cost \(c\). In that convention one maximizes \(\int b\,d\gamma\) over \(\Pi(\mu_1,\dots,\mu_m)\), and a Monge solution is again a graph-induced optimizer. This dual terminology is standard in graph-based and symmetric formulations [2104.09488].

## 2. Twist, splitting sets, and uniqueness of graphical optimizers

The decisive objects are \(c\)-splitting sets. A set \(S\subset X_1\times\cdots\times X_m\) is \(c\)-splitting if there exist functions \(u_i\) such that
\[
\sum_{i=1}^m u_i(x_i)\le c(x_1,\dots,x_m)
\]
everywhere, with equality on \(S\). Optimal plans are supported on such sets, so the Monge problem is governed not only by the global geometry of \(c\) but by the behavior of \(D_{x_1}c\) along splitting sets [1307.6293].

Kim–Pass introduced the condition that \(c\) be twisted on \(c\)-splitting sets: for fixed \(x_1\), the map
\[
(x_2,\dots,x_m)\mapsto D_{x_1}c(x_1,x_2,\dots,x_m)
\]
must be injective on each \(c\)-splitting set contained in \(\{x_1\}\times X_2\times\cdots\times X_m\). Under continuity and semi-concavity of \(c\), and absolute continuity of \(\mu_1\), this condition implies that the optimal plan is unique and concentrated on the graph of a measurable map, hence yields a unique Monge solution [1307.6293].

Pass had earlier obtained a differential criterion of the same flavor. In his theorem, \(c\in C^2\), \((1,m)\)-non-degeneracy, \((1,m)\)-twist, and negative definiteness of a global covariant \(2\)-tensor \(T\), together with the condition that \(\mu_1\) does not charge sets of Hausdorff dimension \(\le n-1\), imply that every Kantorovich optimizer is concentrated on a single graph. Consequently, both the Monge and Kantorovich problems have unique solutions [1007.0424].

Later work weakened the cost-side assumptions by strengthening the marginal regularity assumptions. Pass–Vargas-Jiménez introduced twist on \(c\)-splitting sets with respect to variables \(x_1,x_{k_1},\dots,x_{k_r}\): injectivity of \(D_{x_1}c\) is required only on refined contact sets where the corresponding dual potentials are differentiable. If \(\mu_1,\mu_{k_1},\dots,\mu_{k_r}\) are absolutely continuous, then the optimal plan is again unique and graph-concentrated. When only \(\mu_1\) is assumed absolutely continuous, this condition collapses back to the original twist on splitting sets [2202.06783].

## 3. Cyclical monotonicity, duality, and concentration on several graphs

In the multi-marginal setting, \(c\)-cyclical monotonicity requires that for any finite family \(\{(x_1^{(i)},\dots,x_d^{(i)})\}_{i=1}^n\) in a set \(\Gamma\) and any permutations \(\sigma_2,\dots,\sigma_d\),
\[
\sum_{i=1}^n c\bigl(x_1^{(i)},\dots,x_d^{(i)}\bigr)
\le
\sum_{i=1}^n c\bigl(x_1^{(i)},x_2^{(\sigma_2(i))},\dots,x_d^{(\sigma_d(i))}\bigr).
\]
This is the natural extension of the two-marginal notion: one fixes the first coordinate and independently reshuffles the remaining marginals [1601.05608].

Griessler proved that, under the assumptions that \(c\) is continuous, non-negative, and bounded above by a sum of integrable functions \(f_1+\cdots+f_d\), multi-marginal \(c\)-cyclical monotonicity is sufficient for optimality. The proof passes through \(c\)-splitting sets and dual potentials: every \(c\)-cyclically monotone set is shown to be \(c\)-splitting, and then the dual inequality yields optimality. In the Monge context, this turns \(c\)-cyclical monotonicity of a candidate graph into a direct sufficient criterion for Monge optimality [1601.05608].

A weaker conclusion than a single graph is often the correct one. Moameni’s measure-theoretic approach introduced \(m\)-twist and generalized twist on \(c\)-splitting sets. Under \(m\)-twist, any optimal plan \(\gamma\) admits a decomposition
\[
\gamma=\sum_{i=1}^k \alpha_i\,(\mathrm{Id}\times G_i)_\#\mu_1,
\qquad k\le m,
\]
with measurable weights \(\alpha_i\) summing to \(1\) \(\mu_1\)-a.e.; generalized twist yields an analogous countable decomposition. Thus the support is a finite or countable union of graphs over the first marginal, even when a single Monge graph does not exist [1403.3389].

The several-graph picture was further developed by Moameni in a study of plans concentrated on finitely or countably many graphs. There local differential conditions implying local \(d\)-rectifiability of the support were shown to imply a local \(1\)-twist property, and hence generalized twist under compactness. The paper also analyzed extremality and uniqueness for plans supported on several graphs, showing that multi-graph concentration is a genuine intermediate regime between Monge determinism and fully non-graphical optimal transport [1507.05923].

## 4. Special cost classes: quadratic, symmetric, and graph-structured models

The benchmark positive example is the Gangbo–Świȩch quadratic cost
\[
c(x_0,\dots,x_{m-1})=\sum_{0\le i<j\le m-1}|x_i-x_j|^2.
\]
Under finite second moments and the condition that each \(\mu_i\) vanishes on \((d-1)\)-rectifiable sets, Gangbo–Świȩch proved existence of a unique optimal plan, and that plan is supported on a graph:
\[
\theta=(T_0,\dots,T_{m-1})_\#\mu_0,\qquad T_0=I.
\]
The maps are expressed through convex potentials \(f_i\), and the problem is equivalent to a Wasserstein barycenter problem: the barycenter measure \(\nu\) minimizes \(\sum_i W_2^2(\mu_i,\nu)\), and the maps \(\nabla f_i\) are Brenier maps from \(\mu_i\) to \(\nu\) [1212.1680].

Symmetry introduces a different but equally rigid structure. In symmetric Monge–Kantorovich problems on \(\Omega^m\), with equal marginals \(\mu\) and cyclic invariance, Ghoussoub–Moameni studied costs
\[
c(x_0,\dots,x_{m-1})=\sum_{i=1}^{m-1}\langle u_i(x_0),x_i\rangle
\]
generated by vector fields. Their theorem states that the symmetric optimizer is of Monge type:
\[
\theta=(I,S,S^2,\dots,S^{m-1})_\#\mu,
\]
for a measure-preserving transformation \(S\) satisfying \(S^m=I\) a.e. The accompanying Hamiltonian representation
\[
u_i(x)=\nabla_{i+1}H(x,Sx,S^2x,\dots,S^{m-1}x)
\]
connects symmetric OT, \(m\)-cyclic monotonicity, and polar decomposition of vector fields [1212.1680].

Another structured family is given by graph-type surpluses
\[
b(x_1,\dots,x_m)=\sum_{\{i,j\}\in P} x_i\cdot x_j,
\]
where \(P\) is the edge set of a graph. Pass–Vargas used graph theory to classify when such surpluses admit Monge and unique solutions. Disconnected graphs, or graphs missing a crucial edge from the reference vertex, produce non-Monge and non-unique behavior. By contrast, complete graphs, graphs in which each vertex misses at most one edge, graphs with an inner hub, and certain tree-like gluings of clique-interaction subgraphs admit Monge solutions and uniqueness once \(\mu_1\) and selected additional marginals are absolutely continuous [2104.09488].

## 5. Failure of the Monge ansatz and minimal counterexamples

A common misconception is that benign costs, especially quadratic ones, should always favor Monge solutions. The discrete theory shows otherwise. Friesecke–Vögler constructed a transparent counterexample in the smallest symmetric finite setting: \(N=3\) marginals, \(\ell=3\) sites, and a symmetric pairwise Frenkel–Kontorova cost. The symmetric Kantorovich polytope has \(22\) extreme points, only \(7\) of which are Monge, and the unique minimizer is non-Monge. By superposition, the same mechanism yields a continuous one-dimensional example in which the Monge infimum is not attained and minimizing sequences develop microstructure [1808.04318].

For uniform discrete marginals and the quadratic multi-marginal cost
\[
c(x_1,\dots,x_N)=\sum_{i,j=1}^N |x_i-x_j|^2,
\]
a sharp counterexample was obtained in 2024. For \(N=3\), \(d=2\), \(m=3\), explicit \(3\)-empirical marginals in \(\mathbb{R}^2\) admit an optimal coupling with cost \(68.027\), while the best Monge-type coupling among the \(36\) permutation-based candidates has cost \(68.065\). Hence the Monge ansatz fails. The same paper proves that Monge always holds in three important regimes: \(N=2\) for any \(d,m\), \(d=1\) for any \(N,m\), and \(m=2\) for any \(N,d\). It follows that \((N,d,m)=(3,2,3)\) is the smallest possible triple for failure, and that the set \(E_m\) of \(m\)-empirical measures is not barycentrically convex for \(N\ge 3\), \(d\ge 2\), \(m\ge 3\) [2401.12417].

Symmetry also obstructs uniqueness. In a symmetric multi-marginal problem with equal marginals and permutation-invariant cost, if an optimizer charges a product \(S_1\times S_2\times S_3\) of pairwise disjoint sets with positive mass, then permutation symmetry generates another optimizer with the same cost. In that sense, uniqueness in symmetric multi-marginal transport occurs only under very special circumstances [1507.05923].

## 6. Applications, current extensions, and open structure

The multi-marginal Monge problem appears in several application areas. In economics, costs of matching-for-teams or hedonic pricing type can be written as
\[
c(x_1,\dots,x_m)=\inf_{z\in Z}\sum_{i=1}^m f_i(x_i,z),
\]
and under twist and non-degeneracy hypotheses this yields unique Monge solutions. In the Riemannian quadratic case \(f_i(x_i,z)=t_i d^2(x_i,z)\), the same formalism is equivalent to Wasserstein barycenters, so existence and uniqueness of the barycenter transport structure become a special case of multi-marginal Monge theory [1007.0424][1406.0026].

In density functional theory, the semi-classical limit produces multi-marginal OT with Coulomb cost. The survey literature emphasizes that the two-electron case admits a unique Monge solution under mild assumptions, while for \(m\ge 3\) the picture changes drastically: in one dimension the unique symmetric minimizer is supported on the union of graphs of powers of a map, and in higher dimensions the support can have dimension at least \(2n-2\). A different DFT-derived problem, based on an intermolecular interaction cost
\[
c(\vec x,\vec y)=\sum_{i=1}^{N_\alpha}\sum_{j=1}^{N_\beta} x_i\cdot y_j^*,
\qquad y_j^*=(-2y_j^1,y_j^2,\dots,y_j^d),
\]
admits a unique Monge solution provided \(\rho_1\) and one \(\beta\)-marginal are absolutely continuous; if all \(\beta\)-marginals are Dirac, every coupling is optimal, showing that the extra regularity assumption is necessary for uniqueness [1406.0026][2401.07880].

Recent work has also turned the Monge problem into an explicit computational target. A deep-learning approach based on Hilbert space embeddings and MMD penalties parameterizes the maps \(T_2,\dots,T_N\) by neural networks and solves a penalized multi-marginal Monge problem directly on GPUs; for the pairwise quadratic cost, the paper proves that the MMD penalties enforce the marginal constraints asymptotically [2507.09206]. A complementary dynamical formulation rewrites multi-marginal OT with convex or semi-convex translation-invariant costs as a convex optimization problem on flows of couplings over \(\Omega^k\), with a primal-dual structure and numerically computed quasi-Monge solutions [2509.22494].

Across these developments, a stable dichotomy emerges. For special costs—quadratic barycentric costs, twist-compatible one-dimensional costs, several graph-structured surpluses, and some DFT-derived anisotropic bilinear costs—optimal plans are deterministic or finitely multi-valued. Outside those classes, especially under symmetry or discreteness, optimal transport may concentrate on several graphs, on high-dimensional sets, or on nonattained limits with microstructure. The multi-marginal Monge problem is therefore less a single theorem than a hierarchy of structural regimes, organized by twist, splitting geometry, cyclical monotonicity, marginal regularity, and symmetry.

Source: https://www.emergentmind.com/topics/multi-marginal-monge-problem