The L^p Neumann problem for parabolic operators with coefficients satisfying small Carleson condition

In this paper, we resolve the question of whether the Neumann problem for the parabolic PDE $-\partial_tu + \mathrm{div}(A\nabla u)=0$ on a Lipschitz cylinder $\mathcal O\times\mathbb R$ is solvable for some $p\in (1,\infty)$ under the assumption that the matrix $A$ is elliptic with bounded and measurable coefficients that satisfy a natural Carleson condition (a parabolic analog of the so-called DKP-condition). We prove that for any $1<p<\infty$ the Neumann problem is solvable under the assumption that both the Carleson norm of coefficients and the Lipschitz constant of the domain are sufficiently small (with dependence on $p$). The question of what happens in the "large Carleson norm/large Lipschitz constant" regime remains open, and even for elliptic PDEs this question has only been resolved in two dimensions. This paper complements results from our recent manuscript (by the same authors) in which the parabolic regularity problem has been fully resolved in both the small and large Carleson norm regime. Previously, the Dirichlet problem had been resolved under the same conditions by various authors.