Let $\M,\N$ be $d$-dimensional oriented Riemannian manifolds. Consider for a moment the Dirichlet functionalLet $E:C^{\infty}(M,N) \to \R$, give$\tilde E:C^{\infty}(M,N) \to \R$ be defined by $$E(f)=\int_{\M} | df|^2 \text{Vol}_{\M}. $$ It is a classic fact that a map $f:\M \to \N$ is a critical point of $E$ if and only if $\delta(df)=0$, where $\delta$ is the adjoint of $\nabla^{f^*{\TN}}$. (such a map $f$ is called harmonic).
I am asking whether or not an analogous characterization of critical points holds for this modified functional:
$$\tilde E(f)=\frac{1}{2}\int_{\M} | \bigwedge^2 df|^2 \text{Vol}_{\M}$$$$\tilde E(f)=\frac{1}{2}\int_{\M} | \bigwedge^2 df|^2 \text{Vol}_{\M}.$$ $\tilde E(f)$ measures the mean action of $f$ on 2$D$-cubes. (How $f$ affects areas of surfaces, locally).
I want to know if the set ofQuestion: Does every critical pointspoint of $\tilde E$ coincides with set of maps satisfyingsatisfy $\delta_{\nabla^{\Lambda_2(f^*{\TN})}} \big( \bigwedge^2 df \big) =0$.?
OneOf course, by counting degrees of freedom, this doesn't look reasonable.
The converse direction follows easily from the Euler-Lagrange equation of $\tilde E$: I proved here that the E-L equation is $$h_{f^*\TN} \big( \tr_{\TM}\big( df \otimes \delta_{\nabla^{\Lambda_2(f^*{\TN})}}(\bigwedge^2 df)\big) \bigg)=0,$$
Question: Does every critical point of $\tilde E$ satisfy $\delta_{\nabla^{\Lambda_2(f^*{\TN})}} \big( \bigwedge^2 df \big) =0$?
Edit 1:
Of course, by counting degreesA natural way to produce critical points of freedom, this doesn'ta functional is to look reasonableat its symmetries. $\tilde E$ is conformally-invariant exactly at dimension $4$. So, in dimension $4$ conformal maps are critical. However, I proved that they also satisfy the stronger equation $\delta_{\nabla^{\Lambda_2(f^*{\TN})}} \big( \bigwedge^2 df \big) =0$.
(In fact I proved that a conformal map $\M^d \to \N^d$ satisfies $\delta_{\nabla^{\Lambda_k(f^*{\TN})}} \big( \bigwedge^k df \big) =0$ if and only if $d=2k$ or it is a homothety).
Edit 12:
As explainedI shall prove below (see "Edit 3"), in dimension $d=2$ the answer is positive. (Note that in this case the "number of constraints" is the same). Focusing on the simplest next case, we consider $d=3,\M=\N=\R^3$ with the flat metrics.
$$ \delta \big( df \wedge df \big)(X)=\sum_i \nabla_{e_i}(df \wedge df)(e_i,X)=\sum_i (\nabla_{e_i} df)(e_i) \wedge df(X) + df(e_i) \wedge (\nabla_{e_i} df)(X)$$ $$ =\Delta f \wedge df(X)+\sum_i df(e_i) \wedge (\nabla_{e_i} df)(X).$$ $\nabla_{e_i}$ is the standard component-wise derivative w.r.t the variable $x_i$, or the Levi-Civita connection of $\R^3$, if you prefer.
We can try to look first for harmonic counter-examples, but so far I failed doing even that.
Note that in the second equation we refer to the trace taken w.r.t the Euclidean metric: $$ \tr \big( df \otimes \delta(df \wedge df) \big)=\sum_i df(e_i) \otimes \bigg( \delta \big( df \wedge df \big)(e_i) \bigg)=0.$$
Edit 23: