One method which goes back to works of E. Heinz and F. Sauvigny will be used in the paper at hand: Using isothermal parameters, a minimal surface $X$  and its unit normal mapping $N$  in ${\mathbb{R}}^3$  satisfy the nonlinear elliptic systems $$△X=0,△N=-2KWN.$$  Here, $K$  denotes the Gaussian curvature and $W$  the surface area element of the immersion $X.$  If one can establish certain bounds on the size of $X$  and $N,$  the gradient estimates in [3] ensure upper bounds for the first and second derivatives of $X$  in terms of given a priori-data. A Harnack-type inequality for pseudo-holomorphic functions gives lower bounds for $W.$  Altogether, this gives estimates for the principle curvatures of the immersion, namely in such a way that if the surface growths to become complete, the principle curvatures vanish identically.

This method generalizes to various variational problems if the differential systems satisfy a structure condition of the form $$|△Z|\le \text{const}|\nabla Z{|}^2.$$  For example, the mean curvature system $△X=2HWN$  and the above system for $N$  satisfy this assumption. Some problems arise in higher codimension.

In 1964, Osserman [5] proved that a complete two-dimensional minimal surface in ${\mathbb{R}}^n$  is a plane if all of its normal vectors make a certain positive angle with a fixed axis in space. The assumption that all normals make an acute angle with some fixed direction makes the difference to the Bernstein theorem in ${\mathbb{R}}^3$  where it is sufficient to assume that the normal of the immersion omits a certain neighbourhood of some point of the unit sphere.

Osserman's method is based essentially on results of complex analysis due to the fact that the components of a minimal surface are harmonic functions.

Although it is rather difficult to realize Osserman's condition, it is an important fact that any condition of this type is necessary for curvature estimates.

In this paper we want to examine how to apply the methods developed in [3] and [6] to establish certain a priori-estimates for the principle curvatures of graphs with prescribed mean curvature.

In the first part, we concentrate on some geometric basics of two-dimensional immersion in ${\mathbb{R}}^n$  which are motivated from [2] , and we will derive nonlinear elliptic differential systems to describe the surfaces analytically. The second part contains the mentioned estimates.

1. Surfaces in ${\mathbb{R}}^n$ 

The essential principles of the following are motivated from [1] , [2] , [5] , and [10] .

Let $$B:=\left\{\right(u,v)\in {\mathbb{R}}^2:u^2+v^2<1\}$$  denote the open unit disc, $\overline{B}\subset {\mathbb{R}}^2$  its topological closure. For integer $n\ge 3$  we consider immersions

$$\begin{array}{c}X=X(u,v)\in C^{2+\alpha }(B,{\mathbb{R}}^n)\cap C^0(\overline{B},{\mathbb{R}}^n),\alpha \in (0,1),\end{array}$$

$$\begin{array}{c}X(u,v)=\left(x^1\right(u,v),x^2(u,v),\dots ,x^n(u,v\left)\right),\end{array}$$

with the property $$\text{rank}\partial X(u,v)\equiv \text{rank}\left(\begin{array}{cc}{x}_u^1(u,v)& x_v^1(u,v)\\ ...& ...\\ x_u^n(u,v)& x_v^n(u,v)\end{array}\right)=2\text{for all}(u,v)\in B.$$  This means that $X_u=X_u(u,v)$  and $X_v=X_v(u,v)$  are linearly independent at every point $(u,v)\in B$  and, therefore, span the two-dimensional tangent plane $${\mathcal{T}}_X\left(w\right):=\text{Span}\left\{X_u\right(w),X_v(w\left)\right\}\text{at}w\in B.$$  There exists a normal space ${\mathcal{N}}_X\left(w\right)$  spanned by $n-2$  linearly independent vectors $N_1,N_2,\dots ,N_{n-2}$  which can be assumed to satisfy the orthonormality relations $$N_{\Sigma }\cdot N_{\Theta }^t={\delta }_{\Sigma \Theta }:=\{\begin{array}{c}1\text{if}\Sigma =\Theta \\ 0\text{if}\Sigma \ne \Theta \end{array}\text{for all}\Sigma ,\Theta \in \{1,2,\dots ,n-2\}.$$  Here, by the upper $t$  we denote the transposed vector.

The choice of the normal vectors $N_{\Sigma }\left(w\right),$  $\Sigma =1,\dots ,n-2,$  is not unique.

Example 1 : The unit vectors

$$\begin{array}{ccc}{N}_1& :=& \frac{1}{\sqrt{1+|\nabla {\phi }_1{|}^2}}(-{\phi }_{1,x},-{\phi }_{1,y},1,0,\dots ,0),\end{array}$$

$$\begin{array}{ccc}{N}_2& :=& \frac{1}{\sqrt{1+|\nabla {\phi }_2{|}^2}}(-{\phi }_{2,x},-{\phi }_{2,y},0,1,\dots ,0),\dots \end{array}$$

are normal to the graph $X(x,y)=(x,y,{\phi }_1(x,y),\dots ,{\phi }_n(x,y\left)\right).$ 

Assumption: There exists an orthonormal basis $N_1\left(w\right),N_2\left(w\right),\dots ,N_{n-2}\left(w\right)$  of the normal space ${\mathcal{N}}_X\left(w\right)$  such that $N_{\Sigma }\in C^{1+\alpha }(B,S^{n-1}).$  In the following we will use such regular bases.

1..1 Differential geometric foundations

We consider surfaces $X=X(u,v)$  immersed into the Euclidean space ${\mathbb{R}}^n.$  Inserting $X=X(u,v)$  into the element $d{s}^2:=(d{x}^1{)}^2+\dots +(d{x}^n{)}^2$  yields $$d{s}^2=h_{11}d{u}^2+2h_{12}dudv+h_{22}d{v}^2=h_{ij}d{u}^{i}d{u}^j$$  (note the summation convention in the last term) with the first fundamental form $$I:=(h_{ij}{)}_{i,j=1,2}\equiv (X_{u^i}\cdot X_{u^j}^t{)}_{i,j=1,2}\subset {\mathbb{R}}^{2\times 2}.$$  Using conformal parameters $(u,v)\in B,$  the first fundamental form appears as $$h_{11}=W=h_{22},h_{12}=0\text{in}B$$  with the surface area element $W:=\sqrt{h_{11}{h}_{22}-h_{12}^2}.$  The second fundamental form $I{I}_{\Sigma }=(L_{\Sigma ,ij}{)}_{i,j=1,2}$  w.r.t. the unit normal $N_{\Sigma }$  is defined as $$L_{\Sigma ,ij}:=-X_{u^i}\cdot N_{\Sigma ,u^j}^t=X_{u^i{u}^j}\cdot N_{\Sigma }^t,i,j=1,2.$$  Consider the form $$(L_{\Sigma ,i}^k{)}_{i,k=1,2}:=(L_{\Sigma ,ij}{h}^{jk}{)}_{i,k=1,2}\in {\mathbb{R}}^{2\times 2}$$  with the coefficients $h^{ij}$  of the inverse of the first fundamental form, i.e. $h_{ij}{h}^{jk}={\delta }_i^k$  with the Kronecker symbol ${\delta }_i^k.$ 

Definition 1 : The mean curvature $H_{\Sigma }$  in direction $N_{\Sigma }$  is defined as $$H_{\Sigma }:=\frac{1}{2}\text{trace}(L_{\Sigma ,i}^k{)}_{i,k=1,2}=\frac{L_{\Sigma ,11}{h}_{11}-2L_{\Sigma ,12}{h}_{12}+L_{\Sigma ,22}{h}_{22}}{2(h_{11}{h}_{22}-h_{12}^2)}.$$ 

Besides the mean curvature we introduce the Gaussian curvature w.r.t. to $N_{\Sigma }$  as $$K_{\Sigma }:=\text{det}(L_{\Sigma ,i}^k{)}_{i,k=1,2}=\frac{L_{\Sigma ,11}{L}_{\Sigma ,22}-L_{\Sigma ,12}^2}{h_{11}{h}_{22}-h_{12}^2}.$$  The principle curvatures ${\kappa }_{\Sigma ,1},$  ${\kappa }_{\Sigma ,2}$  w.r.t. $N_{\Sigma }$  are the eigenvalues of the form $(L_{\Sigma ,i}^k{)}_{i,k=1,2},$  that is $$H_{\Sigma }=\frac{{\kappa }_{\Sigma ,1}+{\kappa }_{\Sigma ,2}}{2},K_{\Sigma }={\kappa }_{\Sigma ,1}{\kappa }_{\Sigma ,2}.$$ 

1..2 The differential equations

Let $N_{\Sigma },$  $\Sigma =1,\dots ,n-2,$  be an orthonormal basis of the normal space. We rewrite the derivatives $N_{\Sigma ,u^i}$  and $X_{u^i{u}^j}$  for and $i,j=1,2$  in terms of the moving frame $X_u,X_v,N_1,\dots ,N_{n-2}.$  The following is discussed more detailed in [2] .

Proposition 1 : Let $N_{\Sigma },$  $\Sigma =1,\dots ,n-2,$  be an orthonormal basis of the normal space. Then $$N_{\Sigma ,u^i}=-L_{\Sigma ,ij}{h}^{jk}{X}_{u^k}^t+{\sigma }_{\Sigma ,i}^{\Theta }{N}_{\Theta },i=1,2,\Sigma =1,2,\dots ,n-2,$$  with the torsion coefficients $${\sigma }_{\Sigma ,i}^{\Theta }:=\{\begin{array}{cc}{N}_{\Sigma ,u^i}\cdot N_{\Theta }^t& \text{if}\Sigma \ne \Theta \\ 0& \text{if}\Sigma =\Theta \end{array},$$  as well as $$X_{u^i{u}^j}={\Gamma }_{ij}^k{X}_{u^k}+{\sum }_{\Sigma =1}^{n-2}{L}_{\Sigma ,ij}{N}_{\Sigma }$$  with the Christoffel symbols ${\Gamma }_{ij}^k:=\frac{1}{2}{h}^{kl}(h_{jl,i}+h_{li,j}-h_{ij,l}),$  $h_{ij,k}:=h_{ij,u^k}.$ 

1..3 An elliptic system for $X$ 

Let $(u,v)\in B$  be conformal parameters. Then, there hold ${\Gamma }_{11}^1+{\Gamma }_{22}^1=0,$  ${\Gamma }_{11}^2+{\Gamma }_{22}^2=0.$  Thus, the Gauss equations from Proposition 1 yield $$△X=({\Gamma }_{11}^1+{\Gamma }_{22}^1)X_u+({\Gamma }_{11}^2+{\Gamma }_{22}^2)X_v+{\sum }_{\Sigma =1}^{n-2}(L_{\Sigma ,11}+L_{\Sigma ,22})N_{\Sigma }=2{\sum }_{\Sigma =1}^{n-2}{H}_{\Sigma }W{N}_{\Sigma }.$$  This system is called the mean curvature system.

Let $\left|H_1\right|+\left|H_2\right|+\dots +\left|H_{n-2}\right|\le h_0$  with a real constant $h_0\in (0,+\infty ).$  We estimate $$|△X|\le 2h_0\left|X_u\right|\left|X_v\right|\le h_0|\nabla X{|}^2$$  due to $W=\sqrt{h_{11}{h}_{22}}=\left|X_u\right|\left|X_v\right|.$  Note the quadratic growth of the gradient on the right hand side.

Elliptic systems of this structure are the matter of investigations in [3] .

2. A priori-estimates for graphs in ${\mathbb{R}}^4$ 

In this part we consider immersions

$$\begin{array}{c}X=X(u,v)\in C^{2+\alpha }(B,{\mathbb{R}}^4)\cap C^0(\overline{B},{\mathbb{R}}^4),\alpha \in (0,1),\end{array}$$

$$\begin{array}{c}X(u,v)=\left(x^1\right(u,v),\dots ,x^4(u,v\left)\right).\end{array}$$

Let a function $\overline{\mathcal{\mathscr{H}}}\left(X\right):{\mathbb{R}}^4⟶{\mathbb{R}}^4$  be given. We define a scalar function $\mathcal{\mathscr{H}}(X,Z):{\mathbb{R}}^4\times S^3\to \mathbb{R}$  by $$\mathcal{\mathscr{H}}(X,Z):=\overline{\mathcal{\mathscr{H}}}\left(X\right)\cdot Z^t\text{for}X\in {\mathbb{R}}^4,Z\in S^3.$$  We now consider surfaces of prescribed mean curvature $\mathcal{\mathscr{H}}=\mathcal{\mathscr{H}}(X,Z)$  in the following sense.

Definition 2 : A surface $X=X(u,v)$  is called a (conformally parametrized) surface of prescribed mean curvature $\mathcal{\mathscr{H}}:{\mathbb{R}}^4\times S^3\to \mathbb{R}$  if it satisfies the system

$$\begin{array}{c}△X=2\mathcal{\mathscr{H}}(X,N_1)W{N}_1+2\mathcal{\mathscr{H}}(X,N_2)W{N}_2,\end{array}$$

$$\begin{array}{c}|X_u{|}^2=W=|X_v{|}^2,X_u\cdot X_v^t=0\text{in}B\end{array}$$

for an orthonormal basis $N_1,N_2$  of the normal space.

Remarks:

• 1. By the definition of $\mathcal{\mathscr{H}}$  , the above system is satisfied for all orthonormal bases of the normal space if it is satisfied for at least one orthonormal basis.
• 2. Using conformal parameters the geometric mean curvature $H_{\Sigma }\left(X\right)\equiv H(X,N_{\Sigma })$  of $X$  in direction of the unit normal $N_{\Sigma }$  can be calculated by $$H_{\Sigma }\left(X\right)=\frac{(h_{11}{X}_{uu}-2h_{12}{X}_{uv}+h_{22}{X}_{vv})\cdot N_{\Sigma }^t}{2(h_{11}{h}_{22}-h_{12}^2)}=\frac{△X\cdot N_{\Sigma }^t}{2W}.$$  Thus, by section 1.3, the geometric mean curvature $H_{\Sigma }\left(X\right)$  is equal to $\mathcal{\mathscr{H}}(X,N_{\Sigma })$  if the normal basis is continuously differentiable.

Our main result is the following.

Theorem 1 :

Assumptions: Let $X\in C^{2+\alpha }(B,{\mathbb{R}}^4)\cap C^0(\overline{B},{\mathbb{R}}^4),$  $\alpha \in (0,1),$  be a surface of prescribed mean curvature $\mathcal{\mathscr{H}}(X,Z)\in C^0({\mathbb{R}}^4\times S^3,\mathbb{R})$  with the following properties:

• (A1) The immersion $X=X(u,v)$  is a positively oriented conformal reparametrization of a graph of the form $$X(x,y)=(x,y,\phi (x,y),\psi (x,y\left)\right),(x,y)\in {\overline{B}}_R,$$  where $B_R:=\left\{\right(x,y)\in R^2:x^2+y^2<R^2\},$  and $\phi ,\psi \in C^{2+\alpha }({\overline{B}}_R,\mathbb{R}).$ 
• (A2) The mean curvature $\mathcal{\mathscr{H}}=\mathcal{\mathscr{H}}(X,Z)$  satisfies $$|\mathcal{\mathscr{H}}(X,Z\left)\right|\le \frac{h_0}{2}\text{for all}X\in {\mathbb{R}}^4,Z\in S^3,$$  as well as $$|\mathcal{\mathscr{H}}(X_1,Z_1)-\mathcal{\mathscr{H}}(X_2,Z_2\left)\right|\le h_1|X_1-X_2{|}^{\alpha }+h_2|Z_1-Z_2|\text{for all}{X}_1,X_2\in {\mathbb{R}}^4,Z_1,Z_2\in S^3$$  with constants $h_0,h_1,h_2\in [0,+\infty ).$ 
• (A3) The surface represents a geodesic disc of radius $r>0$  and with center $X(0,0)=(0,0,0,0).$ 
• (A4) With a real constant $d_0>0,$  the Dirichlet energy can be estimated by $$\mathcal{D}\left[X\right]:=\int {\int }_B|\nabla X(u,v){|}^{2}dudv\le d_0{r}^2.$$ 
• (A5) At every point $w\in B,$  each normal vector of the immersion makes an angle of at least $\omega >0$  with the $x_1$  -axis.

Statement: Then, for orthonormal vectors $N_1$  and $N_2$  spanning the normal space at the point $X(0,0),$  there exists a constant $\Theta =\Theta (h_{0}r,h_1{r}^{1+\alpha },h_{2}r,d_0,sin\omega ,\alpha )\in (0,+\infty )$  such that holds $${\kappa }_{\Sigma ,1}(0,0{)}^2+{\kappa }_{\Sigma ,2}(0,0{)}^2\le \frac{1}{r^2}\left\{(h_{0}r{)}^2+\Theta \right\}$$  for the principle curvatures ${\kappa }_{\Sigma ,1}$  and ${\kappa }_{\Sigma ,2}$  in direction $N_{\Sigma }.$ 

Proof: The proof follows [6] , § 2. It holds

$$\begin{array}{ccc}{\kappa }_{\Sigma ,1}(0,0{)}^2+{\kappa }_{\Sigma ,2}(0,0{)}^2& =& 4H_{\Sigma }(0,0{)}^2-2K_{\Sigma }(0,0)\le h_0^2+2|K_{\Sigma }(0,0)|\end{array}$$

$$\begin{array}{ccc}& =& \frac{1}{r^2}\left\{(h_{0}r{)}^2+r^2|K_{\Sigma }(0,0)|\right\}\end{array}$$

for $\Sigma =1,2.$  The desired estimate follows from $$K_{\Sigma }(0,0)=\frac{(X_{uu}\cdot N_{\Sigma })(X_{vv}\cdot N_{\Sigma })-(X_{uv}\cdot N_{\Sigma }{)}^2}{W^2}{|}_{(0,0)}$$  w.r.t. a normal basis $\{N_1,N_2\}.$  Thus, our proof consists of two steps: First, we have to find a lower bound for the surface area element, and then an upper bound for the second derivatives of the immersion.

• 1. In the first part we will prove the estimate $$\frac{W\left(w\right)}{r^2}\ge C_1\text{for}w\in B_{\frac{1}{2}}(0,0)$$  with a constant $C_1=C_1(h_{0}r,d_0,sin\omega )>0$  .
  • 1.1 We fix a normal basis of the surface. Note that the vectors $e_3:=(0,0,1,0)$  and $e_4:=(0,0,0,1)$  are not in any tangent plane of the surface due to the fact that $$\frac{1}{\sqrt{1+|\nabla \phi {|}^2}}(-{\phi }_x,-{\phi }_y,1,0),\frac{1}{\sqrt{1+|\nabla \psi {|}^2}}(-{\psi }_x,-{\psi }_y,0,1)$$  are normal to the surface (in the corresponding surface point). We see that the inner product of the first vector by $e_3$  and the inner product of the second vector by $e_4$  do not vanish. Therefore, the projections $$N_1^{*}:=e_3-\frac{e_3\cdot X_u^t}{|X_u{|}^2}{X}_u-\frac{e_3\cdot X_v^t}{|X_v{|}^2}{X}_v,N_2^{*}:=e_4-\frac{e_4\cdot X_u^t}{|X_u{|}^2}{X}_u-\frac{e_4\cdot X_v^t}{|X_v{|}^2}{X}_v$$  are a basis of the normal space and can be transformed into an orthonormal basis $\{N_1,N_2\}$  of the normal space.
  • 1.2 In the first part of the proof we will work with the above orthonormal basis $\{N_1,N_2\}.$  Using the differential equation (see Remark 1 above) $$△X=2\mathcal{\mathscr{H}}(X,N_1)W{N}_1+2\mathcal{\mathscr{H}}(X,N_2)W{N}_2\text{in}B,$$  Assumption (A1) yields the estimate $$|△X(u,v\left)\right|\le 2h_0|\nabla X(u,v){|}^2\text{in}B.$$  The special structure of this differential inequality the quadratic growth in the gradient on the right hand side allows us to apply the methods of [3] .
    We cite two important consequences of our assumptions.
  • 1.3 From Assumption (A3) we conclude:
    Let $\Gamma \left(B\right)$  denote the set of all continuous and piecewise differentiable curves $\gamma :[0,1]\to \overline{B},$  such that $\gamma \left(0\right)=(0,0)$  and $\gamma \left(1\right)\in \partial B.$  Then, it holds $${inf}_{\gamma \in \Gamma \left(B\right)}{\int }_0^1\left|\frac{d}{dt}X\circ \gamma \left(t\right)\right|dt\ge r.$$  For a proof we refer the reader to [6] .
  • 1.4 Assumption (A5) enables us to interpolate the surface area element $W=W(u,v)$  in terms of $|\nabla x^1{|}^2:$  Namely, it holds $$|\nabla x^1{|}^2\ge W{sin}^2\omega \text{in}B.$$  The proof can be extracted from [5] , Lemma 1.1 (note that this result of Osserman makes only use of the conformal representation of the surface).
    Now, we estimate the surface area element. For this, we define several two-dimensional auxiliary functions and apply Heinz' results for elliptic systems in the plane (cp. [3] ).
  • 1.5 Let $F^{*}(u,v):=\left(x^1\right(u,v),x^2(u,v\left)\right):\overline{B}\to {\mathbb{R}}^2$  denote the plane mapping w.r.t.
    $X=X(u,v).$  Then, for all $w\in B$  there hold

    • (i) $|△F^{*}(w\left)\right|\le \frac{4h_0}{{sin}^2\omega }|\nabla F^{*}(w){|}^2$ 
    • (ii) $|\nabla X(w){|}^2\le \frac{2}{{sin}^2\omega }|\nabla F^{*}\left(w\right){|}^2.$ 

    For the proof we note that the second inequality follows from $W+W=|\nabla X{|}^2$  and from the Osserman inequality in 1.4. The differential inequality $|△X|\le 2h_0|\nabla X{|}^2=4h_{0}W$  (see 1.2) yields $$|△F^{*}|\le |△X|\le 4h_{0}W\le \frac{4h_0}{{sin}^2\omega }|\nabla x^1{|}^2\le \frac{4h_0}{{sin}^2\omega }|\nabla F^{*}{|}^2.$$ 
  • 1.6 Let $w_0\in B$  and $\nu \in (0,1)$  be given such that $B_{2\nu }\left(w_0\right):=\{w\in B:|w-w_0|<2\nu \}\subset B.$  We consider the mapping $$Y\left(w\right):=\frac{1}{r}\left\{X\right(w_0+2\nu w)-X(w_0\left)\right\},w\in \overline{B},$$  as well as the corresponding plane mapping $F\left(w\right):=\left(y^1\right(w),y^2(w\left)\right):\overline{B}\to {\mathbb{R}}^2.$  The immersion $Y$  satisfies $\left|Y_u\right(w){|}^2=\frac{4{\nu }^2}{r^2}W(w_0+2\nu w)=|Y_v\left(w\right){|}^2,$  $Y_u\left(w\right)\cdot Y_v(w{)}^t=0$  and $$|△Y(w\left)\right|\le 2\left(h_{0}r\right)|\nabla Y(w){|}^2\text{in}B.$$ 
  • 1.7 From inequality 1.5 and 1.6 we infer

    $$\begin{array}{ccc}|△F(w\left)\right|& \le & |△Y(w\left)\right|\le 2\left(h_{0}r\right)|\nabla Y(w){|}^2\end{array}$$

    $$\begin{array}{ccc}& =& \frac{16{\nu }^2\left(h_{0}r\right)}{r^2}W(w_0+2\nu w)=\frac{8{\nu }^2\left(h_{0}r\right)}{r^2}|\nabla X(w_0+2\nu w){|}^2\end{array}$$

    $$\begin{array}{ccc}& \le & \frac{16{\nu }^2\left(h_{0}r\right)}{r^2{sin}^2\omega }|\nabla F^{*}(w_0+2\nu w){|}^2=\frac{16{\nu }^2\left(h_{0}r\right)}{r^2{sin}^2\omega }\frac{r^2}{4{\nu }^2}|\nabla F\left(w\right){|}^2\end{array}$$

    $$\begin{array}{ccc}& \le & \frac{4\left(h_{0}r\right)}{{sin}^2\omega }|\nabla F(w){|}^2\end{array}$$

    for all $w\in B.$  Furthermore, by the assumption that the immersion represents a graph of the form $(x,y,\phi (x,y),\psi (x,y\left)\right),$  $(x,y)\in {\overline{B}}_R,$  we conclude that the plane mapping $F=F(u,v)$  is one-to-one and has positive Jacobian $J_F\left(w\right)>0$  for all $w\in B$  (see e.g.
    [9] , [7] ). Additionally, Assumption (A4) gives $$\mathcal{D}\left[F\right]\le \mathcal{D}\left[Y\right]\le \frac{1}{r^2}\mathcal{D}\left[X\right]\le d_0$$  for the Dirichlet integrals. We now apply [3] , Theorem 6, page 254, which gives the following inner gradient estimate: There is a constant $c_1=c_1(h_{0}r,d_0,sin\omega )\in (0,+\infty )$  such that $$|\nabla F(u,v\left)\right|\le c_1(h_{0}r,d_0,sin\omega )\text{for all}(u,v)\in B_{\frac{1}{2}}(0,0).$$ 
  • 1.8 From 1.7 we infer $\frac{4{\nu }^2}{r^2}W(w_0+2\nu w)\le \frac{1}{{sin}^2\omega }|\nabla F(w){|}^2$  for $w\in B.$  Thus, we arrive at the following estimate $$\frac{1}{r^2}W\left(w\right)\le \frac{1}{4{\nu }^2{sin}^2\omega }{c}_1(h_{0}r,d_0,sin\omega )=:c_2(h_{0}r,d_0,sin\omega ,\nu )$$  for all $w\in B_{\nu }\left(w_0\right)$  and for all $w_0\in B$  such that $B_{2\nu }\left(w_0\right)\subset B.$  For $w_0\in B_{\frac{1}{2}}(0,0)$  and $\nu \le \frac{1}{4}$  we get $$\frac{1}{r^2}W\left(w\right)\le c_2(h_{0}r,d_0,sin\omega ,\nu )\text{for all}w\in B_{\frac{1}{2}}(0,0).$$  This estimate will be used in the second part of the proof.
  • 1.9 The properties $J_F\left(w\right)>0$  in $B$  and $\mathcal{D}\left[F\right]\le d_0,$  as well as the structure of the differential inequality for $|△F(w\left)\right|$  (see 1.7) make [3] , Lemma 17, page 255 applicable:
    There exists a constant $c_3=c_3(h_{0}r,d_0,sin\omega )\in (0,+\infty )$  such that $$|\nabla F(w){|}^2\le c_3(h_{0}r,d_0,sin\omega \left)\right|\nabla F(0,0){|}^{\frac{2}{5}}\text{for all}w\in B_{\frac{1}{2}}(0,0).$$  It follows, that (cp. 1.6)

    $$\begin{array}{ccc}\frac{4{\nu }^2}{r^2}W(w_0+2\nu w)& \le & \frac{1}{{sin}^2\omega }|\nabla F(w){|}^2\le \frac{c_3(h_{0}r,d_0,sin\omega )}{{sin}^2\omega }|\nabla F(0,0){|}^{\frac{2}{5}}\end{array}$$

    $$\begin{array}{ccc}& \le & \frac{c_3(h_{0}r,d_0,sin\omega )}{{sin}^2\omega }|\nabla Y(0,0){|}^{\frac{2}{5}}=\frac{c_3(h_{0}r,d_0,sin\omega )}{{sin}^2\omega }{\left[\frac{8{\nu }^2}{r^2}W\left(w_0\right)\right]}^{\frac{1}{5}}.\end{array}$$

    Rearranging yields the important inequality of Harnack type $${\left[\frac{W\left(w_0\right)}{r^2}\right]}^{\frac{1}{5}}\ge \frac{8^{\frac{17}{15}}{\nu }^{\frac{8}{5}}{sin}^2\omega }{c_3(h_{0}r,d_0,sin\omega )}\frac{W(w_0+2\nu w)}{r^2}\text{for all}w\in B_{\frac{1}{2}}(0,0),$$  or, which is equivalent, $$c_4(h_{0}r,d_0,sin\omega ,\nu ){\left[\frac{W\left(w\right)}{r^2}\right]}^5\le \frac{W\left(w_0\right)}{r^2}\text{for all}w\in B_{\nu }\left(w_0\right)$$  with the constant $c_4(h_{0}r,d_0,sin\omega ,\nu ):=\frac{8^{\frac{17}{3}}{\nu }^8{sin}^{10}\omega }{c_3(h_{0}r,d_0,sin\omega {)}^5}\in (0,+\infty ).$ 
• 1.10 Assumption (A4) also ensures that we can estimate the surface area element in at least one point: There exists a $w^{*}\in B_{1-{\nu }_0}(0,0),$  ${\nu }_0:=min(e^{-4\pi d_0},\frac{1}{2})$  such that $$\frac{W\left(w^{*}\right)}{r^2}\ge \frac{1}{4(1-e^{-4\pi d_0})}=:c_5\left(d_0\right)>0.$$  The constant ${\nu }_0$  arises from an application of the Courant-Lebesgue lemma (see [6] ).
• 1.11 We now show an estimate of the surface area element:
  We set $\nu :=\frac{1}{2}{\nu }_0\in (0,\frac{1}{4}]$  and choose an integer $n=n\left(\nu \right)\in \mathbb{N}$  such that $1-2\nu \le \frac{n}{2}\nu \le 1-\nu .$  For an arbitrary $w_0\in B_{1-{\nu }_0}(0,0)$  we define the following points $$w_j:=\frac{j}{n}{w}^{*}+\frac{n-j}{n}{w}_0\text{for}j=0,\dots ,n.$$  Together with 1.10 we have $\left|w_j\right|\le \frac{j}{n}\left|w^{*}\right|+\frac{n-j}{n}\left|w_0\right|<1-{\nu }_0,$  and therefore $B_{2\nu }\left(w_j\right)=B_{{\nu }_0}\left(w_j\right)\subset B.$  Furthermore, we have $$|w_{j+1}-w_j|=|\frac{1}{n}{w}^{*}-\frac{1}{n}{w}_0|\le \frac{1}{n}|w^{*}-w_0|\le \frac{2(1-{\nu }_0)}{n}\le \nu .$$  This implies $w_{j+1}\in B_{\nu }\left(w_j\right)\text{for}j=0,\dots ,n-1.$  We apply the Harnack inequality from 1.9 and obtain $$\frac{W\left(w_0\right)}{r^2}\ge c_4{\left[\frac{W\left(w_1\right)}{r^2}\right]}^5\ge c_4^{1+5}{\left[\frac{W\left(w_2\right)}{r^2}\right]}^{5^2}\ge \dots \ge c_4^{1+5+5^2+\dots +5^{n-1}}{\left[\frac{W\left(w_n\right)}{r^2}\right]}^{5^n}.$$  Recalling $w_n=w^{*}$  , 1.10 gives $$\frac{W\left(w_0\right)}{r^2}\ge c_4^{1+5+5^2+\dots +5^{n-1}}{c}_5(d_0{)}^{5^n}=:C_1(h_{0}r,d_0,sin\omega )>0$$  for all $w_0\in B_{1-{\nu }_0}(0,0)$  . From ${\nu }_0\le \frac{1}{2}$  we now conclude $$\frac{W\left(w\right)}{r^2}\ge C_1(h_{0}r,d_0,sin\omega )\text{for all}w\in B_{\frac{1}{2}}(0,0).$$  This ends the first part of the proof.
• 2. In the second part of the proof we have to estimate the second derivatives of $X$  using the differential equation $$△X=2\mathcal{\mathscr{H}}(X,N_1)W{N}_1+2\mathcal{\mathscr{H}}(X,N_2)W{N}_2.$$  To do so, we have to give Hölder estimates of the right hand side of this equation. Especially, we have to construct an orthonormal basis $\{N_1,N_2\}$  of the normal space at each point, whose Hölder norm can be estimated.
• 2.1 We define the auxiliary function $$Z(u,v)=\frac{1}{r}\left\{X\right(u,v)-X(0,0\left)\right\}=\frac{1}{r}X(u,v),(u,v)\in \overline{B}.$$  Denoting by $W_Z$  the surface area element of $Z,$  we have $|Z_u{|}^2=W_Z=|Z_v{|}^2$  and $Z_u\cdot Z_v^t=0$  in $B.$  It holds $r^2{W}_Z=W_X$  with $W_X:=|X_u{|}^2=|X_v{|}^2.$  We calculate

  $$\begin{array}{ccc}△Z& =& \frac{2}{r}\mathcal{\mathscr{H}}(X,N_1)W_X{N}_1+\frac{2}{r}\mathcal{\mathscr{H}}(X,N_2)W_X{N}_2\end{array}$$

  $$\begin{array}{ccc}& =& 2r\mathcal{\mathscr{H}}(rZ,N_1)W_Z{N}_1+2r\mathcal{\mathscr{H}}(rZ,N_2)W_Z{N}_2.\end{array}$$

• 2.2 Due to 1.8 we have the estimate $$|△Z(w\left)\right|\le 4\left(r{h}_0\right)c_2(h_{0}r,d_0,sin\omega )\text{for all}w\in B_{\frac{1}{2}}(0,0).$$  Furthermore, we get $$\left|Z\right(u,v\left)\right|=\left|Z\right(u,v)-Z(0,0\left)\right|\le \frac{1}{2}|\nabla Z(\tilde{w}\left)\right|\le \sqrt{c_2(h_{0}r,d_0,sin\omega )}\text{in}{B}_{\frac{1}{2}}(0,0).$$  Now, by potential theoretic estimates, there exists a constant $c_6(h_{0}r,d_0,sin\omega ,\alpha )$  such that the Hölder estimate $$\left|Z_{u^i}\right(w_1)-Z_{u^i}(w_2\left)\right|\le c_6(h_{0}r,d_0,sin\omega ,\alpha )|w_1-w_2{|}^{\alpha }\text{for}{w}_1,w_2\in B_{\frac{1}{4}}(0,0),i=1,2$$  holds true. Therefore, for the surface area element it holds $$\left|W_Z\right(w_1)-W_Z(w_2\left)\right|\le c_7(h_{0}r,d_0,sin\omega ,\alpha )|w_1-w_2{|}^{\alpha }\text{for all}{w}_1,w_2\in B_{\frac{1}{4}}(0,0)$$  with the constant $c_7:=4\sqrt{c_2}{c}_6.$ 
• 2.3 Using the mean value theorem we have the Lipschitz estimate $$\left|Z\right(w_1)-Z(w_2\left)\right|\le 2\sqrt{c_2(h_{0}r,d_0,sin\alpha )}|w_1-w_2|\text{for}{w}_1,w_2\in B_{\frac{1}{2}}(0,0).$$  In a certain neighborhood of the origin, we now construct an orthonormal basis $\{N_1,N_2\}$  of the normal space whose Hölder norm can be estimated.
  First we choose vectors ${\overline{N}}_1,{\overline{N}}_2\in {\mathbb{R}}^4$  such that ${\overline{N}}_i\cdot Z_{u^j}(0,0{)}^t=0,$  ${\overline{N}}_i\cdot {\overline{N}}_j^t={\delta }_{ij},$  i.e. $\{{\overline{N}}_1,{\overline{N}}_2\}$  form a orthonormal basis of the normal space at the point $Z(0,0)$  . Due to the Gram-Schmidt orthonormalization, for $k=1,2$  we define vectors $$N_k^{*}\left(w\right):={\overline{N}}_k-\frac{{\overline{N}}_k\cdot Z_u\left(w\right)}{\left|Z_u\right(w){|}^2}{Z}_u\left(w\right)-\frac{{\overline{N}}_k\cdot Z_v\left(w\right)}{\left|Z_v\right(w){|}^2}{Z}_v\left(w\right)\text{in}B.$$  These vectors both belong to the normal space at $Z\left(w\right)$  but they may not be linearly independent. At first, we determine a ${\nu }_1={\nu }_1(h_{0}r,d_0,sin\omega ,\alpha )>0$  such that $$\left|N_k^{*}\right(w){|}^2=1-\frac{[{\overline{N}}_k\cdot Z_u(w){]}^2}{W_Z\left(w\right)}-\frac{[{\overline{N}}_k\cdot Z_v(w){]}^2}{W_Z\left(w\right)}\ge \frac{1}{2}\text{in}{B}_{{\nu }_1}(0,0).$$  This is possible because, first, we calculate

  $$\begin{array}{ccc}|{\overline{N}}_k\cdot Z_{u^{\ell }}(w){|}^2& =& |{\overline{N}}_k\cdot \{Z_{u^{\ell }}\left(w\right)-Z_{u^{\ell }}(0,0)\}{|}^2\le |Z_{u^{\ell }}\left(w\right)-Z_{u^{\ell }}(0,0){|}^2\end{array}$$

  $$\begin{array}{ccc}& \le & c_6(h_{0}r,d_0,sin\omega ,\alpha {)}^2|w{|}^{2\alpha }\end{array}$$

  by 2.2, and next, by 1.11, we have the lower bound $$W_Z\left(w\right)\ge C_1(h_{0}r,d_0,sin\omega )\text{in}{B}_{\frac{1}{2}}(0,0).$$ 
• 2.4 We remark that the vectors $N_k^{*}\left(w\right),$  $k=1,2,$  are Hölder continuous in $B_{{\nu }_1}(0,0)$  and the Hölder estimate $$\left|N_k^{*}\right(w_1)-N_k^{*}(w_2\left)\right|\le c_8(h_{0}r,d_0,sin\omega ,\alpha )|w_1-w_2{|}^{\alpha },w_1,w_2\in B_{{\nu }_1}(0,0),$$  holds true with a constant $c_8(h_{0}r,d_0,sin\omega ,\alpha )$  coming from the Hölder estimate for $Z_{u_j}$  and the lower bound of $W_Z.$  Now, for $k=1,2$  we define $${\tilde{N}}_k\left(w\right):=\frac{N_k^{*}\left(w\right)}{\left|N_k^{*}\right(w\left)\right|}\text{in}{B}_{{\nu }_1}(0,0).$$  These vectors are well defined because of $\left|N_k^{*}\right(w){|}^2\ge \frac{1}{2}$  in $B_{{\nu }_1}(0,0),$  but they are not orthogonal. Note that $$N_1^{*}\cdot N_2^{*}=-\frac{({\overline{N}}_1\cdot Z_u^t)({\overline{N}}_2\cdot Z_u^t)}{W_Z}-\frac{({\overline{N}}_1\cdot Z_v^t)({\overline{N}}_2\cdot Z_v^t)}{W_Z},$$  and, therefore, due to 2.3 $$|{\tilde{N}}_1\cdot {\tilde{N}}_2|=\frac{|N_1^{*}\cdot N_2^{*}|}{\left|N_1^{*}\right|\left|N_2^{*}\right|}\le \frac{2}{C_1}\left(\right|{\overline{N}}_1\cdot Z_u^t\left|\right|{\overline{N}}_2\cdot Z_u^t|+|{\overline{N}}_1\cdot Z_v^t\left|\right|{\overline{N}}_2\cdot Z_v^t\left|\right)\le \frac{4c_6^2}{C_1}|w{|}^{2\alpha }.$$  Thus we can find a ${\nu }_2={\nu }_2(h_{0}r,d_0,sin\omega ,\alpha )$  with $0<{\nu }_2\le {\nu }_1$  such that $$|{\tilde{N}}_1\cdot {\tilde{N}}_2|\le \frac{1}{2}\text{in}{B}_{{\nu }_2}(0,0).$$  Now we define the vectors $$N_1\left(w\right):={\tilde{N}}_1\left(w\right),N_2\left(w\right):=\frac{{\tilde{N}}_2\left(w\right)-(N_1\cdot {\tilde{N}}_2)N_1}{\sqrt{1-(N_1\cdot {\tilde{N}}_2{)}^2}}\text{in}{B}_{{\nu }_2}(0,0).$$  $N_2$  is well defined because for its denominator we have $1-(N_1\cdot {\tilde{N}}_2{)}^2=1-({\tilde{N}}_1\cdot {\tilde{N}}_2{)}^2\ge \frac{3}{4}.$  The vectors $N_1\left(w\right),N_2\left(w\right)$  are an orthonormal basis of the normal space at each point $Z\left(w\right)$  for $w\in B_{{\nu }_2}(0,0)$  . Furthermore, the Hölder estimate $$\left|N_k\right(w_1)-N_k(w_2\left)\right|\le c_9(h_{0}r,d_0,sin\omega ,\alpha )|w_1-w_2{|}^{\alpha }\text{for}{w}_1,w_2\in B_{{\nu }_2}(0,0)$$  holds true with a constant $c_9(h_{0}r,d_0,sin\omega ,\alpha )$  which can be calculated directly using the Hölder estimate for $N_k^{*}$  .
• 2.5 We now use the differential equation (recall Remark 1 to Definition 2) $$△Z=2r\mathcal{\mathscr{H}}(rZ,N_1)W_Z{N}_1+2r\mathcal{\mathscr{H}}(rZ,N_2)W_Z{N}_2\text{in}{B}_{{\nu }_2}(0,0).$$  We already showed the estimate $|△Z(w\left)\right|\le 4\left(h_{0}r\right)c_2(h_{0}r,d_0,sin\omega )$  in $B_{{\nu }_2}(0,0).$  Using the assumptions on $\mathcal{\mathscr{H}}$  we obtain the Hölder estimate

  $$\begin{array}{c}|\mathcal{\mathscr{H}}(rZ\left(w_1\right),N_k\left(w_1\right))-\mathcal{\mathscr{H}}(rZ\left(w_2\right),N_k\left(w_2\right)\left)\right|\end{array}$$

  $$\begin{array}{c}\le h_1{r}^{\alpha }\left|Z\right(w_1)-Z(w_2){|}^{\alpha }+h_2|N_k\left(w_1\right)-N_k\left(w_2\right)|\end{array}$$

  $$\begin{array}{c}\le h_1{2}^{\alpha }{r}^{\alpha }{c}_2^{\frac{\alpha }{2}}|w_1-w_2{|}^{\alpha }+h_2{c}_9|w_1-w_2{|}^{\alpha }.\end{array}$$

  Thus we can find a constant $c_{10}=c_{10}(h_{0}r,h_1{r}^{1+\alpha },h_{2}r,d_0,sin\omega ,\alpha )$  such that $$|△Z(w_1)-△Z(w_2\left)\right|\le c_{10}|w_1-w_2{|}^{\alpha }\text{for}{w}_1,w_2\in B_{{\nu }_2}(0,0).$$  Setting ${\nu }_3:=\frac{1}{2}{\nu }_2$  , the interior Schauder estimates give a constant $C_2\in (0,+\infty )$  such that $$\left|Z_{uu}\right(w\left)\right|,\left|Z_{uv}\right(w\left)\right|,\left|Z_{vv}\right(w\left)\right|\le C_2(h_{0}r,h_1{r}^{1+\alpha },h_{2}r,d_0,sin\omega ,\alpha )\text{in}{B}_{{\nu }_3}(0,0).$$  From the beginning of the proof we recall $${\kappa }_{1,\Sigma }(0,0{)}^2+{\kappa }_{2,\Sigma }(0,0{)}^2\le \frac{1}{r^2}\left\{(h_{0}r{)}^2+\frac{\left|Z_{uu}\right(0,0\left)\right|\left|Z_{vv}\right(0,0\left)\right|+\left|Z_{uv}\right(0,0){|}^2}{W_Z(0,0{)}^2}\right\}.$$  Setting $\Theta (h_{0}r,h_1{r}^{1+\alpha },h_{2}r,d_0,sin\omega ,\alpha ):=\frac{2C_2^2}{C_1^2}$  we arrive at $${\kappa }_{1,\Sigma }(0,0{)}^2+{\kappa }_{2,\Sigma }(0,0{)}^2\le \frac{1}{r^2}\left\{(h_{0}r{)}^2+\Theta \right\}.$$  This proves the statement. $\square$ 

Remarks:

• 1. For the introduction of conformal parameters into the graph we refer to [9] and [7] . In the theorem, we can replace the domain $B_R$  of the graph by a $C^2$  -Jordan domain $\Omega$  using Riemann's mapping theorem (see e.g. [8] ).
• 2. Assumption $\left(A5\right)$  is basically needed to establish the inequality $$|△F^{*}(w\left)\right|\le \frac{4h_0}{{sin}^2\omega }|\nabla F^{*}(w){|}^2$$  for the plane mapping $F^{*}$  . This is the difference to estimates for surfaces in ${\mathbb{R}}^3,$  where $|△F^{*}(w\left)\right|\le 4h_0|\nabla F^{*}(w){|}^2$  holds due to the conformal representation (note, that in this case $|X_u^3{|}^2\le |X_v^1{|}^2+|X_v^2{|}^2$  ).
• 3. The proof does not make use of the special dimension $n=4$  . Thus, it can be carried over to the more general case of dimension $n\ge 4$  .

In the case of minimal surfaces, i.e. $\mathcal{\mathscr{H}}\equiv 0$  , the proof yields a result of Bernstein type since the a priori constant $\Theta$  does not depend on $r$  any more. For $r\to \infty$  we obtain the

Corollary 1 : A complete minimal graph $X(x,y)=(x,y,\phi (x,y),\psi (x,y\left)\right)$  which satisfies assumptions (A4) and (A5) is an affine plane.

Our results are motivated from [5] , where curvature estimates and a Bernstein-type result for minimal immersions in ${\mathbb{R}}^n$  are proved.

Literatur

1. Blaschke, W.; Leichtweiss, K.: Elementare Differentialgeometrie. Grundlehren der mathematischen Wissenschaften 1, Springer, 1973.
2. Brauner, H.: Differentialgeometrie. Friedr. Vieweg & Sohn Verlagsgesellschaft mbH, Braunschweig, 1981.
3. Heinz, E.: On certain nonlinear differential equations and univalent mappings. Journal d'Analyse Math. 5, 197-272, 1956.
4. Lawson, H.B.: Lectures on minimal submanifolds. Vol. I, Mathematics Lecture Series 9, Publish or Perish, Inc., 1980.
5. Osserman, R.: Global properties of minimal surfaces in ${\mathbb{E}}^3$  and ${\mathbb{E}}^n.$  Ann. Math. 80, 340-364, 1964.
6. Sauvigny, F.: A-priori-Abschätzungen der Hauptkrümmungen für Immersionen vom Mittleren-Krümmungs-Typ mittels Uniformisierung und Sätze vom Bernstein-Typ. Habilitationsschrift, Göttingen 1988.
7. Sauvigny, F.: Introduction of isothermal parameters into a Riemannian metric by the continuity method. Analysis 19, No. 3, 235-243, 1999.
8. Sauvigny, F.: Partielle Differentialgleichungen der Geometrie und der Physik. Springer-Verlag, 2004.
9. Schulz, F.: Regularity theory for quasilinear elliptic systems and Monge-Ampère equations in two dimensions. Lecture Notes in Mathematics 1445, Springer-Verlag, 1990.
10. Weyl, H.: Zur Infinitesimalgeometrie: $p$  -dimensionale Flächen im $n$  -dimensionalen Raum. Math. Z. 12, 154-160, 1922.

Matthias Bergner, Steffen Fröhlich Technische Universität Darmstadt Fachbereich Mathematik, AG 4 Schloßgartenstraße 7D-64289 Darmstadt Germany e-mail: bergner@mathematik.tu-darmstadt.de sfroehlich@mathematik.tu-darmstadt.de