On one hand, since the field equations of gravity are generally covariant and of second order derivatives in the metric tensor, one would
naively expect these equations to be derived from an action principle involving metric tensor and its first and second order derivatives [2],
analogous to the situation for many other field theories of physics. In recent years, there have been considerable works for understanding
the role of the higher curvature terms from various points of view, especially with regard to higher-dimensional black hole physics. For
example, thermodynamics and other properties of the static black hole solutions in Gauss–Bonnet gravity have been found by many
authors [3–9]. Also, the Taub–NUT/bolt solutions of higher derivative gravity and their thermodynamics properties have been constructed
[10–13]. Not long ago, M.H. Dehghani introduced two new classes of rotating solutions of second order Lovelock gravity and investigated
their thermodynamics [14,15].