Growth of Accreting Supermassive Black Hole Seeds and Neutrino Radiation

In the framework of microscopic theory of black hole (MTBH), which explores the most important processes of rearrangement of vacuum state and spontaneous breaking of gravitation gauge symmetry at huge energies, we have undertaken a large series of numerical simulations with the goal to trace an evolution of themass assembly history of 377 plausible accreting supermassive black hole seeds in active galactic nuclei (AGNs) to the present time and examine the observable signatures today. Given the redshifts, masses, and luminosities of these black holes at present time collected from the literature, we compute the initial redshifts and masses of the corresponding seed black holes. For the present massesMBH/M⊙ ≃ 1.1 × 10 6 to 1.3 × 1010 of 377 black holes, the computed intermediate seedmasses are ranging fromMSeed BH /M⊙ ≃ 26.4 to 2.9×10 . We also compute the fluxes of ultrahigh energy (UHE) neutrinos produced via simple or modified URCA processes in superdense protomatter nuclei. The AGNs are favored as promising pure UHE neutrino sources, because the computed neutrino fluxes are highly beamed along the plane of accretion disk, peaked at high energies, and collimated in smaller opening angle (θ ≪ 1).


Introduction
With typical bolometric luminosities ∼10 45−48 erg s −1 , the AGNs are amongst the most luminous emitters in the universe, particularly at high energies (gamma-rays) and radio wavelengths. From its historical development, up to current interests, the efforts in the AGN physics have evoked the study of a major unsolved problem of how efficiently such huge energies observed can be generated. This energy scale severely challenges conventional source models. The huge energy release from compact regions of AGN requires extremely high efficiency (typically ≥10 per cent) of conversion of rest mass to other forms of energy. This serves as the main argument in favour of supermassive black holes, with masses of millions to billions of times the mass of the Sun, as central engines of massive AGNs. The astrophysical black holes come in a wide range of masses, from ≥3 ⊙ for stellar mass black holes [1] to ∼10 10 ⊙ for supermassive black holes [2,3]. Demography of local galaxies suggests that most galaxies harbour quiescent supermassive black holes in their nuclei at the present time and that the mass of the hosted black hole is correlated with properties of the host bulge. The visible universe should therefore contain at least 100 billion supermassive black holes. A complex study of evolution of AGNs requires an answer to the key questions such as how did the first black holes form, how did massive black holes get to the galaxy centers, and how did they grow in accreting mass, namely, an understanding of the important phenomenon of mass assembly history of accreting supermassive black hole seeds. The observations support the idea that black holes grow in tandem with their hosts throughout cosmic history, starting from the earliest times. While the exact mechanism for the formation of the first black holes is not currently known, there are several prevailing theories [4]. However, each proposal towards formation and growth of initial seed black holes has its own advantage and limitations in proving the whole view of the issue. In this report we review the mass assembly history of 377 plausible accreting supermassive black hole seeds in AGNs and their neutrino radiation in the framework of gravitation theory, which explores the most important processes of rearrangement of vacuum state and a spontaneous breaking of gravitation gauge symmetry at huge energies. We will proceed according to the following structure. Most observational, theoretical, and computational aspects of the growth of black hole seeds are summarized in Section 2. The other important phenomenon of ultrahigh energy cosmic rays, in relevance to AGNs, is discussed in Section 3. The objectives of suggested approach are outlined in Section 4. In Section 5 we review the spherical accretion on superdense protomatter nuclei, in use. In Section 6 we discuss the growth of the seed black hole at accretion and derive its intermediate mass, initial redshift, and neutrino preradiation time (PRT). Section 7 is devoted to the neutrino radiation produced in superdense protomatter nuclei. The simulation results of the seed black hole intermediate masses, PRTs, seed redshifts, and neutrino fluxes for 377 AGN black holes are brought in Section 8. The concluding remarks are given in Section 9. We will refrain from providing lengthy details of the proposed gravitation theory at huge energies and neutrino flux computations. For these the reader is invited to visit the original papers and appendices of the present paper. In the latter we also complete the spacetime deformation theory, in the model context of gravitation, by new investigation of building up the complex of distortion (DC) of spacetime continuum and showing how it restores the world-deformation tensor, which still has been put in by hand. Finally, note that we regard the considered black holes only as the potential neutrino sources. The obtained results, however, may suffer if not all live black holes at present reside in final stage of their growth driven by the formation of protomatter disk at accretion and they radiate neutrino. We often suppress the indices without notice. Unless otherwise stated, we take geometrized units throughout this paper.

A Breakthrough in Observational and Computational Aspects on Growth of Black Hole Seeds
Significant progress has been made in the last few years in understanding how supermassive black holes form and grow. Given the current masses of 10 6−9 ⊙ , most black hole growth happens in the AGN phase. A significant fraction of the total black hole growth, 60% [6], happens in the most luminous AGN, quasars. In an AGN phase, which lasts ∼10 8 years, the central supermassive black hole can gain up to ∼10 7−8 ⊙ , so even the most massive galaxies will have only a few of these events over their lifetime. Aforesaid gathers support especially from a breakthrough made in recent observational, theoretical, and computational efforts in understanding of evolution of black holes and their host galaxies, particularly through self-regulated growth and feedback from accretionpowered outflows; see, for example, [4,[7][8][9][10][11][12][13][14][15][16][17][18]. Whereas the multiwavelength methods are used to trace the growth of seed BHs, the prospects for future observations are reviewed. The observations provide strong support for the existence of a correlation between supermassive black holes and their hosts out to the highest redshifts. The observations of the quasar luminosity function show that the most supermassive black holes get most of their mass at high redshift, while at low redshift only low mass black holes are still growing [19]. This is observed in both optical [20] and hard Xray luminosity functions [19,21], which indicates that this result is independent of obscuration. Natarajan [13] has reported that the initial black hole seeds form at extremely high redshifts from the direct collapse of pregalactic gas discs. Populating dark matter halos with seeds formed in this fashion and using a Monte-Carlo merger tree approach, he has predicted the black hole mass function at high redshifts and at the present time. The most aspects of the models that describe the growth and accretion history of supermassive black holes and evolution of this scenario have been presented in detail by [9,10]. In these models, at early times the properties of the assembling black hole seeds are more tightly coupled to properties of the dark matter halo as their growth is driven by the merger history of halos. While a clear picture of the history of black hole growth is emerging, significant uncertainties still remain [14], and in spite of recent advances [6,13], the origin of the seed black holes remains an unsolved problem at present. The NuSTAR deep high-energy observations will enable obtaining a nearly complete AGN survey, including heavily obscured Comptonthick sources, up to ∼ 1.5 [22]. A similar mission, ASTRO-H [23], will be launched by Japan in 2014. These observations in combination with observations at longer wavelengths will allow for the detection and identification of most growing supermassive black holes at ∼ 1. The ultradeep X-ray and near-infrared surveys covering at least ∼1 deg 2 are required to constrain the formation of the first black hole seeds. This will likely require the use of the next generation of spacebased observatories such as the James Webb Space Telescope and the International X-ray Observatory. The superb spatial resolution and sensitivity of the Atacama Large Millimeter Array (ALMA) [24] will revolutionize our understanding of galaxy evolution. Combining these new data with existing multiwavelength information will finally allow astrophysicists to pave the way for later efforts by pioneering some of the census of supermassive black hole growth, in use today.

UHE Cosmic-Ray Particles
The galactic sources like supernova remnants (SNRs) or microquasars are thought to accelerate particles at least up to energies of 3 × 10 15 eV. The ultrahigh energy cosmicray (UHECR) particles with even higher energies have since been detected (comprehensive reviews can be found in [25][26][27][28][29]). The accelerated protons or heavier nuclei up to energies exceeding 10 20 eV are firstly observed by [30]. The cosmicray events with the highest energies so far detected have energies of 2 × 10 11 GeV [31] and 3 × 10 11 GeV [32]. These energies are 10 7 times higher than the most energetic manmade accelerator, the LHC at CERN. These highest energies are believed to be reached in extragalactic sources like AGNs or gamma-ray bursts (GRBs). During propagation of such energetic particles through the universe, the threshold for pion photoproduction on the microwave background is ∼2 × 10 10 GeV, and at ∼3 × 10 11 GeV the energy-loss distance is about 20 Mpc. Propagation of cosmic rays over substantially larger distances gives rise to a cutoff in the spectrum at ∼ 10 11 GeV as was first shown by [33,34], the GZK cutoff. The recent confirmation [35,36] of GZK suppression in the Journal of Astrophysics 3 cosmic-ray energy spectrum indicates that the cosmic rays with energies above the GZK cutoff, GZK ∼ 40 EeV, mostly come from relatively close (within the GZK radius, GZK ∼ 100 Mpc) extragalactic sources. However, despite the detailed measurements of the cosmic-ray spectrum, the identification of the sources of the cosmic-ray particles is still an open question as they are deflected in the galactic and extragalactic magnetic fields and hence have lost all information about their origin when reaching Earth. Only at the highest energies beyond ∼10 19.6 GeV cosmic-ray particles may retain enough directional information to locate their sources. The latter must be powerful enough to sustain the energy density in extragalactic cosmic rays of about 3 × 10 −19 erg cm −3 which is equivalent to ∼8 × 10 44 erg Mpc −3 yr −1 . Though it has not been possible up to now to identify the sources of galactic or extragalactic cosmic rays, general considerations allow limiting potential source classes. For example, the existing data on the cosmic-ray spectrum and on the isotropic 100 MeV gammaray background limit significantly the parameter space in which topological defects can generate the flux of the highest energy cosmic rays and rule out models with the standard X-particle mass of 10 16 GeV and higher [37]. Eventually, the neutrinos will serve as unique astronomical messengers, and they will significantly enhance and extend our knowledge on galactic and extragalactic sources of the UHE universe. Indeed, except for oscillations induced by transit in a vacuum Higgs field, neutrinos can penetrate cosmological distances and their trajectories are not deflected by magnetic fields as they are neutral, providing powerful probes of high energy astrophysics in ways which no other particle can. Moreover, the flavor composition of neutrinos originating at astrophysical sources can serve as a probe of new physics in the electroweak sector. Therefore, an appealing possibility among the various hypotheses of the origin of UHECR is so-called Z-burst scenario [38][39][40][41][42][43][44][45][46][47][48][49][50][51]. This suggests that if ZeV astrophysical neutrino beam is sufficiently strong, it can produce a large fraction of observed UHECR particles within 100 Mpc by hitting local light relic neutrinos clustered in dark halos and form UHECR through the hadronic Z (s-channel production) and W-bosons (t-channel production) decays by weak interactions. The discovery of UHE neutrino sources would also clarify the production mechanism of the GeV-TeV gamma rays observed on Earth [43,52,53] as TeV photons are also produced in the up-scattering of photons in reactions to accelerated electrons (inverse-Compton scattering). The direct link between TeV gamma-ray photons and neutrinos through the charged and neutral pion production, which is well known from particle physics, allows for a quite robust prediction of the expected neutrino fluxes provided that the sources are transparent and the observed gamma rays originate from pion decay. The weakest link in the Z-burst hypothesis is probably both unknown boosting mechanism of the primary neutrinos up to huge energies of hundreds ZeV and their large flux required at the resonant energy ] ≃ 2 /(2 ] ) ≃ 4.2 × 10 21 eV (eV/ ] ) well above the GZK cutoff. Such a flux severely challenges conventional source models. Any concomitant photon flux should not violate existing upper limits [37,48,49,54]. The obvious question is then raised: where in the Cosmos are these neutrinos coming from? It turns out that currently, at energies in excess of 10 19 eV, there are only two good candidate source classes for UHE neutrinos: AGNs and GRBs. The AGNs as significant point sources of neutrinos were analyzed in [50,55,56]. While hard to detect, neutrinos have the advantage of representing aforesaid unique fingerprints of hadron interactions and, therefore, of the sources of cosmic rays. Two basic event topologies can be distinguished: track-like patterns of detected Cherenkov light (hits) which originate from muons produced in charged-current interactions of muon neutrinos (muon channel); spherical hit patterns which originate from the hadronic cascade at the vertex of neutrino interactions or the electromagnetic cascade of electrons from chargedcurrent interactions of electron neutrinos (cascade channel).
If the charged-current interaction happens inside the detector or in case of charged-current tau-neutrino interactions, these two topologies overlap which complicates the reconstruction. At the relevant energies, the neutrino is approximately collinear with the muon and, hence, the muon channel is the prime channel for the search for point-like sources of cosmic neutrinos. On the other hand, cascades deposit all of their energy inside the detector and therefore allow for a much better energy reconstruction with a resolution of a few 10%. Finally, numerous reports are available at present in literature on expected discovery potential and sensitivity of experiments to neutrino point-like sources. Currently operating high energy neutrino telescopes attempt to detect UHE neutrinos, such as ANTARES [57,58] which is the most sensitive neutrino telescope in the Northern Hemisphere, Ice-Cube [35,[59][60][61][62][63][64] which is worldwide largest and hence most sensitive neutrino telescope in the Southern Hemisphere, BAIKAL [65], as well as the CR extended experiments of The Telescope Array [66], Pierre Auger Observatory [67,68], and JEM-EUSO mission [69]. The JEM-EUSO mission, which is planned to be launched by a H2B rocket around 2015-2016, is designed to explore the extremes in the universe and fundamental physics through the detection of the extreme energy ( > 10 20 eV) cosmic rays. The possible origins of the soon-to-be famous 28 IceCube neutrino-PeV events [59][60][61] are the first hint for astrophysical neutrino signal. Aartsen et al. have published an observation of two ∼1 PeV neutrinos, with a value 2.8 beyond the hypothesis that these events were atmospherically generated [59]. The analysis revealed an additional 26 neutrino candidates depositing "electromagnetic equivalent energies" ranging from about 30 TeV up to 250 TeV [61]. New results were presented at the IceCube Particle Astrophysics Symposium (IPA 2013) [62][63][64]. If cosmic neutrinos are primarily of extragalactic origin, then the 100 GeV gamma ray flux observed by Fermi-LAT constrains the normalization at PeV energies at injection, which in turn demands a neutrino spectral index Γ < 2.1 [70].

MTBH, Revisited: Preliminaries
For the benefit of the reader, a brief outline of the key ideas behind the microscopic theory of black hole, as a guiding principle, is given in this section to make the rest of the paper understandable. There is a general belief reinforced by statements in textbooks that, according to general relativity (GR), a long-standing standard phenomenological black hole model (PBHM)-namely, the most general Kerr-Newman black hole model, with parameters of mass ( ), angular momentum ( ), and charge ( ), still has to be put in by hand-can describe the growth of accreting supermassive black hole seed. However, such beliefs are suspect and should be critically reexamined. The PBHM cannot be currently accepted as convincing model for addressing the aforementioned problems, because in this framework the very source of gravitational field of the black hole is a kind of curvature singularity at the center of the stationary black hole. A meaningless central singularity develops which is hidden behind the event horizon. The theory breaks down inside the event horizon which is causally disconnected from the exterior world. Either the Kruskal continuation of the Schwarzschild ( = 0, = 0) metric, or the Kerr ( = 0) metric, or the Reissner-Nordstrom ( = 0) metric, shows that the static observers fail to exist inside the horizon. Any object that collapses to form a black hole will go on to collapse to a singularity inside the black hole. Thereby any timelike worldline must strike the central singularity which wholly absorbs the infalling matter. Therefore, the ultimate fate of collapsing matter once it has crossed the black hole surface is unknown. This, in turn, disables any accumulation of matter in the central part and, thus, neither the growth of black holes nor the increase of their massenergy density could occur at accretion of outside matter, or by means of merger processes. As a consequence, the mass and angular momentum of black holes will not change over the lifetime of the universe. But how can one be sure that some hitherto unknown source of pressure does not become important at huge energies and halt the collapse? To fill the void which the standard PBHM presents, one plausible idea to innovate the solution to alluded key problems would appear to be the framework of microscopic theory of black hole. This theory has been originally proposed by [71] and references therein and thoroughly discussed in [72][73][74][75].
Here we recount some of the highlights of the MTBH, which is the extension of PBHM and rather completes it by exploring the most important processes of rearrangement of vacuum state and a spontaneous breaking of gravitation gauge symmetry at huge energies [71,74,76]. We will not be concerned with the actual details of this framework but only use it as a backdrop to validate the theory with some observational tests. For details, the interested reader is invited to consult the original papers. Discussed gravitational theory is consistent with GR up to the limit of neutron stars. But this theory manifests its virtues applied to the physics of internal structure of galactic nuclei. In the latter a significant change of properties of spacetime continuum, so-called inner distortion (ID), arises simultaneously with the strong gravity at huge energies (see Appendix A). Consequently the matter undergoes phase transition of second kind, which supplies a powerful pathway to form a stable superdense protomatter core (SPC) inside the event horizon. Due to this, the stable equilibrium holds in outward layers too and, thus, an accumulation of matter is allowed now around the SPC.
The black hole models presented in phenomenological and microscopic frameworks have been schematically plotted in Figure 1, to guide the eye. A crucial point of the MTBH is that a central singularity cannot occur, which is now replaced by finite though unbelievably extreme conditions held in the SPC, where the static observers existed. The SPC surrounded by the accretion disk presents the microscopic model of AGN. The SPC accommodates the highest energy scale up to hundreds of ZeV in central protomatter core which accounts for the spectral distribution of the resulting radiation of galactic nuclei. External physics of accretion onto the black hole in earlier part of its lifetime is identical to the processes in Schwarzschild's model. However, a strong difference in the model context between the phenomenological black hole and the SPC is arising in the second part of its lifetime (see Section 6). The seed black hole might grow up driven by the accretion of outside matter when it was getting most of its mass. An infalling matter with time forms the protomatter disk around the protomatter core tapering off faster at reaching out the thin edge of the event horizon. At this, metric singularity inevitably disappears (see appendices) and the neutrinos may escape through vista to outside world, even after the neutrino trapping. We study the growth of protomatter disk and derive the intermediate mass and initial redshift of seed black hole and examine luminosities, neutrino surfaces for the disk. In this framework, we have computed the fluxes of UHE neutrinos [75], produced in the medium of the SPC via simple (quark and pionic reactions) or modified URCA processes, even after the neutrino trapping (G. Gamow was inspired to name the process URCA after the name of a casino in Rio de Janeiro, when M. Schenberg remarked to him that "the energy disappears in the nucleus of the supernova as quickly as the money disappeared at that roulette table"). The "trapping" is due to the fact that as the neutrinos are formed in protomatter core at superhigh densities they experience greater difficulty escaping from the protomatter core before being dragged along with the matter; namely, the neutrinos are "trapped" comove with matter. The part of neutrinos annihilates to produce, further, the secondary particles of expected ultrahigh energies. In this model, of course, a key open question is to enlighten the mechanisms that trigger the activity, and how a large amount of matter can be steadily funneled to the central regions to fuel this activity. In high luminosity AGNs the large-scale internal gravitational instabilities drive gas towards the nucleus which trigger big starbursts, and the coeval compact cluster just formed. It seemed they have some connection to the nuclear fueling through mass loss of young stars as well as their tidal disruption and supernovae. Note that we regard the UHECR particles as a signature of existence of superdence protomatter sources in the universe. Since neutrino events are expected to be of sufficient intensity, our estimates can be used to guide investigations of neutrino detectors for the distant future.

Spherical Accretion onto SPC
As alluded to above, the MTBH framework supports the idea of accreting supermassive black holes which link to AGNs. In order to compute the mass accretion ratė, in use, it is necessary to study the accretion onto central supermassive SPC. The main features of spherical accretion can be briefly summed up in the following three idealized models that illustrate some of the associated physics [72].

Free Radial Infall.
We examine the motion of freely moving test particle by exploring the external geometry of the SPC, with the line element (A.7), at = 0. Let us denote the 4vector of velocity of test particle V =̃/̃,̃= (̃,̃,̃,̃), and consider it initially for simplest radial infall V 2 = V 3 = 0. We determine the value of local velocity V̂< 0 of the particle for the moment of crossing the EH sphere, as well as at reaching the surface of central stable SPC. The equation of geodesics is derived from the variational principle ∫ = 0, which is the extremum of the distance along the wordline for the Lagrangian at hand wherė≡̃/ is the -component of 4-momentum and is the affine parameter along the worldline. We are using an affine parametrization (by a rescaling → ( )) such that = const is constant along the curve. A static observer makes measurements with local orthonormal tetrad: The Euler-Lagrange equations for̃,̃, and̃can be derived from the variational principle. A local measurement of the particle's energy made by a static observer in the equatorial plane gives the time component of the 4-momentum as measured in the observer's local orthonormal frame. This is the projection of the 4-momentum along the time basis vector. The Euler-Lagrange equations show that if we orient the coordinate system as initially the particle is moving in the equatorial plane (i.e.,̃= /2,̇= 0), then the particle always remains in this plane. There are two constants of the motion corresponding to the ignorable coordinates̃and̃, namely, the -"energy-at-infinity" and the -angular momentum. We conclude that the free radial infall of a particle from the infinity up to the moment of crossing the EH sphere, as well as at reaching the surface of central body, is absolutely the same as in the Schwarzschild geometry of black hole ( Figure 2(a)). We clear up a general picture of orbits just outside the event horizon by considering the Euler-Lagrange equation for radial component with "effective potential. " The circular orbits are stable if is concave up, namely, at̃> 4̃, wherẽis the mass of SPC. The binding energy per unit mass of a particle in the last stable circular orbit at̃= 4ĩ s̃b ind = ( − )/̃≃ 1 − (27/32) 1/2 . Namely, this is the fraction of rest-mass energy released when test particle originally at rest at infinity spirals slowly toward the SPC to the innermost stable circular orbit and then plunges into it. Thereby one of the important parameters is the capture cross section for particles falling in from infinity: capt = 2 max , where max is the maximum impact parameter of a particle that is captured.

Collisionless Accretion.
The distribution function for a collisionless gas is determined by the collisionless Boltzmann equation or Vlasov equation. For the stationary and spherical flow we obtain theṅ  Figure 2: (a) The free radial infall of a particle from the infinity to EH sphere ( 0 = 1), which is similar to the Schwarzschild geometry of BH. Crossing the EH sphere, a particle continues infall reaching finally the surface ( 0 = 2) of the stable SPC. (b) Approaching the EH sphere ( 0 = 1), the particle concentration increases asymptotically until the threshold value of protomatter. Then, due to the action of cutoff effect, the metric singularity vanishes and the particles well pass EH sphere.
where the particle density ∞ is assumed to be uniform at far from the SPC and the particle speed is V ∞ ≪ 1. During the accretion process the particles approaching the EH become relativistic. Approaching event horizon, the particle concentration increases asymptotically as (Figure 2(b)). Due to the action of cutoff effect, the metric singularity then vanishes and the particles well pass EH sphere ( 0 = 1) and in the sequel form the protomatter disk around the protomatter core.

Hydrodynamic Accretion.
For real dynamical conditions found in considered superdense medium, it is expected that the mean free path for collisions will be much shorter than the characteristic length scale; that is, the accretion of ambient gas onto a stationary, nonrotating compact SPC will be hydrodynamical in nature. For any equation of state obeying the causality constraint the sound speed implies 2 < 1 and the flow must pass through a critical sonic point outside the event horizon. The locally measured particle velocity reads V̂= (1 −̃0 0 / 2 ), where = ∞ / = (̃0 0 /(1 − 2 )) 1/2 and ∞ is the energy at infinity of individual particle of the mass . Thus, the proper flow velocity V̂= → 0 and is subsonic. At̃= /2, the proper velocity equals the speed of light |V̂| = = 1 > and the flow is supersonic. This condition is independent of the magnitude of and is not sufficient by itself to guarantee that the flow passes through a critical point outside EH. For largẽ≥ , it is expected that the particles be nonrelativistic with ≤ ≪ 1 (i.e., ≪ 2 / = 10 13 ), as they were nonrelativistic at infinity ( ∞ ≪ 1). Considering the equation of accretion onto superdense protomatter core, which is an analogue of Bondi equations for spherical, steadystate adiabatic accretion onto the SPC, we determine a mass accretion ratė= where prime ( ) denotes differentiation with respect tõat the point . The gas compression can be estimated as The approximate equality between the sound speed and the mean particle speed implies that the hydrodynamic accretion rate is larger than the collisionless accretion rate by the large factor ≈10 9 .

The Intermediate Mass, PRT, and Initial Redshift of Seed Black Hole
The key objectives of the MTBH framework are then an increase of the mass, Seed BH , gravitational radius, Seed , and of the seed black hole, BH Seed , at accretion of outside matter. Thereby an infalling matter forms protomatter disk around protomatter core tapering off faster at reaching the thin edge of event horizon. So, a practical measure of growth BH Seed → BH may most usefully be the increase of gravitational radius or mass of black hole: where , , and , respectively, are the total mass, density, and volume of protomatter disk. At the valuêB H of gravitational radius, when protomatter disk has finally reached the event horizon of grown-up supermassive black Journal of Astrophysics and in approximation ≪̂B H we set 0 ( 0 ) ≃ 0 ≃ / √ 2.

PRT.
The PRT is referred to as a lapse of time BH from the birth of black hole till neutrino radiation, the earlier part of the lifetime. That is, BH = /̇, wherėis the accretion rate. In approximation at hand ≪ , the PRT reads In case of collisionless accretion, (3) and (11) give BH ≃ 2.6 ⋅ 10 16 cm In case of hydrodynamic accretion, (4) and (11) yield BH ≃ 8.8 ⋅ 10 38 Note that the spherical accretion onto black hole, in general, is not necessarily an efficient mechanism for converting restmass energy into radiation. Accretion onto black hole may be far from spherical accretion, because the accreted gas possesses angular momentum. In this case, the gas will be thrown into circular orbits about the black hole when centrifugal forces will become significant before the gas plunges through the event horizon. Assuming a typical mass-energy conversion efficiency of about ∼ 10%, in approximation ≪ , according to (12) and (13), the resulting relationship of typical PRT versus bolometric luminosity becomes We supplement this by computing neutrino fluxes in the next section.

Redshift of Seed Black Hole.
Interpreting the redshift as a cosmological Doppler effect and that the Hubble law could most easily be understood in terms of expansion of the universe, we are interested in the purely academic question of principle to ask what could be the initial redshift, Seed , of seed black hole if the mass, the luminosity, and the redshift, , of black hole at present time are known. To follow the history of seed black hole to the present time, let us place ourselves at the origin of coordinates = 0 (according to the Cosmological Principle, this is mere convention) and consider a light traveling to us along the − direction, with angular variables fixed. If the light has left a seed black hole, located at , , and , at time , and it has to reach us at a time 0 , then a power series for the redshift as a function of the time of flight is Seed = 0 ( 0 − ) + ⋅ ⋅ ⋅ , where 0 is the present moment and 0 is Hubble's constant. Similar expression, = 0 ( 0 − 1 )+⋅ ⋅ ⋅ , can be written for the current black hole, located at 1 , 1 , and 1 , at time 1 , where 1 = + BH , as seed black hole is an object at early times. Hence, in the first-order approximation by Hubble's constant, we may obtain the following relation between the redshifts of seed and present black holes: Seed ≃ + 0 BH . This relation is in agreement with the scenario of a general recession of distant galaxies away from us in all directions, the furthest naturally being those moving the fastest. This relation, incorporating with (14), for the value 0 = 70 [km]/[s Mpc] favored today yields Seed ≃ + 2.292 × 10 28

UHE Neutrino Fluxes
The Here ( ) = (1 + ) ∫ where = /2 and is the thickness of the protomatter disk at the edge of even horizon. The resulting total UHE neutrino flux of cooling of the SPC can be obtained as URCA ] where the neutrino is radiated in a cone with the beaming angle ∼ ≪ 1, ( ) = (1 + ) ∫ As it is seen, the nucleon modified URCA reactions can contribute efficiently only to extragalactic objects with enough small redshift ≪ 1.

Pionic Reactions.
The pionic reactions, occurring in the superdense protomatter medium of SPC, allow both the distorted energy and momentum to be conserved. This is the analogue of the simple URCA processes: and the two inverse processes. As in the modified URCA reactions, the total rate for all four processes is essentially four times the rate of each reaction alone. The muons are already present when pions appear. The neutrino luminosity of the SPC of given mass,̃, by pionic reactions reads [75] Then, the UHE neutrino total flux is The resulting total energy-loss rate will then be dramatically larger due to the pionic reactions (19) rather than the modified URCA processes.

Quark Reactions.
In the superdense protomatter medium the distorted quark Fermi energies are far below the charmed c-, t-, and b-quark production thresholds. Therefore, only up-, down-, and strange quarks are present. The equilibrium is maintained by reactions like which are decay and its inverse. These reactions constitute simple URCA processes, in which there is a net loss of a ] ] pair at nonzero temperatures. In this application a sufficient accuracy is obtained by assuming -equilibrium and that the neutrinos are not retained in the medium of Λ-like protomatter. The quark reactions (22) and (23) proceed at equal rates in equilibrium, where the participating quarks must reside close to their Fermi surface. Hence, the total energy of flux due to simple URCA processes is rather twice than that of (22) or (23) alone. For example, the spectral fluxes of the UHE antineutrinos and neutrinos for different redshifts from quark reactions are plotted, respectively, in Figures 4 and 5 [75]. The total flux of UHE neutrino can be written as

Simulation
For simulation we use the data of AGN/BH mass and luminosity estimates for 377 black holes presented by [5]. These masses are mostly based on the virial assumption for the broad emission lines, with the broad-line region size determined from either reverberation mapping or optical luminosity. Additional black hole mass estimates based on properties of the host galaxy bulges, either using the observed stellar velocity dispersion or using the fundamental plane relation. Since the aim is to have more than a thousand of realizations, each individual run is simplified, with a use of previous algorithm of the SPC-configurations [71] as a working model, given in Appendix G. Computing the corresponding PRTs, seed black hole intermediate masses, and total neutrino fluxes, a main idea comes to solving an inverse problem. Namely, by the numerous reiterating integrations of the state equations of SPC-configurations we determine those required central values of particle concentratioñ(0) and IDfield (0), for which the integrated total mass of configuration has to be equal to the black hole mass BH given from observations. Along with all integral characteristics, the radius is also computed, which is further used in (10), (14), (15), (18), (21), and (24) for calculating Seed BH , BH , Seed , and ] , respectively. The results are summed up in Tables 1, 2, 3, 4 and 5. Figure 6 gives the intermediate seed masses Seed BH / ⊙ versus the present masses BH / ⊙ of 337 black holes, on logarithmic scales. For the present masses BH / ⊙ ≃ 1.1 × 10 6 to 1.3 × 10 10 , the computed intermediate seed masses  (1) In accordance, the AGNs are favored as promising pure neutrino sources because the computed neutrino fluxes are highly beamed along the plane of accretion disk and peaked at high energies and collimated in smaller opening angle ∼ = /2 ≪ 1. To render our discussion here a bit more transparent and to obtain some feeling for the parameter we may estimate ≃ 1.69×10 −10 , just for example, only, for the suppermassive black hole of typical mass ∼10 9 ⊙ (2 = 5.9 × 10 14 cm), and so ∼ 1 km. But the key problem of fixing the parameter more accurately from experiment would be an important topic for another investigation elsewhere.

Conclusions
The growth of accreting supermassive black hole seeds and their neutrino radiation are found to be common phenomena in the AGNs. In this report, we further expose the assertions made in the framework of microscopic theory of black hole via reviewing the mass assembly history of 377 plausible accreting supermassive black hole seeds. After the numerous reiterating integrations of the state equations of SPC-configurations, we compute their intermediate seed masses, Seed BH , PRTs, initial redshifts, Seed , and neutrino fluxes. All the results are presented in Tables 1-5. Figure 6 gives the intermediate seed masses Seed BH / ⊙ versus the present masses BH / ⊙ of 337 black holes, on logarithmic scales. In accordance, the AGNs are favored as promising pure UHE neutrino sources. Such neutrinos may reveal clues on the puzzle of origin of UHE cosmic rays. We regard the considered black holes only as the potential neutrino sources. The obtained results, however, may suffer and that would be underestimated if not all 377 live black holes in the BH / ⊙ ≃ 1.1 × 10 6 to 1.3 × 10 10 mass range at present reside in final stage of their growth, when the protomatter disk driven by accretion has reached the event horizon.

A. Outline of the Key Points of Proposed Gravitation Theory at Huge Energies
The proposed gravitation theory explores the most important processes of rearrangement of vacuum state and a spontaneous breaking of gravitation gauge symmetry at huge energies. From its historical development, the efforts in gauge treatment of gravity mainly focus on the quantum gravity and microphysics, with the recent interest, for example, in the theory of the quantum superstring or, in the very early universe, in the inflationary model. The papers on the gauge treatment of gravity provide a unified picture of gravity modified models based on several Lie groups. However, currently no single theory has been uniquely accepted as the convincing gauge theory of gravitation which could lead to a consistent quantum theory of gravity. They have evoked the possibility that the treatment of spacetime might involve non-Riemannian features on the scale of the Planck length. This necessitates the study of dynamical theories involving post-Riemannian geometries. It is well known that the notions of space and connections should be separated; see, for example, [78][79][80][81]. The curvature and torsion are in fact properties of a connection, and many different connections are allowed to exist in the same spacetime. Therefore, when considering several connections with different curvature and torsion, one takes spacetime simply as a manifold and connections as additional structures. From this view point in a recent paper [82] we tackle the problem of spacetime deformation. This theory generalizes and, in particular cases, fully recovers the results of the conventional theory. Conceptually and techniquewise this theory is versatile and powerful and manifests its practical and technical virtue in the fact that through a nontrivial choice of explicit form of the world-deformation tensor, which we have at our disposal, in general, we have a way to deform the spacetime which displayed different connections, which may reveal different post-Riemannian spacetime structures as corollary. All the fundamental gravitational structures in fact-the metric as much as the coframes and connections-acquire a spacetime deformation induced theoretical interpretation. There is another line of reasoning which supports the side of this method. We address the theory of teleparallel gravity and construct a consistent Einstein-Cartan (EC) theory with the dynamical torsion.
We show that the equations of the standard EC theory, in which the equation defining torsion is the algebraic type and, in fact, no propagation of torsion is allowed, can be equivalently replaced by the set of modified Einstein-Cartan equations in which the torsion, in general, is dynamical. Moreover, the special physical constraint imposed upon the spacetime deformations yields the short-range propagating spin-spin interaction. For the self-contained arguments in Appendix A.1 and Appendices B and C we complete the   (2) (12)     spacetime deformation theory [82] by new investigation of building up the distortion-complex of spacetime continuum and showing how it restores the world-deformation tensor, which still has been put in by hand. We extend necessary geometrical ideas of spacetime deformation in concise form, without going into the subtleties, as applied to the gravitation theory which underlies the MTBH framework. I have attempted to maintain a balance between being overly detailed and overly schematic. Therefore the text in the appendices should resemble a "hybrid" of a new investigation and some issues of proposed gravitation theory.

A.1. A First Glance at Spacetime Deformation. Consider a smooth deformation map Ω :
4 →M 4 , written in terms of the world-deformation tensor (Ω), the general (M 4 ), and flat ( 4 ) smooth differential 4D-manifolds. The following notational conventions will be used throughout the appendices. All magnitudes related to the space,M 4 , will be denoted by an over "̃". We use the Greek alphabet In the standard theory of gravitation they can be determined from the standard field equations by means of the general linear frames (C.10). However, it should be emphasized that the standard Riemannian space interacting quantum field theory cannot be a satisfactory ground for addressing the most important processes of rearrangement of vacuum state and gauge symmetry breaking in gravity at huge energies. The difficulties arise there because Riemannian geometry, in general, does not admit a group of isometries, and it is impossible to define energy-momentum as Noether local currents related to exact symmetries. This, in turn, posed severe problem of nonuniqueness of the physical vacuum and the associated Fock space. A definition of positive frequency modes cannot, in general, be unambiguously fixed in the past and future, which leads to |in⟩ ̸ = |out⟩, because the state |in⟩ is unstable against decay into many particle |out⟩ states due to interaction processes allowed by lack of Poincaré invariance. A nontrivial Bogolubov transformation between past and future positive frequency modes implies that particles are created from the vacuum and this is one of the reasons for |in⟩ ̸ = |out⟩.

A.2. General Gauge Principle.
Keeping in mind the aforesaid, we develop the alternative framework of the general gauge principle (GGP), which is the distortion gauge induced fiberbundle formulation of gravitation. As this principle was in use as a guide in constructing our theory, we briefly discuss its general implications in Appendix D. The interested reader   is invited to consult the original paper [74] for details. In this, we restrict ourselves to consider only the simplest spacetime deformation map,Ω : . This theory accounts for the gravitation gauge group generated by the hidden local internal symmetry loc .
We assume that a distortion massless gauge field ( ) (≡ ( )) has to act on the external spacetime groups. This field takes values in the Lie algebra of the abelian group loc . We pursue a principle goal of building up the worlddeformation tensor,Ω( ) =̃( )̃( ), where is the  differential form of gauge field = (1/2) ∧ . We connect the structure group , further, to the nonlinear realization of the Lie group of distortion of extended space 6 ( →̃6) (E.1), underlying the 4 . This extension appears to be indispensable for such a realization. In using the 6D language, we will be able to make a necessary reduction to the conventional 4D space. The laws guiding this redaction are given in Appendix E. The nonlinear realization technique or the method of phenomenological Lagrangians [85][86][87][88][89][90][91] provides a way to determine the transformation properties of fields defined on the quotient space. In accordance, we treat the distortion group and its stationary subgroup = (3), respectively, as the dynamical group and its algebraic subgroup. The fundamental field is distortion gauge field (a) and, thus, all the fundamental gravitational structures in fact-the metric as much as the coframes and connectionsacquire a distortion-gauge induced theoretical interpretation. We study the geometrical structure of the space of parameters in terms of Cartan's calculus of exterior forms and derive the Maurer-Cartan structure equations, where the distortion fields (a) are treated as the Goldstone fields. ( )] of (1). This group leads to the renormalizable theory, because gauge invariance gives a conservation of charge, and it also ensures the cancelation of quantum corrections that would otherwise result in infinitely large amplitudes. This has two generators, the third component 3 of isospin ⃗ related to the Pauli spin matrix ⃗ /2, and hypercharge implying = 3 + /2, where is the distortion charge operator assigning the number −1 to particles, but +1 to antiparticles. The group (A.3) entails two neutral gauge bosons of (1), or that coupled to 3 , and of (1) , or that coupled to the hypercharge . Spontaneous symmetry breaking can be achieved by introducing the neutral complex scalar Higgs field. Minimization of the vacuum energy fixes the nonvanishing vacuum expectation value (VEV), which spontaneously breaks the theory, leaving the (1) subgroup intact, that is, leaving one Goldstone boson. Consequently, the left Goldstone boson is gauged away from the scalar sector, but it essentially reappears in the gauge sector providing the longitudinally polarized spin state of one of gauge bosons which acquires mass through its coupling to Higgs scalar. Thus, the two neutral gauge bosons were mixed to form two physical orthogonal states of the massless component of distortion field, ( ) ( = 0), which is responsible for gravitational interactions, and its massive component, ( ) ( ̸ = 0), which is responsible for the ID-regime. Hence, a substantial change of the properties of the spacetime continuum besides the curvature may arise at huge energies. This theory is renormalizable, because gauge invariance gives conservation of charge and also ensures the cancelation of quantum corrections that would otherwise result in infinitely large amplitudes. Without careful thought we expect that in this framework the renormalizability of the theory will not be spoiled in curved space-time too, because the infinities arise from ultraviolet properties of Feynman integrals in momentum space which, in coordinate space, are short distance properties, and locally (over short distances) all the curved spaces look like maximally symmetric (flat) space.
A.4. Model Building: Field Equations. The field equations follow at once from the total gauge invariant Lagrangian in terms of Euler-Lagrange variations, respectively, on both curved and flat spaces. The Lagrangian of distortion gauge field ( ) defined on the flat space is undegenerated Killing form on the Lie algebra of the group loc in adjoint representation, which yields the equation of distortion field (F.1). We are interested in the case of a spherical-symmetric gravitational field 0 ( ) in presence of one-dimensional spacelike ID-field (F.6). In the case at hand, one has the group of motions wherẽ(̃) and̃(̃) (̃∈̃3) are taken to denote the internal pressure and macroscopic density of energy defined in proper frame of reference that is being used. The equations of gravitation ( 0 ) and ID ( ) fields can be given in Feynman gauge [71] as where , ⃗ , and and̃,̃⃗ , and̃are ordinary and distorted energy, momentum, and mass at rest. Hence the matter found in the ID-region ( ̸ = 0) of space-time continuum has undergone phase transition of II-kind; that is, each particle goes off from the mass shell-a shift of mass and energymomentum spectra occurs upwards along the energy scale. The matter in this state is called protomatter with the thermodynamics differing strongly from the thermodynamics of ordinary compressed matter. The resulting deformed metric on 4 in holonomic coordinate basis takes the form As a working model we assume the SPC-configurations given in Appendix G, which are composed of spherical-symmetric distribution of matter in many-phase stratified states. This is quick to estimate the main characteristics of the equilibrium degenerate barionic configurations and will guide us toward first look at some of the associated physics. The simulations confirm in brief the following scenario [71]: the energy density and internal pressure have sharply increased in protomatter core of SPC-configuration (with respect to corresponding central values of neutron star) proportional to gravitational forces of compression. This counteracts the collapse and equilibrium holds even for the masses ∼10 9 ⊙ . This feature can be seen, for example, from Figure 7 where the state equation of the II-class SPC II configuration, with the quark protomatter core, is plotted.

B. A Hard Look at Spacetime Deformation
The holonomic metric onM 4  spacetime, there is, however, a need to introduce the soldering tools, which are the linear frames and forms in tangent fiberbundles to the external smooth differential manifold, whose components are so-called tetrad (vierbein) fields. TheM 4 has at each point a tangent space,̃̃̃4, spanned by the anholonomic orthonormal frame field,̃, as a shorthand for the collection of the 4-tuplet (̃0, . . . ,̃3), wherẽ=̃̃. We use the first half of Latin alphabet ( , , , . . . = 0, 1, 2, 3) to denote the anholonomic indices related to the tangent space. The frame field,̃, then defines a dual vector,̃, of differential forms,̃= (̃0 . . . 3 ), as a shorthand for the collection of thẽ =̃̃, whose values at every point form the dual basis, such that̃⌋̃= , where ⌋ denotes the interior product; namely, this is a ∞ -bilinear map ⌋: Ω 1 → Ω 0 with Ω denoting the ∞ -modulo of differential -forms oñ M 4 . In components, we havẽ̃= . On the manifold, M 4 , the tautological tensor field,̃, of type (1, 1) can be defined which assigns to each tangent space the identity linear transformation. Thus, for any point̃∈M 4 and any vector ∈̃̃M 4 , one has̃(̃) =̃. In terms of the frame field, thẽ give the expression for̃as̃=̃̃=̃0 ⊗̃0 + ⋅ ⋅ ⋅̃3 ⊗̃3, in the sense that both sides yield̃when applied to any tangent vector̃in the domain of definition of the frame field. One can also consider general transformations of the linear group, (4, ), taking any base into any other set of four linearly independent fields. The notation, {̃,̃}, will be used below for general linear frames. Let us introduce so-called first deformation matrices, ( ( ) and̃(̃)) ∈ (4,̃) for all̃, as follows: The holonomic metric can be recast in the form .

(B.15)
The anholonomy objects defined on the tangent space, where the anholonomy coefficients, ∘ , which represent the curls of the base members, are (B.17) The (anholonomic) Levi-Civita (or Christoffel) connection can be written as We have then two choices to recast metric as follows: In the first case, the contribution of the Christoffel symbols constructed by the metric̆=̆̆reads where∇ is the covariant derivative. The contravariant deformed metric,̃] , is defined as the inverse of̃], such that̃]̃] = . Hence, the connection deformation Π acts like a force that deviates the test particles from the geodesic motion in the space, 4 . A metric-affine space (̃4,̃,Γ) is defined to have a metric and a linear connection that need not be dependent on each other. In general, the lifting of the constraints of metric-compatibility and symmetry yields the new geometrical property of the spacetime, which are the nonmetricity 1-form̃and the affine torsion 2-form representing a translational misfit (for a comprehensive discussion see [92][93][94][95]). These, together with the curvature 2-form̃, symbolically can be presented as [96,97]

C. Determination of̃and̃in Standard Theory of Gravitation
Let̃=̃∧̃be the 1-forms of corresponding connections assuming values in the Lorentz Lie algebra. The action for gravitational field can be written in the form where the integral is the usual Einstein action, with the coupling constant relating to the Newton gravitational constant ae = 8 / 4 , is the phenomenological action of the spin-torsion interaction, and ⋆ denotes the Hodge dual. This is a ∞ -linear map ⋆ : Ω → Ω − , which acts on the wedge product monomials of the basis 1-forms as ⋆(̃1 ⋅⋅⋅ ) = 1 ⋅⋅⋅̃+ 1 ⋅⋅⋅ . Here we used the abbreviated notations for the wedge product monomials, where ⋆̃is the dual 3-form relating to the canonical energymomentum tensor,̃, by and ⋆Σ = − ⋆Σ is the dual 3-form corresponding to the canonical spin tensor, which is identical with the dynamical spin tensor̃; namely, The variation of the total action,̃=̃+̃, with respect to thẽ,̃andΦ gives the following field equations: for all̃∈ {Ũ }, where ( ) is the identity element of the group . The fiber̃− 1 at̃∈ 4 is diffeomorphic tõ, wherẽis the fiber space, such that̃− 1 (̃) ≡̃̃≈̃. The action of the structure group oñdefines an isomorphism of the Lie algebrag of onto the Lie algebra of vertical vector fields oñtangent to the fiber at each̃∈̃called fundamental. To involve a drastic revision of the role of gauge fields in the physical concept of the spacetime deformation, we generalize the standard gauge scheme by exploring a new special type of distortion gauge field, ( ), which is assumed to act on the external spacetime groups. Then, we also consider the principle fiber bundle, ( 4 , loc ; ), with the base space 4 , the structure group loc , and the surjection . The matter fields Φ( ) take values in the standard fiber which is the Hilbert vector space where a linear representation ( ) of group loc is given. This space can be regarded as the Lie algebra of the group loc upon which the Lie algebra acts according to the law of the adjoint representation: The GGP accounts for the gravitation gauge group generated by the hidden local internal symmetry loc .  In case of zero curvature, one has = = = ( / ), ‖ ‖ ̸ = 0, where are the inertial coordinates. In this, the conventional gauge theory given on the 4 is restored in both curvilinear and inertial coordinates. Although the distortion gauge field ( ) is a vector field, only the gravitational attraction is presented in the proposed theory of gravitation.

E. A Lie Group of Distortion
The extended space 6 reads 0 where is the Kronecker symbol and * = 1 − . Three spatial = × and three temporal 0 = 0 × components are the basis vectors, respectively, in spaces 3 and 3 , where ± = (1/ √ 2)( 0 ± ), 2 0 = − 2 = 1, ⟨ 0 , ⟩ = 0. The 3D space 3 ± is spanned by the coordinates (± ) . In using this language it is important to consider a reduction to the space 4 which can be achieved in the following way.
(1) In case of free flat space 6 , the subspace 3 is isotropic. And in so far it contributes in line element just only by the square of the moduli = |x 0 |, x 0 ∈ 3 , then, the reduction 6 → 4 = 3 ⊕ 1 can be readily achieved if we use = |x 0 | for conventional time.
(2) In case of curved space, the reduction 6 → 4 can be achieved if we use the projection (̆0) of the temporal component (̆0 ) of basis six-vector(̆,̆0 ) on the given universal direction (̆0 →̆0). By this we choose the time coordinate. Actually, the Lagrangian of physical fields defined on 6 is a function of scalars such that ( ) ( ) = + 0 0 ; then upon the reduction of temporal components of six-vectors 0 0 = 0 ⟨̆0 ,̆0 ⟩ 0 = 0 ⟨̆0,̆0⟩ 0 = 0 0 we may fulfill a reduction to 4 .
A distortion of the basis (E.2) comprises the following two steps. We, at first, consider distortion transformations of the ingredient unit vectors under the distortion gauge field ( ):̆( where Q (=Q ( ) ( )) is an element of the group . This induces the distortion transformations of the ingredient unit vectors , which, in turn, undergo the rotations,̆( ) ( ) = is the element of the group of rotations of planes involving two arbitrary axes around the orthogonal third axis in the given ingredient space 3 . In fact, distortion transformations of basis vectors ( ) and ( ) are not independent but rather are governed by the spontaneous breaking of the distortion symmetry (for more details see [74]). To avoid a further proliferation of indices, hereafter we will use uppercase Latin ( ) in indexing ( ), and so forth. where ≡ at given . The generators (E.4) of the group do not complete the group to the dynamical group , and therefore they cannot be interpreted as the generators of the quotien space / , and the distortion fields cannot be identified directly with the Goldstone fields arising in spontaneous breaking of the distortion symmetry . These objections, however, can be circumvented, because, as it is shown by [74], the distortion group = × (3) can be mapped in a one-to-one manner onto the group = (3) × (3), which is isomorphic to the chiral group (2) × (2), in case of which the method of phenomenological Lagrangians is well known. In aftermath, we arrive at the key relation tan = −ae . (E.6) Given the distortion field , the relation (E.6) uniquely determines six angles of rotations around each of six ( ) axes. In pursuing our goal further, we are necessarily led to extending a whole framework of GGP now for the base 12D smooth differentiable manifold: 12 = 6 ⊕ 6 . (E.7) Here the 6 is related to the spacetime continuum (E.1), but the 6 is displayed as a space of inner degrees of freedom. The

F. Field Equations at Spherical Symmetry
The extended field equations followed at once in terms of Euler-Lagrange variations, respectively, on the spaces 12 and̃1 2 [74]. In accordance, the equation of distortion gauge field = ( ( ) , ( ) ) reads where is the energy-momentum tensor and 0 is the gauge fixing parameter. To render our discussion here more transparent, below we clarify the relation between gravitational and coupling constants. To assist in obtaining actual solutions from the field equations, we may consider the weakfield limit and will envisage that the right-hand side of (F.1) should be in the form (F.5) To obtain a feeling for this point we may consider physical systems which are static as well as spherically symmetrical. We are interested in the case of a spherical-symmetric gravitational field 0 ( ) in presence of one-dimensional space-like ID-field :  where 0 ≡ ae 0 , ≡ ae .

G. SPC-Configurations
The equations describing the equilibrium SPC include the gravitational and ID field equations (A.2), the hydrostatic equilibrium equation, and the state equation specified for each domain of many layered configurations. The resulting stable SPC is formed, which consists of the protomatter core and the outer layers of ordinary matter. A layering of configurations is a consequence of the onset of different regimes in equation of state. In the density range < 4.54 × 10 12 g cm −3 , one uses for both configurations the simple semiempirical formula of state equation given by Harrison and Wheeler, see for example [98]. Above the density > 4.54 × 10 12 g cm −3 , for the simplicity, the Iclass SPC I configuration is thought to be composed of regular n-p-e (neutron-proton-electron) gas (in absence of ID) in intermediate density domain 4.54 × 10 12 g cm −3 ≤ < and of the n-p-e protomatter in presence of ID at > . For the II-class SPC II configuration above the density = 4.09 × 10 14 g cm −3 one considers an onset of melting down of hadrons when nuclear matter consequently turns to quark matter, found in string flip-flop regime. In domain ≤ < , to which the distances 0.4 fm < ≤ 1.6 fm correspond, one has the regular (ID is absent) string flip-flop regime. This is a kind of tunneling effect when the strings joining the quarks stretch themselves violating energy conservation and after touching each other they switch on to the other configuration [71]. The basic technique adopted for calculation of transition matrix element̃is the instanton technique (semiclassical treatment). During the quantum transition from a state 1 of energỹ1 to another one 2 of energỹ2, the lowering of energy of system takes place and the quark matter acquires Δ̃correction to the classical string energy such that the flip-flop energy lowers the energy of quark matter, consequently by lowering the critical density or critical Fermi momentum. If one, for example, looks for the string flip-flop transition amplitude of simple system of where is an imaginary time interval and [̃] is the integration over all the possible string motion. The actioñis proportional to the areãof the surface swept by the strings in the finite region of ID-region of 4 . The strings are initially in the -configuration and finally in the -configuration. The maximal contribution to the path integral comes from the surface 0 of the minimum surface area "instanton". A computation of the transition amplitude is straightforward by summing over all the small vibrations around 0 . In domain ≤ < , one has the string flip-flop regime in presence of ID, at distances 0.25 fm < ≤ 0.4 fm. That is, the system is made of quark protomatter in complete -equilibrium with rearrangement of string connections joining them. In final domain > , the system is made of quarks in one bag in complete -equilibrium at presence of ID. The quarks are under the weak interactions and gluons, including the effects of QCD-perturbative interactions. The QCD vacuum has a complicated structure, which is intimately connected to the gluon-gluon interaction. In most applications, sufficient accuracy is obtained by assuming that all the quarks are almost massless inside a bag. The latter is regarded as noninteracting Fermi gas found in the ID-region of the space-time continuum, at short distances ≤ 0.25 fm. Each configuration is defined by the two free parameters of central values of particle concentratioñ(0) and dimensionless potential of space-like ID-field (0). The interior gravitational potential int 0 ( ) matches the exterior one ext 0 ( ) at the surface of the configuration. The central value of the gravitational potential 0 (0) can be found by reiterating integrations when the sewing condition of the interior and exterior potentials holds. The key question of stability of SPC was studied in [72]. In the relativistic case the total mass-energy of SPC is the extremum in equilibrium for all configurations with the same total number of baryons. While the extrema of̃and occur at the same point in a one-parameter equilibrium sequence, one can look for the extremum of̃=̃2 −̃on equal footing. Minimizing the energy will give the equilibrium configuration, and the second derivative of̃will give stability information. Recall that, for spherical configurations of matter, instantaneously at rest, small radial deviations from equilibrium are governed by a Sturm-Liouville linear eigenvalue equation [98], with the imposition of suitable boundary conditions on normal modes with time dependence ( ⃗ , ) = ( ⃗ ) . A necessary and sufficient condition for stability is that the potential energy be positive defined for all initial data of ( ⃗ , 0), namely, in first order approximation when one does not take into account the rotation and magnetic field, if the square of frequency of normal mode of small perturbations is positive. A relativity tends to destabilize configurations. However, numerical integrations of the stability equations of SPC [72] give for the pressure-averaged value of the adiabatic index Γ 1 = ( lñ/ lñ) the following values: Γ 1 ≈ 2.216 for 28 Journal of Astrophysics the SPC I and Γ 1 ≈ 2.4 for SPC II configurations. This clearly proves the stability of resulting SPC. Note that the SPC is always found inside the event horizon sphere, and therefore it could be observed only in presence of accreting matter.

Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.