On optimal H\"older regularity of solutions to the equation $\mathrm{\Delta }u+b\cdot \mathrm{\nabla }u=0$ in two dimensions

We show that for an $L^2$ drift $b$ in two dimensions, if the Hardy norm of $\text{div }b$ is small, then the weak solutions to $\mathrm{\Delta }u+b\cdot \mathrm{\nabla }u=0$ have the same optimal H"older regularity as in the case of divergence-free drift, that is, $u\in C_{\text{loc}}^{\alpha }$ for all $\alpha \in (0,1)$.