GREEN S FUNCTION FOR SECOND ORDER ELLIPTIC EQUATIONS IN NON-DIVERGENCE FORM arxiv:1712.02160v2 [math.ap] 17 Jul 2018 SUKJUNG HWANG AND SEICK KIM Abstract. We construct the Green function for second order elliptic equations in non-divergence form when the mean oscillations of the coefficients satisfy the Dini condition and the domain has C 1,1 boundary. We also obtain pointwise bounds for the Green functions and its derivatives. 1. Introduction and main results Letbe a bounded C 1,1 domain inr n with n 3. We consider a second-order elliptic operator L in non-divergence form n Lu= a ij (x)d ij u. (1.1) i,j=1 We assume that the coefficient A := (a ij ) is an n n symmetric matrix-valued function defined onr n, which satisfies the uniform ellipticity condition n λ ξ 2 a ij (x)ξ i ξ j Λ ξ 2, ξ=(ξ 1,...,ξ n ) R n, x R n (1.2) i,j=1 for some constants 0<λ Λ. In this article, we are concerned with construction and pointwise estimates for the Green s function G(x, y) of the non-divergent operator L (1.1) in. Unlike the Green s function for uniformly elliptic operators in divergence form, the Green s function for non-divergent elliptic operators do not necessarily enjoy usual pointwise bound G(x, y) c x y 2 n (1.3) even in the case when the coefficient A is uniformly continuous; see [1]. On the other hand, in the case when the coefficient A is Hölder continuous, then it is well known that the Green s function satisfies the pointwise bound (1.3). In this perspective, it is an interesting question to ask what is the minimal regularity condition to ensure the Green s function to have the pointwise bound (1.3). We shall show that if the coefficient A is of Dini mean oscillation, then the Green s function exists and satisfies the pointwise bound (1.3). We shall say that a function is of Dini mean oscillation if its mean oscillation satisfies the Dini condition. Below is a more precise formulation. 2010 Mathematics Subject Classification. Primary 35J08, 35B45 ; Secondary 35J47. S. Hwang is partially supported by the Yonsei University Research Fund (Post Doc. Researcher Supporting Program) of 2017 (project no.: 2017-12-0031). S. Kim is partially supported by National Research Foundation of Korea (NRF) Grant No. 2016R1D1A1B03931680. 1
2 S. HWANG AND S. KIM For x R n and r>0, we denote by B(x, r) the Euclidean ball with radius r centered at x, and denote (x, r) := B(x, r). We shall say that a function g : Ris of Dini mean oscillation inif the mean oscillation functionω g :R + Rdefined by ω g (r) := sup g(y) ḡ (x,r) dy where ḡ (x,r) := g (x,r) (x,r) x satisfies the Dini condition; i.e., 1 ω g (t) dt<+. t 0 It should be clear that if g is Dini continuous, then g is of Dini mean oscillation. However, the Dini mean oscillation condition is strictly weaker than the Dini continuity; see [4] for an example. Also if g is of Dini mean oscillation, then it is routine to verify that g is uniformly continuous in with its modulus of continuity controlled byω g. The formal adjoint operator L is given by n L u= D ij (a ij (x)u). i,j=1 We need to consider the boundary value problem of the form L u=div 2 g+ f in, u= gν ν on, (1.4) Aν ν where g=(g ij ) is an n n matrix-valued function, div 2 g= n i,j=1 D ijg ij, andνis the unit exterior normal vector of. For g L p () and f L p (), where 1<p< and 1 p + 1 p = 1, we say that u in L p () is an adjoint solution of (1.4) if u satisfies ulv= tr(gd 2 v)+ f v (1.5) for any v in W 2,p () W 1,p 0 (). The following lemma is quoted from [7, Lemma 2]. Lemma 1.6. Let 1<p< and assume that g L p () and f L p (). Then there exists a unique adjoint solution u in L p (). Moreover, the following estimates holds u L p () C [ g L p ()+ f L p ()], where a constant C depends on, p, n,λ,λ, and the continuity of A. It is also known that if f L p () with p> n 2, then the adjoint solution of the problem L u= f in, u=0 on, (1.7) is uniformly continuous in; see Theorem 1.8 in [5]. We say is C k,dini if for each point x 0, there exist a constant r > 0 independent of x 0 and a C k,dini function (i.e., C k function whose kth derivatives are uniformly Dini continuous)γ : R n 1 R such that (upon relabeling and
GREEN FUNCTION FOR NON-DIVERGENCE ELLIPTIC EQUATION 3 reorienting the coordinates axes if necessary) in a new coordinate system (x, x n )= (x 1,..., x n 1, x n ), x 0 becomes the origin and B(0, r)={x B(0, r) : x n >γ(x 1,..., x n 1 )}, γ(0 )=0. A few remarks are in order before we state our main theorem. There are many papers in the literature dealing with the existence and estimates of Green s functions or fundamental solutions of non-divergence form elliptic operators with measurable or continuous coefficients. For our best knowledge, the first author who considered Green s function for non-divergence form elliptic operators with measurable coefficients is Bauman [2, 3], who introduced the concept of normalized adjoint solutions; see also Fabes et al. [9]. Fabes and Stroock [8] established L p -integrablity of Green s functions for non-divergence form elliptic operators with measurable coefficients. Krylov [13] showed that the weak uniqueness property holds for solutions of non-divergence form elliptic equations in if and only if there is a unique Green s function in. Escauriaza [6] established bounds for fundamental solution for non-divergence form elliptic operators in terms of nonnegative adjoint solution. We would like to thank Luis Escauriaza for brining our attention to these results in the literature. Now, we state our main theorem. Theorem 1.8. Letbe a bounded C 1,1 domain inr n with n 3. Assume the coefficient A=(a ij ) of non-divergent operator L in (1.1) satisfies the uniform ellipticity condition (1.2) and is of Dini mean oscillation in. Then, there exists a Green s function G(x, y) (for any x, y, x y) of the operator L inand it is unique in the following sense: if u is the unique adjoint solution of the problem (1.7), where f L p () with p> n 2, then u is represented by u(y)= G(x, y) f (x) dx. (1.9) The Green function G(x, y) satisfies the following pointwise estimates: G(x, y) C x y 2 n, (1.10) D x G(x, y) C x y 1 n, (1.11) where C=C(n,λ,Λ,,ω A ). Moreover, if the boundary is C 2,Dini, then we have where C=C(n,λ,Λ,,ω A ). D 2 x G(x, y) C x y n, (1.12) Remark 1.13. In the proof of Theorem 1.8, we will construct the Green s function G (x, y) for the adjoint operator L as a by-product. It is characterized as follows: for q> n 2 and f Lq (), if v W 2,q () W 1,q () is the strong solution of 0 Lv= f in, v=0 on, then, we have the representation formula v(y)= G (x, y) f (x) dx. Also, in the proof of Theorem 1.8, we shall show that G(x, y)=g (y, x), x, y, x y. (1.14) Finally, by the maximum principle, it should be clear that G(x, y) 0.
4 S. HWANG AND S. KIM 2. Proofof Theorem 1.8 2.1. Construction of Green s function. To construct Green s function, we follow the scheme of [12], which in turn is based on [11]. For y andǫ>0, let v ǫ = v ǫ;y W 2,2 () W 1,2 () be a unique strong solution of the problem 0 1 Lv= (y,ǫ) χ (y,ǫ) in, v=0 on. (2.1) Note that A is uniformly continuous in with its modulus of continuity controlled byω A. Therefore, the unique solvability of the problem (2.1) is a consequence of standard W 2,p theory. Also, by the same theory, we see that v ǫ belong to W 2,p () for any p (1, ) and we have an estimate v ǫ W 2,p () Cǫ n+ n p, 0< ǫ<diam, (2.2) where C=C(n,λ,Λ, p,,ω A ). In particular, we see that v ǫ is continuous in. Next, for f C c (), consider the adjoint problem L u= f in, u=0 on. By Lemma 1.6, there exists a unique adjoint solution u in L 2 (), and by (1.5), we have f v ǫ. (2.3) u= (y,ǫ) Let w be the solution of the Dirichlet problem w= f in, w=0 on. (2.4) By the standard L p theory and Sobolev inequality, for 1<p< n 2, we have and for 1<p<n, we have w L q () C(n,, p) f L p (), where 1 q = 1 p 2 n, (2.5) w L q () C(n,, p) f L p (), where 1 q = 1 p 1 n. In particular, for n 2 < p<n, we have by Morrey s theorem that [w] C 0,µ () C(n,, p) f L p (), where µ=1 n q = 2 n p. (2.6) Hereafter, we set g := wi. Note that by (1.7) and (2.4), u L 2 () is an adjoint solution of L u=div 2 g in, u=0 on. (2.7) By Lemma 1.6 and (2.5), we see that u L q () for any 1<q< and satisfies u L q () C f L nq/(n+2q) (), where C=C(n, q,λ,λ,,ω A ). (2.8) Also, by (2.6), we see that g is of Dini mean oscillation inwith ω g (t) C f L p () t 2 n p. Therefore, by Theorem 1.8 of [5], we see that u C(). As a matter of fact, Lemma 2.27 of [5] and Theorem 1.10 of [4] with a scaling argument should reveal that for any x 0 and 0<r<diam, we have ( ) sup u C r n u L 1 ((x 0,r))+ r 2 n p f L p (), (x 0, 1 2 r)
GREEN FUNCTION FOR NON-DIVERGENCE ELLIPTIC EQUATION 5 where C=C(n, p,λ,λ,,ω A ). In particular, if f is supported in(y, r), then by (2.8) and Hölder s inequality, we have ( ) sup u C r n q f L ((y,r))+ r 2 n p f nq/(n+2q) L p ((y,r)) Cr 2 f L ((y,r)). (2.9) (y, 1 2 r) Therefore, if f is supported in(y, r), then it follows from (2.3) and (2.9) that f v ǫ Cr 2 f L ((y,r)), 0< ǫ< 1 2 r. (y,r) By the converse to Hölder s inequality, we obtain v ǫ L 1 ((y,r)) Cr 2, 0< ǫ< 1 2r, 0< r<diam, (2.10) where C=C(n,λ,Λ,,ω A ). We define the approximate Green s function G ǫ (x, y)=v ǫ,y (x)=v ǫ (x). Lemma 2.11. Let x, y with x y. Then where C=C(n,λ,Λ,,ω A ). G ǫ (x, y) C x y 2 n, 0< ǫ< 1 x y, (2.12) 3 Proof. Let r= 2 3 x y. Ifǫ< 1 2 r, then v ǫ satisfies Lv ǫ = 0 in(x, r). Therefore, by the standard elliptic estimate (see [10, Theorem 9.26]) and (2.10), we have v ǫ (x) Cr n v ǫ L1 ((x,r)) Cr n v ǫ L1 ((y,3r)) Cr 2 n. Lemma 2.13. For any y and 0<ǫ<diam, we have G ǫ (x, y) 2n n 2 dx Cr n, r>0, (2.14) \B(y,r) D 2 xg ǫ (x, y) 2 dx Cr n, r>0, (2.15) \B(y,r) where C=C(n,λ,Λ,,ω A ). Proof. We first establish (2.14). In the case when r>3ǫ, we get from (2.12) that G ε (x, y) 2n n 2 dx C x y 2n dx Cr n. \B(y,r) \B(y,r) In the case when r 3ǫ, by (2.2) with p= 2n n+2 and the Sobolev inequality, we have v ǫ L 2n n 2 () C v ǫ W 2, 2n n+2 () Cǫ1 n 2 Cr 1 n 2, (2.16) and thus we still get (2.14). Next, we turn to the proof of (2.15). It is enough to consider the case when r > 2ǫ. Indeed, by (2.2), we have D 2 v ǫ 2 \B(y,r) D 2 v ǫ 2 Cǫ n Cr n. For g C c (\B(y, r)), let u L 2 () be an adjoint solution (2.7) so that we have u= (y,ǫ) tr(gd 2 v ǫ ). (2.17)
6 S. HWANG AND S. KIM Since g=0in(y, r), we see that u is continuous on(y, 1 2r) by [5, Theorem 1.8]. In fact, it should follow from [5, Lemma 2.27] that sup (y, 1 2 r) u Cr n u L 1 ((y,r)). Therefore, by Hölder s inequality and Lemma 1.6, we have sup u Cr n 2 u L 2 ((y,r)) r n 2 g L 2 (). (y, 1 2 r) Since g is supported in\b r (y), by (2.17) and the above estimate, we have tr(gd 2 v ǫ ) Cr n 2 g L 2 (\B(y,r)). \B(y,r) Therefore, (2.15) follows from the converse to the Hölder s inequality. Lemma 2.18. For any y and 0<ǫ<diam, we have where C=C(n,λ,Λ,,ω A ). Proof. We first establish (2.19). Let {x : G ǫ (x, y) >t} Ct n n 2, t>0, (2.19) {x : D 2 xg ǫ (x, y) >t} Ct 1, t>0, (2.20) A t ={x : G ǫ (x, y) >t} and take r=t 1 n 2. Then, by (2.14), we get A t \ B(y, r) t 2n n 2 G ǫ (x, y) 2n n 2 dx Ct 2n n n 2 t n 2 = Ct n A t \B(y,r) Since A t B(y, r) Cr n = Ct n n 2, obviously we thus obtain (2.19). Next, we prove (2.20). Let A t ={x : D 2 x G ǫ(x, y) > t} and take r=t 1 n. Then, by (2.15), we have A t \ B(y, r) t 2 D 2 xg ǫ (x, y) 2 dx Ct 2 t=ct 1. A t \B(y,r) Since A t B(y, r) Cr n = Ct 1, we get (2.20). We are now ready to construct a Green s function. By Lemma 2.13, for any r>0, we have sup G ǫ (, y) W 2,2 (\B(y,r)) <+. 0<ǫ<diam Therefore, by applying a diagonalization process, we see that there exists a sequence of positive numbers{ǫ i } i=1 with lim i ǫ i = 0 and a function G(, y), which belongs to W 2,2 (\B(y, r)) for any r>0, such that G ǫi (, y) G(, y) weakly in W 2,2 (\B(y, r)), r>0. (2.21) Note that by (2.15), we have D 2 xg(x, y) 2 dy Cr n, r>0. (2.22) \B(y,r) n 2.
GREEN FUNCTION FOR NON-DIVERGENCE ELLIPTIC EQUATION 7 Also, by passing to a subsequence if necessary, we derive from (2.14) that G(x, y) 2n n 2 dx Cr n, r>0, \B(y,r) On the other hand, (2.19) implies that for 1<p< n n 2, we have sup G ǫ (, y) L p ()<+. 0<ǫ<diam Therefore, by passing to a subsequence if necessary, we see that G ǫi (, y) G(, y) weakly in L p (), 1< p< n n 2. For f L q () with q> n 2, let u be the unique adjoint solution in Lq () of the problem (1.7). Then by (2.3), we have u(x) dx= (y,ǫ i ) G ǫi (x, y) f (x) dx. By taking the limit, we get the representation formula (1.9), which should yield the uniqueness of the Green s function. Finally, from (2.21) and (2.1), we find that G(, y) belongs to W 2,2 (\B(y, r)) and satisfies LG(, y)=0 in\b(y, r) for all r>0. Since A is uniformly continuous in, by the standard W 2,p theory, we then see that G(, y) is continuous in\ {y}. Moreover, by the same reasoning, we see from Lemma 2.13 that G ǫ (, y) is equicontinuous on\b(y, r) for any r>0. Therefore, we may assume, by passing if necessary to another subsequence, that G ǫi (, y) G(, y) uniformly on \B(y, r) for any r>0. (2.23) In particular, from Lemma 2.11, we see that Therefore, we have shown (1.10). G(x, y) C x y 2 n. 2.2. Construction of Green s function for the adjoint operator. For y and ǫ>0, consider the following adjoint problem: L u= 1 (y,ǫ) χ (y,ǫ) in, u=0 on. (2.24) By Lemma 1.6, there exists a unique adjoint solution u ǫ = u ǫ;y in L 2 () such that u ǫ L 2 () Cǫ n 2, 0< ǫ<diam, where C=C(n,λ,Λ,,ω A ). Next, for f C c (), consider the problem Lv= f in, v=0 on. By the standard L p theory and Sobolev inequality, we have v L q () C f L nq/(n+2q) (), where C=C(n,λ,Λ,, q,ω A ). Also, we have ( ) sup v C r n v L 1 ((x 0,r))+ r 2 f L ((x 0,r)), x0, 0< r<diam. (x 0, 1 2 r)
8 S. HWANG AND S. KIM Therefore, if f is supported in(y, r), then by the above estimates and Hölder s inequality, we get sup (y, 1 2 r) v Cr 2 f L ((y,r)). From the identity we then see that (c.f. (2.10) above) v= (y,ǫ) f u ǫ, (2.25) u ǫ L 1 ((y,r)) Cr 2, 0< ǫ< 1 2r, 0< r<diam, where C=C(n,λ,Λ,,ω A ). We define G ǫ (x, y)=u ǫ,y(x)=u ǫ (x). Then, similar to Lemma 2.11, for x, y with x y, we have G ǫ(x, y) C x y 2 n, 0< ǫ< 1 x y. (2.26) 3 Indeed, if we set r= 2 3 x y, then forǫ< 1 2 r, we have L u ǫ = 0 in(x, r). Since A is of Dini mean oscillation, by [5, Lemma 2.27], we have u ǫ (x) Cr n u ǫ L 1 ((x,r)) Cr n u ǫ L 1 ((y,3r)) Cr 2 n, which should yield (2.26). By using (2.26) and following the proof of Lemma 2.13, we get the following estimate, which is a counterpart of (2.14). For any y and 0<ǫ<diam, we have G 2n ǫ (x, y) n 2 dx Cr n, r>0. (2.27) \B(y,r) 1 Indeed, by taking f= (y,ǫ) χ (y,ǫ) and q= 2n n 2 in (2.8), we have u ǫ L 2n n 2 () C f L 2n n+2 () Cǫ1 n 2, which corresponds to (2.16). Then, by the same proof of Lemma 2.13, we get (2.27). By using (2.27) and proceeding as in the proof of (2.19), we obtain {x : G n ǫ (x, y) >t} Ct n 2, t>0, which in turn implies that for 0<p< n n 2, there exists a constant C p such that G ǫ (x, y) p dx C p, y, 0< ǫ<diam. Therefore, for any 1<p< n n 2, we obtain sup G ǫ (, y) L p ()<+, 0<ǫ<diam and thus, there exists a sequence of positive numbers{ǫ j } j=1 with lim j ǫ j = 0 and a function G (, y) L p () such that G ǫ j (, y) G (, y) weakly in L p (). (2.28) For f L q () with q> n 2, let v W2,q () W 1,q 0 () be the strong solution of Lv= f in, v=0 on.
GREEN FUNCTION FOR NON-DIVERGENCE ELLIPTIC EQUATION 9 Then, we have (c.f. (2.25) above) v(x) dx= (y,ǫ j ) G ǫ j (x, y) f (x) dx, and thus, by taking the limit, we also get the representation formula v(y)= G (x, y) f (x) dx. Finally, by (2.24), (2.28), and Theorem 1.10 of [4], we see that G (, y) is continuous away from its singularity at y. 2.3. Proof of symmetry (1.14). For xy in, choose two sequences{ǫ i } i=1 and {δ j } j=1 such that 0<ǫ i,δ j < 1 3 x y for all i, j andǫ i,δ j 0 as i, j. From the construction of G ǫ (, y) and G ǫ(, x), we observe the following: G ǫi (, y)= (x,δ j ) G δ j (, x). (y,ǫ i ) By (2.28) and continuity of G (, x) away from singularity, we get lim lim i j G δ j (, x)= lim (y,ǫ i ) i G (, x)=g (y, x). (y,ǫ i ) On the other hand, by continuity of G ǫi (, y) and (2.23) we obtain lim lim i j We have thus shown that G ǫi (, y)= lim G ǫi (x, y)=g(x, y). (x,δ j ) i G(x, y)=g (y, x), for any x y. So far, we have seen that there is a subsequence of{ǫ i } tending to zero such that G ǫki (, y) G(, y). However, for x y, we have G ǫ (x, y)= lim j That is, we have Therefore, we find that G ǫ (, y)= lim (x,δ j ) j G δ j (, x)= (y,ǫ) G ǫ (x, y)= G(x, z) dz. (y,ǫ) lim G ǫ (x, y)=g(x, y) for any x y. ǫ 0 G (, x)= (y,ǫ) G(, y). (y,ǫ) 2.4. Proof of estimates (1.11) and (1.12). We now show (1.11) and (1.12). Let v=g(, y) and r= 1 3 x y. Note that we have Lv=0 in (x, 2r), v=0 on (x, 2r). By the standard elliptic theory (e.g. L p theory), we have Therefore, by (1.10), we get from which (1.11) follows. Dv L ((x,r)) Cr 1 n v L 1 ((x,2r)). Dv(x) Cr 1 n v L 1 ((y,4r)) Cr 1 n,
10 S. HWANG AND S. KIM In the case whenhas C 2,Dini boundary, by Theorem 1.5 of [5], we have sup D 2 v Cr n D 2 v L 1 ((x,2r)). (x,r) Therefore, by Hölder s inequality and (2.22), we have D 2 v(x) Cr n 2 D 2 v L 2 ((x,2r)) Cr n 2 D 2 v L 2 (\B(y,r)) Cr n, from which (1.12) follows. References [1] Bauman, Patricia. Equivalence of the Green s functions for diffusion operators inr n : a counterexample. Proc. Amer. Math. Soc. 91 (1984), no. 1, 64 68. [2] Bauman, Patricia. Positive solutions of elliptic equations in nondivergence form and their adjoints. Ark. Mat. 22 (1984), no. 2, 153 173. [3] Bauman, Patricia. A Wiener test for nondivergence structure, second-order elliptic equations. Indiana Univ. Math. J. 34 (1985), no. 4, 825 844. [4] Dong, Hongjie; Kim, Seick. On C 1, C 2, and weak type-(1, 1) estimates for linear elliptic operators. Comm. Partial Differential Equations 42 (2017), no. 3, 417 435. [5] Dong, Hongjie; Escauriaza, Luis; Kim, Seick. On C 1, C 2, and weak type-(1, 1) estimates for linear elliptic operators: part II. Math. Ann. 370 (2018), no. 1-2, 447 489. [6] Escauriaza, Luis. Bounds for the fundamental solution of elliptic and parabolic equations in nondivergence form. Comm. Partial Differential Equations 25 (2000), no. 5-6, 821 845. [7] Escauriaza, Luis; Montaner, Santiago. Some remarks on the L p regularity of second derivatives of solutions to non-divergence elliptic equations and the Dini condition. Rend. Lincei Mat. Appl. 28 (2017) 49 63. [8] Fabes, E. B.; Stroock, D. W. The L p -integrability of Green s functions and fundamental solutions for elliptic and parabolic equations. Duke Math. J. 51 (1984), no. 4, 997 1016. [9] Fabes, E.; Garofalo, N.; Mar?-Malave, S.; Salsa, S. Fatou theorems for some nonlinear elliptic equations. Rev. Mat. Iberoamericana 4 (1988), no. 2, 227 251. [10] Gilbarg, David; Trudinger, Neil S. Elliptic partial differential equations of second order. Reprint of the 1998 edition. Classics in Mathematics. Springer-Verlag, Berlin, 2001. [11] Grüter, Michael; Widman, Kjell-Ove. The Green function for uniformly elliptic equations. Manuscripta Math. 37 (1982), no. 3, 303 342. [12] Hofmann, Steve; Kim, Seick. The Green function estimates for strongly elliptic systems of second order. Manuscripta Math. 124 (2007), no. 2, 139 172. [13] Krylov, N. V. On one-point weak uniqueness for elliptic equations. Comm. Partial Differential Equations 17 (1992), no. 11-12, 1759 1784. (S. Hwang) Department of Mathematics, Yonsei University, 50 Yonsei-ro, Seodaemun-gu, Seoul 03722, Republic of Korea E-mail address: sukjung hwang@yonsei.ac.kr (S. Kim) Department of Mathematics, Yonsei University, 50 Yonsei-ro, Seodaemun-gu, Seoul 03722, Republic of Korea E-mail address: kimseick@yonsei.ac.kr