Higher-Dimensional Vaidya Metric

• The higher-dimensional Vaidya metric is defined as a generalization of the classic spacetime model, incorporating additional dimensions with variable mass functions, electric charge, and matter contributions.
• It employs advanced null and double-null coordinate systems to analyze radiating or accreting black holes, revealing distinct causal structures and horizon behavior.
• Applications span string theory, black hole thermodynamics, and numerical simulations in massive gravity, providing insights into dynamic gravitational collapse and quantum effects.

The higher-dimensional Vaidya metric is the natural generalization of the classic four-dimensional Vaidya spacetime to arbitrary dimension $D = n$ or $D = n+3$, often with inclusion of additional physical effects such as electric charge, cosmological constant, generalized matter sources, and modifications from massive gravity. These metrics model spherically symmetric, radiating or accreting spacetimes, describing gravitational collapse or evaporation processes where the mass function, electric charge, or matter content varies along a privileged null direction. Generalizations to higher dimensions hold particular relevance for research in higher-dimensional gravity, string theory phenomenology, and the study of gravitational collapse and black hole formation in contexts beyond four-dimensional General Relativity.

1. Metric Ansatz and Coordinate Systems

The prototypical higher-dimensional Vaidya metric is constructed on a spacetime that admits a maximally symmetric $(D-2)$-dimensional submanifold (typically, the $S^{D-2}$ sphere) as angular sections. In standard ingoing Eddington–Finkelstein (advanced null) coordinates $(v, r, x^i)$, the line element for $n$-dimensional spacetime with curvature $k$ reads

$$ds^2 = -f(v,r)\,dv^2 + 2\,dv\,dr + r^2\,\gamma_{ij}(x)\,dx^i dx^j,$$

where $\gamma_{ij}$ is the metric on the unit $(n-2)$-sphere or manifold of sectional curvature $k$, and $f(v,r)$ is the lapse function determined by the field equations and matter content. For $D$ spacetime dimensions, the mass term's dependence on $r$ becomes $M(v)/r^{D-3}$, producing distinctive causal structures and horizon behavior that diverge from the four-dimensional case (Han et al., 2023, Chirenti et al., 2012, Hu et al., 2016).

Alternative formulations employ double-null coordinates $(u,v)$ suited for numerical simulations and the analysis of propagation of disturbances: $$ds^2 = -2f(u,v)\,du\,dv + r^2(u,v)\,d\Omega^2_{D-2},$$ where $f(u,v)$ and $r(u,v)$ are smooth functions determined by the generalized Einstein field equations, and $d\Omega^2_{D-2}$ is the metric on the unit $S^{D-2}$ (Chirenti et al., 2012, Chirenti et al., 2011).

2. Mass Functions, Charge, and Matter Content

The classical Vaidya metric models null dust—pure radiation streaming along the $v$-direction—with stress-energy tensor $T_{ab} = \mu(v,r)\, l_a l_b$, where $l_a = \partial_a v$. The generalization to higher dimensions introduces an explicit $r$-dependence in $\mu$: $$\mu(v,r) = \frac{D-2}{16\pi G_D} \frac{dM/dv}{r^{D-2}},$$ where $M(v)$ is the Bondi mass function, now allowed to be an arbitrary monotonic (typically increasing) function of $v$ (Han et al., 2023, Hu et al., 2016, Chirenti et al., 2012).

Further generalizations, such as the charged Vaidya metric, admit time-dependent electric charge $q(v)$ with energy–momentum contributions from the Maxwell sector: $F^{uv} = -F^{vu} = \frac{q(v)}{f\,r^{D-2}},$ and the resulting $r$-dependence in $f(v,r)$ includes terms like $q^2/r^{2(D-3)}$ (Chirenti et al., 2012, Chirenti et al., 2011). Vaidya-like solutions can feature additional nonzero pressure and more generic matter sources, parameterized by supplementary functions $C(v)$ and exponent $o$, which promote $T_{ab}$ from a pure null picture to a mixture including radial pressure and density components (Hu et al., 2016).

A typical form for the lapse function in $D$ dimensions (in the presence of a cosmological constant) is

$$f(v,r) = k + \frac{r^2}{\ell^2} - \frac{16\pi G_D M(v)}{(D-2)\Omega_{D-2}\, r^{D-3}},$$

with extensions for massless or massive gravity, higher curvature, or pressure terms as appropriate (Han et al., 2023, Hu et al., 2016).

3. Field Equations, Energy Conditions, and Horizon Structure

The field equations in higher-dimensional Einstein gravity (possibly with cosmological constant $\Lambda$ and electromagnetic terms) constrain $f(v,r)$ via the Einstein equations with prescribed matter: $$G_{\mu\nu} + \Lambda\,g_{\mu\nu} = 8\pi G_D\, T_{\mu\nu}.$$ In the case of massive gravity, extra interaction terms arising from the dRGT model are included, involving the elementary symmetric polynomials of the square root matrix $\sqrt{g^{-1} f}$ and the singular reference metric $$f_{\mu\nu} = \operatorname{diag}(0,0,c_0^2\gamma_{ij})$$ (Hu et al., 2016).

The admissibility of the Vaidya metric as a solution is controlled by the null and dominant energy conditions. For pure null radiation, $\dot M(v) \ge 0$ ensures the positivity of energy flux. For generalized sources, constraints such as $\mu\ge 0$, $\mu + \rho\ge0$, $P\ge0$ emerge from the respective $T_{ab}$ decomposition (Hu et al., 2016, Chirenti et al., 2012).

Causal structure is governed by the location of marginally trapped surfaces (apparent horizons), whose radii are given by the roots of

$f(v, r_A) = 0,$

or, in double-null form, the fixed-point equation for the outgoing null ray expansion. For charged or cosmological generalizations, this becomes an algebraic equation of higher degree in $r$ that captures the competition between mass, charge, and cosmological repulsion (Chirenti et al., 2012, Shimano et al., 2011). A salient feature at high $D$ is the emergence of new regimes for horizon existence and structure, e.g., suppression of inner horizons unless the charge is sufficiently subdominant relative to mass (Chirenti et al., 2012).

4. Self-Similar Solutions and Trapped Surface Formation

Self-similar higher-dimensional Vaidya spacetimes, defined by homothetic mass functions such as $m(v) = \mu v^n$, possess a conformal Killing vector and model accretion profiles with strict scaling. The existence of closed trapped surfaces, crucial for cosmic censorship, depends on the rise rate of the mass function: a threshold $\frac{dm}{dv} \gtrsim 0.4628$ (numerically) is required for a trapped surface to penetrate into the initially flat region. This generalizes four-dimensional results to arbitrary $D$, and ensures that, provided the threshold is exceeded, no naked singularity is possible, as the entire future domain of dependence inside the trapped surface is causally disconnected from infinity (Shimano et al., 2011).

The maximal possible radius for such a trapped surface asymptotes to the Schwarzschild–Tangherlini horizon radius $r_g = [2M/n]^{1/n}$ as $D \to \infty$, indicating that the quasi-local and global notions of horizon coalesce in the large-$D$ limit (Shimano et al., 2011).

5. Thermodynamics, Misner–Sharp Mass, and Fluctuations

The generalization of the Misner–Sharp quasi-local mass in higher-dimensional Vaidya spacetimes, including massive gravity, yields an effective mass function

$$M_\text{eff}(v,r) = \frac{V_k (D-2)}{16\pi G_D}\, r^{D-3}\left[k + \frac{r^2}{\ell^2} + \text{(higher-order terms)}\right],$$

with corrections determined by the gravity theory and matter content (Hu et al., 2016).

On the apparent horizon, the Clausius relation $\delta Q = T dS$ is exactly satisfied without extra entropy production, evidencing that the system is in thermodynamic equilibrium for all $D$. The unified first law, incorporating potential work terms from pressure, takes the form

$dM_\text{eff} = T\,dS + W\,dV,$

where $W$ is the work density, which typically vanishes for pure null radiation. In massive gravity, these results hold for the full parameter space ensuring real roots to $f=0$, i.e., as long as a horizon exists (Hu et al., 2016).

Metric fluctuations, modeled as small oscillatory components in $M(v)$,

$$M(v) = M_0\left[1 + \epsilon \sin(\omega v)\right],$$

can be treated perturbatively. To $\mathcal{O}(\epsilon^2)$, the horizon radius, area, temperature, and entropy show $D$-dependent correction factors, all of which vanish as $D \to \infty$. The large-$D$ limit thus provides an analytic simplification for the impact of fluctuating horizons on outgoing null rays and for the long-time, averaged black hole thermodynamic variables (Han et al., 2023).

6. Applications: Causal Structure, Shells, and Quasinormal Modes

Higher-dimensional Vaidya spacetimes serve as foundational models for time-dependent black holes in string-theoretic models, brane world cosmology, high-energy astrophysics, and studies of gravitational collapse in higher dimensions. Their analytic tractability applies to:

• Study of dynamical horizons and the process of black hole growth or evaporation under infalling or outgoing massless radiation and/or charged matter (Chirenti et al., 2012, Chirenti et al., 2011).
• Analysis of trapped surfaces and the avoidance of naked singularities in gravitational collapse; threshold results indicate strict bounds on accretion for cosmic censorship (Shimano et al., 2011).
• Investigation of dynamical quasinormal modes and late-time ringdowns; the higher-dimensional form of the wave equation and its coordinate structure allow for numeric and semi-analytic time-domain evolution in nonstationary spacetimes (Chirenti et al., 2011).
• Exploration of quantum gravity effects (as classical proxies), such as fluctuating mass and resulting fluctuating event horizon, allowing computation of induced corrections to black hole thermodynamics and outgoing null geodesics (Han et al., 2023).

7. Extensions: Massive Gravity, Double-Null Formalism, and Open Problems

Recent work extends the higher-dimensional Vaidya metric into the de Rham–Gabadadze–Tolley (dRGT) massive gravity framework. Using a singular reference metric, generalized Vaidya solutions retain the key dynamical features of the standard case, augmented by new polynomial terms with coefficients controlled by the dRGT parameters $c_j$, the reference parameter $c_0$, and the graviton mass $m^2$. The generalized Misner–Sharp mass and first law are explicitly established in this setting, with thermodynamic equilibrium preserved for all allowed parameter choices (Hu et al., 2016).

The double-null formalism provides a powerful framework for both analytic calculations and characteristic numerical evolution of massless and massive field perturbations on dynamical backgrounds, facilitating the computation of apparent and event horizons, tracking null rays, and analyzing stability problems (Chirenti et al., 2012, Chirenti et al., 2011).

Open problems include the full nonlinear stability of higher-dimensional Vaidya spacetimes, the role of higher curvature corrections (Lovelock, Gauss–Bonnet, etc.) in the detailed causal and thermodynamic structure, and the quantum interpretation of time-dependent horizon fluctuations in the context of holography or semiclassical gravity.

---

References:

• S. Hu, X.-M. Wu, and H. Zhang, "Generalized Vaidya Solutions and Misner–Sharp mass for $n$-dimensional massive gravity" (Hu et al., 2016)
• C. Chirenti and A. Saa, "Double null formulation of the general Vaidya metric" (Chirenti et al., 2012)
• C. Chirenti and A. Saa, "Nonstationary regime for quasinormal modes of the charged Vaidya metric" (Chirenti et al., 2011)
• S. Shimano, T. Harada, and M. Tsukamoto, "A trapped surface in the higher-dimensional self-similar Vaidya spacetime" (Shimano et al., 2011)
• T. Koga, H. Maeda, and K. Noda, "Metric fluctuations in higher-dimensional black holes" (Han et al., 2023)