Black holes in the four-dimensional Einstein-Lovelock gravity

A $(3+1)$-dimensional Einstein-Gauss-Bonnet theory of gravity has been recently formulated in [D. Glavan and C. Lin, Phys. Rev. Lett. {\bf 124}, 081301 (2020)] which is different from the pure Einstein theory, i.e., bypasses the Lovelock's theorem and avoids Ostrogradsky instability. The theory was formulated in $D > 4$ dimensions and its action consists of the Einstein-Hilbert term with a cosmological constant, while the Gauss-Bonnet term multiplied by a factor $1/(D-4)$. Then, the four-dimensional theory is defined as the limit $D \to 4$. Here we generalize this approach to the four-dimensional Einstein-Lovelock theory and formulate the most general static $4D$ black-hole solution allowing for a $\Lambda$-term (either positive or negative) and the electric charge $Q$. As metric functions cannot be found in a closed form in the general case, we develop and share publicly the code which constructs the metric functions for every given set of parameters.