Asymptotically mean value harmonic functions in sub-Riemannian and RCD settings

We consider weakly and strongly asymptotically mean value harmonic (amv-harmonic) functions on subriemannian and RCD settings. We demonstrate that, in non-collapsed RCD-spaces with vanishing metric measure boundary, Cheeger harmonic functions are weakly amv-harmonic and that, in Carnot groups, weak amv-harmonicity equivalently characterizes harmonicity in the sense of the sub-Laplacian. In homogeneous Carnot groups of step $2$, we prove a Blaschke-Privaloff-Zaremba type theorem. Similar results are discussed in the settings of Riemannian manifolds and for Alexandrov surfaces.


INTRODUCTION
A key problem in the theory of harmonic functions is to find various characterizations of harmonicity. A particularly strong example of this is the mean value property in the Euclidean setting. Unfortunately, harmonic functions rarely enjoy the mean value property outside Euclidean space. For example, in the setting of Riemannian manifolds it is known to hold on manifolds which are either flat or rank-one symmetric and in dimensions 2-5, see the related discussion on the Lichnerowicz conjecture in Example 4 in [1] and references therein. However a suitable variant beyond the Euclidean setting turns out to be the asymptotic mean value (amv) property In Riemannian manifolds the convergence in (1) is locally uniform for a harmonic function u. However the r-laplacian in (1) can be defined for locally integrable functions on any metric measure space X = (X, d, µ), where u Br(x) stands for the mean-value of u over a ball B r (x) with radius r > 0 centered at x ∈ X. Consequently the amv property can be formulated in a very general setting (note that the mean value (mv) property can be also expressed in the same generality in terms of vanishing r-laplacian, see see Definition 1.1 and also [1,27]). Let us briefly mention that the r-laplacian arises in approximation problems of Riemannian manifolds by graphs [18], in geometric group theory [34], and in relation with functions of bounded variation [22,Theorem 1]. It moreover appears in connection with p-harmonic type equations and stochastic games [40,41,9], as well as in the setting of Carnot groups; the r-laplacian and its relations to the subelliptic harmonic functions have been studied, for instance, in [3,25,26].
For recent work on amv functions on general metric measure spaces we refer the reader to [42,43] and in particular to [2] where we focus on functions with the amv property in doubling metric measure spaces and show their local Hölder regularity and the mv property for their blow-ups. As an important special case, in [2] we consider weighted Euclidean setting and find an elliptic PDE satisfied by amv-harmonic functions.
1.1. Main results. In this paper we continue the investigations of and complement results from [2] by considering functions with the amv property in subriemannian and RCD settings. As in [2] we consider two notions of asymptotically mean value harmonic (amv-harmonic) functions, arising from different ways to interpret the limit in (1). More precisely, we define strong and weak amv-harmonicity as follows. Definition 1.1 (see Definition 1.1 in [2]). Let X = (X, d, µ) be a metric measure space.
(1) A function u ∈ L 1 loc (X) is strongly amv-harmonic, if lim r→0 ∆ r u L ∞ (K) = 0 for any compact set K ⊂ X; (2) A function u ∈ L 2 loc (X) on a metric measure space X is weakly amv-harmonic if lim r→0 X ϕ∆ r udµ = 0 for every compactly supported Lipschitz function ϕ.
To state our first result we refer to Section 3.3 and the references therein for the theory of RCDspaces, which include GH-limits of Riemannian manifolds with uniform Ricci lower bounds and dimension upper bounds (i.e. Ricci limit spaces) as well as Alexandrov spaces equipped with the Hausdorff measure of appropriate dimension.
We refer to Definition 3.2 for the boundary ∂X of an RCD(K, N )-space X = (X, d, µ) and recall X is non-collapsed if µ = H N , the N -Hausdorff measure. Below and in what follows, harmonic functions mean continuous local minimizers of the Cheeger energy, cf. Definition 2.1.
In particular, Theorem 1.1 implies that harmonic functions on Alexandrov spaces without boundary are weakly amv-harmonic. For locally Lipschitz continuous functions, the converse implication in Theorem 1.1 also holds, see Remark 3.3.
It is not known whether harmonic functions are also strongly amv-harmonic in the setting of Theorem 1.1 (see [44,Corollary 3.5] for a related weaker statement). Such a result would need very good control over the rate of change of the H N -measure of balls centered at nearby points: the presence of a weight (e.g. the setting of weakly non-collapsed RCD-spaces, see [24]) rules out the possibility of having even the weak amv property for harmonic functions, see [2,Section 6].
In Carnot groups, homogeneity properties come into play and ensure the good control needed for the strong amv property. Using a second order Taylor expansion our next theorem establishes the equivalence of harmonic and strongly amv-harmonic functions in Carnot groups equipped with homogeneous pseudonorms with sufficiently regular unit balls.
To state the theorem recall that an s-step Carnot group G with Lie algebra g = V 1 ⊕ . . . ⊕ V s can be identified with R m , where m = v 1 + · · · + v s and v i := dim V i , for i = 1, . . . , s (see Section 3.2 for more details on Carnot groups). We denote by z = (z (1) , z (2) , . . . , z (s) ) the coordinates in G, where z (i) ∈ R v i . We consider a pseudonorm ρ on G such that for some F : R × R m−v 1 → R that is C 1 with ∂ 1 F > 0 outside the origin.

Theorem 1.2.
Let Ω ⊂ G be an open domain in a Carnot group G of homogeneous dimension Q equipped with a pseudonorm ρ satisfying (2). If Ω denotes the metric measure space (Ω, ρ, H Q ), where H Q is the Haar measure on G, then the following are equivalent for a function u ∈ HW 1,2 loc (Ω). (i) u is strongly amv-harmonic; (ii) u is weakly amv-harmonic; (iii) u satisfies ∆ G u = 0 in Ω.
In the claim above HW 1,2 loc (Ω) denotes the horizontal Sobolev space and ∆ G is the subelliptic laplacian (13). Theorem 1.2 characterizes amv-harmonicity (in both senses) for a large class of pseudonorms, namely those satisfying (2), and in particular for the Koranyi gauge (12). See [26] for similar results for p-harmonic functions on Carnot groups with the Koranyi gauge. We remark that solutions to the sub-Laplacian equation are known to satisfy certain weighted mean value properties, see [11,Chapter 5.5] and [3].
We close this subsection by restricting ourselves to Carnot groups of step s = 2 where we are able to obtain a Blaschke-Privaloff-Zaremba (BPZ) type theorem. Such results update regularity from continuous to smooth under the assumption that the r-laplacian vanishes pointwise as r → 0. In the theorem below we consider the Koranyi gauge (12) on a Carnot group.

Theorem 1.3.
Let Ω ⊂ G be a domain in the Carnot group G of step 2 equipped with the gauge distance d G and let f : Ω → R be a continuous function satisfying Then f is a sub-elliptic harmonic function, i.e. for all x ∈ Ω it holds that where C depends on v 1 (the dimension of the horizontal layer V 1 ) and the geometry of G.
Similar results are well-known on Riemannian manifolds (cf. Proposition 4.2 for a sketch of proof for the readers convenience). We also obtain BPZ-type results for Alexandrov surfaces (cf. Proposition 4.4), and for Alexandrov spaces (cf. Proposition 4.5). In summary, we show that continuity and the strong amv-harmonicity are sufficient conditions to imply the harmonicity with respect to various natural Laplace operators in some singular spaces.
1.2. Overview and strategy. One of the drawbacks of the r-laplacian ∆ r is its lack of self-adjointness: without further assumptions on the geometry of the underlying space, the difference ∆ r − ∆ * r of the r-laplacian and its formal adjoint (cf. Section 2.2) does not tend to zero in any reasonable sense. In the RCD-setting, where the translation invariance of the measure is not available, a large part of our work is to deal with this issue. We overcome it by introducing the symmetrized r-laplacian of a locally integrable function u defined in a metric measure space X = (X, d, µ) (see Definition 10 for details): In Section 3.1 we study the basic properties of the symmetrized r-laplacian, in particular self-adjointness, which makes it a viable candidate for a yet another laplacian studied beyond the Euclidean setting. Crucially, the symmetrized r-laplacian is connected to the (approximate) Korevaar-Schoen (KS) energy so that, if the KS-energy is a Dirichlet form in the sense of Definition 2.2, then the symmetrized r-laplacians ∆ r u converge in a suitable weak sense, as r → 0, to the Korevaar-Schoen laplacian ∆ KS u of u for any u ∈ D(∆ KS ), see Section 2.3 and Proposition 3.2. It follows that local minimizers of the Korevaar-Schoen energy are characterized as functions whose symmetrized r-laplacian vanishes in the weak sense as r → 0.
Both RCD-spaces and Carnot groups have the property that the Korevaar-Schoen energy is a Dirichlet form. For RCD-spaces, this follows from [28] (in the non-collapsed case also from work of Sturm [52,53]), while for Carnot groups it is implied by [54]. See Propositions 3.6 and 3.3 for precise statements. For Carnot groups (where ∆ r = ∆ r ) this observation leads to the equivalence of weak amv-harmonicity and subelliptic harmonicity. We remark that the novel implication in Theorem 1.2 is that weakly harmonic functions are subelliptic harmonic; indeed, the strong amv property of subelliptic harmonic functions follows with a straightforward calculation as in [26,Lemma 4.1].
In the RCD-setting the equality ∆ r = ∆ r breaks down and, in the presence of boundary, the two rlaplacians do not agree even asymptotically (which can be seen by considering e.g. the standard half plane). To circumvent this problem we use the notion of vanishing metric measure boundary introduced in [32], which has recently been shown in [16] to be equivalent to having no boundary in the sense of [24]. In non-collapsed RCD-spaces the vanishing mm-boundary enables us to obtain the asymptotic vanishing of ∆ r − ∆ r as r → 0 in a suitable weak sense, cf. Lemma 3.8 and Corollary 3.9. This yields the proof of Theorem 1.1. Theorem 1.3, and the BPZ-type theorems in Section 4, are obtained by modifying the original argument in [49] to fit the various nonsmooth settings. The crucial tool is the solvability of the Dirichlet problem with continuous data on the boundary of the metric ball for various notions of Laplacians on singular spaces, and lower estimates for the Laplacian of distance functions. An inspection of the proof of Theorem 1.3 shows that the restriction to step 2 arises from regularity issues of the boundary of the unit ball B G with respect to the Koranyi gauge: in step s > 2 the solvability of the Dirichlet problem with continuous boundary data in B G is an open problem.
Acknowledgements The authors would like to thank Ben Warhurst for a discussion concerning Theorem 1.2.

PRELIMINARIES
2.1. Notation and conventions. Throughout this article, a metric measure space X = (X, d, µ) is a separable metric space (X, d) equipped with a Radon measure µ that is finite and non-trivial on balls, i.e. 0 < µ(B) < ∞ for all balls B ⊂ X.
Given a subset F ⊂ X of a metric space and r > 0, we denote the open and closed r-neighbourhood of F (note that N r (F ) need not be the closure of N r (F ) unless X is a length space). For x ∈ X, we denote by B r (x) := N r ({x}) and B r (x) := N r ({x}), respectively, an open and closed ball centered at x with radius r. The Lipschitz constant of a map f : and the pointwise Lipschitz slope is If x ∈ X is an isolated point, then Lip f (x) = 0. Given u ∈ L 1 loc (X) and a µ-measurable set A ⊂ X with µ(A) > 0, denotes the average of u over A.

Doubling measures and averaging operators.
A measure µ on a separable metric space X is called locally doubling if, for every compact K ⊂ X, there exists r K > 0 and a constant C K > 0, such that for every x ∈ K and 0 < r ≤ r K . Let (X, d, µ) be a metric measure space and r > 0. Given a locally integrable function u ∈ L 1 loc (X), we denote by udµ, x ∈ X, the r-average function of u. Note that A r u(x) = u Br(x) . We will use the two notations interchangeably, depending on whether we want to view the average as a number, or an operator on a function space. Indeed, the function A r u : X → R is measurable, and A r defines a bounded linear operator A r : Moreover, in this case the operator norm satisfies A r L 1 →L 1 = a r L ∞ , see [4,Theorem 3.3]. This is true in particular when µ is a doubling measure. On the other hand, it is true that by the Lebesgue differentiation theorem if µ is infinitesimally doubling, cf. [31,Remark 3.4.29]. If X is doubling as a metric space, then there exists C > 0 so that A r L p →L p ≤ C for every r > 0 and every 1 ≤ p < ∞, cf. [4,Theorem 3.5]. However, A r is not a self-adjoint operator; the formal adjoint A * r of A r is given by for u ∈ L 1 loc (X). Indeed, a direct computation using the Fubini theorem yields that where 1/p + 1/q = 1.
We may express the r-laplacian using the averaging operator as the formal adjoint of the r-laplacian. Note that if A r : Remark 2.1. Identity (5) together with (6) imply the following Green-type formula for u ∈ L p (X), v ∈ L q (X) where 1/p + 1/q = 1: The asymptotic behaviour of ∆ r − ∆ * r as r → 0 plays an important role in Section 3.

Remark 2.2.
While most of our results will be formulated for metric measure spaces, they encompass the case of an open set Ω ⊂ X in the introduction. Indeed, an open subset Ω ⊂ X of a metric measure space can be regarded as a metric measure space Ω = (Ω, d| Ω , µ| Ω ). In particular, if X is locally doubling, then Ω is locally doubling.

Cheeger and Korevaar-Schoen energy, and Dirichlet forms.
In this section we recall some basic notions of calculus in metric measure spaces, needed in the discussion of the harmonic functions according to Cheeger, see Theorems 1.1 and 1.2. Moreover, we discuss the Korevar-Schoen energy and its relation to the Dirichlet forms.
One of the founding stones of the calculus theory in metric measure spaces are the notions of p-weak upper gradient and the minimal p-weak upper gradient, see [31,Chapter 5]. Related is a counterpart of the Sobolev space on a domain Ω ⊂ X, the so-called Newtonian space N 1,p (X) and the local Newtonian space N 1,p loc (Ω), see [31,Chapter 7] and [10] for precise definitions of these spaces and their further properties.
Cheeger's definition guarantees that Ch p is lower semicontinuous with respect to L p -convergence.

Definition 2.1.
Let Ω ⊂ X be a domain in a metric measure space. A continuous function u ∈ N 1,2 loc (Ω) is called harmonic if for any ϕ ∈ LIP c (Ω) there holds We assume continuity of Cheeger harmonic functions in the definition. However, in path-connected doubling metric measure spaces supporting Poincaré inequality, local Hölder (Lipschitz under stronger assumptions) continuity of Cheeger harmonic functions is established, see [36, Theorem 1.1].
Korevaar-Schoen energy. In [35], the authors defined a Sobolev space on a metric measure space X by considering what is now called as the Korevaar-Schoen energy for u ∈ L 2 loc (X): where the approximate energy densities e 2 r (u)dµ are given by For each r > 0, the approximate energy densities arise from a symmetric bi-linear form e r : which satisfies e r (u, u) = e 2 r (u). Furthermore, we define Dirichlet forms. We follow the notation of [12, Part I] and say that a Dirichlet form on L 2 (X) is a lower semicontinuous quadratic functional E : F → [0, ∞] defined on a dense domain F ⊂ L 2 (X), which is a Hilbert space equipped with the norm We denote by E : F × F → R the bi-linear form associated to E and assume the following strong locality property: if f, g ∈ F and g is constant on a neighbourhood of spt f , then E(f, g) = 0.
The form E admits a generator: a self-adjoint operator We call ∆ E the laplacian associated to the Dirichlet form. Moreover, if E is regular, it admits carré du champ, i.e. a positive, symmetric and continuous bi-linear operator Γ : F × F → L 1 (X) with the property that We may localize a given Dirichlet form E on L 2 (X) to a domain Ω ⊂ X. We say that f ∈ F loc (Ω) if there is a sequence (f k ) ⊂ F and open sets Ω k ⊂ Ω such that Ω k ⊂ spt(f k ), Ω = k Ω k , and f k = f on Ω k . By the strong locality, we may define the Carré du Champ on Ω k , where k Ω k = Ω and f = f k and g = g k on Ω k . We can correspondingly localize ∆ E , and we denote by D(∆ E ; Ω) the domain of ∆ E localized to Ω, cf. [12,Section 7].
The Cheeger and Korevaar-Schoen energies do not define Dirichlet forms in the sense defined above in general. Indeed, while the Cheeger energy is lower semicontinuous, it need not be a quadratic form. The Korevaar-Schoen energy is neither quadratic nor lower semicontinuous in general.
However, in RCD spaces and Carnot groups the Cheeger energy is quadratic and therefore constitutes a Dirichlet form. In particular in these classes the Korevaar-Schoen energy is known to be a constant multiple of the Cheeger energy and thus a posteriori defines a Dirichlet form, see [28]. We refer to Section 3 for a more precise discussion.
Spaces with a Korevaar-Schoen Dirichlet form. The Korevaar-Schoen energy has a distinguished role in this paper due to its close connection with the symmetrized r-laplacian, see Section 3.1. For the sake of a unified treatment of weakly amv-harmonic functions in Section 3, we introduce the following terminology for spaces where the Korevaar-Schoen energy defines a Dirichlet form. Definition 2.2. Let X = (X, d, µ) be a locally compact metric measure space. We say that X has a Korevaar-Schoen Dirichlet form (abbreviated KS-Dirichlet form), if E KS is lower semicontinuous with respect to L 2 -convergence and, for all u ∈ KS 1,2 (X), there exists an integrable function e 2 (u) satisfying The function e 2 (u) is called the energy density of u.
By the previous discussion, elaborated in Section 3, RCD spaces and Carnot groups equipped with certain homogeneous norms have a KS-Dirichlet form. Sturm [52,53] introduced the class of M CP (K, N )-spaces, which includes RCD-spaces and showed that an energy closely related to E 2 KS is a Dirichlet form on M CP (K, N ). Note that here we do not consider the existence and lower semicontinuity of the limit in the Korevaar-Schoen energy for maps with metric space targets. This fact is known for RCD and MCP-spaces while for Carnot groups it is -to our knowledge -an open question.
If X has a KS-Dirichlet form, then the domain of E 2 KS is KS 1,2 (X) and carré du champ is given by where the convergence is in the weak sense of measures. We denote by E KS and ∆ KS the bi-linear form, and its generator, associated to E 2 KS .

SYMMETRIZED r-LAPLACIAN AND AMV-HARMONIC FUNCTIONS
In this section we prove Theorems 1.1 and 1.2. We begin by relating the Korevaar-Schoen energy to the symmetrized r-laplacians, thus establishing the general principle that local minimizers correspond to functions with vanishing symmetrized r-laplacian, cf. Proposition 3.2.
In the case of Carnot groups the symmetrized r-laplacian agrees with the r-laplacian. Thus Theorem 1.2 follows directly from Proposition 3.2 and a second order Taylor expansion, once we establish that Carnot groups have a Korevaar-Schoen Dirichlet form, see Proposition 3.3.
Finally we consider non-collapsed RCD-spaces. After reviewing their relevant basic properties we show that having vanishing mm-boundary implies that ∆ r and ∆ r agree asymptotically (Lemma 3.8).
Since by [28] the Korevaar-Schoen energy is a Dirichlet form, Theorem 1.1 follows.
3.1. Symmetrized r-laplacian and the Korevaar-Schoen energy. Throughout this subsection, X = (X, d, µ) denotes a generic metric measure space. To define the symmetrized r-laplacian, we denote by k r : X × X → R the symmetric mean value kernel: Definition 3.1. Let X = (X, d, µ) be a metric measure space, and r > 0. The symmetrized r-laplacian of a function u ∈ L 1 loc (X) is given by Let u ∈ L 1 loc (X) and r > 0. A direct calculation yields the following identity (cf. formula (6)).
Moreover, we have the following elementary identities (recall the definitions of e r and E r in Section 2.3).
Lemma 3.1. Let u, v ∈ L 2 loc (X) and r > 0. Then ∆ r (uv) = u∆ r v + 2e r (u, v) + v∆ r u µ-almost everywhere. In particular, it holds that and whenever u, v ∈ L 2 (X). Here Proof. For µ-a.e. x ∈ X, we have proving the first identity in the claim. Together with Remark 2.1, the identity (11) implies (a): indeed, we have that The identity (b) follows directly from the observation that The characterization of local minimizers of the Korevaar-Schoen energy in terms of the symmetrized r-laplacians -in spaces with a KS-Dirichlet form -is a straightforward consequence of the preceding observations. Proof. Recall that local minimizers of E 2 KS in Ω are characterized by the property that E KS (ϕ, u) = 0 for every ϕ ∈ KS 1,2 (Ω) with compact support in Ω. Recall also, that the carré du champ for E KS is e(u, v). By Lemma 3.1 and the existence of the Korevaar-Schoen energy, this is equivalent with lim r→0 X ϕ ∆ r udµ = 0.

Carnot groups.
A Carnot group G is a connected and simply connected nilpotent Lie group whose Lie algebra g is finite dimensional and admits a stratification of step s ∈ N and m = h s . By using the exponential map, every Carnot group G of step s is isomorphic as a Lie group to (R m , ·), with the group operation given by the Baker-Campbell-Hausdorff formula.
Moreover, the stratification gives rise to a family of dilations δ t : G → G for t > 0, given by A continuous function ρ : G → [0, ∞) is said to be a pseudonorm on G, if ρ satisfies the following conditions: (1) ρ(δ r (x)) = rρ(x) for every r > 0 and x ∈ G, for every x ∈ G. A pseudonorm defines a pseudodistance on a group. Examples of such pseudodistances include the Korányi-Reimann (gauge) distance and the following : where x (i) = (x h i−1 +1 , . . . , x h i ) for i = 1, . . . , s, and x (j) stands for Euclidean norm in R v j . We observe that with the above notation, the coordinates in G can be given a form x = (x (1) , x (2) , . . . , x (s) ).
The number Q := s i=1 i dim(V i ) is called the homogeneous dimension of G. The Lebesgue measure L m =: µ is the Haar measure on G. Moreover, for all t > 0 it holds µ(δ t (E)) = t Q µ(E), for a measurable E ⊂ G. In particular, µ(B(x, r)) = cr Q for any left-invariant metric (c depends on the choice of metric) and the Hausdorff-dimension of G equals Q.
Next we come to Sobolev spaces. Let Ω ⊂ G be a domain and X = {X 1 , . . . , X m } a Jacobian basis of g, see [11,Chapter 1.3]. We define the horizontal Sobolev space HW 1,p (Ω) as space of measurable functions u ∈ L p (Ω) whose distributional horizontal derivatives X j u are L p -functions. We denote the horizontal gradient of u by The norm u p 1,p = u p L p (Ω) + |∇ H u| p L p (Ω) makes HW 1,p (Ω) into a Banach space. We also define the horizontal Laplacian ∆ G u of a C 2 -function u : Ω ⊂ G → R by The various different notions on Sobolev spaces on Ω ⊂ G agree, since G supports a 1-Poincaré inequality, cf. [29]. We mention here that, for u ∈ HW 1,2 (Ω), the function |∇ H u| is a minimal weak upper gradient of u when G is equipped with the Carnot-Carathéodory metric.
Below we will also use the approximate Pansu differential Du x : G → R of a Sobolev function u ∈ HW 1,2 loc (Ω) at almost every x ∈ Ω, which can be seen as a natural derivative in Carnot groups and plays an important role in function and mapping theory (e.g. in quasiconformal analysis). The following property relating the Pansu differential and the horizontal gradient of a function can be found in [45,Remark 3.3], see also [19,Corollary 3.5].

Proposition 3.3.
Let Ω ⊂ G be a domain a Carnot group and ρ be a pseudonorm on G of the form where z ′ = (z (2) , . . . , z (s) ) and F : R × R m−v 1 → R is C 1 outside the origin with ∂ 1 F > 0. Then for every u ∈ HW 1,2 (Ω) we have e 2 (u) = c|∇ H u| 2 almost everywhere, where c is a constant depending only on ρ, the data of G, and e(u) is the KS-energy density of u, cf. Definition 2.2. (14) is given as follows. Consider a norm N on R s with ∂ 1 N > 0, and norms · i on V i for i = 2, . . . , s, which need not necessarily be the Euclidean ones. The function

Remark 3.2. A particular class of homogeneous norms satisfying
is a homogeneous norm. This class includes the Korányi-Reimann gauge (12).

Proof. By the proof of [54, Theorem 4.1], we have that the KS-energy density exists and satisfies
for a.e x ∈ Ω. The last equality follows from the remark on the Pansu differential above. Here B ρ is the unit ball with respect to the pseudonorm ρ.
Since every two pseudonorms are equivalent (see [ This observation together with the assumptions on ρ imply that, for each z ′ ∈ M , where M ⊂ R m−v 1 is a bounded (and closed) set, the function F (·, z ′ ) is invertible and thus there exists a continuous function G : M → R such that for each fixed z ′ Using this and the change of variables, we compute where c is a constant. Indeed, since G is continuous and M is bounded, the integral of G 2+k is finite. Then, by [28,Lemma 4.18] we have that This completes the proof.
We finish the subsection by providing the proof of Theorem 1.2.
Proof of Theorem 1.2. By Proposition 3.3 we have that the energy density e 2 (u) exists and satisfies for every u ∈ HW 1,2 loc (Ω). Since |∇ H u| is the minimal weak upper gradient of u for the metric measure space (Ω, H Q , d cc ) (cf. [29, Proposition 11.6 and Theorem 11.7] for smooth functions from which the claim follows by approximation) it follows that u → E 2 KS (u) = Ω e 2 (u)dH Q is lower semicontinuous and thus the KS-energy is a Dirichlet form in the sense of Definition 2.2. Moreover by [11, Part I, Section 1.5 (A5)] ∆ G is the self-adjoint operator associated to the quadratic form u → Ω |∇ H u| 2 dz. This and (16) imply that ∆ KS is a constant multiple of the sublaplacian ∆ G .
Since ∆ r u = ∆ r u for small enough r in any domain compactly contained in Ω (by the translation invariance of H Q ), Proposition 3.2 implies that u ∈ HW 1,2 loc (Ω) is weakly amv-harmonic if and only if ∆ G u = 0, showing that (ii)⇔(iii). Since (i)⇒ (ii) trivially, it remains to prove (iii)⇒(i).
If ∆ G u = 0 then u ∈ C ∞ (Ω). Now we may argue exactly as in the proof of [26, Lemma 4.1] to obtain u(x) = u Br(x) + Cr 2 ∆ G u + o(r 2 ) (17) using the second order Taylor expansion [26, Lemma 2.7]: the first order terms cancel out since balls B r (x) with respect to ρ are symmetric, and the computation(s) in the proof of [26, Lemma 4.1] remain valid with the sets S(ε, x (2) , . . . , x (s) ) (in the notation of [26]) replaced by the sets { z (1) ≤ G(z ′ )} from the proof of Proposition 3.3. We obtain the constant in (17) where B ρ stands for the unit ball of ρ.
Finally we note that the error term o(r 2 ) in (17) is locally uniform in x (since the same is true for the second order Taylor expansion of a smooth function). This implies that ∆ r u(x) → 0 as r → 0 locally uniformly in x ∈ Ω, establishing (iii)⇒(i).
3.3. RCD-spaces. RCD-spaces arise from the study of synthetic curvature bounds, initially stated in terms of convexity properties of entropy functionals in Wasserstein spaces for Riemannian manifolds, cf. [58,52]. A Riemannian theory of the resulting curvature dimension conditions (CD(K, N )-conditions) is obtained by requiring in addition the quadratic property for the Cheeger energy, see Section 2.3. The notion of RCD-spaces was introduced and substantially developed by Ambrosio-Gigli-Savare, cf. [8]. In terms of the Γ-calculus introduced in Section 2.3 the RCD(K, N )-condition with parameters K ∈ R and N ≥ 1 on a metric measure space (X, d, µ) can be characterized by the Bochner inequality for any f ∈ D(∆) with ∆f ∈ W 1,2 (X) and any g ∈ L ∞ (X) ∩ D(∆).
A cornerstone of RCD-theory is the Bishop-Gromov inequality, due to Sturm [53]. To state it, define s K,N , v K,N : [0, ∞) → R (for K ∈ R and N ≥ 1) by where ω N the the Lebesgue measure of the unit ball in R N .
Note that v K,N (r) = ω N r N + O(r N +2 ) (cf. [32]). However, while N plays the role of dimension upper bound it need not equal Hausdorff dimension of X.
Rectifiability properties of RCD-spaces were established in [44] and improved in [14], where the constancy of dimension of RCD-spaces was proved. Theorem 3.5 (Theorem 3.8 in [14]). Let (X, d, µ) be an RCD(K, N )-space, K ∈ R and N ≥ 1. There exists a natural number n ≤ N such that for µ-almost every x ∈ X, R n (with the Euclidean norm and Lebesgue measure) is the unique tangent cone of X at x.
The number n above is called the essential dimension of X, and X is n-rectifiable, implying that n equals the Hausdorff dimension. The essential dimension appears in a result of Gigli-Tyulenev [28] which states that the Korevaar-Schoen energy is a constant multiple of the Cheeger energy, and thus a Dirichlet form, on RCD-spaces. Proposition 3.6 (cf. Theorem 2.18 and Proposition 4.19 in [28]). Let X be an RCD(K, N )-space, for K ∈ R and N ≥ 1, and let n be the essential dimension of X. Then there is a constant c = c(K, n) so that, if Ω ⊂ X is a domain, then for any u ∈ KS 1,2 (Ω) we have that where |∇u| is the minimal weak upper gradient of u.
We remark that by Proposition 3.6 local minimizers of E 2 KS are exactly the local minimizers of Ch 2 . It is known, cf. [36], that local minimizers of Ch 2 on a domain Ω are locally Lipschitz continuous in Ω and are characterized by the equation ∆u = 0.
Next we concentrate on non-collapsed RCD(K, N ) spaces X = (X, d, H N ). (For an earlier, slightly different, notion of non-collapsed RCD-spaces, see [33].) As a Corollary of the results above, noncollapsed RCD(K, N )-spaces are N -rectifiable and locally Ahlfors N -regular. Every point x ∈ X has a tangent cone (which is an RCD(0, N )-space) and on a set of full measure the tangent cone is unique and isometrically isomorphic to R N . The singular set S, defined as the set of points which have a non-Euclidean tangent cone, can be characterized using the densities as S = {x ∈ X : θ N (x) < 1}. Note that by Theorem 3.4 the limit defining θ N exists and satisfies θ N ≤ 1 everywhere.
The boundary of a non-collapsed RCD-space can defined by using the top dimensional stratum of the singular set. [24] and Definition 4.2 in [32]). The boundary ∂X of an RCD(K, N )space X is the closure of the set S N −1 \ S N −2 = {x ∈ X : the half space R N + is a tangent cone at x}. It has recently been proven that ∂X is (N − 1)-rectifiable, cf. [15]. Moreover it was shown in [16] that if a neighbourhood of a point does not meet the boundary then the metric measure theoretic boundary near that point vanishes. [32]). Let X = (X, d, H N ) be N -rectifiable. We say that X has vanishing mm-boundary (in an open set U ⊂ X), if the signed Radon measures

Definition 3.3 (Definition 1.5 in
are uniformly bounded in the total variation norm (in U ), and converge weakly to zero (in U ) as r → 0.
The following slight reformulation of [16,Theorem 1.2] can be obtained by a suitable re-scaling of the metric. We next employ Theorem 3.7 to obtain an asymptotic connection between the r-laplacian and the symmetrized r-laplacian. For the statement recall that for every ϕ ∈ C c (X \ ∂X).
Proof. By a cut-off argument we may assume that both ϕ and u have compact support in Ω. By Lemma 3.1 we have that The claim now follows from Lemma 3.8.
Proof of Theorem 1.1. Let u : Ω → R be harmonic. In particular, u is locally Lipschitz. By Lemma 3.6 u is a local minimizer of E 2 KS in Ω. Theorem 3.2 implies that lim r→0 Ω ϕ ∆ r udµ = 0 for every u ∈ KS 1,2 (Ω) with compact support in Ω. This and Corollary 3.9 imply that lim r→0 X ϕ∆ r udµ = 0 for every ϕ ∈ LIP c (Ω).

Remark 3.3.
It follows from Corollary 3.9 and the argument above that, if u : Ω → R is a locally Lipschitz continuous weakly amv-harmonic function, then u is Cheeger harmonic.

IMPROVING REGULARITY: THE BPZ-THEOREM FOR AMV-HARMONIC FUNCTIONS
In this section we discuss Blaschke-Privaloff-Zaremba type theorems (called BPZ theorems for the sake of brevity). In its classical version, see e.g. [38, Theorem 2.1.5], the BPZ theorem asserts that given an open set in R n a continuous pointwise amv-harmonic function solves locally the Laplace equation. Thus, the pointwise nullity of the amv-harmonic operator lim r→0 + ∆ r improves the regularity of amv-harmonic functions. Below we show that this is also the case of amv-harmonic functions in the homogeneous Carnot groups of step 2 and on Riemannian manifolds. Moreover, we show a counterpart of the BPZ theorem in Alexandrov spaces and discuss the relationship between amvharmonic functions defined in the weak sense and the Laplacian in Alexandrov spaces.

Carnot groups.
Solutions to the subelliptic Laplace equation in homogeneous Carnot groups are C ∞ , see e.g. the discussion in [11, Section 5.10] and the references therein. Using the notation from Section 3.2 we define B(p 0 , r), an open Korányi-Reimann ball (a gauge ball) centered at point p 0 ∈ G with radius r > 0 as follows Since, by [20,Corollary 1], the gauge balls are the X-NTA domains in step 2 Carnot groups G, they are by the Wiener criterion regular domains and the Dirichlet problem for the continuous boundary data has the classical smooth solutions for the harmonic sub-elliptic equation in B ⊂ G, see the [23, Theorem 3.9] and the discussion on pg. 409 in [20]. However, it is shown in [30], that for H n for n ≥ 2 balls in the Carnot-Carathéodory distance are not regular at the characteristic points (cf. [10,Example 14.4]). Therefore, due to the approach we take in the proof of Theorem 4.1 below, we will restrict our discussion to the case of step 2 Carnot groups and balls with respect to the homogeneous distance d G as in (12). We also remark, that whether gauge balls are X-NTA in general, step s, Carnot groups remains an open problem, see Conjecture 1 on pg. 429 in [20].
For the next theorem recall the distance d G associated to th Koranyi gauge (12) and the horizontal Laplacian ∆ G u (13) from Section 3.2.
Theorem 4.1 (The BPZ theorem for amv-harmonic functions in Carnot groups of step 2). Let Ω ⊂ G be a domain in the Carnot group G of step 2 equipped with the gauge distance d G and let f : Ω → R be a continuous function satisfying lim Then f is a sub-elliptic harmonic function, i.e. for all x ∈ Ω it holds that where C depends on v 1 (the dimension of the horizontal layer V 1 ) and the geometry of G.
The constant C is computed in [26,Lemma 4.1], see also the proof of Theorem 1.2. We further remark that the same lemma shows the Blaschke-Privaloff-Zaremba theorem for C 2 -functions in Carnot groups of step s.
Proof. In the proof we follow the original idea of Privaloff developed for the setting of R n , see [49,Theorem II]. Let p 0 ∈ Ω and B = B(p 0 , R) be a gauge ball centered at p 0 with radius R > 0 such that B ⊂ Ω. The above discussion allows us to infer that the sub-elliptic Dirichlet problem on B with the boundary data f has the (unique) solution, denoted by u, such that u ∈ C(B) and u = f | ∂B . Set φ = f − u. Then φ ∈ C(B) and φ| ∂B ≡ 0. The assertion will be proven if we show that φ ≡ 0 in B. We argue by contradiction. Namely, suppose that there exists q ∈ B such that φ(q) = 0 and without the loss of generality we assume that φ(q) < 0. Let us define the following function on B where (p −1 p 0 ) (1) stands for the Euclidean length in R v 1 of the horizontal part of the point p −1 p 0 ∈ Ω (see (12)). Claim: it holds that F ∈ C(B), F | ∂B ≥ 0 and that F (q) < 0. Indeed, if p ∈ ∂B(p 0 , R), then (p −1 p 0 ) (1) 4 + (p −1 p 0 ) (2) 2 = R 4 for (p −1 p 0 ) (2) ∈ R m−v 1 and so Thus F (p) ≥ 0. Similarly, we find that Hence, by the claim, there is q 1 ∈ B such that F (q 1 ) = min B F (in fact, q 1 ∈ B). Moreover, ∆ r F (q 1 ) ≥ 0 for all r ≤ d(q 1 , H 1 \ B) and by direct computations we verify that Therefore, upon applying the definition of ∆ r , we arrive at the following estimate Then, it holds Recall that, by the classical proposition on the Jacobian basis {X 1 , . . . , X m } (see [11,Corollary 1.3.9]), the subelliptic Laplacian at any point x ∈ G can be written as follows: where a (j) i are the polynomial coefficients. Therefore, in our case the sub-Laplacian reduces to the Laplacian in R v 1 . Namely, we get that Since f is assumed to be pointwise amv-harmonic and u ∈ C 2 (B), the definition of amv-harmonic functions together with [26,Lemma 4.1] imply that upon r → 0 + it holds that ∆ r f (q 1 ) → 0 and ∆ r u(q 1 ) → 0. Moreover, by applying [26,Lemma 4.1] again and by (21) we obtain that since above we assume that φ(q) < 0. In a consequence, we get that 0 ≤ Cv 1 R 2 φ(q) < 0 contradicting our assumption. The proof of the theorem is completed.

Riemannian manifolds.
For the readers convenience, we sketch the proof of the BPZ-theorem in the Riemannian setting. The argument is analogous to that in the proof of Theorem 4.1 Proof. We follow steps of the proof for Theorem 4.1 and argue by contradiction. Fix z 0 ∈ Ω and a ball B = B(z 0 , R) ⊂ Ω for a small enough radius R > 0, which will be determined below. As in the proof of Theorem 4.1, we denote by u the unique solution to the Laplace-Beltrami Dirichlet problem ∆ M u = 0 in B, such that u ∈ C(B) and u = f | ∂B . The existence of such u follows, for instance from the general theory for metric measure spaces, since the Riemannian manifold M with the Ricci curvature bounded from below satisfies the doubling condition and the 1-Poincaré inequality (see [10,Chapter A.7] and [29] for details).
Let us set φ := f − u and notice that φ ∈ C(B) and φ| ∂B ≡ 0. In order to prove the assertion we show that φ ≡ 0 in B. We argue by contradiction and suppose that there exists y ∈ B such that φ(y) = 0 and without the loss of generality we assume that φ(y) < 0.
Similarly as in the proof of Theorem 4.1 we define function F as follows where dist M (·, z 0 ) is a distance function to z 0 ∈ Ω defined in metric g on M . The direct computations for the Beltrami-Laplace operator, see the proof of [37,Theorem 4.1], allow us to conclude that where H z 0 (z) stands for the mean curvature of ∂B(z 0 , r) at z for r = dist M (z, z 0 ). According to the discussion in [47], it holds that H z 0 (z) = n−1 r − 1 3 Ric z 0 (v, v)r + O(r 2 ), asymptotically as r → 0, where v ∈ S n−1 ⊂ T z 0 M (notice that this formula arises as a perturbation of the mean-curvature formula for a sphere in the zero-curvature case). Therefore, by taking R small enough, and hence small enough As in the proof of the previous theorem, we find y m ∈ B such that F (y m ) = min B F (in fact, y m ∈ B). In order to get the contradiction we need to conclude that ∆ r f (y m ) → 0 and ∆ r u(y m ) → 0 as r → 0 + , cf. (20) and the discussion following it. The first convergence follows, as f is amv-harmonic, while the second follows from Claim (1), Chapter 6.32 in [59]. Indeed, by [59] solutions of the Laplace-Beltrami equation are smooth (assumptions of [59, Chapter 6.32] require M to be compact, but we need this result locally). Similarly to the Euclidean argument, based on expanding u in the Taylor expansion up to the second order, one can show that lim r→0 Thus, by analogous reasoning as in the previous theorem we obtain that ≤ cφ(y) 4(n + 2)R 2 < 0.
This contradicts our assumption that φ(y) < 0. Hence, the proof of the theorem is completed.

Alexandrov spaces and surfaces.
The setting of Alexandrov surfaces. The above discussion for Riemannian manifolds can be adapted to the setting of possibly nonsmooth surfaces, namely for closed Alexandrov surfaces with bounded integral curvature, see e.g. [55,56,50,5]. In particular, this together with recent results by Lytchak-Stadler [39, Theorem 1.1] and Petrunin [48] imply also a BPZ theorem for non-collapsed RCD(K, 2) spaces.
An Alexandrov surface (S, d) is a topological 2-manifold with a length metric d inducing the manifold topology which has bounded integral curvature in the sense of [57, Definition 2.1]. Alexandrov surfaces admit a notion of curvature which is a Radon measure, denoted by ω. For this and more details we refer the interested reader to [56,57]. Here we mention that there exist several equivalent definitions based on polyhedral approximations or on the Gauss-Bonnet theorem. Moreover, let us recall that examples of Alexandrov surfaces encompass, for instance, compact polyhedral surfaces (both Euclidean, hyperbolic, spherical or more generally Riemannian), locally CAT (0) surfaces, compact metric surfaces with lower bounded curvature in the sense of Alexandrov (in particular compact Gromov-Hausdorff limits of surfaces with lower bounded Gaussian curvature, see [17,39]), and also surfaces with conical singularities, see [56,57].
Let (S, d) be a closed Alexandrov surface (i.e. compact surface without boundary) equipped with a metric d and the curvature measure ω. We say that (S, d) is without cusps, if ω({x}) < 2π for all x ∈ S. In particular, both the BIC surfaces with non-positive curvature measure and all the smooth surfaces are automatically without cusps.
Following [55,56], if (S, h) is a closed Riemannian surface, then we denote by V (S, h) the space of functions v : S → R such that µ = ∆ h v is a (signed) measure. By using v one may define a new distance, denoted by d h,v , see [56,Proposition 7.1]. Then, the following result relates Alexandrov's surfaces to smooth geometry, see [56,Theorem 7.3]. The proof is similar to the one for Proposition 4.2 and thus we will comment on the key difference only. As in the Riemannian case we argue by contradiction and localize the discussion on a ball B(z 0 , R) ⊂ Ω. Then, one considers functions φ and F similar as in the proof of the above theorem. By the Reschetnyak-Huber theorem we reduce the discussion to the plane with conformally flat metric and by (22) get that ∆ S dist S (z, z 0 ) 2 = e −2u ∆ R 2 dist(z, z 0 ) 2 = 2e −2u > 0. From this point, the reasoning follows strictly the proof of Proposition 4.2.
The setting of Alexandrov spaces. Another setting where a BPZ type theorem can be proven is the setting of Alexandrov spaces. We say that M is an Alexandrov space, if it is locally compact complete length metric space with the lower curvature bound defined via comparison of geodesic triangles in M to geodesic triangles in space form of the appropriate constant sectional curvature, see for instance [13,17,46]. The Hausdorff dimension of an Alexandrov space turns out to be always an integer and it defines the dimension of the given Alexandrov space. Moreover, a point in an n-dimensional Alexandrov space M is called regular, if its tangent cone is isometric to Euclidean n-space with standard metric, see e.g. Section 7.1 in [6]. A point which is not regular is called singular and it turns out that the Hausdorff dimension of the set of all singular points in M does not exceed n − 1. Furthermore, the set of all regular points is contained in a DC 0 Riemannian manifold, see Definitions 3.1 and 2.1 in [6].
Let Ω ⊂ M be a bounded domain in Alexandrov space M . For a semiconcave function u : Ω → R (see Preliminaries in [61]) a natural Laplacian L u defined via the Dirichlet form is a signed Radon measure and its nonsingular part in the Lebesgue decomposition of L u satisfies ∆u(x) = n Σx H x u(ξ, ξ)dξ, for almost all points x ∈ Ω, where H x denotes the Perelman hessian while Σ x stands for the space of directions at x, see [48] and [61,Preliminaries]. One of the key features of hessian H x is that for semiconcave functions in Ω (or more general DC functions, see [6]) one can prove the existence of the second order Taylor expansion based on H x , see [61,Formula (2.6)] and [6, Proposition 7.5].
The above observations allow us to state the following variant of the BPZ theorem in Alexandrov spaces. ∆u(x).
In particular, a semiconcave function u on Ω is harmonic (in the sense of the Laplacian ∆) if and only if lim r→0 + ∆ r u(x) = 0 for every regular point x of u.
The proof of this observation is the immediate consequence of the proof of [61, Lemma 2.3]. We remark that the theorem shows that for a large class of functions the amv-Laplace harmonic functions coincide with those which are harmonic in Alexandrov spaces with respect to Laplacian ∆. However, notice that semiconcave functions are locally Lipschitz. Therefore, in the setting of Alexandrov spaces a characterization of the amv-harmonic functions stated in Proposition 4.5 holds under stronger assumptions than just a continuity as it is in R n , Carnot groups and on manifolds (cf., respectively, [38, Theorem 2.1.5] and Theorem 4.1 and Proposition 4.2). Finally let us comment that, since Alexandrov spaces without boundary are conjectured to have vanishing mm-boundary (see [32,Section 1.3]) and this is known to be the case for two-dimensional Alexandrov spaces (Theorem 1.8 in [32]), the characterization part of Proposition 4.5 for an Alexandrov space M without the boundary is related to our Theorem 1.1.