Padmanabhan’s Method: Emergent Gravity

• Padmanabhan’s method is a family of constructions that reformulate gravitational and cosmological dynamics via boundary data, horizon quantities, and entropy functionals.
• It encompasses emergent cosmic space laws, entropy-functional variational principles, and duality-invariant propagators to derive modified field equations and address zero-point effects.
• The approach has versatile applications—from reproducing Friedmann dynamics and modifying gravity theories to regularizing black hole solutions and refining gauge field formulations.

Searching arXiv for recent and foundational papers on Padmanabhan-related methods. arxiv.search(query="Padmanabhan emergent space holographic equipartition FRW", max_results=5) Padmanabhan’s method denotes a family of constructions associated with T. Padmanabhan in which gravitational, cosmological, or field-theoretic dynamics are reformulated in terms of boundary data, thermodynamic counting, entropy functionals, or duality-modified propagators. In current usage, the expression most often refers to the emergent-space law for Friedmann cosmology, but the same label is also used for an entropy-functional variational principle for gravity, a duality-invariant propagator implementing a zero-point length, a boundary variational principle, a bulk-surface relation for the Einstein–Hilbert action, and a Born-approximation technique for highly damped quasinormal modes (Ai et al., 2013, Hammad, 2014, Cordoba et al., 28 Jun 2026, Koch, 2024, Smoot, 2012, Xi et al., 2010).

1. Range of meanings

As used across the literature, the term does not identify a single algorithm. It instead names several distinct constructions that share a recurring strategy: bulk dynamics is encoded through horizon quantities, boundary terms, effective degrees of freedom, or short-distance modifications.

Usage of the term	Defining structure	Typical domain
Emergent cosmic space	$\frac{dV}{dt}=L_p^2(N_{\rm sur}-N_{\rm bulk})$	FRW cosmology
Entropy-functional gravity	Extremization of $S[\xi]$ or $S[n^a]$	GR and modified gravity
Duality-invariant propagator	$${\cal G}^{\rm(QG)}_D(s)={\cal G}_D(\sqrt{s^2+\ell^2})$$	Minimal-length QFT
Boundary variational principle	Boundary action with gauge-like variation	Maxwell and Yang–Mills
Bulk-surface relation	Einstein–Hilbert bulk term related to a surface term	Gravitational action
QNM method	Born-amplitude pole extraction	Highly damped quasinormal modes

The emergent-space formulation is the most widely used meaning in cosmology (Ai et al., 2013). The entropy-functional and boundary formulations are closer to emergent-gravity and holographic variational ideas (Tuveri et al., 2016, Koch, 2024). The propagator version belongs to a different line of work, where “Padmanabhan’s method” means imposing a duality on worldline length so that a zero-point length enters the propagator (Cordoba et al., 28 Jun 2026). This polysemy is not accidental; it reflects a broader Padmanabhan program in which spacetime dynamics is repeatedly recast in thermodynamic, boundary, or horizon-centered variables.

2. Emergent cosmic space and holographic equipartition

In cosmology, Padmanabhan’s method treats expansion as the emergence of space driven by a mismatch between surface and bulk degrees of freedom. In the flat FRW case, the central postulate is

$$\frac{dV}{dt}=L_p^2\left(N_{\rm sur}-N_{\rm bulk}\right),$$

with

$$V=\frac{4\pi}{3H^3},\qquad A=4\pi H^{-2},\qquad N_{\rm sur}=\frac{A}{L_p^2},\qquad N_{\rm bulk}=\frac{2|E|}{T},$$

where the bulk energy is the Komar energy $E=|(\rho+3p)|V$ and the horizon temperature is $T=H/(2\pi)$. In this setup, the de Sitter condition $N_{\rm sur}=N_{\rm bulk}$ is holographic equipartition, while the non-equilibrium relation above reproduces the acceleration equation and, with the continuity equation, the Friedmann equation (Ai et al., 2013).

A unifying reformulation was later proposed as

$$\frac{\alpha}{H(n-1)}\frac{dN_{\rm sur}}{dt}=N_{\rm sur}-N_{\rm bulk},$$

with $S[\xi]$0. In that formulation, Padmanabhan’s original law, Cai’s higher-dimensional version, and Sheykhi’s apparent-horizon version are treated as special cases obtained by choosing the horizon radius and entropy formula appropriately (Ai et al., 2013).

A large subsequent literature modifies either the surface count or the bulk count. In Einstein–Cartan cosmology, the bulk contribution is reorganized as

$S[\xi]$1

so that

$S[\xi]$2

with $S[\xi]$3; the spin sector is therefore important mainly in the early universe and cannot replace late-time dark energy (Hadi et al., 2016). In Rastall gravity, the relevant replacement is a modified Komar density,

$S[\xi]$4

together with the non-flat apparent-horizon expansion law

$S[\xi]$5

which reproduces the Rastall Friedmann equations once the modified continuity equation is used (Yuan et al., 2016). In a general braneworld embedding, the extrinsic curvature of the brane contributes an independent degree-of-freedom term,

$S[\xi]$6

where $S[\xi]$7 is assigned to the embedding geometry itself (Heydarzade et al., 2015).

Other extensions alter the horizon entropy rather than the bulk source. In entropy-corrected cosmology one starts from

$S[\xi]$8

uses the resulting effective surface count, and recovers the corresponding modified Friedmann equations; the same strategy has been applied to logarithmic and power-law entropy corrections, as well as to RS II and DGP braneworlds (Sheykhi et al., 2013). A flat-FRW holographic model based simultaneously on a holographic-like connection and Padmanabhan’s law yields the single background equation

$S[\xi]$9

with the notable feature that temperature cancels from the combined evolution law, leaving entropy as the controlling quantity (Komatsu, 2023). A fractional-fractal generalization replaces the smooth Hubble horizon by an effective fractal horizon and obtains

$S[n^a]$0

thereby interpreting emergent space in a fractional-quantum-gravity setting (Junior et al., 2023).

The method has also been compared with entropy bounds. For a flat universe, one study argued that the maximum bulk entropy in Padmanabhan’s emergent universe coincides with Bousso’s covariant entropy bound on the null surface defined by the Hubble horizon when the universe is filled only by a cosmological constant or radiation, while matter domination lowers the maximal area and entropy and motivates a D-bound interpretation (Hadi et al., 2016).

3. Entropy-functional emergent gravity

A second major meaning of Padmanabhan’s method is an entropy-functional variational principle in which spacetime is treated as an elastic medium. In one formulation, the basic variable is a displacement field $S[n^a]$1, and the entropy functional is taken to be quadratic in $S[n^a]$2 and $S[n^a]$3: $S[n^a]$4 Extremality for arbitrary deformations yields Einstein gravity in the simplest case and, when generalized, constrains which higher-curvature theories are compatible with the entropy principle (Hammad, 2014).

In the null-vector version used for GR, one varies

$S[n^a]$5

The variation is with respect to $S[n^a]$6 and the Lagrange multiplier $S[n^a]$7, not the metric. Because the resulting equation must hold for arbitrary null vectors, it leads to Einstein’s equations with cosmological constant, with $S[n^a]$8 appearing as an integration constant rather than as a parameter inserted in the action (Tuveri et al., 2016).

That same line of work was given a covariant reinterpretation through Augmented Variational Principles. For Kerr–Schild metrics over a background metric, the AVP action reduces, at first order and up to a total divergence, to the same bulk functional that Padmanabhan varies with respect to null vectors. This clarifies why the background metric can be treated as non-dynamical in the entropy-functional setup and why the null vector carries the effective dynamical content (Tuveri et al., 2016).

The modified-gravity extension is more restrictive than a generic higher-curvature ansatz. At quadratic order in curvature, the entropy-functional method yields

$S[n^a]$9

rather than an arbitrary $${\cal G}^{\rm(QG)}_D(s)={\cal G}_D(\sqrt{s^2+\ell^2})$$0. The corresponding on-shell entropy agrees with previous Noether-charge and Wald-type results, and the paper emphasizes that the formalism narrows the admissible higher-curvature structure rather than generating all possible modifications (Hammad, 2014).

In $${\cal G}^{\rm(QG)}_D(s)={\cal G}_D(\sqrt{s^2+\ell^2})$$1 gravity, a later paper combined Padmanabhan’s entropy tensor with Hammad’s generalized elasticity ansatz to construct an explicit entropy functional without introducing an auxiliary scalar. The resulting functional splits into a GR-like external part and an extra contribution

$${\cal G}^{\rm(QG)}_D(s)={\cal G}_D(\sqrt{s^2+\ell^2})$$2

That work argues that this entropy-functional extra term is incompatible with the conventional “internal entropy”

$${\cal G}^{\rm(QG)}_D(s)={\cal G}_D(\sqrt{s^2+\ell^2})$$3

and therefore questions whether the latter should be interpreted literally as an entropy (Matouš, 2024).

4. Duality-invariant propagator and zero-point length

In a different branch of the literature, “Padmanabhan’s method” refers to a worldline-duality prescription for propagators. The starting point is the relativistic point-particle path integral, modified so that the proper length $${\cal G}^{\rm(QG)}_D(s)={\cal G}_D(\sqrt{s^2+\ell^2})$$4 appears in the self-dual combination

$${\cal G}^{\rm(QG)}_D(s)={\cal G}_D(\sqrt{s^2+\ell^2})$$5

The corresponding path integral is

$${\cal G}^{\rm(QG)}_D(s)={\cal G}_D(\sqrt{s^2+\ell^2})$$6

and the resulting propagator can be written as

$${\cal G}^{\rm(QG)}_D(s)={\cal G}_D(\sqrt{s^2+\ell^2})$$7

In momentum space,

$${\cal G}^{\rm(QG)}_D(s)={\cal G}_D(\sqrt{s^2+\ell^2})$$8

The method suppresses very short paths, implements a zero-point length, softens short-distance singularities, and preserves Euclidean invariance (Cordoba et al., 28 Jun 2026).

A field-theoretic representation of this propagator was constructed to $${\cal G}^{\rm(QG)}_D(s)={\cal G}_D(\sqrt{s^2+\ell^2})$$9 by expressing the duality-modified $$\frac{dV}{dt}=L_p^2\left(N_{\rm sur}-N_{\rm bulk}\right),$$0-dimensional propagator as an infinite sum of ordinary propagators in dimensions $$\frac{dV}{dt}=L_p^2\left(N_{\rm sur}-N_{\rm bulk}\right),$$1. Truncating at $$\frac{dV}{dt}=L_p^2\left(N_{\rm sur}-N_{\rm bulk}\right),$$2 leads to an effective action

$$\frac{dV}{dt}=L_p^2\left(N_{\rm sur}-N_{\rm bulk}\right),$$3

interpreted as a local free massive scalar theory in $$\frac{dV}{dt}=L_p^2\left(N_{\rm sur}-N_{\rm bulk}\right),$$4, with the last two coordinates constrained by

$$\frac{dV}{dt}=L_p^2\left(N_{\rm sur}-N_{\rm bulk}\right),$$5

The extra dimensions are described there as auxiliary rather than fundamental (Cordoba et al., 28 Jun 2026).

In electrodynamics, the same propagator was used to implement T-duality directly at the level of the photon propagator. The gauge-invariant nonlocal action is

$$\frac{dV}{dt}=L_p^2\left(N_{\rm sur}-N_{\rm bulk}\right),$$6

and the static interaction becomes

$$\frac{dV}{dt}=L_p^2\left(N_{\rm sur}-N_{\rm bulk}\right),$$7

The electric field vanishes at the origin, and the same potential is recovered by a gauge-invariant Hamiltonian treatment based on Dirac-dressed states (Gaete et al., 2022).

The propagator has also been used to generate regular charged black-hole geometries. In that construction, the Newtonian and electrostatic kernels $$\frac{dV}{dt}=L_p^2\left(N_{\rm sur}-N_{\rm bulk}\right),$$8 are replaced by $$\frac{dV}{dt}=L_p^2\left(N_{\rm sur}-N_{\rm bulk}\right),$$9, which produces smeared matter and electromagnetic energy densities, a regular mass function, and a static metric that is asymptotically Reissner–Nordström but de Sitter-like at the core. The paper notes that the static solution resembles the Ayón-Beato–García spacetime provided the T-duality scale is redefined as $$V=\frac{4\pi}{3H^3},\qquad A=4\pi H^{-2},\qquad N_{\rm sur}=\frac{A}{L_p^2},\qquad N_{\rm bulk}=\frac{2|E|}{T},$$0 and an additional correction function is set to one (Gaete et al., 2022).

5. Boundary variational principle

A further meaning of Padmanabhan’s method is a boundary variational principle in which bulk field equations are derived from an action composed primarily of a boundary term. The essential feature is the choice of variation: one varies the gauge field in a way that looks like a gauge transformation, but matter fields are held fixed, so the operation is not treated as a true gauge transformation of the full system (Koch, 2024).

For electrodynamics, the starting boundary action is

$$V=\frac{4\pi}{3H^3},\qquad A=4\pi H^{-2},\qquad N_{\rm sur}=\frac{A}{L_p^2},\qquad N_{\rm bulk}=\frac{2|E|}{T},$$1

combined with the interaction term

$$V=\frac{4\pi}{3H^3},\qquad A=4\pi H^{-2},\qquad N_{\rm sur}=\frac{A}{L_p^2},\qquad N_{\rm bulk}=\frac{2|E|}{T},$$2

With the Padmanabhan variation

$$V=\frac{4\pi}{3H^3},\qquad A=4\pi H^{-2},\qquad N_{\rm sur}=\frac{A}{L_p^2},\qquad N_{\rm bulk}=\frac{2|E|}{T},$$3

one has $$V=\frac{4\pi}{3H^3},\qquad A=4\pi H^{-2},\qquad N_{\rm sur}=\frac{A}{L_p^2},\qquad N_{\rm bulk}=\frac{2|E|}{T},$$4, and the total variation reduces to the boundary expression

$$V=\frac{4\pi}{3H^3},\qquad A=4\pi H^{-2},\qquad N_{\rm sur}=\frac{A}{L_p^2},\qquad N_{\rm bulk}=\frac{2|E|}{T},$$5

Because $$V=\frac{4\pi}{3H^3},\qquad A=4\pi H^{-2},\qquad N_{\rm sur}=\frac{A}{L_p^2},\qquad N_{\rm bulk}=\frac{2|E|}{T},$$6 is arbitrary on the boundary, stationarity gives the contracted equation

$$V=\frac{4\pi}{3H^3},\qquad A=4\pi H^{-2},\qquad N_{\rm sur}=\frac{A}{L_p^2},\qquad N_{\rm bulk}=\frac{2|E|}{T},$$7

For sufficiently general boundaries, this implies the full Maxwell equation $$V=\frac{4\pi}{3H^3},\qquad A=4\pi H^{-2},\qquad N_{\rm sur}=\frac{A}{L_p^2},\qquad N_{\rm bulk}=\frac{2|E|}{T},$$8 (Koch, 2024).

The same logic was extended to Yang–Mills theory. There the Padmanabhan variation is the infinitesimal gauge transformation of the gauge field alone,

$$V=\frac{4\pi}{3H^3},\qquad A=4\pi H^{-2},\qquad N_{\rm sur}=\frac{A}{L_p^2},\qquad N_{\rm bulk}=\frac{2|E|}{T},$$9

and the resulting contracted boundary equation is

$E=|(\rho+3p)|V$0

equivalently the normal projection of $E=|(\rho+3p)|V$1 (Koch, 2024).

The paper distinguishes between an “ASC” boundary class, where arbitrary local normals force the standard bulk equations, and an “IF” boundary class, where the contracted equations allow additional structures. In the latter case, the weak-field equations can acquire terms such as

$E=|(\rho+3p)|V$2

and

$E=|(\rho+3p)|V$3

interpreted there as analogs of CP-violating dual couplings and gauge-boson mass terms. The paper draws an explicit analogy with Padmanabhan’s gravitational result that the cosmological constant appears as an integration constant rather than as a fundamental coupling (Koch, 2024).

6. Related reinterpretations, specialized uses, and open issues

Several papers use “Padmanabhan’s method” in more specialized senses. One concerns the differential bulk-surface relation of the Einstein–Hilbert action. In a differential-form gauge-theory reformulation due to Göckeler and Schücker, the Einstein–Hilbert action is split into a bulk term and an exact boundary term, and the paper argues that Padmanabhan’s bulk-surface relation can be rederived naturally in that language (Smoot, 2012). Another reinterpretation constructs Padmanabhan’s surface density of spacetime degrees of freedom as an acceleration-driven phase-space density: using a spacelike foliation along the acceleration direction and projecting the canonical pair along the observer’s velocity, one obtains

$E=|(\rho+3p)|V$4

after integrating over one Euclidean-time period (Hadad, 2016).

The term also appears in black-hole perturbation theory. For highly damped quasinormal modes of a black hole with deficit solid angle and quintessence-like matter, Padmanabhan’s Born-approximation scattering method yields evenly spaced imaginary parts,

$E=|(\rho+3p)|V$5

with spacing determined only by the black-hole horizon surface gravity $E=|(\rho+3p)|V$6 (Xi et al., 2010).

The cosmological literature also contains attempts to reinterpret the driving term $E=|(\rho+3p)|V$7 thermodynamically. In one such proposal, a generalized entropy $E=|(\rho+3p)|V$8 makes entropy proportional to volume and leads to

$E=|(\rho+3p)|V$9

suggesting a bridge between Verlinde’s entropic logic and Padmanabhan’s emergent-space law. The same paper, however, argues that this generalized entropy is incompatible with the unified first law and the standard Friedmann equation unless dark-energy-induced entropy corrections are included (Moradpour, 2016). More generally, several extensions of the emergent-space law are presented as phenomenological or structurally motivated rather than microscopically derived; one unified formula for FRW dynamics explicitly notes that its deeper origin was not found (Ai et al., 2013). A similar caution appears in the boundary variational formulation, where the resulting bulk equations depend on the class of allowed boundaries (Koch, 2024).

Taken together, these usages make Padmanabhan’s method less a single technique than a methodological family. The common thread is a persistent replacement of standard bulk-first dynamics by boundary, horizon, entropy, equipartition, or ultraviolet-duality data. The major points of continuing debate concern how much of this structure is fundamental, which pieces are heuristic, and when thermodynamic terms introduced in extended theories are genuinely entropic rather than effective work or pressure contributions (Matouš, 2024, Moradpour, 2016).