5D Proca-Maxwell System

• The five-dimensional Proca-Maxwell system is defined by a spin-1 theory where a massive Proca field is minimally coupled to a gauge-invariant Maxwell field.
• It employs higher-derivative gravity with an infinite tower of curvature corrections to achieve regular spacetime cores and distinctive frozen state dynamics.
• Numerical and analytical studies reveal charge effects, quasi-horizons, and stability criteria, with implications across flat, AdS, and compactified frameworks.

The five-dimensional Proca-Maxwell system denotes a class of spin-1 field theories in $D=5$ in which a massive vector field and an electromagnetic field are coupled on a five-dimensional background, with the Proca sector typically breaking $U(1)$ gauge invariance while the Maxwell sector remains gauge invariant. In the formulation studied most explicitly in higher-derivative gravity, the matter content is a complex gauged Proca field $B_\mu$ minimally coupled to a Maxwell field $A_\mu$, and the full system is embedded in a five-dimensional gravitational theory containing Einstein-Hilbert gravity plus an infinite tower of generalized quasi-topological curvature corrections (Hao et al., 3 Apr 2026). In the broader five-dimensional literature, the same terminology also covers flat-space Proca and Maxwell dynamics, Kaluza-Klein compactifications, AdS normal and quasinormal modes, dual formulations, and generalized Proca black-hole solutions (Delphenich, 2018, Escalante et al., 2014, Lopes et al., 2024, Lopes et al., 18 May 2026, Dalmazi et al., 2011, Hassaine et al., 18 May 2026).

1. Field content and defining structures

In the higher-derivative gravitating model, spacetime dimension is $D=5$, the metric ansatz is static and spherically symmetric with $S^3$ angular sections, and the matter sector consists of a complex Proca field $B_\mu$ of mass $m$ and gauge charge $q$, together with a Maxwell potential $A_\mu$ coupled through the gauge-covariant derivative $U(1)$0 (Hao et al., 3 Apr 2026). The Proca field strength is

$U(1)$1

while the Maxwell field strength is

$U(1)$2

The specific five-dimensional gravitating configurations studied numerically use

$U(1)$3

together with the Proca ansatz

$U(1)$4

and the Maxwell potential

$U(1)$5

The frequency $U(1)$6 controls the harmonic time dependence of the Proca condensate, and the combination $U(1)$7 enters the reduced field equations (Hao et al., 3 Apr 2026).

Outside the higher-derivative setting, five-dimensional Proca-Maxwell theory is also studied in fixed backgrounds. In flat five-dimensional Minkowski space one may regard the Maxwell equations as the fundamental five-dimensional system and recover four-dimensional Proca dynamics by separating the fifth coordinate, with the identification $U(1)$8 under a five-dimensional null dispersion relation $U(1)$9 (Delphenich, 2018). In AdS$B_\mu$0, Proca and Maxwell fields admit a scalar-type/vector-type decomposition on $B_\mu$1, with $B_\mu$2 physical Proca degrees of freedom and $B_\mu$3 Maxwell degrees of freedom in five dimensions (Lopes et al., 2024, Lopes et al., 18 May 2026). These formulations are mathematically distinct but share the same basic distinction: the Proca sector is massive and non-gauge, whereas the Maxwell sector is massless and gauge-invariant.

2. Action, equations of motion, and higher-curvature sector

The gravitating model is defined by an action split into gravitational and matter parts,

$B_\mu$4

where the higher-curvature terms $B_\mu$5 form an infinite tower of generalized quasi-topological invariants (Hao et al., 3 Apr 2026). The second-order term is the Gauss-Bonnet density,

$B_\mu$6

A convenient coupling choice is

$B_\mu$7

which leads, in static spherical symmetry, to the compact function

$B_\mu$8

with

$B_\mu$9

This $A_\mu$0 encodes the full higher-derivative tower under the spherical ansatz (Hao et al., 3 Apr 2026).

The matter Lagrangians are

$A_\mu$1

The Proca field equation is

$A_\mu$2

and the Maxwell equation with Proca current source is

$A_\mu$3

The cross-coupling is minimal, and there are no additional non-minimal matter-curvature couplings in this model (Hao et al., 3 Apr 2026).

For comparison, the corresponding fixed-background Proca equation in AdS$A_\mu$4 is

$A_\mu$5

which implies $A_\mu$6, while the Maxwell limit is obtained as $A_\mu$7 with $A_\mu$8 (Lopes et al., 2024, Lopes et al., 18 May 2026). In flat-space five-dimensional compactifications, the same mass term appears after reducing the five-dimensional Maxwell action on a compact dimension, yielding a four-dimensional effective theory containing a massive vector zero mode, a tower of massive KK vectors, and a tower of massive scalars $A_\mu$9 (Escalante et al., 2014).

3. Spherical reduction, boundary data, and conserved quantities

Under the static $D=5$0-symmetric ansatz, the full field equations reduce to an ODE system for $D=5$1, $D=5$2, $D=5$3, $D=5$4, and $D=5$5. The first two reduced equations are

$D=5$6

and

$D=5$7

supplemented by one algebraic Proca equation and two second-order matter equations (Hao et al., 3 Apr 2026).

Regularity at the origin requires

$D=5$8

while asymptotic flatness requires

$D=5$9

with $S^3$0 fixed by gauge and set to zero in practice (Hao et al., 3 Apr 2026).

The ADM mass is read from

$S^3$1

and the Proca Noether charge is

$S^3$2

The electric charge is

$S^3$3

and the binding energy is

$S^3$4

Positive $S^3$5 suggests gravitationally bound configurations, while $S^3$6 typically indicates instability (Hao et al., 3 Apr 2026).

The stress tensor is anisotropic. With $S^3$7, $S^3$8, and $S^3$9, one has $B_\mu$0, and the explicit expressions used in numerics depend on $B_\mu$1 and $B_\mu$2 (Hao et al., 3 Apr 2026). Curvature regularity is monitored through the Kretschmann scalar

$B_\mu$3

whose explicit form under the ansatz involves $B_\mu$4, $B_\mu$5, $B_\mu$6, $B_\mu$7, and $B_\mu$8 (Hao et al., 3 Apr 2026).

4. Numerical construction and higher-derivative solution space

The numerical analysis is performed with the dimensionless choices $B_\mu$9 and $m$0, the compactified radial coordinate

$m$1

a finite-element discretization with 10,000 grid points, and Newton-Raphson iteration with relative error tolerance below $m$2 (Hao et al., 3 Apr 2026).

The correction order $m$3 organizes the gravitational sector. The case $m$4 recovers Einstein gravity, $m$5 adds Gauss-Bonnet, and $m$6 resums the full local tower via $m$7 (Hao et al., 3 Apr 2026). The solution space changes qualitatively with $m$8 and $m$9.

In the Einstein baseline, the five-dimensional charged Proca-star branch displays the familiar spiral in $q$0 and $q$1, and for increasing charge $q$2 the available frequency interval shrinks drastically, down to approximately $q$3 for $q$4. Along this charged Einstein branch, the binding energy is negative across the branch, indicating likely instability (Hao et al., 3 Apr 2026).

In the Gauss-Bonnet case, strong coupling such as $q$5 removes the spiral and produces a monotonic branch from $q$6 downwards, but in the neutral low-frequency limit the matter fields concentrate at the center, the Kretschmann scalar diverges, and the interior metric develops a pathology. Gauss-Bonnet alone therefore fails to regularize the core as $q$7 (Hao et al., 3 Apr 2026).

By contrast, higher-order corrections and especially the nonperturbative $q$8 resummation regularize the spacetime. In vacuum, the $q$9 metric is

$A_\mu$0

which behaves near the origin as

$A_\mu$1

corresponding to a smooth de Sitter core rather than a curvature singularity (Hao et al., 3 Apr 2026).

5. Frozen states, quasi-horizons, and the effect of electric charge

The central structural result in the infinite-tower theory is the appearance of a neutral “frozen state” in the supercritical regime. As $A_\mu$2 with $A_\mu$3, the matter fields retract completely inside a critical radius $A_\mu$4, and the metric develops a global minimum approaching zero at $A_\mu$5 without ever forming an actual horizon (Hao et al., 3 Apr 2026). The exterior becomes indistinguishable from an extremal vacuum black-hole geometry.

For the neutral $A_\mu$6 vacuum exterior, one may write

$A_\mu$7

The quasi-horizon location and minimum metric value satisfy

$A_\mu$8

and near extremality,

$A_\mu$9

These relations provide analytic benchmarks for the numerical frozen branch (Hao et al., 3 Apr 2026).

A representative numerical example at $U(1)$00, $U(1)$01, and $U(1)$02 gives $U(1)$03, $U(1)$04, $U(1)$05, and $U(1)$06, very close to $U(1)$07. The energy density is entirely confined within $U(1)$08, the exterior matches the extremal vacuum metric, and curvature invariants remain finite everywhere (Hao et al., 3 Apr 2026).

Introducing electric charge qualitatively changes this low-frequency behavior. The Maxwell sector enters through $U(1)$09 and through the source $U(1)$10, and the paper identifies the physical mechanism as long-range Coulomb repulsion counteracting compression and preventing access to $U(1)$11 (Hao et al., 3 Apr 2026). This “unfreezing” produces a strictly positive lower bound $U(1)$12 that grows monotonically with $U(1)$13.

At $U(1)$14 and $U(1)$15, the charged branch shrinks with increasing $U(1)$16. For $U(1)$17, one finds $U(1)$18, $U(1)$19, and $U(1)$20. For $U(1)$21, one finds $U(1)$22, $U(1)$23, and $U(1)$24. For “hyper-charged” states $U(1)$25, the upper frequency does not reach $U(1)$26, and the binding energy is negative throughout the narrow existence window (Hao et al., 3 Apr 2026).

An important feature of the higher-derivative solutions is that all standard energy conditions were verified numerically across the explored parameter space: $U(1)$27 Even though the tangential pressure $U(1)$28 can be slightly negative in a narrow band, the inequalities remain satisfied (Hao et al., 3 Apr 2026). This distinguishes the construction from regular black-hole mimickers built from exotic matter.

6. Related five-dimensional formulations and research directions

The same five-dimensional Proca-Maxwell nomenclature appears in several adjacent settings. In global AdS$U(1)$29, Proca perturbations can be decomposed into scalar-type master fields $U(1)$30 and a vector-type master field $U(1)$31, each obeying

$U(1)$32

with exact normal-mode frequencies

$U(1)$33

where $U(1)$34, $U(1)$35, and $U(1)$36 (Lopes et al., 2024). In the Maxwell limit $U(1)$37, the scalar-type mode in $U(1)$38 requires Dirichlet-Neumann rather than Dirichlet boundary conditions, yielding

$U(1)$39

while the vector-type Maxwell spectrum is

$U(1)$40

(Lopes et al., 2024).

In five-dimensional Schwarzschild-AdS backgrounds, Proca and Maxwell quasinormal modes obey Schrödinger-type radial equations under ingoing-horizon and Dirichlet-type AdS boundary conditions (Lopes et al., 18 May 2026). A notable $U(1)$41 result is the existence of purely imaginary low-frequency scalar-type Maxwell modes in large black holes,

$U(1)$42

with numerical values $U(1)$43 for $U(1)$44 and $U(1)$45 for $U(1)$46, in agreement with the analytic formula (Lopes et al., 18 May 2026). No unstable modes with $U(1)$47 were found under the stated boundary conditions.

Flat-space and compactified approaches emphasize different structural aspects. Five-dimensional electromagnetism can be formulated pre-metrically through $U(1)$48 and $U(1)$49, or equivalently $U(1)$50, with the five-dimensional null relation $U(1)$51 unifying massless and massive branches (Delphenich, 2018). Compactification of five-dimensional Proca theory on $U(1)$52 yields a four-dimensional theory with one massive vector zero mode, a tower of massive vectors of masses $U(1)$53, and a tower of massive scalars $U(1)$54, while all constraints remain second class and no gauge symmetry survives for $U(1)$55 (Escalante et al., 2014).

Dual descriptions are also available. In $U(1)$56, Maxwell theory is dual to a massless 2-form with action

$U(1)$57

while massive Proca is dual to a massive 3-form and also admits a second-order symmetric-tensor dual description in terms of $U(1)$58, still propagating the correct $U(1)$59 massive spin-1 degrees of freedom (Dalmazi et al., 2011).

A more recent generalized Proca development concerns exact rotating black holes in five dimensions. In that setting, the vector is aligned with a Kerr-Schild null congruence, the on-shell condition $U(1)$60 collapses the functional couplings to constants, and the full nonlinear system reduces to three master equations (Hassaine et al., 18 May 2026). However, that construction explicitly omits the standard Maxwell term and shows that adding a Maxwell sector is incompatible with the aligned Kerr-Schild reduction in the rotating case, except in trivial or stealth-like branches (Hassaine et al., 18 May 2026). This marks an important distinction from the static higher-derivative Proca-Maxwell system, where the Maxwell field is essential and the electric sector drives the charge-induced unfreezing (Hao et al., 3 Apr 2026).

Across these five-dimensional settings, a consistent pattern emerges. The Proca-Maxwell system is not a single model but a family of closely related constructions organized by geometry, boundary conditions, and the role of gravitational couplings. In asymptotically flat higher-derivative gravity, the infinite tower of local curvature corrections yields globally regular, horizonless black-hole mimickers and a charge-sensitive frozen-to-unfrozen transition (Hao et al., 3 Apr 2026). In AdS$U(1)$61, the same spin-1 sectors furnish exact normal modes and distinctive quasinormal spectra (Lopes et al., 2024, Lopes et al., 18 May 2026). In compactified and dual formulations, the theory clarifies how five-dimensional vector dynamics reorganize into lower-dimensional massive and massless sectors (Delphenich, 2018, Escalante et al., 2014, Dalmazi et al., 2011). These lines of work jointly define the modern research landscape of the five-dimensional Proca-Maxwell system.