Schwarzschild–Bonnor–Melvin Geometry

• Schwarzschild–Bonnor–Melvin is a static, axisymmetric solution that immerses a Schwarzschild black hole in a uniform magnetic field via a Harrison transformation.
• The magnetic field distorts the horizon geometry without altering the event horizon radius, area, or thermodynamic quantities, reinforcing key uniqueness properties.
• Extensions include nonlinear electrodynamics, scalar hair, and spinorial rigidity methods, highlighting its role as a benchmark in exploring diverse gravitational models.

Schwarzschild–Bonnor–Melvin denotes the static, axisymmetric spacetime obtained by immersing a Schwarzschild black hole in the Bonnor–Melvin magnetic universe. In the standard Einstein–Maxwell realization, with

$$f(r)=1-\frac{2m}{r},\qquad \Lambda(r,\theta)=1+\frac{B^2}{4}r^2\sin^2\theta,$$

the line element and Maxwell potential are

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$

$$A_\varphi(r,\theta)=\frac{B\,r^2\sin^2\theta}{2\,\Lambda(r,\theta)},\qquad A_t=A_r=A_\theta=0,$$

with $B$ measuring the strength of the asymptotically uniform magnetic field along the $z$-axis. The geometry is therefore not asymptotically flat but approaches the Bonnor–Melvin universe, and in the mathematical relativity literature it appears as the magnetized Schwarzschild solution whose spatial slice is asymptotically Melvin (Astorino, 18 Aug 2025, Masood-ul-Alam et al., 2014).

1. Canonical Einstein–Maxwell form

The standard construction begins from the vacuum Schwarzschild seed written in Schwarzschild-like coordinates $(t,r,\theta,\varphi)$,

$$ds^2_{\rm Schw} = -\Bigl(1-\frac{2m}{r}\Bigr)\,dt^2 + \Bigl(1-\frac{2m}{r}\Bigr)^{-1}dr^2 + r^2\Bigl(d\theta^2+\sin^2\theta\,d\varphi^2\Bigr).$$

In the Ernst formulation one introduces the complex gravitational potential $\mathcal{E}_0=f_0=1-\tfrac{2m}{r}$ and electromagnetic Ernst potential $\Phi_0=0$. A purely magnetic Harrison transformation with real parameter $B$ acts through the function

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$0

and sends

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$1

while leaving the twist $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$2 and the function $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$3 essentially unchanged up to gauge. The resulting line element is precisely the Schwarzschild–Bonnor–Melvin metric written above, and the only nonzero component of the Maxwell four-potential is $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$4 (Astorino, 18 Aug 2025).

This places the solution within the standard Harrison-transform orbit of static axisymmetric Einstein–Maxwell spacetimes. The same source also states that the pure Schwarzschild–Bonnor–Melvin solution is recovered as a limit of the more general Schwarzschild–Bertotti–Robinson–Bonnor–Melvin family by setting

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$5

that is, by turning off the Bertotti–Robinson sector and renaming the Harrison parameter as the physical magnetic field (Astorino, 18 Aug 2025).

2. Horizon data, distortion, and asymptotics

For the standard Schwarzschild–Bonnor–Melvin metric, the Killing horizon of $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$6 is located by $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$7, hence

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$8

The surface gravity remains the Schwarzschild value,

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$9

and the Hawking temperature is correspondingly

$$A_\varphi(r,\theta)=\frac{B\,r^2\sin^2\theta}{2\,\Lambda(r,\theta)},\qquad A_t=A_r=A_\theta=0,$$0

A Komar or ADM mass computed at infinity remains $$A_\varphi(r,\theta)=\frac{B\,r^2\sin^2\theta}{2\,\Lambda(r,\theta)},\qquad A_t=A_r=A_\theta=0,$$1, with the caveat that the normalization of $$A_\varphi(r,\theta)=\frac{B\,r^2\sin^2\theta}{2\,\Lambda(r,\theta)},\qquad A_t=A_r=A_\theta=0,$$2 in a non-flat asymptotic requires care but does not alter the identification $$A_\varphi(r,\theta)=\frac{B\,r^2\sin^2\theta}{2\,\Lambda(r,\theta)},\qquad A_t=A_r=A_\theta=0,$$3 (Astorino, 18 Aug 2025).

The horizon geometry is distorted by the external field. If the two-dimensional horizon surface $$A_\varphi(r,\theta)=\frac{B\,r^2\sin^2\theta}{2\,\Lambda(r,\theta)},\qquad A_t=A_r=A_\theta=0,$$4 is embedded into $$A_\varphi(r,\theta)=\frac{B\,r^2\sin^2\theta}{2\,\Lambda(r,\theta)},\qquad A_t=A_r=A_\theta=0,$$5, the polar radius is larger than the equatorial radius. The equatorial circumference is

$$A_\varphi(r,\theta)=\frac{B\,r^2\sin^2\theta}{2\,\Lambda(r,\theta)},\qquad A_t=A_r=A_\theta=0,$$6

so it is reduced relative to the Schwarzschild value $$A_\varphi(r,\theta)=\frac{B\,r^2\sin^2\theta}{2\,\Lambda(r,\theta)},\qquad A_t=A_r=A_\theta=0,$$7, whereas the polar semicircumference is slightly enlarged. Nevertheless the horizon area remains exactly

$$A_\varphi(r,\theta)=\frac{B\,r^2\sin^2\theta}{2\,\Lambda(r,\theta)},\qquad A_t=A_r=A_\theta=0,$$8

the same as for vacuum Schwarzschild. No conical singularities arise on the axis provided $$A_\varphi(r,\theta)=\frac{B\,r^2\sin^2\theta}{2\,\Lambda(r,\theta)},\qquad A_t=A_r=A_\theta=0,$$9 has its standard period $B$0 and $B$1 everywhere (Astorino, 18 Aug 2025).

The coordinate ranges are the usual exterior ones,

$B$2

Because $B$3 for all $B$4, no new coordinate singularities appear outside $B$5. The only genuine curvature singularity lies at $B$6. As $B$7, $B$8 rather than a constant, so the spacetime is not asymptotically flat but approaches the Bonnor–Melvin universe (Astorino, 18 Aug 2025).

A common misconception is that the external magnetic field should shift the horizon radius or the thermodynamic temperature. In the standard Einstein–Maxwell solution neither effect occurs: the field deforms the horizon intrinsically while leaving $B$9, $z$0, $z$1, and $z$2 unchanged (Astorino, 18 Aug 2025).

3. Asymptotically Melvin initial data and harmonic spinors

On a time-symmetric spatial slice, the relevant geometric setting is an axisymmetric Riemannian $z$3-manifold $z$4 with metric

$z$5

where $z$6 are cylindrical or spherical coordinates at each end, $z$7 is the Killing coordinate, and $z$8 and $z$9 are $(t,r,\theta,\varphi)$0-independent. The Killing field $(t,r,\theta,\varphi)$1 generates closed orbits and its zero set is the symmetry axis. The manifold is asymptotically Melvin with field-strength parameter $(t,r,\theta,\varphi)$2 if on each end $(t,r,\theta,\varphi)$3 one can write

$(t,r,\theta,\varphi)$4

with

$(t,r,\theta,\varphi)$5

and remainders $(t,r,\theta,\varphi)$6, $(t,r,\theta,\varphi)$7, $(t,r,\theta,\varphi)$8. In particular, as $(t,r,\theta,\varphi)$9, $$ds^2_{\rm Schw} = -\Bigl(1-\frac{2m}{r}\Bigr)\,dt^2 + \Bigl(1-\frac{2m}{r}\Bigr)^{-1}dr^2 + r^2\Bigl(d\theta^2+\sin^2\theta\,d\varphi^2\Bigr).$$0, $$ds^2_{\rm Schw} = -\Bigl(1-\frac{2m}{r}\Bigr)\,dt^2 + \Bigl(1-\frac{2m}{r}\Bigr)^{-1}dr^2 + r^2\Bigl(d\theta^2+\sin^2\theta\,d\varphi^2\Bigr).$$1, and $$ds^2_{\rm Schw} = -\Bigl(1-\frac{2m}{r}\Bigr)\,dt^2 + \Bigl(1-\frac{2m}{r}\Bigr)^{-1}dr^2 + r^2\Bigl(d\theta^2+\sin^2\theta\,d\varphi^2\Bigr).$$2 (Masood-ul-Alam et al., 2014).

For $$ds^2_{\rm Schw} = -\Bigl(1-\frac{2m}{r}\Bigr)\,dt^2 + \Bigl(1-\frac{2m}{r}\Bigr)^{-1}dr^2 + r^2\Bigl(d\theta^2+\sin^2\theta\,d\varphi^2\Bigr).$$3-independent spinors $$ds^2_{\rm Schw} = -\Bigl(1-\frac{2m}{r}\Bigr)\,dt^2 + \Bigl(1-\frac{2m}{r}\Bigr)^{-1}dr^2 + r^2\Bigl(d\theta^2+\sin^2\theta\,d\varphi^2\Bigr).$$4, one defines weighted Sobolev norms relative to $$ds^2_{\rm Schw} = -\Bigl(1-\frac{2m}{r}\Bigr)\,dt^2 + \Bigl(1-\frac{2m}{r}\Bigr)^{-1}dr^2 + r^2\Bigl(d\theta^2+\sin^2\theta\,d\varphi^2\Bigr).$$5, and the closed subspace

$$ds^2_{\rm Schw} = -\Bigl(1-\frac{2m}{r}\Bigr)\,dt^2 + \Bigl(1-\frac{2m}{r}\Bigr)^{-1}dr^2 + r^2\Bigl(d\theta^2+\sin^2\theta\,d\varphi^2\Bigr).$$6

Under the assumptions $$ds^2_{\rm Schw} = -\Bigl(1-\frac{2m}{r}\Bigr)\,dt^2 + \Bigl(1-\frac{2m}{r}\Bigr)^{-1}dr^2 + r^2\Bigl(d\theta^2+\sin^2\theta\,d\varphi^2\Bigr).$$7, asymptotically Melvin with parameter $$ds^2_{\rm Schw} = -\Bigl(1-\frac{2m}{r}\Bigr)\,dt^2 + \Bigl(1-\frac{2m}{r}\Bigr)^{-1}dr^2 + r^2\Bigl(d\theta^2+\sin^2\theta\,d\varphi^2\Bigr).$$8, and finitely many ends, the Dirac operator

$$ds^2_{\rm Schw} = -\Bigl(1-\frac{2m}{r}\Bigr)\,dt^2 + \Bigl(1-\frac{2m}{r}\Bigr)^{-1}dr^2 + r^2\Bigl(d\theta^2+\sin^2\theta\,d\varphi^2\Bigr).$$9

is an isomorphism onto its range. In particular there is a unique, up to normalization, nontrivial spinor $\mathcal{E}_0=f_0=1-\tfrac{2m}{r}$0 solving

$\mathcal{E}_0=f_0=1-\tfrac{2m}{r}$1

This is the existence theorem for the harmonic spinor used in the uniqueness analysis of the magnetized Schwarzschild solution (Masood-ul-Alam et al., 2014).

The proof proceeds by introducing an auxiliary asymptotically flat metric

$\mathcal{E}_0=f_0=1-\tfrac{2m}{r}$2

using Bartnik’s semi-Fredholm result for $\mathcal{E}_0=f_0=1-\tfrac{2m}{r}$3, and comparing $\mathcal{E}_0=f_0=1-\tfrac{2m}{r}$4 to $\mathcal{E}_0=f_0=1-\tfrac{2m}{r}$5. With appropriately chosen orthonormal frames,

$\mathcal{E}_0=f_0=1-\tfrac{2m}{r}$6

so $\mathcal{E}_0=f_0=1-\tfrac{2m}{r}$7 is a compact perturbation of a semi-Fredholm operator. Formal self-adjointness, nonnegativity of $\mathcal{E}_0=f_0=1-\tfrac{2m}{r}$8, and the Weitzenböck–Lichnerowicz identity then imply trivial adjoint kernel and hence index zero (Masood-ul-Alam et al., 2014).

In a local orthonormal frame $\mathcal{E}_0=f_0=1-\tfrac{2m}{r}$9, the Dirac operator is

$\Phi_0=0$0

and the Weitzenböck–Lichnerowicz identity is

$\Phi_0=0$1

equivalently

$\Phi_0=0$2

Once $\Phi_0=0$3 satisfies $\Phi_0=0$4 and tends to a constant spinor at infinity, integration over large coordinate balls $\Phi_0=0$5 yields

$\Phi_0=0$6

A direct asymptotic expansion shows

$\Phi_0=0$7

and since the $\Phi_0=0$8-area form on $\Phi_0=0$9 is $B$0,

$B$1

Thus the Melvin parameter is extracted from the spinor energy by

$B$2

This identifies the asymptotic magnetic parameter through a positive-mass-type integral on the initial slice (Masood-ul-Alam et al., 2014).

4. Rigidity and uniqueness of the magnetized Schwarzschild solution

The harmonic-spinor construction is applied to the time-symmetric slice $B$3 of a static Einstein–Maxwell spacetime with two ends and a totally geodesic horizon surface. Two auxiliary metrics are introduced on $B$4,

$B$5

where

$B$6

Here $B$7 is the mass parameter and $B$8 the background Melvin parameter through $B$9 (Masood-ul-Alam et al., 2014).

By conformal-type spinor transformation laws, each $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$00 admits a harmonic spinor pulled back from $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$01, satisfying

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$02

together with zero-mass asymptotics for $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$03 and regularity at the fixed-point set for $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$04. The integrated Weitzenböck identities on $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$05 force both metrics to be flat. Flatness of $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$06 then implies that $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$07 must coincide with the known Schwarzschild–Bonnor–Melvin family. Accordingly, the only static, axisymmetric, asymptotically Melvin solution with a nondegenerate horizon is the magnetized Schwarzschild solution found by Bonnor and Melvin (Masood-ul-Alam et al., 2014).

This rigidity statement is significant because it is formulated in a non-asymptotically-flat setting. The role usually played by asymptotically Euclidean spinorial arguments is replaced by a Melvin-end analysis in weighted Sobolev spaces, together with the extraction of the Melvin parameter from the asymptotics of the harmonic spinor. A plausible implication is that the Schwarzschild–Bonnor–Melvin geometry is distinguished not merely by explicit solution generation but also by a boundary-value characterization intrinsic to asymptotically Melvin initial data.

5. Nonlinear-electrodynamic extension: the ModMax case

In Einstein–ModMax theory, the Schwarzschild–Melvin–Bonnor configuration is generalized by allowing both asymptotic electric and magnetic fields. Denoting the ModMax coupling by $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$08, the black-hole mass by $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$09, and the asymptotic fields by $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$10 and $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$11, one defines

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$12

and writes the metric and gauge potential as

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$13

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$14

In the Maxwell limit $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$15, one recovers the standard Ernst–Schwarzschild–Melvin solution (Barrientos et al., 2024).

The construction is obtained from the ModMax $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$16-metric of two oppositely charged accelerating black holes through a limiting procedure: one moves to a near-Rindler-horizon scaling, sends the acceleration configuration to infinity, and is left with the self-gravitating homogeneous field identified as the ModMax Melvin–Bonnor universe. A subsequent embedding of a Schwarzschild black hole yields the full metric and potential above. The paper attributes the starting ModMax $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$17-metric to Barrientos et al. 2022 and states that for all Melvin–Bonnor–ModMax configurations one has $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$18, so the full field equations reduce to the Maxwell form (Barrientos et al., 2024).

A distinctive feature is a Kerr–Schild representation. With seed fields

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$19

and null, geodesic, shear-free one-form

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$20

the full solution is

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$21

The field equations are satisfied provided $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$22 is null and geodesic with respect to $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$23 (Barrientos et al., 2024).

The global structure remains close to the Maxwellian case. The horizon is at $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$24, no additional Melvin-type horizons arise, the Penrose diagrams coincide with Schwarzschild at fixed $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$25, and no ergoregion or superradiant region appears since the solution is static. The Kretschmann scalar remains finite on and outside the event horizon away from the physical singularity at $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$26, and the ModMax energy–momentum tensor respects the weak, strong, and dominant energy conditions for $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$27. The presence of the asymptotic electric field $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$28 and the screening factor $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$29 modifies the geometry and field lines, while in the limit $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$30 the electric sector is completely suppressed (Barrientos et al., 2024).

6. Scalar, dilatonic, and baryonic extensions

One line of generalization embeds Schwarzschild-like compact objects with scalar hair into Melvin-type fields in Einstein–Maxwell–dilaton theories. Starting from the Just–Fisher–Janis–Robinson–Winicour seed with parameters $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$31 and $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$32,

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$33

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$34

one obtains, by the Harrison–Demiański–Dowker procedure, an exact axisymmetric magnetized solution with a uniform magnetic field asymptotically along the $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$35-axis. In the special case $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$36, the dilaton vanishes and one recovers exactly the standard Schwarzschild–Melvin (Bonnor–Melvin) metric of general relativity. When $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$37, the solution reduces to the JFJW seed. The scalar charge is identified from the asymptotic decay

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$38

For $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$39 there is a regular event horizon at $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$40, but for $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$41 the surface $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$42 remains a Killing horizon while the Ricci and Kretschmann scalars diverge there, so the geometry is a naked singularity rather than a black hole. The same work also gives a conformal map to the Entangled Relativity frame with $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$43, and states that in the limit $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$44 the entangled-frame solution reduces exactly to the analytic compact-star solution used in Arruga and Minazzoli 2021 (Minazzoli et al., 19 Feb 2025).

A second extension introduces scalar dressing and, after a dictionary to the gauged Skyrme–Maxwell–Einstein model, a discrete Baryonic charge. In coordinates $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$45 with

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$46

the magnetized scalar-dressed Schwarzschild–Melvin metric is

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$47

with

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$48

The Harrison-generated Maxwell field is purely magnetic, with

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$49

and one checks directly that $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$50 in this background. The horizon remains at $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$51, the coordinates cover $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$52, and as $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$53 the metric approaches the Melvin universe rather than flat space (Barrientos et al., 27 Jan 2026).

In the Einstein–Scalar–Maxwell seed one introduces

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$54

which satisfies $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$55. After mapping to the Skyrme sector, the entire Baryonic charge is carried by the Callan–Witten term and, after subtracting the pure Melvin background, one finds

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$56

Since $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$57 must be an integer, this relation implicitly quantizes the ratio $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$58 in terms of $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$59, or equivalently renders the black-hole mass $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$60 a function of the discrete Baryon number $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$61 and the magnetic field $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$62. In the large-mass regime,

$$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$63

whereas for moderate or small $$ds^2=\Lambda^2(r,\theta)\Bigl[-f(r)\,dt^2+f(r)^{-1}dr^2+r^2\,d\theta^2\Bigr]+\Lambda^{-2}(r,\theta)\,r^2\sin^2\theta\,d\varphi^2,$$64 sizable deviations from linearity appear. The paper describes this as the first closed-form analytic model of a Schwarzschild black hole carrying a discrete topological Baryonic charge while immersed in a fully backreacting external Melvin magnetic field (Barrientos et al., 27 Jan 2026).

Taken together, these developments show that the Schwarzschild–Bonnor–Melvin geometry serves as a core axisymmetric background across several extensions: pure Einstein–Maxwell magnetization, Melvin-end spinorial rigidity, nonlinear electrodynamics with electric screening, scalar and dilatonic hair, conformal reformulations, and topological charge sectors. This suggests that its importance lies not only in its explicit metric form but also in its role as a reference solution for uniqueness theorems, solution-generating techniques, and matter-coupled deformations.