On Neumann problem for the degenerate Monge–Ampère type equations

In this paper, we study the global C1,1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$C^{1, 1}$\end{document} regularity for viscosity solution of the degenerate Monge–Ampère type equation det[D2u−A(x,Du)]=B(x,u,Du)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\det [D^{2}u-A(x, Du)]=B(x, u, Du)$\end{document} with the Neumann boundary value condition Dνu=φ(x)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$D_{\nu }u=\varphi (x)$\end{document}, where the matrix A is under the regular condition and some structure conditions, and the right-hand term B is nonnegative.


Introduction
In this paper, we study the degenerate Monge-Ampère type equation det D 2 u -A(·, Du) = B(·, u, Du) in , ( 1 ) associated with the Neumann boundary value condition where is a bounded domain, Du and D 2 u denote the gradient and Hessian matrix of second order derivatives of the unknown function u : → R, respectively, A : × R n → R n×n is a symmetric n × n matrix-valued function and A ∈ C 2 (¯ × R n , R n×n ), B :¯ × R × R n → R + ∪ {0} is a nonnegative scalar function, ν denotes the unit inner normal on ∂ and ϕ is a scalar-valued function defined on ∂ . We shall use x, z, and p to denote the points in , R, and R n , respectively. We introduce the definitions of nondegenerate and degenerate Monge-Ampère type equation respectively. The Monge-Ampère type equation (1) is called nondegenerate or degenerate if D 2 u -A is positive definite or nonnegative definite, namely the right-hand side term B ≥ B 0 > 0 or B ≥ 0, respectively, where B 0 is a constant. Accordingly, a solution u of the Neumann boundary value problem (1)-(2) is called elliptic (degenerate elliptic) if D 2 u -A > 0 (≥ 0).
We say that A is strictly regular in , if n i,j,k,l=1 A kl ij (x, p)ξ i ξ j η k η l ≥ c 0 |ξ | 2 |η| 2 (3) holds for all (x, p) ∈ × R n , ξ , η ∈ R n with ξ · η = 0, and some positive constant c 0 , where If c 0 on the right-hand side in (3) is replaced with 0, we say that A is regular in . As usual, the strictly regular condition and regular condition are also said to be the A3 condition and the A3w condition, respectively, see [1,2].
When A = 0, the nondegenerate Monge-Ampère type equation (1) arises in various aspects such as optimal mass transportation problems, geometric optics, and conformal geometry (see, for instance, [1,[9][10][11]). Ma, Trudinger, and Wang [1] obtained a priori interior second order estimates under the A3 condition. The boundary C 2,α estimate for the Dirichlet problem was established by Huang et al. [12] under the A3 condition. Jiang et al. [13] obtained the global smooth solutions under the A3w condition. Furthermore, Leoper [2] proved that the A3w condition is necessary for C 0 regularity of the optimal map in the optimal transportation problem.
For the Neumann boundary value problem, the global regularity of solutions for the classical Monge-Ampère equation was established by Lions, Trudinger, and Urbas in [14]. Subsequently, Jiang, Trudinger, and Xiang [15] proved the second order derivative estimates of solutions for the Monge-Ampère type equation (1) under the A3w condition. However, as far as we know, there is no regularity result for the degenerate Monge-Ampère type equation with Neumann boundary value condition, which is the motivation of this work.
Our global second derivative estimate is obtained by mixing the Bernstein type argument and the barrier argument. In order to construct the global barrier function under the regularity assumption of A, we need to assume the existence of a subsolution (strict subsolution) u to equation (1) satisfying together with the Neumann boundary condition In order to treat the degenerate case for equation (1), we assume that B satisfies the following conditions: and in × R × R n for i, j = 1, . . . , n, where C 1 , C 2 , and C 3 are positive constants, I denotes the unit matrix, B i , B z , and B p i denote the partial derivatives of B with respect to x i , z, and p i , respectively. In fact, we can have the following relaxed version of condition (6) if |B| ≤ C: in × R × R n for i = 1, . . . , n, where C 4 depends on n, sup B, and C 1 .
We define the domain to be uniformly A-convex with respect to for all vectors τ = τ (x) tangent to ∂ and some positive constant δ 0 . For the Neumann problem, the definition of uniform A-convexity of with respect to u is introduced in [15]. We refer the readers to [15] for more descriptions of this definition. Based on (10), we can define a barrier function in the neighborhood of the boundary ∂ by using the defining function of the domain , which will be used in the double normal derivative estimate of u on ∂ in Sect. 3. We now formulate the main result.
We remark that the alternative conditions (a) and (b) for the matrix A are used in the treatment of the third derivative terms in the proof, see the details in Sect. 3.
Since the global second order derivative estimates are independent of the lower bound of B, we have the following existence result by approximation. Before stating the theorem, we introduce the definition of viscosity solution of equation (1). An upper (lower)semicontinuous function u is called a viscosity subsolution (supersolution) of equation (1) if, for any function φ ∈ C 2 ( ) such that uφ has a local maximum (minimum) at some point x 0 ∈ , there holds A continuous function u is a viscosity solution of equation (1) if it is both a viscosity subsolution and a viscosity supersolution of equation (1). Based on the a priori second derivative estimate, we now establish the existence result in the following theorem.

Theorem 1.2
Under the assumption of Theorem 1.1, assume that "B is a nonnegative function" and "a strict subsolution u ∈ C 2 (¯ )" in place of "B is a positive function" and "a subsolution u ∈ C 2 (¯ )", respectively. Assume further that there exists a supersolutionū of equation (1) satisfying D νū ≤ ϕ(x) on ∂ , and A satisfies the quadratic structure condition from below, namely where μ 0 is a positive constant depending on sup |u| and I denotes the n × n identity matrix. Then there exists a degenerate elliptic viscosity solution u ∈ C 1,1 (¯ ) for the Neumann problem (1)-(2).
Remark 1 Note that in Theorem 1.2, if B is a positive function and u ∈ C 2 (¯ ) is merely a subsolution, we can get the existence and uniqueness of the elliptic solution u ∈ C 3,α (¯ ) of problem (1)-(2) for some α ∈ (0, 1). We can prove this by following the proof of Theorem 1.2 in [15] and replacingū in its original proof with u. The reason why we can replacē u with u is that both A in (1) and ϕ in (2) have no u dependence. Indeed, when A in (1) is independent of u, we can obtain Lemma 2.2, which can replace Lemma 2.1 of [15]. Meanwhile, since ϕ in (2) has no u dependence, we are still able to get the key inequality (47) in our proof. Then we can follow the proof of Theorem 1.2 in [15] to establish u ∈ C 3,α (¯ ).
Remark 2 It is clear that Theorem 1.2 can be applied to the Monge-Ampère equation with homogenous right-hand side as follows: where c 0 is a nonnegative constant, α ∈ [2(n -1), +∞) is a constant. Thus, there exists a unique convex solution u ∈ C 1,1 (¯ ) of equation (14) under the Neumann boundary value condition (2). For the Dirichlet problem of the totally degenerate Monge-Ampère equation (14) when c 0 = 0 and A ≡ 0, we refer the reader to [16][17][18].
Remark 3 Condition (a) in Theorems 1.1 and 1.2 can include A = diag(u 2 1 , . . . , u 2 n ) as an example, which satisfies both the A3w condition and the structure condition (13). Both conditions (a) and (b) in Theorems 1.1 and 1.2 can cover the special case when the matrix A is independent of p, so that the standard Monge-Ampère case is already involved. It would be interesting to relax restrictions (a) or (b) of the matrix A in Theorems 1.1 and 1.2.
The paper is organized as follows. In Sect. 2, we introduce two lemmas which will be used to deal with the degenerate term and in the construction of the global barrier function, respectively. In Sect. 3, we give the proof of Theorem 1.1 which constitutes the core of this paper. In this section, we reduce the global estimates to the boundary by using a global auxiliary function, where the computations of the third order derivatives are carried out in detail under the two alternative conditions of Theorem 1.1. For the boundary estimates, we construct a suitable barrier function which can be used to control the degenerate term, originating from the differentiation of equation (1). In this part, we mix the techniques of dealing with the Neumann problems for nondegenerate Monge-Ampère type equations in [14,15] and the Dirichlet problems for degenerate Monge-Ampère type equations in [3,[5][6][7]. In Sect. 4, we complete the proof of Theorem 1.2.

Preliminaries
In this section, we introduce two lemmas. One lemma is used to deal with the degenerate term which arises in differentiating equation (1). The other lemma is used to construct the global barrier function. In equation (1), we first suppose B > 0 in , Then both the matrices {ũ ij } and {ũ ij } are positive definite. We can rewrite equation (1) in the form whereB := log B. By differentiating equation (15) in the direction ξ ∈ R n once and twice respectively, we havẽ where and where . Note that we use the standard summation convention in the context that repeated indices indicate summation from 1 to n unless otherwise specified. Next, we point out several facts that will be used in the proof of Theorem 1.1.
Assume that B satisfies (6), B > 0, andB = log B. Then we have the following properties: holds for i = 1, . . . , n, where the constant C depends on n, C 1 , A, and u C 1 (¯ ) . (ii) If B also satisfies conditions (7)-(8), then holds for i = 1, . . . , n, where the constant C depends on n, C 1 , C 2 , A, and u C 1 (¯ ) , and the constant C depends on C 3 and A.
Proof of Lemma 2.1 By choosing ξ = e i in (18) and usingũ ij = u ij -A ij and (9) (the relaxed version of (6)), we get (20) and finish the proof of conclusion (i).
Next, we turn to proving (ii). By condition (8), we have where δ kl denotes the usual Kronecker delta, the constant C depends on C 3 and A. It follows from (6) and (7) that the sum of the other terms on the right-hand side of (7) has a lower bound -C[1 + max j (|ũ ij |)]B -1 n-1 , where the constant C depends on n, C 1 , C 2 , sup |B|, A, and sup |Du|.
Remark 4 We remark thatB = log B satisfies condition (8) if it is semi-convex in p. The term n k=1B p k D ii u k on the right-hand side of (21) can also be dealt with directly. We shall explain the treatment in the later discussion.
We introduce the following linearized operator of equation (1): We introduce a fundamental lemma, which is important for constructing the global barrier function in the next section.

Lemma 2.2
Let u ∈ C 2 (¯ ) be an elliptic solution of (1), and let u ∈ C 2 (¯ ) be an elliptic holds in for sufficiently large positive constant K and uniform positive constants 1 , C depending on A C 2 , B C 1,1 , , u C 1 , u C 1 , δ, and C 1 in (6), where L is the operator in (23), and T = n iũ ii .
Proof Since A is regular, similar to Lemma 2.1(ii) in [20], we have where D 2 u -A(x, Du) ≥ δI > 0 and condition (6) are used in the second inequality of (25), M 1 = sup x∈ B(x, u, Du), and the constant C changes from line to line. Then estimate (24) follows by further adjusting the constant C.
Remark 5 Note that the barrier inequality in Lemma 2.2 is different from that in [15].
Here the constants 1 and C in estimate (24) of Lemma 2.2 are independent of the positive lower bound of B. Note also that estimate (24) is proved under condition (6), while in [19] a similar estimate is proved under condition (8).

Global second order derivative estimates
In this section, by constructing a global auxiliary function, we can reduce the global estimates to the boundary. On the boundary, in order to get the mixed tangential normal second order derivatives of u, we apply the tangential operator to the boundary value condition (2). A bound of double normal derivative of u on ∂ is obtained by constructing a suitable auxiliary function. We use the key trick of [14] to establish the global second order derivative estimates of solutions for the degenerate Monge-Ampère type equations with the Neumann boundary value condition. For the argument below, we assume that the functions ϕ, ν are smoothly extended to¯ . The constant C in this section changes from line to line.
Mixed tangential normal derivative estimate on ∂ . We introduce the tangential gradient operator δ = (δ 1 , δ 2 , . . . , δ n ), where δ i = n j=1 (δ ijν i ν j )D j for i = 1, . . . , n. Applying this tangential operator to boundary condition (2), we have If τ is a direction tangential to ∂ at any point y ∈ ∂ , we have where the second equality and the last equality are both valid using the fact that τ · ν = 0, and the fifth equality holds by (26). Hence, from (27) we obtain on ∂ , where the constant C depends on sup¯ |Du| and .
Double normal derivative estimate on ∂ . We may consider any boundary point. Without loss of generality, we may take it to be the origin and the x n -axis in the direction of the interior normal. From the uniform A-convexity of and the regularity of A, there exists a C 2 defining function in¯ satisfying together with the inequality in a neighborhood N of ∂ , whenever D ν u ≥ ϕ(x), where δ 0 is a positive constant and I denotes the identity matrix. We employ the auxiliary function where β and d 0 are two positive constants to be determined. By the continuity of D ν φ and D ν φ = -1 on ∂ , there exists a small constant d 0 such that D ν φ ≤ -1 2 in d 0 . It follows from (16) that where T = n i=1ũ ii , the first inequality is obtained by using (20) in Lemma 2.1, the second inequality is obtained by using the fact that By a direct calculation, we obtain where the first and second inequalities are established by using (34) and (6), respectively. Therefore, from (32), (33), and (34), we have We If max¯ d 0 w = w(x) for somex ∈¯ d 0 . The proof of the double normal derivative estimate of u on ∂ splits into three stages according to whetherx ∈ 1 d 0 orx ∈ 2 d 0 , or neither.
We may assume thatũ ij (x) is diagonal. Without loss of generality, we may assume that u 11 (x) ≥ 1. At the pointx, we have whereũ 11 (x) ≥ 1 is used in the last inequality. Atx, we can assume that √ T ≥ 3C δ 0 √ n-1 , otherwise we have already got the global second order derivative estimates. Therefore at x, by (35), we have where we choose β ≥ 3C δ 0 + 1 such that the last inequality holds. It is a contradiction with (36). Therefore,x / ∈ 1 d 0 .
where β is chosen large enough such that the last inequality holds. On ∂ , we have Therefore, for all x ∈¯ d 0 , we have Since w(0) = 0, we have which implies Remark 6 In the above proof, we only use (30) in a neighborhood N of ∂ whenever D ν u ≥ ϕ(x). If D ν u < ϕ(x) in a neighborhood N 0 of some boundary point x 0 ∈ ∂ , from the boundary condition, we directly get an upper bound D νν u(x 0 ) ≤ D ν ϕ(x 0 ). While the lower bound D νν u ≥ n i,j A ij ν i ν j at x 0 can be derived from the ellipticity.
Based on the mixed tangential normal derivative estimate and the double normal derivative estimate on ∂ , we now use the method in [15] to prove Theorem 1.1, which is a modification of the original method in [14].
Proof of Theorem 1.1 We modify the elliptic subsolution u by adding a perturbation function aφ, where a is a small positive constant and φ is the defining function of the domain satisfying (29). Note that the function u X := uaφ is still uniformly elliptic in if a is a sufficiently small. By a direct computation, we have By a direct calculation, we also have We employ the auxiliary function in¯ × R n , where |ξ | = 1, α and κ are positive constants to be determined, v is given by and ξ is given by Assume that G attains its maximum at x 0 ∈ and ξ = ξ 0 . Let then the function H attains its maximum at x 0 ∈ and ξ = ξ 0 . From now on, all the calculations will be made at the point x = x 0 and ξ = ξ 0 unless otherwise specified. At x 0 , we obtain Therefore, at x 0 we have By a direct computation and using (17)- (19), we have where the first inequality is established by (21) in Lemma 2.1, and the second inequality is obtained by using the A3w condition. Note that the term n k=1B p k D ii u k in (21) does not appear on the right-hand side of the first inequality in (57) since it can be subtracted using the definition of the linearized operator L in (23). By a direct computation, we obtain where (20) in Lemma 2.1 is used in the second inequality. It follows from (16) and (20) in Lemma 2.1 that Note that there is no 1 +ũ ii in the coefficient of B -1 n-1 in (59) since the termB p l D l u k is already subtracted in L(u k ).
We shall treat the first two terms on the right-hand side of (62). For this purpose, we will make a more detailed calculations than [15] since the equality D i A ξξ = D ξ A iξ does not hold in general. Because in (62) the vector ξ is not equal to e 1 , the situation here is different from the formula below (24) in [21]. Next, we divide the discussions into the two cases (a) and (b) assumed in the statement of this theorem.
Next, we consider the case x 0 ∈ ∂ , namely the function G in (50) attains its maximum over¯ at x 0 ∈ ∂ and a unit vector ξ . The estimation of the rest of the Hessian D 2 u splits into two stages according to a different direction of ξ .
Case (i). ξ tangential. Since v (x 0 , ξ ) = 0, at x 0 we obtain where c 0 = Ka inf(e K(u-u) ), M = max x∈∂ |D k u(ϕ k -D i uD k ν i )|. The above inequality gives a relationship betweenũ ξξ (x 0 ) and D ν u ξξ (x 0 ) at x 0 , namely On the other hand, by tangential differentiating the boundary condition twice, we obtain Hence, for the tangential direction ξ at x 0 , we have where the double normal derivative estimate (46) on ∂ is used in the second inequality. Thus, at x 0 we have where the mixed tangential normal derivative estimate (28) and the double normal derivative estimate (46) are used to obtain the last inequality, and the constant C changes from line to line. Without loss of generality, we can assume the normal at x 0 to be ν = (0, . . . , 1), and we can assume {ũ ij (x 0 )} i,j<n is diagonal with maximum eigenvalueũ 11 (x 0 ) > 1. Then, at x 0 we obtain (κc 0 -αM)ũ ξξ ≤ C(ũ 11 + 1).
Since G(x 0 , e 1 ) ≤ G(x 0 , ξ ), we havẽ The C 1 estimate is established in [15,22] in a general form. For convenience, we repeat it here.
where μ 0 and δ are nonnegative constants. Then we have the estimate where C depends on μ 0 , δ, , and sup |u|.
With the preparations of a priori estimates (100), (102), and (11), we are able to prove the existence result in Theorem 1.2.
Proof of Theorem 1.2 We shall prove the existence of a viscosity solution u of the Neumann problem of the degenerate Monge-Ampère equation (1)-(2) by approximations.
We consider the following approximating problem: det D 2 u -A(·, Du) = B > 0, in , where B (·, u, Du) := B(·, u, Du) + for small positive constant . If B satisfies conditions (6)- (8), it is obvious that B still satisfies the same conditions. By Remark 1, there exists a unique elliptic solution u ∈ C 3,α of problem (103), where α ∈ (0, 1). Sinceū is still the supersolution of problem (103), and for sufficiently small > 0, u is still the subsolution of problem (103) as well, we get the uniform estimate for sufficiently small > 0. Notice that the gradient estimate in Theorem 4.1 is independent of the lower bound of the right-hand side term of the equation. Therefore, we have where the constant C is independent of . By Theorem 1.1, we obtain the uniform second derivative estimate sup D 2 u ≤ C.
From the remarks in the paragraph below Theorem 1.1 (or from the proof of Theorem 1.1), we know that the constant C in (106) is independent of . With the uniform estimates (104), (105), and (106), we can pass the limit to u . Now letting → 0 (passing to a subsequence if necessary), we get a viscosity solution u of the degenerate Monge-Ampère equation (1) with the Neumann boundary condition (2), which is degenerate elliptic and satisfies the estimate Therefore, the proof of Theorem (1.2) is completed.
Remark 7 In Theorem 1.2, we proved the existence of a viscosity solution u ∈ C 1,1 (¯ ) by smooth approximations. Such uniqueness of viscosity solutions in our A z = 0, B z > 0, and ϕ z = 0 case is not covered in Sect. 4 of [23]. It would be interesting to investigate the uniqueness of the viscosity solutions of the degenerate Monge-Ampère equation (1) satisfying the Neumann boundary condition (2).
Remark 8 In the current paper, we show that "the existence of a subsolution implies the existence of a solution", which is the new point of the existence result to the Neumann problem (1)-(2) of degenerate Monge-Ampère type equations. This idea is also new even for nondegenerate Monge-Ampère type equations. For the nondegenerate case in [15], for a more general problem with A and ϕ depending also on u, det D 2 u -A(·, u, Du) = B(·, u, Du) in , the result is that "the existence of a supersolution implies the existence of a solution". If we consider such a general problem (108) in the degenerate elliptic setting, the supersolution u is only degenerate elliptic and is not enough to construct a barrier function. It needs new ideas to solve the general Neumann problem (108) in the degenerate case.