50 years of energy extraction from rotating black hole: revisiting magnetic Penrose process

Magnetic Penrose process (MPP) is not only the most exciting and fascinating process mining the rotational energy of black hole but it is also the favoured astrophysically viable mechanism for high energy sources and phenomena. It operates in three regimes of efficiency: low, moderate and ultra, depending on the magnetization and charging of spinning black holes in astrophysical setting. In this paper we revisit MPP with a comprehensive discussion of its physics in different regimes, and compare its operation with other competing mechanisms. We show that MPP could in principle foot the bill for powering engine of such phenomena as ultra-high-energy cosmic rays, relativistic jets, fast radio bursts, quasars and AGNs etc. Further it also leads to a number of important observable predictions. All this beautifully bears out the promise of a new vista of energy powerhouse heralded by Roger Penrose half a century ago through this process, and it has today risen in its magnetically empowered version of mid 1980s from a purely thought experiment of academic interest to a realistic powering mechanism for various high-energy astrophysical phenomena.


Introduction
Among the experimental tests successfully passed by Einstein's general theory of relativity are the tests of precession of Mercury's perihelion, deflection of photons by Sun's gravity, measurement of gravitational redshift, orbital decay of binary pulsars [4], direct detection of gravitational waves [5], investigation of properties of the Galactic center supermassive black hole [6], and others. So far, all experimental tests of general relativity at various scales and regimes bear no convincing evidence of any deviation of black holes from the rotating Kerr black hole hypothesis, which pronounces that an astrophysical black hole can be well characterized solely by two parameters, namely its mass M and spin a.
The mass is the most fundamental parameter of a black hole, which in many cases can be measured with relatively high precision through observations of dynamics of nearby objects. For example, the most recent estimates of the mass of order 4.14 × 10 6 M of the currently best known black hole candidate SgrA* located at the center of our Galaxy have been achieved by near infrared observations of S2 star revolving around central black hole [7]. In order to measure spin of a black hole one needs to probe effects occurring in strong gravity regime, as its gravitational contribution has no Newtonian analogue and hence it is very small and hard to measure. Acceleration of interstellar matter floating towards a black hole gets heated up, resulting in X-ray emission from accretion disk or hot spots. Some progress on the spin determination methods has been achieved with observations and modeling of X-ray spectra from both stellar mass and supermassive black holes [8]. In addition, potential detection of gravitational waves from extreme mass ratio inspirals (EMRIs) by future space-based Laser Interferometer Space Antenna (LISA) [9] seems to be a promising avenue for determination of spin of astrophysical black holes. However, since method of observation is model dependent and it cannot be directly measured, the estimated spin values in various models may differ dramatically.
According to the no-hair theorem, there can exist the third black hole parameter, electric charge arising from Einstein-Maxwell equations for rotating charged mass. The charge parameter of a black hole is usually set explicitly equal to zero, which is justified by quick discharge of any net charge of a black hole due to selective accretion of a plasma matter surrounding any astrophysical black hole. However, as black holes are usually embedded into external magnetic field arising due to plasma dynamics, and more specifically twisting of magnetic field lines due to the frame dragging of effect in the vicinity of a rotating black hole induces electric field in both vacuum and plasma cases [10,11]. It posits a net quadrupole charge on the black hole [12,13]. This charge is weak in the same sense as magnetic field, i.e. its stress-energy tensor does not alter spacetime metric. Thus the assumption of the Kerr hypothesis is well founded. However, it would turn out as we will show later that black hole charge would play crucial role in making energy extraction process ultra-efficient, so much so that efficiency could range over 10 10 . In the next section we shall describe the black hole charging mechanisms in astrophysical context and discuss its possible screening by surrounding plasma.
Spacetime around an astrophysical black hole is described by the Kerr metric in the standard form given in the Boyer-Lindquist coordinates, where Σ = r 2 + a 2 cos 2 θ and ∆ = r 2 − 2Mr + a 2 . Physical singularity occurs at the ring r = 0, θ = π/2. The roots of ∆ = 0 define outer and inner horizons located at r ± = M ± (M 2 − a 2 ) 1/2 . It is the outer horizon which is referred as the event horizon r + ≡ r H . The geometry is characterized by existence of the two Killing vectors, timelike, δ/δt and spacelike, δ/δφ, indicating the corresponding conserved quantities, energy E and angular momentum L. One can introduce an observer with timelike four-velocity and zero angular momentum L = u φ = 0, which is infalling into the black hole from rest at infinity. This corresponds to the locally non-rotating frame of reference (LNRF) of the zero angular momentum observers (ZAMO) with the four-velocity given by n α = (n t , 0, 0, n φ ), (n t ) 2 = g φφ Computing the angular velocity of LNRF/ZAMO we get Since Ω/(2Ma) is always positive, ZAMO co-rotates with the black hole being dragged.
One of the most interesting properties of the Kerr black hole geometry is the existence of direct analogy between event horizon area A H of a black hole with thermodynamical entropy S H (see, e.g. [14][15][16]), which implies that a black hole of mass M and spin a has irreducible energy, For extremally rotating black hole this is 71% of its total energy [17,18], while the rest of 29% is the rotational energy and is thus available for extraction. For stellar mass black holes this energy is of order 10 63 eV, while for supermassive black holes of mass M = 10 9 M it is of order 10 74 eV making them the largest energy reservoirs in the Universe. It is therefore most pertinent to tap this enormous source most effectively and ultra efficiently. In this study we shall entirely address to classical black holes and their physics and astrophysics without any reference to quantum effects, such as the Hawking radiation etc. First attempt to tap energy from black hole was made by Roger Penrose in 1969 [1] who pointed out existence of negative energy states of particles orbiting black hole with respect to observer at infinity. Negative energy orbits for neutral particles can exist inside the ergosphere. Static observer with timelike trajectory, having spatial velocity v = 0, needs to satisfy the inequality g tt u t u t < 0 implying g tt < 0. This means static observers cannot exist when g tt turns positive; i.e. they can only exist for r > r stat (θ) ≡ M + (M 2 − a 2 cos 2 θ) 1/2 . The ergosphere is the region bounded by event horizon, r H and static surface r stat (θ), so that r stat (θ) ≥ r H . Energy of a particle with momentum p α measured by an observer of velocity u α (obs) is E = −p α u α (obs) . Since no observer can remain static below r stat (θ), i.e. u α turns spacelike, energy E relative to an observer at infinity can turn negative for some suitable particle parameters. However local energy would now be defined relative to stationary observer -locally non-rotating or zero angular momentum observer, that has radial velocity zero but angular velocity is necessarily non-zero, ω = −g tφ /g φφ , the frame dragging velocity. Thus energy relative to LNRO/ZAMO is conserved and would however be always positive while it could be negative relative to observer at infinity. Thus in ergosphere where g tt > 0, there exist particle orbits of negative energy states relative to infinity. This is the key property that drives Penrose process of energy extraction from a rotating black hole.
Following the original idea of Penrose let us consider the equatorial motion of a freely falling particle 1 which decays inside the ergosphere into two fragments one of which (particle 2) attains negative energy relative to infinity, while the other one (particle 3) escapes to infinity with energy grater than that of incident particle. The efficiency of the process, defined as ratio of extracted to infalling energy, and is given by the relation (see, e.g. [19,20] and derivation in Section 3.3) For split occurring close to horizon, it is then given by It is maximal for extremally rotating black hole (a = M), being η PP = 0.207, or ≈ 21%. While for moderate spins, e.g. a = 0.5M, PP efficiency is < 2%. In addition to low efficiency, PP is inoperable in realistic conditions because threshold relative velocity between two fragments after split required to be greater than half of the speed of light [19,20]. This is the condition for a particle to ride on a negative energy orbit relative to an observer at infinity. There exists no conceivable astrophysical process that could almost instantaneously accelerate neutral particle to such high velocities. Nor there exists any observational evidence to support such a happening. Thus PP though very novel and purely relativistic in nature cannot be astrophysically viable.
In the mid 1980s, PP was revived astrophysically by inclusion of interaction of matter with electromagnetic field surrounding black hole [2,3,21], nicely reviewed in [22]. This was magnetic version of Penrose process (MPP), where the inconvenient relative velocity threshold between two fragments after split could be easily overcome in presence of an external magnetic field in which black hole is immersed. In other words, energy required for a particle to get onto negative energy orbit could now come from particle's interaction with electromagnetic field leaving mechanical velocity completely free. It was then shown that the process turned very efficient and its efficiency could even exceed 100%. For example for electrons around stellar mass black hole, efficiency is greater than 100%, for as low a field as milliGauss [23]. Although energy extraction efficiency 100% was first shown in 1985 [2] for discrete particle accretion and idealized magnetosphere, it is highly gratifying to see that recent fully general relativistic magnetohydrodynamic (GRMHD) simulations [24,25] have wonderfully borne out this most important and interesting feature of the process.
Earlier in 1975, Remo Ruffini and James H. Wilson have considered the process of energy extraction from rotating black hole based on the charge separation in a magnetized plasma accreting into Kerr black hole [26]. This became the first use of gravitationally induced charge of a black hole by the frame-dragging effect for the extraction of rotational energy from black holes. Another similar process of energy extraction from rotating black hole is the Blandford-Znajek mechanism (BZ) [11]. BZ operates on the principle of unipolar generator, similar to classical Faraday disc. Here the role of disc is played by a black hole rotating in magnetic field. As in the case of MPP, black hole's rotation generates electric currents along the horizon surface which convert mechanical spin energy of a black hole into electromagnetic energy to be extracted. In addition, both BZ and MPP act due to existence of quadrupole electric field, being produced by twisting of magnetic field lines. However, BZ requires force-free magnetosphere which can be formed e.g. by a cascade of electron-positron pairs [27]. This requires the threshold magnetic field of order 10 4 G. In addition BZ cannot provide ultra high efficiency, which is the distinguishing feature of MPP [23]. Below we will show that MPP is a general process which includes and approximates to BZ for high magnetic field regime. As mentioned before, MPP operates in three regimes of efficiency: low, moderate and ultra high. BZ is included in the middle -moderate efficiency regime. MPP has thus become one of the leading processes for powering the central engine of high energy astrophysical objects like quasars and AGNs involving black holes.
Following the work of Roger Penrose, several other modifications and variants of the original PP were proposed. These included collisional Penrose process (CPP) [28,29] and its various variants [30,31] based on multiple collisions of particles within the ergosphere, whose energy in the center-of mass could grow arbitrarily high for extremally rotating black hole (see, recent review [32] and references therein). It was however agreed that CPP was unlikely to be relevant in realistic high-energy phenomena, since the efficiency of the process in astrophysical situations (i.e. non-extremal black hole, collision occurring not exactly on horizon, incident particle falls from infinity) was severely constrained with maximum η CPP < 15 [33]. In [34], it has been shown that the efficiency of energy extraction from Kerr naked singularity can reach 157%. In [35], energy extraction from boosted black holes, i.e. black holes moving at relativistic speeds, has been studied. Among similar high-energy processes which have drawn attention in the literature is also the so-called BSW mechanism [36], which can provide arbitrarily high center-of-mass energy for collision again occurring at the horizon of extremally spinning black hole. It has also been generalized to include magnetic field [37] and to many other cases of modified gravity (see, e.g. [38][39][40][41][42][43][44][45], among others). This paper is organized as follows. In Sec. 2 we discuss the mechanisms of magnetization and charging of a rotating black hole in astrophysical context. In Sec. 3 we describe the formalism of MPP, study its operation in different regimes and estimate the efficiency of MPP in several plausible radioactive decay modes. Here we show that the efficiency in some cases can exceed 10 10 in astrophysically reasonable conditions. In Sec. 4 we relate MPP for the explanation of ultra-high-energy cosmic rays (UHECRs) and relativistic jets. We show that the cosmic rays with the highest detected energy in the range 10 18 − 10 20 eV can be explained by neutron beta-decay in ergosphere of SMBH of mass 10 9 M and magnetic field 10 4 G. MPP can also provide possible explanation of the knee of the cosmic ray spectrum as the energy of proton after beta-decay in the ergosphere of Galactic center black hole reaches ∼ 10 15.5 eV. We also discuss how MPP can be related to the explanation of relativistic jets from black holes and AGNs. Employing different processes, such as the charge separation in a plasma, pair production, ionization of accretion disk and chaotic scattering we discuss the jet-like motion of a matter, in particular the possibilities to obtain high Lorentz factors and strong collimation. In Sec. 6 we summarize the main results and give concluding remarks.

Magnetization of black holes
Black holes are indeed embedded into external magnetic fields which can arise due to dynamics of surrounding plasma, e.g. electric currents inside accretion disks, intergalactic or interstellar plasma dynamics. Magnetic field could also be generated in early phases of expansion of the Universe. Strength of magnetic field can vary for each particular black hole candidate, usually being in the range 10 -10 8 G, highly dependent on properties of surrounding plasma. For example, the best known black hole candidate Sgr A*, which is located at the Galactic center is surrounded by highly oriented magnetic field measured at distance of few Schwarzschild radii from the center has field of order 10 − 100 G [46,47]. Recent detection of orbiting hot-spots around the object [48] indicates presence of strong poloidal magnetic field at the ISCO scale.
It is important to note that any astrophysical magnetic field is weak in a sense that its energy-momentum tensor does not modify the background Kerr metric. It is easy to see by comparing magnetic field energy in a given volume with that of black hole's mass energy. This condition for stellar mass black holes [49] is given by Observations of various black hole candidates and astrophysical phenomena occurring in their vicinity indicate that inequality (7) is perfectly satisfied. This implies that an astrophysical black hole is weakly magnetized, hence its effect on neutral test particle dynamics is negligible. On the other hand its effect for motion of charged particles is non-ignorable -rather immense, as ratio of Lorentz to gravitational force would be very large due to large value of charge to mass ratio. Since matter surrounding black hole is usually moving with relativistic velocities and highly ionized, one can characterize relative influence of Lorentz to gravitational force by dimensionless parameter B = |q|BM/(mc 4 ). For electrons close to event horizon of the Galactic center black hole [6], the estimate of this parameter is For protons this ratio is ∼ 2000 times lower. Stellar mass black holes, e.g. in binary systems can attribute magnetic fields of order 10 8 G [50], for which B is of similar order of magnitude as (8). This implies that the effect of magnetic field on dynamics of charged particles is very dominant in realistic astrophysical conditions. In the Kerr geometry, it is natural to assume that external magnetic field would also share symmetries of stationarity and axial symmetry. Using the Killing equation ξ α;β + ξ β;α = 0 one finds solution for electromagnetic field in the form [51] A First solution of Maxwell equations in background rotating black hole spacetime corresponding to black hole embedded in homogeneous magnetic field aligned with the spin axis was obtained by R. Wald in [10], while for magnetic field inclined to the axis of rotation was due to J. Bičák and V. Janiš [52]. Magnetic field of current loop around rotating black hole corresponding to dipole magnetic field configurations was solved by J. Petterson in [53].
As we would show later that orientation of magnetic field configuration plays rather secondary role in the context of MPP operation and its efficiency. The four-potential of asymptotically uniform magnetic field in which the Kerr black hole is immersed reads [10] One can notice that rotation of black hole generates quadrupole electric field given by A t which is the result of twisting of magnetic field lines -the frame-dragging effect. Electric field due to this induced charge is for arbitrary magnetic field configuration in the vicinity of axially symmetric black hole. This induced quadrupolar charge may be referred as black hole charge [13,54] and it is this which is responsible for providing necessary negative energy to particle in the ergosphere. By this way, the inconvenient and unsurmountable condition on relative velocity for the original mechanical PP could be easily overcome. One should note that above is true until the extremality of black hole spin is reached, where the gravitational analogue of Meissner effect comes into play. Extremally rotating black hole expels fields out acting thus like a superconductor [10,55], which happens exactly at a/M = 1. On the other hand, it is generally assumed that the maximum plausible spin of astrophysical black holes is a/M = 0.998 [56] marking equipartition of magnetic and gravitational energy. Thus, the further study and discussion is well in order and well motivated.

Black hole charge
Assumption of electrical neutrality of a black hole in many cases is justified by the presence of a plasma around black hole which can quickly discharge any charge excess. Indeed, one can estimate the discharge timescale of maximally charged black hole using the following arguments. Maximal theoretical value of the charge of a black hole of mass M is given by Q = 2G 1/2 M, written in Gaussian units. For spinning black hole it's of order [13] If such a charge is carried by protons and electrons, this corresponds to the mass M Q = m p,e Q/e, where p and e denote proton and electron, respectively. Luminosity of black hole surrounded by plasma or accretion disk can be derived from infalling matter as L = Ṁ c 2 , whereṀ accretion rate and is the fraction of the rest mass energy radiated away. On the other hand, from the balance of gravitational force and radiation pressure in the vicinity of a black hole one can derive the Eddington luminosity for fully ionized hydrogen plasma surrounding a black hole in the form Defining charged matter accretion rate as the fraction of total accretion rate,Ṁ Q = δ ·Ṁ, we get the neutralization timescale of the maximally charged black hole as which is estimated for positive black hole charge. In case of negative charge, the timescale is ≈ 1835 times lower. Both and δ have values in the range (0 , 1) and in many cases are of similar order of magnitude. This implies that in all astrophysically relevant settings any net charge on black hole would be neutralized relatively quickly unless there is a mechanism preventing the black hole from neutralization. For a spinning black hole immersed in an external stationary magnetic field it is easy to see from (9) and (10), that any local observer within the ergosphere can measure nonzero electric field component if the magnetic field has nonvanishing poloidal component [10,26,57]. Black hole spin contributes to the Faraday induction generating the electrostatic potential A t , which can be associated with induced electric field. The charging process is similar to the classical Faraday's homopolar generator. In case of the uniform external magnetic field, there exists the potential difference between event horizon and infinity taking the form case the black hole acquires net charge Q = 2aMB. The process of charging of a black hole due to spin-induced effect in magnetic field takes place in any magnetic field configuration which shares the symmetry of background Kerr metric spacetime. On the other hand, chosen form of the field configuration may put restrictions on the dynamical timescales of the process of selective accretion due to following charge separation in a plasma. In general, the energy of charged test particle is given by . Difference between electrostatic energy of a charged particle at the event horizon and at infinity is given by E H − E ∞ = qA t|r→r H − qA t|r→∞ ≡ δ. For positive δ, more favourable is the accretion of particles with the sign of charge of q, while for negative δ, it is more likely to accrete particles with −q sign of charge. In both cases this leads to the formation of the net black hole charge with a sign which depends on the orientation of magnetic field with respect to the rotational axis of a black hole [58]. For a magnetic field generated by the dynamics of co-rotating surrounding plasma matter the black hole's charge is more likely to be positive. In [13] it was shown that in realistic cases applied to the Galactic center black hole even a small charge of the black hole can have non-negligible effects on the observed bremsstrahlung emission profile. We would like to emphasize that the black hole charge plays a key role in the black hole energy extraction processes, such as Blandford-Znajek mechanism and magnetic Penrose process [23]. Discharge of induced electric field by oppositely charged accreting matter drives the rotational energy away from the black hole in both of the processes.
However, in order to tap the gravitationally induced electrostatic energy from black hole and support particle acceleration, this electric field should not be screened as usually occurs in the presence of a plasma. In the past the problem has been widely discussed [59,60] resulting with the Komissarov's theorem [61] which states that induced electric field of rotating black hole in magnetic field is not screened at least within ergosphere. Indeed, total screening of black hole's charge may occur only when the following two conditions are satisfied simultaneously where B and D are magnetic and electric field measured LNRO or ZAMO. It is easy to show that the relation B 2 − D 2 which is positive far away from black hole, is negative everywhere inside the ergosphere 1 [61]. Moreover the sign of this relation is independent of the strength of magnetic field being only dependent on the location. In other words, electric field is stronger for stronger magnetic field and within the ergosphere it cannot be screened off. Black hole charge associated with gravitationally induced electric field has been estimated in [13] for the Galactic center supermassive black hole Sgr A*, with an upper limit of 10 15 C. On the other hand, classical estimates of charge, based on the difference between thermal velocities and masses of electrons and protons in the fully ionized plasma around Sgr A*, imply the presence of equilibrium charge of the central body of order 10 8 C. It therefore follows that black holes posses an electric charge in the range 10 2 − 10 12 C per solar mass. This charge is gravitationally weak in a sense that its influence to the spacetime metric can be neglected. For it to be gravitationally significant would require charge of order ∼ 10 20 C per solar mass, and therefore, Kerr black hole hypothesis stands firm and valid. It is important to note that such a charge cannot be measured by imaging of black holes, based on the observations of their shadows. However, even such a small charge associated with the black hole has significant effects on the processes occurring in its neighbourhood, such as acceleration of charged particles to ultra-high energy.

Dynamics of charged particles around black hole
Efficiency of MPP in its original formulation can be derived by applying the conservation law for parameters of the particles before and after decay near black hole. Motion of relativistic charged particles, such as electrons and protons in magnetic field leads to inevitable synchrotron radiation loss which can sufficiently influence particle's energy and angular momentum. Therefore it is important to establish the limits of applicability of the formalism based on the conservation laws. In this subsection we provide general description of the charged particle dynamics in curved spacetime taking into account the influence of radiation-reaction forces and providing corresponding estimates. Detailed studies of charged particle dynamics around Schwarzschild and Kerr black holes in magnetic field neglecting radiation reaction were done by several authors, see, e.g. [58,[62][63][64]. Qualitative and quantitative studies of synchrotron radiation reaction problem in curved background can be found e.g. in [65,66] and in more recent papers [67,68]. Equation of motion for a point charge in its most general form is usually referred as the DeWitt-Brehme equation [69], which contains non-local tail integral term and Ricci term. Detailed discussion of the complete set of equations of motion for radiating charged particles can be found in [70,71]. It was shown, however, that the tail force can be neglected in most of the cases [68], being for electrons around stellar mass black hole 10 −19 times smaller than corresponding gravitational "force" on the horizon scale. Ricci term is irrelevant in vacuum metrics. Thus, after several algebraic manipulations and applying 1 As for energy E, sign of B 2 − D 2 would dependent upon the location of observer whether inside or outside ergosphere. There could exist no static observer in ergosphere, it could at best be stationary with non-zero frame dragging angular velocity.

B (Gauss)
τ e (s) τ p (s) τ Fe (s) 10 12 10 −16 10 −6 10 −5 10 8 10 −8 10 2 10 3 10 4 1 10 10  Landau-Lifshitz method [72], one can simplify the equation of motion for radiating charged test particle around black hole to the following form [68] Du where u α is four-velocity satisfying normalization u α u α = −1 and F µ ν is the Faraday tensor of external electromagnetic field, whose covariant derivative is given by Equation (17) is thus the covariant form of Landau-Lifshitz equation. Solving this equation numerically gives us the rate of energy and angular momentum loss. For our purpose here it is important to find the cooling timescales corresponding to energy and angular momentum loss of charged particles around black hole. Therefore, applied to the Schwarzschild black hole with uniform magnetic field we find the evolution of specific energy E and specific angular momentum L in the form where f = 1 − 2M/r and k = 2q 2 /(3m), and q, m are mass and charge of test particle. For ultrarelativistic charged particle (E 1, or B 1) the leading contribution for energy loss is the first term on the right hand side of Eq. (19). Inserting all constants, the cooling timescale of charged particle is given by In this case the Lorentz force is dominant over gravitational "force" which happens in realistic settings as discussed above Eq. (8). Closer to black hole the cooling timescale increases. In Tab.1 we give typical cooling timescales of electrons, proton and fully ionized iron nuclei for various values of magnetic field. Due to cubic dependence on mass, electrons cool 10 10 times faster than protons. One can compare the cooling timescale with an orbital timescales of particles at ISCO τ orb ≈ 4πr isco /c, which is of order of ∼ 10 −3 s for non-rotating stellar mass black holes or ∼ 10 5 s for supermassive black holes. Thus, the energy loss can be quite relevant especially in case of lighter particles, such as electrons.
If the radiation-reaction force is neglected, two components of generalized four-momentum P α = mu α + qA α are conserved, namely energy and angular momentum of a particle. They can be associated with the Killing vectors in the form In this case the standard approach can be used (see, e.g. [58] and references therein). The motion of charged particle in magnetic field is always bounded, which can be described by introducing the effective potential. In the equatorial motion it takes the form One sees that V eff can attain both positive and negative values depending on the angular momentum, charge and location of the particle. Occurrence of negative energy orbits (NEOs) in the ergosphere is critical requirement for extraction of energy from rotating black holes. Presence of the term −qA t in the effective potential extends the region of existence of NEOs far beyond the ergosphere, while the last two terms in V eff can be negative only within the ergosphere. In fact, the region with possible NEOs and thus, the energy extraction zone for charged test particles extends to infinity [73,74]. Ergosphere is maximal in the equatorial plane, therefore we consider decay of a particle falling onto black hole in the equatorial plane. In addition to energy and angular momentum and their conservation due to the Killing symmetries, the normalization condition u α u α = −k must be satisfied for both charged and uncharged particles, where k = 1 for massive particle and k = 0 for massless particle. In the equatorial motion with the four-velocity u α = u t (1, v, 0, Ω), where v = dr/dt and Ω = dφ/dt we get the angular velocity of a test particle with respect to the asymptotic observer at rest in the form where the sign defines the co or counter rotation with respect to LNRO. The limit of u α tending to a null vector gives the restrictions to the angular velocity of a particle (both charged and uncharged) surrounding black hole in the form [3]

Split of infalling particle
Let us now consider decay of a particle 1, which is not necessarily neutral, into two charged fragments 2 and 3 close to horizon in the equatorial plane. According to conservation of energy and angular momentum after decay, we write where a dot indicates derivative relative to particle's proper time. If energy of particle 2 is negative relative to infinity, then particle 3 attains energy E 3 = E 1 − E 2 > E 1 greater than that of incident particle 1. Infalling negative energy into the black hole results in extraction of its rotational energy.
Using the conservation laws, one can show that the local angular velocities of particles at the point of split satisfy the relation Reminding that u φ = Ω u t = −Ω e/X, where e = (E + qA t )/m and X = g tt + Ω g tφ , Eq.(36) can be rewritten in the form After several algebraic manipulations we get the energy of escaping particle in the form where Ω i = dφ/dτ is the angular velocity of i-th particle, which is given by Eq. (26) and restricted by the limiting values (29). If the parameters are chosen in such a way that at the point of split q 3 A t < 0, then this term plays the dominant role in the energy extraction from black hole. If q 3 > 0, and B and a are positive, it is easy to see that the condition q 3 A t < 0 is perfectly satisfied.

Three regimes of MPP
MPP can operate in three regimes providing low, moderate and ultra high efficiency for the energy extraction from black hole, depending essentially on strength of magnetic field. As in the neutral case, Eq. (5), we define the energy extraction efficiency as the ratio between gain and input energies, i.e. in our notation Using Eqs. (38) and (39) at the point of split the general expression for the efficiency reads where A t is calculated at the point of split. Setting Ω 1 to Ω given by Eq. (26) and Ω 2 = Ω − , Ω 3 = Ω + , which maximizes the efficiency, and reminding that the velocity component u t is related to energy E as mu t = −(E + qA t ), we get which for freely falling massive particle (k = 1) reduces to χ = (1 + 1 + g tt )/2. The expression (41) can also be rewritten as where is the efficiency of original Penrose process leading to Eqs. (5) and (6). In the absence of magnetic field, MPP turns to original Penrose process giving its lower limit with the maximum efficiency η low MPP ≡ η PP = ( √ 2 − 1)/2 ≈ 0.207 or 20.7%, corresponding to extremally rotating black hole. If all particles are charged the efficiency is given by its full form (43). However, as we have shown before in the presence of magnetic field for elementary particles the electromagnetic forces are dominating in the system, implying that | q m A t | |u t |. Relating this to (43) one can simplify the expression to the form This is the moderate regime of MPP, which can operate, as we see, when q 3 > q 1 , thus neutralizing gravitationally induced electric field of black hole. MPP in moderate regime has direct analogy with another famous process, namely Blandford-Znajek mechanism (BZ). In both cases the driving engine is a quadrupole electric field of a black hole which arises due to twisting of magnetic field lines by the frame-dragging effect. It is important to note, that the efficiency of BZ and moderate regime of MPP cannot grow ultra-large due to natural restrictions of global plasma neutrality surrounding black hole. In moderate regime MPP approximates to BZ and that explains why ultra-high efficiency has not been observed in numerical simulations of the process. There could as well exist the third and the most efficient regime of MPP, which requires special attention leading to several important predictions. If the particle 1 is neutral q 1 = 0, with energy E 1 = m 1 , splitting into two charged fragments, it is easy to see that the general expression for efficiency (41) reduces to the form The term χ − 1 on the right hand side of (46) is the efficiency of mechanical PP which depends on purely geometrical factors and its value ranges from 0 for a = 0 to 0.207 for a = 1. One can also see, that A t component of vector potential (10) attains negative values everywhere above horizon, if the spin and magnetic field are co-aligned, i.e. aB > 0. Since for elementary particles, such as electrons and protons, the charge-to-mass ratios q/m are typically very large, the dominant contribution in (46) is due to −q 3 A t /m 1 term. Thus, the expression for the efficiency in this case can be rewritten in the form which can grow to enormous values for elementary particles. If the spin and magnetic field are aligned in opposite directions, i.e. aB < 0, then, in order to get positive η ultra MPP , the charge of escaping particle must be negative, q 3 < 0. This causes MPP to turn ultra-high-efficient regime imparting ultra-high energy to escaping particle.
Note, that all expressions for efficiency in three regimes are quite general and independent of magnetic field configuration. We also do not specify splitting point at this stage. Essentially axial symmetry of spacetime and electromagnetic field is what is required for MPP to operate in three regimes of efficiency depending upon the two parameters, magnetic field strength and charge to mass ratio of particles involved in the process of energy extraction. The other factor that matters is split point should be as close to horizon as possible which means infalling particle is neutral so that it can reach closer to horizon without any hindrance and then it splits or decays into charged fragments of opposite charge. So the ultra-high regime is characterized by magnetic field, particles involved are electron or protons, and it is neutral particle that decays closest to horizon.
The field configuration however matters for direction of escape and final fate of escaping particle, which is expected to move along the magnetic field lines. For example, in case of uniform magnetic field the collimation of escaping particles is maximal, while in case of hypothetical magnetic monopole the particles escape isotropically. In case of closed magnetic field lines such as for dipole magnetic field generated by current loops in the accretion disk of black hole, the escaping zone is concentrated on the polar caps collimating bipolar outflow of matter similar to relativistic jets. On the other hand, if energy of escaping particle is ultra-high powered by ultra-high efficient MPP, particle can get across the field lines and escape in arbitrary direction.

Quantitative estimates
In order to compare different regimes of MPP we choose the external magnetic field to be asymptotically homogeneous [10] with the strength B. Thus the four-vector potential A µ around black hole is given by relations (10) and (11). We start from the expression for efficiency of MPP in ultra-high regime, given by Pair production γ 0 → e − + e + e − < 0 e + 2.6 × 10 6 ultra Table 2. Maximum efficiency of magnetic Penrose process and corresponding energy extraction regime for some typical radioactive decay modes in the vicinity of a black hole of mass 10M with the spin parameter a = 0.8M. Magnetic field is aligned along the rotation axis and has the strength 10 4 G. Initial energy of decaying particle is taken to be of order of its rest mass, m X c 2 . In case of the pair production the energy of photon is taken to be E γ = 2m e c 2 with recoil term neglected for simplicity.
where B is the strength of uniform magnetic field, r S = 2GM/c 2 , and r split is the splitting point. The first term on the right hand side is positive only if the splitting point r H ≤ r split < r S with maximum at r split = r H . This, purely geometrical term is due to the original Penrose process and it varies in the range (0 ; 0.21). The second term gives the contribution due to electromagnetic interaction, which can exceed 1 for electrons with a few milliGauss field and stellar mass black hole [23]. The important point to be noted is that electromagnetic interaction effectively expands ergosphere beyond the geometric bound, r = 2M. That is negative energy orbits are available in much enlarged region which goes on to enhance efficiency and overall working of the process. Efficiency of MPP for escaping proton in ultra-high regime after decay of freely falling neutron above horizon is plotted against magnetic field in Fig.1 for different values of black hole mass and spin. The interesting and important that emerges from this is that MPP does not require rapid rotation of black hole and can as well operate for low spins. In Table 2 we compare MPP efficiency in different regimes for some typical radioactive decay modes occurring above horizon of a black hole with mass 10M , spin a = 0.8M and magnetic field of strength 10 4 G. Magnetic field is taken to be asymptotically uniform.

State of the art
One of the fascinating applications of MPP is the explanation of the origin and production mechanism of ultra-high-energy cosmic rays (UHECR). These are extragalactic charged particles observed at energies above 10 20 eV [75,76]. Detection of particles at such a high energy violates the GZK-cutoff limit [77,78], which requires exotic scenario for its production mechanism. There are indications that UHECRs are produced at AGNs evidenced by recent extrasolar neutrino detection by IceCube collaboration and simultaneous multi-wavelength observations in optical spectra tracing the source to blazar [79,80]. The source TXS 0506+056 is located at the distance of ∼ 1.75Gpc from Earth with relativistic jet directed towards us. Recent observational data lean towards the mixed composition of UHECRs [81,82], although all previous fluorescence measurements of cosmic rays demonstrated dominantly proton flux at high energies [83,84]. As observed by the largest cosmic ray telescopes, such as Pierre Auger Observatory [76] and Telescope Array [85], the highest energy cosmic rays (> 10 18 eV) are originated outside of our Galaxy. Spectrum of cosmic rays demostrate the presence of the knee at ∼ 10 15.5 eV and the ankle at ∼ 10 18.5 eV. Flux of particles above knee decreases considerably, which may indicate the change in the source from Galactic to extragalactic. The origin of cosmic rays above the knee and mechanism of production of highest energy cosmic rays remain under active debate. The main difficulty in understanding of UHECR physics is extremely low flux of these particles. The flux of UHECRs with energy > 10 20 eV is about one particle per km 2 per century. One of the promising acceleration scenarios of UHECRs is the shock acceleration in the plasma of relativistic jets (see, e.g. [86]). However, due to extended acceleration length and interaction looses inside and outside of the jet, realisation of this scenario in the nature is debatable. Below we apply MPP for the explanation of UHECRs which in contrast to jet acceleration models does not require extended acceleration zone. For this we employ the beta-decay of free neutron in the vicinity of supermassive black hole immersed into external magnetic field. We will partially follow the results presented in recent paper [87].

Maximum energy of proton
In Figure 2 (left) we illustrate the Feynman diagram of the beta-decay in the ergosphere of rotating black hole which corresponds to MPP in ultra-efficient regime. Neutron decays into proton, electron and antineutrino. Due to expectation of alignment of accretion flow with a black hole rotation at least near the event horizon, gravitationally induced charge on the black hole is more plausibly positive. This implies that in the neutron beta-decay the escaping particle is proton. One should note that the effect of antineutrino on the energy of proton is negligible.
Applying MPP to the beta-decay of free neutron in vicinity of a standard supermassive black hole (SMBH) with mass M = 10 9 M and magnetic field of strength B = 10 4 G, one can get proton escaping from black hole with energy We plot the dependence of proton's energy on magnetic field in Fig.2 (middle) for various values of black hole masses. Spin is chosen to be a = 0.8M as an example. As we have seen in Fig.1, the spin of the black hole does not change the efficiency of MPP dramatically, although higher the black hole spin, larger is the ergosphere, which enhances probability of operation of MPP in reasiltic conditions.

Propagation of cosmic rays in magnetic field
Depending on the constituents of primary cosmic rays, a significant amount of initial energy will be lost in the photo-pion interactions of cosmic ray particles with low energy photons of cosmic microwave background (CMB). Loss caused by such interactions will grow over cosmological distances if the energy of primary cosmic ray is larger than ∼ 10 19.7 eV. This effect of energy loss by primary cosmic rays is known as the GZK-cutoff [77,78]. Observations of cosmic rays above GZK-cutoff limit indicate the location of sources within < 100Mpc. Indeed, the data obtained by Pierre Auger Observatory point to strong correlation between arrival direction of highest-energy cosmic rays with the location of nearby active galactic nuclei within a distance of 75 − 100Mpc from our Galaxy [88].
However, apart from CMB photons the UHECRs interact with magnetic fields along the propagation distance which can vary from intergalactic and Galactic magnetic fields of order 10 −9 − 10 −5 G up to values exceeding 10 4 G in the source region. Therefore it is important to investigate the synchrotron loss timescales and distances for UHECRs in various magnetic field. Cooling timescale of charged particle in curved spacetime is given by (22) which has cubic dependence on mass of charged particle. Hence, electron cooling timescale is ∼ 10 10 shorter than protons, which makes heavier constituents more realistic cosmic ray candidates (see, Table 1). In Fig.2(right) we demonstrate energy loss of protons with various initial energies over propagation distance due to intereaction with Galactic order magnetic field 10 −5 G. One can also see that in the case of neutron star (sometimes suggested as UHECRs source), whose magnetic field is typically above 10 12 G, the timescale of synchrotron energy loss is extremely short even for protons and ions.

Cosmic rays from Galactic center black hole
Remarkably, MPP gives relatively precise prediction on the maximum proton energy produced by Galactic center black hole, which is the closest and best studied supermassive black hole candidate nowadays [6]. The central black hole Sgr A* has relatively accurate measurements of its mass being of ∼ 4.1 million solar masses based on the observations of surrounding stellar dynamics [89]. Distance measurements are closely related to mass estimates which gives a distance of ∼ 8.1kpc from Solar system to Galactic center. Estimate of spin is model dependent and requires observations of immediate regions of black hole. Based on multi-wavelength, radio, infrared and near-infrared interferometric data current constraints on the spin of Sgr A* is a > 0.4M [90], although one cannot exclude near extremal values [48]. Magnetic field around Sgr A* has equipartition strength of 10G [6], which can reach up to few hundred Gauss on the event horizon scale [47]. Applying MPP to the beta-decay of neutron skirting close to Sgr A* we obtain the high-energy proton at the Galactic center with energy of order Surprisingly, this energy has the same order of magnitude as the knee of the cosmic ray energy spectra, above which flux of particles demonstrate significant suppression [91]. Although many models have been suggested, the explanation of the origin of knee is still in strong debate due to number of uncertainties in Galactic magnetic field structure and lack of direct detections of primary cosmic rays. The sharpness of the knee energy spectra could in fact indicate the existence of a single source at knee energy level [92]. Thus, our model could serve as an alternative and relatively simple explanation of the knee appearing in the cosmic ray spectra. The logical extension of the proposed model can be the search of correlations in the observational data of UHECRs with predictions of our model.

Relativistic jets
Relativistic, collimated ejections of matter named jets have been observed in variety of objects, such as AGNs, X-ray binaries, quasars, etc. It is generally believed that the mechanism of production of jets is strongly connected with accretion processes in combined strong gravitational and electromagnetic fields. Being one of the fundamental problems of the modern relativistic astrophysics, the relativistic jets are under intensive considerations. Large interest is connected with attempts to understand how the energy from accreting matter is converted into the kinetic energy of escaping matter in a narrow cone, so that the jets acquire very large gamma-factors.
Here we describe a simple model of jet-like motion of particles due to the ionization of neutral Keplerian accretion disc in the equatorial plane of weakly magnetized rotating black hole. As we have already seen above, the ionization of neutral particles in the vicinity of black hole can provide ultra-high-energy to one of its charged fragments after the split, so the fragment has enough energy to escape from the black hole. However, the motion of charged matter in the presence of magnetic field is always bounded in the equatorial plane (see, e.g. [58] and refs therein). Therefore the escape of particles from black hole is possible only along the magnetic field lines, which we assume to be aligned with the rotation axis of the black hole. This assumption is well justified due to the fact that any field configuration in the vicinity of black hole generally shares the symmetries of the background spacetime metric. Transition from equatorial motion of charged particles to the linear motion along the rotation axis occurs due to chaotic scattering of high-energy particles in the effective potential. This causes the interchange between the oscillatory and transnational modes of energy E 0 → E z of charged particles. The mathematical technique and numerical modelling of the chaotic scattering of charged particles by black hole is given in [93].

Chaotic scattering of ionized particles
Motion of test particle is limited by the effective potential, so that E = V eff . While bounded in equatorial plane in the presence of magnetic field, the boundaries of motion can be open for escape to infinity along perpendicular direction to the equatorial plane. A particle is able to escape to infinity if its energy is grater than asymptotic value, which in terms of the specific energy (E = E/m) is given by where B is defined by (21). One can derive that where E ∞ is the energy of a particle measured at infinity. Thus the total energy of a particle measured at infinity is composed from the longitudinal E z and transverse E 0 parts. It appears that near the black hole two components of energy E z and E 0 are interchangable [93], while in Minkowski spacetime the energies given by (52) -(54) are integrals of motion and therefore cannot be transferred between two energy modes. In Kerr spacetime metric the conserved quantity is the total covariant time component of momentum i.e., E ∞ . Interchange in two energy modes imply the change of the velocity of charged particles along the magnetic field line, which can increase to extremely large values due to combination of MPP with chaotic scattering effect.

Escape velocity
Asymptotic value of Lorentz gamma factor is related to the energy at infinity in the following form Escape velocity u z = dz/(dτ) or v z = dz/(dt) and the Lorentz gamma in z-direction γ z can be expressed in terms of relation (52) as follows The Lorentz factor is maximal when E 0 is minimal, i.e.
In case of maximal acceleration with positive B, the orbital velocity of particle around black hole vanishes and all energy in equatorial plane transfers into kinetic energy in z direction, so that u φ = 0, for B > 0. In case B < 0, for maximal acceleration the orbital velocity tends to the limit u φ = 2BL. Comparison of trajectories in two cases with negative (left) and positive (rights) magnetic parameter B is represented in Fig. 3, which is reproduced from [93]. The trajectories are found in [93] by numerical integration of the full set of equations of motion. Acceleration is larger in case of positive B. The Fig. 3 represents the cases with relatively low magnetic parameters |B| ∼ 1. For elementary particles around astrophysical black holes the magnetic parameter B is usually very large (see, e.g. (8) for electron around Sgr A* black hole). Therefore, the Lorentz factor of escaping particles after decay following MPP can be extremely large. Energy of charged particles produced in MPP is mainly concentrated in E 0 mode due to considered formalism of decay in equatorial plane. The motion in that case is quasi-circular. Due to chaotic scattering all energy in E 0 mode can be transformed to longitudinal kinetic energy E z in the black hole vicinity. Notable feature of the "transmutation" effect is that it does not require black hole's rotation and can operate also in Schwarzschild spacetime immersed into external magnetic field. However, in order to produce high-energy particles the black hole rotation plays crucial role as it generates electric potential A t providing ultra-high acceleration to the charged particles. Thus, the combination of MPP with chaotic scattering effect for the acceleration of charged particles can serve as a simple model of relativistic jets as the model provides extremely large Lorentz factors (γ z 1) and strong collimation of escaping charged particles along the symmetry axis.

Charge separation in the accretion disk
In the plasma setup, similar mechanism of energy extraction from rotating black holes can described by the mechanism of Ruffini and Wilson [26], based on the charge separation in a magnetized plasma of accretion disk around black hole. In this scenario, the black hole may act as a pulsar [54] with non-zero net charge of the black hole and resulting charged magnetosphere. Plasma of accretion disk generally surrounding black holes is usually considered to be neutral. This assumption is justified by the neutralization of any charged plasma configuration in relatively short timescales. However, taking an analogy of magnetosphere of pulsars, where the plasma acquires net charge density, called Goldreich-Julian (GJ) charge density [94], it is easy to show that similar situation occurs in the relativistically moving plasma around black hole. In particular, given that the matter is ordinary two component fluid of hydrogen plasma, motion of plasma in relativistic speeds v in external magnetic field B induces an electric field component D In the frame co-moving with plasma, electric current density is connected with electric field D by the Ohm's law, which takes the following form for an observer at rest Since electric field in co-moving frame must be neutralized, which usually happens in a very short timescales, electric field with respect to observer at rest does not vanish. Velocity of an arbitrary point of the disk is v = × R, where Ω is angular velocity and R is distance to chosen point of the disk. This implies that net charge density in disk is given by Consequently, charge separation in plasma moving around weakly magnetized rotating black hole can be applied to MPP in a similar way as decay of neutral particles into charged fragments, described in previous sections. General relativistic magnetohydrodynamic simulations (GRMHD) of two component plasma fluid with global charge neutrality revolving around rotating black hole in external magnetic field should be tested as it may lead to energy extraction from black hole in ultra efficient way.

Discussion
The high energy phenomena like quasars, AGNs, GRBs, FRBs, UHECRs and so on ask for two things: one, a very large reservior of energy and two, it has to be harnessed very-very efficiently to give enormous power output. A rotating black hole has 29% its mass in rotational energy which could be mined out. It therefore becomes the natural candidate for the huge reservior we are looking for. Now the question is of very efficient tapping of this rotational energy.
In 1969 Roger Penrose proposed a very novel and purely geometric process [1] of energy extraction from a rotating black hole. The key property driving the process was the existence of negative energy orbits relative to an observer at infinity in the vicinity of black hole horizon, the region called ergosphere. At the root of all this is the property of frame-dragging by which black hole shares its rotation with surrounding space. That is, space around a rotating black hole has inherent rotation so as to give angular velocity even to particle having zero angular momentum. The crucial point is that this property of existence of negative energy orbits in ergosphere is at the root of all processes involved in extracting energy from black hole.
The question then boils down to find the most efficient process. In the original Penrose process, the maximum possible efficiency was 20.7%. More important was the question that energy required to push a particle on NEOs is non-trivial -relative velocity between fragments to be greater 1/2c [19,20]. This requirement was astrophysically unsurrountable because there can be no astrophysical mechanism that could accelerate particles to such relativistic velocity almost instantaneously. Thus Penrose process was very ingenious and interesting but it was astrophysically not viable to power high energy sources. Then various variants of the process were considered [28][29][30][31][32][33][35][36][37] but not of much any consequence or avail.
In 1977 came the BZ [11] which is set with black hole sitting in a magnetic field whose field lines are twisted due to frame dragging producing quadrupolar electric potential difference between pole and equator. The discharge of which drives out energy and angular momentum from black hole. It was astrophysically very exciting and interesting process. However NEOs also play the critical role in working of the process. Further the process requires polarization of vacuum in ergosphere which would require threshold magnetic field of order 10 4 G. Here rotational energy of black hole is extracted electromagnetically -magnetic field plays the role of a catalytic agent.
In 1985 an electromagnetic version of Penrose process -MPP [2,3,21,22] was considered and it was shown that the limit on threshold relative velocity could be easily overcome by particle's interaction with electromagnetic field; i.e. energy required for particle to ride on NEO could now come from −qA t leaving particle's velocity completely free. This was a wonderful revival of the process for astrophysical applications. Not only that it was shown that the process is highly efficient, so much so efficiency could exceed 100% for as low field as milli Gauss. This was for the first demonstration of efficiency exceeding 100%.
Clearly BZ and MPP share the same setting of black hole immersed in magnetic field, and so are the two the same? True, in moderately high magnetic fild range, the two seem to be similar. Recently it has been argued [23] that MPP is a general process which works in all magnetic field range while BZ requires the threshold magnetic field, and for that range the former approximates to the latter. That is BZ is contained in MPP in high field range.
All this was of course in the linear test particle accretion regime which is admittedly an idealized situation. More realistic setting could not be taken up due to lack of computing power at that time. Since around 2010, fully relativistic magnetohydrodynamic flow simulation models have been studied [24,25], it is remarkable that the process continues to remain highly efficient, efficiency ranging to 300%. It turns out that MPP/BZ is the most plausible and important powering mechanism for high energy sources like quasars etc. It is remarkable and very gratifying for one of us that the mid 1980s prediction of efficiency exceeding 100% has been beautifully borne out. Thus MPP is perhaps the key process of mining rotational energy of black hole and powering the high energy sources.
Among astrophysical phenomena which can be directly related to MPP are the ultra-high-energy cosmic rays and relativistic jets. In case of UHECRs, MPP is able to produce protons with energy exceeding 10 20 eV employing neutron beta-decay in the ergosphere of supermassive black hole of mass 10 9 M and magnetic field of 10 4 G. Heavier cosmic ray constituents of similar energy range can also be achieved depending on the decay mode of infalling neutral matter. We summarized some of the potential radioactive decay modes with ultra efficient energy extraction in Table 2. Depending on the energy range and constituents of primary UHECRs, this can provide constraints on the black hole's mass, magnetic field and the distance to the source candidate. Remarkably, the knee energy of the cosmic ray spectrum at ∼ 10 15.5 eV coincides with the maximum proton energy obtained by MPP applied to the Galactic center black hole Sgr A*. Correlations of detected cosmic rays with nearby SMBHs of mass and magnetic field predicted by MPP needs to be explored.
Similarly, MPP can provide jet-like behaviour of charged particles, namely the large Lorentz gamma factor and strong collimation of ejected matter. This can be related to the ionization of accretion disk either by decay of infalling neutral particle into charged fragments or due to charge separation in a relativistic plasma of accretion disk. In the latter scenario black hole may behave like a pulsar [54] with net stable charge of both black hole and surrounding magnetosphere with respect to observer at rest at infinity. Ejection of ionized particles to infinity takes place along magnetic field lines, which can be open to infinity at the polar caps of black holes. Then, charged particles escape to infinity due to chaotic scattering effect caused by the interchange between the oscillatory and translational modes of the total energy of escaping particles E 0 → E z [93]. Thus the combination of MPP with chaotic scattering effect leads to the high Lorentz factors and alignment of trajectories of escaping particles along the rotation axis of a black hole. This is quite promising and exciting that MPP may find further astrophysical applications in the near future, such as to be relevant for the operation of many other high-energy astrophysical phenomena, such as e.g. GRBs, FRBs among others.
Finally, we have come long way down from 1969 -good half a century, and it is remarkable to see what was proposed as a thought experiment had come to stay firmly and beautifully as one of the prime movers in high energy astrophysical phenomena. It is an excellent example of purely geometry driven process aided by magnetic field could wonderfully serve as powering engine for such complex astrophysical systems like quasars, AGNs, UHECR, and so on. In this review we have attempted to recount the historical evolution of the process, and it is fair to say that it has more borne out the promise it held.