
We prove weighted $L_{p,q}$estimates for divergence type higher order
elliptic and parabolic systems with irregular coefficients on Reifenberg flat
domains. In particular, in the parabolic case the coefficients do not have any
regularity assumptions in the time variable. As functions of the spatial
variables, the leading coefficients are permitted to have small mean
oscillations. The weights are in the class of Muckenhoupt weights $A_p$. We
also prove the solvability of the systems in weighted Sobolev spaces.

We study the stationary Stokes system with Dini mean oscillation coefficients
in a domain having $C^{1,\rm{Dini}}$ boundary. We prove that if $(u, p)$ is a
weak solution of the system with zero Dirichlet boundary condition, then
$(Du,p)$ is continuous up to the boundary. We also prove a weak type$(1,1)$
estimate for $(Du, p)$.

We study Green functions for stationary Stokes systems satisfying the
conormal derivative boundary condition. We establish existence, uniqueness, and
various estimates for the Green function under the assumption that weak
solutions of the Stokes system are continuous in the interior of the domain.
Also, we establish the global pointwise bound for the Green function under the
additional assumption that weak solutions of the conormal derivative problem
for the Stokes system are locally bounded up to the boundary. We provide some
examples satisfying such continuity and boundedness properties.

We study the stationary Stokes system in divergence form. The coefficients
are assumed to be merely measurable in one direction and have Dini mean
oscillations in the other directions. We prove that if $(u,p)$ is a weak
solution of the system, then $(Du,p)$ is bounded and its certain linear
combinations are continuous. We also prove a weak type$(1,1)$ estimate for
$(Du,p)$ under a stronger assumption on the $L^1$mean oscillation of the
coefficients. The corresponding results up to the boundary on a half ball are
also established. These results are new even for elliptic equations and
systems.

We prove the solvability in Sobolev spaces of the conormal derivative problem
for the stationary Stokes system with irregular coefficients on bounded
Reifenberg flat domains. The coefficients are assumed to be merely measurable
in one direction, which may differ depending on the local coordinate systems,
and have small mean oscillations in the other directions. In the course of the
proof, we use a local version of the Poincar\'e inequality on Reifenberg flat
domains, the proof of which is of independent interest.

We study the Green function for the stationary Stokes system with bounded
measurable coefficients in a bounded Lipschitz domain $\Omega\subset
\mathbb{R}^n$, $n\ge 3$. We construct the Green function in $\Omega$ under the
condition $(\bf{A1})$ that weak solutions of the system enjoy interior H\"older
continuity. We also prove that $(\bf{A1})$ holds, for example, when the
coefficients are $\mathrm{VMO}$. Moreover, we obtain the global pointwise
estimate for the Green function under the additional assumption $(\bf{A2})$
that weak solutions of Dirichlet problems are locally bounded up to the
boundary of the domain. By proving a priori $L^q$estimates for Stokes systems
with $\mathrm{BMO}$ coefficients on a Reifenberg domain, we verify that
$(\bf{A2})$ is satisfied when the coefficients are $\mathrm{VMO}$ and $\Omega$
is a bounded $C^1$ domain.

We establish the existence and the pointwise bound of the fundamental
solution for the stationary Stokes system with measurable coefficients in the
whole space $\mathbb{R}^d$, $d \ge 3$, under the assumption that weak solutions
of the system are locally H\"older continuous. We also discuss the existence
and the pointwise bound of the Green function for the Stokes system with
measurable coefficients on $\Omega$, where $\Omega$ is an unbounded domain such
that the divergence equation is solvable. Such a domain includes, for example,
half space and an exterior domain.

The aim of this paper is to investigate Green's function for parabolic and
elliptic systems satisfying a possibly nonlocal Robintype boundary condition.
We construct Green's function for parabolic systems with timedependent
coefficients satisfying a possibly nonlocal Robintype boundary condition
assuming that weak solutions of the system are locally H\"{o}lder continuous in
the interior of the domain, and as a corollary we construct Green's function
for elliptic system with a Robintype condition. Also, we obtain Gaussian bound
for Robin Green's function under an additional assumption that weak solutions
of Robin problem are locally bounded up to the boundary. We provide some
examples satisfying such a local boundedness property, and thus have Gaussian
bounds for their Green's functions.

We study the Neumann Green's function for second order parabolic systems in
divergence form with timedependent measurable coefficients in a cylindrical
domain $\mathcal{Q}=\Omega\times (\infty,\infty)$, where $\Omega\subset
\mathbb{R}^n$ is an open connected set such that a multiplicative Sobolev
embedding inequality holds there. Such a domain includes, for example, a
bounded Sobolev extension domain, a special Lipschitz domain, and an unbounded
domain with compact Lipschitz boundary. We construct the Neumann Green's
function in $\mathcal{Q}$ under the assumption that weak solutions of the
systems satisfy an interior H\"older continuity estimate. We also establish
global Gaussian bounds for Neumann Green's function under an additional
assumption that weak solutions with zero Neumann data satisfy a local
boundedness estimate. In the scalar case, such a local boundedness estimate is
a consequence of De GiorgiMoserNash theory holds for equations with bounded
measurable coefficients in Sobolev extension domains, while in the vectorial
case, one may need to impose further regularity assumptions on the coefficients
of the system as well as on the domain to obtain such an estimate. We present a
unified approach valid for both the scalar and vectorial cases and discuss some
applications of our results including the construction of Neumann functions for
second order elliptic systems with measurable coefficients in two dimensional
domains.

We study Neumann functions for divergence form, second order elliptic systems
with bounded measurable coefficients in a bounded Lipschitz domain or a
Lipschitz graph domain. We establish existence, uniqueness, and various
estimates for the Neumann functions under the assumption that weak solutions of
the system enjoy interior H\"older continuity. Also, we establish global
pointwise bounds for the Neumann functions under the assumption that weak
solutions of the system satisfy a certain natural local boundedness estimate.
Moreover, we prove that such a local boundedness estimate for weak solutions of
the system is in fact equivalent to the global pointwise bound for the Neumann
function. We present a unified approach valid for both the scalar and the
vectorial cases.