$\mathrm{L}^p$ bounds for parabolic Riesz transforms with rough coefficients: The case $1<p \leq 2$

We establish the first results on $\mathrm{L}^p$ bounds for Riesz transforms associated with non-autonomous second order parabolic differential operators in divergence form with bounded coefficients that depend measurably on all variables. In the case of complex coefficients, we identify the maximal open range of exponents $1<p \leq2$ through the availability of $\mathrm{L}^p$ resolvent bounds. This open range always contains the lower parabolic Sobolev conjugate of $2$ and the result is sharp in spatial dimension $n \geq 2$. For real coefficients, we prove extrapolation to the full range. Our argument relies on novel space-time off-diagonal bounds based on two complementary geometries: parabolic cubes on small scales and regions modeled after the half-order time derivative of a parabolic Bessel potential on large scales.