
We prove that for any homogeneous, second order, constant complex coefficient
elliptic system $L$, the Dirichlet problem in $\mathbb{R}^{n}_{+}$ with
boundary data in BMO is wellposed in the class of functions $u$ with
$d\mu_u(x',t):=\nabla u(x',t)^2\,t\,dx'dt$ being a Carleson measure. We
establish a Fatou type theorem guaranteeing the existence of the pointwise
nontangential boundary trace for smooth nullsolutions $u$ of such systems
satisfying the said Carleson measure condition. These imply that BMO can be
characterized as the collection of nontangential pointwise traces of smooth
nullsolutions $u$ to the elliptic system $L$ with the property that $\mu_u$ is
a Carleson measure.
We establish a regularity result for the BMODirichlet problem in the
upperhalf space: the nontangential pointwise trace of any given smooth
nullsolutions of $L$ satisfying the above Carleson measure condition belongs
to Sarason's space VMO if and only if $\mu_u$ satsifies a vanishing Carleson
measure condition. Moreover, we obtain the wellposedness of the Dirichlet
problems when the boundary data are prescribed in MorreyCampanato and
solutions are required to satisfy a vanishing Carleson measure condition of
fractional order.
As a consequence, we characterize the space VMO as the closure in BMO of
classes of smooth functions contained in BMO within which uniform continuity
may be suitably quantified (such as the class of smooth functions satisfying a
H\"older or Lipschitz condition), improving Sarason's classical result
describing VMO as the closure in BMO of the space of uniformly continuous
functions with bounded mean oscillations.
Finally, we show that any Calder\'onZygmund operator $T$ satisfying $T(1)=0$
extends as a bounded linear operator in $\mathrm{VMO}$, and characterize the
membership to $\mathrm{VMO}$ via the action of various classes of singular
integral operators.

We prove wellposedness results for the Dirichlet problem in
$\mathbb{R}^{n}_{+}$ for homogeneous, second order, constant complex
coefficient elliptic systems with boundary data in generalized H\"older spaces
$\mathscr{C}^{\omega}(\mathbb{R}^{n1},\mathbb{C}^M)$ and in generalized
MorreyCampanato spaces $\mathscr{E}^{\omega,p}(\mathbb{R}^{n1},\mathbb{C}^M)$
under certain assumptions on the growth function $\omega$. We also identify a
class of growth functions $\omega$ for which
$\mathscr{C}^{\omega}(\mathbb{R}^{n1},\mathbb{C}^M)=\mathscr{E}^{\omega,p}(\mathbb{R}^{n1},\mathbb{C}^M)$
and for which the aforementioned wellposedness results are equivalent, in the
sense that they have the same unique solution, satisfying natural regularity
properties and estimates.

Given a domain above a Lipschitz graph, we establish solvability results for
strongly elliptic secondorder systems in divergenceform, allowed to have
lowerorder (drift) terms, with $L^p$boundary data for $p$ near $2$ (more
precisely, in an interval of the form
$\big(2\varepsilon,\frac{2(n1)}{n2}+\varepsilon\big)$ for some small
$\varepsilon>0$). The main novel aspect of our result is that the coefficients
of the operator do not have to be constant, or have very high regularity,
instead they will satisfy a natural Carleson condition that has appeared first
in the scalar case. A significant example of a system to which our result may
be applied is the Lam\'e operator for isotropic inhomogeneous materials.
Dealing with genuine systems gives rise to substantial new challenges, absent
in the scalar case. Among other things, there is no maximum principle for
general elliptic systems, and the De Giorgi  Nash  Moser theory may also not
apply. We are, nonetheless, successful in establishing estimates for the
squarefunction and the nontangential maximal operator for the solutions of the
elliptic system described earlier, and use these as alternative tools for
proving $L^p$ solvability results for $p$ near $2$.

We show that the boundedness of the HardyLittlewood maximal operator on a
K\"othe function space ${\mathbb{X}}$ and on its K\"othe dual ${\mathbb{X}}'$
is equivalent to the wellposedness of the $\mathbb{X}$Dirichlet and
$\mathbb{X}'$Dirichlet problems in $\mathbb{R}^{n}_{+}$ in the class of all
secondorder, homogeneous, elliptic systems, with constant complex
coefficients. As a consequence, we obtain that the Dirichlet problem for such
systems is wellposed for boundary data in Lebesgue spaces, variable exponent
Lebesgue spaces, Lorentz spaces, Zygmund spaces, as well as their weighted
versions. We also discuss a version of the aforementioned result which
contains, as a particular case, the Dirichlet problem for elliptic systems with
data in the classical Hardy space $H^1$, and the BeurlingHardy space ${\rm
HA}^p$ for $p\in(1,\infty)$. Based on the wellposedness of the $L^p$Dirichlet
problem we then prove the uniqueness of the Poisson kernel associated with such
systems, as well as the fact that they generate a strongly continuous semigroup
in natural settings. Finally, we establish a general Fatou type theorem
guaranteeing the existence of the pointwise nontangential boundary trace for
nullsolutions of such systems.

We investigate closed, symmetric $L^2(\mathbb{R}^n)$realizations $H$ of
Schr\"odingertype operators $( \Delta
+V)\upharpoonright_{C_0^{\infty}(\mathbb{R}^n \setminus \Sigma)}$ whose
potential coefficient $V$ has a countable number of wellseparated
singularities on compact sets $\Sigma_j$, $j \in J$, of $n$dimensional
Lebesgue measure zero, with $J \subseteq \mathbb{N}$ an index set and $\Sigma =
\bigcup_{j \in J} \Sigma_j$. We show that the defect, $\mathrm{def}(H)$, of $H$
can be computed in terms of the individual defects, $\mathrm{def}(H_j)$, of
closed, symmetric $L^2(\mathbb{R}^n)$realizations of $( \Delta +
V_j)\upharpoonright_{C_0^{\infty}(\mathbb{R}^n \setminus \Sigma_j)}$ with
potential coefficient $V_j$ localized around the singularity $\Sigma_j$, $j \in
J$, where $V = \sum_{j \in J} V_j$. In particular, we prove \[ \mathrm{def}(H)
= \sum_{j \in J} \mathrm{def}(H_j), \] including the possibility that one, and
hence both sides equal $\infty$. We first develop an abstract approach to the
question of decoupling of deficiency indices and then apply it to the concrete
case of Schr\"odingertype operators in $L^2(\mathbb{R}^n)$.
Moreover, we also show how operator (and form) bounds for $V$ relative to
$H_0=  \Delta\upharpoonright_{H^2(\mathbb{R}^n)}$ can be estimated in terms of
the operator (and form) bounds of $V_j$, $j \in J$, relative to $H_0$. Again,
we first prove an abstract result and then show its applicability to
Schr\"odingertype operators in $L^2(\mathbb{R}^n)$.
Extensions to secondorder (locally uniformly) elliptic differential
operators on $\mathbb{R}^n$ with a possibly strongly singular potential
coefficient are treated as well.

The principal aim of this paper is to derive an abstract form of the third
Green identity associated with a proper extension $T$ of a symmetric operator
$S$ in a Hilbert space $\mathfrak H$, employing the technique of quasi boundary
triples for $T$. The general results are illustrated with couplings of
Schr\"{o}dinger operators on Lipschitz domains on smooth, boundaryless
Riemannian manifolds.

For an arbitrary open, nonempty, bounded set $\Omega \subset \mathbb{R}^n$,
$n \in \mathbb{N}$, and sufficiently smooth coefficients $a,b,q$, we consider
the closed, strictly positive, higherorder differential operator $A_{\Omega,
2m} (a,b,q)$ in $L^2(\Omega)$ defined on $W_0^{2m,2}(\Omega)$, associated with
the higherorder differential expression $$ \tau_{2m} (a,b,q) :=
\bigg(\sum_{j,k=1}^{n} (i \partial_j  b_j) a_{j,k} (i \partial_k 
b_k)+q\bigg)^m, \quad m \in \mathbb{N}, $$ and its Kreinvon Neumann extension
$A_{K, \Omega, 2m} (a,b,q)$ in $L^2(\Omega)$. Denoting by $N(\lambda; A_{K,
\Omega, 2m} (a,b,q))$, $\lambda > 0$, the eigenvalue counting function
corresponding to the strictly positive eigenvalues of $A_{K, \Omega, 2m}
(a,b,q)$, we derive the bound $$ N(\lambda; A_{K, \Omega, 2m} (a,b,q)) \leq C
v_n (2\pi)^{n} \bigg(1+\frac{2m}{2m+n}\bigg)^{n/(2m)} \lambda^{n/(2m)} , \quad
\lambda > 0, $$ where $C = C(a,b,q,\Omega)>0$ (with $C(I_n,0,0,\Omega) =
\Omega$) is connected to the eigenfunction expansion of the selfadjoint
operator $\widetilde A_{2m} (a,b,q)$ in $L^2(\mathbb{R}^n)$ defined on
$W^{2m,2}(\mathbb{R}^n)$, corresponding to $\tau_{2m} (a,b,q)$. Here $v_n :=
\pi^{n/2}/\Gamma((n+2)/2)$ denotes the (Euclidean) volume of the unit ball in
$\mathbb{R}^n$.
Our method of proof relies on variational considerations exploiting the
fundamental link between the Kreinvon Neumann extension and an underlying
abstract buckling problem, and on the distorted Fourier transform defined in
terms of the eigenfunction transform of $\widetilde A_{2} (a,b,q)$ in
$L^2(\mathbb{R}^n)$.
We also consider the analogous bound for the eigenvalue counting function for
the Friedrichs extension $A_{F,\Omega, 2m} (a,b,q)$ in $L^2(\Omega)$ of
$A_{\Omega, 2m} (a,b,q)$.
No assumptions on the boundary $\partial \Omega$ of $\Omega$ are made.

In this paper we study the selfadjoint Kreinvon Neumann realization $A_K$
of the perturbed Laplacian $\Delta+V$ in a bounded Lipschitz domain
$\Omega\subset\mathbb{R}^n$. We provide an explicit and selfcontained
description of the domain of $A_K$ in terms of Dirichlet and Neumann boundary
traces, and we establish a Weyl asymptotic formula for the eigenvalues of
$A_K$.

We study the infinitesimal generator of the Poisson semigroup in $L^p$
associated with homogeneous, secondorder, strongly elliptic systems with
constant complex coefficients in the upperhalf space, which is proved to be
the DirichlettoNormal mapping in this setting. Also, its domain is identified
as the linear subspace of the $L^p$based Sobolev space of order one on the
boundary of the upperhalf space consisting of functions for which the
Regularity problem is solvable. Moreover, for a class of systems containing the
Lam\'e system, as well as all secondorder, scalar elliptic operators, with
constant complex coefficients, the action of the infinitesimal generator is
explicitly described in terms of singular integral operators whose kernels
involve firstorder derivatives of the canonical fundamental solution of the
given system. Furthermore, arbitrary powers of the infinitesimal generator of
the said Poisson semigroup are also described in terms of higher order Sobolev
spaces and a higher order Regularity problem for the system in question.
Finally, we indicate how our techniques may adapted to treat the case of higher
order systems in graph Lipschitz domains.

This paper has two main goals. First, we are concerned with the
classification of selfadjoint extensions of the Laplacian
$\Delta\big_{C^\infty_0(\Omega)}$ in $L^2(\Omega; d^n x)$. Here, the domain
$\Omega$ belongs to a subclass of bounded Lipschitz domains (which we term
quasiconvex domains), which contain all convex domains, as well as all domains
of class $C^{1,r}$, for $r\in(1/2,1)$. Second, we establish Kreintype formulas
for the resolvents of the various selfadjoint extensions of the Laplacian in
quasiconvex domains and study the properties of the corresponding
WeylTitchmarsh operators (or energydependent DirichlettoNeumann maps).
One significant technical innovation in this paper is an extension of the
classical boundary trace theory for functions in spaces which lack Sobolev
regularity in a traditional sense, but are suitably adapted to the Laplacian.

For an arbitrary nonempty, open set $\Omega \subset \mathbb{R}^n$, $n \in
\mathbb{N}$, of finite (Euclidean) volume, we consider the minimally defined
higherorder Laplacian $( \Delta)^m\big_{C_0^{\infty}(\Omega)}$, $m \in
\mathbb{N}$, and its Kreinvon Neumann extension $A_{K,\Omega,m}$ in
$L^2(\Omega)$. With $N(\lambda,A_{K,\Omega,m})$, $\lambda > 0$, denoting the
eigenvalue counting function corresponding to the strictly positive eigenvalues
of $A_{K,\Omega,m}$, we derive the bound $$ N(\lambda,A_{K,\Omega,m}) \leq (2
\pi)^{n} v_n \Omega \{1 + [2m/(2m+n)]\}^{n/(2m)} \lambda^{n/(2m)}, \quad
\lambda > 0, $$ where $v_n := \pi^{n/2}/\Gamma((n+2)/2)$ denotes the
(Euclidean) volume of the unit ball in $\mathbb{R}^n$.
The proof relies on variational considerations and exploits the fundamental
link between the Kreinvon Neumann extension and an underlying (abstract)
buckling problem.

We identify a large class of constant (complex) coefficient, second order
elliptic systems for which the Dirichlet problem in the upperhalf space with
data in $L^p$based Sobolev spaces, $1<p<\infty$, of arbitrary smoothness
$\ell$, is wellposed in the class of functions whose nontangential maximal
operator of their derivatives up to, and including, order $\ell$ is
$L^p$integrable. This class includes all scalar, complex coefficient elliptic
operators of second order, as well as the Lam\'e system of elasticity, among
others.

We announce a local $T(b)$ theorem, an inductive scheme, and $L^p$
extrapolation results for $L^2$ square function estimates related to the
analysis of integral operators that act on AhlforsDavid regular sets of
arbitrary codimension in ambient quasimetric spaces. The inductive scheme is a
natural application of the local $T(b)$ theorem and it implies the stability of
$L^2$ square function estimates under the socalled big pieces functor. In
particular, this analysis implies $L^p$ and Hardy space square function
estimates for integral operators on uniformly rectifiable subsets of the
Euclidean space.

In the first (and abstract) part of this survey we prove the unitary
equivalence of the inverse of the Kreinvon Neumann extension (on the
orthogonal complement of its kernel) of a densely defined, closed, strictly
positive operator, $S\geq \varepsilon I_{\mathcal{H}}$ for some $\varepsilon
>0$ in a Hilbert space $\mathcal{H}$ to an abstract buckling problem operator.
This establishes the Krein extension as a natural object in elasticity theory
(in analogy to the Friedrichs extension, which found natural applications in
quantum mechanics, elasticity, etc.).
In the second, and principal part of this survey, we study spectral
properties for $H_{K,\Omega}$, the Kreinvon Neumann extension of the
perturbed Laplacian $\Delta+V$ (in short, the perturbed Krein Laplacian)
defined on $C^\infty_0(\Omega)$, where $V$ is measurable, bounded and
nonnegative, in a bounded open set $\Omega\subset\mathbb{R}^n$ belonging to a
class of nonsmooth domains which contains all convex domains, along with all
domains of class $C^{1,r}$, $r>1/2$.

We consider layer potentials associated to elliptic operators $Lu={\rm
div}(A \nabla u)$ acting in the upper halfspace $\mathbb{R}^{n+1}_+$ for
$n\geq 2$, or more generally, in a Lipschitz graph domain, where the
coefficient matrix $A$ is $L^\infty$ and $t$independent, and solutions of
$Lu=0$ satisfy interior estimates of De Giorgi/Nash/Moser type. A
"Calder\'onZygmund" theory is developed for the boundedness of layer
potentials, whereby sharp $L^p$ and endpoint space bounds are deduced from
$L^2$ bounds. Appropriate versions of the classical "jumprelation" formulae
are also derived. The method of layer potentials is then used to establish
wellposedness of boundary value problems for $L$ with data in $L^p$ and
endpoint spaces.

We initiate the study of the finiteness condition
$\int_{\Omega}u(x)^{\beta}\,dx\leq C(\Omega,\beta)<+\infty$ where
$\Omega\subseteq{\mathbb{R}}^n$ is an open set and $u$ is the solution of the
Saint Venant problem $\Delta u=1$ in $\Omega$, $u=0$ on $\partial\Omega$. The
central issue which we address is that of determining the range of values of
the parameter $\beta>0$ for which the aforementioned condition holds under
various hypotheses on the smoothness of $\Omega$ and demands on the nature of
the constant $C(\Omega,\beta)$. Classes of domains for which our analysis
applies include bounded piecewise $C^1$ domains in ${\mathbb{R}}^n$, $n\geq 2$,
with conical singularities (in particular polygonal domains in the plane),
polyhedra in ${\mathbb{R}}^3$, and bounded domains which are locally of class
$C^2$ and which have (finitely many) outwardly pointing cusps. For example, we
show that if $u_N$ is the solution of the Saint Venant problem in the regular
polygon $\Omega_N$ with $N$ sides circumscribed by the unit disc in the plane,
then for each $\beta\in(0,1)$ the following asymptotic formula holds: %
{eqnarray*} \int_{\Omega_N}u_N(x)^{\beta}\,dx=\frac{4^\beta\pi}{1\beta}
+{\mathcal{O}}(N^{\beta1})\quad{as}\,\,N\to\infty. {eqnarray*} % One of the
original motivations for addressing the aforementioned issues was the study of
sublevel set estimates for functions $v$ satisfying $v(0)=0$, $\nabla v(0)=0$
and $\Delta v\geq c>0$.

The principal aim of this short note is to extend a recent result on Gaussian
heat kernel bounds for selfadjoint $L^2(\Om; d^n x)$realizations, $n\in\bbN$,
$n\geq 2$, of divergence form elliptic partial differential expressions $L$
with (nonlocal) Robintype boundary conditions in bounded Lipschitz domains
$\Om \subset \bbR^n$, where $$ Lu =  \sum_{j,k=1}^n\partial_j
a_{j,k}\partial_k u. $$ The (nonlocal) Robintype boundary conditions are then
of the form $$ \nu\cdot A\nabla u + \Theta\big[u\big_{\partial\Om}\big]=0 \,
\text{on} \, \partial\Omega, $$ where $\Theta$ represents an appropriate
operator acting on Sobolev spaces associated with the boundary $\partial \Om$
of $\Om$, and $\nu$ denotes the outward pointing normal unit vector on
$\partial\Om$.

One of the principal topics of this paper concerns the realization of
selfadjoint operators $L_{\Theta, \Om}$ in $L^2(\Om; d^n x)^m$, $m, n \in
\bbN$, associated with divergence form elliptic partial differential
expressions $L$ with (nonlocal) Robintype boundary conditions in bounded
Lipschitz domains $\Om \subset \bbR^n$. In particular, we develop the theory in
the vectorvalued case and hence focus on matrixvalued differential
expressions $L$ which act as $$ Lu = 
\biggl(\sum_{j,k=1}^n\partial_j\bigg(\sum_{\beta = 1}^m
a^{\alpha,\beta}_{j,k}\partial_k u_\beta\bigg) \bigg)_{1\leq\alpha\leq m},
\quad u=(u_1,...,u_m). $$ The (nonlocal) Robintype boundary conditions are
then of the form $$ \nu \cdot A D u + \Theta \big[u\big_{\partial \Om}\big] =
0 \, \text{on $\partial \Om$}, $$ where $\Theta$ represents an appropriate
operator acting on Sobolev spaces associated with the boundary $\partial \Om$
of $\Om$, $\nu$ denotes the outward pointing normal unit vector on
$\partial\Om$, and $Du:=\bigl(\partial_j
u_\alpha\bigr)_{\substack{1\leq\alpha\leq m 1\leq j\leq n}}$.
Assuming $\Theta \geq 0$ in the scalar case $m=1$, we prove Gaussian heat
kernel bounds for $L_{\Theta, \Om}$ by employing positivity preserving
arguments for the associated semigroups and reducing the problem to the
corresponding Gaussian heat kernel bounds for the case of Neumann boundary
conditions on $\partial \Om$. We also discuss additional zeroorder potential
coefficients $V$ and hence operators corresponding to the form sum $L_{\Theta,
\Om} + V$.

We establish square function estimates for integral operators on uniformly
rectifiable sets by proving a local $T(b)$ theorem and applying it to show that
such estimates are stable under the socalled big pieces functor. More
generally, we consider integral operators associated with AhlforsDavid regular
sets of arbitrary codimension in ambient quasimetric spaces. The local $T(b)$
theorem is then used to establish an inductive scheme in which square function
estimates on socalled big pieces of an AhlforsDavid regular set are proved to
be sufficient for square function estimates to hold on the entire set.
Extrapolation results for $L^p$ and Hardy space versions of these estimates are
also established. Moreover, we prove square function estimates for integral
operators associated with variable coefficient kernels, including the Schwartz
kernels of pseudodifferential operators acting between vector bundles on
subdomains with uniformly rectifiable boundaries on manifolds.

We prove that given any positive integer $k$, for each open set $\Omega$ and
any closed subset $D$ of its closure such that $\Omega$ is locally an
(epsilon,delta)domain near points in the boundary of $\Omega$ not contained in
$D$ there exists a linear and bounded extension operator $E$ mapping, for each
$p\in[1,\infty]$, the space $W^{k,p}_D(\Omega)$ into
$W^{k,p}_D({\mathbb{R}}^n)$. Here, with $O$ denoting either $\Omega$ or the
entire ambient, the space $W^{k,p}_D(O)$ is defined as the completion in the
classical Sobolev space $W^{k,p}(O)$ of compactly supported smooth functions
whose supports are disjoint from $D$. In turn, this result is used to develop a
functional analytic theory for the class $W^{k,p}_D(\Omega)$ (including
intrinsic characterizations, boundary traces and extensions results,
interpolation theorems, among other things) which is then employed in the
treatment of mixed boundary value problems formulated in locally
(epsilon,delta)domains.

We consider DirichlettoNeumann maps associated with (not necessarily
selfadjoint) Schrodinger operators in $L^2(\Omega; d^n x)$, $n=2,3$, where
$\Omega$ is an open set with a compact, nonempty boundary satisfying certain
regularity conditions. As an application we describe a reduction of a certain
ratio of modified Fredholm perturbation determinants associated with operators
in $L^2(\Omega; d^n x)$ to modified Fredholm perturbation determinants
associated with operators in $L^2(\partial\Omega; d^{n1}\sigma)$, $n=2,3$.
This leads to a two and threedimensional extension of a variant of a
celebrated formula due to Jost and Pais, which reduces the Fredholm
perturbation determinant associated with a Schrodinger operator on the
halfline $(0,\infty)$ to a simple Wronski determinant of appropriate
distributional solutions of the underlying Schrodinger equation.

We consider DirichlettoNeumann maps associated with (not necessarily
selfadjoint) Schrodinger operators describing nonlocal interactions in
$L^2(\Omega; d^n x)$, $n\geq 2$, where $\Omega$ is an open set with a compact,
nonempty boundary satisfying certain regularity conditions. As an application
we describe a reduction of a certain ratio of Fredholm perturbation
determinants associated with operators in $L^2(\Omega; d^n x)$ to Fredholm
perturbation determinants associated with operators in $L^2(\partial\Omega;
d^{n1}\sigma)$. This leads to an extension of a variant of a celebrated
formula due to Jost and Pais, which reduces the Fredholm perturbation
determinant associated with a Schr\"odinger operator on the halfline
$(0,\infty)$, in the case of local interactions, to a Wronski determinant of
appropriate distributional solutions of the underlying Schrodinger equation.

The aim of this paper is twofold: First, we characterize an essentially
optimal class of boundary operators $\Theta$ which give rise to selfadjoint
Laplacians $\Delta_{\Theta, \Omega}$ in $L^2(\Omega; d^n x)$ with (nonlocal
and local) Robintype boundary conditions on bounded Lipschitz domains
$\Omega\subset\bbR^n$, $n\in\bbN$, $n\geq 2$. Second, we extend Friedlander's
inequalities between Neumann and Dirichlet Laplacian eigenvalues to those
between nonlocal Robin and Dirichlet Laplacian eigenvalues associated with
bounded Lipschitz domains $\Omega$, following an approach introduced by Filonov
for this type of problems.

We provide a systematic study of boundary data maps, that is, 2 \times 2
matrixvalued DirichlettoNeumann and more generally, RobintoRobin maps,
associated with onedimensional Schrodinger operators on a compact interval
[0,R] with separated boundary conditions at 0 and R. Most of our results are
formulated in the nonselfadjoint context.
Our principal results include explicit representations of these boundary data
maps in terms of the resolvent of the underlying Schrodinger operator and the
associated boundary trace maps, Kreintype resolvent formulas relating
Schrodinger operators corresponding to different (separated) boundary
conditions, and a derivation of the Herglotz property of boundary data maps (up
to right multiplication by an appropriate diagonal matrix) in the special
selfadjoint case.

We prove the unitary equivalence of the inverse of the Kreinvon Neumann
extension (on the orthogonal complement of its kernel) of a densely defined,
closed, strictly positive operator, $S\geq \epsilon I_{\mathcal{H}}$ for some
$\epsilon >0$ in a Hilbert space $\mathcal{H}$ to an abstract buckling problem
operator.
In the concrete case where $S=\bar{\Delta_{C_0^\infty(\Omega)}}$ in
$L^2(\Omega; d^n x)$ for $\Omega\subset\mathbb{R}^n$ an open, bounded (and
sufficiently regular) domain, this recovers, as a particular case of a general
result due to G. Grubb, that the eigenvalue problem for the Krein Laplacian
$S_K$ (i.e., the Kreinvon Neumann extension of $S$), \[ S_K v = \lambda v,
\quad \lambda \neq 0, \] is in onetoone correspondence with the problem of
{\em the buckling of a clamped plate}, \[ (\Delta)^2u=\lambda (\Delta) u
\text{in} \Omega, \quad \lambda \neq 0, \quad u\in H_0^2(\Omega), \] where $u$
and $v$ are related via the pair of formulas \[ u = S_F^{1} (\Delta) v, \quad
v = \lambda^{1}(\Delta) u, \] with $S_F$ the Friedrichs extension of $S$.
This establishes the Krein extension as a natural object in elasticity theory
(in analogy to the Friedrichs extension, which found natural applications in
quantum mechanics, elasticity, etc.).