Thermal aspects and particle dynamics of Euler-Heisenberg AdS black hole in 4D Einstein Gauss-Bonnet gravity

1. Problem Statement

The paper addresses the need to understand the interplay between higher-curvature gravity and nonlinear electrodynamics (NLED) in the context of black hole (BH) physics. Specifically, it investigates charged Anti-de Sitter (AdS) black holes within the framework of 4D Einstein-Gauss-Bonnet (EGB) gravity coupled to Euler-Heisenberg (EH) nonlinear electrodynamics.

• Context: Standard General Relativity (GR) and Maxwell electrodynamics often fail to describe extreme conditions (e.g., near singularities or strong fields).
• Motivation:
  • EGB Gravity: A second-order Lovelock theory that avoids ghost instabilities typical of higher-derivative theories. While the Gauss-Bonnet (GB) term is topological in 4D, a regularization scheme (Glavan-Lin) allows it to contribute dynamically to the field equations.
  • Euler-Heisenberg NLED: An effective Lagrangian derived from Quantum Electrodynamics (QED) describing vacuum polarization and pair production in strong fields. It corrects the linear Maxwell theory at high field strengths.
• Goal: To construct exact BH solutions in this combined framework and analyze their horizon structure, thermodynamic properties (including Joule-Thomson expansion), and the dynamics of test particles (geodesics).

2. Methodology

A. Theoretical Framework

The authors start with the action for $D$-dimensional EGB gravity coupled to EH NLED:
$$S = \frac{1}{16\pi} \int d^D x \sqrt{-g} \left[ R - \frac{(D-1)(D-2)}{3}\Lambda - 4L_{EH}(F, G) + \alpha L_{GB} \right]$$

• Lagrangians: $L_{GB}$ is the Gauss-Bonnet term ($R^2 - 4R_{ab}R^{ab} + R_{abcd}R^{abcd}$). $L_{EH}$ is the Euler-Heisenberg Lagrangian, expanded to first order in the coupling parameter $a$: $L_{EH} = -F + \frac{a}{2}F^2 + \frac{7a}{8}G^2$.
• Regularization: To obtain non-trivial 4D solutions, the GB coupling $\alpha$ is rescaled as $\alpha \to \alpha/(D-4)$ before taking the limit $D \to 4$.
• Metric Ansatz: A static, spherically symmetric metric is assumed: $ds^2 = -F(r)dt^2 + F(r)^{-1}dr^2 + r^2 d\Omega^2$.
• Solution: By solving the modified field equations for a purely electric charge $Q$, the authors derive the metric function $F(r)$ (Eq. 18), which includes contributions from mass $M$, charge $Q$, cosmological constant $\Lambda$, GB coupling $\alpha$, and EH parameter $a$.

B. Thermodynamic Analysis

The study is conducted in the extended phase space, where the cosmological constant is treated as thermodynamic pressure ($P = -\Lambda/8\pi$).

• Quantities Derived: Mass (interpreted as Enthalpy $H$), Hawking Temperature ($T_H$), Entropy ($S$), and Heat Capacity ($C_P$).
• Equation of State: The relationship between $P$, $T$, and horizon radius $r_+$ is derived.
• Joule-Thomson (JT) Expansion: The authors calculate the JT coefficient $\mu_J = (\partial T / \partial P)_H$ to determine inversion curves separating heating and cooling regions.

C. Geodesic Dynamics

The motion of test particles is analyzed using the Lagrangian formalism.

• Effective Potential: Derived for massive particles to study orbital stability.
• Key Metrics: Circular orbit radii, Keplerian frequency ($\Omega_K$), epicyclic frequencies ($\Omega_r, \Omega_\theta$), and the Innermost Stable Circular Orbit (ISCO).

3. Key Contributions and Results

A. Horizon Structure

• The metric function admits multiple horizons (event, Cauchy, and cosmological) depending on parameters.
• Branch Selection: The "minus" branch of the solution is identified as the physically meaningful one, reducing to the standard Reissner-Nordström-AdS solution when $\alpha \to 0$ and $a \to 0$.
• Parameter Effects: Both $\alpha$ and $a$ significantly alter the horizon structure, allowing for multiple positive real roots under specific conditions.

B. Thermodynamics

• Entropy: The entropy acquires a logarithmic correction ($S = \pi r_+^2 + 4\alpha \ln r_+$) due to the GB term, deviating from the standard area law.
• Critical Behavior: The system exhibits a liquid-gas phase transition.
  • GB Coupling ($\alpha$): Increasing $\alpha$ raises the critical pressure ($P_c$) and critical temperature ($T_c$) while increasing the critical radius ($r_c$). It dominates the shift in critical points.
  • EH Parameter ($a$): Increasing $a$ slightly decreases $r_c$ and increases $T_c$ and $P_c$, but its effect is comparatively weaker than that of $\alpha$.
• Stability: The heat capacity $C_P$ shows divergence at a specific radius, marking a second-order phase transition between unstable (small) and stable (large) black hole phases.

C. Joule-Thomson Expansion

• Inversion Curves: The study identifies non-trivial inversion curves in the $T-P$ plane, separating heating ($\mu_J < 0$) and cooling ($\mu_J > 0$) regions.
• Parameter Influence:
  • Charge ($Q$): Higher charge increases the inversion temperature and pressure, significantly expanding the cooling region.
  • EH Parameter ($a$): Increasing $a$ shifts the inversion curve to lower temperatures and pressures, reducing the cooling region. This highlights the screening effect of nonlinear electrodynamics.
  • GB Coupling ($\alpha$): Increasing $\alpha$ shifts the inversion curve to higher pressures, indicating that higher-curvature effects modify the equation of state and the balance between gravity and pressure.

D. Particle Dynamics (Geodesics)

• Effective Potential ($V_{eff}$):
  • GB Effect: Increasing $\alpha$ raises the potential barrier and shifts the maximum outward, requiring higher energy for particles to approach the horizon.
  • EH Effect: Increasing $a$ lowers the potential barrier, facilitating inward motion.
• Orbital Stability:
  • Stable Orbits: Increasing $\alpha$ shifts the radius of stable circular orbits outward, expanding the region of stability. Conversely, increasing $a$ reduces the stable orbit radius.
  • ISCO: The ISCO radius decreases with increasing $\alpha$ (allowing stable orbits closer to the horizon) but is minimally affected by $a$.
• Epicyclic Frequencies: The radial epicyclic frequency squared ($\Omega_r^2$) becomes negative in certain regions, defining the ISCO. The GB term pushes the ISCO outward, while the EH term pulls it inward.

4. Significance and Conclusion

This paper provides a comprehensive analysis of how quantum corrections (via EH NLED) and higher-curvature gravity (via EGB) jointly influence black hole physics.

• Physical Insight: The study demonstrates that while both corrections modify thermodynamic and dynamic properties, they often act in competition. The GB term tends to dominate the short-distance behavior and critical thermodynamic shifts, whereas the EH term introduces non-Maxwellian screening effects that alter the cooling/heating dynamics and orbital stability.
• Observational Relevance: The modifications to the ISCO, Keplerian frequency, and phase transition points offer potential observational signatures for distinguishing these theories from standard GR.
• Future Directions: The authors suggest extending this work to rotating solutions, quasinormal modes, and black hole shadows to further probe the strong-field regime of 4D EGB gravity coupled to QED effects.

In summary, the work establishes that the combined framework of 4D EGB gravity and Euler-Heisenberg electrodynamics yields a rich phase structure and distinct dynamical behaviors, offering a viable theoretical laboratory for testing deviations from General Relativity in extreme environments.