THE SCHAUDER ESTIMATES FOR HIGHER-ORDER PARABOLIC SYSTEMS WITH TIME IRREGULAR COEFFICIENTS arxiv: v1 [math.ap] 18 Jul 2013

Transcription

1 THE SCHAUDER ESTIMATES FOR HIGHER-ORDER PARABOLIC SYSTEMS WITH TIME IRREGULAR COEFFICIENTS arxiv: v1 [math.ap] 18 Jul 2013 HONGJIE DONG AND HONG ZHANG Abstract. We prove Schauder estimates for solutions to both divergence and non-divergence type higher-order parabolic systems in the whole space and the half space. We also provide an existence result for divergence type systems in a cylindrical domain. All coefficients are assumed to be only measurable in the time variable and Hölder continuous in the spatial variables. 1. Introduction This paper is devoted to the Schauder estimates for divergence and nondivergence type higher-order parabolic systems. It is well known that the Schauder estimates play an important role in the existence and regularity theories for elliptic and parabolic equations and systems, which have been studied by many authors. For parabolic equations with constant coefficients, the classical approach is to study the fundamental solutions, see e.g., [14, 7]. Moreover, combining the classical approach with the argument of freezing the coefficients, we can deal with general parabolic equations with smooth coefficients, for instance see [14, 17, 1]. For systems it has become customary to use Campanato s technique, first introduced in [4]. Giaquinta [8] provided a comprehensive explanation for the application of this technique to second-order elliptic systems. Schlag [16] applied this technique to second-order parabolic systems. For higher-order systems, we also refer the reader to [12, Chap. 3] and the references therein. The classical Schauder estimates were established under the assumption that coefficients are regular in both space and time. In this paper, we consider the coefficients which are regular only with respect to spatial variables. This type of coefficient has been studied by several authors mostly for second-order equations; see, for instance, [2, 13, 9, 11, 15]. In [9, 10] Lieberman studied interior and boundary Schauder estimates for secondorder parabolic equations with time irregular coefficients using the maximal principle and a Campanato type approach. More recently, Boccia [3] considered higher-order non-divergence type parabolic systems in the whole space. Key words and phrases. higher-order systems, Schauder estimates. H. Dong was partially supported by the NSF under agreement DMS H. Zhang was partially supported by the NSF under agreement DMS

2 2 H. DONG AND H. ZHANG To present our results precisely, let Lu =, β m where m is a positive integer, A αβ D α D β u, Lu =, β m D α = D α 1 1 Dα d d, α = (α 1,,α d ), D α (A αβ D β u), and, for each α and β, A αβ = [A αβ ij ]n i,j=1 is an n n real matrix-valued function. Moreover the leading coefficients satisfy the so-called Legendre Hadamard ellipticity condition (see (2.1)), which is more general than the strong ellipticity condition. The functions used throughout this paper u = (u 1,,u n ) tr, f = (f 1,,f n ) tr, f α = (f 1 α,,f n α) tr are real vector-valued functions. The parabolic systems which we study are u t ( 1) m Lu = f, u t ( 1) m Lu = D α f α in the whole space, or in the half space or cylindrical domains with the Dirichlet boundary condition u = Du =... = D m 1 u = 0, where the first system is in the non-divergence form and the second system is in the divergence form. We assume that all the coefficients and data are Hölder continuous with respect to the spatial variables and merely measurable with respect to the time variable. For the non-divergence form systems, we prove that D 2m u is Hölder continuous in the whole space, and in the half space all the 2mth order derivatives of u are Hölder continuous up to the boundary with the exception of Dd 2mu, where x d is the normal direction of the boundary. For the divergence form systems, we prove that all the mth order derivatives are Hölder continuous up to the boundary. We also prove an existence theorem for the divergence form systems in a cylindrical domain, provided that the boundary of the domain is sufficiently smooth. To our best knowledge, these results are new for higher-order systems and they extend the corresponding results found in Lieberman [9] for second-order scalar equations. In the special case of second-order parabolic systems, compared to [16] our conditions on the coefficients and data are more general. In particular, we do not require the data to be vanish on the lateral boundary of the domain, i.e., the compatibility condition in [16]. For the proof, we use some results in [5], in which Dong and Kim proved the L p estimates for the divergence and non-divergence type higher-order parabolic systems, with the same type of coefficients considered in this paper. Letusgivetheoutlineoftheproofs. Inthedivergencecase, withtheclassical

3 L 2 estimates, we prove D m u (D m u) Qr(X 0 ) 2 dxdt Q r(x 0 ) C( r R )2d2m where u is a solution of Here Q R (X 0 ) SCHAUDER ESTIMATES 3 D m u (D m u) QR (X 0 ) 2 dxdt, r R, u t ( 1) m L 0 u = 0 in Q 2R. L 0 = α = β =m D α (A αβ (t)d β ), and X 0 Q R. The coefficients A αβ which here depend only on t, are called simple coefficients. For a similar inequality corresponding to the boundary estimates, we combine the L p estimates established in [5] and the Sobolev embedding theorem to prove the Hölder continuity and obtain, for instance, the following mean oscillation type estimate for systems with simple coefficients, D Q d m u (Dm d u) Q r (X 0 ) 2 dxdt r (X 0 ) C( r R )2γd2m Dd m u (Dm d u) Q R (X 0) 2 dxdt, r R, Q R (X 0) where X 0 {x d = 0} Q R and γ (0,1) are arbitrary, and u satisfies u t ( 1) m L 0 u = 0 in Q 2R with the Dirichlet boundary condition u = Du =... = D m 1 u = 0 on {x d = 0}. Similar interior and boundary estimates for the non-divergence form systems can be established in the same fashion. In particular, we can always differentiate the system { ut ( 1) m L 0 u = 0 in Q 2R u = 0, D d u = 0,,D m 1 d u = 0 on {x d = 0} Q 2R with respect to tangential direction x. This together with the Wp 1,2m estimates implies that D x D 2m 1 u is in some Hölder space C a 2m,a (Q R ) with a arbitrarily close to 1 from below, which yields the mean oscillation type estimates. For the coefficients which depend on both t and x, we use the standard argument of freezing the coefficients to obtain the Campanato type estimates and then achieve the estimates of the Hölder norms. The paper is organized as follows. In the next section we introduce some notation and state our main results. Section 3 is devoted to some necessary technical lemmas. In Section 4 we make necessary preparations and in

4 4 H. DONG AND H. ZHANG Section 5 prove our main result for divergence type systems, Theorem 2.1. Section 6 deals with non-divergence type systems. 2. main results We first introduce some notation used throughout the paper. A point in R d is denoted by x = (x 1,,x d ). We also denote x = (x,x d ), where x R d 1. A point in R d1 = R R d = {(t,x) : t R,x R d } is denoted by X = (t,x). For T (, ], set O T = (,T) R d, O T = (,T) Rd, where R d = {x = (x 1,,x d ) R d : x d > 0}. In particular, when T =, we use O = R Rd. Denote B r (x 0 ) = {x R d : x x 0 < r}, Q r (t 0,x 0 ) = (t 0 r 2m,t 0 ) B r (x 0 ), B r (x 0) = B r (x 0 ) R d, Q r (t 0,x 0 ) = Q r (t 0,x 0 ) O. We use the abbreviations, for example, B r if x 0 = 0 and Q r if (t 0,x 0 ) = (0,0). The parabolic boundary of Q r (t 0,x 0 ) is defined to be p Q r (t 0,x 0 ) = [t 0,t 0 r 2m ) B r (x 0 ) {t = t 0 r 2m } B r (x 0 ). The parabolic boundary of Q r (t 0,x 0 ) is We denote p Q r (t 0,x 0 ) = ( p Q r (t 0,x 0 ) O ) (Q r (t 0,x 0 ) {x d = 0}). f,g Ω = Ω f tr g = n j=1 For a function f defined on D R d1, we set Ω f j g j. { } f(t,x) f(s,y) [f] a,b,d sup t s a : (t,x),(s,y) D, (t,x) (s,y), x y b where a,b (0,1]. The Hölder semi-norm with respect to t is denoted by { } f(t,x) f(s,x) f a,d sup t s a : (t,x),(s,x) D, t s, where a (0,1]. We will also use the Hölder semi-norm with respect to x { } f(t,x) f(t,y) [f] a,d sup x y a : (t,x),(t,y) D,x y, and denote f a,d = f L (D) [f] a,d, where a (0,1]. The space corresponding to a,d is denoted by Ca (D). For a (0,1], set f a 2m,a,D = f L (D) [f] a 2m,a,D.

5 SCHAUDER ESTIMATES 5 The space corresponding to a 2m,a,D is denoted by C a 2m,a (D). For a (1,2m), not an integer, we define f a 2m,a,D = f L (D) [D α f]a α α <a 2m,{a},D, where {a} = a [a]. We use the following Sobolev space Wp 1,2m ((S,T) Ω) = {u : u t,d α u L p ((S,T) Ω),0 α 2m}. We denote the average of f in D to be (f) D = 1 f(t,x)dxdt = D D D f(t,x)dxdt. Sometimes we take average only with respect to x. For instance, (f) BR (x 0 )(t) = f(t,x)dx. B R (x 0 ) Throughout this paper, we assume that all the coefficients are measurable and bounded, A αβ K. In addition, we impose the Legendre Hadamard ellipticity condition with constant λ > 0 on the leading coefficients, i.e., A αβ ij ξ iξ j η α η β λ ξ 2 η 2m, (2.1) α = β =m whereξ R n, η R d, and η α = η α 1 1 ηα 2 coefficients if α = β = m. Throughout this paper L 0 u = L 0 u = α = β =m 2 ηα d d. Here we call Aαβ ij the leading D α (A αβ (t)d β u), where A αβ are measurable in t and satisfy (2.1). We are ready to state the main results of the paper. The first result is about the Schauder estimate and solvability for divergence type higher-order systems in cylindrical domains. Theorem 2.1. Assume f α C a for α = m and f α L for α < m, where a (0,1). Suppose that the operator L satisfies the Legendre Hadamard condition, i.e., (2.1), and A αβ C a. Let g be a smooth function in R d1 and Ω C m,a. Then u t ( 1) m Lu = D α f α in (0,T) Ω, u = g, Du = Dg,, D m 1 u = D m 1 g on [0,T) Ω, u = g on {0} Ω

6 6 H. DONG AND H. ZHANG has a unique solution u C am 2m,am ([0,T) Ω). Moreover there exists a constant C = C(d,n,m,λ,K, A αβ a,ω,a) such that, where and u am 2m,am,(0,T) Ω C( u L 2 ((0,T) Ω) F G), (2.2) G = F = [f α ] a f α L D α g a D α g L g t L. The following theorem is regarding a priori interior and boundary Schauder estimates for the non-divergence type systems. Theorem 2.2. Suppose that A αβ C a and f C a, where a (0,1). Let u C loc (Rd1 ) be a solution to u t ( 1) m Lu = f, (2.3) where L satisfies the Legendre Hadamard condition. For any R 1, I. Interior case: if (2.3) holds in Q 4R, then there exists a constant C = C(d,n,m,λ,K, A αβ a,r,a) such that u t a,q R D 2m u a 2m,a,Q R C( f a,q 4R u L2 (Q 4R )); II. Boundary case: if (2.3) holds in Q 4R with the Dirichlet boundary condition on {x d = 0} Q 4R, then there exists a constant C = C(d,n,m,λ,K, A αβ a,r,a) such that [D x D 2m 1 u] a C[f] C D γ u 2m,a,Q R a,q 4R L (Q ). 4R γ 2m We note that the boundary estimate in Theorem 2.2 is optimal even for the heat equation, in the sense that near the boundary Dd 2 u might be discontinuous with respect to x if f has jump discontinuities in t. See, for instance, [18, 3] and [9, 16]. In order to estimate the Hölder semi-norm of Dd 2mu, we need to impose more regularity conditions on Aαβ and f. See Remark 6.3 below for a discussion. 3. Technical preparation We need the following version of Campanato s theorem. The proof can be found in [8] and [16]. Lemma 3.1. (i) Let f L 2 (Q 2 ), a (0,1), and assume that f (f) Qr(t 0,x 0 ) 2 dxdt A 2 r 2a, (t 0,x 0 ) Q 1, (3.1) Q r(t 0,x 0 ) and 0 < r 1. Then f C a 2m,a (Q 1 ) and [f] a 2m,a,Q 1 CA,

7 SCHAUDER ESTIMATES 7 with C = C(d). ii) Let f L 2 (Q 2 ), a (0,1), and assume that (3.1) holds for r < x 0d. Moreover f (f) Q r (t 1,x 1 ) 2 dxdt A 2 r 2a, (t 1,x 1 ) {x d = 0} Q 2, Q r (t 1,x 1 ) and 0 < r 1. Then f C a 2m,a (Q 1 ) and with C = C(d). [f] a 2m,a,Q 1 CA, The following technical lemma will be frequently used in this paper. The proof can be found in [9] and [16]. Lemma 3.2. Let Φ be a nonnegative, nondecreasing function on (0,r 0 ] such that Φ(ρ) A( ρ r )a Φ(r)Br b for any 0 < ρ < r r 0, where 0 < b < a are fixed constants. Then with a constant C = C(A,a,b). Φ(r) Cr b (r b 0 Φ(r 0)B) r (0,r 0 ) Here is another technical lemma. Lemma 3.3. Assume that Ω is an open bounded domain in R d. Then there exists a function ξ C0 (Ω) with unit integral such that for any 0 < α m ξ(y)y α dy = 0. Proof. Let M = m l=1 (ld 1 d 1 Ω ). We look for ξ in the form ξ(y) = a β ξ β (y), β m where β is an n-tuple multi-index and ξ β C0 (Ω), β m, are functions to be chosen. The problem is then reduced to the following linear system of {a β }: a β ξ β (y)y α dy = 1 α =0, α m. β m Ω It suffices to choose ξ β such that the (M 1) (M 1) coefficient matrix A := { Ωξ β (y)y α dy } α, β m is invertible. To this end, we use a perturbation argument. It is easily seen that the matrix { y β y α dy } α, β m Ω

8 8 H. DONG AND H. ZHANG is invertible. Since C0 (Ω) is dense in L 2(Ω), for any ε > 0 and β m, there exists ξ β C0 (Ω) such that ξ β y β L2 (Ω) ε. By the Cauchy Schwarz inequality, Ω ξ β (y)y α dy y β (y)y α dy Cε, Ω where C is a positive constant independent of ε. By the continuity of matrix determinant, upon taking ε sufficiently small, A is still invertible. The lemma is proved. 4. Estimates for divergence type systems This section is devoted to the interior and boundary a priori estimates for divergence type higher-order systems Interior estimates. The main result of this subsection is the following proposition. Proposition 4.1. Assume that R 1 and u Cloc (Rd1 ) satisfies u t ( 1) m Lu = D α f α in Q 2R. Suppose that A αβ C a, f α C a if α = m where a (0,1), and f α L if α < m. Then there exists a constant C(d,n,m,λ,K, A αβ a,r,a) such that [D m u] a 2m,a,Q C( D β u R L (Q 2R ) F), (4.1) where β m u 1 2 a 2m,Q R C( β m D β u L (Q 2R ) F), (4.2) F = [f α ] a,q 2R f α L (Q 2R ). For the proof of the proposition, we first consider homogeneous systems with simple coefficients. Lemma 4.2. Assume that u C loc (Rd1 ) and satisfies u t ( 1) m L 0 u = 0 in Q 2R, (4.3) where R (0, ). Then for any X 0 Q R and r < R, There is a constant C = C(d,n,m,λ) such that u (u) Qr(X 0 ) 2 dxdt C( r R )22md u (u) QR (X 0 ) 2 dxdt. Q r(x 0 ) Q R (X 0 ) (4.4)

9 SCHAUDER ESTIMATES 9 Proof. By scaling and translation of the coordinates, without loss of generality, we can assume that R = 1 and X 0 = (0,0). First we consider the case r 1/4. Since Q r u (u) Qr 2 dxdt Cr 22md sup Q r Du 2 Cr 4m2md sup Q r u t 2. By Lemma 4.3 in [5], sup( Du u t ) C u L2 (Q 1 ). Q 1/4 Hence, Q r u (u) Qr 2 dxdt Cr 22md u 2 L 2 (Q 1 ). Clearly, the inequality above holds true in the case when r [1/4,1) as well. Since u (u) Q1 also satisfies (4.3), substituting u by u (u) Q1, the lemma is proved. Lemma 4.3. Assume that r < R 1 and u C loc (Rd1 ) satisfies u t ( 1) m L 0 u = D α f α in Q 2R, and f α C a, if α = m where a (0,1), and f α L if α < m. Then for any X 0 Q R, there exists a constant C = C(n,m,d,λ) such that Q r(x 0 ) D m u (D m u) Qr(X 0 ) 2 dxdt C( r R )22md Q R (X 0 ) D m u (D m u) QR (X 0 ) 2 dxdt C [f α ] 2 a,q 2R R 2a2md C f α 2 L (Q 2R ) R22md. Proof. Let w be the weak solution of the following system w t ( 1) m L 0 w = D α f α in Q 2R, w = 0, Dw = 0,..., D m 1 w = 0 on p Q 2R. Multiply w to both sides of the system and integrate over Q 2R. Due to the homogeneous Dirichlet boundary condition of w, we can replace D α f α by D α (f α (t,x) f α (t,0)) when α = m. By Young s inequality, we get, for any

10 10 H. DONG AND H. ZHANG ε (0,1), D m w 2 dxdt Q 2R C Q 2R ( 1) m (f α (t,x) f α (t,0))d α wdxdt Q 2R ( 1) α f α D α wdxdt ε D m w 2 L 2 (Q 2R ) f α L (Q 2R ) C(ε) f α (t,x) f α (t,0) 2 L 2 (Q 2R ) Q 2R D α w dxdt. (4.5) Since w = 0, Dw = 0,..., D m 1 w = 0 on p Q 2R, then for α < m and R 1 by the Poincaré inequality D α w 2 dxdt CR 2 D m w 2 dxdt, Q 2R Q 2R where C is a universal constant. By the Cauchy Schwarz inequality and the inequality above, we can estimate the last term in (4.5) as follows: f α L (Q 2R ) D α w dxdt Q 2R f α L (Q 2R ) Q 2R 1/2 ( D α w 2 dxdt) 1/2 Q 2R C f α L (Q 2R )R Q 2R 1/2 ( D m w 2 dxdt) 1/2 Q 2R C(ε) f α 2 L (Q 2R ) R2d2m ε D m w 2 dxdt. Q 2R Choosing ε sufficiently small, we obtain D m w 2 dxdt Q 2R C f α (t,x) f α (t,0) 2 L 2 (Q 2R ) C f α 2 L (Q 2R ) R2d2m C [f α ] 2 a,q 2R R 2ad2m C f α 2 L (Q 2R ) R2d2m. (4.6) We now temporarily assume that A αβ and f α are smooth functions. By the classical theory, we know that w is smooth. Consider v = u w, which is a smooth function as well. Then v satisfies v t ( 1) m L 0 v = 0 in Q 2R.

11 SCHAUDER ESTIMATES 11 Since D m v satisfies the same system as v, applying Lemma 4.2 to D m v, we have D m v (D m v) Qr(X 0 ) 2 dxdt Q r(x 0 ) C( r R )2d2m Q R (X 0 ) By the triangle inequality, D m u (D m u) Qr(X 0 ) 2 dxdt Q r(x 0 ) = C Q r(x 0 ) Q r(x 0 ) D m v (D m v) QR (X 0 ) 2 dxdt. (4.7) D m w (D m w) Qr(X 0 ) D m v (D m v) Qr(X 0 ) 2 dxdt D m w 2 dxdtc By (4.6) and (4.7), we get D m u (D m u) Qr(X 0 ) 2 dxdt Q r(x 0 ) C Q r(x 0 ) Q r(x 0 ) D m w 2 dxdtc( r R )2d2m D m v (D m v) Qr(X 0 ) 2 dxdt. Q R (X 0 ) D m v (D m v) QR (X 0 ) 2 dxdt C D m w 2 dxdtc( r Q 2R R )2d2m D m u (D m u) QR (X 0 ) 2 dxdt Q R (X 0 ) C( r R )2d2m D m u (D m u) QR (X 0 ) 2 dxdt Q R (X 0 ) C [f α ] 2 a,q 2R R 2ad2m C f α 2 L (Q 2R ) R2d2m. Thus the lemma is proved under the additional smoothness assumption. By a standard argument of mollification, it is easily seen that we can remove the smoothness assumption. Applying Lemma 3.2 and Lemma 4.3, we get the Campanato type estimate D m u (D m u) Qr(X 0 ) 2 dxdt Q r(x 0 ) C( r R )2ad2m Q R (X 0 ) D m u (D m u) QR (X 0 ) 2 dxdt C( [f α ] 2 a,q 2R f α 2 L (Q 2R ) )r2ad2m

12 12 H. DONG AND H. ZHANG for any r < R 1 and X 0 Q R. Then, 1 r 2ad2m D m u (D m u) Qr(X 0 ) 2 dxdt Q r(x 0 ) C R 2ad2m D m u (D m u) QR (X 0 ) 2 dxdtcf 2 Q R (X 0 ) C R 2a2md D m u 2 dxdtcf 2. Q 2R Thanks to Lemma 3.1, we get [D m u] a 2m,a,Q R C(R (2ad2m) Q 2R D m u 2 dxdtf 2 ) 1 2. (4.8) As for the coefficients which also depend on x, the estimates follow from the standard argument of freezing the coefficients. Specifically, we first consider the operator L that only has the highest-order terms. We fix a y B R, and define L 0y = D α (A αβ (t,y)d β ). Then u t ( 1) m L 0y u = α = β =m D α f α ( 1) m (L L 0y )u = D α fα, where f α = f α ( 1) m1( f α = f α if α < m. β =m ) (A αβ (t,x) A αβ (t,y))d β u if α = m, Applying the method in proving Lemma 4.3, as (4.6) we obtain D m w 2 L 2 (Q 2R ) C( α = β =m Following the proof of (4.8), we reach [A αβ ] 2 a Dm u 2 L (Q 2R ) R2a2md f α 2 L (Q 2R ) R22md [f α ] 2 a,q 2R R 2a2md ). [D m u] a 2m,a,Q R C(R (2ad2m) D m u 2 dxdt Q 2R [A αβ ] 2 a D m u 2 L (Q 2R ) F2 ) 1 2. (4.9) α = β =m

13 SCHAUDER ESTIMATES 13 For systems with lower-order terms, we move all the lower-order terms to the right-hand side to get u t ( 1) m L h u = D α (f α ( 1) m1 A αβ D β u) β <m D α (f α ( 1) m1 β m A αβ D β u), where L h u denotes the sum of the highest-order terms. With the inequality we just got, we only need to substitute ˆf α = f α ( 1) m1 β Aαβ D β u into the estimate (4.9) and notice that [A αβ D β u] a,q 2R C ( A αβ L [D β u] a,q 2R [A αβ ] a Dβ u L (Q 2R )), β <m β <m when α = m. An easy calculation then completes the proof of (4.1). It remains to prove (4.2). We estimate u(t,x 0 ) u(s,x 0 ), where (t,x 0 ), (s,x 0 ) Q R. Let ρ = t s 1 2m andη(x) = 1 ξ( x ρ d ρ ), whereξ(x) is thefunction in Lemma 3.3. We define u(t,y) = u(t,y) T m,x0 u(t,y)u(t,x 0 ), where T m,x0 u(t,y) is the Taylor expansion of u(t,y) in y at x 0 up to mth order. By the triangle inequality, u(t,x 0 ) u(s,x 0 ) u(t,x 0 ) η(x 0 y)u(t,y)dy u(s,x 0 ) η(x 0 y)u(s,y)dy η(x 0 y)u(s,y)dy η(x 0 y)u(t,y)dy. (4.10) The first two terms on the right-hand side of (4.10) can be estimated in a similar fashion: noting that B ρ η(x)dx = 1, u(t,x 0 ) η(x 0 y)u(t,y)dy Since = η(x 0 y)(u(t,x 0 ) u(t,y))dy. u(t,y) u(t,x 0 ) = u(t,y) T m,x0 u(t,y) C[D m u] a 2m,a,Q 2R y x 0 ma, we get u(t,x 0 ) η(x 0 y)u(t,y)dy C[D m u] a 2m,a,Q 2R ρma. (4.11)

14 14 H. DONG AND H. ZHANG For the last term on the right-hand side of (4.10), by the definition of η, η(x 0 y)u(s,y)dy = η(x 0 y)u(s,y)dy. Therefore, = = t s η(x 0 y)u(t,y)dy η(x 0 y)u(t,y)dy η(x 0 y)u t (τ,y)dydτ. η(x 0 y)u(s,y)dy η(x 0 y)u(s,y)dy Because u satisfies the equation, t s t = s η(x 0 y)u t (τ,y)dydτ η(x 0 y)(( 1) m1 D α (A αβ D β u(τ,y)), β m D α f α )dydτ. (4.12) Since η(x 0 y) has compact support in B ρ (x 0 ), we use f α (f α ) Bρ(x 0 ) to substitute f α when α = m. Moreover for α = m, D α (A αβ D β u) = D α( (A αβ (A αβ ) Bρ(x 0 ))D β u ) D α ((A αβ ) Bρ(x 0 )(D β u (D β u) Bρ(x 0 ))).

15 SCHAUDER ESTIMATES 15 We plug all these into (4.12) and integrate by parts. It follows easily that t η(x 0 y)u t (τ,y)dydτ s, β m Q ρ(x 0 ), β m ([A αβ ] a Dβ u L (Q 2R ) [D β u] a,q 2R A αβ L )ρ a D m η(x 0 y) dydt A αβ L D β u L (Q 2R ) f α L (Q 2R ) Q ρ(x 0 ) C( D m u a 2m,a,Q 2R β <m D m η L [f α ] a,q 2R Q ρ ρ a Q ρ(x 0 ) D α η(x 0 y) dydt D α η(x 0 y) dydt D β u a,q 2R )ρ ma C [f α ] a,q 2R ρ ma β m C ( f α L (Q 2R ) A αβ L D β u L (Q 2R ))ρ m1. Here we used D α η L Cρ d α for any α m. Hence, combining with (4.11), u(t,x 0 ) u(s,x 0 ) Cρ ma ( D m u a 2m,a,Q D β u 2R L (Q 2R ) F), β m which together with (4.1) implies (4.2). This completes the proof of the proposition Boundary estimates. As in the previous subsection, we first consider systems with simple coefficients without lower-order terms: u t ( 1) m L 0 u = 0 in Q 2R, (4.13) u = 0, D d u = 0,, D m 1 d u = 0 on {x d = 0} Q 2R. (4.14) Thanks to the L p estimates in the half space obtained in [5], we are able to derive a local Wp 1,2m estimate for solution of (4.13) and (4.14) by the method in Lemma 4.1 in [5]. Namely, u W 1,2m p (Q R ) C(R,p) u L p(q 2R ), (4.15) wherec(r,p) isaconstant. Now wearereadytoprovethefollowing lemma. Lemma 4.4. Assume that 0 < r R <, γ (0,1), and u C loc (Rd1 ) satisfies (4.13) and (4.14). Then there exists a constant C = C(d,n,m,λ,K,γ),

16 16 H. DONG AND H. ZHANG such that Q r (X 0 ) Q r (X 0 ) Q r (X 0 ) u 2 dxdt C( r R )d4m Q R (X 0) D l d u 2 dxdt C( r R )2(m l)d2m D m d u (Dm d u) Q r (X 0 ) 2 dxdt C( r R )2γd2m Q R (X 0) where X 0 {x d = 0} Q R and 0 < l < m. u 2 dxdt, (4.16) Q R (X 0) D l d u 2 dxdt, (4.17) D m d u (Dm d u) Q R (X 0) 2 dxdt, (4.18) Proof. By scaling and translation of the coordinates, without loss of generality, we can assume R = 1 and X 0 = (0,0). From (4.15), for any r 1 < 2, there exists a constant C = C(r 1 ) so that u W 1,2m 2 (Q r 1 ) C u L 2 (Q 2 ). By the Sobolev embedding theorem, u Lp(Q r 1 ) C u W 1,2m 2 (Q r 1 ), where 1/p > 1/2 1/(d1). Using (4.15) again, we get u W 1,2m p (Q r 2 ) C u L p(q r 1 ) C u L 2 (Q 2 ), where r 2 < r 1 and C = C(r 1,r 2 ). Due to a standard argument of bootstrap withasequenceofshrinkingcylindersq rl, forany p > 2thereisaC depends on p such that u W 1,2m p (Q 1 ) C u L 2 (Q 2 ). (4.19) By the Sobolev embedding theorem, we have Wp 1,2m (Q 1 ) C a 2m,a (Q 1 ), u a 2m,a,Q 1 C u W 1,2m p (Q 1 ), (4.20) where a = 2m (d2m)/p. Since p < can be arbitrarily large, a can be arbitrarily close to 2m from below. Let us first prove (4.16). We only need to consider r < 1/2, because otherwise it is obvious. By the Poincaré inequality and (4.20) with some a > m, u 2 dxdt Cr 2m D Q d m u 2 L (Q r ) Q r Crd4m Dd m u 2 L (Q 1/2 ) r Cr d4m u 2 Wp 1,2m (Q 1/2 ) Crd4m u 2 dxdt. Q 1 The last inequality is due to the L p estimate in (4.19) with 1/2 and 1 in place of 1 and 2 respectively. Next, we prove (4.17) and only consider the

17 SCHAUDER ESTIMATES 17 case when r < 1/2 as in proving (4.16). By the Poincaré inequality and (4.20) with some a > m, we have Q r D l d u 2 dxdt r 2(m l) Q r D m d u 2 dxdt Cr 2(m l)d2m sup Dd m u 2 Cr 2(m l)d2m u 2 Q Wp 1,2m (Q 1/2 ) 1/2 Cr 2(m l)d2m u 2 dxdt Cr 2(m l)d2m Dd l u 2 dxdt. Q 1 Q 1 As for (4.18), let v = u xm d m! (Dm d u) Q so that v also satisfies the system as 1 u. Moreover, Dd mv = Dm d u (Dm d u) Q. Then, by (4.20), (4.19), and the 1 Poincaré inequality, for a = mγ (m,m1) and r < 1/2, Q r Dd m v (Dm d v) Q 2 dxdt Cr 2(a m)d2m v 2 a r 2m,a,Q r Cr 2γd2m v 2 dxdt Cr 2γd2m Dd m v 2 dxdt. Q 1 Q 1 The lemma is proved. Similar to Lemma 4.3, we consider and u satisfies (4.14). u t ( 1) m L 0 u = D α f α in Q 2R, (4.21) Lemma 4.5. Assume that u Cloc (Rd1 ) satisfies (4.21) and (4.14), and f α C a for α = m where a (0,1), f α L for α < m. Then for any 0 < r < R 1,γ (0,1) and X 0 {x d = 0} Q R, there exists a constant C = C(d,m,n,λ,K,γ) such that Q r (X 0 ) D m x u 2 dxdt C( r R )4md Q R (X 0) Dx m u 2 dxdtc [f α ] 2 a,q R 2a2md 2R C f α 2 L (Q 2R )R22md, (4.22)

18 18 H. DONG AND H. ZHANG Q r (X 0 ) D m d u (Dm d u) Q r (X 0 ) 2 dxdt C( r R )2γd2m C [f α ] 2 a,q 2R Q R (X 0) D Q d l Dm l x u 2 dxdt r (X 0 ) C( r R )2(m l)d2m C [f α ] 2 where 0 < l < m. a,q 2R D m d u (Dm d u) Q R (X 0) 2 dxdt R 2a2md C Q R (X 0) Dd l Dm l x u 2 dxdt R 2a2md C f α 2 L (Q 2R )R22md, (4.23) f α 2 L (Q 2R )R22md, (4.24) Proof. We follow the proof of Lemma 4.3. Let w be the weak solution of the following system w t ( 1) m L 0 w = D α f α in Q 2R, w = 0, D d w = 0,, D m 1 d w = 0 on p Q 2R. Then the Poincaré inequality and the method in the proof of Lemma 4.3 yield D m w 2 dxdt Q 2R C [f α ] 2 a,q 2R R 2a2md C f α 2 L (Q 2R )R22md. We again use the mollification argument as in Lemma 4.3 so that w is smooth. Let v = u w. Then v satisfies all the conditions in Lemma 4.4, and so does D m x v. Thus we obtain (4.16) with D m x v in place of u. Combining the estimates of v and w, we obtain (4.22). The proofs of (4.23) and (4.24) are similar. Now we are ready to prove the Hölder estimates similar to the interior case. From (4.22), (4.23), and (4.24), taking γ > a and using Lemma 3.2, we get for any X 0 {x d = 0} Q R and r (0,R), Dx m u 2 Dd m u (Dm d u) Q r (X 0 ) 2 Dd l Dm l x u 2 dxdt I Q r (X 0 )

19 where I := C(R (2a2md) SCHAUDER ESTIMATES 19 Q 2R D m u 2 dxdtf 2 )r 2a2md F = [f α ] f a,q α 2R L (Q ). 2R This estimate, together with the interior estimates and Lemma 3.1, implies that [D m u] a C(R (2a2md) D m u 2 dxdtf 2 ) m,a,Q R Q 2R Similar to what we did in proving (4.1), we can apply the argument of freezing the coefficients to deal with the variable coefficients case with lowerorder terms. Here we just state the conclusion: Proposition 4.6. Assume that R 1 and u Cloc (Rd1 ) satisfies u t ( 1) m Lu = D α f α in Q 2R, u = 0, D d u = 0,, D m 1 d u = 0 on {x d = 0} Q 2R. Suppose that L and f α satisfy the conditions in Theorem 2.1. Then there exists a constant C = C(d,n,m,λ,K, A αβ a,r,a) such that u 1 2 a [D m u] 2m,Q a C( D β u R 2m,a,Q R L (Q 2R ) F). β m 5. Proof of Theorem 2.1 Before proving the main theorem, let us show a technical lemma. Lemma 5.1. Assume that a (0,1), u C loc (Rd1 ), Ω is an open set in R d and Ω C 2, T (0, ], then for any ε > 0 u L ((0,T) Ω) ε[u] a 2m,a,(0,T) Ω C(ε) u L2 ((0,T) Ω), where C(ε) is a constant depending on ε. Proof. Since Ω C 2, thereexistsar 0 > 0suchthatforanyX 0 (0,T) Ω, we can choose an arbitrarily small cylinder Q r (s,y) Ω with r (0,r 0 ) so that X 0 Q r (s,y) and u(x 0 ) u(x 0 ) u(t,x)dxdt u(t,x)dxdt Q r(s,y) Q r(s,y) Cr a [u] a 2m,a Cr d2m 2 u L2. Therefore, u L ε[u] a 2m,a C(ε) u L2. The lemma is proved.

20 20 H. DONG AND H. ZHANG Next, we show a global a priori estimate. Proposition 5.2. Let L, Ω, and f α satisfy the conditions in Theorem 2.1. Assume that u C ((0,T) Ω) and satisfies the following system u t ( 1) m Lu = D α f α in (0,T) Ω, u = 0, Du = 0,, D m 1 u = 0 on [0,T) Ω, u = 0 on {0} Ω. Then (2.2) holds with g = 0. (5.1) Proof. Using the standard arguments of partition of unity and flattening the boundary, we combine the interior, boundary estimates, and the estimate in t variable for the divergence form systems to get u am 2m,(0,T) Ω [Dm u] a 2m,a,(0,T) Ω C( β m D β u L ((0,T) Ω) F). By the interpolation inequalities in Hölder spaces, for instance, see [14, Section 8.8], u am 2m,(0,T) Ω [Dm u] a 2m,a,(0,T) Ω C( u L ((0,T) Ω) F). Applying Lemma 5.1 and the interpolation inequalities again, we get u L ((0,T) Ω) C(ε) u L2 ((0,T) Ω)ε[D m u] a 2m,a,(0,T) Ω ε u am 2m,(0,T) Ω. Upon taking ε sufficiently small, we arrive at (2.2). In order to implement the method of continuity, we need the right-hand sideof (2.2)tobeindependentofuandthisleadsustoconsiderthefollowing system, u t ( 1) m Luκu = D α f α in (0,T) Ω, u = 0, Du = 0,, D m 1 u = 0 on [0,T) Ω, u = 0 on {0} Ω, (5.2) where κ is a large constant to be specified later. We rewrite the system as u t ( 1) m L h uκu = D α f α ( 1) m (L h L)u, (5.3)

21 SCHAUDER ESTIMATES 21 where L h u denotes the sum of the highest-order terms. Then multiply u to both sides of (5.3) and integrate over (0,T) Ω. Thus D m u 2 L 2 ((0,T) Ω) κ u 2 L 2 ((0,T) Ω) ( 1) α f α D α udxdt, β m, β <m (0,T) Ω ( 1) m α 1 (0,T) Ω (0,T) Ω A αβ D β ud α udxdt A αβ D β ud α udxdt. (5.4) Take apoint x 0 Ω. Dueto thehomogeneous Dirichlet boundarycondition, if α = m, the factor f α in the firstintegral on the right-hand side above can be replaced by f α (t,x) f α (t,x 0 ). We use the Cauchy Schwarz inequality, Young s inequality, and the interpolation inequality to bound the right-hand side by ( C(n,m,d,λ,K, A αβ a,ε) f α f α (t,x 0 )1 2 dxdt (0,T) Ω (0,T) Ω ) u 2 dxdt ε D m u 2 L 2 ((0,T) Ω). After taking ε sufficiently small to absorb the term ε D m u 2 L 2 ((0,T) Ω) to the left-hand side of (5.4) and choosing κ sufficiently large, we reach u L2 ((0,T) Ω) CF, (5.5) where C = C(d,n,m,λ,K, A αβ a,ω). Combining (5.5) with (2.2), we get the following lemma. Lemma 5.3. Assume that L, Ω, and f α satisfy all the conditions in Theorem 2.1 and u C ((0,T) Ω) satisfies (5.2) with a sufficiently large constant κ depending on n,m,d,λ,k and A αβ a, then where C = C(d,m,n,λ,K, A αβ a,ω,a). Now we are ready to prove Theorem 2.1. u am,am,(0,t) Ω CF, (5.6) 2m Proof of Theorem 2.1. By considering u g instead of u, without loss of generality, we can assume g = 0. It remains to prove the solvability. Let κ be the constant from the previous lemma. Since u satisfies (5.1), the function v := e κt u satisfies v t ( 1) m Lv κv = e κt D α f α.

22 22 H. DONG AND H. ZHANG We then reduce the problem to the solvability of v. By Lemma 5.3, (5.6) holds for v. Consider the following equation v t ( 1) m (slv (1 s) m v)κv = e κt D α f α with the same initial and boundary conditions, where the parameter s [0,1]. Itisknownthatwhens = 0thereisauniquesolutioninC am 2m,am ([0,T) Ω). Then by the method of continuity and the a priori estimate (5.6), we find a solution when s = 1 which gives u. Remark 5.4. Actually the condition f α L for α < m can be relaxed. With slight modification in our proof, f α can be in some Morrey spaces. 6. Estimate for non-divergence type systems In this section, we deal with the non-divergence type systems Interior estimate. First, we consider the interior estimates for the non-divergence form system u t ( 1) m L 0 u = f in Q 2R. For this system, the interior estimates are simple consequences of the corresponding estimates for divergence form systems. We differentiate the system m times with respect to x to get D m u t ( 1) m L 0 D m u = D m f. Thanks to the estimates of the divergence type systems, by (4.8), [D 2m u] a 2m,a,Q R C( D2m u L (Q 2R ) [f] a,q 2R ). Next we deal with the general non-divergence form systems. The coefficients are all in C a. Following exactly the same idea as in handling the divergence form, we first consider L = α = β =m Aαβ D α D β. Fix y B R and let L 0y = A αβ (t,y)d α D β. We rewrite the system as follows α = β =m u t ( 1) m L 0y u = f ( 1) m (L 0y L)u. Set ζ to be an infinitely differentiable function in R d1 such that ζ = 1 in Q 2R, ζ = 0 outside ( (4R) 2m,(4R) 2m ) B 4R. (6.1) FromTheorem2.6in[5], fort = (4R) 2m,letwbetheuniqueW 1,2m 2 -solution of the following system w t ( 1) m L 0y w = ζ(f (f) B4R (t))ζ( 1) m (L 0y L)u

23 SCHAUDER ESTIMATES 23 in ( T,0) R d with the zero initial condition at t = T. It is easy to obtain the estimate for w, D 2m w 2 L 2 (Q 4R ) C f (f) B 4R (t) 2 L 2 (Q 4R ) C (L 0y L)u 2 L 2 (Q 4R ) CR 2a2md [f] 2 a,q 4R CR 2a2md [A αβ ] 2 a D2m u 2 L (Q 4R ). α = β =m We apply the mollification argument so that w is smooth. Let v = u w whichisalso asmoothfunction. Moreover sinceζ = 1inQ 2R, D 2m v satisfies (4.3). Therefore (4.4) holds for D 2m v. Combining the estimate of w and D 2m v, similar to the proof of Proposition 4.1, we get [D 2m u] a 2m,a,Q R C( D2m u L (Q 4R ) [f] a,q 4R ). (6.2) For the operator L with lower-order terms, we rewrite the systems as follows u t ( 1) m L h u = f ( 1) m1 A αβ D α D β u. where L h = α = β =m Aαβ D α D β. Let f = f ( 1) m1 An easy calculation shows that [ f] a,q 4R C ( [f] a,q 4R α β <2m α β <2m [A αβ ] a Dα D β u L (Q 4R )) ). α β <2m A αβ D α D β u. ( A αβ L [D α D β u] a,q 4R From (6.2) with f in place of f, the following estimate holds [D 2m u] a 2m,a,Q C( D γ u R L (Q 4R ) [f] a,q 4R ), γ 2m where C = C(d,m,n,λ,K, A αβ a,r,a). Because it follows immediately that u t = ( 1) m1 Luf, u t a,q R [D 2m u] a 2m,a,Q R C( γ 2m D γ u L (Q 4R ) f a,q 4R ). Implementing a standard interpolation argument, for instance, see [14, Section 8.8] and Lemma 5.1 of [8], we are able to prove the firstpart of Theorem 2.2.

24 24 H. DONG AND H. ZHANG 6.2. Boundary estimate. In the non-divergence case, better estimates than the ones in Lemma 4.4 are necessary. For (4.18), we actually can estimate more normal derivatives up to 2m 1-th order. In order to show this, let us prove the following lemma. A similar result can be found in [6]. Lemma 6.1. Assume that u Cloc (Rd1 ) satisfies (4.13) and (4.14). Let Q(x) be a vector-valued polynomial of order m 1 and P(x) = x m d Q(x). Suppose that (D k Dd m P(x)) Q = (D k Dd m u(t,x)) R Q R, where 0 k m 1. Let v = u P(x), then there exists a constant C = C(d,m,n,λ) such that for any 0 k m 1. D k D m d v L 2 (Q R ) CRm 1 k D m 1 D m d v L 2 (Q R ) Proof. Without loss of generality, let us assume R = 1. Choose ξ(y) C0 (B 1 ) with unit integral. Then let 0 g k (t) = ξ(y)d k Dd m v(t,y)dy, c k = g k (t)dt. B 1 By the Poincaré inequality and Hölder s inequality, the following estimate holds D k Dd m v(t,x) g k(t) 2 dx B 1 = B 1 C C B 1 B 1 Because (D k D m d v) Q 1 B 1 B 1 1 (D k D m d v(t,x) Dk D m d v(t,y))ξ(y)dy 2 dx (D k D m d v(t,x) Dk D m d v(t,y)) 2 dydx D k1 D m d v(t,y) 2 dy. = 0, we have D k D m d v L 2 (Q 1 ) Dk D m d v c k L2 (Q 1 ) D k D m d v g k(t) L2 (Q 1 ) g k(t) c k L2 (Q 1 ) C D k1 D m d v L 2 (Q 1 ) C tg k L2 ( 1,0). (6.3) Since v satisfies the same system as u, by the definition of g k and integrating by parts t g k (t) = ξ(y)d k Dd m tv(t,y)dy B 1 = ( 1) mk1 B 1 D k ξ(y)l 0 D m d v(t,y)dy.

25 SCHAUDER ESTIMATES 25 We integrate by parts again leaving m 1 derivatives on Dd m all the others onto ξ to get t g k (t) C D m 1 Dd m v dy. B 1 v and moving Therefore, t g k (t) 2 C D m 1 Dd m v(t, ) 2 L 2 (B 1 ). (6.4) Combining (6.3) and (6.4), we prove the lemma by induction. Next, we prove an estimate for D 2m 1 u. Lemma 6.2. Assume that u Cloc (Rd1 ) and satisfies (4.13) and (4.14). Then for any 0 < r R < and γ (0,1), there exists a constant C = C(d,n,m,λ,γ) such that for any X 0 {x d = 0} Q R Q r (X 0 ) D 2m 1 u (D 2m 1 u) Q r (X 0 ) 2 dxdt C( r R )2γd2m Q R (X 0) D 2m 1 u (D 2m 1 u) Q R (X 0) 2 dxdt. Proof. By scaling and translation of the coordinates, without loss of generality, we can assume R = 1 and X 0 = (0,0). Moreover, we only need to consider the case when r < 1/2 for the same reason as in the proof of Lemmas 4.2 and 4.4. Recall that in Lemma 4.4 we have u a 2m,a,Q 1 C u W 1,2m p (Q 1 ) C u L 2 (Q 2 ), where a < 2m can be arbitrarily close to 2m. Set a= γ 2m 1 (2m 1,2m). Applying the inequality above with 1/2 and 1 in place of 1 and 2, and the Poincaré inequality, we have Q r D 2m 1 u (D 2m 1 u) Q r 2 dxdt Cr 2γd2m u 2 a 2m,a,Q r Cr2γd2m u 2 W 1,2m p (Q 1/2 ) Cr 2γd2m u 2 L 2 (Q 1 ) Cr2γd2m Dd m u 2 L 2 (Q 1 ). (6.5) Let v be the function in Lemma 6.1. Note that v also satisfies (4.13) and (4.14). Then the inequality (6.5) holds with v in place of u. Since P(x) is a

26 26 H. DONG AND H. ZHANG polynomial of degree 2m 1, using Lemma 6.1, we have D 2m 1 u (D 2m 1 u) Q Q r 2 dxdt = D 2m 1 v (D 2m 1 v) Q r Q r 2 dxdt r Cr 2γd2m Dd m v 2 dxdt Cr 2γd2m D m 1 Dd m v 2 dxdt = Cr 2γd2m Cr 2γd2m Q 1 Q 1 Q 1 The lemma is proved. Q 1 D m 1 Dd m u (Dm 1 Dd m u) Q 2 dxdt 1 D 2m 1 u (D 2m 1 u) Q 2 dxdt. 1 Let us turn to non-homogeneous systems. Consider that u C loc (Rd1 ) and satisfies { ut ( 1) m L 0 u = f in Q 4R, u = 0, D d u = 0,,D m 1 d u = 0 on Q 4R {x d = 0}. We shall show the estimates of all derivatives up to order 2m with the exception of Dd 2mu. Let ζ bethe function in (6.1). For T = (4R)2m, consider the following system w t ( 1) m L 0 w = ζ(f (f) B (t)) in ( T,0) R d, 4R w = 0, D d w = 0,, D m 1 d w = 0 on ( T,0) {x d = 0} w = 0 on {t = T}. WeknowfromTheorem2.6in[5]that thesystem abovehasauniquew 1,2m 2 - solution which satisfies the following estimate D 2m w 2 L 2 (Q 4R ) C f (f) B (t) 2 4R L 2 (Q 4R ) CR2a2md [f] 2. (6.6) a,q 4R We again apply the mollification argument so that w is smooth. Observe that v := u w satisfies { vt ( 1) m L 0 v = (1 ζ)f ζ(f) B (t) in Q 4R 2R, v = 0, D d v = 0,,D m 1 d v = 0 on {x d = 0} Q 2R. Since ζ = 1 in Q 2R, ṽ := D x v satisfies {ṽt ( 1) m L 0 ṽ = 0 in Q 2R, ṽ = 0, D d ṽ = 0,,D m 1 d ṽ = 0 on {x d = 0} Q 2R. Thanks to Lemma 6.2, for any X 0 {x d = 0} Q R we get D x D 2m 1 v (D x D 2m 1 v) Q r (X 0 ) 2 dxdt Q r (X 0 ) C( r R )2γd2m Q R (X 0) (6.7) D x D 2m 1 v (D x D 2m 1 v) Q R (X 0) 2 dxdt. (6.8)

27 SCHAUDER ESTIMATES 27 Therefore, by the triangle inequality, D x D 2m 1 u (D x D 2m 1 u) Q r (X 0 ) 2 dxdt Q r (X 0 ) C( r R )2γd2m C Q 2R Q R (X 0) D 2m w 2 dxdt. D x D 2m 1 u (D x D 2m 1 u) Q R (X 0) 2 dxdt Combining (6.6), we obtain D x D 2m 1 u (D x D 2m 1 u) Q r (X 0 ) 2 dxdt Q r (X 0 ) C( r R )2γd2m Q R (X 0) CR 2ad2m [f] 2. a,q 4R D x D 2m 1 u (D x D 2m 1 u) Q R (X 0) 2 dxdt Taking γ > a, from Lemmas 3.2, 3.1, and the corresponding interior estimates, we obtain the estimate for the Hölder semi-norm, [D x D 2m 1 u] a C( D 2m u 2m,a,Q R L2 (Q 4R ) [f] ). a,q 4R Let us turn to the case that the coefficients are functions of both t and x, and the method of freezing the coefficients can still be implemented in this condition. We first consider the case that L = α = β =m Aαβ D α D β. Fix y B R and rewrite the system as follows u t ( 1) m L 0y u = f ( 1) m (L 0y L)u. We use the same ζ as in (6.1), T = (4R) 2m, and w is the unique W 1,2m 2 - solution of the system w t ( 1) m L 0y w = ζ(f (f) B (t))ζ( 1) m (L 0y L)u 4R in ( T,0) R d with the Dirichlet boundary conditions w = 0, D d w = 0,, D m 1 d w = 0 on ( T,0) {x d = 0} and the zero initial condition at t = T. It is easy to obtain the estimate for w, D 2m w 2 L 2 (Q 4R ) C f (f) B (t) 2 4R L 2 (Q 4R ) C (L 0y L)u 2 L 2 (Q 4R ) CR 2a2md [f] 2 CR 2a2md [A αβ ] 2 a,q a D 2m u 2 4R L (Q ). 4R α = β =m Let v = u w. Then ṽ := D x v satisfies (6.7), which implies (6.8). Consequently, taking γ > a, we get the following inequality [D x D 2m 1 u] a C([f] D 2m u 2m,a,Q R a,q 4R L (Q 4R )). (6.9)

28 28 H. DONG AND H. ZHANG If the operator L has lower-order terms, we can move all the lower-order terms to the right-hand side regarded as a part of f. Hence we proved the second part of Theorem 2.2, i.e., [D x D 2m 1 u] a C([f] 2m,a,Q R a,q D γ u 4R L (Q )). 4R γ 2m The proof of Theorem 2.2 is completed. Remark 6.3. We need more regularity assumptions to get the estimate of [Dd 2mu] a. For example, if assuming A αβ,f C a 2m,a,Q 2m,a, then for R 1, R by using a similar method one can show that [Dd 2m u] a [u 2m,a,Q t ] a C( f R 2m,a,Q a u R 2m,a,Q 8R L2 (Q 8R )). (6.10) We give a sketched proof. First, let us assume that the coefficients are constants and u is smooth and satisfies (4.13) and (4.14). In this case, we can differentiate the equation with respect to t which means u t satisfies the same equation. From Lemma 4.4, we know that for any r < R 1 and X 0 {x d = 0} Q R, there exists a constant C such that u t 2 dxdt C( r R )d4m u t 2 dxdt. (6.11) Q r (X 0 ) Q R (X 0) We are now ready to show (6.10). For simplicity we only consider L which consists of highest-order terms. We use the same ζ as in (6.1), T = (4R) 2m, and fix X 0 Q R. Let w be the W1,2m 2 -solution of the system w t ( 1) m L X0 w = ζ(f (f) Q (t))ζ( 1) m (L X0 L)u 4R in ( T,0) R d with the Dirichlet boundary conditions w = 0, D dw = 0,, D m 1 d w = 0 on ( T,0) {x d = 0} and the zero initial condition at t = T, where L X0 = A αβ (X 0 )D α D β. Let v = u w, which satisfies { vt ( 1) m L X0 v = (f) Q in Q 4R 2R, v = 0, D d v = 0,,D m 1 d v = 0 on {x d = 0} Q 2R. Then (6.11) holds with v in place of u. By the triangle inequality, for any X 1 {x d = 0} Q R and r (0,R), u t 2 dxdt C( r R )d4m u t 2 dxdtc w t 2 dxdt. Q r (X 1 ) Q R (X 1) Q R (X 1) (6.12)

29 SCHAUDER ESTIMATES 29 By Theorem 2.6 in [5], w t 2 L 2 (Q 4R ) C f (f) Q 2 4R L 2 (Q 4R ) (L X 0 L)u 2 L 2 (Q 4R ) CR 2md2a( [f] 2 a ) [A αβ ] 2 2m,a,Q a 4R 2m,a D2m u 2 L (Q 4R ). α = β =m (6.13) We combine (6.12), (6.13), and the corresponding interior estimates to obtain a Hölder estimate, [u t ] a C ( [f] 2m,a,Q a [A αβ ] R 2m,a,Q a 4R 2m,a D 2m u L (Q 4R ) u ) t L2 (Q. 4R ) α = β =m It remains to estimate Dd 2m u. Let γ = (0,0,,m), A γγ D γ D γ u = ( 1) m1 u t A αβ D α D β u( 1) m f. (α,β) (γ,γ) (6.14) Since A γγ is positive definite, it follows easily that [Dd 2m u] a C ( [u t ] 2m,a,Q a D α ) u R 2m,a,Q a [f] R 2m,a,Q a. R 2m,a,Q R α (0,,2m) α =2m Plugging (6.9) and (6.14) into the inequality above and using the interpolation inequalities, we immediately prove (6.10). References [1] Agmon S., Douglis A., Nirenberg L.: Estimates near the boundary for solutions of elliptic partial differential equations satisfying general boundary conditions, I, Comm. Pure Appl. Math. 12, (1959); II, ibid., 17, (1964). [2] Brandt A.: Interior Schauder estimates for parabolic differential-(or difference-) equations via the maximum principle, Israel J. Math. 7 (1969), [3] Boccia S.: Schauder estimates for solutions of high-order parabolic systems, preprint. [4] Campanato S.: Equazioni paraboliche del secondo ordine e spazi L 2,θ (Ω,δ), Ann. Mat. Pura. Appl. 73 (1966), [5] Dong H., Kim D.: On the L p-solvability of higher order parabolic and elliptic systems with BMO coefficients, Arch. Rational Mech. Anal. 199 (2011), no. 3, [6] Dong H., Kim D.: Higher order elliptic system in Sobolev spaces with partial BMO coefficients, J. Funct. Anal. 261 (2011), no. 11, [7] Friedman A.: Partial differential equations of parabolic type, Prentice-Hall, Englewood cliffs, N.J., [8] Giaquinta M.: Introduction to regularity theory for nonlinear elliptic systems, Lectures in Mathematics ETH Zürich. Birkhäuser Verlag, Basel, [9] Lieberman G.: Intermediate Schauder Theory For second order parabolic equations. IV. Time irregularity and regularity, Differential and Integral Equations 5 (1992), no. 6, [10] Lieberman G.: Second order parabolic differential equations, Word Scientific Publishing Co. Pte. Ltd, Singapore-New Jersey-London-Hong Kong, [11] Lorenzi L.: Optimal Schauder estimates for parabolic problems with data measurable with respect to time, SIAM J. Math. Anal. 32 (2000), no. 3,

30 30 H. DONG AND H. ZHANG [12] Lunardi, A.: Analytic semigroups and optimal regularity in parabolic problems, Progress in Nonlinear Differential Equations and their Applications, 16. Birkhäuser Verlag, Basel, [13] Knerr B.: Parabolic interior Schauder estimates by the maximum principle, Arch. Rational Math. Anal. 75 (1980), [14] Krylov N. V.: Lectures on elliptic and parabolic equations in Hölder spaces, American Mathematical Society, Providence, RI, [15] Krylov N. V., Priola E.: Elliptic and parabolic second-order PDEs with growing coefficients, Comm. Partial Differential Equations 35 (2010), no. 1, [16] Schlag W.: Schauder and L p estimates for parabolic system via Campanato spaces, Comm. Partial Differential Equations 21 (1996), no. 7-8, [17] Simon L.: Schauder estimates by scaling, Calc. Var. Partial Differential Equations 5 (1997), no. 5, [18] Sinestrari E., von Wahl W.: On the solutions of the first boundary value problem for the linear parabolic equations, Proc. Roy. Soc. Edinburgh Sect. A 108 (1988), no. 3-4, (H. Dong) Division of Applied Mathematics, Brown University, 182 George Street, Providence, RI 02912, USA address: Hongjie Dong@brown.edu (H. Zhang) Division of Applied Mathematics, Brown University, 182 George Street, Providence, RI 02912, USA address: Hong Zhang@brown.edu