Canonical variational completion of 4D Gauss-Bonnet gravity

A proposal to obtain a finite contribution of second derivative order to the gravitational field equations in (D = 4) dimensions from a renormalized Gauss-Bonnet term in the action has recently received a wave of attention, and triggered a discussion whether the employed renormalization procedure yields a well-defined theory. One of the main criticisms is based on the fact that the resulting field equations cannot be obtained as the Euler-Lagrange equations from a diffeomorphism-invariant action. In this work we use techniques from the inverse calculus of variation as an independent confirmation that the suggested truncated Gauss-Bonnet field equations cannot be variational, in any dimension. For this purpose, we employ canonical variational completion, based on the notion of Vainberg-Tonti Lagrangian, which consists in adding a canonically defined correction term to a given system of equations, so as to make them derivable from an action. We find that in $D>4$ the suggested field equations can be variationally completed, which yields a theory with fourth order field equations. In $D=4$ the variationally completed theory diverges. Our findings are in line with Lovelock's theorem, which states that, in 4 dimensions, the unique second-order Euler-Lagrange equations arising from a scalar density depending on the metric tensors and its derivatives, are the Einstein equations with a cosmological constant.