A massless Fierz-Pauli action is the quadratic flat-background action for a symmetric metric perturbation that yields the vacuum Linearized Einstein equations. One normalization, up to a total derivative, is
Varying the symmetric field with compact support gives the displayed linearized Einstein tensor coefficient. Its linearized coordinate-gauge symmetry is consistent with the contracted Bianchi identity. It is the massless gravitational action, distinct from imposing the massive Fierz-Pauli equations with five polarizations.
Let the symmetric metric perturbation be varied with compact support. Integrating the mixed derivative product by parts makes its integral equal to that of . Thus the supplied massless Fierz-Pauli action has, up to a boundary term, density
For example, the equivalence of the mixed term follows from commuting flat-space partial derivatives after moving one derivative off . Here and .
Integrating each variation by parts, the four terms contribute respectively
The symmetrization in the second term is required because the varied field is symmetric. Comparing this coefficient with the Einstein tensor perturbation above gives
Thus arbitrary compactly supported variations give , exactly the vacuum Linearized Einstein equations. The minus sign and overall normalization of the action do not change those equations; boundary conditions justify the discarded total derivatives.