Maxwell-dilaton-dRGT Massive Gravity

• Maxwell-dilaton-dRGT massive gravity is a theory that couples Einstein gravity with a scalar dilaton and Maxwell field, augmented by dRGT mass terms to break spatial diffeomorphisms.
• The model offers rich black hole and black brane solutions, including Lifshitz branes and spherically symmetric black holes with modified thermodynamics and reverse reentrant phase transitions.
• It produces finite DC transport coefficients and observable effects such as altered black hole shadows, providing a bridge between theoretical predictions and astrophysical observations.

Maxwell-dilaton-dRGT-like Massive Gravity refers to a class of gravitational theories that couple Einstein gravity to both a scalar dilaton and a Maxwell $U(1)$ gauge field, and introduce graviton mass terms of de Rham–Gabadadze–Tolley (dRGT) type. These models are constructed to interpolate ultraviolet (dilaton) and infrared (massive gravity) corrections, allowing both for rich black hole and black brane solutions and for controlled breaking of bulk diffeomorphism invariance, with direct implications for holography, black hole thermodynamics, phase transitions, and observable black hole shadows.

1. Gravitational Action and Field Content

The Maxwell-dilaton-dRGT-like massive gravity model is governed by the action in four dimensions

$$S = \frac{1}{16\pi}\int d^4x \sqrt{-g} \left[ \mathcal{R} - 2(\nabla\varphi)^2 - V(\varphi) - F_{\mu\nu}F^{\mu\nu} + e^{-2\beta\varphi} \, m_g^2 \sum_{i=1}^4 \eta_i u_i(g,h) \right]$$

where:

• $\mathcal{R}$: Ricci scalar of $g_{\mu\nu}$,
• $\varphi$: dilaton field with potential $V(\varphi)$,
• $F_{\mu\nu}$: Maxwell tensor,
• $m_g^2 \sum_i \eta_i u_i$: dRGT-type mass term, $u_i(g,h)$ being the $i$-th symmetric polynomial of the eigenvalues of $$S = \frac{1}{16\pi}\int d^4x \sqrt{-g} \left[ \mathcal{R} - 2(\nabla\varphi)^2 - V(\varphi) - F_{\mu\nu}F^{\mu\nu} + e^{-2\beta\varphi} \, m_g^2 \sum_{i=1}^4 \eta_i u_i(g,h) \right]$$0,
• $$S = \frac{1}{16\pi}\int d^4x \sqrt{-g} \left[ \mathcal{R} - 2(\nabla\varphi)^2 - V(\varphi) - F_{\mu\nu}F^{\mu\nu} + e^{-2\beta\varphi} \, m_g^2 \sum_{i=1}^4 \eta_i u_i(g,h) \right]$$1: fixed (reference) metric,
• $$S = \frac{1}{16\pi}\int d^4x \sqrt{-g} \left[ \mathcal{R} - 2(\nabla\varphi)^2 - V(\varphi) - F_{\mu\nu}F^{\mu\nu} + e^{-2\beta\varphi} \, m_g^2 \sum_{i=1}^4 \eta_i u_i(g,h) \right]$$2: dilaton–massive gravity coupling constant,
• $$S = \frac{1}{16\pi}\int d^4x \sqrt{-g} \left[ \mathcal{R} - 2(\nabla\varphi)^2 - V(\varphi) - F_{\mu\nu}F^{\mu\nu} + e^{-2\beta\varphi} \, m_g^2 \sum_{i=1}^4 \eta_i u_i(g,h) \right]$$3: massive gravity parameters, with in 4D only $$S = \frac{1}{16\pi}\int d^4x \sqrt{-g} \left[ \mathcal{R} - 2(\nabla\varphi)^2 - V(\varphi) - F_{\mu\nu}F^{\mu\nu} + e^{-2\beta\varphi} \, m_g^2 \sum_{i=1}^4 \eta_i u_i(g,h) \right]$$4 and $$S = \frac{1}{16\pi}\int d^4x \sqrt{-g} \left[ \mathcal{R} - 2(\nabla\varphi)^2 - V(\varphi) - F_{\mu\nu}F^{\mu\nu} + e^{-2\beta\varphi} \, m_g^2 \sum_{i=1}^4 \eta_i u_i(g,h) \right]$$5 generally nonzero (Heidari et al., 3 Dec 2025, Yue et al., 2024, Kuang et al., 2017).

By a suitable choice of $$S = \frac{1}{16\pi}\int d^4x \sqrt{-g} \left[ \mathcal{R} - 2(\nabla\varphi)^2 - V(\varphi) - F_{\mu\nu}F^{\mu\nu} + e^{-2\beta\varphi} \, m_g^2 \sum_{i=1}^4 \eta_i u_i(g,h) \right]$$6 (typically diagonal with vanishing $$S = \frac{1}{16\pi}\int d^4x \sqrt{-g} \left[ \mathcal{R} - 2(\nabla\varphi)^2 - V(\varphi) - F_{\mu\nu}F^{\mu\nu} + e^{-2\beta\varphi} \, m_g^2 \sum_{i=1}^4 \eta_i u_i(g,h) \right]$$7 and $$S = \frac{1}{16\pi}\int d^4x \sqrt{-g} \left[ \mathcal{R} - 2(\nabla\varphi)^2 - V(\varphi) - F_{\mu\nu}F^{\mu\nu} + e^{-2\beta\varphi} \, m_g^2 \sum_{i=1}^4 \eta_i u_i(g,h) \right]$$8 components), only spatial diffeomorphisms are broken, enabling controlled IR momentum relaxation in the boundary theory (Kuang et al., 2017, Zhou et al., 2015). The coupling $$S = \frac{1}{16\pi}\int d^4x \sqrt{-g} \left[ \mathcal{R} - 2(\nabla\varphi)^2 - V(\varphi) - F_{\mu\nu}F^{\mu\nu} + e^{-2\beta\varphi} \, m_g^2 \sum_{i=1}^4 \eta_i u_i(g,h) \right]$$9 multiplies the graviton mass term and represents a nonminimal interaction between the dilaton and the massive gravity sector.

2. Solutions: Black Branes and Spherically Symmetric Black Holes

The model admits both planar black brane and spherically symmetric black hole solutions:

• Lifshitz black branes with metric ansatz

$\mathcal{R}$0

support arbitrary dynamical exponent $\mathcal{R}$1, with the dilaton and massive gravity coupling encoded in exponential or two-term ansätze for $\mathcal{R}$2 (Kuang et al., 2017).

• Spherically symmetric black holes with

$\mathcal{R}$3

admit explicit analytic solutions for $\mathcal{R}$4, determined by the parameter $\mathcal{R}$5 (dilaton coupling), and a metric function $\mathcal{R}$6 which consists of a sum of power laws in $\mathcal{R}$7, each corresponding to vacuum, cosmological constant, charge, and dRGT-massive gravity terms (Yue et al., 2024, Heidari et al., 3 Dec 2025).

The gauge field acquires a distinctly dilaton-modulated profile, e.g., for $\mathcal{R}$8,

$\mathcal{R}$9

where $g_{\mu\nu}$0 is the electric charge and $g_{\mu\nu}$1 an integration constant.

3. Thermodynamics and Phase Structure

The model's black holes and branes admit a complete set of thermodynamic quantities derived from geometric and Wald/EUclidean techniques:

• Entropy: $g_{\mu\nu}$2 (area law modulated by dilaton)
• Temperature: $g_{\mu\nu}$3
• Mass, charge, potential, and thermodynamic volume: Computed analytically, with modifications from both Maxwell and massive gravity terms
• Equation of state: $g_{\mu\nu}$4, with $g_{\mu\nu}$5 identified with $g_{\mu\nu}$6

Critical point structure is determined by the conditions

$g_{\mu\nu}$7

allowing for multiple real, positive critical radii.

A novel "reverse reentrant" phase transition is triggered when two critical points are present ($g_{\mu\nu}$8), manifesting as an SBH–LBH–SBH transition sequence (small–large–small black hole order) as temperature increases, opposite to the usual reentrant transition. The phase diagram displays a swallowtail in the $g_{\mu\nu}$9–$\varphi$0 plane and a tricritical point where first- and zeroth-order transition curves meet. These effects are fundamentally controlled by the dilaton–massive gravity coupling $\varphi$1 (Yue et al., 2024).

4. DC Transport, Momentum Dissipation, and Holography

In the context of holography, the Maxwell-dilaton-dRGT-like model describes dual quantum systems with nonrelativistic scaling and explicit momentum relaxation:

• Finite DC conductivities: The graviton mass term (with $\varphi$2 breaking translation invariance) ensures finite DC electric, thermoelectric, and thermal conductivities:

$\varphi$3

and similar expressions for $\varphi$4 and $\varphi$5, where $\varphi$6 sets the effective graviton mass squared (Kuang et al., 2017).

• Wiedemann–Franz law violation: The nontrivial dependence of Lorenz ratios $\varphi$7 on $\varphi$8 and $\varphi$9 signals breakdown of the Wiedemann–Franz law, interpreted as a marker of strong interactions in the dual theory.
• Momentum dissipation: The reference metric $V(\varphi)$0 of the form $V(\varphi)$1 in the mass term breaks spatial diffeomorphism invariance in the bulk and relaxes momentum in the boundary field theory, regularizing otherwise divergent transport coefficients.

5. Stability, Phase Transitions, and Geometrothermodynamics

Thermodynamic stability is governed by local (heat capacity) and global (Gibbs free energy) analyses:

• Local stability: Heat capacity at fixed charge ($V(\varphi)$2) or chemical potential ($V(\varphi)$3) — positive values indicate stability.
• Second-order phase transitions: Divergences of $V(\varphi)$4 correspond to continuous (second-order) SBH–LBH transitions; zeros mark physical-limitation points.
• Geometrothermodynamics: The HPEM metric on the $V(\varphi)$5 space produces a Ricci scalar whose divergences exactly match those of $V(\varphi)$6, providing a robust geometric diagnostic for phase transitions (Heidari et al., 3 Dec 2025).

6. Optical Properties: Photon Sphere and Black Hole Shadows

The structure of unstable null orbits determines observable shadows:

• Photon sphere $V(\varphi)$7: Calculated via extrema of the effective potential $V(\varphi)$8, depending on $V(\varphi)$9.
• Black hole shadow radius:

$F_{\mu\nu}$0

For given EHT measurements (e.g., $F_{\mu\nu}$1), the relation $F_{\mu\nu}$2 constrains the gravity and dilaton parameters (Heidari et al., 3 Dec 2025).

Typically, increasing $F_{\mu\nu}$3 or $F_{\mu\nu}$4 enlarges the shadow, while increasing $F_{\mu\nu}$5 reduces its size. This provides observationally testable constraints on the allowed parameter space of Maxwell-dilaton-dRGT-like massive gravity.

7. Energy Emission and Observational Signatures

The Hawking emission spectrum at high frequency asymptotes to $F_{\mu\nu}$6. The peak of the emission rate

$F_{\mu\nu}$7

is modulated by ($F_{\mu\nu}$8): higher values of $F_{\mu\nu}$9 raise and blue-shift the peak, while higher $m_g^2 \sum_i \eta_i u_i$0 lowers and red-shifts it (Heidari et al., 3 Dec 2025). These features could in principle be mapped to observational data from black hole shadow and spectrum measurements.

---