Harmonic vector fields are critical points of Dirichlet energy function on the set of all unit vector fields on $M$, which is defined as follows:
$$E(X):=\frac{1}{2}\int_M\|dX\|^2\mathrm{dVol_g}\qquad X: (M,g)\to (TM,g_s), \|X\|=1.$$
Where it's Euler-Lagrange equation satisfies the following: $$\bar{\Delta}X=\lambda X,\quad \lambda\in C^\infty(M)$$ where $\bar{\Delta}$ denote rough Laplacian.
$p$-Harmonic vector field: If we use $p$-energy function instead of Dirichlet energy function: $$E_p(X):=\frac{1}{p}\int_M\|\mathrm{d}X\|^p\mathrm{dVol_g},\qquad X: (M,g)\to (TM,g_s),$$ The Euler-Lagrange equation of $E_p(.)$ in all vector fields (not in unit case) is as follows (for function case see this paper): $$\|dX\|^{p-2}\left(\tau(X) + dX\left(\mathrm{grad}(\log \|dX\|^{p-2})\right)\right) =0.$$ where $\tau(X)=\mathrm{div}(dX)$ called tension field of $X$.
My questions are :
Question 1: What is the Euler-Lagrange of $E_p(.)$ on the set of all unit vector fields on $M$?
Question 2: Can one hope to get statements which is generalized and new (not same to exactly) 2-harmonic (harmonic vector field) case?
Thanks.