Maximal regularity for elliptic operators with second-order discontinuous coefficients

We prove maximal regularity for parabolic problems associated to the second-order elliptic operator $$\begin{aligned} L =\Delta +(a-1)\sum _{i,j=1}^N\frac{x_ix_j}{|x|^2}D_{ij}+c\frac{x}{|x|^2}\cdot \nabla -b|x|^{-2} \end{aligned}$$ with $$a>0$$ and $$b,\ c$$ real coefficients.