Axions and axion-like particles are a leading model for the dark matter in the Universe; therefore, dark matter halos may be boson stars in the process of collapsing. We examine a class of static boson stars with a non-minimal coupling to gravity. We modify the gravitational density of the boson field to be proportional to an arbitrary power of the modulus of the field, introducing a non-standard coupling. We find a class of solutions very similar to Newtonian polytropic stars that we denote "quantum polytropes." These quantum polytropes are supported by a non-local quantum pressure and follow an equation very similar to the Lane-Emden equation for classical polytropes. Furthermore, we derive a simple condition on the exponent of the non-linear gravitational coupling, $\alpha>8/3$, beyond which the equilibrium solutions are unstable.

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.

This article is an overview of the contributions numerical relativity has made to our understanding of strong field gravity, to be published in the book "General Relativity and Gravitation: A Centennial Perspective", commemorating the 100th anniversary of general relativity.

We develop a numerical hydrodynamics code using a pseudo-Newtonian formulation that uses the weak field approximation for the geometry, and a generalized source term for the Poisson equation that takes into account relativistic effects. The code was designed to treat moderately relativistic systems such as rapidly rotating neutron stars. The hydrodynamic equations are solved using a finite volume method with High Resolution Shock Capturing (HRSC) techniques. We implement several different slope limiters for second order reconstruction schemes and also investigate higher order reconstructions. We use the method of lines (MoL) to convert the mixed spatial-time partial differential equations into ordinary differential equations (ODEs) that depend only on time. These ODEs are solved using 2nd and 3rd order Runge-Kutta methods. The Poisson equation for the gravitational potential is solved with a multigrid method. In order to confirm the validity of our code, we carry out four different tests including one and two dimensional shock tube tests, stationary star tests of both non-rotating and rotating models and radial oscillation mode tests for spherical stars. In the shock tube tests, the code shows good agreement with analytic solutions which include shocks, rarefaction waves and contact discontinuities. The code is found to be stable and accurate: for example, when solving a stationary stellar model the fractional changes in the maximum density, total mass, and total angular momentum per dynamical time are found to be $3 \times 10^{-6}$, $5 \times 10^{-7}$ and $2 \times 10^{-6}$, respectively. We also find that the frequencies of the radial modes obtained by the numerical simulation of the steady state star agree very well with those obtained by linear analysis.

We present results from numerical solution of the Einstein field equations describing the head-on collision of two solitons boosted to ultra relativistic energies. We show, for the first time, that at sufficiently high energies the collision leads to black hole formation, consistent with hoop conjecture arguments. This implies that the non-linear gravitational interaction between the kinetic energy of the solitons causes gravitational collapse, and that arguments for black hole formation in super-Planck scale particle collisions are robust.

We investigate spherically-symmetric, general relativistic systems of collapsing perfect fluid distributions. We consider neutron star models that are driven to collapse by the addition of an initially "in-going" velocity profile to the nominally static star solution. The neutron star models we use are Tolman-Oppenheimer-Volkoff solutions with an initially isentropic, gamma-law equation of state. The initial values of 1) the amplitude of the velocity profile, and 2) the central density of the star, span a parameter space, and we focus only on that region that gives rise to Type II critical behavior, wherein black holes of arbitrarily small mass can be formed. In contrast to previously published work, we find that--for a specific value of the adiabatic index (Gamma = 2)--the observed Type II critical solution has approximately the same scaling exponent as that calculated for an ultrarelativistic fluid of the same index. Further, we find that the critical solution computed using the ideal-gas equations of state asymptotes to the ultrarelativistic critical solution.

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.

We present a brief synopsis of related work (gr-qc/0007039), describing a study of black hole threshold phenomena for a self-gravitating, massive complex scalar field in spherical symmetry. We construct Type I critical solutions dynamically by tuning a one-parameter family of initial data consisting of a boson star and a massless real scalar field, and numerically evolving this data. The resulting critical solutions appear to correspond to boson stars on the unstable branch, as we show via comparisons between our simulations and perturbation theory. For low-mass critical solutions, we find small ``halos'' of matter in the tails of the solutions, and these distort the profiles which otherwise agree with unstable boson stars. These halos seem to be artifacts of the collisions between the original boson stars and the massless fields, and do not appear to belong to the true critical solutions. From this study, it appears that unstable boson stars are unstable to dispersal (``explosion'') in addition to black hole formation. Given the similarities in macroscopic stability between boson stars and neutron stars, we suggest that similar phenomena could occur in models of neutron stars.

We present a study of black hole threshold phenomena for a self-gravitating, massive complex scalar field in spherical symmetry. We construct Type I critical solutions dynamically by tuning a one-parameter family of initial data composed of a boson star and a massless real scalar field. The real field is used to perturb the boson star via a gravitational interaction which results in a {\em significant} transfer of energy. The resulting critical solutions, which show great similarity with unstable boson stars, persist for a finite time before dispersing or forming a black hole. We extend the stability analysis of Gleiser and Watkins [Nucl. Phys. B319, 733 (1989)], providing a method for calculating the radial dependence of boson star modes of nonzero frequency. We find good agreement between our critical solutions and boson star modes. For critical solutions less than 90% of the maximum boson star mass $M_{\rm max} \simeq 0.633 M_{Pl}^2/m$, a small halo of matter appears in the tail of the solution. This halo appears to be an artifact of the collision between the original boson star and the real field, and does not belong to the true critical solution. It seems that unstable boson stars are unstable to dispersal in addition to black hole formation. Given the similarity in macroscopic stability between boson and neutron stars, we suggest that neutron stars at or beyond the point of instability may also be unstable to explosion.