1. Introduction
Let $\Omega$
be an open subset of $\mathbf {R}^N$
. We consider a general constant-coefficient system of linear partial differential equations acting on functions $u : \Omega \to E$
, that is, a partial differential operator of the form
where $E,\,F$
are finite-dimensional $\mathbf {R}$
-spaces. Here $\alpha = (\alpha _1,\,\ldots,\,\alpha _N) \in (\mathbb {N}_0)^N$
is a multi-index with modulus $|\alpha | = \alpha _1 + \cdots + \alpha _N$
and $\partial ^\alpha = \partial _1^{\alpha _1} \circ \cdots \circ \partial _N^{\alpha _N}$
is the composition of distributional directional partial derivatives. The Fourier transform establishes a one-to-one correspondence between homogeneous operators and their associated principal symbol $\mathbb {A}: \mathbf {R}^N \to \operatorname {Hom}({E,\,F})$
, which in this context is given by the $k$
-homogeneous tensor-valued polynomial
Suppose now that $u \in \mathscr D'(\Omega ; E)$
is a distribution such that $\mathcal {A} u$
is a zero-order distribution, i.e. represented by an $F$
-valued Radon measure. Let $g: \Omega \to F$
be the polar vector of $\mathcal {A} u$
. If, additionally, $\mathcal {A}$
is a constant-rank operator (see e.g. [Reference Arroyo-Rabasa and Simental8, Reference Murat18, Reference Raiţă19]), the De Philippis–Rindler theorem [Reference De Philippis and Rindler11, Theorem 1.1] establishes that $g$
is constrained to take values on the image cone at singular points. More precisely, if $U \subset \Omega$
is a ${\mathscr {L}}_N$
-negligible Borel set, where ${\mathscr {L}}_N$
stands for the $N$
-dimensional Lebesgue measure, then
This property extends a classical result of Alberti [Reference Alberti3], which says that if the distributional gradient $Du$
of a map $u: \mathbf {R}^N \to \mathbf {R}^M$
is represented by an $\mathbf {R}^{M \times N}$
-valued Radon measure, then its polar vector must take values in the cone of rank-one matrices on ${\mathscr {L}}_N$
-negligible sets: if ${\mathscr {L}}_N(U) = 0$
, then
In contrast with this restriction over the polar of $Du$
at singular points, Alberti showed [Reference Alberti1, Theorem 3] that the absolutely continuous part of a gradient measure is fully unconstrained: if $f \in L^1(\mathbf {R}^N;\mathbf {R}^{M \times N})$
, then there exists $u \in BV(\mathbf {R}^N;\mathbf {R}^M)$
such thatFootnote 1
and
Alberti (see [Reference Alberti1, Reference Alberti2]) also established other Lusin-type properties for gradients of arbitrary order. He showed that if $f : \Omega \to \mathbf {R}^{N^k}_\mathrm {sym}$
is continuous, then for any positive measure $\lambda$
on $\Omega$
and any smallness constant $\epsilon > 0$
, one may find a compact set $K \subset \Omega$
and a function $u \in C^k(\Omega )$
such that
and satisfying $L^\infty$
-estimates for $D^k u$
in terms of $f$
. Building on these ideas, Francos [Reference Francos16] (see also [Reference Delladio12]) showed that for any given Borel $f : \Omega \to \mathbf {R}^{N^k}_\mathrm {sym}$
and $\sigma > 0$
, there exists a function $g \in C^{k-1}(\Omega )$
, $k$
-times differentiable almost everywhere, satisfying
Driven by applications to higher-order variational problems where derivatives give rise to surface energies, Fonseca et al. [Reference Fonseca, Leoni and Paroni15, Theorem 1.4] established an analogue of (1.2) for the Hessian operator. More precisely, they showed that if $f : \Omega \to \mathbf {R}^{N \times N}_\mathrm {sym}$
, then there exists $u \in W^{1,1}(\Omega )$
, with $\nabla u \in BV(\Omega ;\mathbf {R}^{N})$
, satisfying
Unlike the Lusin-type properties for higher-order gradients by Francos, this generalization to second-order derivatives is not followed by a straightforward iteration of Alberti's property (1.2). The reason for this drawback is the presence of the symmetry constraints of higher-order curl-free fields.
In light of this digression, the recent developments in the study of fine properties in $BV^\mathcal {A}$
spaces (see e.g. [Reference Ambrosio, Coscia and Dal Maso4, Reference Arroyo-Rabasa6, Reference Arroyo-Rabasa and Skorobogatova9, Reference Breit, Diening and Gmeineder10, Reference Diening and Gmeineder13]), and the De Philippis–Rindler theorem, we are led to ask the following natural question: is the absolutely continuous part of an $\mathcal {A}$
-gradient measure fully unconstrained? Our main results (Theorems 2.1 and 2.2) establish that, at least when $\mathcal {A}$
is an operator of order $k \le 2$
, this is indeed the case; regardless of $\mathcal {A}$
satisfying the constant-rank condition. Our proof hinges on Proposition 2.4, where we show the previous question is equivalent to the validity of (1.2) for Hessians of arbitrary order. This is the main reason why our result is restricted to the case $k = 1,\,2$
.
2. Notation and results
In all that follows, $\Omega$
will denote an open subset of $\mathbf {R}^N$
. We will write ${\mathscr {L}}_N$
to denote the $N$
dimensional Lebesgue measure and ${\mathscr {H}}_{N-1}$
to denote the ($N-1$
)-dimensional Hausdorff outer measure on $\mathbf {R}^N$
. The space $BV(\Omega ;E)$
, of $E$
-valued vector fields with bounded variation on $\Omega$
, consists of all integrable maps $u: \Omega \to E$
whose distributional gradient can be represented by a finite Radon measure taking values on $E \otimes \mathbf {R}^N$
. For such functions, we shall write $\nabla u \in L^1(\Omega ;E \otimes \mathbf {R}^N)$
to denote the absolutely continuous part of $Du$
with respect to ${\mathscr {L}}_N$
. It is well known (see e.g. [Reference Ambrosio, Fusco and Pallara5]) that such maps are Lebesgue continuous outside of a countably ${\mathscr {H}}_{N-1}$
rectifiable set $J_u \subset \Omega$
, with orientation normal $\nu _u$
, called the jump set of $u$
. Moreover, the map $u$
has Lebesgue one-sided limits $u^+(x),\,u^-(x)$
for ${\mathscr {H}}_{N-1}$
almost every $x \in J_u$
with respect to the normal direction $\nu _u(x)$
, and their difference
defines a ${\mathscr {H}}_{N-1}$
-integrable map on $J_u$
. By the Radon–Nikodym theorem, the gradient of a map $u \in BV(\Omega ;E)$
can be decomposed as $Du = \nabla u \, {\mathscr {L}}_N + D^s u$
where $D^s u \perp {\mathscr {L}}_N$
. The singular part $D^s u$
can be further decomposed into mutually singular measures as
where (the Cantor part) $D^c u$
is a measure vanishing on countable unions of ${\mathscr {H}}_{N-1}$
-finite Borel sets, and the so-called jump part
is the restriction of $Du$
to the set of jump discontinuities of $u$
. The subspace $SBV(\Omega ;E)$
of $E$
-valued functions of special bounded variation consists of all functions $u \in BV(\Omega ;E)$
with rectifiable singular gradient, i.e.
We also define
the essential range of the operator $\mathcal {A}$
. It is easy to see that every distribution $\mathcal {A} u$
takes values in this space, which is the smallest space with this property.
With these considerations in mind, we can state our main results. We begin by stating the result for first-order operators:
Theorem 2.1 Let $\mathcal {A}$
be a first-order operator as in (1.1) and let $f : \Omega \to \mathbf {F}_{\mathcal {A}}$
be integrable. Then there exists a map $u \in {SBV}(\Omega ;E)$
satisfying
and
for some constant $C$
that only depends on $\mathcal {A}$
.
A similar statement holds for second-order operators:
Theorem 2.2 Let $\mathcal {A}$
be a second-order operator as in (1.1) and let $f : \Omega \to \mathbf {F}_\mathcal {A}$
be integrable. Then there exists a map $u \in W^{1,1}(\Omega ;E)$
, with $\nabla u \in SBV(\Omega ;E \otimes \mathbf {R}^N)$
, satisfying
and
where $C$
is a constant that only depends on $\mathcal {A}$
.
Remark 2.3 The essential range $\mathbf {F}_\mathcal {A}$
defined in (2.1) coincides with the space $\operatorname {span} \left \{\, \mathbb {A}^k(\xi )[e] \ {:}\ \xi \in \mathbf {R}^N,\, e \in E \,\right \} = \operatorname {span} \, \{\mathcal {I}_{\mathcal {A}}\}$
of all $\mathcal {A}$
-gradient amplitudes in Fourier space (for a proof see [Reference Arroyo-Rabasa, De Philippis and Rindler7, Section 2.5]). For example, if $\mathcal {A} = D^2$
is the Hessian, then $\mathbf {F}_{D^2} = \operatorname {span} \left \{\, \xi \otimes \xi \ {:}\ \xi \in \mathbf {R}^N \,\right \} =: \mathbf {R}^{N \times N}_\mathrm {sym}$
.
Theorems 2.1 and 2.2 will follow directly from the equivalence result below (Proposition 2.4), Alberti's original result for gradients (in the first-order case) and its analogue for Hessians by Fonseca, Leoni and Paroni (in the second-order case). Before stating Proposition 2.4, it will be convenient to introduce some basic notation for symmetric tensors and gradients of higher order.
Symmetric tensors
Let $r \ge 0$
be an integer. For a vector $v \in \mathbf {R}^N$
, we write $v^{\otimes ^r}$
to denote the tensor that results by taking the tensorial product of $v$
with itself $r$
-times (with the convention $v^{\otimes ^0} = 1$
). We consider the subspace
consisting of all symmetric $r$
th order tensors on $\mathbf {R}^N$
. In the following, we shall consider a contraction 
$\langle \,\bullet\,,\,\,\bullet\, \rangle _r : (E \otimes \mathbf {R}^{N^r}_\mathrm {sym} \times \mathbf {R}^N) \to E$ associated with the canonical inner product 
$(\,\bullet\, ,\, \,\bullet\,)$ on $\otimes ^r \mathbf {R}^N$
as
Notice that under our convention for $r=0$
, we have 
$\langle e,\,\,\bullet\, \rangle _0 = e$.
Jump densities for higher-order $BV^k$
-spaces
Let $k \ge 2$
be an integer. We define the space $BV^k(\Omega )$
as the space of integrable functions $u: \Omega \to \mathbf {R}$
whose distributional $k$
th order gradient $D^ku$
can be represented by a Radon measure taking values on $\mathbf {R}^{N^k}_\mathrm {sym}$
; the spaces $BV^k_\mathrm {loc}(\Omega )$
, $BV^k(\Omega ;E)$
are defined accordingly with this definition in the obvious manner. By classical elliptic regularity theory, there is a natural continuous embedding $BV^k_\mathrm {loc}(\Omega ;E) \hookrightarrow W^{k-1,1}_\mathrm {loc}(\Omega ;E)$
. Therefore, the distributional gradient $D^{k-1} u$
of a map $u$
in $BV^k_\mathrm {loc}(\Omega ;E)$
can be represented by an integrable map $\nabla ^{k-1}u : \Omega \to E \otimes \mathbf {R}^{N^{k-1}}_\mathrm {sym}$
, that is,
Moreover, in this case, the tensor field $w:= \nabla ^{k-1} u$
has bounded variation, the jump set $J_w$
of $w$
is a countably ${\mathscr {H}}_{N-1}$
rectifiable set where $w$
has approximate one-sided limits $[w]\in E \otimes \mathbf {R}^{N^{k-1}}_\mathrm {sym}$
, with respect to a fixed orientation $\nu _w$
of $J_w$
, and the difference map $[w]$
is integrable on $J_w$
. Since $D^k u$
is a symmetric-valued tensor measure, it follows that the ${\mathscr {H}}_{N-1}$
density of the jump part of $Dw = D^k u$
is of the form
for some ${\mathscr {H}}_{N-1}$
-integrable map $\lambda _w : J_w \to E$
. Using that $|\nu _w|^2 = 1$
, we can characterize the $E$
-coordinate coefficient $\lambda _w$
directly in terms of $[w]$
and the contraction pairing defined above. Indeed, $\langle [w],\,\nu _w \rangle _{k-1}\cdot |\nu _w|^2 = \langle [w]\otimes \nu _w,\,\nu _w \rangle _{k} = \lambda _w |\nu _\omega |^{2k} = \lambda _w$
. In particular, it follows that
For consistency, we set $\nabla ^0 u := u$
so that, when $k=1$
, we simply get
With these considerations in mind, we are now in the position to state the following equivalence between Lebesgue–Lusin properties for $k$
th order gradients and arbitrary homogeneous operators of order $k$
:
Proposition 2.4 The following statements are equivalent:
(a) Let $\mathcal {A}$
be a $k$th order linear operator as in (1.1) and let $f : \Omega \to \mathbf {F}_{\mathcal {A}}$ be an integrable map. Then there exists $u \in BV^k(\Omega ;E)$, satisfying
\[ \mathcal{A} u \; = \; f\,{\mathscr{L}}_N \; + \; \mathbb{A}^k(\nu_w) \langle [w],\nu_w \rangle_{k-1} \, {\mathscr{H}}_{N-1} \mathbin J_w\,, \]where $w := \nabla ^{k-1} u$
. Moreover, $w \in SBV(\Omega ;E \otimes \mathbf {R}^{k-1}_\mathrm {sym})$ and
\[ \|u\|_{W^{k-1,1}(\Omega)} + \int_\Omega |\nabla w| \, {\rm d}x + \int_{J_w} |[w]| \, {\rm d}{\mathscr{H}}_{N-1} \le C \int_\Omega |f| \, {\rm d}x \]for some constant $C$
depending on $\mathcal {A}$.(b) Let $f : \Omega \to \mathbf {R}^{N^k}_\mathrm {sym}$
be an integrable map. Then there exists $u \in BV^k(\Omega )$ satisfying
\[ D^k u \; = \; f\,{\mathscr{L}}_N \; + \; [w] \otimes \nu_w \; {\mathscr{H}}_{N-1} \, {\mathbin J_w,} \]where $w := \nabla ^{k-1} u$
. Moreover, $w \in SBV(\Omega ;\mathbf {R}^{N^{k-1}}_\mathrm {sym})$ and
\[ \|u\|_{W^{k-1,1}(\Omega)} + \int_\Omega |\nabla w| \, {\rm d}x + \int_{J_w} |[w]| \, {\rm d}{\mathscr{H}}_{N-1} \le C \int_\Omega |f| \, {\rm d}x \]for some constant $C$
depending on $N$ and $k$.
Proof of Proposition 2.4 The implication (a) $\Longrightarrow$
(b) is straightforward from the following observations. Firstly, the principal symbol of the $k$
th order Hessian is given by $D^k(\xi ) = \xi ^{\otimes ^k}$
. Therefore, $D^k$
has the form (1.1) with $E = \mathbf {R}$
and $F = \mathbf {R}^{N^k}_\mathrm {sym}$
, by setting $A_\alpha = M_\alpha$
, where the family $\left \{\, M_\alpha \ {:}\ |\alpha | =k \,\right \}$
is the canonical orthonormal basis of $\mathbf {R}^{N^k}_\mathrm {sym}$
. Hence, from (2.4) and (a) we conclude that $u \in BV^k(\Omega )$
satisfies
with $w = \nabla ^{k-1} u \in SBV(\Omega ;\mathbf {R}^{N^{k-1}}_\mathrm {sym})$
, satisfying
This proves the first implication since $D^k$
implicitly fixes the spatial dimension $N$
.
We now show (b) $\Longrightarrow$
(a). The first step is to use the alternative jet expression $A[D^k u] = \mathcal {A} u$
, where $A : E \otimes \mathbf {R}^{N^k}_\mathrm {sym} \to \mathbf {F}_\mathcal {A}$
is the (unique) linear map satisfying
The existence of $A$
is a direct consequence of the universal property of the tensor product and the $k$
-linearity of the principal symbol on the frequency variable $\xi$
. For the sake of simplicity, let us write $\mathbf {E}_k =E \otimes \mathbf {R}^{N^k}_\mathrm {sym}$
. By construction, we have
Let us consider the Moore–Penrose quasi-inverse $A^{\dagger}$
associated with $A$
. This is an element of $\operatorname {Hom}(\mathbf {F}_\mathcal {A};\mathbf {E}_k)$
satisfying the fundamental identity
In particular, by our assumption on $f$
, $A^{\dagger} f(x)$
is well-defined almost everywhere on $\Omega$
and $A^{\dagger} f \in L^1(\Omega ;\mathbf {E}_k)$
. Now, we make use of the assumption (b) over each $E$
-coordinate to find $u \in BV^k(\Omega ;E)$
satisfying
where $w = \nabla ^{k-1} u\, \in SBV(\Omega ;\mathbf {R}^{N^k}_\mathrm {sym})$
satisfies the estimate
for some constant $C$
depending on $N$
, $k$
and $\dim \mathbf {F}_\mathcal {A}$
. Here, as usual $D^j w = \tilde g\, {\mathscr {H}}_{N-1} \mathbin J_w$
is the jump part of $Dw$
. Notice that in this case $\tilde g = \mathbf 1_{J_w} [w] \otimes \nu _w$
. Pre-composing this expression with $A$
, we get
Since $A$
is a linear map, we may pull in and distribute $A$
into the densities that belong to the sum of the right-hand side, hence concluding that
where $g =\mathbf 1_{J_w} \cdot A([w] \otimes \nu _w)$
. In passing to the last equality in the formula above, we have used the almost everywhere pointwise restriction $f(x) \in \mathbf {F}_\mathcal {A}$
. Now, we use the fact that at jump points (cf. (2.4)), it holds
From this we obtain that $g =\mathbf 1_{J_w}\cdot \mathbb {A}^k(\nu _w)\langle [w],\, \nu _w \rangle _{k-1}$
as desired. Lastly, the bound on the Sobolev norm and total variation of $D^ku$
follow from (2.7), the estimates
and the fact that $\|A^{\dagger} \|_{\mathbf {F}_\mathcal {A} \to \mathbf {E}_k} = \|A^{-1}\|_{\operatorname {Im} A \to \mathbf {F}_\mathcal {A}}$
depends solely on $\mathcal {A}$
.Footnote 2 This finishes the proof.
Proof of Theorems 2.1 and 2.2 The proof of Theorem 2.1 follows directly from the validity of the Lusin property for gradients ([Reference Alberti1, Theorem 3]), the previous proposition (with $k =1$
), and the equivalence in Remark 2.3. The proof of Theorem 2.2, on the other hand, follows from the previous proposition (with $k = 2$
) and [Reference Fonseca, Leoni and Paroni15, Theorem 1.4].
An application for normal currents
Let $m$
be a non-negative integer. An $m$
-dimensional current $T$
on $\Omega$
is an element of the continuous dual $\mathscr D_m(\Omega ) = \mathscr D^m(\Omega )^*$
, where
is the space of smooth and compactly supported $m$
-forms on $\Omega$
. If $m \ge 1$
, the (distributional) boundary operator on currents is defined by duality through the rule
Here, $d$
is the exterior derivative operator acting on $\mathscr D^m(\Omega )$
. The boundary operator then maps $\mathscr D_m(\Omega )$
into $\mathscr D_{m-1}(\Omega )$
. The mass of a current $T \in \mathscr D_m(\Omega )$
is defined as
A current $T$
is called normal if $T$
is representable by a finite Radon measure and either $\partial T$
is representable by a finite Radon measure or $m = 0$
. In particular, every normal $k$
-current is a distribution represented by a measure taking values on the space $\wedge _m \mathbf {R}^N = (\wedge ^m \mathbf {R}^N)^*$
of $m$
-vectors. The space of normal $m$
-dimensional currents on $\Omega$
is denoted by $\mathbf {N}_m(\Omega )$
.
Given an element $\xi \in \mathbf {R}^N \cong \wedge _1 \mathbf {R}^N$
, we write $\xi ^*$
to denote the $1$
-covector
The interior multiplication operator $\llcorner : \wedge _{p} \mathbf {R}^N \times \wedge ^q \mathbf {R}^N \to \wedge _{p-q} \mathbf {R}^N$
is defined as the adjoint of the exterior multiplication (see e.g. [Reference Federer14, Chapter 1.5]):
With this notation in mind, we get the following direct application of Theorem 2.1 for the boundary operator acting on normal currents:
Corollary 2.5 Let $m \in [0,\,N)$
be an integer and let $S : \Omega \to \wedge _m \mathbf {R}^N$
be an integrable $m$
-vector field. Then, there exists a normal current $T \in \mathbf {N}_{m+1}(\Omega ) \cap SBV(\Omega ;\wedge _{m+1}\mathbf {R}^N)$
satisfying
Moreover,
Remark 2.6 A similar result holds for integrable $m$
-forms with $m \in (0,\,N]$
. More precisely, if $\omega \in L^1(\Omega ;\wedge ^m \mathbf {R}^N)$
, then there exists an $(m-1)$
-form $\phi \in SBV(\Omega ;\wedge ^{m-1} \mathbf {R}^N)$
satisfying (see also [Reference Moonens and Pfeffer17, Proposition 2.1])
and
Proof. Since the symbol of the exterior derivative is precisely the exterior multiplication (i.e. $d_m(\xi )\alpha = \xi ^*\wedge \alpha$
), the symbol of the boundary operator $\partial$
on $({m+1})$
-vectors is precisely (see e.g. [Reference Federer14, §4.1.7.])
In particular, the boundary operator $\partial$
defines a constant-coefficient first-order operator from $\mathscr D_{m+1}(\Omega )$
to $\mathscr D_m(\Omega )$
. Since (cf. [Reference Federer14, §1.5.2])
it follows that $\wedge _m \mathbf {R}^N = \operatorname {span}\{ \operatorname {Im} \partial _{m+1}(\xi ) :|\xi | = 1\} = \mathbf {F}_\partial$
. We may thus apply Theorem 1 to find $T \in SBV(\Omega ;\wedge _m \mathbf {R}^N)$
satisfying the desired properties: the fact that $T \in \mathbf {N}_{m+1}(\Omega )$
follows directly from the fact that $\mathbf {M}(T) = \|T\|_{L^1(\Omega )}$
and that $\mathbf {M}(\partial T) \le |DT|(\Omega )$
.
Corollary 2.7 Let $m \in [{1},\,N)$
be an integer and let $S : \Omega \to \wedge _m \mathbf {R}^N$
be an integrable $m$
-vector field. There exists a countably ${\mathscr {H}}_{N-1}$
rectifiable set $J$
with (oriented) normal $\nu$
, and a Borel $(m+1)$
-vector field $g : J \to \wedge _{m+1} \mathbf {R}^N$
satisfying
and
Proof. It is sufficient to use the expression for $\partial T$
in the corollary above and observe that $\partial (\partial T) = 0$
.
The canonical isomorphism $\iota : \mathbf {R}^N \to \wedge _1 \mathbf {R}^N$
, identifying vector fields in $\mathbf {R}^N$
with $1$
-vectors, induces an isometry between vector fields with bounded measure divergence and normal $1$
-currents. In particular, an exciting and direct application of the previous corollary is the following rectifiable completion for systems of vector fields to systems of solenoidal measures.
Corollary 2.8 If $N \ge 2$
, then every vector field ${\vec f} \in L^1(\Omega ;\mathbf {R}^{N})$
is the absolutely continuous part of a solenoidal field measure up to a rectifiable measure. More precisely, there exists a countably ${\mathscr {H}}_{N-1}$
rectifiable set $\Gamma \subset \Omega$
and a Borel vector field $\vec a : \Gamma \to \mathbf {R}^{N}$
satisfying the differential constraint
and also the tangential constraint
Moreover,
for a constant $C$
that depends only on $N$
.
Proof. Applying the assertion of Corollary 2.7 (here we are using that $N \ge 2$
) to the integrable $1$
-vector field $\iota \vec f$
yields the existence of a countably ${\mathscr {H}}_{N-1}$
rectifiable set $\Gamma \subset \Omega$
and an integrable $2$
-vector $g : \Gamma \to \wedge _2 \mathbf {R}^N$
satisfying
and
By the identification discussed before, we conclude that
where $\vec a := \iota ^{-1}(- g \llcorner \nu ^*) :\Gamma \to \mathbf {R}^N$
is integrable on $\Gamma$
and satisfies the asserted $L^1$
bounds on $\Gamma$
. To see the tangential properties of the $\vec a$
, it suffices to observe that $\vec a(x) \cdot \nu (x) = \langle g(x) \llcorner \nu (x)^*,\,\nu (x)^* \rangle = \langle g(x),\, \nu (x)^* \wedge \nu (x)^* \rangle = 0$
for ${\mathscr {H}}_{N-1}$
almost every $x \in \Gamma$
(on Lebesgue points of $g$
and $\nu$
).
Acknowledgements
I would like to thank Guido De Philippis and Giacomo del Nin for very timely insight and comments about the problem discussed in this note. I would also like to thank the referee for their careful corrections and valuable suggestions. Firstly, for bringing to my attention the work of Fonseca, Leoni and Paroni on the Lusin property for functions with bounded Hessians (crucial in establishing Theorem 2). Secondly, for suggesting the identification of divergence-bounded fields and normal $1$
-currents to give a simpler proof of Corollary 6.
