2.3 Nature of the singularity 2 Singularities in AF Spacetimes2.1 Naked singularities and the

2.2 Critical behavior in collapse

2.2.1 Gravitational collapse simulations

For a more detailed discussion of critical behavior see [126Jump To The Next Citation Point In The Article]. Since Gundlach's Living Review article has appeared, the updates in this section will be restricted to results I find especially interesting.

Critical behavior was originally found by Choptuik [78Jump To The Next Citation Point In The Article] in a numerical study of the collapse of a spherically symmetric massless scalar field. For recent reviews see [125, 130]. We note that this is the first completely new phenomenon in general relativity to be discovered by numerical simulation. In collapse of a scalar field, essentially two things can happen: Either a black hole (BH) forms or the scalar waves pass through each other and disperse. Choptuik discovered that for any 1-parameter set of initial data labeled by p, there is a critical value tex2html_wrap_inline1836 such that tex2html_wrap_inline1838 yields a BH. He found

  equation113

where tex2html_wrap_inline1840 is the mass of the eventual BH. The constant tex2html_wrap_inline1842 depends on the parameter of the initial data that is selected but tex2html_wrap_inline1844 is the same for all choices. Furthermore, in terms of logarithmic variables tex2html_wrap_inline1846, tex2html_wrap_inline1848 (tex2html_wrap_inline1850 is the proper time of an observer at r = 0, where r is the radial coordinate, with tex2html_wrap_inline1856 the finite proper time at which the critical evolution concludes, and tex2html_wrap_inline1858 is a constant which scales r), the waveform X repeats (echoes) at intervals of tex2html_wrap_inline1864 in tex2html_wrap_inline1866 if tex2html_wrap_inline1868 is rescaled to tex2html_wrap_inline1870, i.e. tex2html_wrap_inline1872 . The scaling behavior (1Popup Equation) demonstrates that the minimum BH mass (for bosons) is zero. The critical solution itself is a counter-example to cosmic censorship (since the formation of the zero mass BH causes high curvature regions to become visible at tex2html_wrap_inline1874). (See, e.g., the discussion in Hirschmann and Eardley [145Jump To The Next Citation Point In The Article].) The numerical demonstration of this feature of the critical solution was provided by Hamadé and Stewart [135Jump To The Next Citation Point In The Article]. This result caused Hawking to pay off a bet to Preskill and Thorne [61, 170].

 

Click on thumbnail to view image

Figure 3: This figure is the final frame of an animation of Type II critical behavior in Einstein-Yang-Mills collapse. Note the echoing in the near-critical solution. For the entire movie and related references see [77].

Soon after this discovery, scaling and critical phenomena were found in a variety of contexts. Abrahams and Evans [1Jump To The Next Citation Point In The Article] discovered the same phenomenon in axisymmetric gravitational wave collapse with a different value of tex2html_wrap_inline1864 and, to within numerical error, the same value of tex2html_wrap_inline1878 . (Note that the rescaling of r with tex2html_wrap_inline1882 required Choptuik to use adaptive mesh refinement (AMR) to distinguish subsequent echoes. Abrahams and Evans' smaller tex2html_wrap_inline1864 (tex2html_wrap_inline1886) allowed them to see echoing with their 2+1 code without AMR.) Garfinkle [109] confirmed Choptuik's results with a completely different algorithm that does not require AMR. He used Goldwirth and Piran's [119Jump To The Next Citation Point In The Article] method of simulating Christodoulou's [80] formulation of the spherically symmetric scalar field in null coordinates. This formulation allowed the grid to be automatically rescaled by choosing the edge of the grid to be the null ray that just hits the central observer at the end of the critical evolution. (Missing points of null rays that cross the central observer's world line are replaced by interpolation between those that remain.) Hamadé and Stewart [135] have also repeated Choptuik's calculation using null coordinates and AMR. They are able to achieve greater accuracy and find tex2html_wrap_inline1888 .

2.2.2 Critical solutions as an eigenvalue problem

Evans and Coleman [98] realized that self-similar rather than self-periodic collapse might be more tractable both numerically (since ODE's rather than PDE's are involved) and analytically. They discovered that a collapsing radiation fluid had that desirable property. (Note that self-similarity (homothetic motion) is incompatible with AF [93, 108]. However, most of the action occurs in the center so that a match of the self-similar inner region to an outer AF one should always be possible.) In a series of papers, Hirschmann and Eardley [144, 145] developed a (numerical) self-similar solution to the spherically symmetric complex scalar field equations. These are ODE's with too many boundary conditions causing a solution to exist only for certain fixed values of tex2html_wrap_inline1864 . Numerical solution of this eigenvalue problem allows very accurate determination of tex2html_wrap_inline1864 . The self-similarity also allows accurate calculation of tex2html_wrap_inline1878 as follows: The critical tex2html_wrap_inline1896 solution is unstable to a small change in p . At any time t (where t < 0 is increasing toward zero), the amplitude a of the perturbation exhibits power law growth:

  equation148

where tex2html_wrap_inline1906 . At any fixed t, larger a implies larger tex2html_wrap_inline1840 . Equivalently, any fixed amplitude tex2html_wrap_inline1914 will be reached faster for larger eventual tex2html_wrap_inline1840 . Scaling arguments give the dependence of tex2html_wrap_inline1840 on the time at which any fixed amplitude is reached:

  equation156

where

  equation160

Thus

  equation165

Therefore, one need only identify the growth rate of the unstable mode to obtain an accurate value of tex2html_wrap_inline1920 . It is not necessary to undertake the entire dynamical evolution or probe the space of initial data. (While this procedure allowed Hirschmann and Eardley to obtain tex2html_wrap_inline1922 for the complex scalar field solution, they later found [146Jump To The Next Citation Point In The Article] that, in this regime, the complex scalar field has 3 unstable modes. This means [129Jump To The Next Citation Point In The Article, 126] that the Eardley-Hirschmann solution is not a critical solution. A perturbation analysis indicates that the critical solution for the complex scalar field is the discretely self-similar one found for the real scalar field [129Jump To The Next Citation Point In The Article].) Koike et al. [176Jump To The Next Citation Point In The Article] obtain tex2html_wrap_inline1924 for the Evans-Coleman solution. Although the similarities among the critical exponents tex2html_wrap_inline1878 in the collapse computations suggested a universal value, Maison [183] used these same scaling-perturbation methods to show that tex2html_wrap_inline1878 depends on the equation of state tex2html_wrap_inline1930 of the fluid in the Evans-Coleman solution. Gundlach [127] used a similar approach to locate Choptuik's critical solution accurately. This is much harder due to its discrete self-similarity. He reformulates the model as nonlinear hyperbolic boundary value problem with eigenvalue tex2html_wrap_inline1864 and finds tex2html_wrap_inline1934 . As with the self-similar solutions described above, the critical solution is found directly without the need to perform a dynamical evolution or explore the space of initial data. Hara et al. extended the renormalization group approach of [176] and applied it to the continuously-self-similar case [136Jump To The Next Citation Point In The Article]. (For an application of renormalization group methods to cosmology see [162].)

2.2.3 Recent results

Gundlach [129] completed his eigenvalue analysis of the Choptuik solution to find the growth rate of the unstable mode to be tex2html_wrap_inline1936 . He also predicted a periodic ``wiggle'' in the Choptuik mass scaling relation. This was later observed numerically by Hod and Piran [152]. Self-similar critical behavior has been seen in string theory related axion-dilaton models [95, 134] and in the nonlinear tex2html_wrap_inline1938 -model [146Jump To The Next Citation Point In The Article]. Garfinkle and Duncan have shown that subcritical collapse of a spherically symmetric scalar field yields a scaling relation for the maximum curvature observed by the central observer with critical parameters that would be expected on the basis of those found for supercritical collapse [113].

Choptuik et al. [79] have generalized the original Einstein-scalar field calculations to the Einstein-Yang-Mills (EYM) (for SU (2)) case. Here something new was found. Two types of behavior appeared depending on the initial data. In Type I, BH formation had a non-zero mass threshold. The critical solution is a regular, unstable solution to the EYM equations found previously by Bartnik and McKinnon [14]. In Type II collapse, the minimum BH mass is zero with the critical solution similar to that of Choptuik (with a different tex2html_wrap_inline1942, tex2html_wrap_inline1944). Gundlach has also looked at this case with the same results [128]. The Type I behavior arises when the collapsing system has a metastable static solution in addition to the Choptuik critical one [132Jump To The Next Citation Point In The Article].

Brady, Chambers and Gonçalves [71, 52] conjectured that addition of a mass to the scalar field of the original model would break scale invariance and might yield a distinct critical behavior. They found numerically the same Type I and II ``phases'' seen in the EYM case. The Type II solution can be understood as perturbations of Choptuik's original model with a small scalar field mass tex2html_wrap_inline1946 . Here small means that tex2html_wrap_inline1948 where tex2html_wrap_inline1950 is the spatial extent of the original nonzero field region. (The scalar field is well within the Compton wavelength corresponding to tex2html_wrap_inline1946 .) On the other hand, tex2html_wrap_inline1954 yields Type I behavior. The minimum mass critical solution is an unstable soliton of the type found by Seidel and Suen [229]. The massive scalar field can be treated as collapsing dust to yield a criterion for BH formation [120].

The Choptuik solution has also been found to be the critical solution for charged scalar fields [132, 151]. As tex2html_wrap_inline1956, tex2html_wrap_inline1958 for the black hole. Q obeys a power law scaling. Numerical study of the critical collapse of collisionless matter (Einstein-Vlasov equations) has yielded a non-zero minimum BH mass [217, 201]. Bizon and Chmaj [47] have considered the critical collapse of skyrmions.

An astrophysical application of BH critical phenomena has been considered by Niemeyer and Jedamzik [199] and Yokoyama [252]. They consider its implications for primordial BH formation and suggest that it could be important.

2.2.4 Going further

The question is then why these critical phenomena should appear in so many collapsing gravitational systems. The discrete self-similarity of Choptuik's solution may be understood as scaling of a limit cycle [136]. (The self-similarity of other systems may be understood as scaling of a limit point.) Garfinkle [110] originally conjectured that the scale invariance of Einstein's equations might provide an underlying explanation for the self-similarity and discrete-self-similarity found in collapse and proposed a spacetime slicing which might manifestly show this. In fact, he later showed (with Meyer) [116] that, while the originally proposed slicing failed, a foliation based on maximal slicing did make the scaling manifest. These ideas formed the basis of a much more ambitious program by Garfinkle and Gundlach to use underlying actual or approximate symmetries to construct coordinate systems for numerical simulations [115].

An interesting ``toy model'' for general relativity in many contexts has been wave maps, also known as nonlinear tex2html_wrap_inline1938 models. One of these contexts is critical collapse [146]. Recently and independently, Bizon et al. [48] and Liebling et al. [180] evolved wave maps from the base space of 3+1 Minkowski space to the target space tex2html_wrap_inline1964 . They found critical behavior separating singular and nonsingular solutions. For some families of initial data, the critical solution is self-similar and is an intermediate attractor. For others, a static solution separates the singular and nonsingular solutions. However, the static solution has several unstable modes and is therefore not a critical solution in the usual sense. Bizon and Tabor [49] have studied Yang-Mills fields in D + 1 dimensions and found that generic solutions with sufficiently large initial data blow up in a finite time and that the mechanism for blowup depends on D . Husa et al. [158] then considered the collapse of SU (2) nonlinear sigma models coupled to gravity and found a discretely self-similar critical solution for sufficiently large dimensionless coupling constant. They also observe that for sufficiently small coupling constant values, there is a continuously self-similar solution. Interestingly, there is an intermediate range of coupling constant where the discrete self-similarity is intermittent [238].

Until recently, only Abrahams and Evans [1] had ventured beyond spherical symmetry. The first additional departure has been made by Gundlach [131]. He considered spherical and non-spherical perturbations of tex2html_wrap_inline1972 perfect fluid collapse. Only the original (spherical) growing mode survived.

Recently, critical phenomena have been explored in 2+1 gravity. Pretorius and Choptuik [215] numerically evolved circularly symmetric scalar field collapse in 2+1 anti-de Sitter space. They found a continuously self-similar critical solution at the threshold for black hole formation. The BH's which form have BTZ [8] exteriors with strong curvature, spacelike singularities in the interior. Remarkably, Garfinkle obtained an analytic critical solution by assuming continuous self-similarity which agrees quite well with the one obtained numerically [112].



2.3 Nature of the singularity 2 Singularities in AF Spacetimes2.1 Naked singularities and the

image Numerical Approaches to Spacetime Singularities
Beverly K. Berger
http://www.livingreviews.org/lrr-2002-1
© Max-Planck-Gesellschaft. ISSN 1433-8351
Problems/Comments to livrev@aei-potsdam.mpg.de