The idea that spacetime may possess more than the three spatial dimensions we experience is not new. Theories such as string theory naturally require ten or eleven dimensions for mathematical consistency, with the extra dimensions compactified—folded into tiny, closed geometries far below any scale we can probe with ordinary instruments. The Kaluza–Klein framework, dating back nearly a century, showed that when extra dimensions are compactified, gravity leaks into them in a very particular way: rather than spreading in all dimensions, gravitational influence appears to us as the familiar inverse-square law, modified at very short distances by exponentially decaying corrections. Laboratory tests of Newton’s law at submillimeter scales, including torsion-balance experiments probing separations down to tens of micrometers, have so far found no departure from the standard behavior, pushing the compactification scale to increasingly small values. But gravitational waves, which propagate through spacetime itself, offer an entirely independent and complementary probe.

Villani’s approach begins with the fundamental problem of gravitational radiation in such a higher-dimensional setting. In standard general relativity, the emission of gravitational waves from a binary system is described by solving the wave equation for metric perturbations around flat spacetime, using the post-Newtonian expansion—a systematic approximation scheme in powers of v/c, the orbital velocity relative to the speed of light. Decades of work by theorists such as Luc Blanchet and collaborators have refined this framework to remarkably high precision, and the resulting waveforms match the signals detected by the LIGO, Virgo, and KAGRA observatories to extraordinary accuracy. Any modification to this framework would need to preserve these successes while introducing signatures in regimes not yet tested.

In a spacetime with compactified extra dimensions, the situation changes fundamentally because the gravitational field can no longer be described by a single massless mode propagating in three extended spatial dimensions. When the metric is expanded in the compact directions, the graviton—the fundamental quantum of gravitational radiation—gives rise to a tower of additional modes. Villani solves both the homogeneous and the non-homogeneous wave equations in this setting and demonstrates a striking result: the radiated field consists not only of the usual massless mode, which falls off as one over the distance from the source, but also of an infinite sum of so-called pseudo-massive modes. Each of these additional modes behaves as if it had a tiny effective mass associated with its excitation number along the compact dimensions, and each decays exponentially with distance rather than spreading freely through space.

This exponential decay is the key physical feature. A pseudo-massive gravitational mode with effective mass m_d carries a characteristic length scale—the compactification radius of the extra dimension—beyond which its contribution becomes negligible. Near the source, however, these modes are present and they modify the gravitational field in ways that feed back into the orbital dynamics of the binary itself. Villani calculates the metric potentials generated by the binary, derives the resulting equations of motion, and then computes the energy flux carried away by the complete radiated field, including both the massless mode and the tower of pseudo-massive contributions.

One of the most technically significant findings of the study concerns the post-Newtonian structure of the equations of motion. In the standard framework, corrections to Newtonian gravity appear at successive orders: the leading 1PN correction scales as the inverse square of the speed of light, the 2PN correction as the inverse fourth power, and so on. Villani finds that in a spacetime with compactified extra dimensions, a new term appears at the so-called -1PN order—that is, an order below Newtonian in the standard hierarchy, multiplying the Newtonian acceleration itself. Physically, this term arises from the exponential screening of the gravitational interaction by the compactified geometry: the effective Newtonian coupling is multiplied by a factor that depends on the separation between the bodies relative to the compactification scale, so that gravity appears slightly weaker at separations approaching that scale. This unusual ordering, in which an effect appears “before” the Newtonian term in the post-Newtonian counting, reflects the fact that extra-dimensional corrections do not come from relativistic refinements of gravity but from a modification of the underlying force law itself.

The consequences for gravitational-wave astronomy follow directly. Because the energy flux determines how rapidly a binary inspirals—losing orbital energy to radiation and drawing the two objects closer together at an accelerating rate—any modification of the flux changes the predicted waveform. Detectors such as LIGO and Virgo measure not merely the existence of a merger but the detailed phase evolution of the signal as the frequency sweeps upward through the sensitive band. A small additional contribution to the energy flux, or a small alteration of the binding energy, accumulates coherently over hundreds or thousands of orbital cycles, producing a measurable shift in the arrival time of each wave crest. The growing catalog of compact binary coalescences, now numbering in the hundreds of events including the most recent GWTC-4.0 observations, therefore constitutes a dataset of extraordinary sensitivity to exactly the kind of corrections Villani’s calculation predicts.

The study builds on a growing body of work connecting gravitational-wave observations to extra dimensions. Previous research has shown that gravitational-wave detectors could, in principle, act as probes of compactified spacetime geometry, and that the propagation of waves over cosmological distances—so-called standard sirens—can constrain whether gravity “leaks” into large extra dimensions, testing deviations from the inverse-square behavior at astronomical scales. Villani’s contribution is distinct in focusing on the generation rather than the propagation of the waves: he works out how the emission process itself is altered when the gravitational field near the binary is modified by compactified dimensions. This source-side calculation complements propagation-side analyses and could eventually be incorporated into waveform templates used in matched-filter searches, where even minute theoretical mismodeling can reduce detection efficiency or bias parameter estimates.

The mathematical machinery required for the calculation is considerable. Villani expands the metric perturbations in tensor spherical harmonics generalized to higher-dimensional spheres, building on foundational work on symmetric tensor harmonics and their eigenvalue structures on n-dimensional manifolds. The source term of the wave equation—describing the stress-energy of the two orbiting masses—must likewise be decomposed in this higher-dimensional harmonic basis. The resulting expressions for the Green’s functions linking source to radiation involve special functions including exponential integrals, and the author notes that the symbolic computations were carried out with the aid of computer algebra, since the general expressions cannot be manipulated analytically for arbitrary mode numbers. A detailed appendix derives the auxiliary functions that enter the metric potentials, showing how the exponential decay factors of the pseudo-massive modes combine with the source geometry to produce the final flux formula.

Importantly, the framework passes a crucial consistency check. When the compactification radius is taken to infinity’s counterpart—the limit in which the effects of the extra dimensions vanish—and the number of compact dimensions goes to zero, all of Villani’s expressions reduce smoothly to the standard post-Newtonian results found in the canonical literature on compact binary dynamics. The exponentially decaying modes disappear, the -1PN correction vanishes, and the energy flux returns to the well-verified form used in contemporary waveform models. This means the new calculation does not contradict existing gravitational-wave physics; it extends it, embedding the familiar four-dimensional theory inside a larger framework from which it emerges as a limiting case.

What would it take to see the new effects? The exponentially decaying character of the pseudo-massive modes means their reach is governed by the compactification scale, and all current bounds—from submillimeter torsion experiments to astrophysical constraints—indicate that this scale, if extra dimensions exist at all, lies well below the orbital separations of LIGO’s stellar-mass binaries at most frequencies. However, the -1PN correction modifies the effective gravitational constant even at orbital scales through its influence on the binding energy and flux at higher orders, and the analysis makes precise how such corrections would enter the phase of the observed signal. As detector sensitivity improves with future observing runs and next-generation instruments, and as the event count grows, statistical combination of many observations could tighten constraints on any deviation from pure general relativity, indirectly testing the compactified extra-dimensional scenario described here.

The work also carries implications beyond binary mergers. Compactified extra dimensions feature prominently in models of the early universe, where cosmological phase transitions in warped geometries are predicted to generate stochastic backgrounds of gravitational radiation. The same formalism that describes how a binary source radiates in such a spacetime can inform how those primordial signals are generated, shaped, and possibly detected. And on the theoretical side, the demonstration that pseudo-massive modes inevitably accompany gravitational radiation in compactified geometries adds a concrete element to our understanding of how gravity behaves when the dimensionality of spacetime is larger than it appears.

For now, the result is theoretical—no datasets were generated or analyzed, and no observational claim of extra dimensions is made. But it supplies exactly what waveform modelers need: an explicit, controlled calculation of source dynamics and radiation in a well-motivated extension of general relativity, with clear limiting behavior and a novel post-Newtonian structure. As gravitational-wave astronomy matures from a detection science into a precision science, calculations of this kind define the frontier—mapping, signal by signal, how much of the universe’s hidden geometry the ripples of spacetime can reveal.