[guided]We begin with an arbitrary harmonic representative because the theorem asserts that every such representative must vanish. The word harmonic here means harmonic for the Dolbeault Laplacian
\begin{align*}
\Delta_{\bar\partial,E}=\bar\partial_E\bar\partial_E^*+\bar\partial_E^*\bar\partial_E.
\end{align*}
Equivalently, on a compact Hermitian manifold, harmonicity gives the two first-order equations
\begin{align*}
\bar\partial_E\alpha=0
\end{align*}
and
\begin{align*}
\bar\partial_E^*\alpha=0.
\end{align*}
The identity we need is the Bochner-Kodaira-Nakano identity with the sign convention in which the curvature term is exactly
\begin{align*}
\bigl([i\Theta(E,h),\Lambda]\alpha,\alpha\bigr)_{L^2}.
\end{align*}
Here $\nabla'$ is the $(1,0)$-part of the Chern connection on $(E,h)$, and $(\nabla')^*$ is its $L^2$-adjoint. The hypotheses of the Bochner-Kodaira-Nakano identity are met: $(X,\omega)$ is compact Kähler, $(E,h)$ is a Hermitian holomorphic vector bundle, and $\alpha$ is a smooth $E$-valued $(n,q)$-form. Applied to $\alpha$, the identity gives
\begin{align*}
\|\bar\partial_E\alpha\|_{L^2}^2+\|\bar\partial_E^*\alpha\|_{L^2}^2
=
\|\nabla'\alpha\|_{L^2}^2+\|(\nabla')^*\alpha\|_{L^2}^2+\bigl([i\Theta(E,h),\Lambda]\alpha,\alpha\bigr)_{L^2}.
\end{align*}
The reason this identity is useful is that harmonicity annihilates the left-hand side. Indeed, both summands on the left are squared $L^2$ norms of zero forms, so
\begin{align*}
\|\bar\partial_E\alpha\|_{L^2}^2+\|\bar\partial_E^*\alpha\|_{L^2}^2=0.
\end{align*}
Therefore the identity becomes
\begin{align*}
0
=
\|\nabla'\alpha\|_{L^2}^2+\|(\nabla')^*\alpha\|_{L^2}^2+\bigl([i\Theta(E,h),\Lambda]\alpha,\alpha\bigr)_{L^2}.
\end{align*}
This is the central point: the analytic equation for harmonic forms has been converted into a sum of two nonnegative energy terms plus the curvature term appearing in the hypothesis.[/guided]