Regularity of Solutions to the Quaternionic Monge–Ampère Equation

The regularity of solutions to the Dirichlet problem for the quaternionic Monge–Ampère equation is discussed. We prove that the solution to the Dirichlet problem is Hölder continuous under some conditions on the boundary values and the quaternionic Monge–Ampère density from Lp(Ω)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^p(\Omega )$$\end{document} for p>2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p>2$$\end{document}. As a step towards the proof, we provide a refined version of stability for the weak solutions to this equation.


Introduction
The theory of quaternionic plurisubharmonic functions in H n was initiated by Alesker in [1]. He points out there that actually the definition of quaternionic plurisubharmonic function was suggested by Henkin. In particular, the quaternionic Monge-Ampère operator was defined in [1] for continuous functions. In the case of smooth functions, this is the Moore determinant, cf. [17], of the quaternionic Hessian. This will be explained in details in Preliminaries. In [2], the Dirichlet problem for this operator was solved for the right-hand side continuous up to the boundary and a continuous boundary data. It involved the methods based on the approach of Bedford and Taylor [6,7]. Later, the theory was further generalized to the case of hypercomplex manifolds by Alesker Dedicated

Preliminaries
We treat H n = {(q 0 , . . . , q n−1 ) | q i ∈ H} as a right H module, where H = {x 0 + x 1 i + x 2 j + x 3 k | x i ∈ R} and i 2 = j 2 = k 2 = ijk = −1. We often use the following identifications We also treat H n as a hypercomplex manifold with the complex structures I , J , K induced by the right multiplication by i, j, k, respectively (see [20] for details).
A function f in a domain ⊂ H n is called quaternionic plurisubharmonic, qpsh for short, if it is upper-semi-continuous and subharmonic when restricted to any slice of by a right quaternionic line. The set of qpsh functions in will be denoted by QPSH( ). These functions were introduced and studied in [1]. For a smooth function f and the differential operators being qpsh is equivalent to the fact that the matrix ∂ 2 f ∂q k ∂q l k,l is non-negative which means that p k ∂q k ∂q l r l ≥ 0 for all p, r ∈ H n . More details and references for the quaternionic linear algebra may be found in [1].
It is well known, cf. [1], that there is no determinant function defined on the set of all quaternionic matrices having all the properties of determinants of real or complex matrices. This causes the problem since it is tempting to define the quaternionic Monge-Ampère operator of f as a determinant of ∂ 2 f ∂q k ∂q l k,l . One may for instance consider the so-called Deieudonné or Study determinants but these are, among other disadvantages, always non-negative. Luckily, the matrix we obtain is hyperhermitian in the sense For such matrices Moore, cf. [17], has defined the notion of a determinant, we denote it by det M , which reflects enough properties of the determinant function. The original definition of quaternionic Monge-Ampère operator for smooth f was just det M ∂ 2 f ∂q k ∂q l k,l , cf. [1]. It was extended to continuous f by a measure theoretic construction performed originally in the complex case in [6].
In [23] pluripotential theory was adopted to the quaternionic setting using methods from [6,7] and two formal operators d 0 , d 1 . As it was shown in [20] after a suitable choice in [23] they coincide with ∂, the (1, 0) part with respect to the complex structure I of an exterior differential d in H n , and the differential operator ∂ J defined by This operators were studied in the context of the quaternionic Monge-Ampère equation for example in [5]. Because the operators ∂, ∂ J are standard we use them in what follows.
Coming back to the content of [23], for a locally bounded u ∈ QPSH( ) a closed positive current ∂∂ J u was defined. This uses the notion of positivity, introduced by Alesker and Verbitsky in [4,5], for the forms α in H n of a type (2k, 0) with respect to I such that J (α) = α. For us it is sufficient to know that n−1 i=0 dz 2i ∧ dz 2i+1 (corresponding to dz i ∧ dz i from the complex setting) and n = dz 0 ∧ dz 1 ∧ · · · ∧ dz 2n−2 ∧ dz 2n−1 are positive forms. The quaternionic Monge-Ampère operator for locally bounded functions is then defined by (∂∂ J u) n which is a Borel measure in . This definition actually gives a multiple of the operator defined by Alesker since, as checked for example in [4], This also explains the convention to integrate (2n, 0) forms with respect to I in H n , namely for α ∈ 2n,0 The details concerning this approach can be found in [22][23][24].
For the rest of the paper, we denote by a fixed quaternionic strictly pseudoconvex domain i.e. a C 2 −smoothly bounded domain such that there exists ρ ∈ QPSH(U ) ∩ C 2 (U ), ⊂⊂ U , = {ρ < 0}, dρ = 0 on ∂ and (∂∂ J ρ) n ≥ n on U . We will also use the notion of quaternionic capacity, which we denote by cap, introduced in [24]. For a compact set, K it is given by see [24] for details. Let 0 ≤ α ≤ 1, we denote by Li p α ( ) the space of α-Hölder continues functions i.e. f ∈ Li p α ( ) if f ∈ C 0 ( ) and In the next two sections we perform many estimates. We use the notation C(·, •, . . .) for a constant depending only on the quantities ·, •, . . . . In particular C's depending on the same set of variables may vary from line to line.

Stability of Solutions to the Quaternionic Monge-Ampère Equation
In this section, we prove a version of the stability estimate for the quaternionic Monge-Ampère equation. In [20] it was proven that the uniform norm of the difference of solutions is under control of the uniform norm of the difference of boundary data and the L q norm of the difference of Monge-Ampère densities when they belong to L q for q > 2, cf. Proposition 4 in [20]. Our goal here is to prove that the uniform norm of the difference of solutions is under control of the L p norm of that difference for appropriate p, cf. Theorem 1 below.
For the purposes of this section we set the following convention. Whenever u ∈ QPSH( ) is locally bounded then writing (∂∂ J u) n p automatically implies that we assume the Borel measure (∂∂ J u) n has a density with respect to the Lebesgue measure L 4n and this density is in L p ( ). We will need the following comparison of quaternionic capacity and Lebesgue measure. Its proof uses first author's C 0 estimate for the complex Monge-Ampère equation which in turn uses the comparison of the complex capacity and the Lebesgue measure, cf. [15,16].
Remark In the given reference this lemma is shown for compact sets but then the inequality follows easily for any Borel set with the same constant.
Another technical result which we need is the lemma below. It was used implicitly back in [14], since then it turned out to be very useful while performing pluripotential estimates. Its form below is a combination of Lemma 1.5 in [11] and Lemma 2.4 in [9]. In [19] it was attributed to De Giorgi.

Lemma 2 (De Giorgi) Suppose that a non-increasing, right continuous function
for any s ≥ 0, t ∈ [0, 1] and some α, B > 0. There exists s ∞ , depending on α, B and s 0 such that f (s 0 ) ≤ 1 2B , satisfying f (s) = 0 for any s ≥ s ∞ . In fact one can choose s ∞ to be equal s 0 + 2B f (s 0 ) α 1−2 −α . If in addition ( * ) holds for all s, t ≥ 0 then one can take even s ∞ = 2B f (0) α 1−2 −α . Now, we show the main technical fact needed for the proof of Theorem 1. The reasoning we perform in the rest of the section is based on the one presented in [11].
For 0 < α < p−2 np there exists a constant C(c 0 , α, diam( )), depending on c 0 , α and the diameter of , such that for any > 0 1 n for s ≥ 0 and > 0. Firstly, note that for all t, s ≥ 0, > 0 and w ∈ QPSH( ) such that 0 ≤ w ≤ 1 one obtains from inclusions of sets, super-additivity and the comparison principle, cf. Theorem 1.2 in [24], From the Hölder inequality we obtain where q ∈ (1, 2) depends only on p which is the conjugate of p and we choose it so that q p = 1 + nα. Since 1 + nα < 2 p this is always possible. Taking the supremum over all w and n'th root of both sides gives tb (s + t) ≤ C(c 0 , α, diam( ))(b (s)) 1+nα for any s, t ≥ 0 and > 0. One easily checks, as in [19], that the function b satisfies all the assumptions of Lemma 2. This gives From the comparison of volume and capacity (see Lemma 1) it follows that v − u ≤ +C(c 0 , α, diam( ))b (0) nα almost everywhere in . Since u and v are subharmonic we obtain that this holds in , i.e.
We are ready to prove the announced stability of weak solutions.
For any r ≥ 1 and 0 < γ < r r +np + p np p−2 := γ r there exists a constant C(c 0 , γ, diam( )), depending only on c 0 , γ and the diameter of , such that Proof First of all, we may assume that (v − u) + L r ( ) = 0 because otherwise the inequality holds with any constant C(c 0 , γ, diam( )). Arguing as in the beginning of the proof of Proposition 1 we obtain for any > 0 Since on the set {u − v < − } the function v−u + r p is bigger than 1, due to Hölder's inequality we may further estimate Applying Proposition 1 we get Choosing α such that γ = r r +np + p α (which is always possible since when α varies in 0, p−2 np the quantity r r +np + p α varies in (0, γ r )) results in

Hölder Continuity of Weak Solutions
In this section, we proceed to proving regularity of solutions to the Dirichlet problem for the quaternionic Monge-Ampère equation under the conditions on the boundary data and the density. This is the content of Theorem 3. For that goal we first consider a more general situation in Theorem 2 below and then check that under the assumptions of Theorem 3 one can apply Theorem 2.

Theorem 2 Let
⊂ H n be a quaternionic strictly pseudoconvex domain. Suppose p > 2, f ∈ L p ( ) is a non-negative function, φ ∈ C(∂ ) and u is the solution to the Dirichlet problem such that u( ) is finite. If there exists 0 < ν < 1 and b ∈ Li p ν ( ) such that b ≤ u in and b = φ on ∂ then u ∈ Li p α ( ) for any 0 ≤ α < min{ν, 2γ 1 }.

Remark
Let us just emphasize that in the whole section we denote by u, for a locally bounded u, the distributional Laplacian. Since qpsh functions are in particular subharmonic for them u is a positive distribution thus a measure and so it makes sense to write u( ).
We introduce the notation needed for the proof of this theorem. This approach is similar to the one presented in [18]. For a fixed number δ > 0 and a subharmonic u we consider (B(0, 1))δ 4n q− p ≤δ u( p)dL 4n ( p), for q ∈ δ .
The lemma below is a composition of Lemmas 4.2 and 4.3 in [11]. As was noted by Nguyen in [18] proofs given originally in [11] for plurisubharmonic functions used only subharmonicity.

Lemma 3 Let u be a subharmonic function in .
For a fixed 0 < α < 1 the following are equivalent Moreover, there exists a constant c, depending only on the dimension n, such that for all δ > 0 sufficiently small

Lemma 4 Let u and b be as in Theorem 2.
There exists a constant c n such that for all 0 < δ and q ∈ ∂ δ u δ (q) ≤ u(q) + c n δ ν .
Proof Denote by h the harmonic extension of b |∂ to . By the Proposition 2.4 in [18] we know that h ∈ Li p ν ( ). Fix q ∈ ∂ δ and take p, p 0 ∈ H n such that p = p 0 = δ, u δ (q) = u( p + q) and q + p 0 ∈ ∂ . Since b ≤ u ≤ h in , with equalities on ∂ , we obtain the following string of inequalities Proof of Theorem 2 Fix 0 < γ < γ 1 and consider the functioñ u δ (q) = max{û δ (q), u(q) + c n δ ν }, for q ∈ δ u(q) + c n δ ν , for q ∈ \ δ , which by Lemma 4, is qpsh in , asû δ ≤ u δ (cf. Proposition 2.1(4) in [22]); and continuous in . Applying Theorem 1 for u + c n δ ν ,ũ δ and r = 1 we obtain From the construction ofũ δ it follows that the last inequality is equivalent to Estimating further, by Lemma 3 and the trivial estimate (û δ − u − c n δ ν ) + ≤û δ − u, we obtain for sufficiently small δ > 0 This results in for all δ sufficiently small. Since the constant C is independent of δ we obtain, again due to Lemma 3, that u ∈ Li p min{ν,2γ } ( ). Since γ was arbitrary in (0, γ 1 ) this gives our claim.

Theorem 3 Let
⊂ H n be a quaternionic strictly pseudoconvex domain. Suppose p > 2, f ∈ L p ( ) is a non-negative function bounded in a neighborhood of ∂ and φ ∈ C 1,1 (∂ ). Then the Dirichlet problem is solvable and the unique solution is in Li p α ( ) for any 0 ≤ α ≤ 2γ 1 .
Proof The continuous solution u exists and is unique as was shown in [20]. We need to check the assumptions of Theorem 2. For that goal, we construct a function b as in Theorem 2 having in addition the properties that it is subharmonic and the Laplacian of it has finite total mass. This will of course imply that the Laplacian of u has finite total mass since b and u, both subharmonic, will agree on ∂ and b ≤ u in . Take h to be the solution to the Dirichlet problem By the comparison principle, cf. Theorem 1.2 in [24], it is above any v ∈ QPSH( )∩ C( ) such that v |∂ = φ.
We will show that h is Lipschitz in and its Laplacian has finite total mass in . Suppose U is a neighborhood of such that φ is extendable to a functionφ ∈ C 1,1 (U ).
From the properties of h it means that for any q < ĥ (r + q) − C q ≤ h(r ) for all r ∈ .
For r ∈ and q ≤ such that r + q ∈ we thus obtain what results in h(r + q) − h(r ) ≤ C q for r ∈ and q < such that r + q ∈ .
This shows that h is locally Lipschitz continuous in and consequently Lipschitz continuous.
By the assumptions for some M > 0 we have f ≤ M away from a compact K ⊂ . Let u be the continuous solution to our Dirichlet problem. Take B big enough for the function Bρ + h to be below u in a neighborhood of K and such that B n > M. Then by supper-additivity (∂∂ J (Bρ + h)) n ≥ (∂∂ J Bρ) n ≥ f n at least in \ K . The comparison principle, cf. [24], implies that Bρ + h ≤ u in \ K since the inequality holds on the boundary of this set.
We define b = Bρ + h. It is Lipschitz continuous and its Laplacian has finite total mass in since h has these properties.