Search

We present preliminary results from Smooth Particle Hydrodynamics (SPH) simulations of common envelope evolution. We qualitatively compare the interaction between a 0.9M red giant with two different companion masses: a 0.05M brown dwarf and a 0.25M white dwarf companion., Comment: 4 pages, 3 figures, Proceedings of the "Hydrogen Deficient Stars" meeting, Tuebingen, Germany, 2007

The currently-favored model for long-duration gamma-ray bursts (GRBs) invokes explosions from the collapse of a massive star down to a black hole: either directly or through fallback. Those GRBs forming via fallback will produce much less radioactive nickel, and hence it has been argued (without any real calculation) that these systems produce dim supernovae. These fallback black-hole GRBs have been recently been argued as possible progenitors of a newly discovered set of GRBs lacking any associated supernovae. Here we present the first ever radiation-hydrodynamics calculations of the light-curves produced in the hypernova explosion by a delayed-fallback gamma-ray burst. We find that the bolometric light-curve is dominated by shock-deposited energy, not the decay of radioactive elements. As such, observations of such bursts actually probe the density in the progenitor wind more than it does the production of radioactive nickel., Comment: 11 pages (including 3 figures), submitted to ApJ, comments welcome

We use population synthesis to construct distributions of gamma-ray bursts (GRBs) for different proposed progenitor models. We use a description of star formation that takes into account the evolution of metallicity with redshift and galaxy mass, the evolution of galaxy mass with redshift, and the star formation rate with galaxy mass and redshift. We compare predicted distributions with redshift and metallicity to observations of GRB host galaxies and find that the the simple models cannot produce the observed distributions, but that current theoretical models can reproduce the observations within some constraints on the fraction of fallback black holes that produce GRBs., Comment: 17 pages, 9 figures, submitted to the Astrophysical Journal

Theoretical nucleosynthetic yields from supernovae are sensitive to both the details of the progenitor star and the explosion calculation. We attempt to comprehensively identify the sources of uncertainties in these yields. In this paper we concentrate on the variations in yields from a single progenitor arising from common 1-dimensional methods of approximating a supernova explosion. Subsequent papers will examine 3-dimensional effects in the explosion and the progenitor, and trends in mass and composition. For the 1-dimensional explosions we find that both elemental and isotopic yields for Si and heavier elements are a sensitive function of explosion energy. Also, piston-driven and thermal bomb type explosions have different yields for the same explosion energy. Yields derived from 1-dimensional explosions are non-unique., Comment: 87 pages, 22 figures, subitted to the Astrophysical Journal

We study the evolution of first star (Population III) binaries. Under specific conditions, these stars may produce high redshift gamma-ray bursts (GRBs). We demonstrate that the occurrence rate of GRBs does not depend sensitively on evolutionary parameters in the population synthesis models. We show that the first binaries may form a very small group (< 1%) of fast rotating stars through binary tidal interactions that make GRBs. This finding is contrary to the Bromm & Loeb assumption that all stars in close Population III binaries will be spun up by tides and produce a GRB. We find that there is simply not enough fast rotating stars in Population III binaries to expect detection with SWIFT. Predicted detection rates, even with very optimistic assumptions on binary fraction, evolutionary parameters and GRB detection, are very small: 0.1-0.01 per year., Comment: 40 pages, 9 figures, 2 tables: submitted to ApJ

The conditions for the leading r-process site candidate, neutrino-driven winds, can not be reproduced self-consistently in current supernova models. For that reason, we investigate an alternate model involving the mass ejected by fallback in a supernova explosion, through hydrodynamic and nucleosynthesis calculations. The nucleosynthetic products of this ejected material produces r-process elements, including those in the vicinity of the elusive 3rd peak at mass number 195. Trans-iron element production beyond the second peak is made possible by a rapid (<1ms) freezeout of alpha particles which leaves behind a large nucleon (including protons!) to r-process seed ratio. This rapid phase is followed by a relatively long (>15 ms) simmering phase at ~2e9 K, which is the thermodynamic consequence of the hydrodynamic trajectory of the turbulent flows in the fallback outburst. During the slow phase high mass elements beyond the second peak are first made through rapid capture of both protons and neutrons. The flow stays close to valley of stability during this phase. After freeze-out of protons the remaining neutrons cause a shift out to short-lived isotopes as is typical for the r-process. A low electron fraction isn't required in this model, however, the detailed final distribution is sensitive to the electron fraction. Our simulations suggest that supernova fallback is a viable alternative scenario for the r-process., Comment: 12 pages (including 3 figures), submitted to ApJ, comments welcome

We constrain the properties of massive binaries by comparing radial velocity data of Cygnus OB2 with Monte Carlo models. Our comparisons test several popular prescriptions for massive binary parameters. We explore a range of true binary fraction, F, a range of power-law slopes, \alpha, describing the distribution of companion masses, and a range of power-law slopes, \beta, describing the distribution of orbital separations. We also consider distributions of secondary masses described by a Miller-Scalo type initial mass function and by a two-component IMF that includes a substantial ``twin'' population with M_2 ~ M_1. We show that binary fractions F<0.7 are considerably less probable than F>0.8. Thus, the true binary fraction is high. For F=1.0 and a distribution of orbital separations near the canonical Opik's Law distribution (i.e., flat; \beta=0), the power law slope of the mass ratio distribution is \alpha= -0.6 - 0.0. For F~0.8, \alpha is somewhat larger, in the range -0.4 - 1.0. In any case, the secondary star mass function is inconsistent with a Miller-Scalo -like IMF unless the lower end is truncated below ~ 2--4 solar masses. In other words, massive stars preferentially have massive companions. The best fitting models are described by a Salpeter or Miller-Scalo IMF for 60% of secondary star masses with the other 40% of secondaries having M_2 ~ M_1, i.e., ``twins''. These best-fitting model parameters simultaneously predict the fraction of type Ib/c supernovae to be 30-40% of all core-collapse supernovae, in agreement with recent observational estimates., Comment: Replacement for astro-ph/0605069; Substantial changes from original version; 47 pages, 16 figures; Accepted for publication in the Astrophysical Journal

Gamma-ray burst (GRB) afterglow observations have allowed us to significantly constrain the engines producing these energetic explosions. Te redshift and position information provided by these afterglows have already allowed us to limit the progenitors of GRBs to only a few models. The afterglows may also provide another observation that can place further constraints on the GRB progenitor: measurements telling us about the environments surrounding GRBs. Current analyses of GRB afterglows suggest that roughly half of long-duration gamma-ray bursts occur in surroundings with density profiles that are uniform. We study the constraints placed by this observation on both the classic ``collapsar'' massive star progenitor and its relative, the ``helium-merger'' progenitor. We study several aspects of wind mass-loss and find that our modifications to the standard Wolf-Rayet mass-loss paradigm are not sufficient to produce constant density profiles. Although this does not rule out the standard ``collapsar'' progenitor, it does suggest a deficiency with this model. We then focus on the He-merger models and find that such progenitors can fit this particular constraint well. We show how detailed observations can not only determine the correct progenitor for GRBs, but also allow us to study binary evolution physics., Comment: 44 pages including 11 figures

We show that neutrino-driven pulsar kicks can increase the energy of the supernova shock. The observed large velocities of pulsars are believed to originate in the supernova explosion, either from asymmetries in the ejecta or from an anisotropic emission of neutrinos (or other light particles) from the cooling neutron star. In this paper we assume the velocities are caused by anisotropic neutrino emission and study the effects of these neutrino-driven kicks on the supernova explosion. We find that if the collapsed star is marginally unable to produce an explosion, the neutrino-driven mechanisms can drive the convection to make a successful explosion. The resultant explosion is asymmetric, with the strongest ejecta motion roughly in the direction of the neutron star kick. This is in sharp contrast with the ejecta-driven mechanisms, which predict the motion of the ejecta in the opposite direction. This difference can be used to distinguish between the two mechanisms based on the observations of the supernova remnants., Comment: 22 pages including 8 figures, submitted to ApJ, version with high resolution figures can be found at http://qso.lanl.gov/~clf

We compare a suite of 3D explosion calculations and stellar models incorporating advanced physics with observational constraints on the progenitor of Cassiopeia A. We consider binary and single stars from 16 to 40 solar masses with a range of explosion energies and geometries. The parameter space allowed by observations of nitrogen rich high velocity ejecta, ejecta mass, compact remnant mass, and 44Ti and 56Ni abundances individually and as an ensemble is considered. A progenitor of 15-25 solar masses which loses its hydrogen envelope to a binary interaction and undergoes an energetic explosion can match all the observational constraints., Comment: 31 pages, 4 figures, accepted for publication to the Astrophysical Journal

Recent multidimensional hydrodynamic simulations have demonstrated the importance of hydrodynamic motions in the convective boundary and radiative regions of stars to transport of energy, momentum, and composition. The impact of these processes increases with stellar mass. Stellar models which approximate this physics have been tested on several classes of observational problems. In this paper we examine the implications of the improved treatment on supernova progenitors. The improved models predict substantially different interior structures. We present pre-supernova conditions and simple explosion calculations from stellar models with and without the improved mixing treatment at 23 solar masses. The results differ substantially., Comment: 12 pages, 2 figures, accepted for publication in the Astrophysical Journal Letters

An explanation for the large random velocities of pulsars is presented. Like many other models, we propose that the momentum imparted to the star is given at birth. The ultimate source of energy is provided by the intense optically thick neutrino flux that is responsible for radiating the proto-neutron star's gravitational binding energy during the Kelvin-Helmholtz phase. The central feature of the kick mechanism is a radiative-driven magnetoacoustic instability, which we refer to as ``neutrino bubbles.'' Identical in nature to the photon bubble instability, the neutrino bubble instability requires the presence of an equilibrium radiative flux as well as a coherent steady background magnetic field. Over regions of large magnetic flux densities, the neutrino bubble instability is allowed to grow on dynamical timescales ~ 1ms, potentially leading to large luminosity enhancements and density fluctuations. Local luminosity enhancements, which preferentially occur over regions of strong magnetic field, lead to a net global asymmetry in the neutrino emission and the young neutron star is propelled in the direction opposite to these regions. For favorable values of magnetic field structure, size, and strength as well as neutrino bubble saturation amplitude, momentum kicks in excess of 1000 km/s can be achieved. Since the neutrino-powered kick is delivered over the duration of the Kelvin-Helmholtz time ~ a few seconds, one expects spin-kick alignment from this neutrino bubble powered model., Comment: submitted to ApJ

The collapsar model, the now leading model for the engine behind gamma-ray bursts and hypernovae, requires that a star collapses to form a black hole surrounded by an accretion disk of high-angular momentum material. The current best theoretical stellar models, however, do not retain enough angular momentum in the core of the star to make a centrifugally supported disk. In this paper, we present the first calculations of the helium-star/helium-star merger progenitors for the collapsar model. These progenitors invoke the merger of two helium cores during the common envelope inspiral phase of a binary system. We find that, in some cases, the merger can produce cores that are rotating 3-10 times faster than single stars. He-star/He-star gamma-ray burst progenitors have a very different redshift distribution than their single-star gamma-ray burst progenitors and we discuss how gamma-ray burst observations can constrain these progenitors., Comment: 29 pages including 12 figures, submitted to ApJ

The origin and initial mass function of young stellar clusters near the Galactic center are still poorly understood. Two of the more prominent ones, the Arches and Quintuplet clusters, may have formed from a shock-induced burst of star formation, given their similar age and proximity to each other. Their unusual mass distribution, however, may be evidence of a contributing role played by other factors, such as stellar capture from regions outside the clusters themselves. Diffuse X-ray emission from these sources provides us with a valuable, albeit indirect, measure of the stellar mass-loss rate from their constituents. Using recent data acquired with Chandra, we can study the nature and properties of the outflow to not only probe the pertinent physical conditions, such as high metallicity, the magnetic field, and so forth, but also to better constrain the stellar distribution within the clusters, in order to identify their formative history. In this paper, we present a set of three-dimensional smoothed particle hydrodynamics simulations of the wind-wind interactions in both the Arches and Quintuplet clusters. We are guided primarily by the currently known properties of the constituent stars, though we vary the mass-loss rates in order to ascertain the dependence of the measured X-ray flux on the assumed stellar characteristics. Our results are compared with the latest observations of the Arches cluster. Our analysis of the Quintuplet cluster may be used as a basis for comparison with future X-ray observations of this source., Comment: Added comments to sections 2.1 and 2.2, conclusions unchanged; 29 pages, 17 figures, ApJ accepted

The collapse of massive stars not only produces observable outbursts across the entire electromagnetic spectrum but, for Galactic (or near-Galactic) supernovae, detectable signals for ground-based neutrino and gravitational wave detectors. Gravitational waves and neutrinos provide the only means to study the actual engine behind the optical outbursts: the collapsed stellar core. While the neutrinos are most sensitive to details of the equation of state, gravitational waves provide a means to study the mass asymmetries in this central core. We present gravitational wave signals from a series of 3-dimensional core-collapse simulations with asymmetries derived from initial perturbations caused by pre-collapse convection, core rotation, and low-mode convection in the explosion engine itself. A Galactic supernovae will allow us to differentiate these different sources of asymmetry. Combining this signal with other observations of the supernova, from neutrinos to gamma-rays to the compact remnant, dramatically increases the predictive power of the gravitational wave signal. We conclude with a discussion of the gravitational wave signal arising from collapsars, the leading engine for long-duration gamma-ray bursts., Comment: 36 pages including 14 figures, accepted by ApJ

Many neutron stars are observed to be moving with spatial velocities, in excess of 500km/s. A number of mechanisms have been proposed to give neutron stars these high velocities. One of the leading classes of models proposed invokes asymmetries in the core of a massive star just prior to collapse. These asymmetries grow during the collapse, causing the resultant supernova to also be asymmetric. As the ejecta is launched, it pushes off (or ``kicks'') the newly formed neutron star. This paper presents the first 3-dimensional supernova simulations of this process. The ejecta is not the only matter that kicks the newly-formed neutron star. Neutrinos also carry away momentum and the asymmetric collapse leads also to asymmetries in the neutrinos. However, the neutrino asymmetries tend to damp out the neutron star motions and even the most extreme asymmetric collapses presented here do not produce final neutron star velocities above 200km/s., Comment: 7 pages, 4 figures, see http://qso.lanl.gov/~clf/papers/kick.ps.gz for full figures

Recent observations with the Chandra X-ray Observatory have provided us with the capability to discriminate point sources, such as the supermassive black hole Sgr A*, from the diffuse emission within the inner 10'' of the Galaxy. The hot plasma producing the diffuse X-radiation, estimated at ~7.6*10^{31} ergs s^{-1} arcsec^{-2} in the 2-10 keV band, has a RMS electron density ~26 cm^{-3} and a temperature kT ~= 1.3 keV, with a total inferred mass of ~0.1 solar masses. At least some of this gas must be injected into the ISM via stellar winds. In the most recent census, about 25 bright, young stars have been identified as the dominant sources of the overall mass efflux from the Galactic center. In this paper, we use detailed 3-dimensional SPH simulations to study the wind-wind interactions occurring in the inner 3 parsecs of the Galaxy, with a goal of understanding what fraction, if any, of the diffuse X-ray flux measured by Chandra results from the ensuing shock heating of the ambient medium. We conclude that this process alone can account for the entire X-ray flux observed by Chandra in the inner 10'' of the Galaxy. Understanding the X-ray morphology of the environment surrounding Sgr A* will ultimately provide us with a greater precision in modeling the accretion of gas onto this object, which appears to be relatively underluminous compared to its brethren in the nuclei of other galaxies., Comment: 22 pages, 8 figures; new calculations reported, conclusion unchanged; accepted for publication in ApJ

We study the physics behind the neutrino-driven mechanism for gamma-ray bursts and hypernovae, deriving the critical density at which these outbursts occur in the collapsar model. The agreement between this derivation and results from past collapsar simulations (MacFadyen & Woosley 2000) is excellent, implying that we have captured the essential physics. We then use this derivation to study a range of progenitors for collapsar gamma-ray bursts. We derive how much of the star will accrete onto the black hole core before the infall density drops below this critical density, leading to an estimate of the remnant black hole mass for GRBs and hypernovae. We also estimate the time delays between gravity wave or neutrino signals and the onset of the explosion or burst event. This derivation, combined with future observational constraints, provides a physical insight into the structure of the GRB progenitor., Comment: 6 pages total, 3 figures. submitted to ApJ Letters

The new generation of gravitational wave (GW) detectors have the potential to open a novel window onto the violent dynamics of core collapse. Although it is certain that core collapse events generate gravitational radiation, understanding the characteristics of the radiation -- whether it can be measured with these detectors, and the best way to go about doing so -- is a challenging problem. In this chapter we review the promise of GWs as observational probes, including a discussion of the current state of GW detectors, and discuss the status of work to understand the waves generated by stellar core collapse., Comment: 31 pages, 9 figures, to appear in "Stellar Collapse" (Astrophysics and Space Science), edited by C. L. Fryer

Using 220 Gamma-Ray Burst (GRB) redshifts and luminosities derived from the luminosity-variability relationship of Fenimore & Ramirez-Ruiz (2000), we show that there exists a significant correlation between the GRB luminosity and redshift. In particular, we find that the evolution of the average luminosity can be parameterized as L ~ (1+z)^(1.4 +- 0.5), where z is the burst redshift. We discuss the possible reasons behind this evolution and compare it to other known sources that exhibit similar behavior. In addition, we use non-parametric statistical techniques to independently estimate the distributions of the luminosity and redshift of bursts, accounting for the evolution (in contrast to previous studies which have assumed that the luminosity function is independent of redshift). We present these distributions and discuss their implications. Most significantly, we find a co-moving rate density of GRBs that continues to increase to (1+z) ~ 10. From this estimate of the GRB rate density, we then use the population synthesis codes of Fryer et al. (1999) to estimate the star formation rate at high redshifts, based on different progenitor models of GRBs. We find that no matter what the progenitor or population synthesis model, the star formation rate increases or remains constant to very high redshifts (z ~ 10)., Comment: Final version accepted to the Astrophysical Journal

We derive estimates for the characteristics of gravitational radiation from stellar collapse, using recent models of the core-collapse of Chandrasekhar mass white dwarfs (accretion induced collapse), core-collapse supernovae and collapsars, and the collapse of very massive stars (~> 300 Msun). We study gravitational-wave emission mechanisms using several estimation techniques, including two-dimensional numerical computation of quadrupole wave emission, estimates of bar-mode strength, estimates of r-mode emission, and estimates of waves from black hole ringing. We also review the rate at which the relevant collapses are believed to occur, which has a major impact on their relevance as astrophysical sources. Although the latest supernova progenitor simulations produce cores rotating much slower than those used in the past, we find that bar-mode and r-mode instabilities from core-collapse supernovae remain among the leading candidate sources for LIGO-II. Accretion induced collapse (AIC) of a white dwarf could produce gravitational-wave signals similar to those from core-collapse. In the models that we examine, such collapses are not unstable to bar modes; we note that models recently examined by Liu and Lindblom, which have slightly more angular momentum, are certainly unstable to bar formation. Because AIC events are probably 1,000 times less common than core-collapse supernovae, the typical AIC event will be much further away, and thus the observed waves will be much weaker. In the most optimistic circumstances, we find it may be possible to detect gravitational waves from the collapse of 300 Msun Population III stars., Comment: 48 pages including 11 figures, to appear in ApJ. New version incorporates several helpful comments; also corrects some errors in the displayed LIGO noise curve, and several typos

There is growing observational evidence that gamma-ray bursts (GRBs) are powered by black holes accreting rapidly through a disk. The supernova-like outburst that accompanies some gamma-ray bursts suggest that some long-duration GRBs may be driven by the accretion of a rotating stellar core onto a central black hole. Such a system can be produced when a compact remnant spirals into the helium core of its binary companion. During the inspiral, orbital angular momentum is injected into the core. By the time the compact remnant reaches the center of the helium core, it too has gained angular momentum as well as mass, producing a rapidly accreting black hole (or neutron star) at the center of a rotating stellar core. In this paper, we use a 3-dimensional smooth particle hydrodynamics (SPH) code to follow such a merger process and make quantitative estimates of the initial mass and spin of the central compact remnant, as well as the angular momentum in the accreting helium core. From these results, we estimate GRB explosion energies. In all mergers, magnetically driven jets are expected to produce GRB explosions with energies above 10^{51} ergs. For neutrino-annihilation-driven explosions, the GRB energy increases dramatically with helium star mass: the merger of a 2 M_sun compact remnant with a 4 M_sun helium star only produces a 10^{47} ergs explosion in \sim 500 s whereas the merger of a 2 M_sun compact remnant with a 16 M_sun helium star produces a >10^{52} ergs explosion in \sim 65 s. (Abridged), Comment: 28 pages, 12 figures, accepted by ApJ (March 2001 issue)

We derive the theoretical distribution function of black hole masses by studying the formation processes of black holes. We use the results of recent 2D simulations of core-collapse to obtain the relation between remnant and progenitor masses and fold it with an initial mass function for the progenitors. We examine how the calculated black-hole mass distributions are modified by (i) strong wind mass loss at different evolutionary stages of the progenitors, and (ii) the presence of close binary companions to the black-hole progenitors. Thus, we are able to derive the binary black hole mass distribution. The compact remnant distribution is dominated by neutron stars in the mass range 1.2-1.6Msun and falls off exponentially at higher remnant masses. Our results are most sensitive to mass loss from winds which is even more important in close binaries. Wind mass-loss causes the black hole distribution to become flatter and limits the maximum possible black-hole mass (<10-15Msun). We also study the effects of the uncertainties in the explosion and unbinding energies for different progenitors. The distributions are continuous and extend over a broad range. We find no evidence for a gap at low values (3-5Msun) or for a peak at higher values (~7Msun) of black hole masses, but we argue that our black hole mass distribution for binaries is consistent with the current sample of measured black-hole masses in X-ray transients. We discuss possible biases against the detection or formation of X-ray transients with low-mass black holes. We also comment on the possibility of black-hole kicks and their effect on binaries., Comment: 22 pages, submitted to ApJ

A variety of current models for gamma-ray bursts (GRBs) suggest a common engine - a black hole of several solar masses accreting matter from a disk at a rate 0.01 to 10 solar masses per second. Using a numerical model for relativistic disk accretion, we have studied steady-state accretion at these high rates. Inside a radius ~ 10**8 cm, for accretion rates greater than about 0.01 solar masses per second, a global state of balanced power comes to exist between neutrino losses, chiefly pair capture on nucleons, and dissipation. Energy emitted in neutrinos is less, and in the case of low accretion rates, very much less, than the maximum efficiency factor for black hole accretion (0.057 for no rotation; 0.42 for extreme Kerr rotation) times Mdot c**2. The efficiency for producing a pair fireball along the rotational axis by neutrino annihilation is calculated and found to be highly variable and very sensitive to the accretion rate. For some of the higher accretion rates studied, it can be several per cent or more; for accretion rates less than 0.05 solar masses per second, it is essentially zero. The efficiency of the Blandford-Znajek mechanism in extracting rotational energy from the black hole is also estimated. In light of these results, the viability of various gamma-ray burst models is discussed and the sensitivity of the results to disk viscosity, black hole rotation rate, and black hole mass explored. A diverse range of GRB energies seems unavoidable and neutrino annihilation in hyper-accreting black hole systems can explain bursts up to 10**52 erg. Larger energies may be inferred for beaming systems., Comment: 46 pages, includes 9 figures, LaTeX (uses aaspp4.sty), accepted by The Astrophysical Journal. Additional solutions in Tables and Figs. 4 and 5, minor revisions to text, references added

We use a Monte Carlo binary synthesis code to model the formation and evolution of neutron star systems including high-mass X-ray binaries, low-mass X-ray binaries, double neutron star systems and radio pulsars. Our focus is on the signature imprinted on such systems due to natal kicks to neutron stars over and above that imparted by orbital motions. The code incorporates the effect of the galactic potential (including rotation) on the velocities of these systems. A comparison between our models and the observations leads us to infer mean natal kicks between 400-500 km/s. Moreover, to be consistent with all the data, we require a bimodal kick distribution with one peak in the distribution near 0 km/s and the other above 600 km/s., Comment: 41 pages total, 24 text+tables pages, 17 figures, AASTeX, Accepted for publication in ApJ

We study the four double neutron star systems found in the Galactic disk in terms of the orbital characteristics of their immediate progenitors and the natal kicks imparted to neutron stars. Analysis of the effect of the second supernova explosion on the orbital dynamics, combined with recent results from simulations of rapid accretion onto neutron stars lead us to conclude that the observed systems could not have been formed had the explosion been symmetric. Their formation becomes possible if kicks are imparted to the radio-pulsar companions at birth. We identify the constraints imposed on the immediate progenitors of the observed double neutron stars and calculate the ranges within which their binary characteristics (orbital separations and masses of the exploding stars) are restricted. We also study the dependence of these limits on the magnitude of the kick velocity and the time elapsed since the second explosion. For each of the double neutron stars, we derive a minimum kick magnitude required for their formation, and for the two systems in close orbits we find it to exceed 200km/s. Lower limits are also set to the center-of-mass velocities of double neutron stars, and we find them to be consistent with the current proper motion observations., Comment: 25 pages, 6 figs (9 parts), 4 tables, AASTeX, Accepted in ApJ

We present an extensive study of accretion onto neutron stars in which the velocity of the neutron star and structure of the surrounding medium is such that the Bondi-Hoyle accretion exceeds .001 Msun/y. For most cases, hypercritical accretion due to rapid neutrino cooling allows the neutron star to accrete above the Bondi-Hoyle rate as previously pointed out by Chevalier. However, for a subset of simulations which corresponds to evolutionarily common events, convection driven by neutrino heating can lead to explosions by a mechanism similar to that found in core-collapse supernovae. Armed with the results from our calculations, we are in a position to predict the fate of a range of rapid-infall neutron star accretors present in certain low-mass X-ray binaries, common envelope systems, supernova fallbacks and Thorne-Zytkow objects (TZOs). A majority of the common envelope systems that we considered led to explosions expelling the envelope, halting the neutron star's inward spiral, and allowing the formation of close binary systems. Smothered neutron stars produced in collisions also lead to explosions, preventing them from forming millisecond pulsars. For supernovae in which the fallback of material towards the neutron star is large, we find that a black hole is formed within a few seconds. Finally, we argue that the current set of TZO formation scenarios is inadequate and leads instead to hypercritical accretion and black hole formation. Moreover, it appears that many of the current TZ models have structures ill-suited for modeling by mixing length convection. This has prompted us to develop a simple test to determine the viability of this approximation for a variety of convective systems., Comment: 41 pages, 18 figures, for full color version of figure 6 send e-mail to cfryer@as.arizona.edu, to appear in the March 20, 1996 addition of the Astrophysical Journal

Condensed Abstract: We present an extensive study of the inception of supernova explosions by following the evolution of the cores of two massive stars (15 Msun and 25 Msun) in two dimensions. Our calculations begin at the onset of core collapse and stop several 100 ms after the bounce, at which time successful explosions of the appropriate magnitude have been obtained. (...) Guided by our numerical results, we have developed a paradigm for the supernova explosion mechanism. We view a supernova as an open cycle thermodynamic engine in which a reservoir of low-entropy matter (the envelope) is thermally coupled and physically connected to a hot bath (the protoneutron star) by a neutrino flux, and by hydrodynamic instabilities. (...) In essence, a Carnot cycle is established in which convection allows out-of-equilibrium heat transfer mediated by neutrinos to drive low entropy matter to higher entropy and therefore extracts mechanical energy from the heat generated by gravitational collapse. We argue that supernova explosions are nearly guaranteed and self-regulated by the high efficiency of the thermodynamic engine. (...) Convection continues to accumulate energy exterior to the neutron star until a successful explosion has occurred. At this time, the envelope is expelled and therefore uncoupled from the heat source (the neutron star) and the energy input ceases. This paradigm does not invoke new or modified physics over previous treatments, but relies on compellingly straightforward thermodynamic arguments. It provides a robust and self-regulated explosion mechanism to power supernovae which is effective under a wide range of physical parameters., Comment: submitted to ApJ, only full abstract is posted to archive, e-mail mrc@lick.ucsc.edu to obtain complete paper, plain TeX

The high-energy universe has revealed that energetic particles are ubiquitous in the cosmos and play a vital role in the cultivation of cosmic environments on all scales. Our pursuit of more than a century to uncover the origins and fate of these cosmic energetic particles has given rise to some of the most interesting and challenging questions in astrophysics. Within our own galaxy, we have seen that energetic particles engage in a complex interplay with the galactic environment and even drive many of its key characteristics (for more information, see the first white paper in this series). On cosmological scales, the energetic particles supplied by the jets of active galactic nuclei (AGN) are an important source of energy for the intracluster and intergalactic media, providing a mechanism for regulating star formation and black hole growth and cultivating galaxy evolution (AGN feedback). Gamma-ray burst (GRB) afterglows encode information about their circumburst environment, which has implications for massive stellar winds during previous epochs over the stellar lifecycle. As such, GRB afterglows provide a means for studying very high-redshift galaxies since GRBs can be detected even if their host galaxy cannot. It has even been suggest that GRB could be used to measure cosmological distance scales if they could be shown to be standard candles. Though they play a key role in cultivating the cosmological environment and/or enabling our studies of it, there is still much we do not know about AGNs and GRBs, particularly the avenue in which and through which they supply radiation and energetic particles, namely their jets. Despite the enormous progress in particle-in-cell and magnetohydrodynamic simulations, we have yet to pinpoint the processes involved in jet formation and collimation and the conditions under which they can occur. For that matter, we have yet to identify the mechanism(s) through which the jet accelerates energetic particles – is it the commonly invoked diffusive shock acceleration process or is another mechanism, such as magnetic reconnection, required? Do AGNs and GRBs accelerate hadrons, and if so, do they accelerate them to ultra-high energies and are there high-energy neutrinos associated with them? MeV gamma-ray astronomy, enabled by technological advances that will be realized in the coming decade, will provide a unique and indispensable perspective on the persistent mysteries of the energetic universe. This White Paper is the second of a two-part series highlighting the most well-known high-energy cosmic accelerators and contributions that MeV gamma-ray astronomy will bring to understanding their energetic particle phenomena. Specifically, MeV astronomy will: 1. Determine whether AGNs accelerate CRs to ultra-high energies; 2. Provide the missing pieces for the physics of the GRB prompt emission; 3. Measure magnetization in cosmic accelerators and search for acceleration via reconnection.

The Transient Astrophysics Probe (TAP) is a wide-field multi-wavelength transient mission proposed for flight starting in the late 2020s. The mission instruments include unique ``Lobster-eye'' imaging soft X-ray optics that allow an approximately 1600-degrees-squared Field of View (FoV); a high sensitivity, 1-degree-squared FoV soft X-ray telescope; a 1-degree-squared FoV Infrared telescope with bandpass 0.6 to 3 microns; and a set of 8 NaI gamma-ray detectors. TAP's most exciting capability will be the observation of tens per year of X-ray and Infrared counterparts of gravitational waves (GWs) involving stellar-mass black holes and neutron stars detected by LIGO (Laser Interferometer Gravitational-Wave Observatory ) / Virgo / KAGRA (Kamioka (Japan) Gravitational Wave Detector) / LIGO-India, and possibly several per year X-ray counterparts of GWs from supermassive black holes, detected by LISA (Laser Interferometer Space Antenna) and Pulsar Timing Arrays. TAP will also discover hundreds of X-ray transients related to compact objects, including tidal disruption events, supernova shock breakouts, and Gamma-Ray Bursts from the epoch of reionization.