Critical phenomena at the threshold of black hole formation for collisionless matter in spherical symmetry

We perform a numerical study of the critical regime at the threshold of black hole formation in the spherically symmetric, general relativistic collapse of collisionless matter. The coupled Einstein-Vlasov equations are solved using a particle-mesh method in which the evolution of the phase-space distribution function is approximated by a set of particles (or, more precisely, infinitesimally thin shells) moving along geodesics of the spacetime. Individual particles may have non-zero angular momenta, but spherical symmetry dictates that the total angular momentum of the matter distribution vanish. In accord with previous work by Rein et al, our results indicate that the critical behavior in this model is Type I; that is, the smallest black hole in each parametrized family has a finite mass. We present evidence that the critical solutions are characterized by unstable, static spacetimes, with non-trivial distributions of radial momenta for the particles. As expected for Type I solutions, we also find power-law scaling relations for the lifetimes of near-critical configurations as a function of parameter-space distance from criticality.Comment: 32 pages, 10 figure

Numerical evidence for `multi-scalar stars'

We present a class of general relativistic soliton-like solutions composed of multiple minimally coupled, massive, real scalar fields which interact only through the gravitational field. We describe a two-parameter family of solutions we call ``phase-shifted boson stars'' (parameterized by central density rho_0 and phase delta), which are obtained by solving the ordinary differential equations associated with boson stars and then altering the phase between the real and imaginary parts of the field. These solutions are similar to boson stars as well as the oscillating soliton stars found by Seidel and Suen [E. Seidel and W.M. Suen, Phys. Rev. Lett. 66, 1659 (1991)]; in particular, long-time numerical evolutions suggest that phase-shifted boson stars are stable. Our results indicate that scalar soliton-like solutions are perhaps more generic than has been previously thought.Comment: Revtex. 4 pages with 4 figures. Submitted to Phys. Rev.

Evolution of the Schr\"odinger--Newton system for a self--gravitating scalar field

Using numerical techniques, we study the collapse of a scalar field configuration in the Newtonian limit of the spherically symmetric Einstein--Klein--Gordon (EKG) system, which results in the so called Schr\"odinger--Newton (SN) set of equations. We present the numerical code developed to evolve the SN system and topics related, like equilibrium configurations and boundary conditions. Also, we analyze the evolution of different initial configurations and the physical quantities associated to them. In particular, we readdress the issue of the gravitational cooling mechanism for Newtonian systems and find that all systems settle down onto a 0--node equilibrium configuration.Comment: RevTex file, 19 pages, 26 eps figures. Minor changes, matches version to appear in PR

Bondian frames to couple matter with radiation

A study is presented for the non linear evolution of a self gravitating distribution of matter coupled to a massless scalar field. The characteristic formulation for numerical relativity is used to follow the evolution by a sequence of light cones open to the future. Bondian frames are used to endow physical meaning to the matter variables and to the massless scalar field. Asymptotic approaches to the origin and to infinity are achieved; at the boundary surface interior and exterior solutions are matched guaranteeing the Darmois--Lichnerowicz conditions. To show how the scheme works some numerical models are discussed. We exemplify evolving scalar waves on the following fixed backgrounds: A) an atmosphere between the boundary surface of an incompressible mixtured fluid and infinity; B) a polytropic distribution matched to a Schwarzschild exterior; C) a Schwarzschild- Schwarzschild spacetime. The conservation of energy, the Newman--Penrose constant preservation and other expected features are observed.Comment: 20 pages, 6 figures; to appear in General Relativity and Gravitatio

Three-dimensional general relativistic hydrodynamics II: long-term dynamics of single relativistic stars

This is the second in a series of papers on the construction and validation of a three-dimensional code for the solution of the coupled system of the Einstein equations and of the general relativistic hydrodynamic equations, and on the application of this code to problems in general relativistic astrophysics. In particular, we report on the accuracy of our code in the long-term dynamical evolution of relativistic stars and on some new physics results obtained in the process of code testing. The tests involve single non-rotating stars in stable equilibrium, non-rotating stars undergoing radial and quadrupolar oscillations, non-rotating stars on the unstable branch of the equilibrium configurations migrating to the stable branch, non-rotating stars undergoing gravitational collapse to a black hole, and rapidly rotating stars in stable equilibrium and undergoing quasi-radial oscillations. The numerical evolutions have been carried out in full general relativity using different types of polytropic equations of state using either the rest-mass density only, or the rest-mass density and the internal energy as independent variables. New variants of the spacetime evolution and new high resolution shock capturing (HRSC) treatments based on Riemann solvers and slope limiters have been implemented and the results compared with those obtained from previous methods. Finally, we have obtained the first eigenfrequencies of rotating stars in full general relativity and rapid rotation. A long standing problem, such frequencies have not been obtained by other methods. Overall, and to the best of our knowledge, the results presented in this paper represent the most accurate long-term three-dimensional evolutions of relativistic stars available to date.Comment: 19 pages, 17 figure

Classical and Quantum Decay of Oscillatons: Oscillating Self-Gravitating Real Scalar Field Solitons

The oscillating gravitational field of an oscillaton of finite mass M causes it to lose energy by emitting classical scalar field waves, but at a rate that is non-perturbatively tiny for small GMm, where m is the scalar field mass: d(GM)/dt ~ -3797437.776333015 e^[-39.433795197160163/(GMm)]/(GMm)^2. Oscillatons also decay by the quantum process of the annihilation of scalarons into gravitons, which is only perturbatively small in GMm, giving by itself d(GM)/dt ~ - 0.008513223934732692 G m^2 (GMm)^5. Thus the quantum decay is faster than the classical one for Gmm < 39.4338/[ln(1/Gm^2)}-7ln(GMm)+19.9160]. The time for an oscillaton to decay away completely into free scalarons and gravitons is ~ 2/(G^5 m^11) ~ 10^324 yr (1 meV/m)^11. Oscillatons of more than one real scalar field of the same mass generically asymptotically approach a static-geometry U(1) boson star configuration with GMm = GM_0 m, at the rate d(GM/c^3)/dt ~ [(C/(GMm)^4)e^{-alpha/(GMm)}+Q(m/m_{Pl})^2(GMm)^3] [(GMm)^2-(GM_0 m)^2], with GM_0 m depending on the magnitudes and relative phases of the oscillating fields, and with the same constants C, alpha, and Q given numerically above for the single-field case that is equivalent to GM_0 m = 0.Comment: 75 pages, LaTe

General Relativistic MHD Jets

Magnetic fields connecting the immediate environs of rotating black holes to large distances appear to be the most promising mechanism for launching relativistic jets, an idea first developed by Blandford and Znajek in the mid-1970s. To enable an understanding of this process, we provide a brief introduction to dynamics and electromagnetism in the spacetime near black holes. We then present a brief summary of the classical Blandford-Znajek mechanism and its conceptual foundations. Recently, it has become possible to study these effects in much greater detail using numerical simulation. After discussing which aspects of the problem can be handled well by numerical means and which aspects remain beyond the grasp of such methods, we summarize their results so far. Simulations have confirmed that processes akin to the classical Blandford-Znajek mechanism can launch powerful electromagnetically-dominated jets, and have shown how the jet luminosity can be related to black hole spin and concurrent accretion rate. However, they have also shown that the luminosity and variability of jets can depend strongly on magnetic field geometry. We close with a discussion of several important open questions.Comment: 21 pages, 2 figures, To appear in Belloni, T. (ed.): The Jet Paradigm - From Microquasars to Quasars, Lect. Notes Phys. 794 (2009

Theory of magnetically powered jets

The magnetic theory for the production of jets by accreting objects is reviewed with emphasis on outstanding problem areas. An effort is made to show the connections behind the occasionally diverging nomenclature in the literature, to contrast the different points of view about basic mechanisms, and to highlight concepts for interpreting the results of numerical simulations. The role of dissipation of magnetic energy in accelerating the flow is discussed, and its importance for explaining high Lorentz factors. The collimation of jets to the observed narrow angles is discussed, including a critical discussion of the role of `hoop stress'. The transition between disk and outflow is one of the least understood parts of the magnetic theory; its role in setting the mass flux in the wind, in possible modulations of the mass flux, and the uncertainties in treating it realistically are discussed. Current views on most of these problems are still strongly influenced by the restriction to 2 dimensions (axisymmetry) in previous analytical and numerical work; 3-D effects likely to be important are suggested. An interesting problem area is the nature and origin of the strong, preferably highly ordered magnetic fields known to work best for jet production. The observational evidence for such fields and their behavior in numerical simulations is discussed. I argue that the presence or absence of such fields may well be the `second parameter' governing not only the presence of jets but also the X-ray spectra and timing behavior of X-ray binaries.Comment: 29 pages. Publication delays offered the opportunity for further corrections, an expansion of sect 4.2, and one more Fig. To appear in Belloni, T. (ed.): The Jet Paradigm - From Microquasars to Quasars, Lect. Notes Phys. 794 (2009