• 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 well-posed 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 null-solutions $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 null-solutions $u$ to the elliptic system $L$ with the property that $\mu_u$ is a Carleson measure. We establish a regularity result for the BMO-Dirichlet problem in the upper-half space: the nontangential pointwise trace of any given smooth null-solutions 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 well-posedness of the Dirichlet problems when the boundary data are prescribed in Morrey-Campanato 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\'on-Zygmund 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 well-posedness 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}^{n-1},\mathbb{C}^M)$ and in generalized Morrey-Campanato spaces $\mathscr{E}^{\omega,p}(\mathbb{R}^{n-1},\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}^{n-1},\mathbb{C}^M)=\mathscr{E}^{\omega,p}(\mathbb{R}^{n-1},\mathbb{C}^M)$ and for which the aforementioned well-posedness 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 second-order systems in divergence-form, allowed to have lower-order (drift) terms, with $L^p$-boundary data for $p$ near $2$ (more precisely, in an interval of the form $\big(2-\varepsilon,\frac{2(n-1)}{n-2}+\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 square-function 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 Hardy-Littlewood maximal operator on a K\"othe function space ${\mathbb{X}}$ and on its K\"othe dual ${\mathbb{X}}'$ is equivalent to the well-posedness of the $\mathbb{X}$-Dirichlet and $\mathbb{X}'$-Dirichlet problems in $\mathbb{R}^{n}_{+}$ in the class of all second-order, homogeneous, elliptic systems, with constant complex coefficients. As a consequence, we obtain that the Dirichlet problem for such systems is well-posed 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 Beurling-Hardy space ${\rm HA}^p$ for $p\in(1,\infty)$. Based on the well-posedness 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 null-solutions of such systems.
  • We investigate closed, symmetric $L^2(\mathbb{R}^n)$-realizations $H$ of Schr\"odinger-type operators $(- \Delta +V)\upharpoonright_{C_0^{\infty}(\mathbb{R}^n \setminus \Sigma)}$ whose potential coefficient $V$ has a countable number of well-separated 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\"odinger-type 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\"odinger-type operators in $L^2(\mathbb{R}^n)$. Extensions to second-order (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, higher-order differential operator $A_{\Omega, 2m} (a,b,q)$ in $L^2(\Omega)$ defined on $W_0^{2m,2}(\Omega)$, associated with the higher-order 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 Krein--von 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 self-adjoint 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 Krein--von 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 self-adjoint Krein-von 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 self-contained 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, second-order, strongly elliptic systems with constant complex coefficients in the upper-half space, which is proved to be the Dirichlet-to-Normal 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 upper-half 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 second-order, 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 first-order 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 self-adjoint 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 quasi-convex 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 Krein-type formulas for the resolvents of the various self-adjoint extensions of the Laplacian in quasi-convex domains and study the properties of the corresponding Weyl--Titchmarsh operators (or energy-dependent Dirichlet-to-Neumann 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 higher-order Laplacian $(- \Delta)^m\big|_{C_0^{\infty}(\Omega)}$, $m \in \mathbb{N}$, and its Krein--von 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 Krein--von 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 upper-half space with data in $L^p$-based Sobolev spaces, $1<p<\infty$, of arbitrary smoothness $\ell$, is well-posed 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 Ahlfors-David regular sets of arbitrary codimension in ambient quasi-metric 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 so-called 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 Krein--von 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 Krein--von 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 half-space $\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\'on-Zygmund" 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 "jump-relation" formulae are also derived. The method of layer potentials is then used to establish well-posedness 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^{\beta-1})\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 self-adjoint $L^2(\Om; d^n x)$-realizations, $n\in\bbN$, $n\geq 2$, of divergence form elliptic partial differential expressions $L$ with (nonlocal) Robin-type 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) Robin-type 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 self-adjoint 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) Robin-type boundary conditions in bounded Lipschitz domains $\Om \subset \bbR^n$. In particular, we develop the theory in the vector-valued case and hence focus on matrix-valued 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) Robin-type 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 zero-order 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 so-called big pieces functor. More generally, we consider integral operators associated with Ahlfors-David regular sets of arbitrary codimension in ambient quasi-metric spaces. The local $T(b)$ theorem is then used to establish an inductive scheme in which square function estimates on so-called big pieces of an Ahlfors-David 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 Dirichlet-to-Neumann maps associated with (not necessarily self-adjoint) 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^{n-1}\sigma)$, $n=2,3$. This leads to a two- and three-dimensional 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 half-line $(0,\infty)$ to a simple Wronski determinant of appropriate distributional solutions of the underlying Schrodinger equation.
  • We consider Dirichlet-to-Neumann maps associated with (not necessarily self-adjoint) 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^{n-1}\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 half-line $(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 self-adjoint Laplacians $-\Delta_{\Theta, \Omega}$ in $L^2(\Omega; d^n x)$ with (nonlocal and local) Robin-type 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 matrix-valued Dirichlet-to-Neumann and more generally, Robin-to-Robin maps, associated with one-dimensional Schrodinger operators on a compact interval [0,R] with separated boundary conditions at 0 and R. Most of our results are formulated in the non-self-adjoint 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, Krein-type 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 self-adjoint case.
  • We prove the unitary equivalence of the inverse of the Krein--von 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 Krein--von Neumann extension of $S$), \[ S_K v = \lambda v, \quad \lambda \neq 0, \] is in one-to-one 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.).