In https://arxiv.org/abs/1903.03601, on page 13, the propagator of 4d Chern-Simons theory is computed, in the gauge $D^iA_i=0$, where $D^i = (\partial_x,\partial_y,4\partial_z)$.
The gauge-fixed action is $$\begin{equation} S + S_{\text{gf}} = \frac{1}{2\pi}\int_{M}d^4w\, \Bigg(\varepsilon^{ijk}\textrm{Tr}\bigg(A_i\partial_jA_k + \frac{2}{3}A_iA_jA_k\bigg) -\xi^{-1}\textrm{Tr}\Big((D^iA_i)^2\Big)\Bigg)\,, \end{equation}\tag{3.1}$$ where $\xi$ is a gauge-fixing parameter. The propagator $\Delta_{ij}(w,w')$ on $M$ satisfies $$\begin{equation} \frac{1}{\pi}(\varepsilon^{ijk}\partial_j + \xi^{-1}D^iD^k)\Delta_{k\ell}(w,w') = c\,\delta^i_{\,\ell}\,\delta^{(4)}(w-w')\,, \end{equation}\tag{3.2}$$ where $c=t_a\otimes t_b\,(\kappa^{-1})^{ab} \, = t_a\otimes t^a \in \mathfrak{g}^{\otimes 2}$, where the Killing form is $\kappa_{ab}=\textrm{Tr}(t_at_b)$. In Landau gauge ($\xi=0$), the propagator is claimed to be $$\begin{equation} \label{PropagatoronM} \Delta_{ij}(w,w') = -\frac{c}{4\pi}\varepsilon_{ijk}D^k\bigg(\frac{1}{\|w-w'\|^2}\bigg) = -\frac{c}{4\pi}\varepsilon_{ijk}D^k\bigg(\frac{1}{(x-x')^2+(y-y')^2+(z-z')(\bar{z}-\bar{z}')}\bigg) \,. \end{equation}\tag{3.3}$$
Away from $w=w'$, this is obvious from the antisymmetry of the Levi-Civita tensor. However, it is not clear to me why this is the case for $w=w'$. In fact, by computing $D^k\Delta_{kl}$ explicitly, I obtain zero. What am I missing?