2
principle, to investigate the nature of electromagnetism (linear or not), due to different signatures in certain neutrino
phenomena, such as neutrino oscillations, spin-flip and r-processes. The effect of non-linear phenomena on the BH
shadow, BH thermodynamics, deflection angle of light and also wormholes have been investigated too [56–69]. In the
context of primordial physics, instead, NLED, when coupled to a gravitational field, can give the necessary negative
pressure and enhance cosmic inflation [70] and some models also prevent cosmic singularity at the big bang [71–75] and
ensure matter-antimatter asymmetry [76]. The reason to consider NLED in the primordial Universe comes from the
assumption that electromagnetic and gravitational fields were very strong during the evolution of the early universe,
thereby leading to quantum correction and giving birth to NLED [77,78]. Recently, the non-linear electrodynamics
has been also invoked as an available framework for generating the primordial magnetic fields (PMFs) in the Universe
[79,80]. The latter, indeed, is a still open problem of the modern cosmology, and although many mechanisms have been
proposed, this issue is far to be solved. Seed of magnetic fields may arise in different contexts, e.g. string cosmology
[81], inflationary models of the Universe [82,83], non-minimal electromagnetic-gravitational coupling [84,85], gauge
invariance breakdown [83,86], density perturbations [87], gravitational waves in the early Universe [88], Lorentz
violation [89], cosmological defects [90], electroweak anomaly [91], temporary electric charge non-conservation [92],
trace anomaly [93], parity violation of the weak interactions [94]. The current state of art points to an unexplained
physical mechanism that creates large-scale magnetic fields and seems to be present in all astrophysical contexts.
They might be remnants of the early Universe that were amplified later in a pregalactic period, according to one
idea. To create such large-scale fields, super-horizon correlations can only still be created during inflationary epochs.
However, it is still unclear how the electromagnetic conformal symmetry is broken. Different theoretical techniques
have been taken into consideration for this, most notably non-minimal coupling with gravity, which by its very nature
broke conformal symmetry ([95] and reference therein). In a minimal scenario, electromagnetic conformal invariance
can also be overcome. In this instance, the major goal is to modify the electromagnetic Lagrangian to a non-linear
function of F.
= (1/4)Fµν Fµν , as done in [79,80,96].
Since all NLED models significantly depend on scale factors (dimensionless or not), which may cause overlaps
with other physics observables, it is obvious that determining the presence of non-linear phenomena is not free
of uncertainty. Energy extraction from black holes, which is connected to various significant astrophysical events,
including black hole jets and therefore Gamma-ray bursts (GRBs), is one area where NLED effects have not yet been
properly studied [97]. The Blandford-Znajeck (BZ) process [98–103] and the (very recent) magnetic reconnection
mechanism [104,105] are the two different energy extraction techniques used today, along with a revised version of
the original Penrose process [106] called magnetic Penrose process [107–109]. Among them, the BZ mechanism is still
the most widely accepted theory to explain high energy phenomena [110,111] (even if there are still open questions
in certain models or combinations [112–114]). It involves a magnetic field generated by the accretion disk, whose
field lines are accumulated during the accretion process and twisted inside the rotating ergosphere. Charged particles
within the cylinder of twisted lines can be accelerated away from the black hole, composing the jets. A characteristic
feature of this mechanism is that the energy loss rate decays exponentially. This has been confirmed in a good fraction
of observations (X-ray light curves) of GRBs [115]. Furthermore, black holes with brighter accretion disks have more
powerful jets implying a correlation between the two. Even if accretion onto a black hole is the most efficient process
for emitting energy from matter it is not able to reach the energy rate of the GRBs, while other energy extraction ways
such as the Hawking radiation give predictions on temperature, time-scale and energy rate highly in conflict with the
observations [116]. Numerical models of black hole accretion systems have significantly progressed our understanding
of relativistic jets indicating two types of jets, one associated with the disc that is mass-loaded by disc material and
the other associated directly with the black hole [117]. In the first case, however, jets with high Lorentz factors are
not supported. The BZ process, which produces highly relativistic jets by electromagnetically extracting black hole
spin energy, remains the most astrophysically plausible mechanism to do so and is in good agreement with direct
observations [118]. In this sense, understanding the general relativistic magnetohydrodynamic (GRMHD) model of
the bulk flow dynamics near the black hole (where relativistic jets are formed) is essential to study the central engine.
In this paper, in order to determine if non-linear effects may change the rate of energy extraction and the magnetic
field configuration surrounding a (non-charged) black hole encircled by its magnetosphere, we will investigate the
Blandford-Znajek mechanism in the context of the NLED framework.
The layout of the paper is as follows: in Sec. II we derive, for the first time, the general version of energy flux up to
second order in the spin parameter. Sec. III is devoted to computing and solving the magnetohydrodynamic problem
in Kerr-Schild coordinates, searching, in particular, for separated (monopole and paraboloid) solutions. In Sec. IV
we give some estimates of the energy extraction w.r.t. standard BZ mechanism. We study primordial magnetic fields
from (minimally coupled) NLED for different non-linear models in Sec. V, while discussion and conclusions are drawn
in the Sec. VI. In this work, we adopt natural units G=c= 1 and for simplicity set M= 1 in order to handle
adimensional quantities (r,a,...). The negative metric signature (+,−,−,−)is also adopted.