
Many geometric and analytic properties of sets hinge on the properties of
harmonic measure, notoriously missing for sets of higher codimension. The aim
of this manuscript is to develop a version of elliptic theory, associated to a
linear PDE, which ultimately yields a notion analogous to that of the harmonic
measure, for sets of codimension higher than 1.
To this end, we turn to degenerate elliptic equations. Let $\Gamma \subset
\mathbb R^n$ be an Ahlfors regular set of dimension $d<n1$ (not necessarily
integer) and $\Omega = \mathbb R^n \setminus \Gamma$. Let $L =  {\rm div}
A\nabla$ be a degenerate elliptic operator with measurable coefficients such
that the ellipticity constants of the matrix $A$ are bounded from above and
below by a multiple of ${\rm dist}(\cdot, \Gamma)^{d+1n}$. We define weak
solutions; prove trace and extension theorems in suitable weighted Sobolev
spaces; establish the maximum principle, De GiorgiNashMoser estimates, the
Harnack inequality, the H\"older continuity of solutions (inside and at the
boundary). We define the Green function and provide the basic set of pointwise
and/or $L^p$ estimates for the Green function and for its gradient. With this
at hand, we define harmonic measure associated to $L$, establish its doubling
property, nondegeneracy, changeofthepole formulas, and, finally, the
comparison principle for local solutions.
In another article to appear, we will prove that when $\Gamma$ is the graph
of a Lipschitz function with small Lipschitz constant, we can find an elliptic
operator $L$ for which the harmonic measure given here is absolutely continuous
with respect to the $d$Hausdorff measure on $\Gamma$ and vice versa. It thus
extends Dahlberg's theorem to some sets of codimension higher than 1.

In the recent work [DFM1, DFM2] G. David, J. Feneuil, and the first author
have launched a program devoted to an analogue of harmonic measure for
lowerdimensional sets. A relevant class of partial differential equations,
analogous to the class of elliptic PDEs in the classical context, is given by
linear degenerate equations with the degeneracy suitably depending on the
distance to the boundary.
The present paper continues this line of research and focuses on the criteria
of quantitative absolute continuity of the newly defined harmonic measure with
respect to the Hausdorff measure, $\omega\in A_\infty(\sigma)$, in terms of
solvability of boundary value problems. The authors establish, in particular,
square function estimates and solvability of the Dirichlet problem in BMO for
domains with lowerdimensional boundaries under the underlying assumption
$\omega\in A_\infty(\sigma)$. More generally, it is proved that in all domains
with Ahlfors regular boundaries the BMO solvability of the Dirichlet problem is
necessary and sufficient for the absolute continuity of the harmonic measure.

In the present paper we establish sharp exponential decay estimates for
operator and integral kernels of the (not necessarily selfadjoint) operators
$L=(\nablai\mathbf{a})^TA(\nablai\mathbf{a})+V$. The latter class includes,
in particular, the magnetic Schr\"odinger operator
$\left(\nablai\mathbf{a}\right)^2+V$ and the generalized electric
Schr\"odinger operator ${\rm div }A\nabla+V$. Our exponential decay bounds
rest on a generalization of the FeffermanPhong uncertainty principle to the
present context and are governed by the Agmon distance associated to the
corresponding maximal function. In the presence of a scaleinvariant Harnack
inequality, for instance, for the generalized electric Schr\"odinger operator
with real coefficients, we establish both lower and upper estimates for
fundamental solutions, thus demonstrating sharpness of our results. The only
previously known estimates of this type pertain to the classical Schr\"odinger
operator $\Delta +V$.

We consider the Neumann boundary value problem for an elliptic operator
$L=\mbox{div} (A \, \mbox{grad}) + V$ on a bounded biLipschitz domain in
$\mathbf{R}^n$. More generally, we will treat closed manifolds and manifolds
with boundary. The eigenfunctions of $L$ are often localized, as a result of
disorder of the potential $V$, the matrix of coefficients $A$, irregularities
of the boundary, or all of the above. In earlier work, two of us introduced the
function $u$ solving $Lu=1$, and showed numerically that it strongly reflects
this localization. In this paper, we deepen the connection between the
eigenfunctions and this {\em landscape} function $u$ by proving that its
reciprocal $1/u$ acts as an {\em effective potential}. The effective potential
governs the exponential decay of the eigenfunctions of the system and delivers
information on the distribution of eigenvalues near the bottom of the spectrum.
(This version strengthens and simplifies the results by replacing a global
biLipschitz hypothesis on the domain with a local biLipschitz hypothesis.)

The present paper, along with its sequel, establishes the correspondence
between the properties of the solutions of a class of PDEs and the geometry of
sets in Euclidean space. We settle the question of whether (quantitative)
absolute continuity of the elliptic measure with respect to the surface measure
and uniform rectifiability of the boundary are equivalent, in an optimal class
of divergence form elliptic operators satisfying a suitable Carleson measure
condition. The result can be viewed as a quantitative analogue of the Wiener
criterion adapted to the singular $L^p$ data case.
This paper addresses the free boundary problem under the assumption of
smallness of the Carleson measure of the coefficients. Part II of this work
develops an extrapolation argument to bootstrap this result to the general
case. The ideas in Part I constitute a novel application of techniques
developed in geometric measure theory. They highlight the synergy between
several areas. The ideas developed in this paper are well suited to study
singularities arising in variational problems in a geometric setting.

We present here a model of carrier distribution and transport in
semiconductor alloys accounting for quantum localization effects in disordered
materials. This model is based on the recent development of a mathematical
theory of quantum localization which introduces for each type of carrier a
spatial function called \emph{localization landscape}. These landscapes allow
us to predict the localization regions of electron and hole quantum states,
their corresponding energies, and the local densities of states. We show how
the various outputs of these landscapes can be directly implemented into a
driftdiffusion model of carrier transport and into the calculation of
absorption/emission transitions. This creates a new computational model which
accounts for disorder localization effects while also capturing two major
effects of quantum mechanics, namely the reduction of barrier height (tunneling
effect), and the raising of energy ground states (quantum confinement effect),
without having to solve the Schr\"odinger equation. Finally, this model is
applied to several onedimensional structures such as single quantum wells,
ordered and disordered superlattices, or multiquantum wells, where comparisons
with exact Schr\"odinger calculations demonstrate the excellent accuracy of the
approximation provided by the landscape theory.

This paper introduces a novel method to account for quantum disorder effects
into the classical driftdiffusion model of semiconductor transport through the
localization landscape theory. Quantum confinement and quantum tunneling in the
disordered system change dramatically the energy barriers acting on the
perpendicular transport of heterostructures. In addition they lead to
percolative transport through paths of minimal energy in the 2D landscape of
disordered energies of multiple 2D quantum wells. This model solves the carrier
dynamics with quantum effects selfconsistently and provides a computationally
much faster solver when compared with the Schr\"odinger equation resolution.
The theory also provides a good approximation to the density of states for the
disordered system over the full range of energies required to account for
transport at roomtemperature. The currentvoltage characteristics modeled by
3D simulation of a full nitridebased lightemitting diode (LED) structure
with compositional material fluctuations closely match the experimental
behavior of high quality blue LEDs. The model allows also a fine analysis of
the quantum effects involved in carrier transport through such complex
heterostructures. Finally, details of carrier population and recombination in
the different quantum wells are given.

In 1977 the celebrated theorem of B. Dahlberg established that the harmonic
measure is absolutely continuous with respect to the Hausdorff measure on a
Lipschitz graph of dimension $n1$ in $\mathbb R^n$, and later this result has
been extended to more general nontangentially accessible domains and beyond.
In the present paper we prove the first analogue of Dahlberg's theorem in
higher codimension, on a Lipschitz graph $\Gamma$ of dimension $d$ in $\mathbb
R^n$, $d<n1$, with a small Lipschitz constant. We construct a linear
degenerate elliptic operator $L$ such that the corresponding harmonic measure
$\omega_L$ is absolutely continuous with respect to the Hausdorff measure on
$\Gamma$. More generally, we provide sufficient conditions on the matrix of
coefficients of $L$ which guarantee the mutual absolute continuity of
$\omega_L$ and the Hausdorff measure.

In this paper we establish squarefunction estimates on the double and single
layer potentials with rough inputs for divergence form elliptic operators, of
arbitrary even order 2m, with variable tindependent coefficients in the upper
halfspace.

In this paper we establish well posedness of the Neumann problem with
boundary data in $L^2$ or the Sobolev space $\dot W^2_{1}$, in the half space,
for linear elliptic differential operators with coefficients that are constant
in the vertical direction and in addition are self adjoint. This generalizes
the well known wellposedness result of the second order case and is based on a
higher order and one sided version of the classic Rellich identity, and is the
first known well posedness result for a higher order operator with rough
variable coefficients and boundary data in a Lebesgue or Sobolev space.

We show that if $u$ is a solution to a linear elliptic differential equation
of order $2m\geq 2$ in the halfspace with $t$independent coefficients, and if
$u$ satisfies certain area integral estimates, then the Dirichlet and Neumann
boundary values of $u$ exist and lie in a Lebesgue space $L^p(\mathbb{R}^n)$ or
Sobolev space $\dot W^p_{\pm 1}(\mathbb{R}^n)$. Even in the case where $u$ is a
solution to a second order equation, our results are new for certain values
of~$p$.

In this paper, we establish existence, uniqueness, and scaleinvariant
estimates for fundamental solutions of nonhomogeneous second order elliptic
systems with bounded measurable coefficients in $\mathbb{R}^n$ and for the
corresponding Green functions in arbitrary open sets. We impose certain
nonhomogeneous versions of de GiorgiNashMoser bounds on the weak solutions
and investigate in detail the assumptions on the lower order terms sufficient
to guarantee such conditions.

We introduce a new notion of a harmonic measure for a $d$dimensional set in
$\R^n$ with $d<n1$, that is, when the codimension is strictly bigger than 1.
Our measure is associated to a degenerate elliptic PDE, it gives rise to a
comprehensive elliptic theory, and, most notably, it is absolutely continuous
with respect to the $d$dimensional Hausdorff measure on reasonably nice sets.
This note provides general strokes of the proof of the latter statement for
Lipschitz graphs with small Lipschitz constant.

A recent theoretical breakthrough has brought a new tool, called
\emph{localization landscape}, to predict the localization regions of vibration
modes in complex or disordered systems. Here, we report on the first experiment
which measures the localization landscape and demonstrates its predictive
power. Holographic measurement of the static deformation under uniform load of
a thin plate with complex geometry provides direct access to the landscape
function. When put in vibration, this system shows modes precisely confined
within the subregions delineated by the landscape function. Also the maxima of
this function match the measured eigenfrequencies, while the minima of the
valley network gives the frequencies at which modes become extended. This
approach fully characterizes the low frequency spectrum of a complex structure
from a single static measurement. It paves the way to the control and
engineering of eigenmodes in any vibratory system, especially where a
structural or microscopic description is not accessible.

The amplitude of localized quantum states in random or disordered media may
exhibit long range exponential decay. We present here a theory that unveils the
existence of an effective potential which finely governs the confinement of
these states. In this picture, the boundaries of the localization subregions
for low energy eigenfunctions correspond to the barriers of this effective
potential, and the long range exponential decay characteristic of Anderson
localization is explained as the consequence of multiple tunneling in the dense
network of barriers created by this effective potential. Finally, we show that
the Weyl's formula based on this potential turns out to be a remarkable
approximation of the density of states for a large variety of onedimensional
systems, periodic or random.

In the present paper we prove that for any open connected set
$\Omega\subset\mathbb{R}^{n+1}$, $n\geq 1$, and any $E\subset \partial \Omega$
with $\mathcal{H}^n(E)<\infty$, absolute continuity of the harmonic measure
$\omega$ with respect to the Hausdorff measure on $E$ implies that $\omega_E$
is rectifiable. This solves an open problem on harmonic measure which turns out
to be an old conjecture even in the planar case $n=1$.

In the present paper we sketch the proof of the fact that for any open
connected set $\Omega\subset\mathbb{R}^{n+1}$, $n\geq 1$, and any $E\subset
\partial \Omega$ with $0<\mathcal{H}^n(E)<\infty$, absolute continuity of the
harmonic measure $\omega$ with respect to the Hausdorff measure on $E$ implies
that $\omega_E$ is rectifiable.

In this paper we establish squarefunction estimates on the double and single
layer potentials for divergenceform elliptic operators, of arbitrary even
order 2m, with variable tindependent coefficients in the upper halfspace.
This generalizes known results for variablecoefficient secondorder operators,
and also for constantcoefficient higherorder operators.

Recent years have brought significant advances in the theory of higher order
elliptic equations in nonsmooth domains. Sharp pointwise estimates on
derivatives of polyharmonic functions in arbitrary domains were established,
followed by the higher order Wiener test. Certain boundary value problems for
higher order operators with variable nonsmooth coefficients were addressed,
both in divergence form and in composition form, the latter being adapted to
the context of Lipschitz domains. These developments brought new estimates on
the fundamental solutions and the Green function, allowing for the lack of
smoothness of the boundary or of the coefficients of the equation. Building on
our earlier account of history of the subject, this survey presents the current
state of the art, emphasizing the most recent results and emerging open
problems.

The behavior of solutions to the biharmonic equation is wellunderstood in
smooth domains. In the past two decades substantial progress has also been made
for the polyhedral domains and domains with Lipschitz boundaries. However, very
little is known about higher order elliptic equations in the general setting.
In this paper we introduce new integral identities that allow to investigate
the solutions to the biharmonic equation in an arbitrary domain. We establish:
(1) boundedness of the gradient of a solution in any threedimensional domain;
(2) pointwise estimates on the derivatives of the biharmonic Green function;
(3) Wienertype necessary and sufficient conditions for continuity of the
gradient of a solution.

In the present paper we prove that for any open connected set
$\Omega\subset{\mathbb R}^{n+1}$, $n\geq 1$, and any $E\subset \partial\Omega$
with $0<{\mathcal H}^n(E)<\infty$ absolute continuity of the harmonic measure
$\omega$ with respect to the Hausdorff measure on $E$ implies that $\omega_E$
is rectifiable.

Let $L$ be a homogeneous divergence form higher order elliptic operator with
complex bounded measurable coefficients and $(p_(L),\, p_+(L))$ be the maximal
interval of exponents $q\in[1,\,\infty]$ such that the semigroup
$\{e^{tL}\}_{t>0}$ is bounded on $L^q(\mathbb{R}^n)$. In this article, the
authors establish the nontangential maximal function characterizations of the
associated Hardy spaces $H_L^p(\mathbb{R}^n)$ for all $p\in(0,\,p_+(L))$,
which, when $p=1$, answers a question asked by Deng et al. in [J. Funct. Anal.
263 (2012), 604674]. Moreover, the authors characterize $H_L^p(\mathbb{R}^n)$
via various versions of square functions and Lusinarea functions associated to
the operator $L$.

In the present paper we establish the Wiener test for boundary regularity of
the solutions to the polyharmonic operator. We introduce a new notion of
polyharmonic capacity and demonstrate necessary and sufficient conditions on
the capacity of the domain responsible for the regularity of a polyharmonic
function near a boundary point.
In the case of the Laplacian the test for regularity of a boundary point is
the celebrated Wiener criterion of 1924. It was extended to the biharmonic case
in dimension three by [Mayboroda, Maz'ya, Invent. Math. 2009]. As a preliminary
stage of this work, in [Mayboroda, Maz'ya, Invent. Math. 2013] we demonstrated
boundedness of the appropriate derivatives of solutions to the polyharmonic
problem in arbitrary domains, accompanied by sharp estimates on the Green
function. The present work pioneers a new version of capacity and establishes
the Wiener test in the full generality of the polyharmonic equation of
arbitrary order.

The localization subregions of stationary waves in continuous disordered
media have been recently demonstrated to be governed by a hidden landscape that
is the solution of a Dirichlet problem expressed with the wave operator. In
this theory, the strength of Anderson localization confinement is determined by
this landscape, and continuously decreases as the energy increases. However,
this picture has to be changed in discrete lattices in which the eigenmodes
close to the edge of the first Brillouin zone are as localized as the low
energy ones. Here we show that in a 1D discrete lattice, the localization of
low and high energy modes is governed by two different landscapes, the high
energy landscape being the solution of a dual Dirichlet problem deduced from
the low energy one using the symmetries of the Hamiltonian. We illustrate this
feature using the onedimensional tightbinding Hamiltonian with random onsite
potentials as a prototype model. Moreover we show that, besides unveiling the
subregions of Anderson localization, these dual landscapes also provide an
accurate overal estimate of the localization length over the energy spectrum,
especially in the weak disorder regime.

Let $E\subset \mathbb{R}^{n+1}$, $n\ge 2$, be a uniformly rectifiable set of
dimension $n$. Then bounded harmonic functions in $\Omega:=
\mathbb{R}^{n+1}\setminus E$ satisfy Carleson measure estimates, and are
"$\varepsilon$approximable". Our results may be viewed as generalized versions
of the classical F. and M. Riesz theorem, since the estimates that we prove are
equivalent, in more topologically friendly settings, to quantitative mutual
absolute continuity of harmonic measure, and surface measure.