Let $\mathbb{D}^2$ be the closed two-dimensional unit disk, and let $g:\mathbb{D}^2 \to \mathbb{R}$ be a non-constant harmonic function (smooth up to the boundary).
Does there exist a sequence of smooth one-forms $\sigma_n$ on $\mathbb{D}^2$ such that
- $\sigma_n \to dg$ in $L^2$.
- $\sigma_n$ do not vanish on $\mathbb{D}^2$.
- $\|\delta d \sigma_n \|_{L^1} \to 0, \| d \delta \sigma_n \|_{L^1} \to 0$
That is, I want to approximate the one-form $dg$, which may have zeroes in the domain, with nowhere vanishing forms, such that certain second derivatives become negligible in the limit.
I know how to do that without achieving the third condition.
Let me explain a bit further on that condition:
Write $\sigma =dg \in \Omega^1(\mathbb{D}^2)$. Since $g$ is harmonic, $\delta \sigma=0$. Since $\sigma$ is exact, we also have $d\sigma=0$. In fact, for a one-form $\sigma \in \Omega^1(\mathbb{D}^2)$, $d\sigma=\delta \sigma=0$ is equivalent to the existence of a harmonic function $g$ on $\mathbb{D}^2$ such that $\sigma=dg$.
Thus, our desired limit form $\sigma =dg$ satisfies $\|\delta d \sigma \|_{L^1} = \| d \delta \sigma \|_{L^1} =0$.
Hence, if we could approximate $\sigma$ by nowhere-vanishing $\sigma_n$ in $W^{2,1}$, we were done. However, we cannot in general approximate even in $W^{1,1}$ while staying nowhere-vanishing. (See also this comment).
Two more comments:
Trying to approximate (in $L^2$) $g$ by harmonic functions $g_n$ in with non-vanishing differential also doesn't work: This creates a too fast approximation, due to rigidity properties of harmonic functions, again hitting the topological obstruction.
Even trying to approximate (again in $L^2$) $\sigma=dg$ with $\sigma_n$ satisfying $d\delta \sigma_n=0, \delta d \sigma_n=0$ cannot work: If $\sigma_n=f^1_ndx+f^2_ndy$, then these conditions imply that $f^1_n$ and $f^2_n$ are harmonic, which again implies a too fast convergence.