It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

We consider a process where two identical massive particles fall inwards, starting from rest at infinity towards the extremal Kerr black hole, collide outside the event horizon in its vicinity and produce two massless particles. The center of mass energy of collision between the two particles diverges if one of the particle admits a specific critical value of the angular momentum and if the collision takes place at a location arbitrarily close to the event horizon. Assuming the isotropic emission of particles in the center of mass frame we show that one of the massless particles produced has divergent conserved energy comparable to the center of mass energy with probability slightly less than half. This particle enters the black hole event horizon which coincides with the Cauchy horizon, turns back, and emerges through the white hole event horizon into another asymptotic region in the maximal extension of Kerr spacetime. Since conserved energy is preserved in this process, it is perceived as a particle with divergent energy by the observer when it reaches infinity. Thus the extremal white hole appears to be a source of ultra high energy particles. Similar processes wherein collision takes place slightly inside the black hole event horizon or just outside or inside of white hole event horizon also produce high energy particles.

We perform (3+1)-dimensional simulations of primordial black hole (PBH) formation starting from the spheroidal super-horizon perturbations. We investigate how the ellipticity (prolateness or oblateness) affects the threshold of PBH formation in terms of the peak amplitude of curvature perturbation. We find that, in the case of the radiation-dominated universe, the effect of ellipticity on the threshold is negligibly small for large amplitude of perturbations expected for PBH formation.

The production rate of primordial black holes is often calculated by considering a nearly Gaussian distribution of cosmological perturbations, and assuming that black holes will form in regions where the amplitude of such perturbations exceeds a certain threshold. A threshold $\zeta_{\rm th}$ for the curvature perturbation is somewhat inappropriate for this purpose, because it depends significantly on environmental effects, not essential to the local dynamics. By contrast, a threshold $\delta_{\rm th}$ for the density perturbation at horizon crossing seems to provide a more robust criterion. On the other hand, the density perturbation is known to be bounded above by a maximum limit $\delta_{\rm max}$, and given that $\delta_{\rm th}$ is comparable to $\delta_{\rm max}$, the density perturbation will be far from Gaussian near or above the threshold. In this paper, we provide a new plausible estimate for the primordial black hole abundance based on peak theory. In our approach, we assume that the curvature perturbation is given as a random Gaussian field with the power spectrum characterized by a single scale, while an optimized criterion for PBH formation is imposed, based on the locally averaged density perturbation. Both variables are related by the full nonlinear expression derived in the long-wavelength approximation of general relativity. We do not introduce a window function, and the scale of the inhomogeneity is introduced as a random variable in the peak theory. We find that the mass spectrum is shifted to larger mass scales by one order of magnitude or so, compared to a conventional calculation. The abundance of PBHs becomes significantly larger than the conventional one, by many orders of magnitude, mainly due to the optimized criterion for PBH formation and the removal of the suppresion associated with a window function.

In the previous paper [Harada, Cardoso, and Miyata, Phys.\ Rev.\ D {\bf 99} (2019), 044039], it is shown that a hollow transmissive shell collapsing to an ultracompact object of radius very close to its horizon radius generally emits transient Hawking radiation followed by a couple of bursts separated each other by a long time interval. In the current paper, we expand the previous work in two independent directions: changing boundary conditions and specifying the equations of state of the matter. First, we introduce a perfectly reflective surface collapsing to an ultracompact object and find that this model also emits transient Hawking radiation that is followed only by a single burst. Second, we introduce two different collapse dynamics to an ultracompact object and specify the corresponding matter equations of state. We find that transient Hawking radiation is quite commonly seen in early times, while the subsequent bursts strongly depend on the boundary condition and the equation of state or the braking behavior of the surface.

Angular momentum plays very important roles in the formation of PBHs in the matter-dominated phase if it lasts sufficiently long. In fact, most collapsing masses are bounced back due to centrifugal force, since angular momentum significantly grows before collapse. As a consequence, most of the formed PBHs are rapidly rotating near the extreme value $a_{*}=1$, where $a_{*}$ is the nondimensional Kerr parameter at their formation. The smaller the density fluctuation $\sigma_{H}$ at horizon entry is, the stronger the tendency towards the extreme rotation. Combining the effect of angular momentum with that of anisotropy, we estimate the black hole production rate. We find that the production rate suffers from suppression dominantly due to angular momentum for a smaller value of $\sigma_{H}$, while due to anisotrpopy for a larger value of $\sigma_{H}$. We argue that matter domination significantly enhances the production of PBHs despite the suppression. If the matter-dominated phase does not last so long, the effect of the finite duration significantly suppresses PBH formation and weakens the tendency towards large spins. (abridged)

We investigate the effect of inhomogeneity on primordial black hole formation in the matter dominated era. In the gravitational collapse of an inhomogeneous density distribution, a black hole forms if apparent horizon prevents information of the central region of the configuration from leaking. Since information cannot propagate faster than the speed of light, we identify the threshold of the black hole formation by considering the finite speed for propagation of information. We show that the production probability $\beta_{inhom}(\sigma)$ of primordial black holes, where $\sigma$ is density fluctuation at horizon entry, is significantly enhanced from that derived in previous work in which the speed of propagation was effectively regarded as infinite. For $\sigma \ll 1$, we obtain $\beta_{inhom}\simeq 3.70 \sigma^{3/2}$, which is larger by about an order of magnitude than the probability derived in earlier work by assuming instantaneous propagation of information.

A strong gravity naked singular region can give important clues towards understanding classical as well as spontaneous nature of General Relativity. We propose here a model for energy emission from a naked singular region in a self-similar dust spacetime by gluing two self-similar dust solutions at the Cauchy horizon. The energy is defined and evaluated as a surface energy of a null hypersurface, the null shell. Also included are scenarios of spontaneous creation or disappearance of a singularity, end of inflation, black hole formation and bubble nucleation. Our examples investigated here explicitly show that one can model unlimitedly luminous and energetic objects in the framework of General Relativity.

We completely classify Friedmann-Lema\^{i}tre-Robertson-Walker solutions with spatial curvature $K=0,\pm 1$ and equation of state $p=w\rho$, according to their conformal structure, singularities and trapping horizons. We do not assume any energy conditions and allow $\rho < 0$, thereby going beyond the usual well-known solutions. For each spatial curvature, there is an initial spacelike big-bang singularity for $w>-1/3$ and $\rho>0$, while no big-bang singularity for $w<-1$ and $\rho>0$. For $K=0$ or $-1$, $-1<w<-1/3$ and $\rho>0$, there is an initial null big-bang singularity. For each spatial curvature, there is a final spacelike future big-rip singularity for $w<-1$ and $\rho>0$, with null geodesics being future complete for $-5/3\le w<-1$ but incomplete for $w<-5/3$. For $w=-1/3$, the expansion speed is constant. For $-1<w<-1/3$ and $K=1$, the universe contracts from infinity, then bounces and expands back to infinity. For $K=0$, the past boundary consists of timelike infinity and a regular null hypersurface for $-5/3<w<-1$, while it consists of past timelike and past null infinities for $w\le -5/3$. For $w<-1$ and $K=1$, the spacetime contracts from an initial spacelike past big-rip singularity, then bounces and blows up at a final spacelike future big-rip singularity. For $w<-1$ and $K=-1$, the past boundary consists of a regular null hypersurface. The trapping horizons are timelike, null and spacelike for $w\in (-1,1/3)$, $w\in \{1/3, -1\}$ and $w\in (-\infty,-1)\cup (1/3,\infty)$, respectively. A negative energy density ($\rho <0$) is possible only for $K=-1$. In this case, for $w>-1/3$, the universe contracts from infinity, then bounces and expands to infinity; for $-1<w<-1/3$, it starts from a big-bang singularity and contracts to a big-crunch singularity; for $w<-1$, it expands from a regular null hypersurface and contracts to another regular null hypersurface.

We consider a head-on collision of two massive particles that move in the equatorial plane of an extremal Kerr black hole, which results in the production of two massless particles. Focusing on a typical case, where both of the colliding particles have zero angular momenta, we show that a massless particle produced in such a collision can escape to infinity with arbitrarily large energy in the near-horizon limit of the collision point. Furthermore, if we assume that the emission of the produced massless particles is isotropic in the center-of-mass frame but confined to the equatorial plane, the escape probability of the produced massless particle approaches $5/12$ and almost all escaping massless particles have arbitrarily large energy at infinity and an impact parameter approaching $2GM/c^2$, where $M$ is the mass of the black hole.

We investigate primordial black hole formation in the matter-dominated phase of the Universe, where nonspherical effects in gravitational collapse play a crucial role. This is in contrast to the black hole formation in a radiation-dominated era. We apply the Zel'dovich approximation, Thorne's hoop conjecture, and Doroshkevich's probability distribution and subsequently derive the production probability $\beta_{0}$ of primordial black holes. The numerical result obtained is applicable even if the density fluctuation $\sigma$ at horizon entry is of the order of unity. For $\sigma\ll 1$, we find a semi-analytic formula $\beta_{0}\simeq 0.05556 \sigma^{5}$, which is comparable with the Khlopov-Polnarev formula. We find that the production probability in the matter-dominated era is much larger than that in the radiation-dominated era for $\sigma\lesssim 0.05$, while they are comparable with each other for $\sigma\gtrsim 0.05$. We also discuss how $\sigma$ can be written in terms of primordial curvature perturbations.

We construct spacetimes which provide spherical and nonspherical models of black hole formation in the flat Friedmann--Lemaitre--Robertson--Walker (FLRW) universe with the Lemaitre--Tolman--Bondi solution and the Szekeres quasispherical solution, respectively. These dust solutions may contain both shell-crossing and shell-focusing naked singularities. These singularities can be physically regarded as the breakdown of dust description, where strong pressure gradient force plays a role. We adopt the simultaneous big bang condition to extract a growing mode of adiabatic perturbation in the flat FLRW universe. If the density perturbation has a sufficiently homogeneous central region and a sufficiently sharp transition to the background FLRW universe, its central shell-focusing singularity is globally covered. If the density concentration is sufficiently large, no shell-crossing singularity appears and a black hole is formed. If the density concentration is not sufficiently large, a shell-crossing singularity appears. In this case, a large dipole moment significantly advances shell-crossing singularities and they tend to appear before the black hole formation. In contrast, a shell-crossing singularity unavoidably appears in the spherical and nonspherical evolution of cosmological voids. The present analysis is general and applicable to cosmological nonlinear structure formation described by these dust solutions.

In an accretion of fluid, its velocity may transit from subsonic to supersonic. The point at which such transition occurs is called sonic point and often mathematically special. We consider a steady-state and spherically symmetric accretion problem of ideal photon gas in general static spherically symmetric spacetime neglecting back reaction. Our main result is that the EOS of ideal photon gas leads to correspondence between its sonic point and the photon sphere of the spacetime in general situations. Despite of the dependence of the EOS on the dimension of spacetime, this correspondence holds for spacetimes of arbitrary dimensions.

We propose a consistent analytic approach to the efficiency of collisional Penrose process in the vicinity of a maximally rotating Kerr black hole. We focus on a collision with arbitrarily high center-of-mass energy, which occurs if either of the colliding particles has its angular momentum fine-tuned to the critical value to enter the horizon. We show that if the fine-tuned particle is ingoing on the collision, the upper limit of the efficiency is $(2+\sqrt{3})(2-\sqrt{2})\simeq 2.186$, while if the fine-tuned particle is bounced back before the collision, the upper limit is $(2+\sqrt{3})^{2}\simeq 13.93$. Despite earlier claims, the former can be attained for inverse Compton scattering if the fine-tuned particle is massive and starts at rest at infinity, while the latter can be attained for various particle reactions, such as inverse Compton scattering and pair annihilation, if the fine-tuned particle is either massless or highly relativistic at infinity. We discuss the difference between the present and earlier analyses.

The origin of the ultra-high-energy particles we receive on the Earth from the outer space such as EeV cosmic rays and PeV neutrinos remains an enigma. All mechanisms known to us currently make use of electromagnetic interaction to accelerate charged particles. In this paper we propose a mechanism exclusively based on gravity rather than electromagnetic interaction. We show that it is possible to generate ultra-high-energy particles starting from particles with moderate energies using the collisional Penrose process in an overspinning Kerr spacetime transcending the Kerr bound only by an infinitesimal amount, i.e., with the Kerr parameter $a=M(1+\epsilon)$, where we take the limit $\epsilon \rightarrow 0^+$. We consider two massive particles starting from rest at infinity that collide at $r=M$ with divergent center-of-mass energy and produce two massless particles. We show that massless particles produced in the collision can escape to infinity with the ultra-high energies exploiting the collisional Penrose process with the divergent efficiency $\eta \sim {1}/{\sqrt{\epsilon}} \rightarrow \infty$. Assuming the isotropic emission of massless particles in the center-of-mass frame of the colliding particles, we show that half of the particles created in the collisions escape to infinity with the divergent energies. To a distant observer, ultra-high-energy particles appear to originate from a bright spot which is at the angular location $\xi \sim {2M}/{r_{obs}}$ with respect to the singularity on the side which is rotating towards the observer. We show that the anisotropy in emission in the center-of-mass frame, which is dictated by the differential cross-section of underlying particle physics process, leaves a district signature on the spectrum of ultra-high-energy massless particles. Thus, it provides a unique probe into fundamental particle physics.

Primordial black holes (PBHs) are those which may have formed in the early Universe and affected the subsequent evolution of the Universe through their Hawking radiation and gravitational field. To constrain the early Universe from the observational constraint on the abundance of PBHs, it is essential to determine the formation threshold for primordial cosmological fluctuations, which are naturally described by cosmological long-wavelength solutions. I will briefly review our recent analytical and numerical results on the PBH formation.

The center-of-mass energy of two particles can become arbitrarily large if they collide near the event horizon of an extremal Kerr black hole, which is called the Ba$\rm \tilde n$ados-Silk-West (BSW) effect. We consider such a high-energy collision of two particles which started from infinity and follow geodesics in the equatorial plane and investigate the energy extraction from such a high-energy particle collision and the production of particles in the equatorial plane. We analytically show that, on the one hand, if the produced particles are as massive as the colliding particles, the energy-extraction efficiency is bounded by $2.19$ approximately. On the other hand, if a very massive particle is to be produced as a result of the high-energy collision, which has negative energy and necessarily falls into the black hole, the upper limit of the energy-extraction efficiency is increased to $(2+\sqrt{3})^2 \simeq 13.9$. Thus, higher efficiency of the energy extraction, which is typically as large as 10, provides strong evidence for the production of a heavy particle.

The effect of the Gauss-Bonnet term on the existence and dynamical stability of thin-shell wormholes as negative tension branes is studied in the arbitrary dimensional spherically, planar, and hyperbolically symmetric spacetimes. We consider radial perturbations against the shell for the solutions which have the Z${}_2$ symmetry and admit the general relativistic limit. It is shown that the Gauss-Bonnet term shrinks the parameter region admitting static wormholes. The effect of the Gauss-Bonnet term on the stability depends on the spacetime symmetry. For planar symmetric wormholes, the Gauss-Bonnet term does not affect their stability. If the coupling constant is positive but small, the Gauss-Bonnet term tends to destabilize spherically symmetric wormholes, while it stabilizes hypebolically symmetric wormholes. The Gauss-Bonnet term can destabilize hypebolically symmetric wormholes as a non-perturbative effect, however, spherically symmetric wormholes cannot be stable.

We make a critical comparison between ultra-high energy particle collisions around an extremal Kerr black hole and that around an over-spinning Kerr singularity, mainly focusing on the issue of the timescale of collisions. We show that the time required for two massive particles with the proton mass or two massless particles of GeV energies to collide around the Kerr black hole with Planck energy is several orders of magnitude longer than the age of the Universe for astro-physically relevant masses of black holes, whereas time required in the over-spinning case is of the order of ten million years which is much shorter than the age of the Universe. Thus from the point of view of observation of Planck scale collisions, the over-spinning Kerr geometry, subject to their occurrence, has distinct advantage over their black hole counterparts.

We construct cosmological long-wavelength solutions without symmetry in general gauge conditions compatible with the long-wavelength scheme. We then specify the relationship among the solutions in different time slicings. Applying this general framework to spherical symmetry, we derive the correspondence relation between long-wavelength solutions in the constant mean curvature slicing with conformally flat spatial coordinates and asymptotic quasihomogeneous solutions in the comoving gauge and compare the numerical results of PBH formation in these two different approaches. To discuss the PBH formation, it is convenient and conventional to use $\tilde{\delta}_{c}$, the value which the averaged density perturbation at threshold in the comoving slicing would take at horizon entry in the lowest-order long-wavelength expansion. We numerically find that within compensated models, the sharper the transition from the overdense region to the FRW universe is, the larger the $\tilde{\delta}_{c}$ becomes. We suggest that, for the equation of state $p=(\Gamma-1)\rho$, we can apply the analytic formulas for the minimum $\tilde{\delta}_{c, {\rm min}}\simeq [3\Gamma/(3\Gamma+2)]\sin^{2}\left[\pi\sqrt{\Gamma-1}/(3\Gamma-2)\right]$ and the maximum $\tilde{\delta}_{c, {\rm max}}\simeq 3\Gamma/(3\Gamma+2)$. As for the threshold peak value of the curvature variable $\psi_{0,c}$, we find that the sharper the transition is, the smaller the $\psi_{0,c}$ becomes. We analytically explain this feature. Using simplified models, we also analytically deduce an environmental effect that $\psi_{0,c}$ can be significantly larger (smaller) if the underlying density perturbation of much longer wavelength is positive (negative).

We investigate negative tension branes as stable thin shell wormholes in Reissner-Nordstrom-(anti) de Sitter spacetimes in $d$ dimensional Einstein gravity. Imposing Z2 symmetry, we construct and classify traversable static thin shell wormholes in spherical, planar (or cylindrical) and hyperbolic symmetries. In spherical geometry, we find the higher dimensional counterpart of Barcelo and Visser's wormholes, which are stable against spherically symmetric perturbations. We also find the classes of thin shell wormholes in planar and hyperbolic symmetries with a negative cosmological constant, which are stable against perturbations preserving symmetries. In most cases, stable wormholes are found with the combination of an electric charge and a negative cosmological constant. However, as special cases, we find stable wormholes even with vanishing cosmological constant in spherical symmetry and with vanishing electric charge in hyperbolic symmetry.

Rapidly rotating Kerr black holes can accelerate particles to arbitrarily high energy if the angular momentum of the particle is fine-tuned to some critical value. This phenomenon is robust as it is founded on the basic properties of geodesic orbits around a near-extremal Kerr black hole. On the other hand, the maximum energy of the acceleration is subjected to several physical effects. There is convincing evidence that the particle acceleration to arbitrarily high energy is one of the universal properties of general near-extremal black holes. We also discuss gravitational particle acceleration in more general context. This article is intended to provide a pedagogical introduction to and a brief overview of this topic for non-specialists.

The claim that an overdense (positive curvature) region in the early universe cannot extend beyond some maximum scale and remain part of our universe, first made 40 years ago, has recently been questioned by Kopp et al. Their analysis is elucidating and demonstrates that one cannot constrain the form of primordial density perturbations using this argument. However, the notion of a separate-universe scale still applies and it places an important upper limit on the mass of primordial black holes forming at any epoch. We calculate this scale for equations of state of the form $p = k \rho c^2$ with $-1 <k < \infty$, refining earlier calculations on account of the Kopp et al. criticisms. For $-1/3 < k < \infty$, the scale is always of order the cosmological particle horizon size, with a numerical factor depending on $k$. This confirms the earlier claim that a primordial black hole cannot be much larger than the particle horizon at formation. For $-1 < k< -1/3$, as expected for some periods in the history of the universe, the situation changes radically, in that a sufficiently large positive-curvature region produces a baby universe rather than a black hole. There is still a separate-universe scale but the interpretation of these solutions requires care.

Primordial black holes (PBHs) are an important tool in cosmology to probe the primordial spectrum of small-scale curvature perturbations that reenter the cosmological horizon during radiation domination epoch. We numerically solve the evolution of spherically symmetric highly perturbed configurations to clarify the criteria of PBHs formation using an extremely wide class of curvature profiles characterized by five parameters, (in contrast to only two parameters used in all previous papers) which specify the curvature profiles not only at the central region but also at the outer boundary of configurations. It is shown that formation or non-formation of PBHs is determined entirely by only two master parameters one of which can be presented as an integral of curvature over initial configurations and the other is presented in terms of the position of the boundary and the edge of the core.

We study high energy particle collisions around the higher dimensional black holes. It is shown that the center of mass energy can be arbitrarily large in the vicinity of the event horizon like the Banados, Silk and West (BSW) process in four dimensions if the following two conditions are satisfied: (i) the horizon is extremal and (ii) the parameters of either of the two colliding particles are fine-tuned, which is called a critical particle condition. We also show that a test particle which satisfies the critical particle condition can reach the event horizon from a distant region for a simple case. Finally, we discuss the relation between the BSW process and the linear instability of test fields around the extremal black holes, the latter of which has been recently found by Aretakis.

Based on a physical argument, we derive a new analytic formula for the amplitude of density perturbation at the threshold of primordial black hole formation in the universe dominated by a perfect fluid with the equation of state $p=w\rho c^{2}$ for $w\ge 0$. The formula gives $\delta^{\rm UH}_{H c}=\sin^{2}[\pi \sqrt{w}/(1+3w)]$ and $\tilde{\delta}_{c}=[3(1+w)/(5+3w)]\sin^{2}[\pi\sqrt{w}/(1+3w)]$, where $\delta^{\rm UH}_{H c}$ and $\tilde{\delta}_{c}$ are the amplitude of the density perturbation at the horizon crossing time in the uniform Hubble slice and the amplitude measure used in numerical simulations, respectively, while the conventional one gives $\delta^{\rm UH}_{H c}=w$ and $\tilde{\delta}_{c}=3w(1+w)/(5+3w)$. Our formula shows a much better agreement with the result of recent numerical simulations both qualitatively and quantitatively than the conventional formula. For a radiation fluid, our formula gives $\delta^{\rm UH}_{H c}=\sin^{2}(\sqrt{3}\pi/6)\simeq 0.6203$ and $\tilde{\delta}_{c}=(2/3)\sin^{2}(\sqrt{3}\pi/6)\simeq 0.4135$. We also discuss the maximum amplitude and the cosmological implications of the present result.

For particle collision near rapidly rotating Kerr black holes, the center-of-mass energy can be arbitrarily high if the angular momentum of either of the colliding particles is fine-tuned. Recently, it has been shown that particles which are produced by such a particle collision and escape to infinity cannot be very massive nor very energetic. For electrically charged black holes there is a similar phenomenon, where the center-of-mass energy for the collision of charged particles near the horizon can be arbitrarily high. One might expect that there would exist a similar bound on the energy and mass of particles that are produced by such a particle collision and escape to infinity. In this paper, however, we see that this expectation is not the case. We explicitly show that superheavy and highly energetic charged particles produced by the collision near maximally charged black holes can escape to infinity at least within classical theory if the backreaction and self-force of the particle can be neglected.

We show a numerical scheme to solve the moment equations of the radiative transfer, i.e., M1 model which follows the evolution of the energy density, $ E $, and the energy flux, $ \mbox{\boldmath$F$} $. In our scheme we reconstruct the intensity from $ E $ and $ \mbox{\boldmath$F$} $ so that it is consistent with the closure relation, relation, $ \chi = (3 + 4 f ^2)/(5 + 2 \sqrt{4 - 3 f ^2}) $. Here the symbols, $ \chi $, $ f = |\mbox{\boldmath$F$}|/(cE) $, and $ c $, denote the Eddington factor, the reduced flux, and the speed of light, respectively. We evaluate the numerical flux across the cell surface from the kinetically reconstructed intensity. It is an explicit function of $ E $ and $ \mbox{\boldmath$F$} $ in the neighboring cells across the surface considered. We include absorption and reemission within a numerical cell in the evaluation of the numerical flux. The numerical flux approaches to the diffusion approximation when the numerical cell itself is optically thick. Our numerical flux gives a stable solution even when some regions computed are very optically thick. We show the advantages of the numerical flux with examples. They include flash of beamed photons and irradiated protoplanetary disks.

The possible violation of the weak cosmic censorship hypothesis in astrophysical phenomena can provide us with the information about trans-Planckian physics through observations. We present negative evidence, however, that one should not expect such a possibility at least when the deviation from spherical symmetry is negligible and the parameter values of collapse are astrophysically reasonable. Taking the Lema\^itre-Tolman-Bondi solution as the model most likely to counter the weak hypothesis, we show that the mass ($ \gtrsim 1.5 M_\odot $) and density ($ \gtrsim 1.5 \times 10^{15} \; {\rm g/cm^3}$) of the collapsing object produce a gravitational field strong enough to capture any null rays soon after emanating from the singularity.

It is well known that the sum of signed magnifications is invariant for mass lens systems. In this paper, we discuss the signed magnification sums of general spherical lens models including the singular isothermal sphere, the Schwarzschild lens and the Ellis wormhole, the last of which is an example of the traversable wormholes of the Morris-Thorne class. We show that the signed magnication sums are a very useful tool to distinguish exotic lens objects. For example, we show that one can distinguish the Ellis wormholes from the Schwarzschild lens with the signed magnification sums.

For the last decade, the gravitational lensing in the strong gravitational field has been studied eagerly. It is well known that, for the lensing by a black hole, infinite number of Einstein rings are formed by the light rays which wind around the black hole nearly on the photon sphere, which are called relativistic Einstein rings. This is also the case for the lensing by a wormhole. In this paper, we study the Einstein ring and relativistic Einstein rings for the Schwarzschild black hole and the Ellis wormhole, the latter of which is an example of traversable wormholes of the Morris-Thorne class. Given the configuration of the gravitational lensing and the radii of the Einstein ring and relativistic Einstein rings, we can distinguish between a black hole and a wormhole in principle. We conclude that we can detect the relativistic Einstein rings by wormholes which have the radii of the throat $a\simeq 0.5$pc at a galactic center with the distance 10Mpc and which have $a\simeq 10$AU in our galaxy using by the most powerful modern instruments which have the resolution of $10^{-2}$arcsecond such as a 10-meter optical-infrared telescope. The black holes which make the Einstein rings of the same size as the ones by the wormholes are galactic supermassive black holes and the relativistic Einstein rings by the black holes are too small to measure at this moment. We may test some hypotheses of astrophysical wormholes by using the Einstein ring and relativistic Einstein rings in the future.

The center-of-mass energy of two particles colliding near the horizon of a maximally rotating black hole can be arbitrarily high if the angular momentum of either of the incident particles is fine-tuned, which we call a critical particle. We study particle emission from such high-energy collision and reaction in the equatorial plane fully analytically. We show that the unconditional upper limit of the energy of the emitted particle is given by 218.6% of that of the injected critical particle, irrespective of the details of the reaction and this upper limit can be realized for massless particle emission. The upper limit of the energy extraction efficiency for this emission as a collisional Penrose process is given by 146.6%, which can be realized in the collision of two massive particles with optimized mass ratio. Moreover, we analyze perfectly elastic collision, Compton scattering, and pair annihilation and show that net positive energy extraction is really possible for these three reactions. The Compton scattering is most efficient among them and the efficiency can reach 137.2%. On the other hand, our result is qualitatively consistent with the earlier claim that the mass and energy of the emitted particle are at most of order the total energy of the injected particles and hence we can observe neither super-heavy nor super-energetic particles.

We study high energy charged particle collisions near the horizon in an electromagnetic field around a rotating black hole and reveal the condition of the fine-tuning to obtain arbitrarily large center-of-mass (CM) energy. We demonstrate that the CM energy can be arbitrarily large as the uniformly magnetized rotating black hole arbitrarily approaches maximal rotation under the situation that a charged particle plunges from the innermost stable circular orbit (ISCO) and collides with another particle near the horizon. Recently, Frolov [Phys. Rev. D 85, 024020 (2012)] proposed that the CM energy can be arbitrarily high if the magnetic field is arbitrarily strong, when a particle collides with a charged particle orbiting the ISCO with finite energy near the horizon of a uniformly magnetized Schwarzschild black hole. We show that the charged particle orbiting the ISCO around a spinning black hole needs arbitrarily high energy in the strong field limit. This suggests that Frolov's process is unstable against the black hole spin. Nevertheless, we see that magnetic fields may substantially promote the capability of rotating black holes as particle accelerators in astrophysical situations.

An inspiralling object of mass $\mu$ around a Kerr black hole of mass $M (\gg \mu)$ experiences a continuous transition near the innermost stable circular orbit from adiabatic inspiral to plunge into the horizon as gravitational radiation extracts its energy and angular momentum. We investigate the collision of such an object with a generic counterpart around a Kerr black hole. We find that the angular momentum of the object is fine-tuned through gravitational radiation and that the high-velocity collision of the object with a generic counterpart naturally occurs around a nearly maximally rotating black hole. We also find that the centre-of-mass energy can be far beyond the Planck energy for dark matter particles colliding around a stellar mass black hole and as high as $10^{58}$ erg for stellar mass compact objects colliding around a supermassive black hole, where the present transition formalism is well justified. Therefore, rapidly rotating black holes can accelerate objects inspiralling around them to energy high enough to be of great physical interest.

We obtain an explicit expression for the center-of-mass (CM) energy of two colliding general geodesic massive and massless particles at any spacetime point around a Kerr black hole. Applying this, we show that the CM energy can be arbitrarily high only in the limit to the horizon and then derive a formula for the CM energy of two general geodesic particles colliding near the horizon in terms of the conserved quantities of each particle and the polar angle. We present the necessary and sufficient condition for the CM energy to be arbitrarily high in terms of the conserved quantities of each particle. To have an arbitrarily high CM energy, the angular momentum of either of the two particles must be fine-tuned to the critical value $L_{i}=\Omega_{H}^{-1}E_{i}$, where $\Omega_{H}$ is the angular velocity of the horizon and $E_{i}$ and $L_{i}$ are the energy and angular momentum of particle $i$ ($=1,2$), respectively. We show that, in the direct collision scenario, the collision with an arbitrarily high CM energy can occur near the horizon of maximally rotating black holes not only at the equator but also on a belt centered at the equator. This belt lies between latitudes $\pm acos(\sqrt{3}-1)\simeq \pm 42.94^{\circ}$. This is also true in the scenario through the collision of a last stable orbit particle.

We derive a general formula for the center-of-mass (CM) energy for the near-horizon collision of two particles of the same rest mass on the equatorial plane around a Kerr black hole. We then apply this formula to a particle which plunges from the innermost stable circular orbit (ISCO) and collides with another particle near the horizon. It is found that the maximum value of the CM energy $E_{\rm cm}$ is given by $E_{\rm cm}/(2m_{0})\simeq 1.40/\sqrt[4]{1-a_{*}^{2}}$ for a nearly maximally rotating black hole, where $m_{0}$ is the rest mass of each particle and $a_{*}$ is the nondimensional Kerr parameter. This coincides with the known upper bound for a particle which begins at rest at infinity within a factor of 2. Moreover, we also consider the collision of a particle orbiting the ISCO with another particle on the ISCO and find that the maximum CM energy is then given by $E_{\rm cm}/(2m_{0})\simeq 1.77/\sqrt[6]{1-a_{*}^{2}}$. In view of the astrophysical significance of the ISCO, this result implies that particles can collide around a rotating black hole with an arbitrarily high CM energy without any artificial fine-tuning in an astrophysical context if we can take the maximal limit of the black hole spin or $a_{*}\to 1$. On the other hand, even if we take Thorne's bound on the spin parameter into account, highly or moderately relativistic collisions are expected to occur quite naturally, for $E_{\rm cm}/(2m_{0})$ takes 6.95 (maximum) and 3.86 (generic) near the horizon and 4.11 (maximum) and 2.43 (generic) on the ISCO for $a_{*}=0.998$. This implies that high-velocity collisions of compact objects are naturally expected around a rapidly rotating supermassive black hole. Implications to accretion flows onto a rapidly rotating black hole are also discussed.

We investigate static spherically symmetric vacuum solutions in the IR limit of projectable nonrelativistic quantum gravity, including the renormalisable quantum gravity recently proposed by Ho\v{r}ava. It is found that the projectability condition plays an important role. Without the cosmological constant, the spacetime is uniquely given by the Schwarzschild solution. With the cosmological constant, the spacetime is uniquely given by the Kottler (Schwarzschild-(anti) de Sitter) solution for the entirely vacuum spacetime. However, in addition to the Kottler solution, the static spherical and hyperbolic universes are uniquely admissible for the locally empty region, for the positive and negative cosmological constants, respectively, if its nonvanishing contribution to the global Hamiltonian constraint can be compensated by that from the nonempty or nonstatic region. This implies that static spherically symmetric entirely vacuum solutions would not admit the freedom to reproduce the observed flat rotation curves of galaxies. On the other hand, the result for locally empty regions implies that the IR limit of nonrelativistic quantum gravity theories does not simply recover general relativity but includes it.

In loop quantum cosmology, the Hamiltonian reduces to a finite difference operator. We study the initial singularity and the large volume limit against the ambiguities in the discretisation and the operator ordering within a homogeneous, isotropic and spatially flat model with the cosmological constant. We find that the absence of the singularity strongly depends on the choice of the operator ordering and the requirement for the absence singles out a very small class of orderings. Moreover we find a general ordering rule required for the absence of the singularity. We also find that the large volume limit naturally recovers a smooth wave function in the discretisation where each step corresponds to a fixed volume increment but not in the one where each step corresponds to a fixed area increment. If loop quantum cosmology is to be a phenomenological realisation of full loop quantum gravity, these results are important to fix the theoretical ambiguities.

In this paper we continue our study on the accretion process onto super-spinning Kerr objects with no event horizon (super-spinars). We discuss the counterpart of the Bondi accretion onto black holes. We first report the results of our numerical simulations. We found a quasi steady-state configuration for any choice of the parameters of our model. The most interesting feature is the presence of hot outflows. Unlike jets and outflows produced around black holes, which are thought to be powered by magnetic fields and emitted from the poles, here the outflows are produced by the repulsive gravitational force at a small distance from the super-spinar and are ejected around the equatorial plane. In some circumstances, the amount of matter in the outflow is considerable, which can indeed significantly reduce the gas mass accretion rate. Finally, we discuss a possible scenario of the accretion process in more realistic situations, which cannot be simulated by our code.

This review addresses the issue of whether there are physically realistic self-similar solutions in which a primordial black hole is attached to an exact or asymptotically Friedmann model for an equation of state of the form $p=(\gamma-1)\rho c^2$. In the positive pressure case ($1 < \gamma < 2$), there is no such solution when the black hole is attached to an exact Friedmann background via a sonic point. However, it has been claimed that there is a one-parameter family of asymptotically Friedmann black hole solutions providing the ratio of the black hole size to the cosmological horizon size is in a narrow range above some critical value. There are also "universal" black holes in which the black hole has an apparent horizon but no event horizon. It turns out that both these types of solution are only asymptotically {\it quasi}-Friedmann, because they contain a solid angle deficit at large distances, but they are not necessarily excluded observationally. We also consider the possibility of self-similar black hole solutions in a universe dominated by a scalar field. If the field is massless, the situation resembles the stiff fluid case, so any black hole solution is again contrived, although there may still be universal black hole solutions. The situation is less clear if the scalar field is rolling down a potential and therefore massive, as in the quintessence scenario. Although no explicit asymptotically Friedmann black hole solutions of this kind are known, they are not excluded and comparison with the $0 < \gamma < 2/3$ perfect fluid case suggests that they should exist if the black hole is not too large. This implies that a black hole might grow as fast as the cosmological horizon in a quintessence-dominated universe in some circumstances, supporting the proposal that accretion onto primordial black holes may have played a role in the production of the supermassive black holes in galactic nuclei.

The specific angular momentum of a Kerr black hole must not be larger than its mass. The observational confirmation of this bound which we call a Kerr bound directly suggests the existence of a black hole. In order to investigate observational testability of this bound by using the X-ray energy spectrum of black hole candidates, we calculate energy spectra for a super-spinning object (or a naked singularity) which is described by a Kerr metric but whose specific angular momentum is larger than its mass, and then compare the spectra of this object with those of a black hole. We assume an optically thick and geometrically thin disc around the super-spinning object and calculate its thermal energy spectrum seen by a distant observer by solving general relativistic radiative transfer equations including usual special and general relativistic effects such as Doppler boosting, gravitational redshift, light bending and frame-dragging. Surprisingly, for a given black hole, we can always find its super-spinning counterpart with its spin $a_*$ in the range $5/3<a_*<8\sqrt{6}/3$ whose observed spectrum is very similar to and practically indistinguishable from that of the black hole. As a result, we conclude that to confirm the Kerr bound we need more than the X-ray thermal spectrum of the black hole candidates.

The accretion process onto spinning objects in Kerr spacetimes is studied with numerical simulations. Our results show that accretion onto compact objects with Kerr parameter (characterizing the spin) $|a| < M$ and $|a| > M$ is very different. In the super-spinning case, for $|a|$ moderately larger than $M$, the accretion onto the central object is extremely suppressed due to a repulsive force at short distance. The accreting matter cannot reach the central object, but instead is accumulated around it, forming a high density cloud that continues to grow. The radiation emitted in the accretion process will be harder and more intense than the one coming from standard black holes; e.g. $\gamma$-rays could be produced as seen in some observations. Gravitational collapse of this cloud might even give rise to violent bursts. As $|a|$ increases, a larger amount of accreting matter reaches the central object and the growth of the cloud becomes less efficient. Our simulations find that a quasi-steady state of the accretion process exists for $|a|/M \gtrsim 1.4$, independently of the mass accretion rate at large radii. For such high values of the Kerr parameter, the accreting matter forms a thin disk at very small radii. We provide some analytical arguments to strengthen the numerical results; in particular, we estimate the radius where the gravitational force changes from attractive to repulsive and the critical value $|a|/M \approx 1.4$ separating the two qualitatively different regimes of accretion. We briefly discuss the observational signatures which could be used to look for such exotic objects in the Galaxy and/or in the Universe.

Motivated by the cosmological wormhole solutions obtained from our recent numerical investigations, we provide a definition of a wormhole which applies to dynamical situations. Our numerical solutions do not have timelike trapping horizons but they are wormholes in the sense that they connect two or more asymptotic regions. Although the null energy condition must be violated for static wormholes, we find that it can still be satisfied in the dynamical context. Two analytic solutions for a cosmological wormhole connecting two Friedmann universes without trapping horizons are presented.

We investigate self-similar solutions which are asymptotic to the Friedmann universe at spatial infinity and contain a scalar field with potential. The potential is required to be exponential by self-similarity. It is found that there are two distinct one-parameter families of asymptotic solutions,one is asymptotic to the proper Friedmann universe, while the other is asymptotic to the quasi-Friedmann universe, i.e., the Friedmann universe with anomalous solid angle. The asymptotically proper Friedmann solution is possible only if the universe is accelerated or the potential is negative. If the potential is positive, the density perturbation in the asymptotically proper Friedmann solution rapidly falls off at spatial infinity, while the mass perturbation is compensated. In the asymptotically quasi-Friedmann solution, the density perturbation falls off only in proportion to the inverse square of the areal radius and the relative mass perturbation approaches a nonzero constant at spatial infinity. The present result shows that a necessary condition holds in order that a self-gravitating body grows self-similarly due to the constant accretion of quintessence in an accelerating universe.

Based on the asymptotic analysis of ordinary differential equations, we classify all spherically symmetric self-similar solutions to the Einstein equations which are asymptotically Friedmann at large distances and contain a perfect fluid with equation of state $p=(\gamma -1)\mu$ with $0<\gamma<2/3$. This corresponds to a ``dark energy'' fluid and the Friedmann solution is accelerated in this case due to anti-gravity. This extends the previous analysis of spherically symmetric self-similar solutions for fluids with positive pressure ($\gamma>1$). However, in the latter case there is an additional parameter associated with the weak discontinuity at the sonic point and the solutions are only asymptotically ``quasi-Friedmann'', in the sense that they exhibit an angle deficit at large distances. In the $0<\gamma<2/3$ case, there is no sonic point and there exists a one-parameter family of solutions which are {\it genuinely} asymptotically Friedmann at large distances. We find eight classes of asymptotic behavior: Friedmann or quasi-Friedmann or quasi-static or constant-velocity at large distances, quasi-Friedmann or positive-mass singular or negative-mass singular at small distances, and quasi-Kantowski-Sachs at intermediate distances. The self-similar asymptotically quasi-static and quasi-Kantowski-Sachs solutions are analytically extendible and of great cosmological interest. We also investigate their conformal diagrams. The results of the present analysis are utilized in an accompanying paper to obtain and physically interpret numerical solutions.

We use a combination of numerical and analytical methods, exploiting the equations derived in a preceding paper, to classify all spherically symmetric self-similar solutions which are asymptotically Friedmann at large distances and contain a perfect fluid with equation of state $p=(\gamma -1)\mu$ with $0<\gamma<2/3$. The expansion of the Friedmann universe is accelerated in this case. We find a one-parameter family of self-similar solutions representing a black hole embedded in a Friedmann background. This suggests that, in contrast to the positive pressure case, black holes in a universe with dark energy can grow as fast as the Hubble horizon if they are not too large. There are also self-similar solutions which contain a central naked singularity with negative mass and solutions which represent a Friedmann universe connected to either another Friedmann universe or some other cosmological model. The latter are interpreted as self-similar cosmological white hole or wormhole solutions. The throats of these wormholes are defined as two-dimensional spheres with minimal area on a spacelike hypersurface and they are all non-traversable because of the absence of a past null infinity.

We calculate the quantum radiation power of black holes which are asymptotic to the Einstein-de Sitter universe at spatial and null infinities. We consider two limiting mass accretion scenarios, no accretion and significant accretion. We find that the radiation power strongly depends on not only the asymptotic condition but also the mass accretion scenario. For the no accretion case, we consider the Einstein-Straus solution, where a black hole of constant mass resides in the dust Friedmann universe. We find negative cosmological correction besides the expected redshift factor. This is given in terms of the cubic root of ratio in size of the black hole to the cosmological horizon, so that it is currently of order $10^{-5} (M/10^{6}M_{\odot})^{1/3} (t/14 {Gyr})^{-1/3}$ but could have been significant at the formation epoch of primordial black holes. Due to the cosmological effects, this black hole has not settled down to an equilibrium state. This cosmological correction may be interpreted in an analogy with the radiation from a moving mirror in a flat spacetime. For the significant accretion case, we consider the Sultana-Dyer solution, where a black hole tends to increase its mass in proportion to the cosmological scale factor. In this model, we find that the radiation power is apparently the same as the Hawking radiation from the Schwarzschild black hole of which mass is that of the growing mass at each moment. Hence, the energy loss rate decreases and tends to vanish as time proceeds. Consequently, the energy loss due to evaporation is insignificant compared to huge mass accretion onto the black hole. Based on this model, we propose a definition of quasi-equilibrium temperature for general conformal stationary black holes.

Applying the first and generalised second laws of thermodynamics for a realistic process of near critical black hole formation, we derive an entropy bound, which is identical to Bekenstein's one for radiation. Relying upon this bound, we derive an absolute minimum mass $\sim0.04 \sqrt{g_{*}}m_{\rm Pl}$, where $g_{*}$ and $m_{\rm Pl}$ is the effective degrees of freedom for the initial temparature and the Planck mass, respectively. Since this minimum mass coincides with the lower bound on masses of which black holes can be regarded as classical against the Hawking evaporation, the thermodynamical argument will not prohibit the formation of the smallest classical black hole. For more general situations, we derive a minimum mass, which may depend on the initial value for entropy per particle. For primordial black holes, however, we show that this minimum mass can not be much greater than the Planck mass at any formation epoch of the Universe, as long as $g_{*}$ is within a reasonable range. We also derive a size-independent upper bound on the entropy density of a stiff fluid in terms of the energy density.

We consider the possible existence of self-similar solutions containing black holes in a Friedmann background with a stiff fluid or a scalar field. We carefully study the relationship between the self-similar equations in these two cases and emphasize the crucial role of the similarity horizon. We show that there is no self-similar black hole solution surrounded by an exact or asymptotically flat Friedmann background containing a massless scalar field. This result also applies for a scalar field with a potential, providing the universe is decelerating. However, if there is a potential and the universe is accelerating (as in the quintessence scenario), the result only applies for an exact Friedmann background. This extends the result previously found in the stiff fluid case and strongly suggests that accretion onto primordial black holes is ineffective even during scalar field domination. It also contradicts recent claims that such black holes can grow appreciably by accreting quintessence. Appreciable growth might be possible with very special matter fields but this requires {\it ad hoc} and probably unphysical conditions.

We discuss a new class of solutions to the Einstein equations which describe a primordial black hole (PBH) in a flat Friedmann background. Such solutions arise if a Schwarzschild black hole is patched onto a Friedmann background via a transition region. They are possible providing the black hole event horizon is larger than the cosmological apparent horizon. Such solutions have a number of strange features. In particular, one has to define the black hole and cosmological horizons carefully and one then finds that the mass contained within the black hole event horizon decreases when it is larger than the Friedmann cosmological apparent horizon, although its area always increases. These solutions involve two distinct future null infinities and are interpreted as the conversion of a white hole into a black hole. Although such solutions may not form from gravitational collapse in the same way as standard PBHs, there is nothing unphysical about them, since all energy and causality conditions are satisfied. Their conformal diagram is a natural amalgamation of the Kruskal diagram for the extended Schwarzschild solution and the conformal diagram for a black hole in a flat Friedmann background. In this paper, such solutions are obtained numerically for a spherically symmetric universe containing a massless scalar field, but it is likely that they exist for more general matter fields and less symmetric systems.

We provide precise constraints on the size of any black holes forming in the early Universe for a variety of formation scenarios. In particular, we prove that the size of the apparent horizon of a primordial black hole formed by causal processes in a flat Friedmann universe is considerably smaller than the cosmological apparent horizon size for an equation of state $p=k\rho$ ($1/3<k<1$). This also applies for a stiff equation of state ($k=1$) or for a massless scalar field. The apparent horizon of a primordial black hole formed through hydrodynamical processes is also considerably smaller than the cosmological apparent horizon for $0<k\le 1$. We derive an expression for the maximum size which an overdense region can have without being a separate closed universe rather than part of our own. Newtonian argument shows that a black hole smaller than the cosmological horizon can never accrete much.

The evolution of primordial black holes in a flat Friedmann universe with a massless scalar field is investigated in fully general relativistic numerical relativity. A primordial black hole is expected to form with a scale comparable to the cosmological apparent horizon, in which case it may go through an initial phase with significant accretion. However, if it is very close to the cosmological apparent horizon size, the accretion is suppressed due to general relativistic effects. In any case, it soon gets smaller than the cosmological horizon and thereafter it can be approximated as an isolated vacuum solution with decaying mass accretion. In this situation the dynamical and inhomogeneous scalar field is typically equivalent to a perfect fluid with a stiff equation of state $p=\rho$. The black hole mass never increases by more than a factor of two, despite recent claims that primordial black holes might grow substantially through accreting quintessence. It is found that the gravitational memory scenario, proposed for primordial black holes in Brans-Dicke and scalar-tensor theories of gravity, is highly unphysical.

There are many solutions to the Einstein field equations that demonstrate naked singularity (NS) formation after regular evolution. It is possible, however, that such a quantum effect as particle creation prevents NSs from forming. We investigate the relation between the curvature strength and the quantum effects of NSs in a very wide class of spherical dust collapse. Through a perturbative calculation, we find that if the NS is very strong, the quantum particle creation diverges as the Cauchy horizon is approached, while if the NS is very weak, the creation should be finite. In the context of cosmic censorship, strong NSs will be subjected to the backreaction of quantum effects and may disappear or be hidden behind horizons, while weak NSs will not.

The evolution of weak discontinuity is investigated in the flat FRW universe with a single scalar field and with multiple scalar fields. We consider both massless scalar fields and scalar fields with exponential potentials. Then we find that a new type of instability, i.e. kink instability develops in the flat FRW universe with massless scalar fields. The kink instability develops with scalar fields with considerably steep exponential potentials, while less steep exponential potentials do not suffer from kink instability. In particular, assisted inflation with multiple scalar fields does not suffer from kink instability. The stability of general spherically symmetric self-similar solutions is also discussed.

We present the results of indirect searches for Weakly Interacting Massive Particles (WIMPs) with 1679.6 live days of data from the Super-Kamiokande detector using neutrino-induced upward through-going muons. The search is performed by looking for an excess of high energy muon neutrinos from WIMP annihilations in the Sun, the core of the Earth, and the Galactic Center, as compared to the number expected from the atmospheric neutrino background. No statistically significant excess was seen. We calculate flux limits in various angular cones around each of the above celestial objects. We obtain conservative model-independent upper limits on WIMP-nucleon cross-section as a function of WIMP mass and compare these results with the corresponding results from direct dark matter detection experiments.

The gravitational interaction is scale-free in both Newtonian gravity and general theory of relativity. The concept of self-similarity arises from this nature. Self-similar solutions reproduce themselves as the scale changes. This property results in great simplification of the governing partial differential equations. In addition, some self-similar solutions can describe the asymptotic behaviors of more general solutions. Newtonian gravity contains only one dimensional constant, the gravitational constant, while the general relativity contains another dimensional constant, the speed of light, besides the gravitational constant. Due to this crucial difference, incomplete similarity can be more interesting in general relativity than in Newtonian gravity. Kinematic self-similarity has been defined and studied as an example of incomplete similarity in general relativity, in an effort to pursue a wider application of self-similarity in general relativity. We review the mathematical and physical aspects of kinematic self-similar solutions in general relativity.

We consider the possibilities for obtaining information about the equation of state for quark matter by using future direct observational data on gravitational waves. We study the nonradial oscillations of both fluid and spacetime modes of pure quark stars. If we observe the $f$ and the lowest $w_{\rm II}$ modes from quark stars, by using the simultaneously obtained radiation radius we can constrain the bag constant $B$ with reasonable accuracy, independently of the $s$ quark mass.

Generic spherically symmetric self-similar collapse results in strong naked-singularity formation. In this paper we are concerned with particle creation during a naked-singularity formation in spherically symmetric self-similar collapse without specifying the collapsing matter. In the generic case, the power of particle emission is found to be proportional to the inverse square of the remaining time to the Cauchy horizon (CH). The constant of proportion can be arbitrarily large in the limit to marginally naked singularity. Therefore, the unbounded power is especially striking in the case that an event horizon is very close to the CH because the emitted energy can be arbitrarily large in spite of a cutoff expected from quantum gravity. Above results suggest the instability of the CH in spherically symmetric self-similar spacetime from quantum field theory and seem to support the existence of a semiclassical cosmic censor. The divergence of redshifts and blueshifts of emitted particles is found to cause the divergence of power to positive or negative infinity, depending on the coupling manner of scalar fields to gravity. On the other hand, it is found that there is a special class of self-similar spacetimes in which the semiclassical instability of the CH is not efficient. The analyses in this paper are based on the geometric optics approximation, which is justified in two dimensions but needs justification in four dimensions.

A stability criterion is derived in general relativity for self-similar solutions with a scalar field and those with a stiff fluid, which is a perfect fluid with the equation of state $P=\rho$. A wide class of self-similar solutions turn out to be unstable against kink mode perturbation. According to the criterion, the Evans-Coleman stiff-fluid solution is unstable and cannot be a critical solution for the spherical collapse of a stiff fluid if we allow sufficiently small discontinuity in the density gradient field in the initial data sets. The self-similar scalar-field solution, which was recently found numerically by Brady {\it et al.} (2002 {\it Class. Quantum. Grav.} {\bf 19} 6359), is also unstable. Both the flat Friedmann universe with a scalar field and that with a stiff fluid suffer from kink instability at the particle horizon scale.

Nonspherical stellar collapse to a black hole is one of the most promising gravitational wave sources for gravitational wave detectors. We numerically study gravitational waves from a slightly nonspherical stellar collapse to a black hole in linearized Einstein theory. We adopt a spherically collapsing star as the zeroth-order solution and gravitational waves are computed using perturbation theory on the spherical background. In this paper we focus on the perturbation of odd-parity modes. Using the polytropic equations of state with polytropic indices $n_p=1$ and 3, we qualitatively study gravitational waves emitted during the collapse of neutron stars and supermassive stars to black holes from a marginally stable equilibrium configuration. Since the matter perturbation profiles can be chosen arbitrarily, we provide a few types for them. For $n_p=1$, the gravitational waveforms are mainly characterized by a black hole quasinormal mode ringing, irrespective of perturbation profiles given initially. However, for $n_p=3$, the waveforms depend strongly on the initial perturbation profiles. In other words, the gravitational waveforms strongly depend on the stellar configuration and, in turn, on the ad hoc choice of the functional form of the perturbation in the case of supermassive stars.

We study through numerical simulation the spherical collapse of isothermal gas in Newtonian gravity. We observe a critical behavior which occurs at the threshold of gravitational instability leading to core formation. For a given initial density profile, we find a critical temperature, which is of the same order as the virial temperature of the initial configuration. For the exact critical temperature, the collapse converges to a self-similar form, the first member in Hunter's family of self-similar solutions. For a temperature close to the critical value, the collapse first approaches this critical solution. Later on, in the supercritical case, the collapse converges to another self-similar solution, which is called the Larson-Penston solution. In the subcritical case, the gas bounces and disperses to infinity. We find two scaling laws: one for the collapsed mass in the supercritical case and the other for the maximum density reached before dispersal in the subcritical case. The value of the critical exponent is measured to be $\simeq 0.11$ in the supercritical case, which agrees well with the predicted value $\simeq 0.10567$. These critical properties are quite similar to those observed in the collapse of a radiation fluid in general relativity. We study the response of the system to temperature fluctuation and discuss astrophysical implications for the insterstellar medium structure and for the star formation process. Newtonian critical behavior is important not only because it provides a simple model for general relativity but also because it is relevant for astrophysical systems such as molecular clouds.

In cosmological models with a varying gravitational constant, it is not clear whether primordial black holes preserve the value of $G$ at their formation epoch. We investigate this question by using the Tolman-Bondi model to study the evolution of a background scalar field when a black hole forms from the collapse of dust in a flat Friedmann universe. Providing the back reaction of the scalar field on the metric can be neglected, we find that the value of the scalar field at the event horizon very quickly assumes the background cosmological value. This suggests that there is very little gravitational memory.

We investigate a particle velocity in the $\kappa$-Minkowski space-time, which is one of the realization of a noncommutative space-time. We emphasize that arrival time analyses by high-energy $\gamma$-rays or neutrinos, which have been considered as powerful tools to restrict the violation of Lorentz invariance, are not effective to detect space-time noncommutativity. In contrast with these examples, we point out a possibility that {\it low-energy massive particles} play an important role to detect it.

Gravitational collapse is one of the most fruitful subjects in gravitational physics. It is well known that singularity formation is inevitable in complete gravitational collapse. It was conjectured that such a singularity should be hidden by horizons if it is formed from generic initial data with physically reasonable matter fields. Many possible counterexamples to this conjecture have been proposed over the past three decades, although none of them has proved to be sufficiently generic. In these examples, there appears a singularity that is not hidden by horizons. This singularity is called a `naked singularity.' The appearance of a naked singularity represents the formation of an observable high-curvature, strong-gravity region. In this paper we review examples of naked singularity formation and recent progress in research of observable physical processes - gravitational radiation and quantum particle creation - from a forming naked singularity.

We discuss astrophysical implications of $\kappa$-Minkowski space-time, in which there appears space-time noncommutativity. We first derive a velocity formula for particles based on the motion of a wave packet. The result is that a massless particle moves at a constant speed as in the usual Minkowski space-time, which implies that an arrival time analysis by $\gamma$-rays from Markarian (Mk) 421 does not exclude space-time noncommutativity. Based on this observation, we analyze reaction processes in $\kappa$-Minkowski space-time which are related to the puzzling detections of extremely high-energy cosmic rays above the Greisen-Zatsepin-Kuzmin cutoff and of high-energy ($\sim$20 TeV) $\gamma$-rays from Mk 501.

A stability criterion is derived for self-similar solutions with perfect fluids which obey the equation of state $P=k\rho$ in general relativity. A wide class of self-similar solutions turn out to be unstable against the so-called kink mode. The criterion is directly related to the classification of sonic points. The criterion gives a sufficient condition for instability of the solution. For a transonic point in collapse, all primary-direction nodal-point solutions are unstable, while all secondary-direction nodal-point solutions and saddle-point ones are stable against the kink mode. The situation is reversed in expansion. Applications are the following: the expanding flat Friedmann solution for $1/3 \le k < 1$ and the collapsing one for $0< k \le 1/3$ are unstable; the static self-similar solution is unstable; nonanalytic self-similar collapse solutions are unstable; the Larson-Penston (attractor) solution is stable for this mode for $0<k\alt 0.036$, while it is unstable for $0.036\alt k $; the Evans-Coleman (critical) solution is stable for this mode for $0<k\alt 0.89$, while it is unstable for $0.89\alt k$. The last application suggests that the Evans-Coleman solution for $0.89\alt k $ is {\em not critical} because it has at least two unstable modes.

We investigate the stability of self-similar solutions for a gravitationally collapsing isothermal sphere in Newtonian gravity by means of a normal mode analysis. It is found that the Hunter series of solutions are highly unstable, while neither the Larson-Penston solution nor the homogeneous collapse one have an analytic unstable mode. Since the homogeneous collapse solution is known to suffer the kink instability, the present result and recent numerical simulations strongly support a proposition that the Larson-Penston solution will be realized in astrophysical situations. It is also found that the Hunter (A) solution has a single unstable mode, which implies that it is a critical solution associated with some critical phenomena which are analogous to those in general relativity. The critical exponent $\gamma$ is calculated as $\gamma\simeq 0.10567$. In contrast to the general relativistic case, the order parameter will be the collapsed mass. In order to obtain a complete picture of the Newtonian critical phenomena, full numerical simulations will be needed.

A direct method is developed to reconstruct the equation of state for high-density nuclear matter from the relationship between any two properties of neutron stars, such as masses, radii, moments of inertia, baryonic masses, binding energies, gravitational redshifts, and their combinations.

It has been pointed out that the value of the gravitational ``constant'' in the early universe may be different from that at present. In that case, it was conjectured that primordial black holes may ``remember'' the value of the gravitational constant in the early universe. The present analysis shows that this is not the case, at least in certain contexts.

We study the spherical collapse of a perfect fluid with an equation of state $P=k\rho$ by full general relativistic numerical simulations. For $0<k\alt 0.036$, it has been known that there exists a general relativistic counterpart of the Larson-Penston self-similar Newtonian solution. The numerical simulations strongly suggest that, in the neighborhood of the center, generic collapse converges to this solution in an approach to a singularity and that self-similar solutions other than this solution, including a ``critical solution'' in the black hole critical behavior, are relevant only when the parameters which parametrize initial data are fine-tuned. This result is supported by a mode analysis on the pertinent self-similar solutions. Since a naked singularity forms in the general relativistic Larson-Penston solution for $0<k\alt0.0105$, this will be the most serious known counterexample against cosmic censorship. It also provides strong evidence for the self-similarity hypothesis in general relativistic gravitational collapse. The direct consequence is that critical phenomena will be observed in the collapse of isothermal gas in Newton gravity, and the critical exponent $\gamma$ will be given by $\gamma\approx 0.11$, though the order parameter cannot be the black hole mass.

Naked singularity occurs in the gravitational collapse of an inhomogeneous dust ball from an initial density profile which is physically reasonable. We show that explosive radiation is emitted during the formation process of the naked singularity. The energy flux is proportional to $(t_{\rm CH}-t)^{-3/2}$ for a minimally coupled massless scalar field, while is proportional to $(t_{\rm CH}-t)^{-1}$ for a conformally coupled massless scalar field, where $t_{\rm CH}-t$ is the `remained time' until the distant observer could observe the singularity if the naked singularity was formed. As a consequence, the radiated energy grows unboundedly for both scalar fields. The amount of the power and the energy depends on parameters which characterize the initial density profile but do not depend on the gravitational mass of the cloud. In particular, there is characteristic frequency $\nu_{s}$ of singularity above which the divergent energy is radiated. The energy flux is dominated by particles of which the wave length is about $t_{\rm CH}-t$ at each moment. The observed total spectrum is nonthermal, i.e., $\nu dN/d\nu \sim (\nu/\nu_{s})^{-1}$ for $\nu>\nu_{s}$. If the naked singularity formation could continue until a considerable fraction of the total energy of the dust cloud is radiated, the radiated energy would reach about $10^{54}(M/M_{\odot})$ erg. The calculations are based on the geometrical optics approximation which turns out to be consistent as a rough order estimate. The analysis does not depend on whether or not the naked singularity occurs in its exact meaning. This phenomenon may provide a new candidate for a source of ultra high energy cosmic rays or a central engine of gamma ray bursts.

It is known that the gravitational collapse of a dust ball results in naked singularity formation from an initial density profile which is physically reasonable. In this paper, we show that explosive radiation is emitted during the formation process of the naked singularity.

We clarify the causal structure of an inflating magnetic monopole. The spacetime diagram shows explicitly that this model is free from ``graceful exit'' problem, while the monopole itself undergoes ``eternal inflation''. We also discuss general nature of inflationary spacetimes.

It was recently shown that the metric functions which describe a spherically symmetric space-time with vanishing radial pressure can be explicitly integrated. We investigate the nakedness and curvature strength of the shell-focusing singularity in that space-time. If the singularity is naked, the relation between the circumferential radius and the Misner-Sharp mass is given by $R\approx 2y_{0} m^{\beta}$ with $ 1/3<\beta\le 1$ along the first radial null geodesic from the singularity. The $\beta$ is closely related to the curvature strength of the naked singularity. For example, for the outgoing or ingoing null geodesic, if the strong curvature condition (SCC) by Tipler holds, then $\beta$ must be equal to 1. We define the ``gravity dominance condition'' (GDC) for a geodesic. If GDC is satisfied for the null geodesic, both SCC and the limiting focusing condition (LFC) by Kr\'olak hold for $\beta=1$ and $y_{0}\ne 1$, not SCC but only LFC holds for $1/2\le \beta <1$, and neither holds for $1/3<\beta <1/2$, for the null geodesic. On the other hand, if GDC is satisfied for the timelike geodesic $r=0$, both SCC and LFC are satisfied for the timelike geodesic, irrespective of the value of $\beta$. Several examples are also discussed.

The final fate of the spherically symmetric collapse of a perfect fluid which follows the $\gamma$-law equation of state and adiabatic condition is investigated. Full general relativistic hydrodynamics is solved numerically using a retarded time coordinate, the so-called observer time coordinate. Thanks to this coordinate, the causal structure of the resultant space-time is automatically constructed. Then, it is found that a globally naked, shell-focusing singularity can occur at the center from relativistically high-density, isentropic and time symmetric initial data if $\gamma \alt 1.01$ within the numerical accuracy. The result is free from the assumption of self-similarity. The upper limit of $\gamma$ with which a naked singularity can occur from generic initial data is consistent with the result of Ori and Piran based on the assumption of self-similarity.

We investigate collapse of a spherical cloud of counter-rotating particles. An explicit solution is given using an elliptic integral. If the specific angular momentum $L(r)=O(r^2)$ at $r\to 0$, no central singularity occurs. With $L(r)$ like that, there is a finite region around the center that bounces. On the other hand, if the order of $L(r)$ is higher than that, a central singularity occurs. In marginally bound collapse with $L(r)=4F(r)$, a naked singularity occurs, where $F(r)$ is the Misner-Sharp mass. The solution for this case is expressed by elementary functions. For $ 4 <L/F<\infty$ at $r\to0$, there is a finite region around the center that bounces and a naked singularity occurs. For $ 0 \le L/F< 4$ at $r\to0$, there is no such region. The results suggests that rotation may play a crucial role on the final fate of collapse.

We investigate neutron stars in scalar-tensor theories. We examine their secular stability against spherically symmetric perturbations by use of a turning point method. For some choices of the coupling function contained in the theories, the number of the stable equilibrium solutions changes and the realized equilibrium may change discontinuously as the asymptotic value of the scalar field or total baryon number is changed continuously. The behaviour of the stable equilibrium solutions is explained by fold and cusp catastrophes. Whether the cusp catastrophe appears or not depends on the choices of the coupling function. These types of the catastrophes are structurally stable. Recently discovered spontaneous scalarization, which is non-perturbative strong-field phenomenon due to the presence of the gravitational scalar field, is well described in terms of the cusp catastrophe.

A stability analysis of a spherically symmetric star in scalar-tensor theories of gravity is given in terms of the frequencies of quasi-normal modes. The scalar-tensor theories have a scalar field which is related to gravitation. There is an arbitrary function, the so-called coupling function, which determines the strength of the coupling between the gravitational scalar field and matter. Instability is induced by the scalar field for some ranges of the value of the first derivative of the coupling function. This instability leads to significant discrepancies with the results of binary-pulsar-timing experiments and hence, by the stability analysis, we can exclude the ranges of the first derivative of the coupling function in which the instability sets in. In this article, the constraint on the first derivative of the coupling function from the stability of relativistic stars is found. Analysis in terms of the quasi-normal mode frequencies accounts for the parameter dependence of the wave form of the scalar gravitational waves emitted from the Oppenheimer-Snyder collapse. The spontaneous scalarization is also discussed.

Unlike general relativity, scalar-tensor theories of gravity predict scalar gravitational waves even from a spherically symmetric gravitational collapse. We solve numerically the generation and propagation of the scalar gravitational wave from a spherically symmetric and homogeneous dust collapse under the approximation that we can neglect the back reaction of the scalar wave on the space-time, and examine how the amplitude, characteristic frequency and wave form of the observed scalar gravitational wave depend on the initial radius and mass of the dust and parameters contained in the theory. In the Brans-Dicke theory, through the observation of the scalar gravitational wave, it is possible to determine the initial radius and mass and a parameter contained in the theory. In the scalar-tensor theories, it would be possible to get the information of the first derivative of the coupling function contained in the theory because the wave form of the scalar gravitational wave greatly depends on it.