Skip to content
BY 4.0 license Open Access Published by De Gruyter March 12, 2021

Sharp Hardy Identities and Inequalities on Carnot Groups

  • Joshua Flynn , Nguyen Lam and Guozhen Lu EMAIL logo

Abstract

In this paper we establish general weighted Hardy identities for several subelliptic settings including Hardy identities on the Heisenberg group, Carnot groups with respect to a homogeneous gauge and Carnot–Carathéodory metric, general nilpotent groups, and certain families of Hörmander vector fields. We also introduce new weighted uncertainty principles in these settings. This is done by continuing the program initiated by [N. Lam, G. Lu and L. Zhang, Factorizations and Hardy’s-type identities and inequalities on upper half spaces, Calc. Var. Partial Differential Equations 58 2019, 6, Paper No. 183; N. Lam, G. Lu and L. Zhang, Geometric Hardy’s inequalities with general distance functions, J. Funct. Anal. 279 2020, 8, Article ID 108673] of using the Bessel pairs introduced by [N. Ghoussoub and A. Moradifam, Functional Inequalities: New Perspectives and New Applications, Math. Surveys Monogr. 187, American Mathematical Society, Providence, 2013] to obtain Hardy identities. Using these identities, we are able to improve significantly existing Hardy inequalities in the literature in the aforementioned subelliptic settings. In particular, we establish the Hardy identities and inequalities in the spirit of [H. Brezis and J. L. Vázquez, Blow-up solutions of some nonlinear elliptic problems, Rev. Mat. Univ. Complut. Madrid 10 1997, 443–469] and [H. Brezis and M. Marcus, Hardy’s inequalities revisited. Dedicated to Ennio De Giorgi, Ann. Sc. Norm. Super. Pisa Cl. Sci. (4) 25 1997, 1–2, 217–237] in these settings.

1 Introduction

Let G denote a Carnot group with sub-Laplacian =j=1mXj2, horizontal gradient H=(X1,,Xm), homogeneous dimension Q, and -gauge d (see preliminaries for precise definitions). In this paper, we employ factorization of differential operators and Bessel pairs to obtain weighted Hardy identities on Carnot groups G of the form

(1.1) B R ( 0 ) V ( d ) | H ( u φ ( d ) ) | 2 φ 2 ( d ) 𝑑 x = B R ( 0 ) V ( d ) | H u | 2 𝑑 x - B R ( 0 ) W ( d ) u 2 | H d | 2 𝑑 x

for uC0(BR(0){0}), where (V,W) is a pair of positive C1-functions on (0,R) which satisfy the condition that there exists a positive C1-solution φ:(0,R) to the differential equation

( V φ ) + W φ = 0 on  ( 0 , R ) .

Such a pair of functions is called a Bessel pair. Here, dx is a Haar measure on G, and BR(0) is the gauge ball determined by d with radius R and center being the origin 0. This result is contained in Theorem 2 below. As a particular important application, we obtain weighted Hardy identities on the Heisenberg group HN=2N+1, the simplest and most important nonabelian Carnot group. It is important to note that, for an appropriately chosen Bessel pair, our results recover the Hardy identity corresponding to known Hardy inequalities on the Heisenberg group and more general Carnot groups, thereby strongly improving these known inequalities by finding exact remainder terms.

In addition to these identities, we also obtain Hardy identities of the form

(1.2)

0 < d K < R V ( d K ) | H ( u φ ( d K ) ) | 2 φ ( d K ) 2 𝑑 x = 0 < d K < R V ( d K ) | H u | 2 𝑑 x - 0 < d K < R W ( d K ) u 2 𝑑 x
- 0 < d K < R V ( d K ) | u | 2 φ - 1 ( d K ) φ ( d K ) [ d K - α - 1 d K ] 𝑑 x ,

where H is the horizontal gradient with respect to certain smooth vector fields X1,,Xm, and dK is the Carnot–Carathéodory distance (also known CC distance) to a closed subset K with respect to this system of vector fields. A key point here is that, by [24], the CC distance function dK satisfies the horizontal eikonal equation |HdK|=1, unlike a general homogeneous gauge d for which |Hd|=1 is typically false. This is the content of Theorem 16.

We recall that, on the Heisenberg group HN=2N+1, the sharp Hardy inequality

(1.3) H N | H u | 2 𝑑 z 𝑑 t ( Q - 2 ) 2 4 H N | u | 2 d 2 | H d | 2 𝑑 z 𝑑 t , u C 0 ( H N { 0 } ) ,

was established by Garofalo and Lanconelli in [14]. The Lp-Hardy inequality on the Heisenberg group for 1<p<Q=2N+2 was established by Niu, Zhang and Wang [26] by using the Picone’s identity for the p-subLaplacian (see also Zhang and Niu [28])

H N | H u | p 𝑑 z 𝑑 t ( Q - p p ) p H N | u | p d p | H d | p 𝑑 z 𝑑 t , u C 0 ( H N { 0 } ) .

When d is replaced by the Carnot–Carathéodory metric dcc, then there holds

(1.4) H N | H u | 2 𝑑 z 𝑑 t c H N | u | 2 d cc 2 𝑑 z 𝑑 t , u C 0 ( H N { 0 } ) ,

for some constant c>0. See for example [13, 22, 27]. Using the fact that dcc and d are equivalent metrics, one may also obtain a Hardy inequality of the form (1.4) with dcc replaced by d, though the best constant appears to be unknown in either case. Nevertheless, while (1.3) is sharp, there are no extremizers, and so it is natural to try to improve this inequality by adding nontrivial nonnegative terms to the right-hand side of (1.3). As mentioned above, we in fact determine exact remainder terms so that (1.3) becomes

H N | H u | 2 𝑑 z 𝑑 t = ( Q - 2 ) 2 4 H N | u | 2 d 2 | H d | 2 𝑑 z 𝑑 t + remainder rerms .

Therefore, we obtain a strict improvement of the known Hardy inequality (1.3), and the counterparts for Carnot groups. This is the content of Theorem 6.

The approach used here is inspired by the recent work of Lam, Lu and Zhang in [20] where they used factorization of differential operators and Bessel pairs to prove, among other things, various Hardy identities on the half-space

+ N = { ( x 1 , , x N ) N : x N > 0 } .

Their main identity theorem is stated as follows.

Theorem A.

Let 0<R, and let V and W be positive C1-functions on (0,R). If (rN+1V,rN+1W) is a Bessel pair on (0,R) with positive solution φ, then there holds

B R ( 0 ) + N V ( | x | ) | u | 2 𝑑 x - B R ( 0 ) + N [ W ( | x | ) - V ( | x | ) | x | ] | u | 2 𝑑 x = B R ( 0 ) + N V ( | x | ) | ( u φ 1 x N ) | φ 2 x N 2 𝑑 x

for uC0(BR(0)R+N),

In addition to this theorem, Lam, Lu and Zhang also proved the following Hardy identities on strips of the form A2R=N-1×(0,2R) which imply in particular Hardy identities on the half-space +N with the Euclidean norm |x| replaced by the distance function d(x)=xN, the height from the boundary. Their result is stated as follows.

Theorem B.

Let 0<R, let V and W be positive C1-functions on (0,R), and let dR=min{xN,2R-xN} be a distance function to the boundary of A2R. If (V,W) is a Bessel pair on (0,R) with positive solution φ, then for uC0(A2R) there holds

A 2 R V ( d R ( x ) ) | u | 2 𝑑 x - A 2 R W ( d R ( x ) ) | u | 2 𝑑 x = A 2 R V ( d R ( x ) ) | ( u φ ( d R ( x ) ) ) | 2 φ 2 ( d R ( x ) ) 𝑑 x

and

A 2 R V ( d R ( x ) ) | u x N | 2 𝑑 x - A 2 R W ( d R ( x ) ) | u | 2 𝑑 x = A 2 R V ( d R ( x ) ) | x N ( u φ ( d R ( x ) ) ) | 2 φ 2 ( d R ( x ) ) 𝑑 x .

Actually, similar half-space Hardy inequalities have been established for the Heisenberg group. We recall the following result of Luan and Yang given in [23].

Theorem C.

Let HN+={(z,t)HN:t>0} denote the upper half-space of HN. Then there holds for uC0(HN+),

(1.5) H N + | H u | 2 𝑑 z 𝑑 t H N + | z | 2 t 2 | u | 2 𝑑 z 𝑑 t

and this inequality is sharp and without extremizers.

See [18] for similar Hardy inequalities on different Heisenberg half spaces.

In the spirit of Theorems A, B and 1.5, we establish the Hardy identities on strips

A 2 R V ( d R ) | H ( u φ ( d R ) ) | 2 φ 2 ( d R ) 𝑑 z 𝑑 t = A 2 R V ( d R ) | u | 2 𝑑 z 𝑑 t - 1 4 A 2 R W ( d R ) | z | 2 | u | 2 𝑑 z 𝑑 t ,
A 2 R V φ 2 ( d R ) | t ( u φ ( d R ) ) | 2 𝑑 z 𝑑 t = A 2 R V ( d R ) | u t | 2 𝑑 z 𝑑 t - A 2 R W ( d R ) | u | 2 𝑑 z 𝑑 t ,

and half-gauge balls

B R ( 0 ) H N + V | ( u φ 1 t ) | 2 φ 2 t 2 𝑑 z 𝑑 t = B R ( 0 ) H N + V | u | 2 𝑑 z 𝑑 t - B R ( 0 ) H N + [ W - V d ] | u | 2 | d | 2 𝑑 z 𝑑 t .

Here dR(z,t)=min{t,2R-t} and d(z,t)=t is the vertical height function. The identity on strips implies in particular improvements of (1.5). These identities are the contents of Theorems 8 and 12.

We note that, in [23], (1.5) has an interesting improved version by adding to the right side a nonnegative nontrivial term with weight given in terms of the homogeneous gauge d; however, they do not obtain exact remainder terms. Moreover, this inequality is not coordinate independent, and so the best constant appearing in front of |x|2t-2|u|2𝑑x depends on the coordinates chosen for the left invariant vector fields on HN. In particular, our choice of coordinates are different and the best constant turns out to be 116 for our choice of coordinates.

We also mention that, our Hardy identities of the form (1.2) for the general distance function dK are in the spirit of the following Hardy identities for general distance functions (for example, the distance to a subset) on Euclidean space obtained by Lam, Lu and Zhang in [21].

Theorem D.

Let 0<R, let V and W be positive C1-functions on (0,R), and let d be a distance function. Assume that for some αR, Δd(x)-α-1d(x) exists on {0<d(x)<R} in the sense of distributions and (rα-1V,rα-1W) is a Bessel pair on (0,R) with positive solution φ. Then, for uC0({0<d(x)<R}),

0 < d ( x ) < R V ( d ( x ) ) | u ( x ) | 2 𝑑 x - 0 < d ( x ) < R W ( d ( x ) ) | u ( x ) | 2 𝑑 x
= 0 < d ( x ) < R V ( d ( x ) ) φ 2 ( d ( x ) ) | ( u ( x ) φ ( d ( x ) ) ) | 2 𝑑 x - 0 < d ( x ) < R V ( d ( x ) ) | u ( x ) | 2 [ Δ d ( x ) - α - 1 d ( x ) ] φ ( d ( x ) ) φ ( d ( x ) ) 𝑑 x .

Among the particularly interesting applications of the Hardy identity (1.1) are Hardy identities and inequalities which improve Hardy inequalities of the form

B R ( 0 ) | u | 2 𝑑 x B R ( 0 ) | u | 2 d 2 | d | 2 𝑑 x

by adding positive nontrivial L2-terms to the right-hand side. In fact, Brezis and Vázquez established the remarkable inequality

Ω | u | 2 𝑑 x - ( N - 2 2 ) 2 Ω | u | 2 | x | 2 𝑑 x z 0 2 ω N 2 N | Ω | 2 N Ω | u | 2 𝑑 x ,

where ΩN, N>2, is any bounded domain, uH01(Ω), and z0 is the first zero of the Bessel function J0. (See also [4] for the distance function the boundary.) We introduce similar inequalities on Carnot groups (see Theorem 6) by establishing identities of the form

B R ( 0 ) | H u | 2 𝑑 x - ( N - 2 2 ) 2 B R ( 0 ) | u | 2 d 2 𝑑 x = z 0 2 R 2 B R ( 0 ) | u | 2 𝑑 x + B R ( 0 ) | H ( u d Q 2 J 0 ; R ( d ) ) | 2 J 0 ; R ( d ) 2 d Q 𝑑 x ,

where J0;R(r)=J0(z0rR). This is done by choosing an appropriate Bessel pair (V,W) and applying (1.1). We note that such Brezis–Vázquez-type inequalities [5] on Carnot groups were observed first by Kombe in [19].

We take the opportunity to mention some forthcoming works which continue the program of establishing Hardy identities by use of Bessel pairs, including the work in this paper and the work of Lam, Lu and Zhang already mentioned. First, we mention that Duy, Lam and Lu establish recently in [7] the following Lp-Hardy identities by use of p-Bessel pairs.

Theorem E.

Let N1, p>1, 0<R, and let A and W be positive C1-functions on (0,R). If (rN-1V,rN-1W) is a p-Bessel pair on (0,R), then for all uC0(BR(0){0}),

B R ( 0 ) V ( | x | ) | u | p 𝑑 x - B R ( 0 ) W ( | x | ) | u | p 𝑑 x = B R ( 0 ) V ( | x | ) C p ( u , φ ( u φ ) ) 𝑑 x

and

B R ( 0 ) V ( | x | ) | u | p 𝑑 x - B R ( 0 ) W ( | x | ) | u | p 𝑑 x = B R ( 0 ) V ( | x | ) C p ( u , φ ( u φ ) ) 𝑑 x .

Here φ is the positive solution of

( r N - 1 V ( r ) | y | p - 2 y ) + r N - 1 W ( r ) | y | p - 2 y = 0

on (0,R). Moreover, R is the radial derivative, and Cp is defined by

C p ( x , y ) = | x | p - | x - y | p - p | x - y | p - 2 ( x - y ) y .

We also point out here that we have established in [10] the Lp-Hardy identities for general distance functions in the Euclidean setting in which we generalize Theorem D to the Lp-setting. We have also established in [12] the Lp-Hardy identities for domains in Euclidean space. Our Lp-Hardy identities from [10] for general distance functions is stated as follows.

Theorem F.

Let p>1, 0<R, let d be a distance function on RN, and let V,WC1(0,R) be positive. Assume that for some αR, Δd(x)-α-1d(x) exists on {0<d(x)<R} in the sense of distributions, and that (rα-1V,rα-1W) is a p-Bessel pair on (0,R) with positive solution φ. Then, for uC0({0<d(x)<}), there holds

0 < d ( x ) < R V ( d ( x ) ) | u | p 𝑑 x - 0 < d ( x ) < R W ( d ( x ) ) | u | p 𝑑 x
= 0 < d ( x ) < R V ( d ( x ) ) C p ( u , φ ( d ( x ) ) ( u φ ( d ( x ) ) ) ) 𝑑 x
- 0 < d ( x ) < R V ( d ( x ) ) | u | p φ ( d ( x ) ) 1 - p | φ ( d ( x ) ) | p - 2 φ ( d ( x ) ) [ Δ d ( x ) - α - 1 d ( x ) ] 𝑑 x

and

0 < d ( x ) < R V ( d ( x ) ) | d u | p 𝑑 x - 0 < d ( x ) < R W ( d ( x ) ) | u | p 𝑑 x
= 0 < d ( x ) < R V ( d ( x ) ) C p ( d u , φ ( d ( x ) ) d ( u φ ( d ( x ) ) ) ) 𝑑 x
- 0 < d ( x ) < R V ( d ( x ) ) | u | p φ ( d ( x ) ) 1 - p | φ ( d ( x ) ) | p - 2 φ ( d ( x ) ) [ Δ d ( x ) - α - 1 d ( x ) ] 𝑑 x .

Another important geometric inequality is the Hardy uncertainty principle given by

(1.6) | Q + λ - 2 | 2 G u 2 | H d | 2 𝑑 x ( G u 2 d 2 | H d | 2 𝑑 x ) 1 2 ( G | H u | 2 𝑑 x ) 1 2 .

The uncertainty principle for the Heisenberg group may be found in [14]. Using factorization, we also obtain a new family of inequalities which generalize (1.6) to a family of uncertainty principles whose weights are given in terms of Bessel pairs. These inequalities are of the form

1 2 B R ( 0 ) u 2 | H d | 2 [ W + V ( φ φ ) 2 ] 𝑑 x ( B R ( 0 ) u 2 | H d | 2 V ( φ φ ) 2 𝑑 x ) 1 2 ( B R ( 0 ) V | H u | 2 𝑑 x ) 1 2

and are the content of Theorem 13. These inequalities appear to be entirely new even for Euclidean space, and they give a new interpretation of the classical Hardy uncertainty principle.

We recall that Ghoussoub and Moradifam introduced in [15] (see also their book [16]) the notion of Bessel pairs to enhance and unify several improved Hardy-type inequalities in the literature for the Euclidean setting. They established the following important characterization theorem.

Theorem G.

Assume N1, let 0<R, and let V and W be C1-functions on (0,R) such that

0 R 1 r N - 1 V ( r ) 𝑑 r = 𝑎𝑛𝑑 0 R r N - 1 V ( r ) 𝑑 r < .

If (rN-1V,rN-1W) is a Bessel pair on (0,R), then for all uC0(BR(0)) there holds

(1.7) B R ( 0 ) V ( | x | ) | u | 2 B R ( 0 ) W ( | x | ) | u | 2 𝑑 x .

Also, if (1.7) holds for all uC0(BR(0)), then (rN-1V,rN-1cW) is a Bessel pair on (0,R) for some c>0.

It is worth mentioning that these results extend straightforwardly to more general settings. We describe two such generalizations now. The first is the case of simply connected nilpotent Lie groups G where we specify a system of left invariant vector fields X={X1,,Xm} which bracket generate the Lie algebra of G. Then, as was shown by Nagel, Ricci and Stein in [25], the sub-Laplacian =j=1mXj2 has a global fundamental solution ΓX(x,y). Then (1.1) holds if we take H=(X1,,Xm), d(x)=dX(x)=ΓX(x,0)2-q for some q>2, and BR(0)={dX(x)<R}.

The second is the case that X={X1,,Xm} is a system of Hörmander vector fields homogeneous of degree -1 (see the preliminaries for precise details). In fact, recently it was shown in [1] that, for these kinds of vector fields, the Hörmander operator =j=1mXj2 admits a global fundamental solution ΓX(x,y), and one can, similar to the nilpotent case, take d(x)=ΓX(x,0)2-q, where q is the homogeneous dimension arising from the family of dilations used to define X. Analogous Hardy identities may be formulated just as in the nilpotent case just mentioned. These identities are the contents of Theorem 20.

Finally, we mention that the Lp-weighted Hardy’s identities and inequalities (for p2) with the p-Bessel pairs on Carnot groups and for Hörmander vector fields have been established in [11]. This is inspired by the earlier works of the authors on Lp-weighted Hardy’s identities and inequalities (for p2) with the p-Bessel pairs in Euclidean spaces [7, 10, 12].

Organization.

The paper is outlined as follows. We first recall in the Preliminaries section facts about Hörmander vector fields, the Heisenberg group, and general Carnot groups. In the following section, we first state and prove the main factorization lemma. From here, we establish the main Hardy identities on Carnot groups, and then apply this result by choosing specific Bessel pairs. Next, we establish another factorization lemma and use it to establish Hardy identities on the Heisenberg half-space and half-gauge ball. Next, we use the factorization lemma to state and prove general weighted uncertainty principles as well as some applications. Then we establish Hardy identities for general distance functions dK, where dK is the CC distance to a closed set K. Lastly, we state and prove Hardy identities and specific inequalities for more general settings including certain vector fields and simply connected nilpotent Lie groups.

For convenience of the reader, we point out explicitly the main results of the paper. Firstly, in Theorem 2 we prove the main Hardy identity for Carnot groups. Secondly, in Theorem 6, we apply Theorem 2 to obtain a Hardy identity in the spirit of Brezis–Vázquez. Thirdly, in Theorem 8 and Theorem 12, we prove Hardy identities related to the Heisenberg half space. Fourthly, in Theorem 13, we introduce a new class of uncertainty principles for groups. Lastly, in Theorems 16 and 20, we prove Hardy identities for the more general settings with respect to the CC metric and global fundamental solution, respectively.

Remark 1.

We employ the following notational conventions. Firstly, if a theorem is the result of a referenced paper, then that theorem is labeled by a Latin letter (e.g., such as Theorem A). Secondly, for a general Carnot group, its measure will generally be denoted by dx; however, in the case of the Heisenberg group, we opt to use dzdt when it is more natural. Secondly, general points shall be denoted by x, and the Lebesgue measure shall be denoted by dx. In particular, a point in the Heisenberg group may be denoted by (x1,,xN,y1,,yN,t) or (z,t) or x, and its measure will always be denoted by dx or dzdt. Lastly, in general, d shall denote a homogeneous gauge, dX shall denote a power of the fundamental solution corresponding to a system of vector fields X, dcc shall denote the Carnot–Carathéodory (CC) metric with respect to a given system of vector fields, and dK shall denote the CC distance to a closed subset K.

2 Preliminaries

We begin by recalling facts about general vector fields, and then proceed to recalling facts and properties about the Heisenberg group and more general Carnot groups. Thus, let X={X1,,Xm} denote a collection of smooth vector fields on N, i.e.,

X j = i = 1 N a i j ( x ) x i , j = 1 , , m ,

with aijC(N). A subunit curve for X is a Lipschitz continuous curve γ:[0,T]N, T>0, such that there exist measurable functions h1,,hm:[0,T] satisfying

γ ( t ) = j = 1 m h j ( t ) X j ( γ ( t ) ) , j = 1 m h j 2 ( t ) 1

for almost every t[0,T]. The Carnot–Carathéodory distance between two points x,yN is then defined to be

d cc ( x , y ) = inf { T 0 : there exists a subunit curve  γ : [ 0 , T ] N  with  γ ( 0 ) = x , γ ( T ) = y } .

In case the vector fields X1,,Xm are such that the Lie algebra generated by X is N-dimensional at every point in N, then X is said to satisfy Hörmander’s condition. The Chow–Rashevsky theorem asserts that, if X satisfies Hörmander’s condition, then dcc(x,y)< for all x,yN. In particular, dcc becomes a metric on N.

Associated with X is the Hörmander sum of squares operator X=j=1mXj2, and, in case X satisfies the Hörmander condition, it is a well-known result of Hörmander that X is hypoelliptic. We also define the horizontal gradient H and divergence as follows:

H u = ( X 1 u , , X m u ) , div H ( a 1 , , a m ) = j = 1 m X j * a j ,

where a=(a1,,am):Nm is any sufficiently regular function, and Xj* is the formal adjoint of Xj.

Now, suppose for convenience that X has a global fundamental solution ΓX(x,y). Then motivated by the fact that, for N3, the function Γeu(x,y)=1N(N-2)ωN|x-y|2-N, where ωN is the volume of the Euclidean unit ball, is the fundamental solution for the Euclidean Laplacian, and |x|=Γeu(x,0)12-Q defines the Euclidean norm, we may define the function dX(x)=ΓX(x,0)12-N and use it to define the corresponding ball

B X , R ( 0 ) = { x N : d X ( x ) < R } .

An important family of Hörmander vector fields are the left-invariant vector fields of the Heisenberg group HN=2N+1. Recall that H may be defined as 2N+1 equipped with the group law

( a , b , t ) ( a , b , t ) = ( a + a , b + b , t + t + 1 2 ( a , b - a , b ) ) ,

where a,b,a,bN are horizontal coordinates, and t,t are central coordinates. At a given point (z,t)=(x1,,xN,y1,,yN,t)HN, the corresponding left invariant vector fields are given by

(2.1) X i = x i - y i 2 t , Y i = y i + x i 2 t , T = t

for i=1,,N. These vector fields satisfy the Heisenberg relation

[ X i , Y j ] = δ i j T , i , j = 1 , , N ,

and so the Lie algebra 𝔥 of H is a graded nilpotent Lie algebra of step 2. The corresponding sub-Laplacian is given by

(2.2) = j = 1 m ( X j 2 + Y j 2 ) = Δ 2 N + 1 4 | z | 2 2 t 2 - j = 1 m ( x j y j - y j x j ) t ,

and the horizontal gradient is given as above.

The Heisenberg group is a homogeneous group in the sense that it has a family of anisotropic dilations δλ:HH given by δλ(z,t)=(λz,λ2t) for λ>0. The corresponding homogeneous dimension of H is given by Q=2N+2. Moreover, H has a homogeneous gauge given by

(2.3) d ( z , t ) = ( | z | 4 + 16 t 2 ) 1 4 ,

and this norm is such that d2-Q is a scalar multiple of the fundamental solution of . In particular, for an appropriate c>0, the balls BX,cR(0) are exactly the balls gauge balls BR(0)={xHN:d(x)<R}. The Haar measure, which will be denoted by dx or dzdt, on H is given by the Lebesgue measure dx=dzdt on 2N+1. The Haar measure and homogeneous gauge enjoy the following homogeneity properties:

d ( δ λ ( z , t ) ) = λ Q d z d t , d ( δ λ ( z , t ) ) = λ d ( z , t ) .

The Heisenberg group belongs to a large family of nilpotent groups, called the Carnot groups, which are important in many fields of analysis and have a rich geometric structure. We recall Carnot groups now. First, recall that a stratified Lie algebra 𝔤 of step r is a Lie algebra with subspaces V1,,Vr satisfying

𝔤 = V 1 V r and [ V 1 , V i ] = V i + 1 , i = 1 , , r - 1 , [ V 1 , V r ] = 0 .

If G is a simply connected Lie group with a stratified Lie algebra 𝔤, then G is called a Carnot group. In fact, if Nj=dimVj, and using the exponential map, we may identify G with N1××Nr so that each point gG is identified with a point x=(x(1),,x(r)) in N1××Nr so that x(j)Nj; such an identification will always be made. Just as for the Heisenberg group, the Carnot groups have a natural family of anisotropic dilations δλ:GG given by

δ λ ( x ) = ( λ x ( 1 ) , λ 2 x ( 2 ) , , λ r x ( r ) ) .

The homogeneous dimension Q of the group is defined to be Q=j=1mjdimVj. Let dx denote a Haar measure on G, which can be taken to be the Lebesgue measure on N1×Nr. The left-invariant vector fields of the first layer on G may be written in the following way:

X j = x j + h = 2 r k = 1 N h a j , k ( h ) ( x ( 1 ) , , x ( h - 1 ) ) x k ( h ) , j = 1 , , N 1 ,

where aj,k are degree h-1 polynomials which are homogeneous with respect to δλ. These vector fields generate the Lie algebra of G and so every left invariant vector on G may be obtained from taking linear combinations of brackets of the Xj. The horizontal gradient, horizontal divergence, and sub-Laplacian of G are defined as before.

We now recall the notions of symmetric homogeneous norms and gauges in the following definitions.

Definition 1.

A nonnegative continuous function d on G is called a symmetric homogeneous norm provided

  1. d ( δ λ ( x ) ) = λ d ( x ) for all λ>0, xG,

  2. d ( x ) > 0 iff x0,

  3. d ( x - 1 ) = d ( x ) for every xG.

Definition 2.

Suppose d is a symmetric homogeneous norm on G which is smooth away from the origin. If d satisfies d2-Q=0 away from the origin, then d is called an -gauge. An -radial function is any function u:G{0} satisfying u=fd for a suitable f:(0,) and a given -gauge d.

Note that the Heisenberg homogeneous gauge (2.3) satisfies this definition. Moreover, it is well known that (e.g., see [3]), if Γ(x,y) denotes the fundamental solution of , then any -gauge d is a scalar multiple of Γ12-Q(x,0) away from the origin (see for example [3]). This fact allows for a nice formula for acting on radial functions, and we recall this useful formula in the following lemma.

Lemma A.

If f(d) is a smooth radial function on G{0}, then

( f ( d ) ) = | H d | 2 ( f ′′ ( d ) + Q - 1 d f ( d ) ) .

The proof may be found in [3, Proposition 5.4.3].

In fact, these properties of G and its fundamental solution hold in a more general setting. Indeed, if {Y1,,Ym} is a system of C(n) homogeneous linearly independent vector fields with homogeneous degree 1 with respect to an anisotropic dilation

δ λ : ( x 1 , , x n ) ( λ σ 1 x 1 , , λ σ n x n ) ,

where 1=σ1σn, and if this system satisfies Hörmander’s finite rank condition at x=0, then, by [1], their sum of squares operator Y=Y12++Ym2 admits a global fundamental solution ΓY satisfying

  1. Y , y Γ Y ( x , y ) = δ x ,

  2. Γ Y is strictly positive,

  3. Γ ( δ λ ( x ) ; δ λ ( y ) ) = λ 2 - q Γ ( x ; y ) for all xy and λ>0,

  4. Γ ( x ; y ) = Γ ( y ; x ) for xy.

Here, q=j=1nσj is the homogeneous dimension of the system, and by being homogeneous of degree 1 it is understood that

Y j ( f δ λ ) = λ ( Y j f ) δ λ

for all λ>0, fC(N), and j=1,,m. Following the proof of [3, Proposition 5.4.3], it is easy to see that dY=Γ12-q satisfies

Y d Y = ( q - 1 ) | Y d Y | 2 d Y ,

where Y=(Y1,,Ym). Moreover, if f(dY) is a radial function with respect to dY, then there holds

Y ( f ( d Y ) ) = | Y d | 2 ( f ′′ ( d Y ) + q - 1 d Y f ( d Y ) ) .

3 Main Results and Proofs

3.1 Hardy Identities on Carnot Groups

We begin by introducing our main factorization lemma which our main results rely on. While this lemma is stated for , H and divH on G, it is worth mentioning that this lemma holds for more general operators M, M*, the formal adjoint M*, and M*M replacing H, divH, and , respectively.

Lemma 1.

Let ΩG be a domain, let V,φ:ΩR>0 be C1(Ω), and let sR. Define

T s = V H - s V φ - 1 H φ , T s + = - div H ( V ) - s V φ - 1 H φ ,

where Ts+ is the formal adjoint of Ts, and Ts acts on compactly supported smooth functions supported in Ω. Then

T s + T s u = - div H ( V H u ) + s ( u φ - 1 H V H φ - V u φ - 2 H φ H φ + V u φ - 1 φ ) + s 2 V u φ - 2 H φ H φ .

Moreover, if T=T1, then

(3.1)

T + T u = - div H ( V H u ) + u φ - 1 H V H φ + V u φ - 1 φ ,
| T u | 2 = V φ 2 | H ( u φ - 1 ) | 2 .

Proof.

First compute

- div H ( V ( V H u ) ) = - div H ( V H u ) ,
- div H ( s V φ - 1 u H φ ) = - s ( φ - 1 u H φ H V - V u φ - 2 H φ H φ + V φ - 1 H u H φ + V φ - 1 u φ ) ,
- s V φ - 1 H φ ( V H u ) = - s V φ - 1 H φ H u ,
- s V φ - 1 H φ ( - s V φ - 1 H φ ) = s 2 V φ - 2 H φ H φ .

Collecting these identities together, we get

T s + T s u = - div H ( V H u ) + s ( u φ - 1 H V H φ - V u φ - 2 H φ H φ + V u φ - 1 φ ) + s 2 V u φ - 2 H φ H φ ,

as desired. Evaluating at s=1 provided the desired identity (3.1).

Next, by expanding

| T u | 2 = ( V H u - u V φ - 1 H φ ) ( V H u - u V φ - 1 H φ )
= V | H u | 2 - 2 V φ - 1 u H φ H u + u 2 V φ - 2 H φ H φ

we see that

V φ 2 | H ( u φ - 1 ) | 2 = V φ 2 ( φ - 1 H u - φ - 2 u H φ ) ( φ - 1 H u - φ - 2 u H φ )
= V | H u | 2 - 2 φ - 1 u H φ H u + φ - 2 u 2 H φ φ H φ
= | T u | 2 ,

which is the desired identity. ∎

Now, if there is a third positive C1-function W such that (V,W) is a Bessel pair with positive solution φ on an interval (0,R), then Lemma 1 may be applied to the functions V(d(x)) and φ(d(x)) to establish the following Hardy identities. We remark that similar results were obtained in [17, 26] where they assume the weights satisfy a general differential inequality; however, our approach with Bessel pairs provides a fairly simple way to produce many interesting weights.

Theorem 2.

Let G be a Carnot group with symmetric homogeneous gauge d, and let BR(0) denote the d-gauge ball of radius R for a given 0<R. Let V and W be positive C1-functions on (0,R). If (rQ-1V,rQ-1W) is a Bessel pair on (0,R) with positive solution φ, then for uC0(BR(0){0}) there holds

B R ( 0 ) V ( d ) | H ( u φ ( d ) ) | 2 φ ( d ) 2 𝑑 x = B R ( 0 ) V ( d ) | H u | 2 𝑑 x - B R ( 0 ) W ( d ) u 2 | H d | 2 𝑑 x .

In particular, the following Hardy’s inequality holds:

B R ( 0 ) V ( d ) | H u | 2 𝑑 x B R ( 0 ) W ( d ) | u | 2 | H d | 2 𝑑 x .

Proof.

First observe that

H V ( d ) = V ( d ) H d , H φ ( d ) = φ ( d ) H d .

Next, by using the factorization lemma (Lemma 1), we have

B R ( 0 ) V φ 2 | H ( u φ - 1 ) | 2 𝑑 x = B R ( 0 ) | T u | 2 𝑑 x
= B R ( 0 ) u T + T u 𝑑 x
= - B R ( 0 ) u div H ( V H u ) 𝑑 x + B R ( 0 ) u 2 φ - 1 H V H φ d x + B R ( 0 ) V u 2 φ - 1 φ 𝑑 x .

Next, we use the divergence theorem and Lemma A to obtain

B R ( 0 ) V φ 2 | H ( u φ - 1 ) | 2 𝑑 x = B R ( 0 ) V | H u | 2 𝑑 x + B R ( 0 ) u 2 V φ - 1 φ | H d | 2 𝑑 x
+ B R ( 0 ) V u 2 φ - 1 ( φ ′′ + Q - 1 d φ ) | H d | 2 𝑑 x
= B R ( 0 ) V | H u | 2 𝑑 x + B R ( 0 ) u 2 φ - 1 ( V φ ′′ + Q - 1 d V φ + V φ ) | H d | 2 𝑑 x .

Then, using that (rQ-1V,rQ-1W) is a Bessel pair, we have that

Q - 1 r V ( r ) φ ( r ) + V ( r ) φ ( r ) + V ( r ) φ ′′ ( r ) = - W ( r ) φ ( r ) ,

and so

B R ( 0 ) V φ 2 | H ( u φ - 1 ) | 2 𝑑 x = B R ( 0 ) V | H u | 2 𝑑 x - B R ( 0 ) u 2 φ - 1 W | H d | 2 𝑑 x ,

which is the desired identity. From

B R ( 0 ) V ( d ) | ( u φ ( d ) ) | 2 φ ( d ) 2 𝑑 x 0

the last inequality of the theorem follows. ∎

Remark 3.

One of the main strengths of the identities given in Theorem 2 is that it allows one to check for the existence of extremizers of

B R ( 0 ) V ( d ) | H u | 2 𝑑 x B R ( 0 ) W ( d ) u 2 | H d | 2 𝑑 x

provided that this inequality is in fact sharp. Indeed, if

B R ( 0 ) V ( d ) | H u | 2 𝑑 x = B R ( 0 ) W ( d ) u 2 | H d | 2 𝑑 x ,

then

B R ( 0 ) V ( d ) | H ( u φ ( d ) ) | 2 φ ( d ) 2 𝑑 x = 0 .

This implies that

u = φ ( d ) ,

and so the existence of extremizers is dictated by the integrability of φ in the appropriate weighted Sobolev space.

Remark 4.

Throughout, we assume that uC0(G{0}). This assumption may be weakened to assuming uC0(G) provided that the origin is a removable singularity with respect to the weighted Sobolev norm

u 2 = B R ( 0 ) V ( d ) | H u | 2 𝑑 x + B R ( 0 ) W ( d ) | u | 2 | H d | 2 𝑑 x .

On the other hand, if either V or W is too singular, then uC0(G{0}) may be necessary for integrability.

As mentioned above, by choosing specific Bessel pairs and using Theorem 2, we may improve geometric inequalities in the sub-elliptic setting. We now state several such applications and include their respective proofs which depict the appropriate choice of Bessel pair.

The first identity we present is the Hardy identity with power weights which imply the classical weighted Hardy inequalities on Carnot groups.

Corollary 5.

For all uC0(G{0}), and all λR, there holds

G d - λ | H u | 2 d x - ( Q - 2 - λ ) 2 4 G d - λ - 2 | u | 2 | H d | 2 d x = G d 2 - Q | H ( d Q - 2 - λ 2 u ) | 2 d x .

Proof.

The identity follows from applying Theorem 2 to the Bessel pair

( r Q - 1 V , r Q - 1 W ) = ( r Q - 1 - λ , ( Q - 2 - λ ) 2 4 r Q - 3 - λ ) ,

which has as a positive solution φ(r)=r-Q-2-λ2. ∎

The following corollary establishes the Hardy identity which implies the Brezis–Vázquez-type inequality in our setting.

Theorem 6.

Let G be a Carnot group with symmetric homogeneous gauge d, and let BR(0) denote the d-gauge ball of radius R for a given 0<R<. Then, for uC0(BR(0){0}), there holds

B R ( 0 ) | H u | 2 𝑑 x - ( Q - 2 2 ) 2 B R ( 0 ) | u | 2 d 2 𝑑 x = z 0 2 R 2 B R ( 0 ) | u | 2 𝑑 x + B R ( 0 ) | H u d Q 2 J 0 ; R ( d ) | 2 J 0 ; R ( d ) 2 d Q 𝑑 x .

Proof.

The identity follows from applying Theorem 2 to the Bessel pair

( r Q - 1 V , r Q - 1 W ) = ( r Q - 1 , r Q - 1 ( ( N - 2 2 ) 2 r - 2 + z 0 2 R 2 ) )

which has as a positive solution

φ ( r ) = r Q 2 J 0 ( r z 0 R ) = r - Q 2 J 0 ; R ( r ) .

3.2 Hardy Identities on the Heisenberg Half Space

Define the Heisenberg half-space HN+={(z,t)HN:t>0} and the half-gauge ball BR+(0)=BR(0)HN+. In this subsection, we introduce general Hardy identities in terms of Bessel pairs on the half-space HN+ and half-gauge ball BN+, and, for specific choices of the Bessel pair, we recover and find exact remainder terms for the half-space inequalities given by Luan and Yang in [23]. To begin, first define the strip

A 2 R = { ( z , t ) H N : t ( 0 , 2 R ) } ,

and define the vertical height function dR(x,t)=min{t,2R-t}, and set d(z,t)=t. Using Lemma 1, we will obtain Hardy inequalities on A2R in this section.

To obtain Hardy identities in terms of the vertical derivative t, the following factorization lemma specialized to the vertical derivative T=t on HN is needed. In fact, the following factorization is effectively used in [20], but for completion and to state the factorization in the same generality as Lemma 1, we state and prove the following lemma.

Lemma 7.

Let

S = V t - V φ - 1 φ t , S + = - t ( V ) - V φ - 1 φ t ,

where S+ is the formal adjoint of S, and S acts on compactly supported functions with appropriate domains. Then

S + S u = - t ( V u t ) + φ - 1 V t φ t u + V φ - 1 2 φ t 2 u , | S u | 2 = V φ 2 | t ( u φ ) | 2 .

Proof.

First compute

t ( V φ - 1 φ t u ) = φ - 1 V t φ t u - φ - 2 V ( φ t ) 2 u + V φ - 1 2 φ t 2 u + V φ - 1 φ t u t ,

from which we obtain

S + S u = - t ( V u t ) + φ - 1 V t φ t u - φ - 2 V ( φ t ) 2 u + V φ - 1 2 φ t 2 u + V φ - 1 φ t u t
- V φ - 1 φ t u t + V φ - 2 ( φ t ) 2 u
= - t ( V u t ) + φ - 1 V t φ t u + V φ - 1 2 φ t 2 u ,

as desired. Lastly, computing

| S u | 2 = V ( u t ) 2 + V φ - 2 ( φ t ) 2 ( u t ) 2 - 2 V φ - 1 φ t u

and

| t ( u φ ) | 2 = φ - 2 ( u t ) 2 + φ - 4 u 2 ( φ t ) 2 - 2 φ - 3 u φ t u t ,

we find

| S u | 2 = V φ 2 | t ( u φ ) | 2 .

We may now state the main Hardy identity for strips.

Theorem 8.

Let 0<R, and let V and W be positive C1-functions on (0,R). If (rQ+1V,rQ+1W) is a Bessel pair on (0,R) with positive solution φ, then for uC0(A2R), there holds

A 2 R V ( d R ) | H ( u φ ( d R ) ) | 2 φ 2 ( d R ) 𝑑 z 𝑑 t = A 2 R V ( d R ) | u | 2 𝑑 z 𝑑 t - 1 4 A 2 R W ( d R ) | z | 2 | u | 2 𝑑 z 𝑑 t

and

A 2 R V φ 2 | t ( u φ ) | 2 𝑑 z 𝑑 t = A 2 R V ( d R ) | u t | 2 𝑑 z 𝑑 t - A 2 R W ( d K ) | u | 2 𝑑 z 𝑑 t ,

Proof.

First observe that

φ ( d R ) = φ ′′ ( d R ) | z | 2 , | H d R | 2 = 1 4 | z | 2 ,

and so, by Lemma 1, we have that

T + T u = - div H ( V H u ) + 1 4 u φ ( d R ) φ ( d R ) | z | 2 V ( d R ) + 1 4 u V ( d R ) φ ′′ ( d R ) φ ( d R ) | z | 2
= - div H ( V H u ) + 1 4 u | z | 2 1 φ ( d R ) [ φ ( d R ) V ( d R ) + φ ′′ ( d R ) V ( d R ) ] .

Therefore, using that (V,W) is a Bessel pair and the divergence theorem, we find

A 2 R V ( d R ) | H ( u φ ( d R ) ) | 2 φ 2 ( d R ) 𝑑 z 𝑑 t = A 2 R | T u | 2 𝑑 z 𝑑 t
= A 2 R V ( d R ) | H u | 2 𝑑 z 𝑑 t - 1 4 A 2 R W ( d R ) | z | 2 | u | 2 𝑑 z 𝑑 t ,

as desired.

To obtain the second identity, one may simply appeal to [20, Theorem 1.3]; we include the following short proof for completeness. Observe

V t φ t = V ( d K ) φ ( d K ) , 2 φ t 2 = φ ′′ ( d K ) ,

and so, using Lemma 7, we find

S + S u = - t ( V u t ) + u φ - 1 [ V φ + V φ ′′ ] .

Therefore, using that (V,W) is a Bessel pair and integration by parts, we find

A 2 R V φ 2 | t ( u φ ) | 2 𝑑 z 𝑑 t = A 2 R | S u | 2 𝑑 z 𝑑 t = A 2 R V ( d R ) | u t | 2 𝑑 z 𝑑 t - A 2 R W ( d K ) | u | 2 𝑑 z 𝑑 t ,

as desired. ∎

Using Theorem 8, we may obtain the following Hardy identity which improves and recovers (1.5).

Corollary 9.

For uC0(HN+), there holds

H N + | H t α 2 u | 2 t 𝑑 z 𝑑 t = H N + t α + 1 | H u | 2 𝑑 z 𝑑 t - α 2 16 H N + | z | 2 t α - 1 | u | 2 𝑑 z 𝑑 t

and

H N + | t ( t α 2 u ) | 2 t 𝑑 z 𝑑 t = H N + t α + 1 | u t | 2 𝑑 z 𝑑 t - α 2 4 H N + t α - 1 | u | 2 𝑑 z 𝑑 t .

In particular, we have the algebraically interesting Hardy inequalities

H N + t α + 1 | T u | 2 𝑑 z 𝑑 t = H N + t α + 1 | [ X j , Y j ] u | 2 𝑑 z 𝑑 t α 2 4 H N + t α - 1 | u | 2 𝑑 z 𝑑 t .

Proof.

The identity follows from applying Theorem 8 to the Bessel pair

( V , W ) = ( r α + 1 , α 2 4 r α - 1 )

with positive solution φ=r-α2 on (0,). ∎

Remark 10.

We point out that to obtain Hardy identity from [23], one needs to take α=-1, and the change of variable t4t is needed since they use the vector fields

X j = x j - 2 y j t , Y j = y j + 2 x j t ,

compare this with (2.1).

We now state and prove the Hardy identities for the half ball BR+(0). We will again need a new factorization lemma.

Lemma 11.

Let V,φ,Y be sufficiently smooth and define

U = V H - V φ - 1 H φ + V Y - 1 H Y ,
U + = - div H ( V ) - V φ - 1 H φ + V Y - 1 H Y .

where U+ is the formal adjoint of U, and U acts on compactly supported smooth functions with appropriate domains. Then

U + U u = - div H ( V H u ) + u φ - 1 H V H φ + V u φ - 1 Δ φ + 2 V u Y - 2 | H Y | 2 - V u Y - 1 Δ Y
- u Y - 1 H V H Y - 2 V u φ - 1 Y - 1 H φ H Y ,
| U u | 2 = V φ 2 Y - 2 | H ( u φ - 1 Y ) | 2 .

Proof.

If T denotes the operator in Lemma 1, then

U = T + V Y - 1 H Y , U + = T + + V Y - 1 H Y .

Consequently, there holds

(3.2)

U + U u = ( T + + V Y - 1 H Y ) ( T + V Y - 1 H Y ) u
= T + T u + T + ( u V Y - 1 H Y ) + V Y - 1 H Y ( T u ) + u V Y - 2 | H Y | 2 .

Since we already have T+T from Lemma 1, we need only compute the middle two terms

(3.3)

T + ( u V Y - 1 H Y ) = - div H ( V u Y - 1 H Y ) - u V Y - 1 φ - 1 H φ H Y
= u Y - 1 H V H Y - V Y - 1 H u H Y + V u Y - 2 | H Y | 2 - V u Y - 1 Y
- u V Y - 1 φ - 1 H φ H Y ,
V Y - 1 H Y ( T u ) = V Y - 1 H Y H u - u V φ - 1 Y - 1 H Y H φ .

Therefore, using that

T + T u = - div H ( V H u ) + u φ - 1 H V H φ + V u φ - 1 φ ,

we may combine (3.2) and (3.3) to obtain the desired identity for U+U.

The computation for |Uu|2 is similar to the computation in Lemma 1. ∎

We remark that Lemma 1 may be obtained from Lemma 7 by taking Y=0.

With this factorization lemma established, we are now prepared to state and prove the following Hardy identity on the half-gauge ball BR+(0).

Theorem 12.

Let 0<R, and let V and W be positive C1-functions on (0,R). If (rQ+1V,rQ+1W) is a Bessel pair on (0,R), then for uC0(BRHN+) there holds

B R ( 0 ) H N + V | ( u φ 1 t ) | 2 φ 2 t 2 𝑑 z 𝑑 t = B R ( 0 ) H N + V | u | 2 𝑑 z 𝑑 t - B R ( 0 ) H N + [ W - V d ] | u | 2 | d | 2 𝑑 z 𝑑 t .

Proof.

This theorem will follow by taking Y=1t in Lemma 11. To see this, first compute

1 t = | z | 2 2 t 3 , H d H 1 t = - 1 t 2 H d H t = - 1 t | H d | 2 .

We then compute

u φ - 1 H V H φ = u φ - 1 φ V | H d | 2 ,
V u φ - 1 φ = V u φ - 1 | H d | 2 ( φ ′′ + Q - 1 d φ ) ,
2 V u t 2 | H 1 t | 2 = 1 2 V u | z | 2 t 2 ,
V u t t = 1 2 V u | z | 2 t 2 ,
u t H V H 1 t = - u V | H d | 2 ,
2 V u φ - 1 t H φ H 1 t = - 2 V u φ - 1 φ | H d | 2 .

Collecting everything together, we have

U + U u = - div H ( V H u ) + u V | H d | 2 + u φ - 1 | H d | 2 [ V φ ′′ + V φ + Q + 1 d V φ ]
= - div H ( V H u ) + u V | H d | 2 + u W | H d | 2 ,

where, to get the last equality, we used the fact that (rQ+1V,rQ+1W) is a Bessel pair. At last,

B R ( 0 ) H N + V | ( u φ 1 t ) | 2 φ 2 t 2 𝑑 z 𝑑 t = B R ( 0 ) H N + | U u | 2 𝑑 z 𝑑 t
= B R ( 0 ) H N + V | u | 2 𝑑 z 𝑑 t - B R ( 0 ) H N + [ W - V d ] | u | 2 | d | 2 𝑑 z 𝑑 t ,

as desired. ∎

3.3 Weighted Uncertainty Principles

Lemma 1 may also be used to obtain Caffarelli–Kohn–Nirenberg-type inequalities with more general weights than simple power weights. The following results are new even for the Euclidean setting.

Theorem 13.

If (rQ-1V,rQ-1W) is a Bessel pair on (0,R) with solution φ, then for uC0(BR(0){0}) there holds

1 2 B R ( 0 ) u 2 | H d | 2 [ W ( d ) + V ( d ) ( φ ( d ) φ ( d ) ) 2 ] 𝑑 x ( B R ( 0 ) u 2 | H d | 2 V ( d ) ( φ ( d ) φ ( d ) ) 2 𝑑 x ) 1 2
× ( B R ( 0 ) V ( d ) | H u | 2 𝑑 x ) 1 2 .

Proof.

The following proof is motivated by Costa’s work [6]. (See also [9, 8].) Using Lemma 1, we find

B R ( 0 ) | T s u | 2 𝑑 x = - B R ( 0 ) u ¯ div H ( V H u ) 𝑑 x + s 2 B R ( 0 ) | u | 2 V ( φ φ ) 2 | H d | 2 𝑑 x
+ s B R ( 0 ) | u | 2 [ V φ φ - V ( φ φ ) 2 + V 1 φ ( φ ′′ + Q - 1 d φ ) ] | H d | 2 𝑑 x .

Then, using that (rQ-1V,rQ-1W) is a Bessel pair, we find

B R ( 0 ) | T s u | 2 𝑑 x = B R ( 0 ) | H u | 2 V 𝑑 x + s 2 B R ( 0 ) | u | 2 V ( φ φ ) 2 | H d | 2 𝑑 x - s B R ( 0 ) | u | 2 [ W + V ( φ φ ) 2 ] | H d | 2 𝑑 x .

At last, by using that As2+Bs+C0 for all s implies 12|B|A12C12, with

A = B R ( 0 ) | u | 2 V ( φ φ ) 2 | H d | 2 𝑑 x ,
B = - B R ( 0 ) u 2 | H d | 2 [ W + V ( φ φ ) 2 ] 𝑑 x ,
C = B R ( 0 ) V | H u | 2 𝑑 x ,

we obtain

1 2 B R ( 0 ) u 2 | H d | 2 [ W + V ( φ φ ) 2 ] 𝑑 x ( B R ( 0 ) u 2 | H d | 2 V ( φ φ ) 2 𝑑 x ) 1 2 ( B R ( 0 ) V | H u | 2 𝑑 x ) 1 2 ,

as desired. ∎

By taking

( V , W ) = ( r λ , | Q + λ - 1 | r λ - 2 ) ,

then (rQ-1V,rQ-1W) is a Bessel pair with positive solution φ=r1-λ-Q on (0,). Using this Bessel pair in Theorem 13, we recover the weighted Hardy uncertainty principles as a corollary

| Q + λ - 2 | 2 G u 2 d λ | H d | 2 𝑑 x ( G u 2 d λ + 2 | H d | 2 𝑑 x ) 1 2 ( G | H u | 2 d λ 𝑑 x ) 1 2

for all λ and uC0(G{0}). In fact, one may show that if

V = r - 2 b , V ( φ φ ) 2 = α r - 2 a

were to hold for some a,b,α, then W changes sign on (0,) unless a=b+1. Therefore, it is reasonable to interpret Corollary 14 as a general family of uncertainty principles.

Now, using that W>0 and V(φφ)20, Theorem 13 may be used to obtain general weighted Hardy-type inequalities of the following form.

Corollary 14.

Suppose V>0 and φ0 are C1-functions on (0,R) such that

( r Q - 1 V φ ) < 0 on  ( 0 , R ) .

Then, for all uC0(BR(0){0}), there holds

1 4 B R ( 0 ) | u | 2 V ( d ) ( φ ( d ) φ ( d ) ) 2 | H d | 2 𝑑 x B R ( 0 ) | H u | 2 V ( d ) 𝑑 x .

Proof.

If (rQ-1Vφ)<0 for a given pair (V,φ), then taking

W = - 1 r Q - 1 φ ( r Q - 1 V φ ) ,

we obtain a Bessel pair (rQ-1V,rQ-1W). Since W>0 and V(φφ)20, there holds

B R ( 0 ) u 2 | H d | 2 [ W + V ( φ φ ) 2 ] 𝑑 x B R ( 0 ) u 2 | H d | 2 V ( φ φ ) 2 𝑑 x ,

and so Theorem 13 implies

1 2 B R ( 0 ) u 2 | H d | 2 V ( φ φ ) 2 𝑑 x ( B R ( 0 ) u 2 | H d | 2 V ( φ φ ) 2 𝑑 x ) 1 2 ( B R ( 0 ) V | H u | 2 𝑑 x ) 1 2 .

Rearranging terms gives the desired inequality. ∎

In the case that V=1, Corollary 14 may be used to obtain the following sufficient condition for a certain Hardy-type inequality to hold. For simplicity, we state it for Euclidean space only.

Corollary 15.

Let w(r) be a nonnegative C1-function on (0,R). Then, for

1 4 B R ( 0 ) | u | 2 w ( | x | ) 𝑑 x B R ( 0 ) | u | 2 𝑑 x

to hold for uC0(BR(0){0}), it is sufficient for w to satisfy the differential inequality

(3.4) ( r N - 1 w 1 2 e w 1 2 𝑑 r ) < 0 .

This inequality is equivalent to

(3.5) N - 1 r w + 1 2 w + w 3 2 < 0 .

Proof.

Setting w=(φφ)2, we see that

φ = C e w 1 2 𝑑 r , φ = C w 1 2 e w 1 2 𝑑 r

for some C>0. The result follows from Corollary 14.

To see why (3.5) is equivalent, first observe that

φ = w 1 2 φ , φ ′′ = 1 2 w - 1 2 φ + w 1 2 φ = 1 2 w - 1 2 φ + w φ .

Consequently, (3.4) holds if and only if

0 > ( r N - 1 φ ) = ( N - 1 ) r N - 2 φ + r N - 1 φ ′′ = ( N - 1 ) r N - 2 w 1 2 φ + 1 2 w - 1 2 φ + w φ ,

and whence (3.5). ∎

It is interesting to note that the Hardy potential w(r)=(N-2)2r-2 does not satisfy (3.5). However, the class of weights satisfying (3.5) is not small. Indeed, for every Bessel pair (V,W) with positive solution φ, one may simply take w=(φ/φ)2. It is also worth mentioning that the more singular w is at the origin, the more likely (3.5) is to hold.

3.4 General Distance Functions in Carnot–Carathéodory Spaces

We first recall a result for a general class of vector fields given in [24] by Monti and Serra-Cassano. In preparation, we state three cases in which their result (and hence our results) apply. Note that Case A is satisfied by any Carnot group, and hence the Heisenberg group.

Case A.

Let X1,,XmC(N;N), m<N, satisfy Hörmander’s condition and be of the form

X j = x j + i = m + 1 N a i j ( x ) x i , j = 1 , , m ,

where aijC(N).

Case B.

Let X1,,XNC(N;N) be of the form

X 1 = x 1 , X 2 = p 2 ( x 1 ) x 2 , , X N = p N ( x 1 , , x N - 1 ) x N ,

where pjC(j-1), j=1,,N, are functions vanishing on a set of null (j-1)-dimensional Lebesgue measure.

Case C.

Let X1,,XmC(N;N) and these vector fields span N at every point.

Their theorem may now be stated.

Theorem H.

Let (RN,d) be a CC complete space induced by X1,,Xm belonging to any one of the cases A, B or C. Let KRN be a closed set, and let dK denote the CC distance from K. Then dK satisfies the horizontal eikonal equation

| H d K ( x ) | = 1

for almost every xRNK.

In particular, if K={0} is the origin, and if dcc=d{0} is the CC-distance to the origin, then dcc satisfies the eikonal equation |Hdcc|=1, in contrast with the homogeneous gauges considered above. Using this fact and Lemma 1, we may establish the following Hardy identity theorem for general distance functions arising from cases A, B and C above.

Theorem 16.

Let (RN,d) be a CC complete space induced by X1,,Xm satisfying any one of the cases A, B or C. Let KRN be a closed set, and let dK denote the CC distance from K. Let ΩR={xH:0<dK(x)<R} for a given 0<R, and let V and W be positive C1-functions on (0,R). If (rα-1V,rα-1W) is a Bessel pair on (0,R) with positive solution φ for some αR, then for uC0(ΩR{0}) there holds

Ω R V ( d K ) | H ( u φ ( d K ) ) | 2 φ ( d K ) 2 𝑑 x = Ω R V ( d K ) | H u | 2 𝑑 x - Ω R W ( d K ) u 2 𝑑 x
- Ω R V ( d K ) | u | 2 φ - 1 ( d K ) φ ( d K ) [ d K - α - 1 d K ] 𝑑 x .

In particular, if

(3.6) φ ( d K ) [ d K - α - 1 d K ] 0

in the sense of distributions, then

Ω R V ( d K ) | H u | 2 𝑑 x Ω R W ( d K ) | u | 2 𝑑 x .

Proof.

First observe that

(3.7)

H V ( d K ) = V ( d K ) H d K ,
H φ ( d K ) = φ ( d K ) H d K ,
H V ( d K ) H φ ( K ) = V ( d K ) φ ( d K ) ,
φ ( d K ) = j = 1 m X j 2 φ ( d K ) = φ ′′ ( d K ) + φ ( d K ) d K ,

where we have invoked Theorem H to use |HdK|=1. Now, by using the factorization lemma (Lemma 1), we have

Ω R V φ 2 | H ( u φ - 1 ) | 2 𝑑 x = Ω R | T u | 2 𝑑 x
= Ω R u T + T u 𝑑 x
= - Ω R u div H ( V H u ) 𝑑 x + Ω R u 2 φ - 1 H V H φ d x + Ω R V u 2 φ - 1 φ 𝑑 x .

Next, we use the divergence theorem and (3.7) to obtain

Ω R V φ 2 | H ( u φ - 1 ) | 2 𝑑 x = Ω R V | H u | 2 𝑑 x + Ω R u 2 V φ - 1 φ 𝑑 x + Ω R V u 2 φ - 1 ( φ ′′ + φ d K ) 𝑑 x
= Ω R V | H u | 2 𝑑 x + Ω R u 2 φ - 1 ( V φ ′′ + α - 1 d V φ - α - 1 d V φ + φ d K + V φ ) 𝑑 x .

Then, using that (rα-1V,rα-1W) is a Bessel pair, we have that

α - 1 r V ( r ) φ ( r ) + V ( r ) φ ( r ) + V ( r ) φ ′′ ( r ) = - W ( r ) φ ( r ) ,

and so

Ω R V ( d K ) φ 2 | H ( u φ - 1 ( d K ) ) | 2 𝑑 x = Ω R V ( d K ) | H u | 2 𝑑 x - Ω R u 2 W ( d K ) 𝑑 x
- Ω R V ( d K ) | u | 2 φ - 1 ( d K ) φ ( d K ) [ d K - α - 1 d K ] 𝑑 x .

which is the desired identity. ∎

It is also worth mentioning that verifying if (3.6) holds is difficult in general. So, for sake of concreteness, we state applications for Theorem 16 for the particular case when the vector fields are the left-invariant vector fields on the Heisenberg group HN, and the distance is the CC distance to the center Z of HN (this is shown below). We note that, in this case, ΩR is the infinite cylinder of radius R with axis being Z. First we establish the general Hardy identity for this setting, and then its applications.

Corollary 17.

Let Z={(z,t)H:z=0} denote the center of the Heisenberg group H=R2N+1, and let dZ(x) denote the CC distance from a point xH to Z. Let ΩR={xH:dZ(x)<R} for a given 0<R, and let V and W be positive C1-functions on (0,R). If (r2N-1V,r2N-1W) is a Bessel pair on (0,R) with positive solution φ for some αR, then for uC0(ΩRZ) there holds

(3.8) Ω R V ( | z | ) | H ( u φ ( | z | ) ) | 2 φ ( | z | ) 2 𝑑 x = Ω R V ( | z | ) | H u | 2 𝑑 x - Ω R W ( | z | ) u 2 𝑑 x .

Proof.

First, the CC distance dZ(x) from a point xHN to the Heisenberg group center Z is computed. To begin, let deu(x,y) denote the Euclidean distance between two points x,y2N+1 when considered as an Euclidean space, and let

d eu , Z ( x ) = inf { d eu ( x , y ) : y Z }

denote the Euclidean distance from a point x2N+1 to the Heisenberg group center Z. It is well known that the CC distance from a point (z,0)HN to the origin (0,0) is given by

d cc ( ( z , 0 ) , ( 0 , 0 ) ) = | z | .

Due to left-invariance, the distance from any point (z,t) to (0,t) is thus |z|.

Now, it is also well known that, for two points x,y2N+1, there holds deu(x,y)dcc(x,y), and so deu,Z(x)dZ(x). Given that

d cc ( ( z , t ) , ( 0 , t ) ) = | z | = d eu , Z ( ( z , t ) , ( 0 , t ) ) d Z ( z , t ) ,

it follows that dZ(z,t)=|z| for all (z,t)HN. In particular, ΩR is the infinite cylinder of radius R with axis being Z.

Next, using (2.2), we compute

d K ( z , t ) = Δ 2 N | z | = 2 N - 1 | z | = 2 N - 1 d cc ( z , t ) .

It is also easy to verify directly that |HdZ|=1. With these observations, we find that

- Ω R V ( d cc ) | u | 2 φ - 1 ( d cc ) φ ( d cc ) [ d cc - 2 N - 1 d cc ] 𝑑 x = 0 .

The identity in (3.8) then follows by applying Theorem 16 with α=2N. ∎

We use this theorem to establish the following weighted Hardy identities for dZ.

Corollary 18.

For all uC0(HNZ), and all λR, there holds

H N d Z - λ | H u | 2 𝑑 x - ( 2 N - 2 - λ ) 2 4 H N d Z - λ - 2 | u | 2 𝑑 x = H N d Z 2 - 2 N | H ( d Z 2 N - 2 - λ 2 u ) | 2 𝑑 x .

Proof.

The identity follows from applying Corollary 18 to the Bessel pair

( r 2 N - 1 V , r 2 N - 1 W ) = ( r 2 N - 1 - λ , ( 2 N - 2 - λ ) 2 4 r Q - 3 - λ ) ,

which has as a positive solution φ(r)=r-2N-2-λ2. ∎

Recalling that on N, the sharp Hardy inequality takes the form

N | u | 2 𝑑 x ( N - 2 ) 2 4 N | u | 2 | x | 2 𝑑 x ,

and so it is apparent that the Hardy inequality on 2 reduces to the trivial inequality

N | u | 2 𝑑 x 0 .

Regardless of this, one may obtain critical Hardy inequalities of the form

Ω R 1 | x | α - 2 | u | 2 𝑑 x - 1 4 Ω R | u ( x ) | 2 | x | 2 | log R | x | | 2 𝑑 x = Ω R 1 | x | α - 2 log R | x | | ( u ( x ) log R | x | ) | 2 d x ;

see [21] and the references therein. Similar to the Euclidean case, if we take λ=2N-2 in Corollary 18, then we obtain trivial identities. Therefore, we establish analogous critical Hardy identities for the distance function dK to the center of HK.

Corollary 19.

For all uC0(ΩRZ), and all λR, there holds

Ω R d K 2 - 2 N | H u ( x ) | 2 𝑑 x - 1 4 Ω R | u ( x ) | 2 d Z ( x ) 2 N | log R d K ( x ) | 2 𝑑 x = Ω R 1 d K ( x ) 2 N - 2 | ( u ( x ) log R | x | ) | 2 log R d K ( x ) d x .

Proof.

The identity follows from applying Corollary 18 to the Bessel pair

( r 2 N - 1 V , r 2 N - 1 W ) = ( r , 1 4 r | log r R | 2 ) ,

i.e., with

( V , W ) = ( r 2 - 2 N , 1 4 r 2 N | log r R | 2 )

with solution φ=|logrR|. ∎

3.5 Hardy Identities for Nilpotent Groups and Hörmander Vector fields

In this subsection we state the analogous results for Theorem 2 in the more general setting of simply connected nilpotent groups and homogeneous Hörmander vector fields. The main point here is that the appropriate sub-Laplacians in these settings have a global fundamental solutions whose powers generalize the homogeneous gauges of Carnot groups. Thus, first consider a simply connected nilpotent Lie group G, where we specify a system of left invariant vector fields X1,,Xm which bracket generate the Lie algebra of G. Indeed, as was shown by Nagel, Ricci and Stein in [25], the sub-Laplacian =j=1mXj2 has a global fundamental solution Γ(x,y). Then, by taking dX(x)=Γ(x,0)2-q for some q, we find

( d X ) = | H d X | 2 q - 1 d X ,

where H=(X1,,Xm) is the horizontal gradient corresponding to the system X1,,Xm.

Similarly, if X={X1,,Xm} is a system of Hörmander vector fields homogeneous of degree -1, then, by [2], the Hörmander operator =j=1mXj2 admits a fundamental solution ΓX(x,y). Moreover, we also define dX(x)=ΓX(x,0)12-q and observe

( d X ) = | H d X | 2 q - 1 d X

also holds.

As was defined for general Hörmander vector fields in the preliminaries, we shall take for a suitable domain the sets

B X , R ( 0 ) = { x N : d X ( x ) < R } .

We now state the two main Hardy identities in this setting. Of course the Bessel pairs used in the previous section may be used in the following theorem, and so the analogous results also hold in this setting.

Theorem 20.

Let X1,,Xm denote either a system of homogeneous Hörmander vector fields on RN, or a bracket generating collection of left invariant vector fields on a simply connected nilpotent group. Next, let dX(x)=ΓX(x,0)12-q, where ΓX(x,y) is the fundamental solution of L=j=1mXj2, and let BX,R(0)={dX<R}. Lastly, let V and W be positive C1-functions on (0,R). If (rQ-1V,rQ-1W) is a Bessel pair on (0,R) with positive solution φ, then for uC0(BX,R{0}) there holds

(3.9) B X , R ( 0 ) V ( d X ) | H ( u φ ( d X ) ) | 2 φ ( d X ) 2 𝑑 x = B X , R ( 0 ) V ( d X ) | H u | 2 𝑑 x - B X , R ( 0 ) W ( d X ) u 2 | H d X | 2 𝑑 x .

In particular,

B X , R ( 0 ) V ( d X ) | H u | 2 𝑑 x B X , R ( 0 ) W ( d X ) | u | 2 | H d | 2 𝑑 x .

Proof.

The proof is identical to the proof of Theorem 2. ∎

Funding source: Simons Foundation

Award Identifier / Grant number: 519099

Funding statement: J. Flynn and G. Lu were partially supported by a grant from the Simons Foundation.

References

[1] S. Biagi and A. Bonfiglioli, The existence of a global fundamental solution for homogeneous Hörmander operators via a global lifting method, Proc. Lond. Math. Soc. (3) 114 (2017), no. 5, 855–889. 10.1112/plms.12024Search in Google Scholar

[2] S. Biagi, A. Bonfiglioli and M. Bramanti, Global estimates for the fundamental solution of homogeneous Hörmander operators, preprint (2019), https://arxiv.org/abs/1906.07836. 10.1007/s10231-021-01183-6Search in Google Scholar

[3] 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

[4] H. Brezis and M. Marcus, Hardy’s inequalities revisited. Dedicated to Ennio De Giorgi, Ann. Sc. Norm. Super. Pisa Cl. Sci. (4) 25 (1997), no. 1–2, 217–237. Search in Google Scholar

[5] H. Brezis and J. L. Vázquez, Blow-up solutions of some nonlinear elliptic problems, Rev. Mat. Univ. Complut. Madrid 10 (1997), 443–469. 10.5209/rev_REMA.1997.v10.n2.17459Search in Google Scholar

[6] D. G. Costa, Some new and short proofs for a class of Caffarelli–Kohn–Nirenberg type inequalities, J. Math. Anal. Appl. 337 (2008), no. 1, 311–317. 10.1016/j.jmaa.2007.03.062Search in Google Scholar

[7] N. T. Duy, N. Lam and G. Lu, p-Bessel pairs, Hardy’s identities and inequalities and Hardy–Sobolev inequalities with monomial weights, preprint. 10.1007/s12220-021-00847-2Search in Google Scholar

[8] J. Flynn, Sharp Caffarelli–Kohn–Nirenberg-type inequalities on Carnot groups, Adv. Nonlinear Stud. 20 (2020), no. 1, 95–111. 10.1515/ans-2019-2065Search in Google Scholar

[9] J. Flynn, Sharp L2-Caffarelli–Kohn–Nirenberg inequalities for Grushin vector fields, Nonlinear Anal. 199 (2020), Article ID 111961. 10.1016/j.na.2020.111961Search in Google Scholar

[10] J. Flynn, N. Lam and G. Lu, Lp-Hardy identities and inequalities for general distance functions, preprint. Search in Google Scholar

[11] J. Flynn, N. Lam and G. Lu, Lp-Hardy identities and inequalities on Carnot groups and for Hörmander’s vector fields, preprint. Search in Google Scholar

[12] J. Flynn, N. Lam and G. Lu, Lp Hardy identities and inequalities with respect to the distance and mean distance to the boundary, preprint. Search in Google Scholar

[13] V. Franceschi and D. Prandi, Hardy-type inequalities for the Carnot–Carathéodory distance in the heisenberg group, J. Geom. Anal. 31 (2021), 2455–2480. 10.1007/s12220-020-00360-ySearch in Google Scholar

[14] N. Garofalo and E. Lanconelli, Frequency functions on the Heisenberg group, the uncertainty principle and unique continuation, Ann. Inst. Fourier (Grenoble) 40 (1990), no. 2, 313–356. 10.5802/aif.1215Search in Google Scholar

[15] N. Ghoussoub and A. Moradifam, Bessel pairs and optimal Hardy and Hardy–Rellich inequalities, Math. Ann. 349 (2011), no. 1, 1–57. 10.1007/s00208-010-0510-xSearch in Google Scholar

[16] N. Ghoussoub and A. Moradifam, Functional Inequalities: New Perspectives and New Applications, Math. Surveys Monogr. 187, American Mathematical Society, Providence, 2013. 10.1090/surv/187Search in Google Scholar

[17] J. A. Goldstein, I. Kombe and A. Yener, A unified approach to weighted Hardy type inequalities on Carnot groups, Discrete Contin. Dyn. Syst. 37 (2017), no. 4, 2009–2021. 10.3934/dcds.2017085Search in Google Scholar

[18] J. Han, P. Niu and W. Qin, Hardy inequalities in half spaces of the Heisenberg group, Bull. Korean Math. Soc. 45 (2008), no. 3, 405–417. 10.4134/BKMS.2008.45.3.405Search in Google Scholar

[19] I. Kombe, Hardy, Rellich and uncertainty principle inequalities on Carnot groups, preprint (2006), https://arxiv.org/abs/math/0611850. Search in Google Scholar

[20] N. Lam, G. Lu and L. Zhang, Factorizations and Hardy’s type identities and inequalities on upper half spaces, Calc. Var. Partial Differential Equations 58 (2019), no. 6, Paper No. 183. 10.1007/s00526-019-1633-xSearch in Google Scholar

[21] N. Lam, G. Lu and L. Zhang, Geometric Hardy’s inequalities with general distance functions, J. Funct. Anal. 279 (2020), no. 8, Article ID 108673. 10.1016/j.jfa.2020.108673Search in Google Scholar

[22] J. Lehrbäck, Hardy inequalities and Assouad dimensions, J. Anal. Math. 131 (2017), 367–398. 10.1007/s11854-017-0013-8Search in Google Scholar

[23] J.-W. Luan and Q.-H. Yang, A Hardy type inequality in the half-space on n and Heisenberg group, J. Math. Anal. Appl. 347 (2008), no. 2, 645–651. 10.1016/j.jmaa.2008.06.048Search in Google Scholar

[24] R. Monti and F. Serra Cassano, Surface measures in Carnot-Carathéodory spaces, Calc. Var. Partial Differential Equations 13 (2001), no. 3, 339–376. 10.1007/s005260000076Search in Google Scholar

[25] A. Nagel, F. Ricci and E. M. Stein, Harmonic analysis and fundamental solutions on nilpotent Lie groups, Analysis and Partial Differential Equations, Lecture Notes Pure Appl. Math. 122, Dekker, New York (1990), 249–275. Search in Google Scholar

[26] P. Niu, H. Zhang and Y. Wang, Hardy type and Rellich type inequalities on the Heisenberg group, Proc. Amer. Math. Soc. 129 (2001), no. 12, 3623–3630. 10.1090/S0002-9939-01-06011-7Search in Google Scholar

[27] Q.-H. Yang, Hardy type inequalities related to Carnot–Carathéodory distance on the Heisenberg group, Proc. Amer. Math. Soc. 141 (2013), no. 1, 351–362. 10.1090/S0002-9939-2012-11322-XSearch in Google Scholar

[28] H. Zhang and P. Niu, Hardy-type inequalities and Pohozaev-type identities for a class of p-degenerate subelliptic operators and applications, Nonlinear Anal. 54 (2003), no. 1, 165–186. 10.1016/S0362-546X(03)00062-2Search in Google Scholar

Published Online: 2021-03-12
Published in Print: 2021-05-01

© 2021 Walter de Gruyter GmbH, Berlin/Boston

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

Downloaded on 19.4.2024 from https://www.degruyter.com/document/doi/10.1515/ans-2021-2123/html
Scroll to top button