Einstein-Maxwell-Neutral Scalar Field System

• The Einstein–Maxwell–neutral scalar field system is defined by coupling a Lorentzian metric, a Maxwell field, and an electrically neutral scalar, outlining both minimal and non-minimal interactions.
• It employs analytical methods in static, plane-wave, and higher-dimensional settings to generate exact solutions including scalarized black holes and integrable sigma-model structures.
• The framework reveals practical implications for cosmic censorship, no-hair theorems, and the evolution of charged compact objects through diverse coupling scenarios.

The Einstein–Maxwell–neutral scalar field system is the class of gravitating field theories in which a Lorentzian metric $g_{\mu\nu}$, a Maxwell field $F_{\mu\nu}$, and an electrically neutral scalar field are coupled either minimally or through non-minimal functions of the scalar and curvature. In the minimal massless case, a standard $n$-dimensional action is

$$S[g_{\mu\nu}, A_\mu, \phi] = \int d^n x \sqrt{-g}\left( \frac{1}{2\kappa_n} R - \frac{1}{4} F_{\mu\nu} F^{\mu\nu} - \frac{1}{2} (\nabla \phi)^2 \right),$$

with field equations $$G_{\mu\nu}=\kappa_n(T^{(\mathrm{em})}_{\mu\nu}+T^{(\phi)}_{\mu\nu})$$, $\nabla_\nu F^{\mu\nu}=0$, and $\Box\phi=0$ (Maeda et al., 2019). In four dimensions, the same minimal system appears with a cosmological constant in the $\alpha=0$ sector of Einstein–Maxwell–Scalar–Gauss–Bonnet theory, where the action reduces to Einstein–Hilbert–Maxwell theory with a minimally coupled neutral massless scalar field (Sang et al., 2022). Much of the modern literature concerns extensions in which the scalar is still neutral in the gauge sense—there is no $D_\mu=\nabla_\mu-iqA_\mu$—but couples non-minimally either to the Maxwell invariant $F_{\mu\nu}F^{\mu\nu}$ or to curvature invariants; these extensions support exact plane-wave models, higher-dimensional static families, scalarized black holes, reflecting-star configurations, strong-cosmic-censorship analyses, and global nonlinear evolution results (Fernandes, 2020).

1. Core formulations and model classes

At the level of field content, “neutral” means that the scalar does not carry gauge charge and does not couple minimally to the Maxwell potential. In the minimally coupled system, the scalar stress tensor is

$F_{\mu\nu}$0

while the electromagnetic stress tensor is

$F_{\mu\nu}$1

This structure underlies the higher-dimensional static reductions and solution-generating transformations developed for the minimally coupled massless system (Maeda et al., 2019).

A widely used non-minimal extension is the four-dimensional Einstein–Maxwell–scalar action

$F_{\mu\nu}$2

with scalar-dependent gauge kinetic function $F_{\mu\nu}$3 and, in some models, scalar potential $F_{\mu\nu}$4 (Fernandes, 2020). Another recurrent variant uses a coupling $F_{\mu\nu}$5, with $F_{\mu\nu}$6, for a neutral scalar around charged horizonless stars; in the linearized regime one choice is $F_{\mu\nu}$7, which induces an effective mass

$F_{\mu\nu}$8

in a Reissner–Nordström exterior (Peng, 2019). A further extension couples the scalar to curvature through the Gauss–Bonnet invariant,

$F_{\mu\nu}$9

which reduces to the standard Einstein–Maxwell–neutral scalar field system when $n$0 (Sang et al., 2022).

A distinct minimal model allows the scalar kinetic term to have either sign,

$n$1

with $n$2 for a canonical scalar and $n$3 for a phantom scalar (Fabris et al., 2022).

Model class	Characteristic coupling	Representative phenomenon
Minimal EMS	$n$4	static higher-dimensional classifications
Maxwell-coupled EMS	$n$5	scalarization, hairy black holes, reflecting stars
Curvature-coupled EMS	$n$6	SCC modified by $n$7 regularity

These formulations already indicate the main bifurcation in the subject. Minimal systems are strongly constrained by no-hair and no-soliton results, whereas non-minimal couplings alter both the effective scalar dynamics and the admissible regularity class of the coupled equations.

2. Exact sectors, integrable backgrounds, and solution generation

One explicit realization of the Einstein–Maxwell–neutral scalar field system appears in a plane-wave setting where the Einstein–Maxwell sector is solved exactly and the neutral scalar is treated as a test field. In light-cone coordinates $n$8, $n$9, an electromagnetic plane wave with transverse potential $$S[g_{\mu\nu}, A_\mu, \phi] = \int d^n x \sqrt{-g}\left( \frac{1}{2\kappa_n} R - \frac{1}{4} F_{\mu\nu} F^{\mu\nu} - \frac{1}{2} (\nabla \phi)^2 \right),$$0 sources a Brinkmann pp-wave metric

$$S[g_{\mu\nu}, A_\mu, \phi] = \int d^n x \sqrt{-g}\left( \frac{1}{2\kappa_n} R - \frac{1}{4} F_{\mu\nu} F^{\mu\nu} - \frac{1}{2} (\nabla \phi)^2 \right),$$1

with

$$S[g_{\mu\nu}, A_\mu, \phi] = \int d^n x \sqrt{-g}\left( \frac{1}{2\kappa_n} R - \frac{1}{4} F_{\mu\nu} F^{\mu\nu} - \frac{1}{2} (\nabla \phi)^2 \right),$$2

For constant $$S[g_{\mu\nu}, A_\mu, \phi] = \int d^n x \sqrt{-g}\left( \frac{1}{2\kappa_n} R - \frac{1}{4} F_{\mu\nu} F^{\mu\nu} - \frac{1}{2} (\nabla \phi)^2 \right),$$3, neutral geodesics satisfy a transverse oscillator equation,

$$S[g_{\mu\nu}, A_\mu, \phi] = \int d^n x \sqrt{-g}\left( \frac{1}{2\kappa_n} R - \frac{1}{4} F_{\mu\nu} F^{\mu\nu} - \frac{1}{2} (\nabla \phi)^2 \right),$$4

and a neutral massive scalar satisfies a Klein–Gordon equation that reduces to a two-dimensional harmonic-oscillator problem with discrete transverse spectrum

$$S[g_{\mu\nu}, A_\mu, \phi] = \int d^n x \sqrt{-g}\left( \frac{1}{2\kappa_n} R - \frac{1}{4} F_{\mu\nu} F^{\mu\nu} - \frac{1}{2} (\nabla \phi)^2 \right),$$5

The same background gives analogous mode structure for Majorana spinors and abelian vector fields (Holten, 2022).

In the static sector, the minimally coupled higher-dimensional theory admits a reduced nonlinear sigma-model structure. For

$$S[g_{\mu\nu}, A_\mu, \phi] = \int d^n x \sqrt{-g}\left( \frac{1}{2\kappa_n} R - \frac{1}{4} F_{\mu\nu} F^{\mu\nu} - \frac{1}{2} (\nabla \phi)^2 \right),$$6

the reduced action can be written in terms of three scalars $$S[g_{\mu\nu}, A_\mu, \phi] = \int d^n x \sqrt{-g}\left( \frac{1}{2\kappa_n} R - \frac{1}{4} F_{\mu\nu} F^{\mu\nu} - \frac{1}{2} (\nabla \phi)^2 \right),$$7, and this structure yields two solution-generating maps. The higher-dimensional Buchdahl transformation takes a static vacuum seed and adds a massless neutral scalar, while the higher-dimensional JRW transformation maps a neutral scalar solution to a charged one with the same scalar field (Maeda et al., 2019). Applied to Schwarzschild–Tangherlini-type seeds, these maps generate generalized JNW solutions and their charged relatives; applied to multi-center and cylindrically symmetric seeds, they generate ghost-scalar multi-center configurations and charged scalar-hairy cylindrical solutions (Maeda et al., 2019).

A complementary exact classification treats all static, purely electric solutions with Einstein base manifold in arbitrary dimension. For $$S[g_{\mu\nu}, A_\mu, \phi] = \int d^n x \sqrt{-g}\left( \frac{1}{2\kappa_n} R - \frac{1}{4} F_{\mu\nu} F^{\mu\nu} - \frac{1}{2} (\nabla \phi)^2 \right),$$8, the warped metric

$$S[g_{\mu\nu}, A_\mu, \phi] = \int d^n x \sqrt{-g}\left( \frac{1}{2\kappa_n} R - \frac{1}{4} F_{\mu\nu} F^{\mu\nu} - \frac{1}{2} (\nabla \phi)^2 \right),$$9

reduces the field equations to

$$G_{\mu\nu}=\kappa_n(T^{(\mathrm{em})}_{\mu\nu}+T^{(\phi)}_{\mu\nu})$$0

together with a master equation for $$G_{\mu\nu}=\kappa_n(T^{(\mathrm{em})}_{\mu\nu}+T^{(\phi)}_{\mu\nu})$$1. The resulting classification with a non-constant real scalar field consists of nine solutions for $$G_{\mu\nu}=\kappa_n(T^{(\mathrm{em})}_{\mu\nu}+T^{(\phi)}_{\mu\nu})$$2 and three solutions for $$G_{\mu\nu}=\kappa_n(T^{(\mathrm{em})}_{\mu\nu}+T^{(\phi)}_{\mu\nu})$$3, all expressible in elementary functions. Two features are singled out for $$G_{\mu\nu}=\kappa_n(T^{(\mathrm{em})}_{\mu\nu}+T^{(\phi)}_{\mu\nu})$$4: asymptotically flat solution is not unique when both Maxwell and scalar fields are nontrivial, and some branches asymptotically approach Bertotti–Robinson spacetime. In accordance with the no-hair theorem, none of the solutions are endowed of a Killing horizon (Maeda et al., 2016).

Taken together, these results show that the system has two analytically tractable limits. One is the plane-wave regime, where neutral matter propagates on the “gravity of light” background as an exactly solvable oscillator problem. The other is the static warped-product regime, where the field equations collapse to integrable ordinary differential systems and hidden target-space symmetries organize solution generation.

3. Hair formation, compact objects, and rotating black holes

The most developed phenomenology concerns non-minimally coupled Einstein–Maxwell–scalar models with asymptotically flat black holes. In the “scalarized” class $$G_{\mu\nu}=\kappa_n(T^{(\mathrm{em})}_{\mu\nu}+T^{(\phi)}_{\mu\nu})$$5, the Reissner–Nordström solution is admitted because $$G_{\mu\nu}=\kappa_n(T^{(\mathrm{em})}_{\mu\nu}+T^{(\phi)}_{\mu\nu})$$6, and scalarization is controlled by the effective mass

$$G_{\mu\nu}=\kappa_n(T^{(\mathrm{em})}_{\mu\nu}+T^{(\phi)}_{\mu\nu})$$7

for linearized perturbations around RN. Adding a scalar mass $$G_{\mu\nu}=\kappa_n(T^{(\mathrm{em})}_{\mu\nu}+T^{(\phi)}_{\mu\nu})$$8 shifts the existence line for scalarized solutions to larger $$G_{\mu\nu}=\kappa_n(T^{(\mathrm{em})}_{\mu\nu}+T^{(\phi)}_{\mu\nu})$$9 and shrinks the domain of existence, while quartic self-interaction $\nabla_\nu F^{\mu\nu}=0$0 does not affect the onset of scalarization but enlarges the nonlinear branch. The paper finds the scalarized solutions thermodynamically preferred over RN wherever both exist, and all studied scalarized and dilatonic black holes are stable against spherical perturbations (Fernandes, 2020).

A separate branch appears for the positive-coupling model

$\nabla_\nu F^{\mu\nu}=0$1

Here the RN background is not driven unstable by a tachyonic mode; instead the scalar-hairy branch is a nonlinear effect. The scalar field grows monotonically with the radial coordinate and asymptotically approaches a finite constant, and the existence domain lies entirely within the RN parameter region. Increasing the charge $\nabla_\nu F^{\mu\nu}=0$2 drives the hairy solutions further from RN, while excessive charging returns the system to the hairless branch. Strengthening the coupling $\nabla_\nu F^{\mu\nu}=0$3 increases the deviation from the trivial state but compresses the existence domain. The computed quasinormal modes show these solutions stable against linearized scalar perturbations (Guo et al., 3 Jun 2025).

Neutral-scalar hair also appears around charged horizonless reflecting stars. In a Reissner–Nordström exterior with reflecting boundary condition $\nabla_\nu F^{\mu\nu}=0$4 and asymptotic decay $\nabla_\nu F^{\mu\nu}=0$5, the linearized radial equation

$\nabla_\nu F^{\mu\nu}=0$6

implies a no-hair theorem for sufficiently small $\nabla_\nu F^{\mu\nu}=0$7, whereas for larger coupling the paper finds regular scalar hairy configurations supported by the reflecting star (Peng, 2019).

Slow rotation introduces another layer of structure. In a two-parameter Einstein–Maxwell–scalar theory with gauge kinetic function $\nabla_\nu F^{\mu\nu}=0$8, the static black hole interpolates between Reissner–Nordström and Einstein–Maxwell–dilaton behavior, and its slowly rotating extension modifies gyromagnetic ratio, frame dragging, circular-orbit periods, ISCO radius, and radiative efficiency. The gyromagnetic ratio is always less than $\nabla_\nu F^{\mu\nu}=0$9 in the Einstein–Maxwell–dilaton limit $\Box\phi=0$0, but in the broader theory it can reach the Kerr–Newman value $\Box\phi=0$1 by increasing $\Box\phi=0$2 and $\Box\phi=0$3 simultaneously (Qiu, 2021).

The common theme is that neutral-scalar hair in Einstein–Maxwell systems is strongly model-dependent. Minimal neutral scalars are highly constrained, but couplings of the form $\Box\phi=0$4 can generate scalarized or disconnected hairy branches, alter thermodynamic ordering, and change the ringdown and orbital properties of the resulting compact objects.

4. Black-hole interiors, cosmic censorship, and singularity structure

For the spherically symmetric Einstein–Maxwell–scalar field system with a neutral scalar and nonzero charge, rigorous nonlinear evolution results connect exterior decay to interior blow-up. The system can be reduced in double-null coordinates to equations for the area radius $\Box\phi=0$5, renormalized Hawking mass $\Box\phi=0$6, and scalar $\Box\phi=0$7, with

$\Box\phi=0$8

For compactly supported future-admissible data, one obtains decay for higher-order derivatives of the scalar in the exterior, sharp late-time tails of Price-law type, and, when combined with earlier lower bounds, generic mass inflation at the Cauchy horizon. The same framework implies $\Box\phi=0$9-future-inextendibility for a generic class of data, and also supports the global construction of two-ended black holes containing both null and spacelike singularities (Gautam, 2024).

In de Sitter settings, the neutral scalar is frequently used as a probe of strong cosmic censorship. In Einstein–Maxwell–Scalar–Gauss–Bonnet theory, the scalar equation

$\alpha=0$0

and the higher-derivative structure of $\alpha=0$1 tighten the weak-extendibility criterion from the familiar $\alpha=0$2 condition to $\alpha=0$3. Near the Cauchy horizon this yields the threshold

$\alpha=0$4

for extendibility. The quasinormal-mode analysis shows that, for $\alpha=0$5, SCC is restored in the interval $\alpha=0$6, whereas for $\alpha=0$7 and $\alpha=0$8 SCC remains violated (Sang et al., 2022).

A higher-dimensional test-field version of the same question appears in Einstein–Maxwell–Gauss–Bonnet–de Sitter black holes. There the neutral massless scalar obeys $\alpha=0$9 on the fixed EMGB background, and the criterion is again stated in terms of

$D_\mu=\nabla_\mu-iqA_\mu$0

with Christodoulou SCC violated for $D_\mu=\nabla_\mu-iqA_\mu$1 and the $D_\mu=\nabla_\mu-iqA_\mu$2 version violated for $D_\mu=\nabla_\mu-iqA_\mu$3. The behavior is dimension-dependent: in $D_\mu=\nabla_\mu-iqA_\mu$4, the SCC-violation region grows as the Gauss–Bonnet coupling increases and $D_\mu=\nabla_\mu-iqA_\mu$5 SCC can be violated near extremality; in $D_\mu=\nabla_\mu-iqA_\mu$6, the violation region first grows and then shrinks, and the $D_\mu=\nabla_\mu-iqA_\mu$7 version is respected (Gan et al., 2019).

A different interior and global-structure classification arises in the minimal massless system with kinetic sign $D_\mu=\nabla_\mu-iqA_\mu$8. For $D_\mu=\nabla_\mu-iqA_\mu$9, static spherically symmetric solutions yield only naked singularities. For $F_{\mu\nu}F^{\mu\nu}$0, the same Einstein–Maxwell–neutral scalar system admits two families of charged black holes with infinite-area, zero-temperature horizons and three families of traversable wormholes; the distinction is entirely controlled by the sign of the scalar kinetic term (Fabris et al., 2022).

These results make the interior problem central to the subject. Neutral scalars may be test fields or fully coupled matter, but in both roles they diagnose whether charged black-hole interiors are weakly extendible, mass-inflating, or converted into null or spacelike singular boundaries.

5. No-hair theorems, no-soliton theorems, and structural obstructions

The existence of hairy black holes in non-minimal theories is accompanied by strong no-go results for regular horizonless configurations. For the generic Einstein–Maxwell–scalar model

$F_{\mu\nu}F^{\mu\nu}$1

three non-existence theorems are established. First, in the consistent truncation $F_{\mu\nu}F^{\mu\nu}$2, a Derrick-type virial identity excludes stationary and axisymmetric self-gravitating scalar solitons unless the scalar potential is somewhere negative; rotation alone cannot support such solitons. Second, for static spacetimes and constant-sign $F_{\mu\nu}F^{\mu\nu}$3, a Heusler-type argument shows that no self-gravitating electromagnetic-scalar solitons exist. Third, a Lichnerowicz-type argument extends the non-existence result to strictly stationary spacetimes with non-negative scalar potential and dominant energy condition (Herdeiro et al., 2019).

These obstructions align with the higher-dimensional static classification, where the complete set of non-constant real scalar solutions has no Killing horizon and displays curvature singularities whenever $F_{\mu\nu}F^{\mu\nu}$4 (Maeda et al., 2016). They also align with the minimal four-dimensional classification in which canonical scalars yield only naked singularities, and only the phantom sign opens black-hole or wormhole branches (Fabris et al., 2022).

A plausible implication is that regular, asymptotically flat, horizonless objects in Einstein–Maxwell–neutral scalar theories require at least one of the assumptions behind the no-go theorems to fail. The data surveyed here point to three such failure modes: a negative region in the scalar potential, a sign-changing or singular gauge kinetic function $F_{\mu\nu}F^{\mu\nu}$5, or a phantom kinetic term. Each route appears explicitly in the literature, but each also changes the energy-condition content of the model (Herdeiro et al., 2019).

6. Analytical structures and representative mathematical techniques

The subject is held together by a set of recurring analytical reductions. In the static sector, the minimally coupled higher-dimensional system admits a reduced action

$F_{\mu\nu}F^{\mu\nu}$6

which can be rewritten as a nonlinear sigma model with target-space metric

$F_{\mu\nu}F^{\mu\nu}$7

Buchdahl and JRW transformations are then understood as target-space symmetries acting on $F_{\mu\nu}F^{\mu\nu}$8 (Maeda et al., 2019).

In the dynamical spherical problem, large-data exterior decay and mass inflation are derived through a vector-field hierarchy built from double-null geometry, the Hawking mass, and commutator fields $F_{\mu\nu}F^{\mu\nu}$9, $F_{\mu\nu}$00, and $F_{\mu\nu}$01, where $F_{\mu\nu}$02 acts as a scaling vector field near the horizon. The strategy combines energy boundedness, integrated local energy decay, $F_{\mu\nu}$03-weighted estimates, and commutator identities for the reduced wave operator (Gautam, 2024). This suggests that, even in the fully coupled Einstein–Maxwell–neutral scalar system, the decisive analytical objects are not only the metric and the scalar itself, but also renormalized geometric quantities such as $F_{\mu\nu}$04, $F_{\mu\nu}$05, and the adapted null derivatives.

Finally, exact solvability in the plane-wave Einstein–Maxwell background shows in especially transparent form how a neutral field interacts with electromagnetism only through curvature. The minimally coupled scalar, a Majorana spinor component, and the transverse modes of an abelian vector all satisfy the same effective oscillator equation in the constant-flux plane-wave geometry, with discrete transverse spectra set by the electromagnetic energy density through the induced pp-wave curvature (Holten, 2022).

Across these settings, the Einstein–Maxwell–neutral scalar field system is best understood not as a single model but as a structured family of theories. Minimal couplings emphasize classification, singularity formation, and no-hair constraints; non-minimal Maxwell or curvature couplings open scalarized branches and modify cosmic censorship; exact plane-wave and higher-dimensional reductions provide solvable laboratories; and nonlinear PDE results connect exterior decay to interior blow-up. The literature therefore presents the system as a junction of exact relativity, geometric analysis, and compact-object phenomenology rather than as a single canonical equation set.