Home Mathematics A continuous model of transportation in the Heisenberg group
Article Open Access

A continuous model of transportation in the Heisenberg group

  • Michele Circelli ORCID logo EMAIL logo and Albert Clop ORCID logo
Published/Copyright: August 29, 2025

Abstract

We present a minimization problem with a horizontal divergence-type constraint in the Heisenberg group. Our study explores its dual formulation and examines its relationship with the congested optimal transport problem, for 1 < p < + , as well as the Monge–Kantorovich problem, in the limit case p = 1 .

MSC 2020: 49Q22

1 Introduction

The optimal transport problem, originally posed by G. Monge in [37] in the late 18th century, has become a highly active research area due to its numerous applications, as highlighted in [29] and [43]. In the Monge–Kantorovich formulation, which dates back to 1942, we are given two probability measures μ and ν on a metric space ( X , d ) and aim to solve the following minimization problem:

(1.1) inf γ Π ( μ , ν ) X × X c ( x , y ) 𝑑 γ ( x , y ) ,

where the set Π ( μ , ν ) consists of transport plans, i.e. probability measures γ on the product space X × X , with prescribed marginals μ and ν. The cost function c 0 is typically a power of the distance d. Solutions exist under quite general assumptions, such as when X is a Polish space and c is lower semicontinuous and bounded from below.

Initially, the existence of solutions to (1.1) induced by measurable maps was extensively studied within Euclidean and Riemannian settings by many researchers (see, for instance, [3, 45, 30] and references therein for comprehensive overviews). Later, the research focus shifted towards sub-Riemannian frameworks, especially Carnot groups, which arise as limits of Riemannian manifolds in a very precise sense (see [36]).

The Heisenberg group n is the simplest model of non-commutative Carnot group: it is characterized by the selection of 2 n vector fields, X 1 , , X 2 n , in 2 n + 1 , known as horizontal vector fields, and a metric on the sub-bundle H n they generate on the tangent bundle T n . Since the number of vector fields is less than the dimension of the space, the metric is non-Riemannian at every point. Moreover, displacement occurs only along integral curves of horizontal vector fields, also called horizontal curves. The presence of this metric naturally leads to the definition of an intrinsic sub-Riemannian distance d SR .

The first significant work on the Monge–Kantorovich problem in n

(1.2) inf γ Π ( μ , ν ) n × n d SR ( x , y ) α 𝑑 γ ( x , y )

is [5], where L. Ambrosio and S. Rigot proved the existence and the uniqueness of the solution to (1.2) for α = 2 , and that this solution is induced by a map; subsequently, this result was extended to general sub-Riemannian manifolds in [2, 31], and to essentially non-branching metric measure spaces in [20]. The first work on the problem (1.2) for α = 1 is [27], where the authors proved the existence of solutions induced by maps. Subsequently, the analogous result was proven in more general measure metric spaces, see for instance [9] and [19]. One of the main obstacles in transitioning from the Euclidean to the sub-Riemannian setting is the loss of uniqueness of geodesics, even locally, and hence, addressing the problem in this new framework required new and very different approaches due to the unique geometry of the space.

In this direction we continue, in this paper, the study initiated in [23] on the congested optimal transport problem in n : this study is motivated by the proposal of Petitot and Tondut in [42] to model the functional geometry of the visual cortex (see also [25], where this model was expressed using sub-Riemannian geometry). The visual signal propagates along the connectivity, leading to congestion phenomena that can be modeled with the congested optimal transport problem in n (see [23, Introduction] and references therein).

The congested optimal transport problem was first introduced in the Euclidean setting by Carlier, Jimenez and Santambrogio in [18]. In the Heisenberg group it reads as follows: given two probability measures μ and ν over the closure of a bounded domain Ω, and a cost function G : Ω × + + , which is convex and exhibits p-growth in the second variable, it consists of solving

(1.3) min Q 𝒬 H p ( μ , ν ) Ω G ( x , i Q ( x ) ) 𝑑 x .

Here, the set 𝒬 H p ( μ , ν ) consists of probability measures Q over the space H of horizontal curves, with prescribed projections μ and ν and for which the associated horizontal traffic intensity i Q

Ω ¯ φ ( x ) 𝑑 i Q ( x ) := H ( 0 1 φ ( σ ( t ) ) | σ ˙ ( t ) | H 𝑑 t ) 𝑑 Q ( σ ) for all  φ C ( Ω ¯ )

is a p-summable function; | | H denotes a norm on the sub-bundle H n . In [23] the authors proved that solutions to (1.3) are equilibrium configurations of Wardrop type. Although (1.3) is a convex minimization problem, and hence the existence of solutions follows from the direct method in the calculus of variations, it is also a minimization problem over a set of measures over curves, and hence it has two layers of infinite dimensions. Therefore, in the spirit of [12] and [44], we investigate the relation between (1.3) and the following more tractable problem:

(1.4) min 𝒘 Ω 𝒢 ( x , 𝒘 ( x ) ) 𝑑 x ,

where

𝒢 ( x , 𝒘 ) = G ( x , | 𝒘 | H ) ,

and the minimum is taken among all p-summable horizontal vector fields 𝒘 : Ω H n , with prescribed horizontal divergence

(1.5) div H 𝒘 = μ - ν .

This problem is the Heisenberg version of the well-known continuous model of transportation introduced in the Euclidean setting by Beckmann in [7]. Following the approach [12] and [44], which relies on the Dacorogna–Moser flow argument, together with a standard regularization procedure, we prove that the two problems (1.4) and (1.3) are actually equivalent. Moreover, in some sense, minimizers to (1.3) and (1.4) correspond to each other, see Proposition 5.13 and Theorem 5.14. The main difference here with the Euclidean case is that, in order to apply the Dacorogna–Moser flow argument, we must pass through a Riemannian approximation of the Heisenberg group.

Using standard duality arguments, (1.4) can be formulated as a classical problem in the calculus of variations, and hence it is strongly related to the regularity theory for PDEs in the Heisenberg group. More precisely, in the spirit of [46], we show in Theorem 5.10 that (1.4) is equivalent to the following problem:

(1.6) max φ H W 1 , q ( Ω ) { - Ω φ d ( μ - ν ) - Ω 𝒢 * ( x , H φ ( x ) ) 𝑑 x } ,

where 𝒢 * denotes the Legendre transform of 𝒢 in the second variable.

Here, H φ denotes the vector of derivatives of φ in horizontal directions, and H W 1 , q ( Ω ) is its associated horizontal Sobolev space, where the exponent q is the dual of p (see Section 2 for its definition). Moreover, if ( 𝒘 0 , φ 0 ) is a pair of minimizer for (1.4) and maximizer for (1.6), then, if 𝒢 * is differentiable, we have

𝒘 0 ( x ) = D 𝒢 * ( x , H φ 0 ( x ) ) a.e. in  Ω ,

where D 𝒢 * is the Euclidean gradient of 𝒢 * in 2 n + 1 , and φ 0 solves the associated Euler–Lagrange equation of horizontal divergence-type:

(1.7) div H ( D 𝒢 * ( x , H φ ( x ) ) ) = μ - ν .

In the Euclidean setting, the equivalence between (1.3), (1.4) and (1.6) was first investigated in [13] for the autonomous case. More recently, this equivalence was extended to the non-autonomous case in [14], where the authors also assume that the right-hand side in (1.5) belongs to the dual of some Sobolev space, whose elements are not measures in general.

In the limiting case p = 1 , the simplest version of the problem (1.4) reads as

(1.8) inf 𝒘 L 1 { Ω | 𝒘 ( x ) | H 𝑑 x : div H 𝒘 = μ - ν } .

As observed in [45, Subsection 4.2.2], this problem is not well-posed a priori: due to the non-reflexivity of L 1 , there may not exist any L 1 horizontal vector field minimizing the L 1 norm under a horizontal divergence constraint. This is why one states the problem (1.8) in the larger space of compactly supported horizontal vector measures, whose horizontal divergence is the signed Radon measure μ - ν :

(1.9) min { 𝒘 : div H 𝒘 = μ - ν } ,

where 𝒘 is the total variation of the vector measure 𝒘 . In the spirit of [45, Section 4], we show in Theorem 4.4 that (1.9) equals the value of the Monge–Kantorovich problem (1.2), with α = 1 . A key role in the proof of this equivalence is played by the transport density, which measures how much transport is taking place in any region of n . We delve deeper into the study of the properties of the transport density in n , initiated in [23], explicitly describing its relation with Kantorovich potentials. To do so, we adapt some results on the differentiability of the Kantorovich potentials present in [30], which were specific to the Riemannian case. While adapting these results, we strongly rely on results from [5] and [21], which were based on the explicit computation of geodesics, the smoothness of the pointed sub-Riemannian distance, and the measure contraction property.

As byproducts of the previous equivalence, we obtain two important facts: first, if either μ or ν is absolutely continuous with respect to the Lebesgue measure (which is a Haar measure for the group structure), then (1.8) admits a solution. Second, one can conclude that the problem (1.9) also admits a dual formulation (as it does for its congested version (1.4)), which is given by the well-known Kantorovich duality theorem:

max u H W 1 , { u d ( μ - ν ) : H u 1 } .

In addition, in Theorem 4.7 we prove that the Lagrangian reformulation of (1.9) is

min Q H SR ( σ ) 𝑑 Q ( σ ) = min Q 𝑑 i Q ( x ) ,

where SR ( σ ) denotes the sub-Riemannian length (see Section 2 for its definition). We emphasize that in the proof, the geometry of the space plays no role, so the result holds in the more general setting of Polish metric spaces.

The paper is organized as follows: In Section 2, we introduce the Heisenberg group and recall some well-known results on optimal transport theory in this setting. In Section 3, we prove some results about the transport set and the Kantorovich potential, which will be useful later. Section 4 deals with the aforementioned three problems in the limit case p = 1 , while the last section addresses the equivalence of the three problems for 1 < p < + . First, we introduce both (1.4) and (1.6) in their natural functional analytic setting; second, we prove the equivalence between (1.4) and the congested optimal transport problem (1.3).

2 Preliminaries

2.1 The Heisenberg group n

Let n 1 . The n-th Heisenberg group n is the connected and simply connected Lie group, whose Lie algebra 𝔥 n is stratified of step 2: i.e. it is the direct sum of two linear subspaces

𝔥 n = 𝔥 1 n 𝔥 2 n ,

where 𝔥 1 = span { X 1 , , X n , X n + 1 , , X 2 n } is the horizontal layer, 𝔥 2 = span { X 2 n + 1 } and the only non-trivial bracket-relation between the vector fields X 1 , , X n , X n + 1 , , X 2 n , X 2 n + 1 is

[ X j , X n + j ] = X 2 n + 1 for all  j = 1 , , n .

The homogeneous dimension of n is

N := i = 1 2 i dim ( 𝔥 i n ) = 2 n + 2 .

The exponential map exp : 𝔥 n n is a diffeomorphism, hence it induces a global diffeomorphism between n and 2 n + 1 :

n exp ( x 1 X 1 + x 2 X 2 + + x 2 n X 2 n + x 2 n + 1 X 2 n + 1 ) ( x 1 , , x 2 n + 1 ) 2 n + 1 .

In this system of coordinates the group law will be given, through the Baker–Campbell–Hausdorff formula which defines a group structure on 𝔥 n , by

x y := ( x 1 + y 1 , , x 2 n + y 2 n , x 2 n + 1 + y 2 n + 1 + 1 2 j = 1 2 n ( x j y n + j - x n + j y j ) ) ,

the unit element e n is 0 2 n + 1 , the center of the group is

L := { ( 0 , , 0 , x 2 n + 1 ) n : x 2 n + 1 }

and the vector fields X 1 , , X 2 n + 1 , left invariant with respect to (2.1), in coordinates read as

{ X j := x j - x n + j 2 x 2 n + 1 , j = 1 , , n , X n + j := x n + j + x j 2 x 2 n + 1 , j = 1 , , n , X 2 n + 1 := x 2 n + 1 .

According to the two step-stratification of the algebra, n is naturally endowed with a family of intrinsic non-isotropic dilations that reads as

δ τ ( ( x 1 , , x 2 n + 1 ) ) := ( τ x 1 , , τ x 2 n , τ 2 x 2 n + 1 ) for all  x n  an all  τ > 0 .

See [11, Section 2] for more details about this topic.

The horizontal layer 𝔥 1 n of the algebra defines a sub-bundle H n of the tangent bundle T n , whose fibre at a point x n is

H x n = span { X 1 ( x ) , , X n ( x ) , X n + 1 ( x ) , , X 2 n ( x ) } .

One can equip n with a left-invariant sub-Riemannian metric as follows: we fix an inner product , H on 𝔥 1 n that makes { X 1 , , X n , X n + 1 , , X 2 n } an orthonormal basis and we denote by | | H the norm associated with such a scalar product. We keep the same notation both for the corresponding scalar product and norm on each fibre H x n , x n .

An absolutely continuous curve σ AC ( [ 0 , 1 ] , 2 n + 1 ) is said to be horizontal if its velocity vector σ ˙ ( t ) belongs to H σ ( t ) n at almost every t [ 0 , 1 ] . Its sub-Riemannian length is

(2.1) SR ( σ ) := 0 1 | σ ˙ ( t ) | H 𝑑 t .

Let us just remark that SR is invariant under reparametrizations of σ.

The Rashevsky–Chow’s theorem guarantees that the sub-Riemannian distance

d SR ( x , y ) := inf { SR ( σ ) : σ  is horizontal , σ ( 0 ) = x , σ ( 1 ) = y } for all  x , y n

is well-defined and it induces the Euclidean topology, see [1, Theorem 3.31]. In particular, ( n , d SR ) is a Polish space.

If x n and r > 0 , we denote by

B ( x , r ) := { y n : d SR ( x , y ) < r }

the sub-Riemannian ball centered in x with radius r. The Lebesgue measure 2 n + 1 is a Haar measure on n 2 n + 1 ; we will often denote it by dx.

For any open set Ω n and any measurable function φ : Ω , we denote by

H φ = i = 1 2 n X i φ X i

its horizontal gradient, in the sense of distributions.

We say that a measurable map φ : Ω is Pansu differentiable at x n if there exists an homogeneous group homomorphism L : n such that

lim y x φ ( y ) - φ ( x ) - L ( x - 1 y ) d SR ( x , y ) = 0 ,

where a group homomorphism L : n is homogeneous if L ( δ τ ( x ) ) = τ L ( x ) , for any x n and any τ > 0 .

If the map L exists, it is unique and will be denoted by D H φ ( x ) .

If φ : Ω is Pansu differentiable at x Ω then φ has directional derivatives at x along all the directions X j for all j = 1 , , 2 n . Moreover, if y = ( y 1 , , y 2 n + 1 ) n , then

(2.2) D H φ ( x ) ( y ) = j = 1 2 n y j X j φ ( x ) = H φ ( x ) , π x ( y ) H ,

where the map y π x ( y ) is the smooth section of H n defined as

(2.3) π x ( y ) := i = 1 2 n y i X i ( y ) ,

see [33, Proposition 5.6]. If φ C and σ is an horizontal curve, a simple computation based on the classical chain rule implies that

(2.4) d d t φ ( σ ( t ) ) | t = s = H φ ( σ ( s ) ) , σ ˙ ( s ) H ,

for a.e. s [ 0 , 1 ] .

For every 1 q , the space

(2.5) H W 1 , q ( Ω ) := { φ : Ω  measurable : φ L q ( Ω ) , H φ L q ( Ω , H Ω ) }

is a Banach space equipped with the norm

φ H W 1 , q ( Ω ) := φ L q ( Ω ) + H φ L q ( Ω , H Ω ) .

2.1.1 Geodesics in n

We call minimizing horizontal curve any curve σ AC ( [ 0 , 1 ] , 2 n + 1 ) , which is horizontal and such that

SR ( σ ) = d SR ( σ ( 0 ) , σ ( 1 ) ) ,

and geodesic any minimizing horizontal curve parametrized proportionally to the arc-length.

The space ( n , d SR ) is a complete and locally compact metric length space, therefore it is geodesic [15, Theorem 2.5.28].

Usually, in this kind of spaces, geodesics between two points may not be unique: in the case of the Heisenberg group geodesics can be computed explicitly and it is possible to detect a set in which these are unique.

We denote by

(2.6) E := { ( x , y ) n × n : x - 1 y L } ,

then it holds the following characterization for geodesics parametrized on [ 0 , 1 ] .

Theorem 2.1.

A non-trivial geodesics parametrized on [ 0 , 1 ] and starting from 0 is the restriction to [ 0 , 1 ] of the curve σ χ , θ ( t ) = ( x 1 ( t ) , , x 2 n + 1 ( t ) ) either of the form

(2.7) x j ( t ) = χ j sin ( θ s ) - χ n + j ( 1 - cos ( θ s ) ) θ , j = 1 , , n ,
(2.8) x n + j ( t ) = χ n + j sin ( θ s ) + χ j ( 1 - cos ( θ s ) ) θ , j = 1 , , n ,
(2.9) x 2 n + 1 ( t ) = | χ | 2 2 θ 2 ( θ s - sin ( θ s ) )

for some χ R 2 n { 0 } and θ [ - 2 π , 2 π ] { 0 } , or of the form

( x 1 ( t ) , , x 2 n + 1 ( t ) ) = ( χ 1 t , , χ 2 n t , 0 )

for some χ R 2 n { 0 } and θ = 0 . In particular, it holds

| χ | 2 n = | σ ˙ χ , θ | H = d SR ( x , σ χ , θ ( 1 ) ) .

Moreover, it holds:

  1. For all ( x , y ) E , there is a unique geodesic x σ χ , θ parametrized on [ 0 , 1 ] between x and y , for some χ 2 n { 0 } and some θ ( - 2 π , 2 π ) .

  2. If ( x , y ) E , then x - 1 y = ( 0 , , 0 , z 2 n + 1 ) for some z 2 n + 1 { 0 } , there are infinitely many geodesics parametrized on [ 0 , 1 ] between x and y . These curves are all curves of the form x σ χ , 2 π , if z 2 n + 1 > 0 , x σ χ , - 2 π , if z 2 n + 1 < 0 , for all χ 2 n such that | χ | = 4 π | z 2 n + 1 | .

From the description of geodesics given in Theorem 2.1 it follows that ( n , d SR ) is a non-branching metric space: any two geodesics which coincide on a non-trivial interval coincide on the whole intersection of their intervals of definition.

We denote by Geo ( n ) the space of geodesics parametrized on [ 0 , 1 ] :

Geo ( n ) := { σ AC ( [ 0 , 1 ] , 2 n + 1 ) : σ  is a geodesic } .

Let us just remark that there exists a map

(2.10) S : n × n Geo ( n ) ,

which is γ-measurable for any positive Borel measure γ on n × n , and that associates with any pair of points ( x , y ) n × n a geodesic S ( x , y ) := σ x , y Geo ( n ) , between x and y. The existence of such a map follows from the theory of Souslin sets and general theorems about measurable selections, see for instance [10, Theorem 6.9.2 and Theorem 7.4.1]. Moreover, S is jointly continuous on E. If e t is the evaluation map at time t [ 0 , 1 ] , then the map

S t := e t S : n × n n

associates to any two points x , y n the point S t ( x , y ) := σ x , y ( t ) of n at distance t d SR ( x , y ) from x on the selected geodesic σ x , y between x and y.

2.1.2 H Ω -valued finite Radon measures

Let Ω n be a set, and denote by H Ω the restriction of the horizontal bundle H n to the set Ω. We denote by C c ( Ω , H Ω ) the class of continuous horizontal vector fields with compact support in Ω, and by C 0 ( Ω , H Ω ) its completion with respect to the uniform norm

ϕ = sup { | ϕ ( x ) | H : x Ω } , ϕ C ( Ω , H Ω ) .

The normed space ( C 0 ( Ω , H Ω ) , ) is a Banach space.

An H Ω -valued finite Radon measure 𝒘 on Ω is a continuous linear functional on C 0 ( Ω , H Ω )

𝒘 C 0 ( Ω , H Ω ) ,

and its action against a horizontal vector field ϕ C 0 ( Ω , H Ω ) is denoted by

Ω ϕ 𝑑 𝒘 .

We denote by

( Ω , H Ω ) := C 0 ( Ω , H Ω )

the space of all H Ω -valued finite Radon measures on Ω and we equip it with the dual norm

𝒘 ( Ω , H Ω ) := 𝒘 ( C 0 ( Ω , H Ω ) ) = sup { Ω ϕ 𝑑 𝒘 : ϕ C 0 ( Ω , H Ω ) , ϕ 1 } .

Since the horizontal bundle H Ω has a global trivialization, it follows from the Riesz theorem (see for instance [4]) that ( Ω , H Ω ) can be identified with the space ( Ω ) 2 n of 2 n -valued finite Radon measure through the identification

Ω ϕ 𝑑 𝒘 = i = 1 2 n Ω ϕ i 𝑑 𝒘 i for all  ϕ = i = 1 2 n ϕ i X i C 0 ( Ω , H Ω ) ,

for some unique ( 𝒘 1 , , 𝒘 2 n ) ( Ω ) 2 n . Moreover, it holds

𝒘 ( Ω , H Ω ) = | 𝒘 | ( Ω ) ,

where | 𝒘 | + ( Ω ) denotes the total variation of the vector measure ( 𝒘 1 , , 𝒘 2 n ) .

Finally, if Ω is compact, then C c ( Ω , H Ω ) = C 0 ( Ω , H Ω ) = C ( Ω , H Ω ) and

( Ω , H Ω ) = ( C ( Ω , H Ω ) ) .

2.1.3 Mollification in n

Let us consider a mollifier for the group structure, i.e. a function ρ C c ( n ) such that ρ 0 , supp ( ρ ) B ( 0 , 1 ) and n ρ ( x ) 𝑑 x = 1 ; for any ε > 0 we denote by

ρ ε ( x ) := ε - N ρ ( δ 1 ε ( x ) ) .

If φ L loc 1 , then

ρ ε φ ( x ) := n ρ ε ( x y - 1 ) φ ( y ) 𝑑 y

is smooth and it enjoys many standard properties, see for instance [32].

More generally, one may also mollify Radon measures: given a finite Radon measure λ ( n ) , resp. a finite vector Radon measure 𝒘 = ( 𝒘 1 , , 𝒘 2 n ) ( n , H n ) , then the mollified functions λ ε C ( n ) , resp. 𝒘 ε C ( n , H n ) , are defined as

λ ε ( x ) := ρ ε λ ( x ) = n ρ ε ( x y - 1 ) 𝑑 λ ( y )

and

𝒘 ε := j = 1 2 n 𝒘 j ε X j ,

where 𝒘 j ε ( x ) := ρ ε 𝒘 j ( x ) , for any j = 1 , , 2 n .

2.1.4 Riemannian approximation of n

The Heisenberg group ( n , d SR ) arises as the pointed Hausdorff–Gromov limit of Riemannian manifolds, defined as follows, in which the non-horizontal direction is increasingly penalized. One can equip n 2 n + 1 with a left-invariant Riemannian metric g ϵ , with ϵ > 0 , that has

X 1 , , X 2 n , ϵ X 2 n + 1

as orthonormal frame for T n . We denote by | | ϵ := g ϵ ( , ) and by d ϵ the associated Riemannian distance: it turns out to be left-invariant with respect to (2.1), since X 1 , , X 2 n , ϵ X 2 n + 1 are themselves left invariant. It follows from the definition that, for fixed x , y n , the function d ϵ ( x , y ) is decreasing in ϵ and for any ϵ ( 0 , ϵ ¯ ) ,

d ϵ ( x , y ) d SR ( x , y ) .

Moreover, for any x , y n it holds that

d SR ( x , y ) = lim ϵ 0 d ϵ ( x , y ) ,

see [36], and the convergence is uniform on compact subsets, see [17]. In [35, Theorem 1.1] the author proved that ( n , d ϵ ) converges to ( n , d SR ) in the pointed Gromov–Hausdorff sense, as ϵ goes to 0.

In the end, if φ C ( n ) , we denote by

ϵ φ := i = 1 2 n X i φ X i + ϵ 2 X 2 n + 1 φ X 2 n + 1 .

2.2 Optimal transport theory in n

This subsection contains some well-known results about the Monge–Kantorovich problem in n with cost function equal to the Carnot–Carathéodory distance.

Let ( M 1 , d 1 ) and ( M 2 , d 2 ) be two Polish spaces and let F : M 1 M 2 be a Borel map. Given a finite Radon measure λ ( M 1 ) , we denote by F # λ ( M 2 ) the push-forward measure, that is the finite Radon measure defined as F # λ ( A ) := λ ( F - 1 ( A ) ) for all A Borel set in M 2 . By 𝒫 ( M ) we denote the subset of ( M ) consisting of probability measures on the metric space M, and by 𝒫 c ( M ) the subset of compactly supported ones.

Now, given μ , ν 𝒫 c ( n ) , we denote by

Π ( μ , ν ) = { γ 𝒫 ( n × n ) : ( π 1 ) # γ = μ , ( π 2 ) # γ = ν }

the set of transport plans between μ and ν, where π 1 and π 2 are the projection on the first and second factor, respectively; this set is compact with respect to the weak convergence of measures.

It is well known that the Monge–Kantorovich problem between μ and ν, associated with the Carnot–Carathéodory distance,

(2.11) inf γ Π ( μ , ν ) n × n d SR ( x , y ) 𝑑 γ ( x , y ) ,

admits at least one solution. See for instance [45, Theorem 1.7]. We denote by

Π 1 ( μ , ν ) := { γ Π ( μ , ν ) : γ  solves (2.11) }

the set of optimal transport plans. It is a closed subset of the compact set Π ( μ , ν ) , with respect to the weak convergence of measures.

Let us denote by

(2.12) Lip 1 ( n , d SR ) := { u : n : | u ( x ) - u ( y ) | d SR ( x , y )  for all  x , y n } .

It holds the following important theorem, see for instance [45, Proposition 3.1].

Theorem 2.2 (Kantorovich duality theorem).

There exists u ¯ Lip 1 ( H n , d SR ) that solves

(2.13) sup { n u d ( μ - ν ) : u Lip 1 ( n , d SR ) } .

Moreover, it holds that

min γ Π ( μ , ν ) n × n d SR ( x , y ) 𝑑 γ ( x , y ) = n u ¯ d ( μ - ν ) ,

and γ Π ( μ , ν ) is optimal if, and only if,

u ¯ ( x ) - u ¯ ( y ) = d SR ( x , y ) γ -a.e. in  n × n .

Any u Lip 1 ( n , d SR ) that solves (2.13) is a Kantorovich potential.

If at least one between μ and ν is absolutely continuous with respect to the Haar measure, then the explicit representation of geodesics implies that any optimal transport plan is concentrated on the set E n × n of pairs of points connected by a unique geodesic, see (2.6) for its definition.

Lemma 2.3.

Let us suppose that either μ L 2 n + 1 , or ν L 2 n + 1 , and γ Π 1 ( μ , ν ) . Then for γ-a.e. ( x , y ) in H n × H n there exists a unique geodesic between x and y. That is,

γ ( { ( x , y ) n × n : x - 1 y L } ) = 0 .

For the proof of this result we refer to [27, Lemma 4.1].

3 Properties of transport rays and Kantorovich potentials in n

In the spirit of [45, Chapter 3] and [30, Section 2] we extend to the Heisenberg setting some properties of transport rays and of transport set.

3.1 Transport rays and transport set

Let us consider μ , ν 𝒫 c ( n ) and a Kantorovich potential u Lip 1 ( n , d SR ) .

Let us just remark that if x supp ( μ ) and y supp ( ν ) are such that

(3.1) u ( x ) - u ( y ) = d SR ( x , y ) ,

then

u ( σ ( t ) ) = u ( x ) - d SR ( x , σ ( t ) ) for all  t [ 0 , 1 ] ,

where σ is any geodesic between x and y.

Definition 3.1.

A transport ray is a non-trivial geodesic σ : [ 0 , 1 ] n such that

  1. σ ( 0 ) supp ( μ ) and σ ( 1 ) supp ( ν ) ,

  2. u ( σ ( 0 ) ) - u ( σ ( 1 ) ) = d SR ( σ ( 0 ) , σ ( 1 ) ) .

We denote by 𝒯 1 the set of all points which lie on transport rays

𝒯 1 := { x σ ( [ 0 , 1 ] ) : σ  is transport ray } ,

and by 𝒯 0 the complementary set of rays of length zero

𝒯 0 := { x supp ( μ ) supp ( ν ) : | u ( x ) - u ( z ) | < d SR ( x , z )  for all  z supp ( μ ) supp ( ν ) , x z } .

We will call transport set the set

(3.2) 𝒯 := 𝒯 1 𝒯 0 .

We observe that

(3.3) supp ( μ ) supp ( ν ) 𝒯 .

Moreover, as in [30, Lemma 8] the following result holds.

Theorem 3.2.

The transport set T is compact.

Proof.

Thanks to Hopf–Rinow theorem it is enough to prove that 𝒯 is a closed and bounded set. Let us consider the function v : n × n , v ( x , y ) = u ( x ) - u ( y ) . Since v is continuous, it attains a maximum L < on supp ( μ ) × supp ( ν ) , which is a compact set. Let us prove that L 0 . If supp ( μ ) supp ( ν ) , then for all x supp ( μ ) supp ( ν ) we have that ( x , x ) supp ( μ ) × supp ( ν ) and v ( x , x ) = 0 . Otherwise, from (3.3) it follows that 𝒯 1 , hence there exists at least one transport ray σ. If x = σ ( 0 ) and y = σ ( 1 ) , then v ( x , y ) = d SR ( x , y ) > 0 . Hence L 0 .

We can suppose that A := 𝒯 supp ( μ ) supp ( ν ) ; otherwise, from the previous theorem it follows that 𝒯 = supp ( μ ) supp ( ν ) , which is compact. Hence, any z A lies on a transport ray σ z . Let us denote by a = σ z ( 0 ) and b = σ z ( 1 ) . Then

d SR ( a , z ) + d SR ( b , z ) = d SR ( a , b ) = v ( a , b ) L .

Hence, A lies in the union of the L-neighborhoods of the compact sets supp ( μ ) and supp ( ν ) , thus 𝒯 is bounded.

Let us prove that 𝒯 is closed. Let us consider ( z n ) n 𝒯 , converging to some z, we prove that z 𝒯 . If there exists a subsequence ( z n k ) k supp ( μ ) supp ( ν ) , then z supp ( μ ) supp ( ν ) by the compactness of supp ( μ ) and supp ( ν ) . Let us suppose that z n A for all n : there exists a transport ray σ n , whose endpoints we denote by a n = σ n ( 0 ) and b n := σ n ( 1 ) . There exist two subsequences

(3.4) a n j a supp ( μ ) and b n j b supp ( ν ) when  j .

Since

d SR ( z n j , a n j ) + d SR ( z n j , b n j ) = d SR ( a n j , b n j ) = u ( a n j ) - u ( b n j ) for all  n j ,

we have that

(3.5) d SR ( z , a ) + d SR ( z , b ) = d SR ( a , b ) = u ( a ) - u ( b ) .

Hence there are two possibilities

  1. a = b z = a = b 𝒯 by (3.5),

  2. a b z 𝒯 1 by (3.4) and (3.5). ∎

3.2 Differentiability of Kantorovich potential

Following [30, Lemma 10], one can prove that any Kantorovich potential u is Pansu differentiable in the interior of transport rays.

Proposition 3.3.

Let u Lip 1 ( H n , d SR ) and x , y H n , x y such that u ( x ) - u ( y ) = d SR ( x , y ) . Let σ : [ 0 , 1 ] H n be a geodesic between x and y, starting from x. Then u is Pansu differentiable at σ ( t ) , for all t ] 0 , 1 [ , and

(3.6) H u ( σ ( t ) ) = - σ ˙ ( t ) | σ ˙ ( t ) | H .

For the proof of this result we need the following two lemmas. The first one collects a differentiability property of the Carnot–Carathéodory distance function from a fixed point y n . Let us denote L y := y L .

Lemma 3.4.

Let y H n . Then the function d y ( ) := d SR ( , y ) is of class C in the euclidean sense on H n L y . In particular, d y is Pansu differentiable on H n L y .

Moreover, if x H n L y and σ : [ 0 , 1 ] H n is the geodesic between x and y, starting from y, then

H d y ( x ) = σ ˙ ( 1 ) | σ ˙ ( 1 ) | H .

Proof.

The core of the proof is contained in [5, Lemma 3.11].

If we denote by Φ ( χ , φ ) := y σ χ , θ ( 1 ) , where σ χ , θ : [ 0 , 1 ] n is any geodesic starting from 0 as in Theorem 2.1, then this map is a C -diffeomorphism from 2 n { 0 } × ( - 2 π , 2 π ) onto n L y in the usual euclidean sense. If z = Φ ( χ , θ ) n L y , then d y ( z ) = | χ | . Therefore d y is C in n L y , thus it Pansu differentiable in n L y . In addition, the condition d y Lip 1 ( n , d SR ) implies that | H d y ( z ) | 1 for all z n L y .

Now let x n L y and σ : [ 0 , 1 ] n be the geodesic between x and y. Since d y ( σ ( t ) ) = t d SR ( x , y ) for all t [ 0 , 1 ] , and σ ( t ) n L y for all t ] 0 , 1 ] , we can differentiate with respect to t:

(3.7) d SR ( x , y ) = d d t d y ( σ ( t ) ) = j = 1 2 n X j ( d y ( σ ( t ) ) ) σ ˙ j ( t ) | H d y ( σ ( t ) ) | H | σ ˙ ( t ) | H d SR ( x , y ) ,

where we used the fact that σ is a geodesic, then | σ ˙ ( t ) | H = d SR ( x , y ) . Hence, all the inequalities in (3.7) are equalities: in particular, we get that

H d y ( σ ( t ) ) = σ ˙ ( t ) | σ ˙ ( t ) | H for all  t ] 0 , 1 ] .

Lemma 3.5.

Let f , g , h : H n R three functions such that f ( x ) g ( x ) h ( x ) for all x H n . Let y H n such that f ( y ) = g ( y ) = h ( y ) and f , h are Pansu differentiable at y with H f ( y ) = H h ( y ) . Then g is Pansu differentiable in y and H g ( y ) = H f ( y ) = H h ( y ) .

Proof.

Let y n as in the hypothesis and x be any point in n . From (2.2) we have that

D H f ( y ) ( x ) = H f ( y ) , π y ( x ) H = H h ( y ) , π y ( x ) H = D H h ( y ) ( x ) .

Therefore,

f ( x ) - f ( y ) - D H f ( y ) ( y - 1 x ) g ( x ) - g ( y ) - D H f ( y ) ( y - 1 x ) = D H h ( y ) ( y - 1 x )
h ( x ) - h ( y ) - D H h ( y ) ( y - 1 x ) .

If we divide the previous inequalities by d SR ( x , y ) and we let x tend to y, by using Pansu differentiability of f and g and [33, Proposition 5.6] again, we get the thesis. ∎

Proof of Proposition 3.3.

Let t 0 ] 0 , 1 [ and a , b σ ( [ 0 , 1 ] ) such that u ( a ) > u ( σ ( t 0 ) ) > u ( b ) .

The description of geodesics given in Theorem 2.1 implies that σ ( s ) - 1 σ ( t ) L , for any s , t ] 0 , 1 [ . Therefore, it follows from Lemma 3.4 that the functions d a ( ) and d b ( ) are smooth in a neighborhood of σ ( t 0 ) . Since u Lip 1 ( n , d SR ) , it holds that

u ( a ) - u ( z ) d a ( z ) for all  z n .

Moreover, a and b lie on the geodesic between x and y, then u ( a ) = u ( b ) + d SR ( a , b ) , and hence

d b ( z ) u ( z ) - u ( b ) d SR ( a , b ) - d a ( z ) for all  z n ,

where equalities hold if z = σ ( t 0 ) , since σ is a geodesic. Moreover, from Lemma 3.4 again, it follows that

H d b ( σ ( t 0 ) ) = - H d a ( σ ( t 0 ) ) = - σ ˙ ( t 0 ) | σ ˙ ( t 0 ) | H .

Hence, from Lemma 3.5, it follows that u is Pansu differentiable in σ ( t 0 ) and

H u ( σ ( t 0 ) ) = - σ ˙ ( t 0 ) | σ ˙ ( t 0 ) | H .

4 The Beckmann problem for p = 1

4.1 Vector transport density in n

In the spirit of [45, Chapter3] we define a vector version of the horizontal transport density, introduced [23].

From now on, we fix a selection of geodesics

(4.1) S : n × n Geo ( n ) , ( x , y ) S ( x , y ) = σ x , y

that is γ-measurable for any γ Π ( μ , ν ) , see (2.10). Note that throughout the rest of this section, all results depend a priori on the measurable selection S of geodesics. However, this is not a serious issue for our purposes.

Hence, one may associate with any γ Π 1 ( μ , ν ) a positive and finite Radon measure a γ + ( n ) , called horizontal transport density, defined as

n φ 𝑑 a γ := n × n ( 0 1 φ ( σ x , y ( t ) ) | σ ˙ x , y ( t ) | H 𝑑 t ) 𝑑 γ ( x , y )

for any φ C ( n ) . If either μ 2 n + 1 , or μ 2 n + 1 , then γ ( n × n E ) = 0 and therefore the definition of a γ is independent of the choice of S.

One may also associate with any γ Π 1 ( μ , ν ) a vector version of the transport density a γ .

Definition 4.1 (Vector horizontal transport density).

For any γ Π 1 ( μ , ν ) , the vector horizontal transport density associated with γ is the H n -valued finite Radon measure 𝒘 γ ( n , H n ) defined as

(4.2) n ϕ 𝑑 𝒘 γ := n × n ( 0 1 ϕ ( σ x , y ( t ) ) , σ ˙ x , y ( t ) H 𝑑 t ) 𝑑 γ ( x , y )

for any ϕ C ( n , H n ) .

We defined the scalar transport density, resp. the vector one, in duality with C ( n ) , resp. C ( n , H n ) because they are both supported on the transport set 𝒯 , which is compact, see Theorem 3.2.

In particular, if γ ¯ Π 1 ( μ , ν ) , we have

(4.3) | 𝒘 γ ¯ | a γ ¯ ,

as measures. Therefore, it follows from the very definition of a γ ¯ , and the fact that γ ¯ is an optimal transport plan, that

(4.4) 𝒘 γ ¯ ( n , H n ) min γ Π ( μ , ν ) n × n d SR ( x , y ) 𝑑 γ ( x , y ) .

Proposition 4.2.

If u Lip ( H n , d SR ) is a Kantorovich potential, then for any γ Π 1 ( μ , ν ) , the measure | 𝐰 γ | is absolutely continuous with respect to a γ , with density - H u

𝒘 γ = - ( H u ) a γ .

Proof.

Let γ Π 1 ( μ , ν ) . Then, Proposition 3.3 implies that

σ ˙ x , y ( t ) = | σ ˙ x , y ( t ) | H σ ˙ x , y ( t ) | σ ˙ x , y ( t ) | H = - d SR ( x , y ) H u ( σ x , y ( t ) )

for every t ] 0 , 1 [ and for γ-almost every ( x , y ) (with x y otherwise both expressions vanish), where σ x , y = S ( x , y ) is the geodesic between x and y, selected by the map S in (4.1). Hence, it follows that

n ϕ 𝑑 𝒘 γ = n × n ( 0 1 - d SR ( x , y ) H u ( σ x , y ( t ) ) , ϕ ( σ x , y ( t ) ) H d t ) 𝑑 γ ( x , y )
= - 0 1 ( n × n H u ( σ x , y ( t ) ) , ϕ ( σ x , y ( t ) ) H d SR ( x , y ) 𝑑 γ ( x , y ) ) 𝑑 t
= - 0 1 ( n H u ( z ) , ϕ ( z ) H d ( S t ) # ( d SR γ ) ( z ) ) 𝑑 t

for every ϕ C ( n , H n ) , where S t = e t S . Analogously

n φ 𝑑 a γ = 0 1 ( n φ ( z ) d ( S t ) # ( d SR γ ) ( z ) ) 𝑑 t

for every φ C ( n ) . Since, by definition, 𝒘 γ and a γ are concentrated on the set of differentiability points of the Kantorovich potential u, it follows that

n ϕ 𝑑 𝒘 γ = - n H u , ϕ H 𝑑 a γ = - n ϕ d ( H u ) a γ

for every ϕ C ( n , H n ) ; hence - H u is the density of the measure | 𝒘 γ | with respect to a γ , and we write

𝒘 γ = - ( H u ) a γ .

From (4.3) and [23, Theorem 3.3 and Theorem 3.6] the next theorem easily follows.

Theorem 4.3 (Absolute continuity and summability of vector horizontal transport densities).

If one has μ L 2 n + 1 , either ν L 2 n + 1 , then there exists γ Π 1 ( μ , ν ) such that 𝐰 γ L 1 ( H n , H H n ) . Moreover, if μ L p for some p [ 1 , ] , it holds: if p < 2 n + 3 2 n + 2 , then there exists γ Π 1 ( μ , ν ) such that 𝐰 γ L p ( H n , H H n ) ; otherwise there exists γ Π 1 ( μ , ν ) such that 𝐰 γ L s ( H n , H H n ) for s < 2 n + 3 2 n + 2 .

4.2 The problem

Given any optimal transport plan γ Π 1 ( μ , ν ) we can test the vector horizontal transport density 𝒘 γ against ϕ = H φ , for any φ C ( n ) . From (2.4) it follows that

(4.5)

n H φ d 𝒘 γ = n × n ( 0 1 H φ ( σ x , y ( t ) ) , σ ˙ x , y ( t ) H 𝑑 t ) 𝑑 γ ( x , y )
= n × n ( 0 1 d d t [ φ ( σ x , y ( t ) ) ] 𝑑 t ) 𝑑 γ ( x , y )
= n × n ( φ ( y ) - φ ( x ) ) 𝑑 γ ( x , y ) = - n φ d ( μ - ν ) .

Now, given an H n -valued compactly supported Radon measure 𝒘 c ( n , H n ) we can define its horizontal divergence div H 𝒘 in duality with the smooth functions

div H 𝒘 , φ := - n H φ ( x ) 𝑑 𝒘 for all  φ C ( n ) .

With this definition in mind we can rewrite (4.5) in the following way:

div H 𝒘 γ = μ - ν for all  γ Π 1 ( μ , ν ) ,

i.e. the horizontal divergence of the measure 𝒘 γ is the signed Radon measure μ - ν , for all γ Π 1 ( μ , ν ) . This says that the following infimum is finite:

(BP) inf { 𝒘 ( n , H n ) : 𝒘 c ( n , H n ) , div H 𝒘 = μ - ν } .

This is the Heisenberg version of the well-known Beckmann problem and it will turn out to be the limit case of the wider class of problems introduced in Section 5, as explained in the introduction of the paper.

In the spirit of [45, Subsection 4.2], we show that this problem is equivalent to the Monge–Kantorovich problem (2.11).

Theorem 4.4.

The problem (BP) admits a solution. Moreover,

(4.6) ( B P ) = ( 2.11 ) ,

where

( 2.11 ) = min γ Π ( μ , ν ) n × n d SR ( x , y ) 𝑑 γ ( x , y ) ,

and a solution to (BP) can be built from a solution to (2.11), as in (4.2).

Proof.

First we prove the equality between (2.11) and (BP). We start by proving that (BP) (2.11). Take an arbitrary function φ C Lip 1 ( n , d SR ) . Pansu’s Theorem implies that H φ 1 , see [41], hence for any 𝒘 admissible we get by definition that

𝒘 ( n , H n ) n ( - H φ ) 𝑑 𝒘 = n φ d ( μ - ν ) .

Now we take φ ε = ρ ε u , where u is a Kantorovich potential and ρ ε is a mollifier for the group structure. It follows that φ ε C Lip 1 ( n , d SR ) and it converges uniformly on compact sets to the Kantorovich potential u, see [26, Proposition 2.14]. By the previous inequality and letting ε 0 we get that

𝒘 ( n , H n ) n u d ( μ - ν ) = ( 2.11 ) ,

where we used Theorem 2.2 in the last equality. Since the previous inequality holds for any admissible 𝒘 , we may take the minimum in the left hand-side and get

( B P ) ( 2.11 ) .

Now we will prove the converse inequality: for any γ Π 1 ( μ , ν ) , we know that the vector Radon measure 𝒘 γ in H n , defined in (4.2), is a compactly supported measure, which satisfies the divergence constraint thanks to (4.5). Moreover, from (4.4) we know that

𝒘 γ ( n , H n ) ( 2.11 ) .

Hence

( B P ) 𝒘 γ ( n , H n ) ( 2.11 ) .

Corollary 4.5.

If either μ L 2 n + 1 , or ν L 2 n + 1 , then there exists 𝐰 L 1 ( H n , H H n ) that solves (BP).

Proof.

The thesis follows from Theorem 4.3 and Theorem 4.4. ∎

As a consequence of Theorem 4.4, one also obtains the following dual formulation.

Corollary 4.6.

It holds that

(DP) max { n u d ( μ - ν ) : u H W 1 , , H u 1 } = ( B P ) .

Proof.

The thesis follows from Theorem 2.2 and Theorem 4.4. ∎

For any γ Π 1 ( μ , ν ) and any Kantorovich potential u Lip 1 ( n , d SR ) , the pair ( a γ , u ) solves the Monge–Kantorovich system

{ div H ( ( H u ) a γ ) = μ - ν , | H u | H = 1 a γ -a.e. , | H u | H 1 .

The second condition holds because | H u | H = 1 in the relative interior of transport rays and a γ is supported on the transport set.

4.3 Lagrangian formulation

Following [3, Lecture 9, section 3], we now introduce a Lagrangian formulation of (BP) in n . Let us consider the space C ( [ 0 , 1 ] , 2 n + 1 ) equipped with the topology of uniform convergence and let us denote by 𝒫 ( C ( [ 0 , 1 ] , 2 n + 1 ) ) the set of Borel probability measures over this space. We are interested in the problem

(4.7) inf { C ( [ 0 , 1 ] , 2 n + 1 ) ( σ ) 𝑑 Q ( σ ) : Q 𝒫 ( C ( [ 0 , 1 ] , 2 n + 1 ) ) , ( e 0 ) # Q = μ , ( e 1 ) # Q = ν } ,

where the functional is the length induced by the sub-Riemannian distance,

(4.8) ( σ ) := sup { i = 1 k d SR ( σ ( t i - 1 ) , σ ( t i ) ) : k ,  0 = t 0 < t 1 < < t k - 1 < t k = 1 } .

This definition of length makes sense even for continuous curves and it extends the sub-Riemannian one, in the sense that if σ is horizontal, then ( σ ) = SR ( σ ) . See for instance [38, Theorem 1.3.5]. It immediately follows that the functional C ( [ 0 , 1 ] , 2 n + 1 ) σ ( σ ) [ 0 , + ] is lower semicontinuous with respect to the topology of uniform convergence, hence Borel.

Hence, the following theorem holds.

Theorem 4.7.

If μ , ν P c ( H n ) , then (4.7) admits a solution.

Moreover, it holds

( 4.7 ) = min γ Π ( μ , ν ) { n × n d SR ( x , y ) 𝑑 γ ( x , y ) } < +

and Q is optimal if and only if Q is supported on the set of minimizing horizontal curves and ( e 0 , e 1 ) # Q Π 1 ( μ , ν ) .

Proof.

Let us first observe that, for any γ Π ( μ , ν ) , the measure

Q γ := S # γ 𝒫 ( C ( [ 0 , 1 ] , 2 n + 1 ) )

satisfies

( e 0 ) # Q = μ , ( e 1 ) # Q = ν ,

by definition. This ensures that the set of admissible measures Q for (4.7) is nonempty. Moreover, Q γ is supported on Geo( n ) by definition.

Let Q 𝒫 ( C ( [ 0 , 1 ] , 2 n + 1 ) ) be admissible for (4.7). Then

(4.9)

C ( [ 0 , 1 ] , 2 n + 1 ) ( σ ) 𝑑 Q ( σ ) C ( [ 0 , 1 ] , 2 n + 1 ) d SR ( σ ( 0 ) , σ ( 1 ) ) 𝑑 Q ( σ )
= n × n d SR ( x , y ) d ( e 0 , e 1 ) # Q ( x , y )
min γ Π ( μ , ν ) { n × n d SR ( x , y ) 𝑑 γ ( x , y ) } ,

where the first inequality follows from the definition of (4.8), the second equality follows from the definition of push-forward, and the last inequality follows from the fact that ( e 0 , e 1 ) # Q Π ( μ , ν ) .

To prove the converse inequality, let us consider an optimal transport plan γ ¯ Π 1 ( μ , ν ) . From the discussion at the beginning of the proof, the measure Q γ ¯ := S # γ ¯ 𝒫 ( C ( [ 0 , 1 ] , 2 n + 1 ) ) is admissible for (4.7) and supported on Geo( n ). Hence,

(4.10)

min γ Π ( μ , ν ) { n × n d SR ( x , y ) 𝑑 γ ( x , y ) } = n × n d SR ( x , y ) 𝑑 γ ¯ ( x , y )
= Geo ( n ) d SR ( σ ( 0 ) , σ ( 1 ) ) 𝑑 Q γ ¯ ( σ )
= Geo ( n ) ( σ ) 𝑑 Q γ ¯ ( σ ) .

Thus, the equality of optimization problems follows. Moreover, if Q 𝒫 ( C ( [ 0 , 1 ] , 2 n + 1 ) ) , then Q is an optimizer for (4.7) if and only if (4.9) holds with equalities, which occurs if and only if Q is supported on the set of minimizing horizontal curves and ( e 0 , e 1 ) # Q Π 1 ( μ , ν ) . ∎

This theorem is independent on the compactness of supports and on the ambient space: it holds in more general Polish metric spaces ( X , d ) with μ , ν 𝒫 ( X ) such that

X d ( x , 0 ) 𝑑 μ ( x ) + X d ( 0 , y ) ( y ) 𝑑 ν ( y ) < + ,

and is the length induced by the distance d as in (4.8).

5 The Beckmann problem for 1 < p < +

From now on we suppose Ω is a bounded domain, with regular C 1 , 1 boundary (in the Euclidean sense). Such a Ω support a q-Poincaré inequality, i.e. there exists c = c ( n , q , Ω ) > 0 such that

(5.1) Ω | φ ( x ) - φ Ω | q 𝑑 x c Ω | H φ ( x ) | q 𝑑 x for all  φ H W 1 , q ( Ω )

for any 1 q < + , where φ Ω = 1 2 n + 1 ( Ω ) Ω φ 𝑑 x . See [34, (1.26) and Theorem 1.21].

We consider the space

H := { σ AC ( [ 0 , 1 ] , Ω ¯ ) : σ  is horizontal } ,

viewed as subset of C ( [ 0 , 1 ] , Ω ¯ ) equipped with the topology of uniform convergence. We stress the fact that the regularity assumption on the boundary of Ω guarantees that any two points in Ω ¯ can be connected by an element of H, see [23, Lemma 4.1]. This assumption replaces the stronger convexity assumption typically used in the Euclidean setting, since there are no non-trivial geodesically convex set in n , see [39].

5.1 Vector traffic intensities in n

In this subsection we will see how also the congested optimal transport problem is related with some minimization problems, over a set of vector fields with divergence constraint.

Let μ , ν 𝒫 ( Ω ¯ ) be two probability measures over the closure of Ω and let us consider the space C ( [ 0 , 1 ] , Ω ¯ ) equipped with the topology of uniform convergence. A horizontal traffic plan between μ and ν is a Borel probability measure Q 𝒫 ( C ( [ 0 , 1 ] , Ω ¯ ) ) such that

  1. C ( [ 0 , 1 ] , Ω ¯ ) ( σ ) 𝑑 Q ( σ ) < + ,

  2. Q is concentrated on the set H,

  3. ( e 0 ) # Q = μ and ( e 1 ) # Q = ν .

We denote by

𝒬 H ( μ , ν ) := { Q 𝒫 ( C ( [ 0 , 1 ] , Ω ¯ ) ) : Q  horizontal traffic plan between  μ  and  ν } .

As for the transport plans, one can also associate with any horizontal traffic plan Q 𝒬 H ( μ , ν ) both a scalar and a vector measures. The scalar measure is the horizontal traffic intensity i Q + ( Ω ¯ ) defined as

Ω ¯ φ ( x ) 𝑑 i Q ( x ) := C ( [ 0 , 1 ] , Ω ¯ ) ( 0 1 φ ( σ ( t ) ) | σ ˙ ( t ) | H 𝑑 t ) 𝑑 Q ( σ ) for all  φ C ( Ω ¯ ) .

Example 5.1.

Let us consider the two measures

μ := δ x and ν := δ y for some  x n  and some  y n L x ,

and let σ : [ 0 , 1 ] n be the unique geodesic between x and y, according to Theorem 2.1. Let us consider a bounded domain Ω, with regular C 1 , 1 boundary, such that x , y Ω and σ ( [ 0 , 1 ] ) Ω . Then the measure

Q := δ σ 𝒬 H ( μ , ν ) ,

and its horizontal traffic intensity is precisely

i Q = 1 σ ( [ 0 , 1 ] ) .

As for the vector measure, we introduce now the following definition.

Definition 5.2 (Vector horizontal traffic intensity).

Let Q 𝒬 H ( μ , ν ) be a horizontal traffic plan admissible between μ and ν. One can associate with Q an H Ω -valued finite Radon measure 𝒘 Q ( Ω ¯ , H Ω ¯ ) defined as

Ω ¯ ϕ ( x ) 𝑑 𝒘 Q := C ( [ 0 , 1 ] , Ω ¯ ) ( 0 1 ϕ ( σ ( t ) ) , σ ˙ ( t ) H 𝑑 t ) 𝑑 Q ( σ )

for any continuous horizontal vector field ϕ C ( Ω ¯ , H Ω ¯ ) . We will call this measure vector horizontal traffic intensity induced by Q .

The following lemma holds.

Lemma 5.3.

Let Q Q H ( μ , ν ) , then

  1. 𝒘 Q ( Ω ¯ , H Ω ¯ ) C ( [ 0 , 1 ] , Ω ¯ ) ( σ ) 𝑑 Q ( σ ) < + ,

  2. | 𝒘 Q | i Q as measures,

  3. div H 𝒘 Q = μ - ν .

Proof.

(1) This follows from the definition of total variation norm and of horizontal traffic plan, while (2) follows from the following inequality:

| Ω ¯ ϕ 𝑑 𝒘 | Ω ¯ | ϕ | H 𝑑 i Q for all  ϕ C ( Ω ¯ , H Ω ¯ ) .

(3) By definition, 𝒘 Q is compactly supported. Then its horizontal divergence is well defined: let φ C ( n ) , it holds

Ω ¯ H φ d 𝒘 Q = C ( [ 0 , 1 ] , Ω ¯ ) ( 0 1 H φ ( σ ( t ) ) , σ ˙ ( t ) H 𝑑 t ) 𝑑 Q ( σ )
= C ( [ 0 , 1 ] , Ω ¯ ) ( 0 1 d d t ( φ ( σ ( t ) ) ) 𝑑 t ) 𝑑 Q ( σ )
= C ( [ 0 , 1 ] , Ω ¯ ) [ φ ( σ ( 1 ) ) - φ ( σ ( 0 ) ) ] 𝑑 Q ( σ )
= Ω ¯ φ d ( ( e 1 ) # Q - ( e 0 ) # Q ) = - Ω ¯ φ d ( μ - ν ) ,

and the thesis follows. ∎

Example 5.4.

Following [14, Example 4.1], we consider x n , y n L x , σ and Ω as in Example 5.1. Moreover, we consider the two discrete measures

μ = ν = δ x + δ y 2

and the horizontal traffic plan

Q := 1 2 δ σ + 1 2 δ σ ~ 𝒬 H ( μ , ν ) ,

where σ ~ ( t ) := σ ( 1 - t ) . Again

i Q = 1 σ ( [ 0 , 1 ] )

but

𝒘 Q 0 .

This shows that, in the previous Lemma, (2) is not an equality in general, because the curves of Q may produce some cancellations: this is due to the fact that 𝒘 Q takes into account the orientation of the curves, while i Q does not.

Lemma 5.5.

Let 1 < p < + . If

𝒬 H p ( μ , ν ) := { Q 𝒬 H ( μ , ν ) : i Q L p ( Ω ) } ,

then

{ 𝒘 L p ( Ω , H Ω ) : div H 𝒘 = μ - ν } .

Moreover, μ - ν ( H W 1 , q ( Ω ) ) , with q = p p - 1 .

Proof.

Let Q 𝒬 H p ( μ , ν ) . Lemma 5.3 implies that 𝒘 Q L p ( Ω , H Ω ) and div H 𝒘 Q = μ - ν . Moreover,

+ > 𝒘 Q L p ( Ω , H Ω ) = sup { Ω ϕ , 𝒘 Q H 𝑑 x : ϕ L q ( Ω , H Ω ) , ϕ L q ( Ω , H Ω ) 1 }
sup { Ω H φ , 𝒘 Q H 𝑑 x : φ H W 1 , q ( Ω ) , φ H W 1 , q ( Ω ) 1 }
= sup { Ω φ d ( μ - ν ) : φ H W 1 , q ( Ω ) , φ H W 1 , q ( Ω ) 1 } = μ - ν ( H W 1 , q ( Ω ) ) .

As a consequence, we see that if μ , ν are Dirac masses then 𝒬 H p ( μ , ν ) = whenever p N N - 1 . Indeed, this is so because μ - ν ( H W 1 , q ( Ω ) ) only for q > N , since otherwise H W 1 , q ( Ω ) contains non-continuous elements.

5.2 Well-posedness of the problem

In the beginning of this subsection, we want to determine the image of the horizontal divergence operator div H on L p ( Ω , H Ω ) . To do this, we start with the following fact.

Lemma 5.6.

If 1 < p < , and q = p p - 1 , then the horizontal divergence induces a continuous linear operator div H : L p ( Ω , H Ω ) ( H W 1 , q ( Ω ) ) with norm 1 .

Proof.

Any given 𝒘 L p ( Ω , H Ω ) induces in a natural way an H n -valued compactly supported Radon measure 𝒘 c ( n , H n ) . Hence the horizontal divergence div H 𝒘 is well defined,

div H 𝒘 , φ = - Ω H φ d 𝒘 for all  φ C ( n ) .

Moreover, since 𝒘 L p ( Ω , H Ω ) , we easily see that

| div H 𝒘 , φ | 𝒘 L p ( Ω , H Ω ) φ H W 1 , q ( Ω ) for all  φ H W 1 , q ( Ω ) .

This allows us to extend the action of div H 𝒘 on any testing function φ H W 1 , q ( Ω ) , still having the above inequality. In particular, it tells us that div H has norm 1 as claimed. ∎

The operator div H : L p ( Ω , H Ω ) ( H W 1 , q ( Ω ) ) is not surjective. Indeed, for any 𝒘 L p ( Ω , H Ω ) we see that

div H 𝒘 , 1 = 0 ,

while there are many other elements in ( H W 1 , q ( Ω ) ) whose action against 1 is not 0 (for instance, any strictly positive function in L p ( Ω ) ). This motivates us to introduce a particular closed subspace ( H W 1 , q ) ( Ω ) of ( H W 1 , q ( Ω ) ) ,

( H W 1 , q ) ( Ω ) := { f ( H W 1 , q ( Ω ) ) : f , 1 = 0 } .

We endow ( H W 1 , q ) ( Ω ) with the environmental dual norm ( H W 1 , q ( Ω ) ) . Following [46, Section 2] one can characterize such space.

Proposition 5.7.

It holds that

( H W 1 , q ) ( Ω ) = div H ( L p ( Ω , H Ω ) ) .

Moreover, there is a constant c = c ( n , q , Ω ) 1 such that for every f ( H W 1 , q ) ( Ω ) , there is a horizontal vector field 𝐰 L p ( Ω , H Ω ) such that

div H 𝒘 = f ,

and 𝐰 L p ( Ω , H Ω ) c f ( H W 1 , q ( Ω ) ) .

Proof.

Let us consider a horizontal vector field 𝒘 L p ( Ω , H Ω ) . We know that its weak divergence defines a distribution div H 𝒘 ( H W 1 , q ( Ω ) ) and it is trivial to prove that div H 𝒘 , 1 = 0 . Therefore,

div H 𝒘 ( H W 1 , q ) ( Ω ) .

Conversely, let us consider the following minimization problem

min φ H W 1 , q ( Ω ) ( φ ) ,

where

φ ( φ ) := 1 q Ω | H φ | H q 𝑑 x + f , φ .

This problem admits at least one solution. Indeed, we can restrict the minimization to the subspace

(5.2) H W 1 , q ( Ω ) := { φ H W 1 , q ( Ω ) : Ω φ ( x ) 𝑑 x = 0 } ,

which is a convex subset and closed with respect to weak topology of H W 1 , q ( Ω ) . The functional is lower semicontinuous with respect to the weak topology of H W 1 , q ( Ω ) and it is coercive on H W 1 , q ( Ω ) : indeed, if φ H W 1 , q ( Ω ) , then from (5.1) it follows the existence of a constant c = c ( n , q , Ω ) > 0 such that

| ( φ ) | 1 q H φ L q ( Ω , H Ω ) q - ( c 1 q + 1 ) f ( H W 1 , q ( Ω ) ) H φ L q ( Ω , H Ω )
= H φ L q ( Ω , H Ω ) ( 1 q H φ L q ( Ω , H Ω ) q - 1 - ( c 1 q + 1 ) f ( H W 1 , q ( Ω ) ) ) .

Hence, the existence of minimizers follows again from (5.1).

Computing Euler–Lagrange equations, a solution φ of the problem above turns out to satisfy

- Ω | H φ | q - 2 H φ , H ψ H 𝑑 x = f , ψ for all  ψ H W 1 , q .

This means that there exists 𝒘 = | H φ | q - 2 H φ L p ( Ω , H Ω ) such that div H 𝒘 = f . Moreover, testing f against φ we get

𝒘 L p ( Ω , H Ω ) p = Ω | 𝒘 | H p 𝑑 x = Ω | H φ | H q 𝑑 x = - f , φ
f ( H W 1 , q ( Ω ) ) φ H W 1 , q ( Ω ) ( c 1 q + 1 ) f ( H W 1 , q ( Ω ) ) H φ L q ( Ω , H Ω )
= ( c 1 q + 1 ) f ( H W 1 , q ( Ω ) ) 𝒘 L p ( Ω , H Ω ) p - 1 ,

as desired. ∎

Let now consider a function 𝒢 : Ω ¯ × 2 n such that

x 𝒢 ( x , w ) is Lebesgue-measurable for any  w 2 n .

Let us suppose in addition that 𝒢 ( x , 0 ) = 0 for any x Ω ¯ ,

(5.3) w 𝒢 ( x , w ) is convex for any  x Ω ¯

and it satisfies

(5.4) a p | w | p - h 0 ( x ) 𝒢 ( x , w ) b p | w | p + h 1 ( x ) for any  ( x , w ) Ω ¯ × 2 n

for some h 0 , h 1 L 1 ( Ω ) , h 0 , h 1 0 , some constants a , b > 0 and p ( 1 , ) .

Following [46], we are interested in the Beckmann-type problem

(${\mathcal{B}}$) inf 𝒘 L p ( Ω , H Ω ) { Ω 𝒢 ( x , 𝒘 ( x ) ) 𝑑 x : div H 𝒘 = f } .

Since div H 𝒘 ( H W 1 , q ) ( Ω ) for any 𝒘 L p ( Ω , H Ω ) , the natural functional setting for the above problem is f ( H W 1 , q ) ( Ω ) .

Theorem 5.8.

If f ( H W 1 , q ) ( Ω ) , then ((${\mathcal{B}}$)) admits a minimizer. If G is strictly convex, then the minimizer is unique.

Proof.

Let f ( H W 1 , q ) ( Ω ) , then Proposition 5.7 implies that there exists at least one vector field 𝒘 L p ( Ω , H Ω ) such that div H 𝒘 = f : hence (5.4) implies that ( ) < + . Let ( 𝒘 n ) n L p ( Ω , H Ω ) be a minimizing sequence. As we are dealing with a minimizing sequence, from (5.4) it follows that this sequence is bounded in L p ( Ω , H Ω ) : hence, up to a subsequence, it is weakly convergent to some 𝒘 ~ L p ( Ω , H Ω ) . This vector field is still admissible: indeed, by the weak convergence it follows that

- Ω 𝒘 ~ , H φ H 𝑑 x = - lim n Ω 𝒘 n , H φ H 𝑑 x = f , φ for all  φ H W 1 , q ( Ω ) .

Since 𝒢 is convex in the second variable, the integral functional is convex, as well, and lower semicontinuous; this implies the lower semicontinuity of the integral functional with respect to the weak convergence in L p ( Ω , H Ω ) , hence

Ω 𝒢 ( x , 𝒘 ~ ( x ) ) 𝑑 x lim inf n + Ω 𝒢 ( x , 𝒘 n ( x ) ) 𝑑 x = ( ) < + .

This proves that 𝒘 ~ is a minimum. ∎

5.3 Dual formulation

The aim of this subsection is to show that the Beckmann-type problem ((${\mathcal{B}}$)) admits the dual formulation

(${\mathcal{D}}$) sup φ H W 1 , q ( Ω ) { - f , φ - Ω 𝒢 * ( x , H φ ( x ) ) 𝑑 x } ,

where f ( H W 1 , q ) ( Ω ) , 𝒢 * is the Legendre transform of the map w 𝒢 ( x , w ) (see Theorem 5.9), and again we are identifying H φ = i = 1 2 n X i φ X i with ( X 1 φ , , X 2 n φ ) : Ω 2 n .

By definition the function

x 𝒢 * ( x , z ) is Lebesgue-measurable for any  z 2 n

and

(5.5) 2 n z 𝒢 * ( x , z ) is convex, for any  x Ω ¯ .

The condition (5.4) implies the following q-growth condition on 𝒢 * :

(5.6) 1 q b q - 1 | z | q - h 1 ( x ) 𝒢 * ( x , z ) 1 q a q - 1 | z | q + h 0 ( x ) for any  ( x , z ) Ω ¯ × 2 n

where q = p p - 1 .

We follow the proof developed in [13], which is based on the following classical theorem of convex analysis.

Theorem 5.9 (Convex duality, [28, Proposition 5]).

Let F : Y R a convex and lower semicontinuous functional on a reflexive Banach space Y. Let X be another reflexive Banach space and A : X Y a bounded linear operator, with adjoint operator A : Y X . Then we have

(5.7) sup x X x , x - ( A x ) = inf y Y { * ( y ) : A y = x } , x X ,

where F * : Y R denotes the Legendre–Fenchel transform of F .

Moreover, if the supremum in (5.7) is attained at some x 0 X , then infimum in (5.7) is attained at some y 0 Y such that

(5.8) y 0 ( A x 0 ) ,

where F denotes the subdifferential

( y ) := { y Y : ( y ) + * ( y ) = y , y } , y Y .

We are now ready to show that ((${\mathcal{B}}$)) = ((${\mathcal{D}}$)).

Theorem 5.10.

If f ( H W 1 , q ) ( Ω ) , then the problem ((${\mathcal{D}}$)) admits a solution and

min 𝒘 L p ( Ω , H Ω ) { Ω 𝒢 ( x , 𝒘 ( x ) ) 𝑑 x : div H 𝒘 = f } = max φ H W 1 , q ( Ω ) { - Ω 𝒢 * ( x , H φ ( x ) ) 𝑑 x - f , φ } .

Moreover, if 𝐰 0 L p ( Ω , H Ω ) in a minimizer for ((${\mathcal{B}}$)) and φ 0 H W 1 , q ( Ω ) is a maximizer for ((${\mathcal{D}}$)), then we have the following primal-dual optimality condition

(5.9) 𝒘 0 𝒢 * ( x , H φ 0 ( x ) ) a.e. in  Ω .

Proof.

First we observe that the existence of a minimizer for the problem ((${\mathcal{B}}$)) follows from Theorem 5.8. Thanks to the Direct method in calculus of variations also ((${\mathcal{D}}$)) admits at least a solution belonging the space H W 1 , q ( Ω ) of Sobolev functions with zero mean. Indeed, (5.5) implies that the functional ~ : H W 1 , q ( Ω )

~ ( φ ) = Ω 𝒢 * ( x , H φ ( x ) ) 𝑑 x + f , φ

is lower semicontinuous with respect to the weak topology of H W 1 , q ( Ω ) ; moreover it is coercive on H W 1 , q ( Ω ) : let φ H W 1 , q ( Ω ) from (5.6) and (5.1)

~ ( φ ) 1 q b q - 1 H φ L q ( Ω , H Ω ) q - f ( H W 1 , q ( Ω ) ) φ H W 1 , q ( Ω ) - d
1 q b q - 1 H φ L q ( Ω , H Ω ) q - ( c 1 q + 1 ) f ( H W 1 , q ( Ω ) ) H φ L q ( Ω , H Ω ) - d
H φ L q ( Ω , H Ω ) ( 1 q b q - 1 H φ L q ( Ω , H Ω ) q - 1 - ( c 1 q + 1 ) f ( H W 1 , q ( Ω ) ) ) - d .

where d > 0 and c = c ( n , q , Ω ) is the constant appearing in (5.1). The existence of solutions follows from again from (5.1).

Let now denote by X = H W 1 , q ( Ω ) , Y = L q ( Ω , H Ω ) and let us consider the operator

A : X Y , A ( φ ) = H φ for all  φ X .

This operator is linear; moreover, it is bounded since

A ( φ ) Y = H φ L q ( Ω , H Ω ) φ H W 1 , q ( Ω ) .

We denote by : Y the functional

( ϕ ) := Ω 𝒢 * ( x , ϕ ( x ) ) 𝑑 x .

By the convexity and the lower semicontinuity of 𝒢 it follows that * : L p ( Ω , H Ω )

* ( 𝒘 ) = Ω 𝒢 * * ( x , 𝒘 ( x ) ) 𝑑 x = Ω 𝒢 ( x , 𝒘 ( x ) ) 𝑑 x .

Moreover, it is a straightforward consequence of the definition of the continuous linear operator

div H : L p ( Ω , H Ω ) ( H W 1 , q ( Ω ) )

that A = - div H .

The rest of the thesis follows from the primal-dual optimality condition (5.8). ∎

Remark 5.11.

In [46] the author provides a different proof of the equality

( m a t h c a l B ) = ( m a t h c a l D ) ,

which works also in the Heisenberg group. For the sake of completeness we just highlight the main steps. Let : ( H W 1 , q ( Ω ) ) be the functional

( k ) := min 𝒘 L p ( Ω , H Ω ) { Ω 𝒢 ( x , 𝒘 ( x ) ) 𝑑 x : div H 𝒘 = f + k } .

Lemma 5.7 and (5.4) imply that the functional ( k ) is well-defined and real-valued if and only if k ( H W 1 , q ) ( Ω ) . A computation implies that its Legendre transform * : H W 1 , q ( Ω ) is precisely

* ( φ ) = - f , φ + Ω 𝒢 * ( x , - H φ ( x ) ) 𝑑 x .

By definition, * * ( 0 ) = sup φ H W 1 , q ( Ω ) { - * ( φ ) } ; again, since one can restrict the minimization to H W 1 , q ( Ω ) , hence sup - * < + . By taking the sup on - φ instead of φ we also have

* * ( 0 ) = sup φ H W 1 , q ( Ω ) { - f , φ - Ω 𝒢 * ( x , H φ ( x ) ) 𝑑 x } .

One can prove that is convex and l.s.c., then * * ( 0 ) = ( 0 ) and the thesis follows.

Let us remark that if the function w 𝒢 ( x , w ) is strictly convex, then the map z 𝒢 * ( x , z ) is C 1 (see [47, Proposition 4.26]). Therefore, its subdifferential reduces to the singleton D 𝒢 * and, according to (5.9) in Theorem 5.10, the unique solution to ((${\mathcal{B}}$)) can be written as

𝒘 0 = D 𝒢 * ( x , H φ 0 ( x ) ) a.e. in  Ω ,

where φ 0 is a weak solution to the Euler–Lagrange equation

(5.10) div H ( D 𝒢 * ( x , H φ ( x ) ) ) = f .

Example 5.12.

If 𝒢 ( w ) = 1 p | w | p , then 𝒢 * ( z ) = 1 q | z | q and (5.10) becomes the degenerate q-Laplace equation

(5.11) div H ( | H φ | H q - 2 H φ ) = f in  Ω .

Based on the results in [16] and [48, Theorem 3.1], it can be shown that if f , X 2 n + 1 f L loc 2 ( Ω ) , and φ 0 H W 1 , q ( Ω ) is a weak solution to (5.11), then

F := | H φ 0 | q - 2 2 H φ 0 H W loc 1 , 2 ( Ω , H Ω ) .

Consequently, following [13], the solution 𝒘 0 to ((${\mathcal{B}}$)) satisfies

𝒘 0 = D 𝒢 * ( H φ 0 ) = | F | q - 2 q F H W loc 1 , r ( Ω , H Ω ) ,

where

r := { 2 if  q = 2 , N q ( N - 1 ) q + 2 - N if  2 < q < 4 .

On the other hand, the optimal regularity for solutions is the Hölder continuity of the horizontal derivatives: for 2 q < , it has been established in [40], with f L loc p ( Ω ) and p > N .

See also [24] for Lipschitz regularity for solutions to (5.10), with a different choice of 𝒢 and f = 0 , and [22] for its parabolic counterpart.

5.4 Lagrangian formulation

The congested optimal transport problem, adapted to the Heisenberg group in [23], consists of the following minimization problem:

(${\mathcal{W}}$) inf Q 𝒬 H p ( μ , ν ) Ω G ( x , i Q ( x ) ) 𝑑 x ,

where the function G : Ω ¯ × satisfies:

(5.12) x G ( x , i ) is Lebesgue-measurable for any  i ,

and

(5.13) G ( x , 0 ) = 0 for any  x Ω ¯ .

Additionally, we assume:

(5.14) i G ( x , i )  is convex and non-decreasing for any  x Ω ¯ ,

and it exhibits p-growth in the second variable:

(5.15) a p i p - h 0 ( x ) G ( x , i ) b p i p + h 1 ( x ) for any  ( x , i ) Ω ¯ × ,

for some p ( 1 , ) , where a , b > 0 and h 0 , h 1 L p ( Ω , + ) . It turns out that solutions to ((${\mathcal{W}}$)) can, in some cases, exhibit the nature of Wardrop equilibria. For further details, we refer the interested reader to [23].

In this subsection, we investigate the equivalence between ((${\mathcal{W}}$)) and ((${\mathcal{B}}$)), where f = μ - ν , and the function 𝒢 ( x , ) = G ( x , | | H ) for any x Ω ¯ :

inf 𝒘 L p ( Ω , H Ω ) { Ω G ( x , | 𝒘 ( x ) | H ) 𝑑 x : div H 𝒘 = μ - ν } .

Moreover, we will demonstrate how to pass from a solution of one problem to a solution of the other, and vice versa.

First, we state the following result.

Proposition 5.13.

If Q H p ( μ , ν ) , then μ - ν ( H W 1 , q ) ( Ω ) . More precisely, for every Q Q H p ( μ , ν ) , 𝐰 Q is admissible for ((${\mathcal{B}}$)) and | 𝐰 Q | H i Q . Therefore,

(5.16) ( ) ( 𝒲 ) .

Proof.

The result follows from Lemma 5.5, Proposition 5.7 and the monotonicity of G. ∎

Hence, the equivalence between the two problems is reduced to the converse inequality.

To prove this, we follow [12, Theorem 2.1] and [44], which are based on the Dacorogna–Moser scheme combined with a regularization procedure. In the Euclidean setting, this regularization procedure is unnecessary if μ and ν are sufficiently smooth, because every solution 𝒘 to ((${\mathcal{B}}$)) belongs to some Sobolev space. Hence, one can define the flow in the weaker sense of DiPerna-Lions, as seen in [13]. The absence of a sub-Riemannian DiPerna–Lions theory, as discussed in [6], is one of the reasons to force a regularization in our setting.

Theorem 5.14.

If μ - ν ( H W 1 , q ) ( Ω ) , then Q H p ( μ , ν ) . More precisely, for every 𝐰 admissible for ((${\mathcal{B}}$)), there is Q 𝐰 Q H p ( μ , ν ) such that i Q 𝐰 | 𝐰 | H . Therefore,

( ) ( 𝒲 ) .

Proof.

Let 𝒘 = j = 1 2 n 𝒘 j X j L p ( Ω , H Ω ) be a solution to ((${\mathcal{B}}$)). We extend it by 0 outside Ω and we consider the horizontal vector field 𝒘 ε = j = 1 2 n 𝒘 j ε X j , where

𝒘 j ε ( x ) := ρ ε 𝒘 j ( x ) for all  j = 1 , , 2 n ,

ρ ε is a mollifier for the group structure of n and Ω ε := B ( 0 , ε ) Ω . We notice that, for any j = 1 , , 2 n , 𝒘 j ε C c ( n ) with support in Ω ε . A simple computation implies that

div H 𝒘 ε = μ ε - ν ε ,

where

μ ε = ρ ε μ + ε , ν ε = ρ ε ν + ε C ,

and we are supposing that both μ and ν are extended by 0 outside Ω. Let us introduce the non-autonomous horizontal vector field

𝒘 ^ ε ( t , x ) := 𝒘 ε ( x ) μ ¯ ε ( t , x ) for all  ( t , x ) [ 0 , 1 ] × Ω ¯ ε ,

where

(5.17) μ ¯ ε ( t , x ) := ( 1 - t ) μ ε ( x ) + t ν ε ( x ) > ε > 0 , ( t , x ) [ 0 , 1 ] × Ω ¯ ε .

Let us notice that 𝒘 ^ ε ( t , ) C c ( n , H n ) with support in Ω ε , for any t [ 0 , 1 ] . Moreover, for any ϵ > 0 it holds that

μ ε - ν ε = div H 𝒘 ε = div ϵ 𝒘 ε ,

where div ϵ denotes the Riemannian ϵ-divergence in n 2 n + 1 , namely

φ div ϵ 𝒘 ε d x = - g ϵ ( 𝒘 ε , ϵ φ ) 𝑑 x

for any φ C ( n ) , see also Section 2.1.4. In particular, by noticing that 𝒘 ^ ε ( t , ) C c ( 2 n + 1 , 2 n + 1 ) , for any t [ 0 , 1 ] , the curve μ ¯ ε ( t , ) satisfies the following initial value problem for the Riemannian continuity equation

(5.18) { t λ + div ϵ ( 𝒘 ^ ε λ ) = 0 , λ ( 0 , ) = μ ε .

Since 𝒘 ^ ε ( t , ) is smooth for any t [ 0 , 1 ] , the unique solution to (5.18) is given by ( Ψ ε ( t , ) ) # μ ε , where

Ψ ε : [ 0 , 1 ] × 2 n + 1 2 n + 1

is the flow of the vector field 𝒘 ^ ε in the Riemannian sense,

{ d d t Ψ ε ( t , x ) = 𝒘 ^ ε ( t , Ψ ε ( t , x ) ) , Ψ ε ( 0 , x ) = x .

Let us remark that, since 𝒘 ^ ε ( t , ) has compact support in Ω ε , for any t [ 0 , 1 ] , then Ψ ε ( t , x ) Ω ε if x Ω ε , and Ψ ε ( t , x ) = x otherwise.

Hence

( Ψ ε ( t , ) ) # μ ε = μ ¯ ε ( t , ) for all  t [ 0 , 1 ] .

If we denote by

Φ ε : Ω ¯ ε C ( [ 0 , 1 ] , Ω ¯ ε ) ,

the trajectory map, i.e. Φ ε ( x ) := Ψ ε ( , x ) , then the measure

Q ε := ( Φ ε ) # μ ε + ( C ( [ 0 , 1 ] , Ω ¯ ε ) )

satisfies

( e t ) # Q ε = ( e t Φ ε ) # μ ε = ( Ψ ε ( t , ) ) # μ ε = μ ¯ ε ( t , ) for all  t [ 0 , 1 ] ;

in particular,

(5.19) ( e 0 ) # Q ε = μ ε , ( e 1 ) # Q ε = ν ε .

By construction, Q ε is supported on the integral curves of the horizontal vector field 𝒘 ^ ε and the definition of Q ε = Φ # ε μ ε and a change of variables imply

C ( [ 0 , 1 ] , Ω ¯ ε ) SR ( σ ) 𝑑 Q ε ( σ ) = Ω ¯ ε ( 0 1 | 𝒘 ^ ε ( t , Ψ ε ( t , x ) ) | H 𝑑 t ) μ ε ( x ) 𝑑 x
= 0 1 ( Ω ¯ ε | 𝒘 ^ ε ( t , y ) | H 𝑑 Ψ ε ( t , ) # μ ε ( y ) ) 𝑑 t = 0 1 ( Ω ¯ ε | 𝒘 ^ ε ( t , y ) | H μ ¯ ε ( t , y ) 𝑑 y ) 𝑑 t
= Ω ¯ ε | 𝒘 ε ( y ) | H 𝑑 y 𝒘 L 1 ( Ω , H Ω ) 2 n + 1 ( Ω ) 1 q 𝒘 L p ( Ω , H Ω ) < + ,

because 𝒘 L p ( Ω , H Ω ) .

Hence, one can define the horizontal traffic intensity i Q ε associated with Q ε as

Ω ¯ ε φ ( x ) 𝑑 i Q ε ( x ) := C ( [ 0 , 1 ] , Ω ¯ ε ) ( 0 1 φ ( σ ( t ) ) | σ ˙ ( t ) | H 𝑑 t ) 𝑑 Q ε ( σ )

for any φ C ( Ω ¯ ε ) . Moreover, the definition of Q ε and a change of variables imply

(5.20)

Ω ¯ ε φ ( x ) 𝑑 i Q ε ( x ) = Ω ¯ ε ( 0 1 φ ( Ψ ε ( t , x ) ) | 𝒘 ^ ε ( t , Ψ ε ( t , x ) ) | H 𝑑 t ) μ ε ( x ) 𝑑 x
= 0 1 ( Ω ¯ ε φ ( y ) | 𝒘 ^ ε ( t , y ) | H μ ¯ ε ( t , y ) 𝑑 y ) 𝑑 t
= Ω ¯ ε φ ( y ) | 𝒘 ε ( y ) | H 𝑑 y

for any φ C ( Ω ¯ ε ) ; hence i Q ε = | 𝒘 ε | H C c ( n ) , with support in Ω ε , and again

(5.21) i Q ε L 1 ( Ω ε ) = 𝒘 ε L 1 ( Ω ε , H Ω ε ) 2 n + 1 ( Ω ) 1 q 𝒘 L p ( Ω , H Ω ) < + .

Analogously we define and compute the vector traffic intensity 𝒘 Q ε associated with Q ε : given ϕ C ( Ω ¯ ε , H Ω ¯ ε )

Ω ¯ ε ϕ ( x ) 𝑑 𝒘 Q ε ( x ) := Ω ¯ ε ( 0 1 ϕ ( Ψ ε ( t , x ) ) , 𝒘 ^ ε ( t , Ψ ε ( t , x ) ) H 𝑑 t ) μ ε ( x ) 𝑑 x = Ω ¯ ε ϕ ( y ) , 𝒘 ε ( y ) H 𝑑 y ,

hence 𝒘 Q ε = 𝒘 ε C c ( n , H n ) , with support in Ω ε .

For any 0 < ε < 1 , we denote by Q ~ ε := Q ε 1 + ε and we fix a compact set Ω n such that Ω Ω ε Ω . Moreover, it holds that

Q ~ ε ( C ( [ 0 , 1 ] , Ω ) ) = Q ε ( C ( [ 0 , 1 ] , Ω ) ) 1 + ε = 1 1 + ε Ω μ ε ( x ) 𝑑 x = 1 1 + ε Ω ε μ ε ( x ) 𝑑 x = 1 ,

therefore Q ~ ε 𝒫 ( C ( [ 0 , 1 ] , Ω ) ) . Since both i Q ~ ε and 𝒘 Q ~ ε are invariant by reparametrization, we may suppose that Q ~ ε is supported on curves parametrized with constant speed, for any 0 < ε < 1 . The sequence ( Q ~ ε ) 0 < ε < 1 𝒫 ( C ( [ 0 , 1 ] , Ω ) ) is tight: indeed, if we denote by

H ( [ 0 , 1 ] , Ω ) := { σ : [ 0 , 1 ] Ω : σ  is horizontal }

and we consider K > 0 , then the set { σ H ( [ 0 , 1 ] , Ω ) : | σ ˙ | H K } is compact by the Ascoli–Arzelà Theorem. Moreover,

Q ~ ε ( { σ H ( [ 0 , 1 ] , Ω ) : | σ ˙ | H > K } ) = Q ~ ε ( { σ H ( [ 0 , 1 ] , Ω ) : SR ( σ ) > K } )
1 ( 1 + ε ) K Ω i Q ε ( x ) 𝑑 x c K 𝒘 L p ( Ω , H Ω ) ,

with c > 0 independent of ε thanks to (5.21). Therefore, ( Q ~ ε ) 0 < ε < 1 is tight, and thus admits a (not relabeled) subsequence that weakly converges to some Q 𝒫 ( C ( [ 0 , 1 ] , Ω ) ) . It is obvious that Q is concentrated on curves valued in Ω ¯ , therefore Q 𝒫 ( C ( [ 0 , 1 ] , Ω ¯ ) ) . Moreover, following the same proof of [18, Lemma 2.8], one can prove that Q is concentrated on H.

From the fact that Q ~ ε Q , as ε 0 , and (5.19) we can conclude that

Q 𝒬 H ( μ , ν ) .

In the end, i Q ~ ε = | 𝒘 ε | H 1 + ε converges to | 𝒘 | H in L p ( Ω ) , hence it converges to the same limit in + ( Ω ) . Then, for any φ C ( Ω ) , φ 0 ,

Ω φ ( x ) | 𝒘 ( x ) | H 𝑑 x = lim ε 0 1 1 + ε Ω φ ( x ) | 𝒘 ε ( x ) | H 𝑑 x
= lim ε 0 C ( [ 0 , 1 ] , Ω ) ( 0 1 φ ( σ ( t ) ) | σ ˙ ( t ) | H 𝑑 t ) 𝑑 Q ~ ε ( σ )
C ( [ 0 , 1 ] , Ω ) ( 0 1 φ ( σ ( t ) ) | σ ˙ ( t ) | H 𝑑 t ) 𝑑 Q ( σ ) = Ω φ ( x ) 𝑑 i Q ( x ) ,

by the lower semicontinuity of the map Q ( φ ( σ ( t ) ) | σ ˙ ( t ) | H 𝑑 t ) 𝑑 Q ( σ ) with respect to the weak convergence of measures. In particular, this means that

i Q | 𝒘 | H L p ( Ω ) ,

hence

Q 𝒬 H p ( μ , ν ) .

By using the monotonicity of G,

( 𝒲 ) Ω G ( x , i Q ( x ) ) 𝑑 x Ω 𝒢 ( x , 𝒘 ( x ) ) 𝑑 x = ( ) .

An alternative proof of the previous theorem can be found in the literature, specifically in [13, Theorem 3.2], and it is also applicable in n . This proof relies on the well-known superposition principle and requires additional assumptions on the measures μ and ν, namely that they are p-summable and bounded below. This approach does not require going through the regularization procedure of solutions to ((${\mathcal{B}}$)). More precisely, let us suppose that

μ , ν L p ( Ω )  and bounded by below;

take a minimizer 𝒘 of ((${\mathcal{B}}$)) and consider the time-depending horizontal vector field

(5.22) 𝒘 ^ ( t , x ) = 𝒘 ( x ) ( 1 - t ) μ ( x ) + t ν ( x ) .

We know that

μ - ν = div H 𝒘 = div ϵ 𝒘

for any ϵ > 0 : hence, the linear interpolating curve μ ¯ ( t , x ) = ( 1 - t ) μ ( x ) + t ν ( x ) is a positive measure-valued distributional solution of the continuity equation (5.18) in 2 n + 1 , with initial datum λ ( 0 , ) = μ . Since 𝒘 is a solution to ((${\mathcal{B}}$)), then

0 1 Ω | 𝒘 ^ ( t , x ) | ϵ 𝑑 μ ¯ ( t , x ) 𝑑 t = Ω | 𝒘 ( x ) | H 𝑑 x < + ,

and we can apply [8, Theorem 5.8]: it follows that there exists a probability measure Q 𝒫 ( C ( [ 0 , 1 ] , Ω ¯ ) ) such that Q ( AC ( [ 0 , 1 ] , Ω ¯ ) ) = 1 ,

μ ¯ ( t , ) = ( e t ) # Q

and Q-a.e. curve σ is an integral curve of 𝒘 ^ ( t , x ) . Then Q-a.e. σ is horizontal curve because the vector field 𝒘 ( t , x ) is horizontal, and hence Q 𝒬 H ( μ , ν ) .

Following the same computation as (5.20), one can prove that

Ω ¯ φ ( x ) 𝑑 i Q ( x ) = Ω φ ( x ) | 𝒘 ( x ) | H 𝑑 x

for every φ C ( Ω ¯ ) . This implies that i Q = | 𝒘 | H and thus Q 𝒬 H p ( μ , ν ) and it solves ((${\mathcal{W}}$)), concluding the proof.

Proposition 5.13 and Theorem 5.14, combined with the existence of solutions for ((${\mathcal{B}}$)) as established by Theorem 5.8, lead to the following two straightforward consequences.

Corollary 5.15.

We have Q H p ( μ , ν ) if, and only if, μ - ν ( H W 1 , q ) ( Ω ) .

Corollary 5.16.

If Q H p ( μ , ν ) , or equivalently μ - ν ( H W 1 , q ) ( Ω ) , then:

  1. ( (${\mathcal{B}}$) ) = ( (${\mathcal{W}}$) ) and both ( (${\mathcal{B}}$) ) and ( (${\mathcal{W}}$) ) admit solutions,

  2. for every Q solving ( (${\mathcal{W}}$) ), 𝒘 Q solves ( (${\mathcal{B}}$) ),

  3. for every 𝒘 solving ( (${\mathcal{B}}$) ), there is Q 𝒘 that solves ( (${\mathcal{W}}$) ).

Remark 5.17.

If the function i G ( x , i ) is strictly increasing, for any x Ω ¯ , then for any Q 𝒬 H p ( μ , ν ) , Q solves ((${\mathcal{W}}$)) if and only if 𝒘 Q solves ((${\mathcal{B}}$)) and | 𝒘 Q | H = i Q .

Indeed, Corollary 5.16, Lemma 5.3 and the monotonicity of G imply

( 𝒲 ) = ( ) Ω G ( x , | 𝒘 Q ( x ) | H ) 𝑑 x Ω G ( x , i Q ( x ) ) 𝑑 x for all  Q 𝒬 H p ( μ , ν ) .

Therefore the result follows.

On the other hand, the convexity assumption on i G ( x , i ) , for any x Ω ¯ , ensures the existence of minimizers for both problems ((${\mathcal{B}}$)) and ((${\mathcal{W}}$)). However, this convexity plays no role in the proof of the equivalence ((${\mathcal{B}}$)) = ((${\mathcal{W}}$)).


Communicated by Zoltan Balogh


Funding statement: Michele Circelli is supported by the project PRIN 2022 F4F2LH – CUP J53D23003760006, Regularity problems in sub-Riemannian structures. Albert Clop acknowledges partial support from grant PID2020-112881GB-I00, PID2021-125021NAI00, PID2024-156429NB-I00 (Spanish Government) and 2021-SGR-00071 (Catalan Government). Albert Clop is also supported by the Spanish State Research Agency, through the Severo Ochoa and María de Maeztu Program for Centers and Units of Excellence in R&D (CEX2020-001084-M), as well as the CERCA Programme/Generalitat de Catalunya.

Acknowledgements

The authors would like to thank the anonymous referee for the careful reading of the draft and for the valuable suggestions, which significantly contributed to the improvement of this paper.

References

[1] A. Agrachev, D. Barilari and U. Boscain, A Comprehensive Introduction to Sub-Riemannian Geometry, Cambridge Stud. Adv. Math. 181, Cambridge University, Cambridge, 2019. 10.1017/9781108677325Search in Google Scholar

[2] A. Agrachev and P. Lee, Optimal transportation under nonholonomic constraints, Trans. Amer. Math. Soc. 361 (2009), no. 11, 6019–6047. 10.1090/S0002-9947-09-04813-2Search in Google Scholar

[3] L. Ambrosio, E. Brué and D. Semola, Lectures on Optimal Transport, Unitext 130, Springer, Cham, 2021. 10.1007/978-3-030-72162-6Search in Google Scholar

[4] L. Ambrosio, N. Fusco and D. Pallara, Functions of Bounded Variation and Free Discontinuity Problems, Oxford Math. Monogr., Oxford University, New York, 2000. 10.1093/oso/9780198502456.001.0001Search in Google Scholar

[5] L. Ambrosio and S. Rigot, Optimal mass transportation in the Heisenberg group, J. Funct. Anal. 208 (2004), no. 2, 261–301. 10.1016/S0022-1236(03)00019-3Search in Google Scholar

[6] L. Ambrosio and D. Trevisan, Lecture notes on the DiPerna–Lions theory in abstract measure spaces, Ann. Fac. Sci. Toulouse Math. (6) 26 (2017), no. 4, 729–766. 10.5802/afst.1551Search in Google Scholar

[7] M. Beckmann, A continuous model of transportation, Econometrica 20 (1952), 643–660. 10.2307/1907646Search in Google Scholar

[8] P. Bernard, Young measures, superposition and transport, Indiana Univ. Math. J. 57 (2008), no. 1, 247–275. 10.1512/iumj.2008.57.3163Search in Google Scholar

[9] S. Bianchini and F. Cavalletti, The Monge problem for distance cost in geodesic spaces, Comm. Math. Phys. 318 (2013), no. 3, 615–673. 10.1007/s00220-013-1663-8Search in Google Scholar

[10] V. I. Bogachev, Measure Theory. Vol. I, II, Springer, Berlin, 2007. 10.1007/978-3-540-34514-5Search in Google Scholar

[11] A. Bonfiglioli, E. Lanconelli and F. Uguzzoni, Stratified Lie Groups and Potential Theory for Their Sub-Laplacians, Springer Monogr. Math., Springer, Berlin, 2007. Search in Google Scholar

[12] L. Brasco and G. Carlier, Congested traffic equilibria and degenerate anisotropic PDEs, Dyn. Games Appl. 3 (2013), no. 4, 508–522. 10.1007/s13235-013-0081-zSearch in Google Scholar

[13] L. Brasco, G. Carlier and F. Santambrogio, Congested traffic dynamics, weak flows and very degenerate elliptic equations [corrected version of mr2584740], J. Math. Pures Appl. (9) 93 (2010), no. 6, 652–671. 10.1016/j.matpur.2010.03.010Search in Google Scholar

[14] L. Brasco and M. Petrache, A continuous model of transportation revisited, Zap. Nauchn. Sem. S.-Peterburg. Otdel. Mat. Inst. Steklov. (POMI) 411 (2013), 5–37. Search in Google Scholar

[15] D. Burago, Y. Burago and S. Ivanov, A Course in Metric Geometry, Grad. Stud. Math. 33, American Mathematical Society, Providence, 2001. 10.1090/gsm/033Search in Google Scholar

[16] L. Capogna, Regularity of quasi-linear equations in the Heisenberg group, Comm. Pure Appl. Math. 50 (1997), no. 9, 867–889. 10.1002/(SICI)1097-0312(199709)50:9<867::AID-CPA3>3.0.CO;2-3Search in Google Scholar

[17] L. Capogna and G. Citti, Regularity for subelliptic PDE through uniform estimates in multi-scale geometries, Bull. Math. Sci. 6 (2016), no. 2, 173–230. 10.1007/s13373-015-0076-8Search in Google Scholar

[18] G. Carlier, C. Jimenez and F. Santambrogio, Optimal transportation with traffic congestion and Wardrop equilibria, SIAM J. Control Optim. 47 (2008), no. 3, 1330–1350. 10.1137/060672832Search in Google Scholar

[19] F. Cavalletti, An overview of L 1 optimal transportation on metric measure spaces, Measure Theory in Non-Smooth Spaces, Partial Differ. Equ. Meas. Theory, De Gruyter, Warsaw (2017), 98–144. 10.1515/9783110550832-003Search in Google Scholar

[20] F. Cavalletti and A. Mondino, Optimal maps in essentially non-branching spaces, Commun. Contemp. Math. 19 (2017), no. 6, Article ID 1750007. 10.1142/S0219199717500079Search in Google Scholar

[21] T. Champion and L. De Pascale, The Monge problem for strictly convex norms in d , J. Eur. Math. Soc. (JEMS) 12 (2010), no. 6, 1355–1369. 10.4171/jems/234Search in Google Scholar

[22] M. Circelli, Gradient estimates for an orthotropic nonlinear diffusion equation in the Heisenberg group, preprint (2025), https://arxiv.org/abs/2505.01354. Search in Google Scholar

[23] M. Circelli and G. Citti, Transport densities and congested optimal transport problem in the Heisenberg group, preprint (2024), https://arxiv.org/abs/2303.16052. Search in Google Scholar

[24] M. Circelli, G. Citti and A. Clop, Lipschitz regularity for solutions to an orthotropic q-laplacian-type equation in the Heisenberg group, preprint (2025), https://arxiv.org/abs/2407.07548. Search in Google Scholar

[25] G. Citti and A. Sarti, A cortical based model of perceptual completion in the roto-translation space, J. Math. Imaging Vision 24 (2006), no. 3, 307–326. 10.1007/s10851-005-3630-2Search in Google Scholar

[26] G. E. Comi and V. Magnani, The Gauss–Green theorem in stratified groups, Adv. Math. 360 (2020), Article ID 106916. 10.1016/j.aim.2019.106916Search in Google Scholar

[27] L. De Pascale and S. Rigot, Monge’s transport problem in the Heisenberg group, Adv. Calc. Var. 4 (2011), no. 2, 195–227. 10.1515/acv.2010.026Search in Google Scholar

[28] I. Ekeland, Convexity Methods in Hamiltonian Mechanics, Ergeb. Math. Grenzgeb. (3) 19, Springer, Berlin, 1990. 10.1007/978-3-642-74331-3Search in Google Scholar

[29] L. C. Evans, Partial differential equations and Monge–Kantorovich mass transfer, Current Developments in Mathematics (Cambridge 1997), International Press, Boston (1999), 65–126. 10.4310/CDM.1997.v1997.n1.a2Search in Google Scholar

[30] M. Feldman and R. J. McCann, Monge’s transport problem on a Riemannian manifold, Trans. Amer. Math. Soc. 354 (2002), no. 4, 1667–1697. 10.1090/S0002-9947-01-02930-0Search in Google Scholar

[31] A. Figalli and L. Rifford, Mass transportation on sub-Riemannian manifolds, Geom. Funct. Anal. 20 (2010), no. 1, 124–159. 10.1007/s00039-010-0053-zSearch in Google Scholar

[32] G. B. Folland and E. M. Stein, Hardy Spaces on Homogeneous Groups, Math. Notes 28, Princeton University, Princeton, 1982. 10.1515/9780691222455Search in Google Scholar

[33] B. Franchi, R. Serapioni and F. Serra Cassano, Rectifiability and perimeter in the Heisenberg group, Math. Ann. 321 (2001), no. 3, 479–531. 10.1007/s002080100228Search in Google Scholar

[34] N. Garofalo and D.-M. Nhieu, Isoperimetric and Sobolev inequalities for Carnot–Carathéodory spaces and the existence of minimal surfaces, Comm. Pure Appl. Math. 49 (1996), no. 10, 1081–1144. 10.1002/(SICI)1097-0312(199610)49:10<1081::AID-CPA3>3.0.CO;2-ASearch in Google Scholar

[35] Z. Ge, Collapsing Riemannian metrics to Carnot–Carathéodory metrics and Laplacians to sub-Laplacians, Canad. J. Math. 45 (1993), no. 3, 537–553. 10.4153/CJM-1993-028-6Search in Google Scholar

[36] M. Gromov, Structures métriques pour les variétés riemanniennes, Textes Math. 1, CEDIC, Paris, 1981. Search in Google Scholar

[37] G. Monge, Mémoire sur la théorie des déblais et des remblais, Mem. Math. Phys. Acad. Roy. Sci. (1781), 666–704. Search in Google Scholar

[38] R. Monti, Distances, boundaries and surface measures in Carnot–Carathéodory spaces, UTM PhDTS, Department of Mathematics, University of Trento, 2001. Search in Google Scholar

[39] R. Monti and M. Rickly, Geodetically convex sets in the Heisenberg group, J. Convex Anal. 12 (2005), no. 1, 187–196. Search in Google Scholar

[40] S. Mukherjee and Y. Sire, Regularity of inhomogeneous quasi-linear equations on the Heisenberg group, Anal. Theory Appl. 37 (2021), no. 4, 520–540. 10.4208/ata.OA-2020-0040Search in Google Scholar

[41] P. Pansu, Métriques de Carnot–Carathéodory et quasiisométries des espaces symétriques de rang un, Ann. of Math. (2) 129 (1989), no. 1, 1–60. 10.2307/1971484Search in Google Scholar

[42] J. Petitot and Y. Tondut, Vers une neurogéométrie. Fibrations corticales, structures de contact et contours subjectifs modaux, Math. Inform. Sci. Humaines (1999), no. 145, 5–101. 10.4000/msh.2809Search in Google Scholar

[43] S. T. Rachev and L. Rüschendorf, Mass Transportation Problems. Vol. II, Probab. Appl. (New York), Springer, New York, 1998. Search in Google Scholar

[44] F. Santambrogio, A Dacorogna–Moser approach to flow decomposition and minimal flow problems, Congrès SMAI 2013, ESAIM Proc. Surveys 45, EDP Sciences, Les Ulis (2014), 265–274. 10.1051/proc/201445027Search in Google Scholar

[45] F. Santambrogio, Optimal Transport for Applied Mathematicians, Progr. Nonlinear Differential Equations Appl. 87, Birkhäuser/Springer, Cham, 2015. 10.1007/978-3-319-20828-2Search in Google Scholar

[46] F. Santambrogio, Regularity via duality in calculus of variations and degenerate elliptic PDEs, J. Math. Anal. Appl. 457 (2018), no. 2, 1649–1674. 10.1016/j.jmaa.2017.01.030Search in Google Scholar

[47] F. Santambrogio, A Course in the Calculus of Variations—Optimization, Regularity, and Modeling, Universitext, Springer, Cham, 2023. 10.1007/978-3-031-45036-5Search in Google Scholar

[48] X. Zhong, Regularity for variational problems in the Heisenberg group, preprint (2017), https://arxiv.org/abs/1711.03284. Search in Google Scholar

Received: 2024-11-19
Accepted: 2025-06-12
Published Online: 2025-08-29
Published in Print: 2025-10-01

© 2025 Walter de Gruyter GmbH, Berlin/Boston

This work is licensed under the Creative Commons Attribution 4.0 International License.

Downloaded on 31.1.2026 from https://www.degruyterbrill.com/document/doi/10.1515/acv-2024-0119/html
Scroll to top button