Gravity with higher-curvature terms and second-order field equations: f(R) meets Gauss-Bonnet

General relativity is expected to break down in the high-curvature regime. Beyond an effective field theory treatment with higher-order operators, it is important to identify consistent theories with highercurvature terms at the nonperturbative level. Two well-studied examples are f(R) gravity and Einsteindilaton- Gauss-Bonnet (EdGB) gravity. The former shares the same vacuum solutions as general relativity, including black holes, while the latter suffers from well-posedness issues due to quadratic curvature terms in the strong-coupling regime. We show that combining these two theories leads to genuinely new phenomena beyond their simple superposition. The resulting framework falls outside Horndeski’s class, as it can be recast as a gravitational theory involving two nonminimally coupled scalar fields with nontrivial mutual interactions. This construction naturally extends EdGB gravity to include arbitrary higher-curvature terms, providing a versatile setting to address fundamental questions. Focusing on quadratic and quartic corrections, we find that: (i) black holes are modified by f(R) terms, unlike the case without Gauss-Bonnet interactions; (ii) the resulting solutions retain the qualitative nonperturbative features of EdGB black holes with certain couplings, such as a minimum mass and multiple branches; (iii) a nontrivial mechanism suppresses the divergence of the Ricci scalar in the black-hole interior; (iv) still, even with quartic corrections, the singularity structure and elliptic regions inside the horizon remain similar to those of pure EdGB gravity. This suggests that, at the nonperturbative level, the theory’s ill-posedness cannot be resolved by adding individual higher-order terms. This conjecture could be tested by studying the nonlinear dynamics, which remains governed by second-order field equations.