  Given a principal bundle $P\to M$ over a Riemannian manifold with compact
structure group $G$, let us consider a stationary Yang-Mills connection $A$
with energy $\int_M |F_A|^2\le \Lambda$. If we consider a sequence of such
connections $A_i$, then it is understood that up to subsequence we can converge
$A_i\to A$ to a singular limit connection such that the energy measures
converge $|F_{A_i}|^2 dv_g\to |F_A|^2dv_g +\nu$, where $\nu=e(x)d\lambda^{n-4}$
is the $n-4$ rectifiable defect measure. Our main result is to show, without
additional assumptions, that for $n-4$ a.e. point the energy density $e(x)$ may
be computed explicitly as the sum of the bubble energies arising from blow ups
at $x$. Each of these bubbles may be realized as a Yang Mills connection over
$S^4$ itself.
  Our second main theorem is to show this hessian bound holds automatically.
Precisely, given a connection $A$ as above we have the apriori estimate $\int_M
|\nabla^2 F_A| < C(\Lambda,\dim G,M)$ for the curvature. It is important to
note this result is proved in tandem with the energy quantization, and not
before it. Indeed, we will in fact prove an effective version of the energy
identity, and it is this effective version which will lead to both the $L^1$
hessian bound and the classical energy quantization results. In the course of
the proof we will provide a quantitative version of the bubble tree
decomposition which hold in all dimensions with effective estimates for a fixed
stationary connections. To produce to strongest estimates in the paper we
introduce an $\epsilon$-gauge condition, which generalizes the usual Coulomb
gauge and which will exist, with effective control, even over singular regions.
On these $\epsilon$-gauges we will provide a new superconvexity estimate which
will be a key tool in analyzing higher dimensional annular regions.
