Adaptive Mesh Refinement for Coupled Elliptic-Hyperbolic Systems

We present a modification to the Berger and Oliger adaptive mesh refinement algorithm designed to solve systems of coupled, non-linear, hyperbolic and elliptic partial differential equations. Such systems typically arise during constrained evolution of the field equations of general relativity. The novel aspect of this algorithm is a technique of "extrapolation and delayed solution" used to deal with the non-local nature of the solution of the elliptic equations, driven by dynamical sources, within the usual Berger and Oliger time-stepping framework. We show empirical results demonstrating the effectiveness of this technique in axisymmetric gravitational collapse simulations. We also describe several other details of the code, including truncation error estimation using a self-shadow hierarchy, and the refinement-boundary interpolation operators that are used to help suppress spurious high-frequency solution components ("noise").Comment: 31 pages, 15 figures; replaced with published versio

Gravitational collapse in anti de Sitter space

A numerical and analytic treatment is presented here of the evolution of initial data of the kind that was conjectured by Hertog, Horowitz and Maeda to lead to a violation of cosmic censorship. That initial data is essentially a thick domain wall connecting two regions of anti de Sitter space. The evolution results in no violation of cosmic censorship, but rather the formation of a small black hole.Comment: 9 pages, 13 figure

Driven neutron star collapse: Type~I critical phenomena and the initial black hole mass distribution

We study the general relativistic collapse of neutron star (NS) models in spherical symmetry. Our initially stable models are driven to collapse by the addition of one of two things: an initially in-going velocity profile, or a shell of minimally coupled, massless scalar field that falls onto the star. Tolman-Oppenheimer-Volkoff (TOV) solutions with an initially isentropic, gamma-law equation of state serve as our NS models. The initial values of the velocity profile's amplitude and the star's central density span a parameter space which we have surveyed extensively and which we find provides a rich picture of the possible end states of NS collapse. This parameter space survey elucidates the boundary between Type I and Type II critical behavior in perfect fluids which coincides, on the subcritical side, with the boundary between dispersed and bound end states. For our particular model, initial velocity amplitudes greater than 0.3c are needed to probe the regime where arbitrarily small black holes can form. In addition, we investigate Type I behavior in our system by varying the initial amplitude of the initially imploding scalar field. In this case we find that the Type I critical solutions resemble TOV solutions on the 1-mode unstable branch of equilibrium solutions, and that the critical solutions' frequencies agree well with the fundamental mode frequencies of the unstable equilibria. Additionally, the critical solution's scaling exponent is shown to be well approximated by a linear function of the initial star's central density.Comment: Submitted to Phys. Rev. D., 24 pages, 25 monochrome figures. arXiv admin note: substantial text overlap with arXiv:gr-qc/031011

Black Hole Criticality in the Brans-Dicke Model

We study the collapse of a free scalar field in the Brans-Dicke model of gravity. At the critical point of black hole formation, the model admits two distinctive solutions dependent on the value of the coupling parameter. We find one solution to be discretely self-similar and the other to exhibit continuous self-similarity.Comment: 4 pages, REVTeX 3.0, 5 figures include

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

Nonminimally coupled topological-defect boson stars: Static solutions

We consider spherically symmetric static composite structures consisting of a boson star and a global monopole, minimally or non-minimally coupled to the general relativistic gravitational field. In the non-minimally coupled case, Marunovic and Murkovic have shown that these objects, so-called boson D-stars, can be sufficiently gravitationally compact so as to potentially mimic black holes. Here, we present the results of an extensive numerical parameter space survey which reveals additional new and unexpected phenomenology in the model. In particular, focusing on families of boson D-stars which are parameterized by the central amplitude of the boson field, we find configurations for both the minimally and non-minimally coupled cases that contain one or more shells of bosonic matter located far from the origin. In parameter space, each shell spontaneously appears as one tunes through some critical central amplitude of the boson field. In some cases the shells apparently materialize at spatial infinity: in these instances their areal radii are observed to obey a universal scaling law in the vicinity of the critical amplitude. We derive this law from the equations of motion and the asymptotic behavior of the fields.Comment: 17 pages, 24 figure

Non-Universality of Critical Behaviour in Spherically Symmetric Gravitational Collapse

The aim of the present letter is to explain the `critical behaviour' observed in numerical studies of spherically symmetric gravitational collaps of a perfect fluid. A simple expression results for the critical index $\gamma$ of the black hole mass considered as an order parameter. $\gamma$ turns out to vary strongly with the parameter $k$ of the assumed equation of state $p=k\rho$.Comment: 6