# On the Dyadosphere of Black Holes

## I Introduction

The title of the original talk I was supposed to give was”on the critical mass”. I was planning to cover the latest developments since my early work on the role of the critical mass in the physics and astrophysics of neutron stars and Black Holes. To discuss as well, how after my propoposal of identifying Cygnus X1 as a Black Hole approximately ten binary systems with very similar characteristic have been found in our own galaxy.

I was also going to stress that very little progress has been made in the intervening years about the determination of the critical mass of neutron stars due to the formidable theoretical difficulties in understanding the behaviour of matter at supranuclear densities and that indeed the theorem about the absolute upper limit to a neutron star mass, established by C. Rhoades and myself, stands still as the pivotal point in the paradigm of identification of Black Holes versus Neutron Stars.

It was my intention to show how, surprisingly, progress is being made recently in gaining a new understanding of the white dwarf critical mass: a classic problem considered by many people well settled, which is extremely attractive from a theoretical point of view. The white dwarfs offer the possibility of testing the concept of critical mass in a regime where all physical theories work in the exact range of their tested validities: there are not unknown due to extrapolation of physical theories to yet unexplored domains as in the matter at supranuclear densities. The new effects can be computed by generalizing to a relativistic treatment the classical results of the Feynmann-Metropolis- Teller work based on the Fermi-Thomas model of the atom to relativistic regimes. The consequences on the estimate of the critical mass can be of paramount importance in the understanding the observed values of the neutron stars in binary systems as well the process of formation of such systems.

However, the recent observations of the Beppo-SAX satellite, the discovery of a very regular afterglow to GRB’s, the fact their x-ray flux varies regularly with time according to precise power laws (Costa et al.; Van Paradijs et al.), and especially the optical and infrared identification of these sources (Kulkarni et al., Halpern et al. and Ramaprakash et al. ) which has established their correct cosmological setting and determined their formidable energy requirements convinced me to reconsider some theoretical results on black hole astrophysics (see e.g. Ruffini et al., Fang, Ruffini and Sato) and to develop a detailed model for a direct confrontation with the observational results.

For these reasons I have deceided, with the consensus of the organizers, to modify slightly the title of my talk to a new talk and deceided to insert this material on the critical mass in a different publication. This will allow me to introduce today a most urgent new concept: the Dyadosphere of a Black Hole and offer this new concept to Humitaka Sato for his sixtyfifth birthday.

## II The basic results on EMBH

It is now generally accepted that Black Holes are uniquely characterized by their total mass-energy $`E`$, charge $`Q`$, and angular momentum $`L`$. The uniqueness theorem has finally been proved after some twenty five years of mathematical work by numerous authors recently reviewed in ref..

I assume that the process of gravitational collapse for a core larger than the neutron star critical mass will generally lead to a black hole characterized by all three fundamental parameters of mass-energy $`M`$, angular momentum $`L`$, and charge $`Q`$ (see Ruffini and Wheeler ). While the first two parameters $`M`$ and $`L`$ evolve on very long time scales, the charge $`Q`$ will dissipate on much shorter time scales. In particular I show that if the electromagnetic field of an EMBH is overcritical in the sense of Heisenberg and Euler, it will dissipate on an extremely short time scale of approximately $`10^7`$ sec. through a new process involving the “dyadosphere” of a black hole. In a very short time $`O(\frac{\mathrm{}}{mc^2})`$, a very large number of pairs is created there and reaches thermodynamic equilibrium with a photon gas. In the ensuing enormous PEM-pulse, a large fraction of the extractable energy of the EMBH in the sense of Christodoulou-Ruffini will be carried away. The PEM-pulse will interact with some of the baryonic matter of the uncollapsed material and the associated emission is closely related to the observed properties of GRB sources. The electromagnetic field of the remnant will further dissipate in the acceleration of cosmic rays or in the propulsion of jets on much longer time scales.

The Christodoulou-Ruffini energy-mass formula for black holes gives

$`E^2`$ $`=`$ $`M^2c^4=\left(M_{\mathrm{ir}}c^2+{\displaystyle \frac{Q^2}{2\rho _+}}\right)^2+{\displaystyle \frac{L^2c^2}{\rho _+^2}},`$ (1)   
$`S`$ $`=`$ $`4\pi \rho _+^2=4\pi (r_+^2+{\displaystyle \frac{L^2}{c^2M^2}})=16\pi \left({\displaystyle \frac{G^2}{c^4}}\right)M_{\mathrm{ir}}^2,`$ (2)   
with

$$\frac{1}{\rho _+^4}\left(\frac{G^2}{c^8}\right)\left(Q^4+\frac{L^2c^2}{4}\right)1,$$
(3)   
where $`M_{\mathrm{ir}}`$ is the irreducible mass, $`r_+`$ is the horizon radius and $`\rho _+`$ is the quasi spheroidal cylindical coordinate of the horizon evaluated at the equatorial plane. (I use c.g.s. units). From eqs. (1) and (3) it follows that up to 29$`\%`$ of the mass-energy of an extreme rotating black hole with $`L_{\mathrm{max}}=2r_+^2c^3/G`$ can be stored as rotational energy and gedanken experiments have been conceived to extract such energy (see e. g. Ruffini and Wheeler and ref.). In the case of black holes endowed with an electromagnetic field, it follows from the same equations that up to 50$`\%`$ of the mass energy of an extreme EMBH with $`Q_{\mathrm{max}}=r_+c^2/\sqrt{G}`$ can be stored in its electromagnetic field. It is appropriate to recall that even in the case of an extreme EMBH the charge to mass ratio is $`10^{18}`$ smaller then the typical charge to mass ratio found in nuclear matter, owing to the different strength and range of the nuclear and gravitational interactions. In other words it is enough to have a difference of one quantum of charge per $`10^{18}`$ nucleons in the collapsing matter for an EMBH to be extreme. By applying the classic work of Heisenberg and Euler to EMBH’s, as reformulated in a relativistic-invariant form by Schwinger, Damour and Ruffini showed that a large fraction of the energy of an EMBH can be extracted by pair creation. This energy extraction process only works for EMBH black holes with $`M_{\mathrm{ir}}<10^6M_{}`$. They also claimed that such an energy source might lead to a natural explanation for GRB’s and for ultra high energy cosmic rays.

The general considerations presented in Damour and Ruffini are correct. However, that work has an underlying assumption which only surfaces in the very last formula: that the pair created in the process of vacuum polarization is absorbed by the EMBH. That view is now fundamentally modified here by the introduction of the novel concept of the Dyadosphere of an EMBH and by the considerations which follow from its introduction.

## III The extension of the Dyadosphere

For simplicity we use the nonrotating Reissner-Nordstrom EMBH to illustrate the basic gravitational-electrodynamical process. The rotating Kerr-Newmann EMBH case will be considered elsewhere.

By introducing the dimensionless mass and charge parameters $`\mu =\frac{M}{M_{}}>3.2`$, $`\xi =\frac{Q}{Q_{\mathrm{max}}}1`$, the horizon radius may be expressed as

$$r_+=\frac{GM}{c^2}\left[1+\sqrt{1\frac{Q^2}{GM^2}}\right]=1.4710^5\mu (1+\sqrt{1\xi ^2})\mathrm{cm}.$$
(4)   
Outside the horizon the electromagnetic field measured in the orthonormal tetrad of an observer at rest at a given point $`r`$ in the usual Boyer-Lindquist coordinates (see e.g. ) has only one nonvanishing component $`𝐄=\frac{Q}{r^2}\widehat{r}`$ along the radial direction. We can evaluate the radius $`r_{\mathrm{ds}}`$ at which the electric field strength reaches the critical value $`_\mathrm{c}=\frac{m^2c^3}{\mathrm{}e}`$ introduced by Heisenberg and Euler, where $`m`$ and $`e`$ are the mass and charge of the electron. This defines the outer radius of the Dyadosphere, which extends down to the horizon and within which the electric field strength exceeds the critical value. Using the Planck charge $`q_\mathrm{c}=(\mathrm{}c)^{\frac{1}{2}}`$ and the Planck mass $`m_\mathrm{p}=(\frac{\mathrm{}c}{G})^{\frac{1}{2}}`$, I express this outer radius in the form

$$r_{\mathrm{ds}}=\left(\frac{\mathrm{}}{mc}\right)^{\frac{1}{2}}\left(\frac{GM}{c^2}\right)^{\frac{1}{2}}\left(\frac{m_\mathrm{p}}{m}\right)^{\frac{1}{2}}\left(\frac{e}{q_\mathrm{p}}\right)^{\frac{1}{2}}\left(\frac{Q}{\sqrt{G}M}\right)^{\frac{1}{2}}=1.1210^8\sqrt{\mu \xi }\mathrm{cm},$$
(5)   
which clearly shows the hybrid gravitational and quantum nature of this quantity. The radial interval $`r_+rr_{\mathrm{ds}}`$ describing the Dyadosphere as a function of the mass is illustrated in Fig. for selected values of the charge paramter $`\xi `$. It is important to note that the Dyadosphere radius is maximized for the extreme case $`\xi =1`$ and that the region exists for EMBH’s with mass larger than the upper limit for neutron stars, namely $`3.2M_{}`$ all the way up to a maximum mass of $`610^5M_{}`$. Correspondingly smaller values of the maximum mass are obtained for values of $`\xi =0.1,0.01`$ as indicated in this figure. For EMBH’s with mass larger than the maximum value stated above, the electromagnetic field (whose strength decreases inversely with the mass) never becomes critical.

## IV The inner structure of Dyadosphere

I turn now to the crucial issue of the number and energy densities of pairs created in the Dyadosphere. If we consider a shell of proper thickness $`\delta \frac{MG}{c^2}`$, the electric field is aproximately constant. We can then at each value of the radius $`r`$ model the electric field as created by a capacitor of width $`\delta `$ and surface charge density

$$\sigma (r)=\frac{Q}{4\pi r^2},$$
(6)   
If we now turn to the process of pair creation in order to apply the QED results, we assume $`\delta =\frac{\mathrm{}}{mc}`$ and we can express the rate of pair creation at a given radius $`r`$ by 

$$\frac{dN}{\sqrt{g}d^4x}=\frac{1}{4\pi c}\left(\frac{eE}{\pi \mathrm{}}\right)^2e^{\frac{\pi E_c}{E}}=\frac{1}{4\pi c}\left(\frac{4e\sigma }{\mathrm{}}\right)^2e^{\frac{\pi \sigma _c}{\sigma }},$$
(7)   
where $`E=4\pi \sigma `$, $`\sigma _c=\frac{1}{4\pi }_\mathrm{c}`$ is the critical surface charge density and $`\sqrt{g}d^4x`$ is the invariant four volume. We have for each value of the radius $`r`$, the rate of pair creation per unit proper time

$$\frac{dN}{d\tau }=\frac{1}{4\pi c}\left(\frac{4e\sigma }{\mathrm{}}\right)^2e^{\frac{\pi \sigma _c}{\sigma }}4\pi r^2\left(\frac{\mathrm{}}{mc}\right).$$
(8)   
The pair creation process will continue until a value of the field $`EE_c`$ is reached, or correspondingly, in the capacitor language, till the surface charge density reaches the critical value $`\sigma _c`$. For $`E<E_c`$ or $`\sigma <\sigma _c`$ the pair creation process is exponentially supressed. We then have

$$\sigma \sigma _c=\frac{e}{4\pi r^2}\frac{\mathrm{\Delta }N}{\mathrm{\Delta }\tau }\mathrm{\Delta }\tau .$$
(9)   
Based on eq.(8), one can get,

$$\mathrm{\Delta }\tau =\frac{\sigma \sigma _c}{\frac{e}{4\pi c}\left(\frac{4e\sigma }{\mathrm{}}\right)^2e^{\frac{\pi \sigma _c}{\sigma }}\left(\frac{\mathrm{}}{mc}\right)}<1.99\left(\frac{\mathrm{}}{mc^2\alpha }\right)=1.7610^{19}\mathrm{sec}.,$$
(10)   
where $`\alpha =\frac{e^2}{4\pi \mathrm{}c}`$ is the fine structure constant. ( Details are given in ref..) The time given by eq.(10) is so short, that the light travel time is smaller or aproximately equal to the width $`\delta `$. Under these circumstances the correlation between shells can be approximately neglected, thus we can justify the aproximation of describing the pair creation process shell by shell.

If we now turn to the Dyadosphere its extension goes from the horizon $`r_+`$ all the way to a radius $`r_{ds}`$ where $`\frac{Q}{4\pi r_{ds}}=E_c`$. The inner layer of the first shell of width $`\delta `$, consisting of charges opposite to the one of EMBH’s, will be captured by the horizon r+, leading to a new EMBH with charge $`Q_c=4\pi r_+^2E_c`$. The outer layer, of oppositely charged particles, will enter the Dyadosphere. The remaining $`(r_{ds}r_+)/\frac{\mathrm{}}{mc}`$ shells will also contribute to the plasma constituting the Dyadosphere. It is noteworty to stress tha the number of pair created is not, as nively expected $`N_{e^+e^{}}=\frac{QQ_c}{e}`$, but the number is tremendously amplified. In the limit $`r_{ds}\frac{GM}{c^2}`$, we have

$$N_{e^+e^{}}\frac{QQ_c}{e}\left[1+\frac{(r_{ds}r_+)}{\frac{\mathrm{}}{mc}}\right].$$
(11)   
## V Density and energy density of pairs in the Dyadosphere

We can now estimate from eq.(9) the number of pairs created at a given radius $`r`$ in a shell volume of proper thickness $`\delta =\frac{\mathrm{}}{mc}`$

$$\mathrm{\Delta }N=\frac{4\pi r^2}{e}(\sigma \sigma _c)=\frac{Q}{e}\left[1\left(\frac{r}{r_{ds}}\right)^2\right],$$
(12)   
and correspondingly the number density of pairs created as a function of the radial coordinate, details in ref. and 

$$n_{e^+e^{}}(r)=\frac{\mathrm{\Delta }N}{4\pi r^2\left(\frac{\mathrm{}}{mc}\right)}=\frac{Q}{e4\pi r^2\left(\frac{\mathrm{}}{mc}\right)}\left[1\left(\frac{r}{r_{ds}}\right)^2\right].$$
(13)   
In Figs. and , we plot the density of pairs for $`10M_{}`$ and $`10^5M_{}`$ EMBH for selected values of $`\xi `$. These two values of the mass were chosen to be representative of objects typical of the galactic population or for the nuclei of galaxies compatible with our upper limit of the maximum mass of $`610^5M_{}`$.

Knowing the initial electrostatic energy density, and the final electric field value $`E_\mathrm{c}(r)=\frac{Q_\mathrm{c}}{r^2}`$ holding after the evolution of the Dyadosphere, we compute the energy density of created pairs as a function of the radial coordinate by evaluating difference between the initial and final field configurations in the Dyadosphere

$$\rho _{e^+e^{}}=\frac{1}{8\pi }(E^2(r)E_c^2(r)).$$
(14)   
Their total energy is then

$$E_{e^+e^{}}^{\mathrm{tot}}=\frac{1}{2}\frac{Q^2}{r_+}(1\frac{r_+}{r_{\mathrm{ds}}})(1\left(\frac{r_+}{r_{\mathrm{ds}}}\right)^2).$$
(15)   
By taking the ratio of eq.(14) and eq.(13), we obtain the average energy per pair at each value of the radial coordinate. Fig. plots this quantity for an EMBH of $`10M_{}`$ and Fig. for an EMBH of $`10^5M_{}`$. In the first case the energy of pairs near the horizon can reach $`10`$GeV, while in the second case it never goes over a few MeV. These two values of the mass were chosen to be representative of objects typical of the galactic population or for the nuclei of galaxies compatible with our upper limit of the maximum mass of $`610^5M_{}`$.

Due to the very large pair density given by eq. and to the sizes of the cross-sections for the process $`e^+e^{}\gamma +\gamma `$, the system is expected to thermalize to a plasma configuration for which

$$N_{e^+}=N_e^{}=N_\gamma =N_{\mathrm{pair}}$$
(16)   
and reach an average temperature

$$kT_{}=\frac{E_{e^+e^{}}^{\mathrm{tot}}}{3N_{\mathrm{pair}}2.7},$$
(17)   
where $`k`$ is Boltzmann’s constant. The average energy per pair $`\frac{E_{e^+e^{}}^{\mathrm{tot}}}{N_{\mathrm{pair}}}`$ is shown as a function of the EMBH mass for selected values of the charge parameter $`\xi `$ in Fig..

## VI The energy extraction from EMBH and the Gamma-Ray bursts afterglow

Finally we can now estimate the total energy extracted by the pair creation process in EMBH’s of different masses for selected values of the charge parameter $`\xi `$ and compare and contrast these values with the maximum extractable energy given by the mass formula for Black Holes (see eqs. (1) and (3). This comparison is summarized in Fig.. The efficiency of energy extraction by pair creation sharply decreases as the maximum value of the EMBH mass for which vacuum polarization occurs is reached. In the opposite limit the efficiency approaches $`100\%`$.

As shown by Ruffini et al. the further evolution of this plasma leads to a relativistic expansion, $`e^+e^{}`$ annihilation and an enormous pair-electromagnetic-pulse “PEM-pulse”.

The evolution of this PEM-pulse has been studied by introducing a variety of models based on relativistic hydrodynamical equations, on quasi analitical and numerical treatments and comparing and constracting the results with the ones obtained by numerical simulation in the super-computer of Livermore in California . It was possible to follow the evolution of the PEM-pulse and its reaching relativistic velocities. A variety of models has been considered both in the absence and presence of baryonic matter, three examples are given in Fig.8.

The PEM-pulse in its evolution leads to a variety of emmision profiles, as a function of the EMBH mass, charge and of the amount of baryonic matter surrounding EMBH’s. A brief example is given in Fig.9 where most explicitly it is shown how the observed power law the afterglow of the Gamma-rays burst is recovered with their characteristic time scale of $`10^7`$ sec. Details for the computations and other example are given in ref..

## VII Conclusions

If this basic scenario is confirmed by observations, other fundamental questions should be investigated to understand the origin of the dyadosphere. Some preliminary work along these lines has already been done by Wilson  who has shown how relativistic magneto-hydrodynamical processes occurring in the accreting material around an already formed black hole lead naturally to very effective charge separation and to reaching the critical value of the charge given by the limiting value of Eq. (3) on a timescale of $`10^2`$$`10^3GM/c^3`$. These studies show the very clear tendency of powerful processes of charge separation and of magnetosphere formation to occur in accretion processes. They cannot, however, be simply extrapolated to the formation of the dyadosphere. The time scale of the dyadosphere discharge is of the order of $`10^{19}`$ sec (for a detailed discussion including relativistic effects, see Jantzen and Ruffini, ). Such a time scale is much shorter than the characteristic magnetohydrodynamical time scales. We expect that the formation of the dyadosphere should only occur during the gravitational collapse itself and in the process of formation of the EMBH, with the formation of a charge depleted region with an electric field sufficient to polarize the vacuum. A complementary aspect dealing with the extremization and equipartition of the electromagnetic energy in a gravitating rotating object advanced in Ruffini and Treves  should also now be considered again. Modeling such a problem seems to be extremely difficult at the present time, also in absence of detailed observations narrowing down the values of the fundamental parameters involved. Confirmation of the basic predictions of the dyadosphere model by observations of gamma ray bursts and their afterglow would provide motivation and essential information to attack such a difficult problem.

The model of GRB’s, based on the Dyadosphere, is very different from the ones debated in the last ten years in the scientific literature, it can overcome some of their basic difficulties and on some specific aspects it can also have important analogies:

* the EMBH drastically differs from the most popular binary neutron stars model  debated in the last 16 years. It presents over those models the distinct advantage of not having the recognized difficulty of explaining the energetic of GRB’s , and especially offers the possibility of computing explicitly and uniquely the details of the initial conditions by the energetic of the Dyadosphere as given in the above Fig.4 and Fig.5.
* the further evolution of the Dyadosphere (the PEM-pulse), has been analyzed by a variety of theoretical models based on relativistic hydrodynamics equations both on semi- analytical-numerical idealized models as well as with the help of a fully general relativistic hydrodynamical code. The combined results are at variance with many of the considerations in the current published papers dealing with fireballs and relativistic winds , though for a few aspects they present striking analogies. Details are given in ref..
* the new treatment lead uniquely to the computation of the further evolution of the PEM Pulse in the afterglow era and again these specific analysis which can be uniquely derived from the above assumptions are at variance with the ones in the existing literature  though they present some interesting analogies and they appear to be consistent with the latest observations .

It goes without saying that the Heisenberg-Euler process of vacuum polarization considered here and in Damour and Ruffini  has nothing to do with the evaporation of black holes considered by Hawking  either from a qualitative or quantitative point of view. The effective temperature of black hole evaporation given by Hawking for a $`10M_{}`$ black hole is $`T6.210^9K^{}`$, which implies a black hole lifetime of $`\tau 10^{66}`$ years and an energy flux of $`10^{24}`$ergs/sec, while a $`10^5M_{}`$ black hole has a Hawking temperature of $`T6.210^{13}K^{}`$ and lifetime of $`\tau 10^{78}`$ years with an energy flux of $`10^{32}`$ergs/sec.

From a fundamental point of view, the process occuring in the Dyadosphere represents the first mechanism capable of extracting large amounts of energy from a Black Hole with an extremely high efficiency.