---
title: Regular Magnetically Charged BHs (NRCBHs)
url: https://www.emergentmind.com/topics/regular-magnetically-charged-black-holes-nrcbhs
type: topic
---

# Regular Magnetically Charged BHs (NRCBHs)

Searching arXiv for recent and foundational papers on regular magnetically charged black holes, nonlinear electrodynamics, and related perturbative/observational studies.
Regular magnetically charged black holes (NRCBHs) are magnetically charged black-hole geometries sourced most often by nonlinear electrodynamics (NED) and constructed so that the Reissner–Nordström central singularity is replaced either by a regular core or, in black-bounce realizations, by a minimal-area bounce. In the standard spherically symmetric sector the metric is typically written as \(ds^2=-f(r)dt^2+f(r)^{-1}dr^2+r^2d\Omega^2\), while the magnetic field is a monopole with \(F_{\theta\phi}=Q_m\sin\theta\); depending on convention, the electromagnetic invariant is written either as \(F=2Q_m^2/r^4\) or \(\mathcal{F}=q^2/(2r^4)\). The literature includes minimally coupled GR constructions, Fan–Wang families, Bronnikov- and Kruglov-type models, rational and power-Maxwell theories, and extensions in 4D Einstein–Gauss–Bonnet (EGB), Einstein cubic gravity (ECG), and \(f(R,T)\) gravity [1409.5988][2210.00295][1509.05580][2005.00360][2312.16614][2506.13676].

## 1. Core definitions and source theories

A common starting point is Einstein gravity minimally coupled to NED,
\[
S=\frac{1}{16\pi}\int d^4x\sqrt{-g}\,[R-L(F)],
\]
or closely related sign-convention variants. In this framework, Bronnikov-type magnetic solutions use
\[
\mathcal{L}(F)=F\cosh^{-2}\!\left(a\left|\frac{F}{2}\right|^{1/4}\right),
\]
while Fan and Wang introduced a one-parameter family
\[
L(F)=\frac{4\mu}{\alpha}(\alpha F)^{(\nu+3)/4}\left[1+(\alpha F)^{\nu/4}\right]^{-(\mu+\nu)/\nu},
\]
with \(\mu=3\) and \(\nu=1,2,3\) corresponding to Maxwellian, Bardeen-like, and Hayward-like subclasses, respectively. In these models the magnetic potential is \(A=Q_m\cos\theta\,d\phi\), the field invariant is \(F(r)=2Q_m^2/r^4\), and the magnetic parameter is related by \(Q_m=q^2/\sqrt{2\alpha}\) [1409.5988][2210.00295].

Other minimally coupled NED models generate regular magnetic geometries with different ultraviolet behavior. Examples include Kruglov’s exponential NED,
\[
\mathcal{L}_{NED}=-\mathcal{F}\exp(-\beta\mathcal{F}),
\]
the rational model
\[
{\cal L}=-\frac{\cal F}{1+2\beta{\cal F}},
\]
the cosh-type model
\[
\mathcal{L}(\mathcal{F})=-\frac{\mathcal{F}}{\cosh u},
\qquad
u=(|\beta\mathcal{F}|)^{1/4},
\]
and the deformed birefringent model
\[
\mathcal{L}=\frac{\mathcal{F}}{1+2\beta\mathcal{F}-\frac{\gamma}{2}\mathcal{G}^2}.
\]
Across these papers, the weak-field Maxwell limit is built in, while magnetic charge rather than electric charge is the preferred route to regularity in the minimally coupled setting [2005.00360][2108.07175][1803.02191][1509.05580].

The same theme extends beyond GR. Exact or perturbative NRCBHs have been constructed in 4D EGB gravity, in ECG coupled to NLE, and in \(f(R,T)=R+2\beta T\) gravity coupled to NLED. These modified-gravity versions preserve the magnetic monopole sector but alter the branch structure, the field equations, and the thermodynamic and perturbative spectra [2005.00360][2305.18254][2312.16614].

## 2. Regularity, core geometries, and energy conditions

The defining issue is whether the center is genuinely regular. Different NRCBH families realize this in different ways. The Bronnikov magnetic black hole has
\[
f(r)=1-\frac{2M}{r}\left[1-\tanh\!\left(\frac{q^2}{2Mr}\right)\right],
\]
for which \(f(r)\to 1\) as \(r\to 0\), the NED Lagrangian and energy density decay exponentially, and the curvature invariants satisfy
\[
R\to 0,\qquad
R_{\mu\nu}R^{\mu\nu}\to 0,\qquad
R_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}\to 0,
\]
yielding a Minkowski-like core. A related cosh-NED model likewise has \(f(0)=1\), exponentially vanishing energy density, and curvature invariants that vanish at the center [1409.5988][1803.02191].

By contrast, the Fan–Wang families with \(\mu=3\) approach
\[
f(r)\approx 1-\frac{2M}{q^3}r^2,
\]
so the core is de Sitter with
\[
\Lambda_{\rm eff}=\frac{6M}{q^3},
\]
and the curvature scalars remain finite:
\[
R\to 4\Lambda_{\rm eff},\qquad
R_{\mu\nu}R^{\mu\nu}\to 4\Lambda_{\rm eff}^2,\qquad
K\to \frac{8}{3}\Lambda_{\rm eff}^2.
\]
The modified Hayward geometry sourced by rational NLED also has
\[
f(r)=1-\frac{r^2}{l^2}+\cdots,
\]
so its core is de Sitter with \(\Lambda_{\rm eff}=3/l^2\), and the Ricci and Kretschmann invariants are bounded; the paper explicitly states that the limiting curvature conjecture takes place. The \(f(R,T)\) constructions likewise admit a regular center, with \(A(r)\approx 1-\text{const}\times r^2\) for \(a=3\), and finite \(R\) and \(K\) for \(a\ge 3\) [2210.00295][2108.07175][2312.16614].

Energy conditions are model dependent. The deformed birefringent NED model of Ma is built to satisfy the weak energy condition, and the paper verifies \(\rho\ge 0\), \(\rho+p_r=0\), and \(\rho+p_t\ge 0\) everywhere for the magnetic solution. In \(f(R,T)\)+NLED, the null, weak, and strong energy conditions can all be arranged to hold outside the outer event horizon by appropriate parameter choices. By contrast, some papers do not analyze the weak, null, or dominant energy conditions at all; this is explicit in the 4D EGB exponential-NED work and in the rational-NED modified Hayward model [1509.05580][2312.16614][2005.00360][2108.07175].

The label “regular” is therefore not uniform across the literature. The 4D EGB magnetic solution of Hegde, Karananas, and collaborators is discussed in the context of regular black holes, but the paper does not compute \(R\), \(R_{\mu\nu}R^{\mu\nu}\), or \(R_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}\), and does not provide a near-origin expansion proving regularity at \(r=0\). In that case, the central regularity claim is contextual rather than invariant-based [2005.00360].

## 3. Horizons, extremality, and parameter space

In essentially all NRCBH models, horizons are defined by the zeros of \(f(r)\) or \(A(r)\), and the same qualitative trichotomy recurs: two horizons, one degenerate extremal horizon, or no horizon. The Bronnikov solution has subextremal, extremal, and horizonless sectors; for \(M=1\) the paper gives
\[
q_{\rm ext}\simeq 1.0554,\qquad r_{\rm ext}\simeq 0.8712,
\]
obtained from \(f(r_0)=0=f'(r_0)\). Fan–Wang solutions show analogous behavior, with allowed charge intervals
\[
q/M\in[0,0.2963]\quad (\nu=1),\qquad
q/M\in[0,0.7698]\quad (\nu=2),\qquad
q/M\in[0,1.0583]\quad (\nu=3),
\]
for black-hole horizons when \(\mu=3\) [1409.5988][2210.00295].

The horizon structure is particularly sensitive to modified-gravity couplings. In 4D EGB with exponential NED, the physically viable branch is the negative branch of the square-root solution, and the paper reports explicit critical values for horizon disappearance. For \(M=1\), \(r_+=1\), and \(\beta=0.1\), increasing \(q=\{0.1,0.2,0.3,0.4,0.5\}\) yields
\[
\alpha_c=\{0.990007,0.960109,0.910537,0.841623,0.753734\},
\]
and for \(\alpha>\alpha_c\) one has \(f(r_+)>0\) and thus no horizon. In ECG, sufficiently large cubic coupling can also eliminate horizons and produce naked singularities; the paper gives the neutral examples \(G^2\lambda=16(GM)^4\) as nakedly singular and \(G^2\lambda=8(GM)^4\) as still black-hole-like [2005.00360][2305.18254].

Some recent observationally oriented NRCBH metrics exhibit larger extremal magnetic charges than Reissner–Nordström. One 2025 model reports
\[
q_{\rm ext}=2.537862\,M
\]
from \(f(r_e)=0=f'(r_e)\), while also noting that its own table and figure captions contain inconsistent values around \(q\approx 3.79\) for \(M=1\). That paper therefore provides a concrete example in which the existence of an extremal bound is clear, but its precise numerical value is internally unsettled [2507.05145].

## 4. Thermodynamics and phase structure

The thermodynamic analysis is usually built on the standard surface-gravity formula
\[
T_H=\frac{f'(r_h)}{4\pi},
\]
with entropy given either by the area law \(S=A/4\) in GR and minimally coupled NED, or by Wald’s entropy in higher-curvature gravity. In many NRCBH families, \(T_H\) vanishes at extremality, develops a maximum at intermediate radius, and is accompanied by heat-capacity divergences interpreted as second-order phase transitions. This pattern appears in Bronnikov-type, cosh-NED, rational-NED, Fan–Wang-type, and ECG-corrected models [1409.5988][1803.02191][2108.07175][2005.00360][2305.18254].

The detailed phase structure depends strongly on the matter model. In the cosh-NED magnetic solution,
\[
T_H=0
\]
at the extremal point \(x_h\simeq 0.93\), \(T_H\) reaches a maximum at \(x_+\simeq 1.82\), and the fixed-charge heat capacity diverges there, signaling a second-order phase transition; the outer branch is locally stable for \(x_+<1.82\) and unstable for \(x_+>1.82\). In Ma’s birefringent rational-NED model, the fixed-charge heat capacity similarly changes sign at the temperature maximum, and the paper states that the smaller black hole is more stable in the usual fixed-charge ensemble [1803.02191][1509.05580].

A separate line of work emphasizes remnants. In the modified Hayward geometry supported by rational NLED, the temperature and heat capacity vanish at a nonzero radius \(x_h=x_1\), which the paper interprets as a black-hole remnant, while a second-order phase transition occurs at a larger \(x_h=x_2\) where \(C_q\to\pm\infty\). In ECG with small charge and small cubic coupling, the Hawking temperature is regularized and leaves a thermodynamically stable remnant for small \(r_h\neq 0\). This suggests that remnant formation is not tied to a single NED model, although the mechanism is model specific [2108.07175][2305.18254].

Modified frameworks add further layers. In 4D EGB with exponential NED, negative Hawking temperature is used as a diagnostic of parameter regions without physical black holes. In one recent NRCBH model, the tunneling temperature is corrected by a generalized uncertainty principle to
\[
T_{\rm GUP}=\frac{f'(r_h)}{4\pi}\sqrt{1-2\beta m_p^2},
\]
and the same paper studies a Joule–Thomson coefficient \(\mu_J\), finding that cooling behavior is favored for larger charge and larger horizon radius [2005.00360][2507.05145].

## 5. Perturbations, quasinormal modes, and stability

Linear perturbations do not produce a single universal picture across NRCBHs. For Bronnikov’s regular magnetic black hole, electromagnetic perturbations in NED obey a Schrödinger-like equation with an effective potential containing \(\bar{\mathcal{L}}_F\), \(\bar{\mathcal{L}}_F'\), and \(\bar{\mathcal{L}}_F''\). In that model the potential nearly coincides with the Reissner–Nordström electromagnetic potential away from the central region, and the quasinormal spectra are correspondingly very close to RN. The paper finds positive-definite effective barriers outside the event horizon, \(\omega_I>0\) in the convention \(\omega=\omega_R-i\omega_I\), and damped time-domain ringdown, concluding stability against NED electromagnetic perturbations [1409.5988].

The coupled gravito-electromagnetic problem can behave quite differently. In the Fan–Wang families analyzed through Chandrasekhar’s procedure, odd-parity gravitational perturbations couple to even-parity electromagnetic perturbations, and the master system is matrix valued,
\[
\frac{d^2Y}{dr_*^2}+[\omega^2 I - V(r)]Y=0.
\]
For \(\mu=3\) and \(\nu=1,2,3\), all computed modes satisfy \(\operatorname{Im}\omega<0\), so the spacetimes are linearly stable in the explored sector, but the spectra depend strongly on \(q\) and on the family label \(\nu\), and are reported to differ markedly from Reissner–Nordström [2210.00295].

Scalar perturbations in modified gravity preserve the same broad stability conclusion. For the regular magnetic solutions in \(f(R,T)\)+NLED, massive scalar quasinormal modes were computed with sixth-order WKB plus Padé averaging, all imaginary parts were found to be negative, and an eikonal analysis was also given. In that case, increasing the trace coupling \(\beta\) lowers both \(\omega_R\) and \(|\omega_I|\), while increasing the magnetic parameter raises \(\omega_R\) and tends to reduce \(|\omega_I|\) [2312.16614].

Taken together, these studies show that “NRCBH ringdown” is not a single template. Some regular magnetic geometries are almost exterior-isospectral to RN, while others acquire substantial NED-induced gravito-electromagnetic mixing and branch structure. A plausible implication is that observational distinguishability depends more on the specific NED or modified-gravity realization than on regularity alone.

## 6. Optical appearance, lensing, and orbital dynamics

Weak-field and strong-field observables are likewise model dependent. For the Bronnikov regular magnetic black hole, the Gauss–Bonnet and geodesic analyses agree that the weak deflection angle is smaller than in the singular RN case with the same mass and magnetic charge, and a cold non-magnetized plasma increases the bending angle. In the shadow problem, the effective geometry induced by NLED must be used: for both Euler–Heisenberg and Bronnikov models, the shadow was computed from the effective metric rather than from the background metric, and comparison with the EHT image of M87* produced upper limits on magnetic charge, including
\[
Q_m<0.7\,M \quad (1\sigma),\qquad Q_m<1.0\,M \quad (2\sigma)
\]
for the Bronnikov case [2101.08409][1912.08231].

A more recent magnetically charged regular black-hole metric of the form
\[
A(r)=1-\frac{2M}{r}\left[1-\tanh\!\left(\frac{q}{r}\right)\right]
\]
was studied in the thin-disk and EHT context. There the parameter \(q\) decreases the photon-sphere radius and shadow size, lowers \(r_{\rm ISCO}\) and \(r_{\rm mb}\), modestly increases the Novikov–Thorne flux and disk temperature, and changes the widths of direct, lensed, and photon-ring impact-parameter intervals. Using the EHT shadow radii of M87* and Sgr A*, the paper reports
\[
q\lesssim 0.43 \quad \text{for M87* at }1\sigma,
\qquad
q\lesssim 0.32 \quad \text{for Sgr A* at }1\sigma
\]
in units with \(M=1\) [2406.07300].

Not all regular metrics are observationally distinct from RN outside the horizon. For the Balart–Vagenas Fermi–Dirac-based regular solution, the paper finds that photon effective potentials, photon circular orbits, shadow radii, stable circular orbits, and even the regular-versus-chaotic dynamics of charged particles in external magnetic fields are almost indistinguishable from RN outside the event horizon. The paper attributes this to an extremely small exterior difference between the two metric functions. This directly cautions against treating regularity as an automatically observable property [2209.08439].

One 2025 NRCBH study extends the observational program to GUP-corrected thermodynamics, Gauss–Bonnet weak lensing, plasma refraction, orbital frequencies, quasi-periodic-oscillation proxies, and Joule–Thomson expansion. It finds that sufficiently large magnetic charge can drive the weak-deflection series negative, that plasma modifies the deflection through
\[
n(r)=\sqrt{1-\omega_p^2(r)f(r)/\omega_0^2},
\]
and that the Keplerian frequency \(\Omega(r)\) develops charge-sensitive maxima. The same paper explicitly warns that its extremal-charge numerics are internally inconsistent, so its observational trends are more robust than its quoted extremal threshold [2507.05145].

## 7. Extensions, caveats, and no-go results

The NRCBH literature now includes rotating and nonstandard global extensions. Toshmatov, Stuchlík, and Ahmedov constructed rotating regular black holes in GR coupled to NED using a Kerr-like metric with a radius-dependent mass function \(\rho(r)\). Their Maxwellian class approaches Maxwell electrodynamics in the weak-field limit, is regular for \(M=0\) and \(\mu\ge 3\), and yields epicyclic frequencies smaller than in Kerr at the same spin. The same work shows that the weak and strong energy conditions are violated inside the Cauchy horizon of the rotating geometries [1704.07300].

A different extension is the black-bounce realization of magnetically charged regular black holes in Einstein–NLED coupled to a self-interacting phantom scalar. In those solutions the physical spacetime has a minimum areal radius \(R_{\min}=|a|\), the non-Lorentzian region \(r<|a|\) is excised, and the causal structure interpolates among black bounces, one-way wormholes, and traversable wormholes according to the relation between the horizon scale and the bounce scale. These geometries are regular on the physical domain but rely on phantom matter and explicit weak-energy-condition violation in the total stress tensor [2202.02303].

Two major caveats have become central. First, in 4D EGB gravity the original dimensional-regularization proposal is controversial; the magnetic exponential-NED paper explicitly notes that its solution should remain valid only for spherically symmetric spacetimes in the regularized theories, and not beyond spherical symmetry. Second, regularity of the background metric does not guarantee regularity of effective light propagation. The 2025 power-Maxwell analysis shows that the effective photon metric can become singular in electric sectors, and that in particular magnetic configurations nonradial photon trajectories can even become spacelike in the background spacetime. This suggests that causal and optical regularity are logically distinct notions [2005.00360][2506.13676].

The strongest recent limitation is a no-go result for nonminimally coupled electromagnetic fields. In theories containing \(R F_{ab}F^{ab}\), \(R_{ab}F^{a}{}_{c}F^{bc}\), or \(R_{abcd}F^{ab}F^{cd}\), static spherically symmetric magnetically charged regular black holes are generically excluded under mild smoothness assumptions. The obstruction arises from \(r^{-6}\) terms in the trace equation near the center. A combined coupling
\[
\alpha R F+\beta R_{ab}F^{a}{}_{c}F^{bc}+\gamma R_{abcd}F^{ab}F^{cd}
\]
admits a necessary fine-tuning condition
\[
12\alpha+3\beta+2\gamma=0
\]
to cancel the leading magnetic divergence, but this is not shown to be sufficient for a regular solution. A plausible implication is that minimally coupled NED remains the comparatively robust arena for NRCBH model building, whereas nonminimal electromagnetic-curvature couplings sharply constrain the existence of physically realistic regular magnetic cores [2606.13773].

Source: https://www.emergentmind.com/topics/regular-magnetically-charged-black-holes-nrcbhs