According to Lovelock's theorem, the Hilbert-Einstein and the Lovelock actions are indistinguishable from their field equations. However, they have different scalar-tensor counterparts, which correspond to the Brans-Dicke and the \emph{Lovelock-Brans-Dicke} (LBD) gravities, respectively. In this paper the LBD model of alternative gravity with the Lagrangian density $\mathscr{L}_{\text{LBD}}=\frac{1}{16\pi}\left[\phi\left( R +\frac{a}{\sqrt{-g}}{}^*RR + b\mathcal{G}\right)-\frac{\omega_{\text{L}}}{\phi}\nabla_\alpha \phi \nabla^\alpha\phi \right]$ is developed, where ${}^*RR$ and $\mathcal{G}$ respectively denote the topological Chern-Pontryagin and Gauss-Bonnet invariants. The field equation, the kinematical and dynamical wave equations, and the constraint from energy-momentum conservation are all derived. It is shown that, the LBD gravity reduces to general relativity in the limit $\omega_{\text{L}}\to\infty$ unless the "topological balance condition" holds, and in vacuum it can be conformally transformed into the dynamical Chern-Simons gravity and the generalized Gauss-Bonnet dark energy with Horndeski-like or Galileon-like kinetics. Moreover, the LBD gravity allows for the late-time cosmic acceleration without dark energy. Finally, the LBD gravity is generalized into the Lovelock-scalar-tensor gravity, and its equivalence to fourth-order modified gravities is established. It is also emphasized that the standard expressions for the contributions of generalized Gauss-Bonnet dependence can be further simplified.

The thermodynamics of the Universe is restudied by requiring its compatibility with the holographic-style gravitational equations which govern the dynamics of both the cosmological apparent horizon and the entire Universe, and possible solutions are proposed to the existent confusions regarding the apparent-horizon temperature and the cosmic entropy evolution. We start from the generic Lambda Cold Dark Matter ($\Lambda$CDM) cosmology of general relativity (GR) to establish a framework for the gravitational thermodynamics. The Cai-Kim Clausius equation for the isochoric process of an instantaneous apparent horizon indicates that, the Universe and its horizon entropies encode the \emph{positive heat out} thermodynamic sign convention, which encourages us to adjust the traditional positive-heat-in Gibbs equation into the positive-heat-out version $dE_m=-T_mdS_m-P_mdV$. It turns out that the standard and the generalized second laws (GSLs) of nondecreasing entropies are always respected by the event-horizon system as long as the expanding Universe is dominated by nonexotic matter $-1\leq w_m\leq 1$, while for the apparent-horizon simple open system the two second laws hold if $-1\leq w_m<-1/3$; also, the artificial local equilibrium assumption is abandoned in the GSL. All constraints regarding entropy evolution are expressed by the equation of state parameter, which show that from a thermodynamic perspective the phantom dark energy is less favored than the cosmological constant and the quintessence. Finally, the whole framework is extended from GR and $\Lambda$CDM to modified gravities with field equations $R_{\mu\nu}-Rg_{\mu\nu}/2=8\pi G_{\text{eff}} T_{\mu\nu}^{\text{(eff)}}$. Furthermore, this paper argues that the Cai-Kim temperature is more suitable than Hayward, and the Bekenstein-Hawking and Wald entropies cannot unconditionally apply to the event and particle horizons.

Inspired by the Wald-Kodama entropy $S=A/(4G_{\text{eff}})$ where $A$ is the horizon area and $G_{\text{eff}}$ is the effective gravitational coupling strength in modified gravity with field equation $R_{\mu\nu}-Rg_{\mu\nu}/2=$ $8\pi G_{\text{eff}} T_{\mu\nu}^{\text{(eff)}}$, we develop a unified and compact formulation in which the Friedmann equations can be derived from thermodynamics of the Universe. The Hawking and Misner-Sharp masses are generalized by replacing Newton's constant $G$ with $G_{\text{eff}}$, and the unified first law of equilibrium thermodynamics is supplemented by a nonequilibrium energy dissipation term $\mathcal{E}$ which arises from the revised continuity equation of the perfect-fluid effective matter content and is related to the evolution of $G_{\text{eff}}$. By identifying the mass as the total internal energy, the unified first law for the interior and its smooth transit to the apparent horizon yield both Friedmann equations, while the nonequilibrium Clausius relation with entropy production for an isochoric process provides an alternative derivation on the horizon. We also analyze the equilibrium situation $G_{\text{eff}}=G=\text{constant}$, provide a viability test of the generalized geometric masses, and discuss the continuity/conservation equation. Finally, the general formulation is applied to the FRW cosmology of minimally coupled $f(R)$, generalized Brans-Dicke, scalar-tensor-chameleon, quadratic, $f(R,\mathcal{G})$ generalized Gauss-Bonnet and dynamical Chern-Simons gravity. In these theories we also analyze the $f(R)$-Brans-Dicke equivalence, find that the chameleon effect causes extra energy dissipation and entropy production, geometrically reconstruct the mass $\rho_m V$ for the physical matter content, and show the self-inconsistency of $f(R,\mathcal{G})$ gravity in problems involving $G_{\text{eff}}$.

This paper studies a generic fourth-order theory of gravity with Lagrangian density $f(R,R_c^2,R_m^2, \mathscr{L}_m)$. By considering explicit $R^2$ dependence and imposing the "coherence condition" $f_{R^2}\!=\!f_{R_m^2}\!=\! -f_{R_c^2}/4$, the field equations of $f(R,R^2,R_c^2,R_m^2, \mathscr{L}_m)$ gravity can be smoothly reduced to that of $f(R,\mathcal{G},\mathscr{L}_m)$ generalized Gauss-Bonnet gravity. We use Noether's conservation law to study the $f(\mathcal{R}_1,\mathcal{R}_2\ldots,\mathcal{R}_n,\mathscr{L}_m)$ model with nonminimal coupling between $\mathscr{L}_m$ and Riemannian invariants $ \mathcal{R}_i$, and conjecture that the gradient of nonminimal gravitational coupling strength $\nabla^\mu f_{\!\mathscr{L}_m}$ is the only source for energy-momentum non-conservation. This conjecture is applied to the $f(R,R_c^2,R_m^2, \mathscr{L}_m)$ model, and the equations of continuity and non-geodesic motion of different matter contents are investigated. Finally, the field equation for Lagrangians including the traceless-Ricci square and traceless-Riemann (Weyl) square invariants is derived, the $f(R,R_c^2,R_m^2, \mathscr{L}_m)$ model is compared with the $f(R,R_c^2,R_m^2,T)+2\kappa \mathscr{L}_m$ model, and consequences of nonminimal coupling for black hole and wormhole physics are considered.