---
title: Sub-Riemannian Transport
url: https://www.emergentmind.com/topics/sub-riemannian-transport
type: topic
---

# Sub-Riemannian Transport

Sub-Riemannian transport refers to the theory and analysis of optimal mass transportation where the geometry of the underlying space is sub-Riemannian rather than Riemannian. In contrast with classical (Euclidean or Riemannian) optimal transport, sub-Riemannian transport involves additional constraints on admissible curves: only horizontal trajectories (i.e., those tangent to a prescribed distribution that is bracket-generating but not full rank everywhere) are permitted. This leads to profound differences in the structure of geodesics, transport maps, regularity, and metric inequalities, and increasingly finds applications across control theory, stochastic processes, thermodynamics, data-driven modeling, and physical systems with nonholonomic constraints [2408.14707].

## 1. Sub-Riemannian Geometry and Foundations of Optimal Transport

A sub-Riemannian manifold $(M, \Delta, g)$ consists of a smooth connected manifold $M$, a smooth bracket-generating distribution $\Delta \subset TM$, and a smooth inner product $g$ defined on $\Delta$. A Lipschitz curve $\gamma: [0,1] \to M$ is horizontal iff $\dot\gamma(t) \in\Delta_{\gamma(t)}$ a.e.; its length is $l(\gamma) = \int_0^1 \|\dot{\gamma}(t)\|_g \,dt$. The sub-Riemannian (Carnot–Carathéodory) distance is
$$
d_{SR}(x,y) = \inf \{ l(\gamma) : \gamma(0)=x,\,\gamma(1)=y,\,\dot{\gamma} \in \Delta \}.
$$
The bracket-generating (Hörmander) condition ensures that $d_{SR}$ is finite and defines the topology of $M$.

Optimal mass transport (OMT) in this setting involves finding, for $\mu_0, \mu_1 \in \mathcal{P}_c(M)$, a minimizer to
$$
T: M \to M, \quad T_\#\mu_0 = \mu_1, \quad \int_M c(x,T(x))\,d\mu_0(x) \to \min,
$$
with $c(x,y) = d_{SR}(x,y)^2$. The Monge and Kantorovich formulations (relaxing to transport plans) generalize from Riemannian to sub-Riemannian geometry, but with new subtleties regarding geodesics, regularity, and the existence of abnormal minimizers [1706.07308][1705.05380][2507.20959].

## 2. Geodesics, Horizontal Distributions, and Sub-Riemannian Distance

Sub-Riemannian geodesics correspond to curves solving Hamiltonian equations governed by the sub-Riemannian Hamiltonian $H(x,p) = \frac{1}{2}g_x^*(p|_{\Delta_x},p|_{\Delta_x})$. Normal geodesics are projections of flows with $H>0$; abnormal geodesics arise as critical points of the endpoint map and may have $H=0$. The presence or absence of nontrivial abnormal generators ('ideal' structures) is decisive for the uniqueness and regularity of optimal maps: if abnormal minimizers are absent, the squared distance $d_{SR}^2$ is semiconvex away from the cut locus and standard regularity results (e.g., uniqueness of optimal transport, Jacobian estimates) hold [1705.05380].

For example, on the Heisenberg group $\mathbb{H}^n$ (as a step-2 Carnot group), the horizontal distribution is spanned by $2n$ left-invariant vector fields, and geodesics exhibit characteristic sub-Riemannian behavior: motion is 'fast' along horizontal directions and 'slow' in the directions generated via their brackets. This anisotropy is formalized in the asymptotics of the Carnot–Carathéodory distance, which locally scales linearly in horizontal and sub-linearly in “vertical” directions [2303.16052][2511.01515].

## 3. Sub-Riemannian Optimal Transport Theory: Formulations and Structure

Both static (Kantorovich) and dynamic (Benamou--Brenier) formulations admit sub-Riemannian generalizations. The Benamou--Brenier framework seeks curves $(\mu_t,v_t)$ subject to a sub-Riemannian continuity equation
$$
\partial_t \mu_t + \operatorname{div}(\mu_t v_t) = 0,
$$
with $v_t(x) \in \Delta_x$ and cost $\int_0^1\int_M \|v_t(x)\|^2_g\,d\mu_t(x)\,dt$. The equivalence between dynamic and static formulations holds under completeness and absence of abnormal minimizers: every optimal plan arises from displacement interpolation along normal geodesics [2507.20959]. The geodesic equation in probability space inherits sub-Riemannian structure, and the tangent space at a density $\rho$ is characterized by the sub-elliptic continuity equation $\dot{\rho} = -\nabla\cdot(\rho a a^\top \nabla\phi)$, with metric
$$
g_{sr}(\dot{\rho}_1, \dot{\rho}_2) = \int_M (a^\top\nabla\phi_1)\cdot(a^\top\nabla\phi_2) \rho\,dx,
$$
mirroring the Otto formalism but using the horizontal Laplacian [1910.07480].

Uniqueness and regularity of Monge minimizers can be proven in settings where the singular set (points joined via singular minimizers) is negligible, leveraging the C$^2$-semiconvexity of Kantorovich potentials and geometric growth properties of the distribution [1706.07308].

## 4. Holonomy, Mixing, and Labeled-Particle Sub-Riemannian OMT

A major development is the analysis of 'labeled-particle' optimal transport, where the sub-Riemannian structure arises in the space of covariance factors in the finite-dimensional Gaussian case. The configuration manifold is $GL^+(n)$, with the projection $\pi: GL^+(n) \to \operatorname{Sym}^+(n)$ a principal $SO(n,\Sigma_{\text{ref}})$-bundle. The horizontal distribution is determined by the Ehresmann connection:
\[
\operatorname{Hor}_\Phi = \{ \dot{\Phi}: \dot{\Phi}\Phi^{-1} \in \operatorname{Sym}(n)\}.
\]
The corresponding sub-Riemannian metric is $\operatorname{Tr}(A\Sigma_{\text{ref}}A)$ with $\dot{\Phi}=A\Phi$, $A=A^\top$. Normal geodesics satisfy a Hamiltonian system with no abnormal extremals.

A novel phenomenon emerges: along closed loops in the base (covariances), horizontal lifts yield a holonomy $\Theta\in SO(n,\Sigma_{\text{ref}})$ which measures the 'mixing' or permutation of labels. In classical (Riemannian) OMT, such mixing is generic and unconstrained, but in the sub-Riemannian model, the holonomy group is generated by non-integrability, and control protocols (e.g., isoparallel-mass-transport, IMT) can be designed to guarantee $\Theta=I$ ('no mixing'): tracer particles return to their original positions after traversing closed curves. This structure is fully developed in [2408.14707].

## 5. Interpolation Inequalities, Distortion Coefficients, and Metric Geometry

On ideal sub-Riemannian manifolds, distortion coefficients $\beta_t(x,y)$, defined via the Jacobian of the sub-Riemannian exponential map (using Jacobi fields), underpin sharp Brunn–Minkowski, Borell–Brascamp–Lieb, and interpolation inequalities. For $t \in [0,1]$,
\[
\rho_t(\gamma_{x,y}(t))^{-1/n} \geq \tau_{1-t}(x,y)\rho_0(x)^{-1/n} + \tau_t(x,y)\rho_1(y)^{-1/n},
\]
with $\tau_t(x,y) = \beta_t(x,y)^{1/n}$, and $\gamma_{x,y}(t)$ the unique normal geodesic from $x$ to $y$ at time $t$. The geodesic ‘dimension’ $N$ (generally $N\geq n$) quantifies the mass-spreading rate, deviating from the Riemannian case due to horizontal-vertical splitting and non-integrability [1705.05380]. Explicit exponent calculations are available for the Heisenberg group, generalized H-type Carnot groups, and the Grushin plane; for instance, the Brunn–Minkowski formula on $\mathbb{H}^3$ involves exponent $5/3$ rather than $1$.

## 6. Sub-Riemannian Transport in Applied and Stochastic Settings

Sub-Riemannian transport structures naturally occur in complex systems with nonholonomic constraints, such as:

- **Ocean neutral transport**: Here, water parcels move along (locally defined) neutral planes—contact distributions $D_x = \ker \eta(x)$, with $\eta$ built from climatological gradients. The horizontal (neutral) constraint and non-integrability (helicity) yield global accessibility (by the Chow–Rashevskii theorem) and highly anisotropic geodesics. Stochastic analogues involve hypoelliptic Brownian diffusions tangent to $D$; transition densities reflect the sub-Riemannian distance, and mixing times for dianeutral dispersion are estimated at $\sim$centuries [2511.01515].

- **Congested transport**: For the Heisenberg group, the sub-Riemannian constraint restricts admissible paths to horizontal curves in the congested optimal transport problem, leading to modified Eulerian and Lagrangian formulations of traffic intensity and equilibrium (Wardrop) states [2303.16052].

- **Transport on path spaces and inequalities**: Talagrand-type transport inequalities can be generalized to horizontal Brownian motion on step-2 Carnot groups, but projection arguments used in the Euclidean setting for path-space measures fail due to the non-commutativity and vertical blows-up. Riemannian approximations recover some inequalities for $p<2$, but not at $p=2$ [2602.06646].

## 7. Regularity, Cut Locus, and Advanced Metric Properties

The regularity theory for sub-Riemannian transport is complex: the squared distance function $d_{SR}^2/2$ is semiconvex away from the cut locus, which is larger and subtler than in the Riemannian case. The cut locus is characterized as the set where semiconvexity fails. The analysis of Jacobi fields and sub-Riemannian index forms provides comparison theorems (e.g., for Sasakian manifolds) and controls on the Hessian of the distance. In applied settings, these results enable comparison principles and coupling estimates for hypoelliptic diffusions and underpin entropy dissipation and $\Gamma$-calculus for kinetic Fokker-Planck equations on sub-Riemannian density manifolds [2212.07715][1910.07480].

---

**References:**  
- [2408.14707] Sub-Riemannian Geometry, Mixing, and the Holonomy of Optimal Mass Transport  
- [1706.07308] Mass Transportation on sub-Riemannian structures of rank two in dimension four  
- [2303.16052] Transport densities and congested optimal transport problem in the Heisenberg group  
- [2602.06646] Talagrand-type transport inequalities for path spaces over Carnot groups  
- [2507.20959] Benamou-Brenier and Kantorovich are equivalent on sub-Riemannian manifolds with no abnormal geodesics  
- [1910.07480] Entropy dissipation for degenerate stochastic differential equations via sub-Riemannian density manifold  
- [2212.07715] Variations of the sub-Riemannian distance on Sasakian manifolds with applications to coupling  
- [1705.05380] Sub-Riemannian interpolation inequalities  
- [2511.01515] Ocean neutral transport: sub-Riemannian geometry and hypoelliptic diffusion

Source: https://www.emergentmind.com/topics/sub-riemannian-transport