Investigations of the dynamic modes of the Poincare gauge theory of gravity found only two good propagating torsion modes; they are effectively a scalar and a pseudoscalar. Cosmology affords a natural situation where one might see observational effects of these modes. Here we consider only the ``scalar torsion'' mode. This mode has certain distinctive and interesting qualities. In particular this type of torsion does not interact directly with any known matter and it allows a critical non-zero value for the affine scalar curvature. Via numerical evolution of the coupled nonlinear equations we show that this mode can contribute an oscillating aspect to the expansion rate of the Universe. From the examination of specific cases of the parameters and initial conditions we show that for suitable ranges of the parameters the dynamic ``scalar torsion'' model can display features similar to those of the presently observed accelerating universe.

For the Poincare gauge theory of gravity we consider the dynamical scalar torsion mode in a cosmological context. We explore in particular the possibility of using dynamical torsion to explain the current state of the accelerating Universe. With certain suitable sets of chosen parameters, this model can give a (qualitatively) proper description of the current universe without a cosmological constant, and the universe described is oscillating with a period of the Hubble time.

We present a formalism for constructing quasi-equilibrium binary black hole initial data suitable for numerical evolution. We construct quasi-equilibrium models by imposing an approximate helical Killing symmetry appropriate for quasi-circular orbits. We use the sum of two Kerr-Schild metrics as our background metric, thereby improving on conformally flat backgrounds that do not accommodate rotating black holes and providing a horizon-penetrating lapse, convenient for implementing black hole excision. We set inner boundary conditions at an excision radius well inside the apparent horizon and construct these boundary conditions to incorporate the quasi-equilibrium condition and recover the solution for isolated black holes in the limit of large separation. We use our formalism both to generate initial data for binary black hole evolutions and to construct a crude quasi-equilibrium, inspiral sequence for binary black holes of fixed irreducible mass.

We present a numerical code designed to study astrophysical phenomena involving dynamical spacetimes containing black holes in the presence of relativistic hydrodynamic matter. We present evolutions of the collapse of a fluid star from the onset of collapse to the settling of the resulting black hole to a final stationary state. In order to evolve stably after the black hole forms, we excise a region inside the hole before a singularity is encountered. This excision region is introduced after the appearance of an apparent horizon, but while a significant amount of matter remains outside the hole. We test our code by evolving accurately a vacuum Schwarzschild black hole, a relativistic Bondi accretion flow onto a black hole, Oppenheimer-Snyder dust collapse, and the collapse of nonrotating and rotating stars. These systems are tracked reliably for hundreds of M following excision, where M is the mass of the black hole. We perform these tests both in axisymmetry and in full 3+1 dimensions. We then apply our code to study the effect of the stellar spin parameter J/M^2 on the final outcome of gravitational collapse of rapidly rotating n = 1 polytropes. We find that a black hole forms only if J/M^2<1, in agreement with previous simulations. When J/M^2>1, the collapsing star forms a torus which fragments into nonaxisymmetric clumps, capable of generating appreciable ``splash'' gravitational radiation.

We experiment with modifications of the BSSN form of the Einstein field equations (a reformulation of the ADM equations) and demonstrate how these modifications affect the stability of numerical black hole evolution calculations. We use excision to evolve both non-rotating and rotating Kerr-Schild black holes in octant and equatorial symmetry, and without any symmetry assumptions, and obtain accurate and stable simulations for specific angular momenta J/M of up to about 0.9M.

We evolve a scalar field in a fixed Kerr-Schild background geometry to test simple $(3+1)$-dimensional algorithms for singularity excision. We compare both centered and upwind schemes for handling the shift (advection) terms, as well as different approaches for implementing the excision boundary conditions, for both static and boosted black holes. By first determining the scalar field evolution in a static frame with a $(1+1)$-dimensional code, we obtain the solution to very high precision. This solution then provides a useful testbed for simulations in full $(3+1)$ dimensions. We show that some algorithms which are stable for non-boosted black holes become unstable when the boost velocity becomes high.

A quasi-equilibrium (QE) computational scheme was recently developed in general relativity to calculate the complete gravitational wavetrain emitted during the inspiral phase of compact binaries. The QE method exploits the fact that the the gravitational radiation inspiral timescale is much longer than the orbital period everywhere outside the ISCO. Here we demonstrate the validity and advantages of the QE scheme by solving a model problem in relativistic scalar gravitation theory. By adopting scalar gravitation, we are able to numerically track without approximation the damping of a simple, quasi-periodic radiating system (an oscillating spherical matter shell) to final equilibrium, and then use the exact numerical results to calibrate the QE approximation method. In particular, we calculate the emitted gravitational wavetrain three different ways: by integrating the exact coupled dynamical field and matter equations, by using the scalar-wave monopole approximation formula (corresponding to the quadrupole formula in general relativity), and by adopting the QE scheme. We find that the monopole formula works well for weak field cases, but fails when the fields become even moderately strong. By contrast, the QE scheme remains quite reliable for moderately strong fields, and begins to breakdown only for ultra-strong fields. The QE scheme thus provides a promising technique to construct the complete wavetrain from binary inspiral outside the ISCO, where the gravitational fields are strong, but where the computational resources required to follow the system for more than a few orbits by direct numerical integration of the exact equations are prohibitive.