Green's function for elliptic systems: existence and Delmotte-Deuschel bounds

We prove that for an open domain $D \subset \mathbb{R}^d$ with $d \geq 2$ , for every (measurable) uniformly elliptic tensor field $a$ and for almost every point $y \in D$ , there exists a unique Green's function centred in $y$ associated to the vectorial operator $-\nabla \cdot a\nabla$ in D. In particular, when $d > 2$ this result also implies the existence of the fundamental solution for elliptic systems, i.e. the Green function for $-\nabla \cdot a\nabla$ in $\mathbb{R}^d$. Moreover, introducing an ensemble $\langle\cdot \rangle$ over the set of uniformly elliptic tensor fields, under the assumption of stationarity we infer for the fundamental solution $G$ some pointwise bounds for $\langle |G(\cdot; x,y)|\rangle$, $\langle|\nabla_x G(\cdot; x,y)|\rangle$ and $\langle |\nabla_x\nabla_y G(\cdot; x,y)|\rangle$. These estimates scale optimally in space and provide a generalization to systems of the bounds obtained by Delmotte and Deuschel for the scalar case.