Extension Theory and Kreĭn-type Resolvent Formulas for Nonsmooth Boundary Value Problems

## Abstract

The theory of selfadjoint extensions of symmetric operators, and more generally the theory of extensions of dual pairs, was implemented some years ago for boundary value problems for elliptic operators on smooth bounded domains. Recently, the questions have been taken up again for nonsmooth domains, with results first on -domains for symmetric or smooth second-order operators, and next on quasi-convex Lipschitz domains for the selfadjoint realizations of the Laplacian. In the present work we show that pseudodifferential methods can be used to obtain a full characterization, including Kreĭn resolvent formulas, of the realizations of nonselfadjoint second-order operators on domains; more precisely, we treat domains with -smoothness and operators with -coefficients, for suitable and . The advantage of the pseudodifferential boundary operator calculus is that the operators are represented by a principal part and a lower-order remainder, leading to regularity results; in particular we analyze resolvents, Poisson solution operators and Dirichlet-to-Neumann operators in this way, also in Sobolev spaces of negative order. Some unbounded domains are allowed.

Key words: Extension theory; Krein resolvent formula; elliptic boundary value problems; pseudodifferential boundary operators; symbol smoothing; M-functions; nonsmooth domains; nonsmooth coefficients
MSC (2000): 35J25, 35P05, 35S15, 46E35, 47A10, 47A20, 47G30

## 1 Introduction

The systematic theory of selfadjoint extensions of a symmetric operator in a Hilbert space , or more generally, adjoint pairs of extensions of a given dual pair of operators in , has its origin in fundamental works of Kreĭn [42], Vishik [59] and Birman [18]. There have been several lines of development since then. For one thing, there are the early works of Grubb [29][31] completing and extending the theories and giving an implementation for results for boundary value problems for elliptic PDEs. Another line has been the development by, among others, Kochubei [40], Gorbachuk–Gorbachuk [28], Derkach–Malamud [23], Malamud–Mogilevskii [47], where the tendency has been to incorporate the problems into studies of relations (generalizing operators), with applications to (operator valued) ODEs; keywords in this development are boundary triples, Weyl-Titchmarsh -functions. More recently this has been applied to PDEs (e.g., Amrein-Pearson [9], Behrndt and coauthors [11, 12, 13, 14, 16], Brown-Marletta-Naboko-Wood [21], Kopachevskiĭ-Kreĭn [41], Malamud [46], Ryzhov [53]). Further references are given in Brown-Grubb-Wood [20], where a connection between the two lines of development is worked out.

One of the interesting aims is to establish Kreĭn resolvent formulas, linking the resolvent of a general operator with the resolvent of a fixed reference operator by expressing the difference in terms of operators connected to boundary conditions, encoding spectral information.

In the applications to elliptic PDEs, Kreĭn-type resolvent formulas are by now well-established in the case of operators with smooth coefficients on smooth domains, but there remain challenging questions about the validity in nonsmooth cases, and their applications.

One difficulty in implementing the extension theory in nonsmooth cases lies in the fact that one needs mapping properties of direct and inverse operators not only in the most usual Sobolev spaces, but also in spaces of low order, even of negative order over the boundary. Another difficulty is to arrive at a theory where ellipticity considerations are still applicable, in the way that the operators are defined from principal symbols plus lower-order error terms. This is important for regularity questions, as well as for questions of spectral estimates.

Gesztesy and Mitrea have addressed the extension problem for the Laplacian on Lipschitz domains, showing Kreĭn-type resolvent formulas in [24, 25, 26] involving Robin problems under the hypothesis that the boundary is of Hölder class . More recently, they have described the selfadjoint realizations of the Laplacian in [27] (based on the abstract theory of [29]), under a more general hypothesis of quasi-convexity, which includes convex domains and necessitates nonstandard boundary value spaces. Posilicano and Raimondi gave in [51] an analysis of selfadjoint realizations of second-order problems on -domains. Grubb treated nonselfadjoint realizations on -domains in [34], including Neumann-type boundary conditions

 χu=Cγ0u, (1.1)

with a differential operator of order 1, where the other mentioned works mainly treat cases (1.1) with of order or nonlocal. ([34] can be considered as a pilot project for the present paper.) It should also be mentioned that Behrndt and Micheler [15] recently have shown how a parametrization of the selfadjoint realizations of the Laplace operator on a Lipschitz domain can be obtained by use of the theory of quasi-boundary triples due to Behrndt et al. (cf. e.g. [12]). Compared with our results less regularity of the boundary is needed in the analysis. But the results are restricted to the Laplacian, while in the following we work with general second order elliptic operators. Moreover, in order to deal with Lipschitz boundary, where the usual results on elliptic regularity might fail, Behrndt and Micheler work with suitable more abstractly defined function spaces. In the case that the boundary is of class for some , their function spaces coincide with the classical ones, which we use in the following, cf. [15, Theorem 4.10] and the discussion below.

Our aim in this paper is to set up a construction of general extensions and resolvents that works in Sobolev spaces when the regularity of is in a scale of function spaces larger than , the coefficients of the elliptic operator in another larger scale, yet allowing the use of pseudodifferential calculi that can take ellipticity of boundary conditions into account and give precise information on the principal parts of the operators. We here choose to work with operators having coefficients in scales of Sobolev spaces and their generalizations to Besov and Bessel-potential spaces, since this allows rather precise multiplication properties, and convenient trace mapping results; then Hölder space properties can be read off using the well-known embedding theorems. The resulting hypothesis on is that it can be parametrized by functions in the Besov space for some . We note that this assumption is equivalent to , where is the regularity number of the Besov space , which is the scaling exponent of the highest order parts of the norms under dilations of functions. It is the most relevant number for Sobolev embeddings, estimates of nonlinearities and applications to nonlinear partial differential equations. We note that (locally) is inbetween and for and any . But the regularity number of is the same as the one of and can be much smaller than .

The theory of pseudodifferential boundary value problems (originating in Boutet de Monvel [19] and further developed e.g. in the book of Grubb [33]; introductory material is given in [35]) is well-established for operators with -coefficients on domains. It has been extended to nonsmooth cases by Abels [2], along the lines of the extension of pseudodifferential operators on open sets in Kumano-Go and Nagase [44], Marshall [49], Taylor [55], [56]. These results have been applied to studies of the Stokes operator in Abels [3] and Abels and Terasawa [5], which in particular imply optimal regularity results for the instationary Stokes system, cf. Abels [4]. For applications to quasi-linear differential equations and free boundary value problems non-smooth coefficients are essential, cf. e.g. Abels [1] and Abels and Terasawa [6]. The present paper builds on [2] and ideas of [5] and develops additional material.

Our final results will be formulated for operators acting between Sobolev spaces, but along the way we also need -based variants with for the operator- and domain-coefficients. Here the integral exponent will be called when we describe the domain and its boundary , and when we describe the given partial differential operators and boundary operators and their rules of calculus. There is then an optimal choice of how to link and , together with the dimension and the smoothness parameters of the spaces where the operators act; this is expressed in Assumption 2.18.

The results in the paper have been applied in Grubb [36] to show spectral asymptotic estimates for the boundary term in Kreĭn formulas established here.

We originally intended to include -order operators with , but the coefficients in Green’s formula needed an extra, lengthy development of symbol classes that made us postpone this to a future publication.

Plan of the paper. In Section 2, we recall the facts on Besov and Bessel-potential function spaces that we shall need, define the domains with boundary in these smoothness classes, and establish a useful diffeomorphism property. Nonsmooth pseudodifferential operators are recalled, with mapping- and composition-properties, and Green’s formula for second-order nonsmooth elliptic operators on appropriate nonsmooth domains is established. The Appendix gives further information on pseudodifferential boundary operators (dbo’s) with nonsmooth coefficients, extending some results of [2] to classes. Section 3 recalls the abstract extension theory of [29], [31], [20]. In Section 4 we use the dbo calculus to construct the resolvent and Poisson solution operator for the Dirichlet problem in the nonsmooth situation, by localization and parameter-dependent estimates. The construction shows that the principal part of the resolvent belongs to the class of non-smooth pseudodifferential boundary operators, which is essential for the subsequent analysis. Section 5 gives an extension of Green’s formula to low-order spaces, and provides an analysis of and the associated Dirichlet-to-Neumann operator , needed for the interpretation of the abstract theory. In particular it is shown that the operators coincide with operators of the pseudodifferential calculi up to lower order operators, which is one of the central results of the paper. Finally, the interpretation is worked out in Section 6, leading to a full validity of the characterization of the closed realizations of in terms of boundary conditions, and including Kreĭn-type resolvent formulas for all closed realizations . Section 7 gives a special analysis of the Neumann-type boundary conditions (1.1) entering in the theory, showing in particular that regularity of solutions holds when is elliptic.

## 2 Basics on function spaces and operators on nonsmooth domains

### 2.1 Function spaces on nonsmooth domains

For convenience we here recall the definitions and properties of function spaces that will be used throughout this paper. Proofs can be found e.g. in Triebel [57] and Bergh and Löfström [17]. All spaces are Banach spaces, some -based spaces are also Hilbert spaces.

The usual multi-index notation for differential operators with , , and , , will be employed.

For the spaces defined over , the Fourier transform is used to define operators such as (also called ), for suitable functions . In particular, with , stands for . denotes the Schwartz space of smooth, rapidly decreasing functions and its dual space, the space of tempered distributions.

Function spaces. The Bessel potential space in of order is defined for by

 Hsp(Rn)={f∈S′(Rn):⟨Dx⟩sf∈Lp(Rn)},

normed by . For , a non-negative integer, equals the space of -functions with derivatives up to order in , also denoted . In the case , we omit the lower index and simply write instead of . We denote the sesquilinear duality pairing of with by (linear in , conjugate linear in ).

To describe the regularity, both of domains and of operator-coefficients, we shall also need Besov spaces , where . These are defined by , where

 ∥f∥Bsp,q(Rn) = (∞∑j=02sjq∥φj(Dx)f∥qLp(Rn))1qif q<∞, ∥f∥Bsp,∞(Rn) = supj∈N02sj∥φj(Dx)f∥Lp(Rn).

Here, , , is a partition of unity on such that and if , chosen such that for all , .

The parameter indicates the smoothness of the functions. The second parameter is called the integration exponent. The third parameter is called the summation exponent; it measures smoothness on a finer scale than , which can be seen by the following simple relations:

 Bsp,1(Rn)↪Bsp,q1(Rn) ↪Bsp,q2(Rn)↪Bsp,∞(Rn) if 1≤q1≤q2≤∞, (2.1) Bs+εp,∞(Rn) ↪Bsp,1(Rn), (2.2)

where , , and are arbitrary. (The sign indicates continuous embedding.) The embeddings follow directly from the definition and the fact that if . Here, is the space of sequences such that in the case and if , provided with the hereby defined norm.

We recall that for and , equals the Sobolev-Slobodetskiĭ space , whereas for , it is that equals . (In the following, the - and -notation will be used for clarity; these scales of spaces have the best interpolation properties.) In the case , all three spaces coincide, for general :

 Hs(Rn)=Hs2(Rn)=Bs2,2(Rn)=Ws2(Rn). (2.3)

The spaces , also denoted when (Hölder-Zygmund spaces), play a special role. For , can be identified with the Hölder space , defined for , and , and also denoted when . For , there are sharp inclusions

 Ckb(Rn)↪Ck−1,1(Rn)↪Ck(Rn);

here, is the usual space of bounded continuous functions with bounded continuous derivatives up to order .

At this point, let us recall some interpolation results: Denoting the real and complex interpolation functors by and , respectively, we have that if with , , and , , then

 (Bs0p,q0(Rn),Bs1p,q1(Rn))θ,r=Bsp,r(Rn). (2.4)

If additionally for some and , then

 (Bs0p0,q0(Rn),Bs1p1,q1(Rn))[θ]=Bsp,q(Rn), (2.5)

cf. [17, Theorem 6.4.5] or [57, Section 2.4.1 Theorem]. Using the same notation, we have in particular for the Bessel potential spaces

 (Hs0p(Rn),Hs1p(Rn))θ,r =Bsp,r(Rn), (2.6) (Hs0p0(Rn),Hs1p1(Rn))[θ] =Hsp(Rn)

(cf. [17, Theorem 6.4.5]).

General embedding properties. For any , we have the following embeddings between Besov spaces and Bessel potential spaces:

 Bs+εp,q1(Rn) ↪Hsp(Rn)↪Bs−εp,q2(Rn)for all ε>0,1≤q1,q2≤∞,s∈R, Bsp,min(2,p)(Rn) ↪Hsp(Rn)↪Bsp,max(2,p)(Rn)for all s∈R, (2.7)

cf. e.g. [17, Theorem 6.4.4].

There are the following Sobolev embeddings for Bessel potential spaces:

 Hs1p1(Rn) ↪Hs0p0(Rn)if s1≥s0,s1−np1≥s0−np0, (2.8) Hsp(Rn) ↪Cα(Rn)if%  α=s−np>0, (2.9)

provided that , . In particular, for .

For the Besov spaces a Sobolev-type embedding is given by

 Bs1p1,q(Rn) ↪Bs0p0,q(Rn) if s1≥s0 and s1−np1≥s0−np0, (2.10)

for any . In particular, combining this with (2.1), we get

 Bs1p1,q(Rn)↪Bs0∞,q(Rn)↪Bs0∞,∞(Rn)=Cs0(Rn) (2.11)

whenever . In the opposite direction, we have from (2.1) and (2.2)

 Cα(Rn)=Bα∞,∞(Rn)↪Bα−ε∞,2(Rn) (2.12)

when , . We also note that

 Hs1p1(Rn)∪Bs1p1,p1(Rn)↪Hs0p0(Rn)∩Bs0p0,p0(Rn) (2.13)

if and ; this can be found in [57, Section 2.8.1, equation (17)].

Function spaces over subsets of . The Bessel potential and Besov spaces are defined on a domain with -boundary (see Definition 2.4 below) simply by restriction:

 Hsp(Ω) = {f∈D′(Ω):f=f′|Ω,f′∈Hsp(Rn)}, (2.14) Bsp,q(Ω) = {f∈D′(Ω):f=f′|Ω,f′∈Bsp,q(Rn)},

for and . Here is defined by for all , embedded in by extension by zero. The spaces are equipped with the quotient norms, e.g.,

 ∥f∥Bsp,q(Ω) = inff′∈Bsp,q(Rn):f′|Ω=f∥f′∥Bsp,q(Rn). (2.15)

In particular, is for and equal to the usual Sobolev space of -functions with derivatives up to order in . We recall that there is an extension operator which is a bounded linear operator , for all , , and satisfies for all . This holds when is merely a Lipschitz domain, cf. e.g. Stein [54, Chapter VI, Section 3.2] and trivially carries over to for . Moreover, in view of the fact that is a retract of , one has that all interpolation and Sobolev embedding results for are inherited by the spaces on .

We shall also need the spaces

 Hs0(¯¯¯¯Ω)={u∈Hs(Rn):suppu⊂¯¯¯¯Ω}.

Here, identifies in a natural way with the dual space of , for all , cf. [50, Theorem 3.30]. For integer , equals the closure of in and is usually denoted (see also [50, Theorem 3.33]).

Traces. Next, let us recall the well-known trace theorems: The trace map from to , defined on smooth functions with bounded support, extends by continuity to continuous maps for , , ,

 γ0:Hsp(Rn+) → Bs−1pp,p(Rn−1), γ0:Bsp,q(Rn+) → Bs−1pp,q(Rn−1).

All of these maps are surjective and have continuous right inverses.

Vector-valued Besov and Bessel potential spaces. In the following let be a Banach space. Then , , is defined as the space of strongly measurable functions with

 ∥f∥Lp(Rn;X):=(∫Rn∥f(x)∥pXdx)1p<∞

and is the space of all strongly measurable and essentially bounded functions. Similarly, , , denotes the -valued variant of .

Furthermore, let be the space of smooth rapidly decreasing functions and let denote the space of tempered -valued distributions. Then the -valued variants of the Bessel potential and Besov spaces of order are defined as

 Hsp(Rn;X) :={f∈S′(Rn;X):⟨Dx⟩sf∈Lp(Rn;X)}if 1

where . Here, is defined by

 ⟨p(Dx)f,¯¯¯¯φ⟩=⟨f,¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯p(Dx)φ⟩for all φ∈S(Rn).

We will also make use of the Banach space

 BUC(Rn;X)={f∈L∞(Rn;X):f is uniformly continuous}

with the supremum norm.

The properties of the function spaces discussed above for the scalar case, carry over to the vector-valued case. For details, we refer to e.g. [10, 38] (for the Bochner integral and its properties) and to [8] (for vector-valued function spaces).

In the following we will use some special anisotropic Sobolev spaces.

###### Definition 2.1

Let or , and let , with coordinates . For , , set

 Wk(2,p)(Ω)={f∈L2(I;Lp(Rn−1)):∂αxf∈L2(I;Lp(Rn−1)),|α|≤k}.
###### Lemma 2.2

One has for that

 Wk(2,p)(Ω) ↪ BUC(I;Bk−12p,2(Rn−1)). (2.16)

Here, the trace mapping is surjective from to . Namely, when , then is in with (where is the semigroup generated by ). extends to a function .

###### Proof.

First of all,

 Wk(2,p)(Ω) ↪ L2(I;Hkp(Rn−1))∩Hk(I;Lp(Rn−1)).

To obtain

 L2(I;Hkp(Rn−1))∩Hk(I;Lp(Rn−1))↪BUC(I;Bk−12p,2(Rn−1))

one can apply [17, Corollary 3.12.3] with , , a result from Lions’ trace method of real interpolation, to obtain

 u|xn=0∈(Lp(Rn−1),Hkp(Rn−1))1−12k,2=Bk−12p,2(Rn−1), (2.17)

for every ; the identity follows from (2.6). Next, this is combined with the strong continuity of the translations , , in as in the proof of [8, Chapter III, Theorem 4.10.2].

For the last assertion, let , ; here . Then

 Hkp(Rn−1)={f∈Lp(Rn−1):⟨Dx′⟩kf∈Lp(Rn−1)}=:D(Ak)

for all , , as explained in [17, Theorem 6.2.3]. Now when is given, let for . Since

 Wk(2,p)(Rn+)=⋂0≤j≤kHj(0,∞;Hk−jp(Rn−1)),

we have by [22, Corollary 3.5.6, Theorem 3.4.2] that , so . Here, we use that , as noted above in (2.17). is extended to a function in by a standard “reflection” in , as explained e.g. in [45, Th. I 2.2]. ∎

By use of this lemma we derive the following product estimate, which is essential for the low boundary regularity that we shall allow:

###### Lemma 2.3

Let and be as in Definition 2.1. Then for every and such that there is some such that

 ∥fg∥Hk(Ω)≤Ck,p∥f∥Hk(Ω)∥g∥Wk(2,p)(Ω),

for all , . Moreover, if and , then

 ∥f∂xjg∥L2(Ω)≤Cp∥f∥H1−τ(Ω)∥g∥W1(2,p)(Ω) (2.18)

uniformly with respect to .

###### Proof.

First of all

 Wk(2,p)(Ω)↪BUC(I;Bk−12p,2(Rn−1))↪L∞(Ω)

by (2.16); the second embedding follows from (2.11) since . Furthermore,

Here we apply (2.16) with for the first embedding, and for the second embedding we use (2.8), noting that is equivalent to . Next we observe that for all

 ∂αx(fg)=∑0≤β≤α,β≠α(αβ)∂α−βxf∂βxg+f∂αxg. (2.19)

Since and , where , we have . For the other terms where , we note that

 ∂α−βxf∈L2(I;Hk−|α|+|β|(Rn−1)),∂βxg∈L∞(I;Bk−12−|β|p,2(Rn−1)). (2.20)

One has in general

 ∥uv∥L2(Rn−1)≤CM,M′∥u∥HM(Rn−1)∥v∥BM′p,2(Rn−1), (2.21)

provided that and . This estimate easily follows from the Sobolev-type embedding theorems: If , then from (2.11) we have and the statement is trivial. If , then (2.13) implies with , and implies by (2.8) that with . If , one can choose some such that and apply the preceding case. The estimate is also a consequence of Hanouzet [37, Théorème 3].

Using (2.21) for products of functions as in (2.20), we obtain altogether that for all .

Finally, if , then we have that

 H1−τ(Ω)↪BUC(I;H12−τ(Rn−1))↪BUC(I;Lr(Rn−1)),1r=12−1p.

Therefore

 ∥f∂xjg∥L2(Ω)≤∥f∥L∞(I;Lr)∥∂xjg∥L2(I,Lp)≤C∥f∥H1−τ(Ω)∥g∥W1(