Photospheric emission may originate from relativistic outflows in two qualitatively different regimes:\ last scattering of photons inside the outflow at the photospheric radius, or radiative diffusion to the boundary of the outflow. In this work the measurement of temperature and flux of the thermal component in the early afterglows of several gamma-ray bursts (GRBs) along with the total flux in the prompt phase are used to determine initial radii of the outflow as well as its Lorentz factors. Results indicate that in some cases the outflow has relatively low Lorentz factors $\Gamma<10$, favouring cocoon interpretation, while in other cases Lorentz factors are larger $\Gamma>10$, indicating diffusive photospheric origin of the thermal component, associated with an ultrarelativistic outflow.

The phenomenon of Bose-Einstein condensation is traditionally associated with and experimentally verified for low temperatures: either of nano-Kelvin scale for alkali atoms [1-3] or room temperatures for quasi-particles [4,5] or photons in two dimensions [6]. Here we demonstrate out of first principles that for certain initial conditions non-equilibrium plasma at relativistic temperatures of billions of Kelvin undergoes condensation, predicted by Zeldovich and Levich in their seminal work [7]. We determine the necessary conditions for the onset of condensation and discuss the possibilities to observe such a phenomenon in laboratory and astrophysical conditions.

Within the binary-driven hypernova I (BdHN I) scenario, the gamma-ray burst GRB190114C originates in a binary system composed of a massive carbon-oxygen core (CO$_{core}$), and a binary neutron star (NS) companion. As the CO$_{core}$ undergoes a supernova explosion with the creation of a new neutron star ($\nu$NS), hypercritical accretion occurs onto the companion binary neutron star until it exceeds the critical mass for gravitational collapse. The formation of a black hole (BH) captures $10^{57}$ baryons by enclosing them within its horizon, and thus a cavity of approximately $10^{11}$ cm is formed around it with initial density $10^{-7}$ g/cm$^3$. A further depletion of baryons in the cavity originates from the expansion of the electron-positron-photon ($e^{+}e^{-}\gamma$) plasma formed at the collapse, reaching a density of $10^{-14}$ g/cm$^3$ by the end of the interaction. It is demonstrated here using an analytical model complemented by a hydrodynamical numerical simulation that part of the $e^{+}e^{-}\gamma$ plasma is reflected off the walls of the cavity. The consequent outflow and its observed properties are shown to coincide with the featureless emission occurring in a time interval of duration $t_{rf}$, measured in the rest frame of the source, between $11$ and $20$ s of the GBM observation. Moreover, similar features of the GRB light curve were previously observed in GRB 090926A and GRB 130427A, all belonging to the BdHN I class. This interpretation supports the general conceptual framework presented in Ruffini et al. (2019) and guarantees that a low baryon density is reached in the cavity, a necessary condition for the operation of the "inner engine" of the GRB presented in an accompanying article (Ruffini & Moradi 2019).

We motivate how the most recent progress in the understanding the nature of the GeV radiation in most energetic gamma-ray bursts (GRBs), the binary-driven hypernovae (BdHNe), has led to the solution of a forty years unsolved problem in relativistic astrophysics: how to extract the rotational energy from a Kerr black hole for powering synchrotron emission and ultra high-energy cosmic rays. The "inner engine" is identified in the proper use of a classical solution introduced by Wald in 1974 duly extended to the most extreme conditions found around the newborn black hole in a BdHN. The energy extraction process occurs in a sequence impulsive processes each accelerating protons to $10^{21}$ eV in a timescale of $10^{-6}$ s and in presence of an external magnetic field of $10^{14}$ G. Specific example is given for a black hole of initial angular momentum $J=0.3\,M^2$ and mass $M\approx 3\,M_\odot$ leading to the GeV radiation of $10^{49}$ erg$\cdot$s$^{-1}$. The process can energetically continue for thousands of years.

We analyze GRB 151027A within the binary-driven hypernova (BdHN) approach, with progenitor a carbon-oxygen core on the verge of a supernova (SN) explosion and a binary companion neutron star (NS). The hypercritical accretion of the SN ejecta onto the NS leads to its gravitational collapse into a black hole (BH), to the emission of the GRB and to a copious $e^+e^-$ plasma. The impact of this $e^+e^-$ plasma on the SN ejecta explains {the} early soft X-ray flare observed in long GRBs. We here apply this approach to the UPE and to the hard X-ray flares. We use GRB 151027A as a prototype. From the time-integrated and the time-resolved analysis we identify a double component in the UPE and confirm its ultra-relativistic nature. We confirm the mildly-relativistic nature of the soft X-ray flare, of the hard X-ray flare and of the ETE. We show that the ETE identifies the transition from a SN to the HN. We then address the theoretical justification of these observations by integrating the hydrodynamical propagation equations of the $e^+ e^-$ into the SN ejecta, the latter independently obtained from 3D smoothed-particle-hydrodynamics simulations. We conclude that the UPE, the hard X-ray flare and the soft X-ray flare do not form a causally connected sequence: Within our model they are the manifestation of \textbf{the same} physical process of the BH formation as seen through different viewing angles, implied by the morphology and the $\sim 300$~s rotation period of the HN ejecta.

The non-equilibrium electron-positron-photon plasma thermalization process is studied using relativistic Boltzmann solver, taking into account quantum corrections both in non-relativistic and relativistic cases. Collision integrals are computed from exact QED matrix elements for all binary and triple interactions in the plasma. It is shown that in non-relativistic case (temperatures $k_B T\leq 0.3 m_e c^2$) binary interaction rates dominate over triple ones, resulting in establishment of the kinetic equilibrium prior to final relaxation towards the thermal equilibrium, in agreement with the previous studies. On the contrary, in relativistic case (final temperatures $k_B T\geq 0.3 m_e c^2$) triple interaction rates are fast enough to prevent the establishment of kinetic equilibrium. It is shown that thermalization process strongly depends on quantum degeneracy in initial state, but does not depend on plasma composition.

We present a new efficient method to compute Uehling-Uhlenbeck collision integral for all two-particle interactions in relativistic plasma with drastic improvement in computation time with respect to existing methods. Plasma is assumed isotropic in momentum space. The set of reactions consists of: Moeller and Bhabha scattering, Compton scattering, two-photon pair annihilation, and two-photon pair production, which are described by QED matrix elements. In our method exact energy and particle number conservation laws are satisfied. Reaction rates are compared, where possible, with the corresponding analytical expressions and convergence of numerical rates is demonstrated.

Photon-photon scattering of gamma-rays on the cosmic microwave background has been studied using the low energy approximation of the total cross section by Zdziarski & Svensson (1989); Svensson & Zdziarski (1990). Here, the cosmic horizon due to photon-photon scattering is accurately determined using the exact cross section and we find that photon-photon scattering dominates over the pair production at energies smaller than 1.68 GeV and at redshifts larger than 180.

We analyze the early X-ray flares in the GRB "flare-plateau-afterglow" (FPA) phase observed by Swift-XRT. The FPA occurs only in one of the seven GRB subclasses: the binary-driven hypernovae (BdHNe). This subclass consists of long GRBs with a carbon-oxygen core and a neutron star (NS) binary companion as progenitors. The hypercritical accretion of the supernova (SN) ejecta onto the NS can lead to the gravitational collapse of the NS into a black hole. Consequently, one can observe a GRB emission with isotropic energy $E_{iso}\gtrsim10^{52}$~erg, as well as the associated GeV emission and the FPA phase. Previous work had shown that gamma-ray spikes in the prompt emission occur at $\sim 10^{15}$--$10^{17}$~cm with Lorentz gamma factor $\Gamma\sim10^{2}$--$10^{3}$. Using a novel data analysis we show that the time of occurrence, duration, luminosity and total energy of the X-ray flares correlate with $E_{iso}$. A crucial feature is the observation of thermal emission in the X-ray flares that we show occurs at radii $\sim10^{12}$~cm with $\Gamma\lesssim 4$. These model independent observations cannot be explained by the "fireball" model, which postulates synchrotron and inverse Compton radiation from a single ultra relativistic jetted emission extending from the prompt to the late afterglow and GeV emission phases. We show that in BdHNe a collision between the GRB and the SN ejecta occurs at $\simeq10^{10}$~cm reaching transparency at $\sim10^{12}$~cm with $\Gamma\lesssim4$. The agreement between the thermal emission observations and these theoretically derived values validates our model and opens the possibility of testing each BdHN episode with the corresponding Lorentz gamma factor.

This paper summarizes the limits on propagation of ultra high energy particles in the Universe, set up by their interactions with cosmic background of photons and neutrinos. By taking into account cosmic evolution of these backgrounds and considering appropriate interactions we derive the mean free path for ultra high energy photons, protons and neutrinos. For photons the relevant processes are the Breit-Wheeler process as well as the double pair production process. For protons the relevant reactions are the photopion production and the Bethe-Heitler process. We discuss the interplay between the energy loss length and mean free path for the Bethe-Heitler process. Neutrino opacity is determined by its scattering off the cosmic background neutrino. We compute for the first time the high energy neutrino horizon as a function of its energy.

The problem of the expansion of a relativistic plasma generated when a large amount of energy is released in a small volume has been considered by many authors. We use the analytical solution of Bisnovatyi-Kogan and Murzina (1995) for the spherically symmetric relativistic expansion. The light-curves and the spectra from transparency of an electron-positron-photon plasma are obtained. We compare our results with the work of Goodman (1986).

Recent observations, especially by the Fermi satellite, point out the importance of the thermal component in GRB spectra. This fact revives strong interest in photospheric emission from relativistic outflows. Early studies already suggested that the observed spectrum of photospheric emission from relativistically moving objects differs in shape from the Planck spectrum. However, this component appears to be subdominant in many GRBs and the origin of the dominant component is still unclear. One of the popular ideas is that energy dissipation near the photosphere may produce a non-thermal spectrum and account for such emission. Before considering such models, though, one has to determine precise spectral and timing characteristics of the photospheric emission in the simplest possible case. Hence this paper focuses on various physical effects which make the photospheric emission spectrum different from the black body spectrum and quantifies them.

Optically thick energy dominated plasma created in the source of Gamma-Ray Bursts (GRBs) expands radially with acceleration and forms a shell with constant width measured in the laboratory frame. When strong Lorentz factor gradients are present within the shell it is supposed to spread at sufficiently large radii. There are two possible mechanisms of spreading: hydrodynamical and thermal ones. We consider both mechanisms evaluating the amount of spreading that occurs during expansion up to the moment when the expanding shell becomes transparent for photons. We compute the hydrodynamical spreading of an ultrarelativistically expanding shell. In the case of thermal spreading we compute the velocity spread as a function of two parameters: comoving temperature and bulk Lorentz factor of relativistic Maxwellian distribution. Based on this result we determine the value of thermal spreading of relativistically expanding shell. We found that thermal spreading is negligible for typical GRB parameters. Instead hydrodynamical spreading appears to be significant, with the shell width reaching $\sim10^{10}$ cm for total energy $E=10^{54}$ erg and baryonic loading $B=10^{-2}$. Within the fireshell model such spreading will result in the duration of Proper Gamma-Ray Bursts up to several seconds.

We consider the formation of photon spectrum at the photosphere of ultrarelativistically expanding outflow. We use the Fokker-Planck approximation to the Boltzmann equation, and obtain the generalized Kompaneets equation which takes into account anisotropic distribution of photons developed near the photosphere. This equation is solved numerically for relativistic steady wind and the observed spectrum is found in agreement with previous studies. We also study the photospheric emission for different temperature dependences on radius in such outflows. In particular, we found that for $T\propto r^{-2}$ the Band low energy photon index of the observed spectrum is $\alpha\simeq -1$, as typically observed in Gamma Ray Bursts.

Electron-positron plasma is believed to play imporant role both in the early Universe and in sources of Gamma-Ray Bursts (GRBs). We focus on analogy and difference between physical conditions of electron-positron plasma in the early Universe and in sources of GRBs. We discuss a) dynamical differences, namely thermal acceleration of the outflow in GRB sources vs cosmological deceleration; b) nuclear composition differences as synthesis of light elements in the early Universe and possible destruction of heavy elements in GRB plasma; c) different physical conditions during last scattering of photons by electrons. Only during the acceleration phase of the optically thick electron-positron plasma comoving observer may find it similar to the early Universe. This similarity breaks down during the coasting phase. Reprocessing of nuclear abundances may likely take place in GRB sources. Heavy nuclear elements are then destroyed, resulting mainly in protons with small admixture of helium. Unlike the primordial plasma which recombines to form neutral hydrogen, and emits the Cosmic Microwave Background Radiation, GRB plasma does not cool down enough to recombine.

We studied the decoupling of photons from ultra-relativistic spherically symmetric outflows expanding with constant velocity by means of Monte-Carlo (MC) simulation. For outflows with finite width we confirm the existence of two regimes: photon thick and photon thin introduced recently by Ruffini, Siutsou, Vereshchagin (2011), hereafter RSV. The probability density function of photon last scattering is shown to be very different in these two cases. We also obtained spectra as well as light curves. In photon thick case, the time integrated spectrum is much broader than the Planck function and its shape is well described by the fuzzy photosphere approximation introduced by RSV. In the photon thin case we confirm the crucial role of photon diffusion, hence the probability density of decoupling has a maximum near the diffusion radius, well below the photosphere. Its spectrum has Band shape. It is produced when the outflow is optically thick and its peak is formed at diffusion radius.

We study the process of energy conversion from overcritical electric field into electron-positron-photon plasma. We solve numerically Vlasov-Boltzmann equations for pairs and photons assuming the system to be homogeneous and anisotropic. All the 2-particle QED interactions between pairs and photons are described by collision terms. We evidence several epochs of this energy conversion, each of them associated to a specific physical process. Firstly pair creation occurs, secondly back reaction results in plasma oscillations. Thirdly photons are produced by electron-positron annihilation. Finally particle interactions lead to completely equilibrated thermal electron-positron-photon plasma.

(Shortened) GRB080319B, with an isotropic energy E_{iso}=1.32x10^{54}erg, and GRB050904, with E_{iso}=1.04x10^{54}erg, offer the possibility of studying the spectral properties of the prompt radiation of two of the most energetic Gamma-Ray Bursts (GRBs). This allows us to probe the validity of the fireshell model for GRBs beyond 10^{54}erg, well outside the energy range where it has been successfully tested up to now (10^{49}-10^{53}erg). We find that in the low energy region, the prompt emission spectra observed by Swift BAT reveals more power than theoretically predicted. The opportunities offered by these observations to improve the fireshell model are outlined. One of the distinguishing features of the fireshell model is that it relates the observed spectra to the spectrum in the comoving frame of the fireshell. Originally, a fully radiative condition and a comoving thermal spectrum were adopted. An additional power-law in the comoving thermal spectrum is required [...] in the fireshell model for GRBs 080319B and 050904. A new phenomenological parameter \alpha is correspondingly introduced in the model. We perform numerical simulations of the prompt emission in the Swift BAT bandpass by assuming different values of \alpha [...]. We compare them with the GRB080319B and GRB050904 observed time-resolved spectra, as well as with their time-integrated spectra and light curves. Although GRB080319B and GRB050904 are at very different redshifts (z=0.937 and z=6.29 respectively), a value of \alpha=-1.8 leads for both of them to a good agreement between the numerical simulations and the observed BAT light curves, time-resolved and time-integrated spectra. Such a modified spectrum is also consistent with the observations of previously analyzed less energetic GRBs and reasons for this additional agreement are given. Perspectives for future low energy missions are outlined.

In this paper we reexamine the optical depth of ultrarelativistic spherically symmetric outflows and reevaluate the photospheric radius for each model during both the acceleration and coasting phases. It is shown that for both the wind and the shell models there are two asymptotic solutions for the optical depth during the coasting phase of the outflow. In particular we show that quite counterintuitively a geometrically thin shell may appear as a thick wind for photons propagating inside it. For this reason we introduce notions of photon thick and photon thin outflows, which appear more general and better physically motivated with respect to winds and shells. Photosphere of relativistic outflow is a dynamic surface. We study its geometry and find that the photosphere of photon thin outflow has always a convex shape, while in the photon thick one it is initially convex (there is always a photon thin layer in any outflow) and then it becomes concave asymptotically approaching the photosphere of an infinitely long wind. We find that both instantaneous and time integrated observed spectra are very close to the thermal one for photon thick outflows, in line with existing studies. It is our main finding that the photospheric emission from the photon thin outflow produces non thermal time integrated spectra, which may be described by the Band function well known in the GRB literature. We find that energetic GRBs should produce photon thin outflows with photospheric emission lasting less than one second for the total energy $E_0\leq10^{54}$ erg and baryonic loading parameter $B\leq10^{-2}$. It means that only time integrated spectra may be observed from such GRBs.

We study the frequency of the plasma oscillations of electron-positron pairs created by the vacuum polarization in an uniform electric field with strength E in the range 0.2 Ec < E < 10 Ec. Following the approach adopted in [1] we work out one second order ordinary differential equation for a variable related to the velocity from which we can recover the classical plasma oscillation equation when E -> 0. Thereby, we focus our attention on its evolution in time studying how this oscillation frequency approaches the plasma frequency. The time-scale needed to approach to the plasma frequency and the power spectrum of these oscillations are computed. The characteristic frequency of the power spectrum is determined uniquely from the initial value of the electric field strength. The effects of plasma degeneracy and pair annihilation are discussed.

By numerically solving the relativistic Boltzmann equations, we compute the time scale for relaxation to thermal equilibrium for an optically thick electron-positron plasma with baryon loading. We focus on the time scales of electromagnetic interactions. The collisional integrals are obtained directly from the corresponding QED matrix elements. Thermalization time scales are computed for a wide range of values of both the total energy density (over 10 orders of magnitude) and of the baryonic loading parameter (over 6 orders of magnitude). This also allows us to study such interesting limiting cases as the almost purely electron-positron plasma or electron-proton plasma as well as intermediate cases. These results appear to be important both for laboratory experiments aimed at generating optically thick pair plasmas as well as for astrophysical models in which electron-positron pair plasmas play a relevant role.

From the interaction of physics and astrophysics we are witnessing in these years a splendid synthesis of theoretical, experimental and observational results originating from three fundametal physical processes. They were originally proposed by Dirac, by Breit and Wheeler and by Sauter, Heisenberg, Euler and Schwinger. The vacuum polarization process in strong electromagnetic field, pioneered by Sauter, Heisenberg, Euler and Schwinger, introduced the concept of critical electric field. It has been searched without success for more than forty years by heavy-ion collisions in many of the leading particle accelerators worldwide. The novel situation today is that these same processes can be studied on a much more grandiose scale during the gravitational collapse leading to the formation of a black hole being observed in Gamma Ray Bursts. This report is dedicated to the scientific race in act. The theoretical and experimental work developed in Earth-based laboratories is confronted with the theoretical interpretation of space-based observations of phenomena originating on cosmological scales. What has become clear in the last ten years is that all the three above mentioned processes, duly extended in the general relativistic framework, are necessary for the understanding of the physics of the gravitational collapse to a black hole. Vice versa, the natural arena where these processes can be observed in mutual interaction and on an unprecedented scale, is indeed the realm of relativistic astrophysics.

(Shortened) [...] After recalling the basic features of the "fireshell model", we emphasize the following novel results: 1) the interpretation of the X-ray flares in GRB afterglows as due to the interaction of the optically thin fireshell with isolated clouds in the CircumBurst Medium (CBM); 2) an interpretation as "fake - disguised" short GRBs of the GRBs belonging to the class identified by Norris & Bonnell [...] consistent with an origin from the final coalescence of a binary system in the halo of their host galaxies with particularly low CBM density [...]; 3) the first attempt to study a genuine short GRB with the analysis of GRB 050509B, that reveals indeed still an open question; 4) the interpretation of the GRB-SN association in the case of GRB 060218 via the "induced gravitational collapse" process; 5) a first attempt to understand the nature of the "Amati relation", a phenomenological correlation between the isotropic-equivalent radiated energy of the prompt emission E_{iso} with the cosmological rest-frame \nu F_{\nu} spectrum peak energy E_{p,i}. In addition, recent progress on the thermalization of the electron-positron plasma close to their formation phase, as well as the structure of the electrodynamics of Kerr-Newman Black Holes are presented. An outlook for possible explanation of high-energy phenomena in GRBs to be expected from the AGILE and the Fermi satellites are discussed. As an example of high energy process, the work by Enrico Fermi dealing with ultrarelativistic collisions is examined. It is clear that all the GRB physics points to the existence of overcritical electrodynamical fields. In this sense we present some progresses on a unified approach to heavy nuclei and neutron stars cores, which leads to the existence of overcritical fields under the neutron star crust.

In the recent Letter [1] we considered the approach of nonequilibrium pair plasma towards thermal equilibrium state adopting a kinetic treatment and solving numerically the relativistic Boltzmann equations. It was shown that plasma in the energy range 0.1-10 MeV first reaches kinetic equilibrium, on a timescale t_{k}<10^{-14} sec, with detailed balance between binary interactions such as Compton, Bhabha and Moller scattering, and pair production and annihilation. Later the electron-positron-photon plasma approaches thermal equilibrium on a timescale t_{th}<10^{-12} sec, with detailed balance for all direct and inverse reactions. In the present paper we systematically present details of the computational scheme used in [1], as well as generalize our treatment, considering proton loading of the pair plasma. When proton loading is large, protons thermalize first by proton-proton scattering, and then with the electron-positron-photon plasma by proton-electron scattering. In the opposite case of small proton loading proton-electron scattering dominates over proton-proton one. Thus in all cases the plasma, even with proton admixture, reaches thermal equilibrium configuration on a timescale t_{th}<10^{-11} sec. We show that it is crucial to account for not only binary but also triple direct and inverse interactions between electrons, positrons, photons and protons. Several explicit examples are given and the corresponding timescales for reaching kinetic and thermal equilibria are determined.

The fireshell model for Gamma-Ray Bursts (GRBs) naturally leads to a canonical GRB composed of a proper-GRB (P-GRB) and an afterglow. P-GRBs, introduced by us in 2001, are sometimes considered "precursors" of the main GRB event in the current literature. We show in this paper how the fireshell model leads to the understanding of the structure of GRBs, with precise estimates of the time sequence and intensities of the P-GRB and the of the afterglow. It leads as well to a natural classification of the canonical GRBs which overcomes the traditional one in short and long GRBs.

The Gamma-Ray Bursts (GRBs) offer the unprecedented opportunity to observe for the first time the blackholic energy extracted by the vacuum polarization during the process of gravitational collapse to a black hole leading to the formation of an electron-positron plasma. The uniqueness of the Kerr-Newman black hole implies that very different processes originating from the gravitational collapse a) of a single star in a binary system induced by the companion, or b) of two neutron stars, or c) of a neutron star and a white dwarf, do lead to the same structure for the observed GRB. The recent progress of the numerical integration of the relativistic Boltzmann equations with collision integrals including 2-body and 3-body interactions between the particles offer a powerful conceptual tool in order to differentiate the traditional "fireball" picture, an expanding hot cavity considered by Cavallo and Rees, as opposed to the "fireshell" model, composed of an internally cold shell of relativistically expanding electron-positron-baryon plasma. The analysis of the fireshell naturally leads to a canonical GRB composed of a proper-GRB and an extended afterglow. By recalling the three interpretational paradigms for GRBs we show how the fireshell model leads to an understanding of the GRB structure and to an alternative classification of short and long GRBs.

(Shortened) We show by example how the uncoding of Gamma-Ray Bursts (GRBs) offers unprecedented possibilities to foster new knowledge in fundamental physics and in astrophysics. After recalling some of the classic work on vacuum polarization in uniform electric fields by Klein, Sauter, Heisenberg, Euler and Schwinger, we summarize some of the efforts to observe these effects in heavy ions and high energy ion collisions. We then turn to the theory of vacuum polarization around a Kerr-Newman black hole, leading to the extraction of the blackholic energy, to the concept of dyadosphere and dyadotorus, and to the creation of an electron-positron-photon plasma. We then present a new theoretical approach encompassing the physics of neutron stars and heavy nuclei. It is shown that configurations of nuclear matter in bulk with global charge neutrality can exist on macroscopic scales and with electric fields close to the critical value near their surfaces. These configurations may represent an initial condition for the process of gravitational collapse, leading to the creation of an electron-positron-photon plasma: the basic self-accelerating system explaining both the energetics and the high energy Lorentz factor observed in GRBs. We then turn to recall the two basic interpretational paradigms of our GRB model. [...] We then turn to the special role of the baryon loading in discriminating between "genuine" short and long or "fake" short GRBs [...] We finally turn to the GRB-Supernova Time Sequence (GSTS) paradigm: the concept of induced gravitational collapse. [...] We then present some general conclusions.

We discuss the temporal evolution of the pair plasma created in Gamma-Ray Burst sources. A particular attention is paid to the relaxation of the plasma into thermal equilibrium. We also discuss the connection between the dynamics of expansion and the spatial geometry of the plasma. The role of the baryonic loading parameter is emphasized.

Starting from a nonequilibrium configuration we analyse the essential role of the direct and the inverse binary and triple interactions in reaching an asymptotic thermal equilibrium in a homogeneous isotropic electron-positron-photon plasma. We focus on energies in the range 0.1--10 MeV. We numerically integrate the integro-partial differential relativistic Boltzmann equation with the exact QED collisional integrals taking into account all binary and triple interactions in the plasma. We show that first, when detailed balance is reached for all binary interactions on a timescale $t_{k}\lesssim10^{-14}$sec, photons and electron-positron pairs establish kinetic equilibrium. Successively, when triple interactions fulfill the detailed balance on a timescale $t_{eq}\lesssim10^{-12}$sec, the plasma reaches thermal equilibrium. It is shown that neglecting the inverse triple interactions prevents reaching thermal equilibrium. Our results obtained in the theoretical physics domain also find application in astrophysics and cosmology.

Analysis of cosmological solutions of GX dark energy models shows the presence of a separatrix in the phase space of solutions which divides them into two classes: Friedmannian-like with initial singularity and non-Friedmannian solutions. The invariants are found, that reveal the intrinsic symmetry in GX models.

We derive luminosity distance equation in Gurzadyan-Xue cosmological models and compared it with available supernovae and radio galaxies data sets. We found that the luminosity distance does not depend explicitly the speed of light and the gravitation constant, and depends only on the matter parameter (GX-invariant) and curvature.

We evidence the existence of plasma oscillations of electrons-positron pairs created by the vacuum polarization in an uniform electric field with E < Ec. Our general treatment, encompassing also the traditional, well studied case of E > Ec, shows the existence in both cases of a maximum Lorentz factor acquired by electrons and positrons and allows determination of the a maximal length of oscillation. We quantitatively estimate how plasma oscillations reduce the rate of pair creation and increase the time scale of the pair production. These results are particularly relevant in view of the experimental progress in approaching the field strengths E < Ec.

We outline the main results of our GRB model, based on the three interpretation paradigms we proposed in July 2001, comparing and contrasting them with the ones in the current literature. Thanks to the observations by Swift and by VLT, this analysis points to a "canonical GRB" originating from markedly different astrophysical scenarios. The communality is that they are all emitted in the formation of a black hole with small or null angular momentum. The following sequence appears to be canonical: the vacuum polarization process creating an optically thick self accelerating electron-positron plasma; the engulfment of baryonic mass during the plasma expansion; the adiabatic expansion of the optically thick "fireshell" up to the transparency; the interaction of the remaining accelerated baryons with the interstellar medium (ISM). This leads to the canonical GRB composed of a proper GRB (P-GRB), emitted at the moment of transparency, followed by an extended afterglow. The parameters are the plasma total energy, the fireshell baryon loading and the ISM filamentary distribution around the source. In the limit of no baryon loading the total energy is radiated in the P-GRB. In this limit, the canonical GRBs explain as well the short GRBs.

We compare and contrast the different approaches to the optically thick adiabatic phase of GRB all the way to the transparency. Special attention is given to the role of the rate equation to be self consistently solved with the relativistic hydrodynamic equations. The works of Shemi and Piran (1990), Piran, Shemi and Narayan (1993), Meszaros, Laguna and Rees (1993) and Ruffini, Salmonson, Wilson and Xue (1999,2000) are compared and contrasted. The role of the baryonic loading in these three treatments is pointed out. Constraints on initial conditions for the fireball produced by electro-magnetic black hole are obtained.

The luminosity and the spectral distribution of the afterglow of GRB 031203 have been presented within our theoretical framework, which envisages the GRB structure as composed by a proper-GRB, emitted at the transparency of an electron-positron plasma with suitable baryon loading, and an afterglow comprising the "prompt emission" as due to external shocks. In addition to the GRB emission, there appears to be a prolonged soft X-Ray emission lasting for 10^6-10^7 seconds followed by an exponential decay. This additional source has been called by us URCA-3. It is urgent to establish if this component is related to the GRB or to the Supernova (SN). In this second case, there are two possibilities: either the interaction of the SN ejecta with the interstellar medium or, possibly, the cooling of a young neutron star formed in the SN 2003lw process. The analogies and the differences between this triptych GRB 031203 / SN 2003lw / URCA-3 and the corresponding ones GRB 980425 / SN 1998bw / URCA-1 and GRB 030329 / SN 2003dh / URCA-2, as well as GRB 060218 / SN 2006aj are discussed.

Perturbation theory within Newtonian approximation is presented for cosmological models with varying physical constants. Analytical solutions for perturbations dynamics are obtained for each Gurzadyan-Xue model with pressureless matter and radiation. We found that perturbations grow during entire expansion within GX models, including curvature- and vacuum-dominated stage when they cease to grow in the standard cosmological model.

Non-perturbative quantum geometric effects in Loop Quantum Cosmology predict a $\rho^2$ modification to the Friedmann equation at high energies. The quadratic term is negative definite and can lead to generic bounces when the matter energy density becomes equal to a critical value of the order of the Planck density. The non-singular bounce is achieved for arbitrary matter without violation of positive energy conditions. By performing a qualitative analysis we explore the nature of the bounce for inflationary and Cyclic model potentials. For the former we show that inflationary trajectories are attractors of the dynamics after the bounce implying that inflation can be harmoniously embedded in LQC. For the latter difficulties associated with singularities in cyclic models can be overcome. We show that non-singular cyclic models can be constructed with a small variation in the original Cyclic model potential by making it slightly positive in the regime where scalar field is negative.

The formula for dark energy density derived by Gurzadyan and Xue provides a value of density parameter of dark energy in remarkable agreement with current cosmological datasets, unlike numerous phenomenological dark energy scenarios where the corresponding value is postulated. This formula suggests the possibility of variation of physical constants such as the speed of light and the gravitational constant. Considering several cosmological models based on that formula and deriving the cosmological equations for each case, we show that, in all models source terms appear in the continuity equation. So, one one hand, GX models make up a rich set covering a lot of currently proposed models of dark energy, on the other hand, they reveal hidden symmetries, with a particular role of the separatrix $\Omega_m=2/3$, and link with the issue of the content of physical constants.

The expansion of the electron-positron plasma in the GRB phenomenon is compared and contrasted in the treatments of Meszaros, Laguna and Rees, of Shemi, Piran and Narayan, and of Ruffini et al. The role of the correct numerical integration of the hydrodynamical equations, as well as of the rate equation for the electron-positron plasma loaded with a baryonic mass, are outlined and confronted for crucial differences.

In the previous papers we derived cosmological equations for models with Gurzadyan-Xue dark energy, have performed their qualitative analysis and, particularly, have revealed a remarkable hidden invariance in the models with respect to the separatrix $\Omega_{sep}$ in their phase portraits. Now, new analytic solutions for these models are obtained, showing additional symmetries at various curvatures. The likelihood analysis with supernovae and radio galaxies data and their characteristics (age, the deceleration parameter), again demonstrate the crucial role of the $\Omega_{sep}$ for all models, in spite of the diversity of both, the initial equations and their solutions.

The dark energy formula derived by Gurzadyan and Xue which leads to a value fitting the SN data, provides a scaling relation between the physical constants and cosmological parameters and defines a set of cosmological models. In previous works we have considered several of those models and derived the cosmological equations for each case. In this letter, we present the phase portrait analysis of those models. Surprisingly we found, first, that the separatrix in the phase space which determines the character of solutions depends solely on the value of the current matter density. Namely, at $\Omega_m>2/3$ the equations describe Friedmannian Universe with the classical singularity at the beginning. While at $\Omega_m<2/3$ all solutions for all models start with zero density and non vanishing scale factor. Secondly, more remarkable, the value $\Omega_{sep}=2/3$ defining the separatrix is the same for all models, which reveales an underlying invariance hidden in the models, possibly, due to the basic nature of the GX-scaling.

We consider cosmological implications of the formula for the dark energy density derived by Gurzadyan and Xue which predicts a value fitting the observational one. Cosmological models with varying by time physical constants, namely, speed of light and gravitational constant and/or their combinations, are considered. In one of the models, for example, vacuum energy density induces effective negative curvature, while another one has an unusual asymptotic. This analysis also explicitely rises the issue of the meaning and content of physical units and constants in cosmological context.

We use the Wilkinson Microwave Anisotropy Probe (WMAP) data on the spectrum of cosmic microwave background anisotropies to put constraints on the present amount of lepton asymmetry L, parameterized by the dimensionless chemical potential (also called degeneracy parameter) xi and on the effective number of relativistic particle species. We assume a flat cosmological model with three thermally distributed neutrino species having all the same mass and chemical potential, plus an additional amount of effectively massless exotic particle species. The extra energy density associated to these species is parameterized through an effective number of additional species DeltaN_oth$. We find that 0<|xi|<1.1 and correspondingly 0<|L|<0.9 at 2sigma, so that WMAP data alone cannot firmly rule out scenarios with a large lepton number; moreover, a small preference for this kind of scenarios is actually found. We also discuss the effect of the asymmetry on the estimation of other parameters and in particular of the neutrino mass. In the case of perfect lepton symmetry, we obtain the standard results. When the amount of asymmetry is left free, we find \sum m_nu < 3.6 eV at 2sigma. Finally we study how the determination of |L| is affected by the assumptions on DeltaN_oth. We find that lower values of the extra energy density allow for larger values of the lepton asymmetry, effectively ruling out, at 2sigma level, lepton symmetric models with DeltaN_oth=0.

In this paper we provide approximate analytical analysis of stability of nonsingular inflationary chaotic-type cosmological models. Initial conditions for nonsingular solutions at the bounce correspond to dominance of potential part of the energy density of the scalar field over its kinetic part both within general relativity and gauge theories of gravity. Moreover, scalar field at the bounce exceeds the planckian value and on expansion stage these models correspond to chaotic inflation. Such solutions can be well approximated by explicitly solvable model with constant effective potential (cosmological term) and massless scalar field during the bounce and on stages of quasi-exponential contraction and expansion. Perturbative analysis shows that nonsingular inflationary solutions are exponentially unstable during contraction stage. This result is compared with numerical calculations.

New type of nonsingular oscillating solutions for the Universe described by cosmological equations of gauge theories of gravity is presented. Advantages of these solutions with respect to existing nonsingular solutions within framework of general relativity and gauge gravity are discussed. It is shown in particular that these solutions have nonzero measure and stable on contraction stage unlike usual nonsingular solutions.

Solutions of gravitational equations of gauge theories of gravity in homogeneous isotropic world with massive scalar field are investigated in the case of flat cosmological models. Special attention is dedicated to general behavior of solutions on contraction stage. It is shown, that on expansion stage inflationary solutions are generic feature of the model. At the same time on contraction stage an effective equation of state is similar to the case of massless scalar field $p=\rho$ at the beginning of evolution and tends to the equation of ultrarelativistic gas $p=\rho/3$. The Hubble parameter tends to some negative value on contraction stage, depending on the mass of the scalar field. Nonsingular solutions in this model are unstable.

In addition to the problem of galaxy formation, one of the greatest open questions of cosmology is represented by the existence of an asymmetry between matter and antimatter in the baryonic component of the Universe. We believe that a net lepton number for the three neutrino species can be used to understand this asymmetry. This also implies an asymmetry in the matter-antimatter component of the leptons. The existence of a nonnull lepton number for the neutrinos can easily explain a cosmological abundance of neutrinos consistent with the one needed to explain both the rotation curves of galaxies and the flatness of the Universe. Some propedeutic results are presented in order to attack this problem.