Critical metrics for the Dirac-Yamabe eigenvalue ratio and Einstein-Dirac equations

Let $(M^n,g)$ be a closed spin manifold of dimension $n\geq 3$ with positive Yamabe constant. We prove that critical metrics for the mixed eigenvalue ratio $\mathcal H(g)=\frac{λ_1(D_g)^2}{λ_1(L_g)},$ where $D_g$ is the Dirac operator and $L_g$ is the conformal Laplacian are characterized by a harmo