We present results from general relativistic calculations of nuclear ignition in white dwarf stars triggered by near encounters with rotating intermediate mass black holes with different spin and alignment parameters. These encounters create thermonuclear environments characteristic of Type Ia supernovae capable of producing both calcium and iron group elements in arbitrary ratios, depending primarily on the proximity of the interaction which acts as a strong moderator of nucleosynthesis. We explore the effects of black hole spin and spin-orbital alignment on burn product synthesis to determine whether they might also be capable of moderating reactive flows. When normalized to equivalent impact penetration, accounting for frame dragging corrections, the influence of spin is weak, no more than 25% as measured by nuclear energy release and mass of burn products, even for near maximally rotating black holes. Stars on prograde trajectories approach closer to the black hole and produce significantly more unbound debris and iron group elements than is possible by encounters with nonrotating black holes or by retrograde orbits, at more than 50% mass conversion efficiency. The debris contains several radioisotopes, most notably Ni56, made in amounts that produce sub-luminous (but still observable) light curves compared to branch-normal SNe Ia.

Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.

We present results from general relativistic calculations of the tidal disruption of white dwarf stars from near encounters with intermediate mass black holes. We follow the evolution of 0.2 and $0.6 M_\odot$ stars on parabolic trajectories that approach $10^3$ - $10^4 M_\odot$ black holes as close as a few Schwarzschild radii at periapsis, paying particular attention to the effect tidal disruption has on thermonuclear reactions and the synthesis of intermediate to heavy ion elements. These encounters create diverse thermonuclear environments characteristic of Type I supernovae and capable of producing both intermediate and heavy mass elements in arbitrary ratios, depending on the strength (or proximity) of the interaction. Nuclear ignition is triggered in all of our calculations, even at weak tidal strengths $\beta \sim 2.6$ and large periapsis radius $R_P \sim 28$ Schwarzschild radii. A strong inverse correlation exists between the mass ratio of calcium to iron group elements and tidal strength, with $\beta \lesssim 5$ producing predominately calcium-rich debris. At these moderate to weak interactions, nucleosynthesis is not especially efficient, limiting the total mass and outflows of calcium group elements to $< 15$\% of available nuclear fuel. Iron group elements however continue to be produced in greater quantity and ratio with increasing tidal strength, peaking at $\sim 60$\% mass conversion efficiency in our closest encounter cases. These events generate short bursts of gravitational waves with characteristic frequencies 0.1-0.7 Hz and strain amplitudes $0.5\times10^{-22}$ - $3.5\times10^{-22}$ at 10 Mpc source distance.

We present results from two-dimensional, general relativistic, viscous, radiation hydrodynamic numerical simulations of Shakura-Sunyaev thin disks accreting onto stellar mass Schwarzschild black holes. We consider cases on both the gas- and radiation-pressure-dominated branches of the thermal equilibrium curve, with mass accretion rates spanning the range from $\dot{M} = 0.01 L_\mathrm{Edd}/c^2$ to $10 L_\mathrm{Edd}/c^2$. The simulations directly test the stability of this standard disk model on the different branches. We find clear evidence of thermal instability for all radiation-pressure-dominated disks, resulting universally in the vertical collapse of the disks, which in some cases then settle onto the stable, gas-pressure-dominated branch. Although these results are consistent with decades-old theoretical predictions, they appear to be in conflict with available observational data from black hole X-ray binaries. We also find evidence for a radiation-pressure-driven instability that breaks the unstable disks up into alternating rings of high and low surface density on a timescale comparable to the thermal collapse. Since radiation is included self-consistently in the simulations, we are able to calculate lightcurves and power density spectra (PDS). For the most part, we measure radiative efficiencies (ratio of luminosity to mass accretion rate) close to 6%, as expected for a non-rotating black hole. The PDS appear as broken power laws, with a break typically around 100 Hz. There is no evidence of significant excess power at any frequencies, i.e. no quasi-periodic oscillations are observed.

We present results of three-dimensional, multi-physics simulations of an AGN jet colliding with an intergalactic cloud. The purpose of these simulations is to assess the degree of "positive feedback," i.e. jet-induced star formation, that results. We have specifically tailored our simulation parameters to facilitate comparison with recent observations of Minkowski's Object (M.O.), a stellar nursery located at the termination point of a radio jet coming from galaxy NGC 541. As shown in our simulations, such a collision triggers shocks which propagate around and through the cloud. These shocks condense the gas and under the right circumstances may trigger cooling instabilities, creating runaway increases in density, to the point that individual clumps can become Jeans unstable. Our simulations provide information about the expected star formation rate, total mass converted to \ion{H}{1}, H$_2$, and stars, and the relative velocity of the stars and gas. Our results confirm the possibility of jet-induced star formation, and agree well with the observations of M.O.

We have extended Cosmos++, a multi-dimensional unstructured adaptive mesh code for solving the covariant Newtonian and general relativistic radiation magnetohydrodynamic (MHD) equations, to accommodate both discrete finite volume and arbitrarily high order finite element structures. The new finite element implementation, called CosmosDG, is based on a discontinuous Galerkin (DG) formulation, using both entropy-based artificial viscosity and slope limiting procedures for regularization of shocks. High order multi-stage forward Euler and strong stability preserving Runge-Kutta time integration options complement high order spatial discretization. We have also added flexibility in the code infrastructure allowing for both adaptive mesh and adaptive basis order refinement to be performed separately or simultaneously in a local (cell-by-cell) manner. We discuss in this report the DG formulation and present tests demonstrating the robustness, accuracy, and convergence of our numerical methods applied to special and general relativistic MHD, though we note an equivalent capability currently also exists in CosmosDG for Newtonian systems.

We present the Grackle chemistry and cooling library for astrophysical simulations and models. Grackle provides a treatment of non-equilibrium primordial chemistry and cooling for H, D, and He species, including H2 formation on dust grains; tabulated primordial and metal cooling; multiple UV background models; and support for radiation transfer and arbitrary heat sources. The library has an easily implementable interface for simulation codes written in C, C++, and Fortran as well as a Python interface with added convenience functions for semi-analytical models. As an open-source project, Grackle provides a community resource for accessing and disseminating astrochemical data and numerical methods. We present the full details of the core functionality, the simulation and Python interfaces, testing infrastructure, performance, and range of applicability. Grackle is a fully open-source project and new contributions are welcome.

In this paper we report on recent upgrades to our general relativistic radiation magnetohydrodynamics code, Cosmos++, including the development of a new primitive inversion scheme and a hybrid implicit-explicit solver with a more general closure relation for the radiation equations. The new hybrid solver helps stabilize the treatment of the radiation source terms, while the new closure allows for a much broader range of optical depths to be considered. These changes allow us to expand by orders of magnitude the range of temperatures, opacities, and mass accretion rates, and move a step closer toward our goal of performing global simulations of radiation-pressure-dominated black hole accretion disks. In this work we test and validate the new method against an array of problems. We also demonstrate its ability to handle super-Eddington, quasi-spherical accretion. Even with just a single proof-of-principle simulation, we already see tantalizing hints of the interesting phenomenology associated with the coupling of radiation and gas in super-Eddington accretion flows.

Using three-dimensional, moving-mesh simulations, we investigate the future evolution of the recently discovered gas cloud G2 traveling through the galactic center. We consider the case of a spherical cloud initially in pressure equilibrium with the background. Our suite of simulations explores the following parameters: the equation of state, radial profiles of the background gas, and start times for the evolution. Our primary focus is on how the fate of this cloud will affect the future activity of Sgr A*. From our simulations we expect an average feeding rate in the range of $5-19 \times 10^{-8}$ solar masses per year beginning in 2013 and lasting for at least 7 years (our simulations stop in year 2020). The accretion varies by less than a factor of three on timescales <1 month, and shows no more than a factor of 10 difference between the maximum and minimum. These rates are comparable to the current estimated accretion rate in the immediate vicinity of Sgr A*, although they represent only a small (<5%) increase over the current expected feeding rate at the effective inner boundary of our simulations. Therefore, the break up of cloud G2 may have only a minimal effect on the brightness and variability of Sgr A* over the next decade. This is because current models of the galactic center predict that most of the gas will be caught up in outflows. However, if the accreted G2 material can remain cold, it may not mix well with the hot, diffuse background gas, and instead accrete efficiently onto Sgr A*. Further observations of G2 will give us an unprecedented opportunity to test this idea. The break up of the cloud itself may also be observable. By tracking the amount of cloud energy that is dissipated during our simulations, we are able to get a rough estimate of the luminosity associated with its tidal disruption; we find values of a few $10^{36}$ erg/s.

Using three-dimensional, moving-mesh simulations, we investigate the future evolution of the recently discovered gas cloud G2 traveling through the galactic center. From our simulations we expect an average feeding rate onto Sgr A* in the range of $(5-19) \times 10^{-8} M_\odot\mathrm{~yr}^{-1}$ beginning in 2014. This accretion varies by less than a factor of three on timescales of about 1 month, and shows no more than a factor of 10 difference between the maximum and minimum observed rates within any given model. These rates are comparable to the current estimated accretion rate in the immediate vicinity of Sgr A*, although they represent only a small (<10%) increase over the current expected feeding rate at the effective inner boundary of our simulations $(r_\mathrm{acc} = 750 R_S \sim 10^{15} \mathrm{cm})$. We also explore multiple possible equations of state to describe the gas. In examining the Br-$\gamma$ light curves produced from our simulations, we find that all of our isothermal models predict significant (factor of 10) enhancements in the luminosity of G2 as it approaches pericenter, in conflict with observations. Models that instead allow the cloud to heat as it is compressed do better at matching observations.

Modeling the radiation generated by accreting matter is an important step towards realistic simulations of black hole accretion disks, especially at high accretion rates. To this end, we have recently added radiation transport to the existing general relativistic magnetohydrodynamic code, Cosmos++. However, before attempting to model radiative accretion disks, we have tested the new code using a series of shock tube and Bondi (spherical inflow) problems. The four radiative shock tube tests, first presented by Farris et al. (2008), have known analytic solutions, allowing us to calculate errors and convergence rates for our code. The Bondi problem only has an analytic solution when radiative processes are ignored, but is pertinent because it is closer to the physics we ultimately want to study. In our simulations, we include Thomson scattering and thermal bremsstrahlung in the opacity, focusing exclusively on the super-Eddington regime. Unlike accretion onto bodies with solid surfaces, super-Eddington accretion onto black holes does not produce super-Eddington luminosity. In our examples, despite accreting at up to 300 times the Eddington rate, our measured luminosity is always several orders of magnitude below Eddington.

This paper presents a numerical study over a wide parameter space of the likelihood of the dynamical bar-mode instability in differentially rotating magnetized neutron stars. The innovative aspect of this study is the incorporation of magnetic fields in such a context, which have thus far been neglected in the purely hydrodynamical simulations available in the literature. The investigation uses the Cosmos++ code which allows us to perform three dimensional simulations on a cylindrical grid at high resolution. A sample of Newtonian magneto-hydrodynamical simulations starting from a set of models previously analyzed by other authors without magnetic fields has been performed, providing estimates of the effects of magnetic fields on the dynamical bar-mode deformation of rotating neutron stars. Overall, our results suggest that the effect of magnetic fields are not likely to be very significant in realistic configurations. Only in the most extreme cases are the magnetic fields able to suppress growth of the bar mode.

In recent work we presented the first results of global general relativistic magnetohydrodynamic (GRMHD) simulations of tilted (or misaligned) accretion disks around rotating black holes. The simulated tilted disks showed dramatic differences from comparable untilted disks, such as asymmetrical accretion onto the hole through opposing "plunging streams" and global precession of the disk powered by a torque provided by the black hole. However, those simulations used a traditional spherical-polar grid that was purposefully underresolved along the pole, which prevented us from assessing the behavior of any jets that may have been associated with the tilted disks. To address this shortcoming we have added a block-structured "cubed-sphere" grid option to the Cosmos++ GRMHD code, which will allow us to simultaneously resolve the disk and polar regions. Here we present our implementation of this grid and the results of a small suite of validation tests intended to demonstrate that the new grid performs as expected. The most important test in this work is a comparison of identical tilted disks, one evolved using our spherical-polar grid and the other with the cubed-sphere grid. We also demonstrate an interesting dependence of the early-time evolution of our disks on their orientation with respect to the grid alignment. This dependence arises from the differing treatment of current sheets within the disks, especially whether they are aligned with symmetry planes of the grid or not.

This paper presents a continuation of our efforts to numerically study accretion disks that are misaligned (tilted) with respect to the rotation axis of a Kerr black hole. Here we present results of a global numerical simulation which fully incorporates the effects of the black hole spacetime as well as magnetorotational turbulence that is the primary source of angular momentum transport in the flow. This simulation shows dramatic differences from comparable simulations of untilted disks. Accretion onto the hole occurs predominantly through two opposing plunging streams that start from high latitudes with respect to both the black-hole and disk midplanes. This is due to the aspherical nature of the gravitational spacetime around the rotating black hole. These plunging streams start from a larger radius than would be expected for an untilted disk. In this regard the tilted black hole effectively acts like an untilted black hole of lesser spin. Throughout the duration of the simulation, the main body of the disk remains tilted with respect to the symmetry plane of the black hole; thus there is no indication of a Bardeen-Petterson effect in the disk at large. The torque of the black hole instead principally causes a global precession of the main disk body. In this simulation the precession has a frequency of $3 (M_\odot/M)$ Hz, a value consistent with many observed low-frequency quasi-periodic oscillations. However, this value is strongly dependent on the size of the disk, so this frequency may be expected to vary over a large range.

We posit that accreting compact objects, including stellar mass black holes and neutron stars as well as supermassive black holes, may undergo extended periods of accretion during which the angular momentum of the disk at large scales is misaligned with that of the compact object. In such a scenario, Lense-Thirring precession caused by the rotating compact object can dramatically affect the disk. In this presentation we describe results from a three-dimensional relativistic magnetohydrodynamic simulation of an MRI turbulent disk accreting onto a tilted rapidly rotating black hole. For this case, the disk does not achieve the commonly described Bardeen-Petterson configuration; rather, it remains nearly planar, undergoing a slow global precession. Accretion from the disk onto the hole occurs predominantly through two opposing plunging streams that start from high latitudes with respect to both the black-hole and disk midplanes. This is a consequence of the non-sphericity of the gravitational spacetime of the black hole.

We use the relativistic hydrodynamics code Cosmos++ to model the evolution of the radio nebula triggered by the Dec. 27, 2004 giant flare event of soft gamma repeater 1806-20. We primarily focus on the rebrightening and centroid motion occurring subsequent to day 20 following the flare event. We model this period as a mildly relativistic (gamma ~ 1.07 - 1.67) jetted outflow expanding into the interstellar medium (ISM). We demonstrate that a jet with total energy ~ 10^46 ergs confined to a half opening angle ~ 20 degrees fits the key observables of this event, e.g. the flux lightcurve, emission map centroid position, and aspect ratio. In particular, we find excellent agreement with observations if the rebrightening is due to the jet, moving at 0.5c and inclined ~ 0 - 40 degrees toward the observer, colliding with a density discontinuity in the ISM at a radius of several 10^16 cm. We also find that a jet with a higher velocity, >~ 0.7c, and larger inclination, >~ 70 degrees, moving into a uniform ISM can fit the observations in general, but tends to miss the details of rebrightening. The latter, uniform ISM model predicts an ISM density more than 100 times lower than that of the former model, and thus suggests an independent test which might discriminate between the two. One of the strongest constraints of both models is that the data seems to require a non-uniform jet in order to be well fit.

A new code and methodology are introduced for solving the general relativistic magnetohydrodynamic (GRMHD) equations in fixed background spacetimes using time-explicit, finite-volume discretization. The code has options for solving the GRMHD equations using traditional artificial-viscosity (AV) or non-oscillatory central difference (NOCD) methods, or a new extended AV (eAV) scheme using artificial-viscosity together with a dual energy-flux-conserving formulation. The dual energy approach allows for accurate modeling of highly relativistic flows at boost factors well beyond what has been achieved to date by standard artificial viscosity methods. It provides the benefit of Godunov methods in capturing high Lorentz boosted flows but without complicated Riemann solvers, and the advantages of traditional artificial viscosity methods in their speed and flexibility. Additionally, the GRMHD equations are solved on an unstructured grid that supports local adaptive mesh refinement using a fully threaded oct-tree (in three dimensions) network to traverse the grid hierarchy across levels and immediate neighbors. A number of tests are presented to demonstrate robustness of the numerical algorithms and adaptive mesh framework over a wide spectrum of problems, boosts, and astrophysical applications, including relativistic shock tubes, shock collisions, magnetosonic shocks, Alfven wave propagation, blast waves, magnetized Bondi flow, and the magneto-rotational instability in Kerr black hole spacetimes.

We present results from fully general relativistic three-dimensional numerical studies of thick-disk accretion onto a rapidly-rotating (Kerr) black hole with a spin axis that is tilted (not aligned) with the angular momentum vector of the disk. We initialize the problem with the solution for an aligned, constant angular momentum, accreting thick disk, which is then allowed to respond to the Lense-Thirring precession of the tilted black hole. The precession causes the disk to warp, beginning at the inner edge and moving out on roughly the Lense-Thirring precession timescale. The propagation of the warp stops at a radius in the disk at which other dynamical timescales, primarily the azimuthal sound-crossing time, become shorter than the precession time. At this point, the warp effectively freezes into the disk and the evolution becomes quasi-static, except in cases where the sound-crossing time in the bulk of the disk is shorter than the local precession timescale. We see evidence that such disks undergo near solid-body precession after the initial warping has frozen in. Simultaneous to the warping of the disk, there is also a tendency for the midplane to align with the symmetry plane of the black hole due to the preferential accretion of the most tilted disk gas. This alignment is not as pronounced, however, as it would be if more efficient angular momentum transport (e.g. from viscosity or magneto-rotational instability) were considered.

We present results from two-dimensional numerical simulations of the interactions between magnetized shocks and radiative clouds. Our primary goal is to characterize the dynamical evolution of the shocked clouds. We perform runs in both the strong and weak magnetic field limits and consider three different field orientations. For the geometries considered, we generally find that magnetic fields external to, but concentrated near, the surface of the cloud suppress the growth of destructive hydrodynamic instabilities. External fields also increase the compression of the cloud by effectively acting as a confinement mechanism driven by the interstellar flow and local field stretching. This can have a dramatic effect on both the efficiency of radiative cooling, which tends to increase with increasing magnetic field strength, and on the size and distribution of condensed cooled fragments. In contrast, fields acting predominately internally to the cloud tend to resist compression, thereby inhibiting cooling. We observe that, even at modest strengths, internal fields can completely suppress low-temperature cooling.

We discuss obbservations and numerical simulations which show that radiative shocks in jet-cloud collisions can trigger the collapse of intergalactic clouds and subsequent star formation in low luminosity, 'FR-I' type, radio galaxies.

Jets from radio galaxies can have dramatic effects on the medium through which they propagate. We review observational evidence for jet-induced star formation in low ('FR-I') and high ('FR-II') luminosity radio galaxies, at low and high redshifts respectively. We then discuss numerical simulations which are aimed to explain a jet-induced starburst ('Minkowski's Object') in the nearby FR-I type radio galaxy NGC 541. We conclude that jets can induce star formation in moderately dense (10 cm-3), warm (10,000 K) gas; that this may be more common in the dense environments of forming, active galaxies; and that this may provide a mechanism for 'positive' feedback from AGN in the galaxy formation process.

Accumulating observational evidence for a number of radio galaxies suggests an association between their jets and regions of active star formation. The standard picture is that shocks generated by the jet propagate through an inhomogeneous medium and trigger the collapse of overdense clouds, which then become active star-forming regions. In this contribution, we report on recent hydrodynamic simulations of radiative shock-cloud interactions using two different cooling models: an equilibrium cooling-curve model assuming solar metallicities and a non-equilibrium chemistry model appropriate for primordial gas clouds. We consider a range of initial cloud densities and shock speeds in order to quantify the role of cooling in the evolution. Our results indicate that for moderate cloud densities (>1 cm^{-3}) and shock Mach numbers (<20), cooling processes can be highly efficient and result in more than 50% of the initial cloud mass cooling to below 100 K. We also use our results to estimate the final H_2 mass fraction for the simulations that use the non-equilibrium chemistry package. This is an important measurement, since H_2 is the dominant coolant for a primordial gas cloud. We find peak H_2 mass fractions of >0.01 and total H_2 mass fractions of >10^{-5} for the cloud gas. Finally, we compare our results with the observations of jet-induced star formation in ``Minkowski's Object.'' We conclude that its morphology, star formation rate (~ 0.3M_solar/yr) and stellar mass (~ 1.2 x 10^7 M_solar) can be explained by the interaction of a 90,000 km/s jet with an ensemble of moderately dense (~ 10 cm^{-3}), warm (10^4 K) intergalactic clouds in the vicinity of its associated radio galaxy at the center of the galaxy cluster.

(Abridged) Many dwarf galaxies exhibit sub-Solar metallicities, with some star-to-star variation, despite often containing multiple generations of stars. The total metal content in these systems is much less than expected from the heavy element production of massive stars in each episode of star formation. Such a deficiency implies that a substantial fraction of the enriched material has been lost from these small galaxies. Mass ejection from dwarf galaxies may have important consequences for the evolution of the intergalactic medium and for the evolution of massive galaxies, which themselves may have formed via the merger of smaller systems. We report here the results of three-dimensional simulations of the evolution of supernova-enriched gas within dwarf spheroidal galaxies (dSph's), with the aim of determining the retention efficiency of supernova ejecta. We consider two galaxy models, selected to represent opposite ends of the dSph sequence. For each model galaxy we investigate a number of scenarios, ranging from a single supernova in smooth gas distributions to more complex multiple supernovae in highly disturbed gas distributions. The results of these investigations suggest that, for low star-formation efficiencies, it is difficult to completely expel the enriched material from the galaxy. Most of the enriched gas is, however, lost from the core of the galaxy following multiple supernovae, especially if the interstellar medium is already highly disturbed by processes such as photo-ionization and stellar winds. If subsequent star formation occurs predominantly within the core where most of the residual gas is concentrated, then these results could explain the poor self-enrichment efficiency observed in dwarf galaxies.

We have developed a new massively-parallel radiation-hydrodynamics code (Cosmos) for Newtonian and relativistic astrophysical problems that also includes radiative cooling, self-gravity, and non-equilibrium, multi-species chemistry. Several numerical methods are implemented for the hydrodynamics, including options for both internal and total energy conserving schemes. Radiation is treated using flux-limited diffusion. The chemistry incorporates 27 reactions, including both collisional and radiative processes for atomic hydrogen and helium gases, and molecular hydrogen chains. In this paper we discuss the equations and present results from test problems carried out to verify the robustness and accuracy of our code in the Newtonian regime. An earlier paper presented tests of the relativistic capabilities of Cosmos.

Recent results from linear perturbation theory suggest that first-order cosmological quark-hadron phase transitions occurring as deflagrations may be ``borderline'' unstable, and those occurring as detonations may give rise to growing modes behind the interface boundary. However, since nonlinear effects can play important roles in the development of perturbations, unstable behavior cannot be asserted entirely by linear analysis, and the uncertainty of these recent studies is compounded further by nonlinearities in the hydrodynamics and self-interaction fields. In this paper we investigate the growth of perturbations and the stability of quark-hadron phase transitions in the early Universe by solving numerically the fully nonlinear relativistic hydrodynamics equations coupled to a scalar field with a quartic self-interaction potential regulating the transitions. We consider single, perturbed, phase transitions propagating either by detonation or deflagration, as well as multiple phase and shock front interactions in 1+2 dimensional spacetimes.

High resolution, non-oscillatory, central difference (NOCD) numerical schemes are introduced as alternatives to more traditional artificial viscosity (AV) and Godunov methods for solving the fully general relativistic hydrodynamics equations. These new approaches provide the advantages of Godunov methods in capturing ultra-relativistic flows without the cost and complication of Riemann solvers, and the advantages of AV methods in their speed, ease of implementation, and general applicability without explicitly using artificial viscosity for shock capturing. Shock tube, wall shock, and dust accretion tests, all with adiabatic equations of state, are presented and compared against equivalent solutions from both AV and Godunov based codes. In the process we address the accuracy of time-explicit NOCD and AV methods over a wide range of Lorentz factors.

We present results from numerical evolutions of single black holes distorted by axisymmetric, but equatorially asymmetric, gravitational (Brill) waves. Net radiated energies, apparent horizon embeddings, and recoil velocities are shown for a range of Brill wave parameters, including both even and odd parity distortions of Schwarzschild black holes. We find that a wave packet initially concentrated on the black hole throat, a likely model also for highly asymmetric stellar collapse and late stage binary mergers, can generate a maximum recoil velocity of about 150 (23) km/sec for even (odd) parity perturbations, significantly less than that required to eject black holes from galactic cores.

We investigate the properties of the Lyman alpha forest as predicted by numerical simulations for a range of currently viable cosmological models. This is done in order to understand the dependencies of the forest on cosmological parameters. Focusing on the redshift range from two to four, we show that: (1) most of the evolution in the distributions of optical depth, flux and column density can be understood by simple scaling relations, (2) the shape of optical depth distribution is a sensitive probe of the amplitude of density fluctuations on scales of a few hundred kpc, (3) the mean of the b distribution (a measure of the width of the absorption lines) is also very sensitive to fluctuations on these scales, and decreases as they increase. We perform a preliminary comparison to observations, where available. A number of other properties are also examined, including the evolution in the number of lines, the two-point flux distribution and the HeII opacity.

The comparison of quasar absorption spectra with numerically simulated spectra from hierarchical cosmological models of structure formation promises to be a valuable tool to discriminate among these models. We present simulation results for the column density, Doppler b parameter, and optical depth probability distributions for five popular cosmological models.

In this paper we critically examine predictions of the Lya forest within the standard cold dark matter (SCDM) model, paying particular attention to the low end of the column-density distribution. We show in particular that the width of these lines, typically measured by the b-parameter of a Voigt profile, is sensitive to spatial resolution in numerical simulations and has previously been overestimated. The new result, which predicts a distribution with a median b of around 20-22 km/s at z=3, is substantially below that observed. We examine a number of possible causes of this discrepancy and argue that it is unlikely to be rectified by an increase in the thermal broadening of the absorbing gas, but is instead telling us something about the distribution of matter on these scales. Although the median differs, the shape of the b-parameter distribution agrees quite well with that observed, and the high-end tail is naturally produced by the filamentary nature of gravitational collapse in these models. In particular, we demonstrate that lines of sight which obliquely intersect a filament or sheet tend to produce absorption lines with larger b parameters. We also examine the physical nature of the gas which is responsible for the forest, showing that for lines with neutral column densities below N_HI ~ 10^14 cm^-2 (for this model at z=3), the peculiar infall velocity is actually slower than the Hubble flow, while larger lines have, on average, turned around and are collapsing.

In this paper we review the importance of the Lyman alpha forest as a probe of structure formation in the universe. We first discuss the statistics used to describe the Lyman alpha forest and the numerical techniques used to produce simulated spectra of the forest from a given cosmological model. We then discuss the physical picture of the absorbing structures that emerges from these numerical simulations. Finally, we comment on how two of the statistics, the slope of the column density distribution and the b parameter distribution, may be used to constrain competing cosmologies.

We present results from the first fully nonlinear numerical calculations of the head--on collision of two unequal mass black holes. Selected waveforms of the most dominant l=2, 3 and 4 quasinormal modes are shown, as are the total radiated energies and recoil velocities for a range of mass ratios and initial separations. Our results validate the close and distant separation limit perturbation studies, and suggest that the head--on collision scenario is not likely to produce an astrophysically significant recoil effect.

We examine the structure of baryonic wakes in the cosmological fluid which would form behind GUT-scale cosmic strings at early times (redshifts z > 100) in a neutrino-dominated universe. We show, using simple analytical arguments as well as 1- and 2-dimensional hydrodynamical simulations, that these wakes will NOT be able to form interesting cosmological objects before the neutrino component collapses. The width of the baryonic wakes (< 10 kpc comoving) is smaller than the scale of wiggles on the strings and are probably not enhanced by the wiggliness of the string network.

We discuss the origin and physical nature of the Lya forest absorption systems as found in hydrodynamical simulations of the Intergalactic Medium in a standard Cold Dark Matter cosmology (Omega=1, H0=50 km/s/Mpc, sigma_8=0.7). The structures of the systems that give rise to the Lya forest span a wide range in morphologies, depending on the density contrast. The absorption properties of the forest clouds correlate strongly with those of the underlying physical systems from which they arise. The highest column density systems (log N_HI > 15), correspond to highly overdense spheroidal structures, moderate column density systems (13 < log N_HI < 14), correspond to filaments, and the lowest density absorption systems originate from discrete fluctuations within underdense regions a few megaparsecs across, cosmic minivoids. Most of the intergalactic HeII opacity arises from these underdense regions. We are able to account for the distribution of optical depths in our synthesized spectra entirely by absorption due to discrete systems. We find that virtually all the baryons in the simulation fragment into structures that we can identify with discrete absorption lines, with at most 5% remaining in a smoothly distributed component (the Gunn-Peterson effect). We compare our results with the cloud ionization parameters inferred from Keck HIRES measurements of carbon and silicon in the Lya forest. Combining with constraints imposed by measurements of the mean intergalactic HI opacity, we find 0.03 < Omega_b < 0.08 and 0.3 < (Gamma_HI/ 10^-12 s^-1) < 1 at z=3-3.5.

We investigate the possibility of very early formation of primordial star clusters from high-\sigma perturbations in cold dark matter dominated structure formation scenarios. For this we have developed a powerful 2-level hierarchical cosmological code with a realistic and robust treatment of multi-species primordial gas chemistry, paying special attention to the formation and destruction of hydrogen molecules, non-equilibrium ionization, and cooling processes. We performed 3-D simulations at small scales and at high redshifts and find that, analogous to simulations of large scale structure, a complex system of filaments, sheets, and spherical knots at the intersections of filaments form. On the mass scales covered by our simulations (5x10^5 - 1x10^9\Ms) that collapse at redshifts z>25, we find that only at the spherical knots can enough H2 be formed (n_{H_2}/n_H > 5x10^-4) to cool the gas appreciably. Quantities such as the time dependence of the formation of H2 molecules, the final H2 fraction, and central densities from the simulations are compared to the theoretical predictions of Abel (1995) and Tegmark et al. (1997) and found to agree remarkably well. Comparing the 3-D results to an isobaric collapse model we further discuss the possible implications of the extensive merging of small structure that is inherent in hierarchical models. Typically only 5-8% percent of the total baryonic mass in the collapsing structures is found to cool significanlty. Assuming the Padoan (1995) model for star formation our results would predict the first stellar systems to be as small as ~30\Ms. Some implications for primordial globular cluster formation scenarios are also discussed.

We simulate the Lya forest in a standard CDM universe using a 2-level hierarchical grid code to evolve the dark and baryonic matter components self-consistently. We solve the time-dependent ionization equations for hydrogen and helium, adopting the Haardt & Madau (1996) estimate for the metagalactic UV radiation background. We compare our simulation results with the measured properties of the Lya forest by constructing synthetic spectra which we analyze in a manner similar to the method used by observers. The HI column density and Doppler parameter distributions we obtain agree closely with those measured by the Keck HIRES and earlier high spectral resolution observations over the column density range 10^{12} < NHI < 10^{16}. In particular, we find a power-law column density distribution persists to the lowest values reported from the HIRES measurements, in agreement with the incompleteness-corrected observations. We find evolution in the cloud number density and opacity comparable to the observed evolution. We also compare with measured values of the intergalactic HeII opacity. Our results require a HeII ionizing background lower than the Haardt & Madau estimate by a factor of 4, corresponding to a soft intrinsic QSO spectrum of alpha_Q = 1.8 - 2.

We developed a three-dimensional 2-level hierarchical cosmological code with a realistic and robust treatment of multi-species non-equilibrium ionization and radiative cooling processes, and use it to investigate primordial star forming regions that originate from high-\sigma perturbations in a standard CDM dominated cosmology. We find it is possible to produce gravitationally bound and cooled structures at very high redshift (z ~ 40) with baryonic masses as small as ~1000Ms. The molecular hydrogen formation in these small scale structures follows very well the analytical predictions of Abel (1995) and Tegmark et al. (1996). We also discuss the minimum mass that cosmological structures must have in order to be able to cool and collapse.

We have reviewed the chemistry and cooling behaviour of low-density (n<10^4 cm^-3) primordial gas and devised a cooling model wich involves 19 collisional and 9 radiative processes and is applicable for temperatures in the range (1 K < T < 10^8 K). We derived new fits of rate coefficients for the photo-attachment of neutral hydrogen, the formation of molecular hydrogen via H-, charge exchange between H2 and H+, electron detachment of H- by neutral hydrogen, dissociative recombination of H2 with slow electrons, photodissociation of H2+, and photodissociation of H2. Further it was found that the molecular hydrogen produced through the gas-phase processes, H2+ + H -> H2 + H+, and H- + H -> H2 + e-, is likely to be converted into its para configuration on a faster time scale than the formation time scale. We have tested the model extensively and shown it to agree well with former studies. We further studied the chemical kinetics in great detail and devised a minimal model which is substantially simpler than the full reaction network but predicts correct abundances. This minimal model shows convincingly that 12 collisional processes are sufficient to model the H, He, H+, H-, He+, He++, and H2 abundances in low density primordial gas for applications with no radiation fields.

We have developed a method of solving for multi-species chemical reaction flows in non--equilibrium and self--consistently with the hydrodynamic equations in an expanding FLRW universe. The method is based on a backward differencing scheme for the required stability when solving stiff sets of equations and is designed to be efficient for three-dimensional calculations without sacrificing accuracy. In all, 28 kinetic reactions are solved including both collisional and radiative processes for the following nine separate species: H, H+, He, He+, He++, H-, H2+, H2, and e-. The method identifies those reactions (involving H- and H2+) ocurring on the shortest time scales, decoupling them from the rest of the network and imposing equilibrium concentrations to good accuracy over typical cosmological dynamical times. Several tests of our code are presented, including radiative shock waves, cosmological sheets, conservation constraints, and fully three-dimensional simulations of CDM cosmological evolutions in which we compare our method to results obtained when the packaged routine LSODAR is substituted for our algorithms.

We perform a double line of sight (DLOS) analysis of the Lyman-alpha forest structures that form and evolve in cosmological N-body/hydrodynamic simulations. Pairs of simulated spectra, extracted from lines of sight separated by distances from D=12.5kpc up to 800kpc, and a ``control sample'' of unrelated lines of sight, are analyzed at redshifts 3, 2, and 1. Coincident line samples are defined for HI column density thresholds of Nco = 10^{12.5}, 10^{13}, and 10^{14} per square cm. We find that: 1) Under the assumption of a single structure size, a Bayesian analysis yields sizes that are larger for smaller Nco, and at fixed Nco the size decreases with decreasing redshift. However, these derived sizes are found to increase with increasing D indicating that the assumption of a single structure size is invalid. 2) The column densities of coincident pairs are highly correlated for small D, with increasing scatter as D is increased, consistent with structures that have a centrally peaked N(HI) that decreases gradually with radius. 3) The velocity difference distribution for coincident lines is very narrow for small D, and widens as D is increased to meet the expectation for chance coincidences in unrelated lines of sight. This behavior is indicative of organized motion within the structures. 4) For small D, the distribution of anticoincident line column densities, Nac, falls steeply as Nac increases from the cutoff value, but has a significant tail at large values which is inconsistent with a population of spherical absorbers with sharp edges, and consistent with a flattened geometry. The conclusions reached on the basis of the DLOS analysis are validated by an examination of the three-dimensional structures and velocity flows in the simulation data.

Soon after the first QSO was identified, Gunn & Peterson searched for the expected characteristic absorption trough on the blueward side of Lya in the spectrum of the QSO due to an Intergalactic Medium (IGM). They failed to find it, placing a constraint on the density of neutral hydrogen in the IGM that was less than 1 part in 10^5 of that residing in galaxies, and concluded that the IGM must be highly ionized. Soon afterwards absorption by Lya was detected in the IGM; not by a diffuse component, but by a clumpy component of intergalactic gas clouds, the Lya forest. The goal then became to search for `excess' absorption beyond that expected from the forest. No clear absorption by a diffuse HI component, however, was ever detected. The results of recent numerical hydrodynamics computations of the formation of the Lya forest appear to indicate that the division between a clumpy component and a diffuse one may be inappropriate. We find that in a CDM-dominated cosmology, no substantial diffuse medium should be expected. Instead, the medium condenses into a network of complex structures that reveal themselves as discrete absorption systems in the spectra of QSOs. The lowest column density lines arise from the fine structure in minivoids -- small regions with densities below the cosmic mean. These results suggest that the long sought for diffuse HI Gunn-Peterson effect may not exist.

We investigate the imprint of nonlinear matter condensations on the Cosmic Microwave Background (CMB) in $\Omega_{0}<1$ cold dark matter (CDM) model universes. We consider simulation domains ranging from $120h^{-1}$ Mpc to $360h^{-1}$ Mpc in size. We concentrate on the secondary temperature anisotropies induced by time varying gravitational potentials occurring after decoupling. Specifically, we investigate the importance of the Rees-Sciama effect due to: (1) intrinsic changes in the gravitational potential of forming, nonlinear structures, (2) proper motion of nonlinear structures, and (3) late time decay of gravitational potential perturbations in open universes. CMB temperature anisotropies are obtained by numerically evolving matter inhomogeneities and CMB photons from an early, linear epoch ($z=100$) to the present, nonlinear epoch $(z=0)$. We test the dependence and relative importance of these secondary temperature anisotropies as a function of the scale of the underlying matter (voids, superclusters) and as a function of $\Omega_{0}$. The results of the $\Omega_{0}<1$ models are compared to a similarly executed $\Omega_{0}=1.0$ simulation. We find that in low density models all three sources of anisotropy could be relevant and reach levels of $\Delta T/T \sim 10^{-6}$. In particular, we find that for $\Omega_{0}<1$ at large scales, secondary temperature anisotropies are dominated by the decaying potential.

We have performed a multi-species hydrodynamical simulation of the formation and evolution of Lyman alpha clouds in a flat CDM dominated universe with an external flux of ionizing radiation. We solve the fully coupled non-equilibrium rate equations for the following species: H, H^+, H^-, H_2, H_2^+, He, He^+, He^{++}, and e^-. The statistical properties of the distribution and evolution of both hydrogen and helium absorption lines are extracted and compared to observed data. We find excellent agreement for the following neutral hydrogen data: the distribution of column densities is fit well by a power law with exponent beta=1.55 with a possible deficiency of lines above column density 10^{15} cm^{-2}; the integrated distribution matches observed data over a broad range of column densities 10^{13} to 10^{17} cm^{-2}; a Gaussian statistical fit to the Doppler parameter distribution yields a median of 35.6 km s^{-1}; the evolution of the number of clouds with column densities larger than 10^{14} cm^{-2} follows a power law with exponent gamma=2.22. Analogous calculations are presented for HeII absorption lines and we find the ratio of Doppler parameters b_{HeII}/b_{HI} = 0.87. Our data also suggests that Ly$\alpha$ clouds may belong to two morphologically different groups: small column density clouds which tend to reside in sheets or filamental structures and are very elongated and/or flattened, and the large column density clouds which are typically found at the intersections of these caustic structures and are quasi-spherical.

We trace the evolution of cosmic microwave background photons propagating through a reionized model universe. The reionization of the intergalactic medium is achieved by UV photons emitted from the decaying `hot' dark matter neutrinos. The model universe is represented by a comoving cube through which we propagate photons and subject them to the Sachs-Wolfe effect and Thomson scattering off the ionized fraction of the IGM, beginning at z=1900 until the present z=0. The simulation follows the evolving matter inhomogeneities and microwave background photons into the present nonlinear regime and yields temperature maps of the microwave sky as a comoving observer would detect. Our temperature maps display anisotropies on scales ranging from 2 to 8 degrees with rms amplitudes ranging from DT/T_{rms} =2.8 to 3.4 times 10^{-5} without reionization. If we allow the neutrinos to decay and subsequently reionize the IGM, the level of the temperature anisotropies is reduced to DT/T=0.8 to 1.1 times 10^{-5}. This includes the Doppler shifts suffered by the CMB photons as they scatter off the moving ionized fraction of the IGM. If, however, reionization is caused by a late and sudden energy input into the IGM, we find that some reduction in $\Delta T/T$ occurs, but not enough to make this scenario consistent with observations. We find that the decaying neutrino hypothesis is more efficient at reducing primordial anisotropies and generates smaller Doppler anisotropies than a model reionized via a late and sudden energy input into the IGM.