We study the origin of Na I absorbing gas in ultraluminous infrared galaxies motivated by the recent observations by Martin of extremely superthermal linewidths in this cool gas. We model the effects of repeated supernova explosions driving supershells in the central regions of molecular disks with M_d=10^10 M_\sun, using cylindrically symmetric gas dynamical simulations run with ZEUS-3D. The shocked swept-up shells quickly cool and fragment by Rayleigh-Taylor instability as they accelerate out of the dense, stratified disks. The numerical resolution of the cooling and compression at the shock fronts determines the peak shell density, and so the speed of Rayleigh-Taylor fragmentation. We identify cooled shells and shell fragments as Na I absorbing gas and study its kinematics. We find that simulations with a numerical resolution of \le 0.2 pc produce multiple Rayleigh-Taylor fragmented shells in a given line of sight. We suggest that the observed wide Na I absorption lines, <v> = 320 \pm 120 km s^-1 are produced by these multiple fragmented shells traveling at different velocities. We also suggest that some shell fragments can be accelerated above the observed average terminal velocity of 750 km s^-1 by the same energy-driven wind with an instantaneous starburst of \sim 10^9 M_\sun. The bulk of mass is traveling with the observed average shell velocity 330 \pm 100 km s^-1. Our results show that an energy-driven bubble causing Rayleigh-Taylor instabilities can explain the kinematics of cool gas seen in the Na I observations without invoking additional physics relying primarily on momentum conservation, such as entrainment of gas by Kelvin-Helmholtz instabilities, ram pressure driving of cold clouds by a hot wind, or radiation pressure acting on dust. (abridged)

Gravitational wave emission from the gravitational collapse of massive stars has been studied for more than three decades. Current state of the art numerical investigations of collapse include those that use progenitors with realistic angular momentum profiles, properly treat microphysics issues, account for general relativity, and examine non--axisymmetric effects in three dimensions. Such simulations predict that gravitational waves from various phenomena associated with gravitational collapse could be detectable with advanced ground--based and future space--based interferometric observatories.

One interesting class of gravitational radiation sources includes rapidly rotating astrophysical objects that encounter dynamical instabilities. We have carried out a set of simulations of rotationally induced instabilities in differentially rotating polytropes. An $n$=1.5 polytrope with the Maclaurin rotation law will encounter the $m$=2 bar instability at $T/|W| \gtrsim 0.27$. Our results indicate that the remnant of this instability is a persistent bar-like structure that emits a long-lived gravitational radiation signal. Furthermore, dynamical instability is shown to occur in $n$=3.33 polytropes with the $j$-constant rotation law at $T/|W| \gtrsim 0.14$. In this case, the dominant mode of instability is $m$=1. Such instability may allow a centrifugally-hung core to begin collapsing to neutron star densities on a dynamical timescale. If it occurs in a supermassive star, it may produce gravitational radiation detectable by LISA.

Dynamical instability is shown to occur in differentially rotating polytropes with N = 3.33 and $T/|W| \gtrsim 0.14$. This instability has a strong m=1 mode, although the m=2, 3, and 4 modes also appear. Such instability may allow a centrifugally-hung core to begin collapsing to neutron star densities on a dynamical timescale. The gravitational radiation emitted by such unstable cores may be detectable with advanced ground-based detectors, such as LIGO II. If the instability occurs in a supermassive star, it may produce gravitational radiation detectable by the space-based detector LISA.

Thermal emission from a rotating, supermassive star will cause the configuration to contract slowly and spin up. If internal viscosity and magnetic fields are sufficiently weak, the contracting star will rotate differentially. For each of the six initial angular momentum distributions considered, a cooling supermassive star will likely encounter the dynamical bar mode instability, which may trigger the growth of nonaxisymmetric bars. This scenario will lead to the generation of long wavelength gravitational waves, which could be detectable by future space-based laser interferometers like LISA.

The existence of supermassive black holes is supported by a growing body of observations. Supermassive black holes and their formation events are likely candidates for detection by proposed long-wavelength, space-based gravitational wave interferometers like LISA. However, the nature of the progenitors of supermassive black holes is rather uncertain. Supermassive black hole formation scenarios that involve either the stellar dynamical evolution of dense clusters or the hydrodynamical evolution of supermassive stars have been proposed. Each of these formation scenarios is reviewed and the evolution of supermassive stars is then examined in some detail. Supermassive stars that rotate uniformly during their secular cooling phase will spin up to the mass-shedding limit and eventually contract to the point of relativistic collapse. Supermassive stars that rotate differentially as they cool will likely encounter the dynamical bar mode instability prior to the onset of relativistic collapse. A supermassive star that undergoes this bar distortion, prior to or during collapse, may be a strong source of quasiperiodic, long-wavelength gravitational radiation.

Adaptive techniques are crucial for successful numerical modeling of gravitational waves from astrophysical sources such as coalescing compact binaries, since the radiation typically has wavelengths much larger than the scale of the sources. We have carried out an important step toward this goal, the evolution of weak gravitational waves using adaptive mesh refinement in the Einstein equations. The 2-level adaptive simulation is compared with unigrid runs at coarse and fine resolution, and is shown to track closely the features of the fine grid run.

Compact astrophysical objects that rotate rapidly may encounter the dynamical ``bar instability.'' The bar-like deformation induced by this rotational instability causes the object to become a potentially strong source of gravitational radiation. We have carried out a set of long-duration simulations of the bar instability with two Eulerian hydrodynamics codes. Our results indicate that the remnant of this instability is a persistent bar-like structure that emits a long-lived gravitational radiation signal.

The 3-level leapfrog time integration algorithm is an attractive choice for numerical relativity simulations since it is time-symmetric and avoids non-physical damping. In Newtonian problems without velocity dependent forces, this method enjoys the advantage of long term stability. However, for more general differential equations, whether ordinary or partial, delayed onset numerical instabilities can arise and destroy the solution. A known cure for such instabilities appears to have been overlooked in many application areas. We give an improved cure ("deloused leapfrog") that both reduces memory demands (important for 3+1 dimensional wave equations) and allows for the use of adaptive timesteps without a loss in accuracy. We show both that the instability arises and that the cure we propose works in highly relativistic problems such as tightly bound geodesics, spatially homogeneous spacetimes, and strong gravitational waves. In the gravitational wave test case (polarized waves in a Gowdy spacetime) the deloused leapfrog method was five to eight times less CPU costly at various accuracies than the implicit Crank-Nicholson method, which is not subject to this instability.

Warped H I gas layers in the outer regions of spiral galaxies usually display a noticeably twisted structure. This structure almost certainly arises primarily as a result of differential precession in the H I disk as it settles toward a preferred orientation in an underlying dark halo potential well that is not spherically symmetric. In an attempt to better understand the structure and evolution of these twisted, warped disk structures, we have adopted the "twist-equation" formalism originally developed by Petterson (1977) to study accretion onto compact objects. Utilizing more recent treatments of this formalism, we have generalized the twist-equation to allow for the treatment of non-Keplerian disks and from it have derived a steady-state structure of twisted disks that develops from free precession in a nonspherical, logarithmic halo potential. We have used this steady-state solution to produce H I maps of five galaxies (M83, NGC 300, NGC 2841, NGC 5033, NGC5055), which match the general features of the observed maps of these galaxies quite well. In addition, the model provides an avenue through which the kinematical viscosity of the H I disk and the quadrupole distortion of the dark halo in each galaxy can be quantified. This generalized equation can also be used to examine the time-evolutionary behavior of warped galaxy disks.

The orbital separation of compact binary stars will shrink with time due to the emission of gravitational radiation. This inspiralling phase of a binary system's evolution generally will be very long compared to the system's orbital period, but the final coalescence may be dynamical and driven to a large degree by hydrodynamic effects, particularly if there is a critical separation at which the system becomes dynamically unstable toward merger. Indeed, if weakly relativistic systems (such as white dwarf-white dwarf binaries) encounter a point of dynamical instability at some critically close separation, coalescence may be entirely a classical, hydrodynamic process. Therefore, a proper investigation of this stage of binary evolution must include three-dimensional hydrodynamic simulations. We have constructed equilibrium sequences of synchronously rotating, equal-mass binaries in circular orbit with a single parameter - the binary separation - varying along each sequence. Sequences have been constructed with various polytropic as well as realistic white dwarf and neutron star equations of state. Using a Newtonian, finite-difference hydrodynamics code, we have examined the dynamical stability of individual models along these equilibrium sequences. Our simulations indicate that no points of instability exist on the sequences we analyzed that had relatively soft equations of state (polytropic sequences with polytropic index $n=1.0$ and 1.5 and two white dwarf sequences). However, we did identify dynamically unstable binary models on sequences with stiffer equations of state ($n=0.5$ polytropic sequence and two neutron star sequences). We thus infer that binary systems with soft equations of state are not driven to merger by a dynamical instability.

Laboratory searches for the detection of gravitational waves have focused on the detection of burst signals emitted during a supernova explosion, but have not resulted in any confirmed detections. An alternative approach has been to search for continuous wave (CW) gravitational radiation from the Crab pulsar. In this paper, we examine the possibility of detecting CW gravitational radiation from pulsars and show that nearby millisecond pulsars are generally much better candidates. We show that the minimum strain h_c ~ 10E-26 that can be detected by tuning an antenna to the frequency of the milli- second pulsar PSR 1957+20, with presently available detector technology, is orders of magnitude better than what has been accomplished so far by observing the Crab pulsar, and within an order of magnitude of the maximum strain that may be produced by it. In addition, we point out that there is likely to be a population of rapidly rotating neutron stars (not necessarily radio pulsars) in the solar neighborhood whose spindown evolution is driven by gravitational radiation. We argue that the projected sensitivity of modern resonant detectors is sufficient to detect the subset of this population that lies within 0.1 kpc of the sun.