The discovery of Jupiter-mass planets in close orbits about their parent stars has challenged models of planet formation. Recent observations have shown that a number of these planets have highly ...inclined, sometimes retrograde orbits about their parent stars, prompting much speculation as to their origin. It is known that migration alone cannot account for the observed population of these misaligned hot Jupiters, which suggests that dynamical processes after the gas disk dissipates play a substantial role in yielding the observed inclination and eccentricity distributions. One particularly promising candidate is planet-planet scattering, which is not very well understood in the nonlinear regime of tides. Through three-dimensional hydrodynamical simulations of multi-orbit encounters, we show that planets that are scattered into an orbit about their parent stars with closest approach distance being less than approximately three times the tidal radius are either destroyed or completely ejected from the system. We find that as few as 9 and as many as 12 of the currently known hot Jupiters have a maximum initial apastron for scattering that lies well within the ice line, implying that these planets must have migrated either before or after the scattering event that brought them to their current positions. If stellar tides are unimportant (Q * 107), disk migration is required to explain the existence of the hot Jupiters present in these systems. Additionally, we find that the disruption and/or ejection of Jupiter-mass planets deposits a Sun's worth of angular momentum onto the host star. For systems in which planet-planet scattering is common, we predict that planetary hosts have up to a 35% chance of possessing an obliquity relative to the invariable plane of greater than 90?.
We show that the luminosity of a star-forming galaxy is capped by the production and subsequent expulsion of cosmic rays from its interstellar medium. By defining an Eddington luminosity in cosmic ...rays, we show that the star formation rate of a given galaxy is limited by its mass content and the cosmic-ray mean free path. When the cosmic-ray luminosity and pressure reach a critical value as a result of vigorous star formation, hydrostatic balance is lost, a galactic-scale cosmic-ray-driven wind develops, and star formation is choked off. Cosmic-ray pressure driven winds are likely to produce wind velocities in proportion to and significantly in excess of the galactic escape velocity. It is possible that cosmic-ray feedback results in the Faber-Jackson relation for a plausible set of input parameters that describe cosmic-ray production and transport, which are calibrated by observations of the Milky Way's interstellar cosmic rays as well as other nearby galaxies.
ABSTRACT In this paper, we model the observable signatures of tidal disruptions of white dwarf (WD) stars using massive black holes (MBHs) of moderate mass, 103-105 M . When the WD passes deep enough ...within the MBH's tidal field, these signatures include thermonuclear transients from burning during maximum compression. We combine a hydrodynamic simulation that includes nuclear burning of the disruption of a 0.6 M C/O WD with a Monte Carlo radiative transfer calculation to synthesize the properties of a representative transient. The transient's emission emerges in the optical, with light curves and spectra reminiscent of Type I supernovae. The properties are strongly viewing angle dependent, and key spectral signatures are 10,000 km s−1 doppler shifts, due to the orbital motion of the unbound ejecta. Disruptions of He WDs likely produce large quantities of intermediate-mass elements, offering a possible production mechanism for Ca-rich transients. Accompanying multi-wavelength transients are fueled by accretion and arise from the nascent accretion disk and relativistic jet. If MBHs of moderate mass exist with number densities similar to those of supermassive BHs, both high-energy wide-field monitors and upcoming optical surveys should detect tens to hundreds of WD tidal disruptions per year. The current best strategy for their detection may therefore be deep optical follow-up of high-energy transients of unusually long duration. The detection rate or the nondetection of these transients by current and upcoming surveys can thus be used to place meaningful constraints on the extrapolation of the MBH mass function to moderate masses.
We present three-dimensional local hydrodynamic simulations of flows around objects embedded within stellar envelopes using a "wind tunnel" formalism. Our simulations model the common envelope ...dynamical inspiral phase in binary star systems in terms of dimensionless flow characteristics. We present suites of simulations that study the effects of varying the binary mass ratio, stellar structure, equation of state, relative Mach number of the object's motion through the gas, and density gradients across the gravitational focusing scale. For each model, we measure coefficients of accretion and drag experienced by the embedded object. These coefficients regulate the coupled evolution of the object's masses and orbital tightening during the dynamical inspiral phase of the common envelope. We extrapolate our simulation results to accreting black holes with masses comparable to that of the population of LIGO black holes. We demonstrate that the mass and spin accrued by these black holes per unit orbital tightening are directly related to the ratio of accretion to drag coefficients. We thus infer that the mass and dimensionless spin of initially nonrotating black holes change by of order 1% and 0.05, respectively, in a typical example scenario. Our prediction that the masses and spins of black holes remain largely unmodified by a common envelope phase aids in the interpretation of the properties of the growing observed population of merging binary black holes. Even if these black holes passed through a common envelope phase during their assembly, features of mass and spin imparted by previous evolutionary epochs should be preserved.
During a common envelope (CE) episode in a binary system, the engulfed companion spirals to tighter orbital separations under the influence of drag from the surrounding envelope material. As this ...object sweeps through material with a steep radial gradient of density, net angular momentum is introduced into the flow, potentially leading to the formation of an accretion disk. The presence of a disk would have dramatic consequences for the outcome of the interaction because accretion might be accompanied by strong, polar outflows with enough energy to unbind the entire envelope. Without a detailed understanding of the necessary conditions for disk formation during CE, therefore, it is difficult to accurately predict the population of merging compact binaries. This paper examines the conditions for disk formation around objects embedded within CEs using the "wind tunnel" formalism developed by MacLeod et al. We find that the formation of disks is highly dependent on the compressibility of the envelope material. Disks form only in the most compressible of stellar envelope gas, found in envelopes' outer layers in zones of partial ionization. These zones are largest in low-mass stellar envelopes, but comprise small portions of the envelope mass and radius in all cases. We conclude that disk formation and associated accretion feedback in CE is rare, and if it occurs, transitory. The implication for LIGO black hole binary assembly is that by avoiding strong accretion feedback, CE interactions should still result in the substantial orbital tightening needed to produce merging binaries.
We use a simple framework to calculate the time evolution of the composition of the fallback material onto a supermassive black hole arising from the tidal disruption of main-sequence stars. We study ...stars with masses between 0.8 and 3.0 M , at evolutionary stages from zero-age main sequence to terminal-age main sequence, built using the Modules for Experiments in Stellar Astrophysics code. We show that most stars develop enhancements in nitrogen (14N) and depletions in carbon (12C) and oxygen (16O) over their lifetimes, and that these features are more pronounced for higher mass stars. We find that, in an accretion-powered tidal disruption flare, these features become prominent only after the time of peak of the fallback rate and appear at earlier times for stars of increasing mass. We postulate that no severe compositional changes resulting from the fallback material should be expected near peak for a wide range of stellar masses and, as such, are unable to explain the extreme helium-to-hydrogen line ratios observed in some TDEs. On the other hand, the resulting compositional changes could help explain the presence of nitrogen-rich features, which are currently only detected after peak. When combined with the shape of the light curve, the time evolution of the composition of the fallback material provides a clear method to help constrain the nature of the disrupted star. This will enable a better characterization of the event by helping break the degeneracy between the mass of the star and the mass of the black hole when fitting tidal disruption light curves.
A neo-neutron star is a hot neutron star that has just become transparent to neutrinos. In a core-collapse supernova or accretion-induced collapse of a white dwarf, the neo-neutron star phase ...directly follows the proto-neutron star phase, about 30-60 s after the initial collapse. It will also be present in a binary neutron star merger in the case where the "born-again" hot massive compact star does not immediately collapse into a black hole. Eddington or even super-Eddington luminosities are present for some time. A neo-neutron star produced in a core-collapse supernova is not directly observable, but the one produced by a binary merger, likely associated with an off-axis short gamma-ray burst, may be observable for some time as well as when produced in the accretion-induced collapse of a white dwarf. We present a first step in the study of this neo-neutron star phase in a spherically symmetric configuration, thus ignoring fast rotation and also ignoring the effect of strong magnetic fields. We put particular emphasis on determining how long the star can sustain a near-Eddington luminosity and also show the importance of positrons and contraction energy during the neo-neutron star phase. We finally discuss the observational prospects for neutron star mergers triggered by LIGO and for accretion-induced collapse transients.
The frequency of short gamma-ray bursts (GRBs) in galaxies with distinct star formation histories can be used to constrain the lifetime of the progenitor systems. As an Illustration, we consider here ...the constraints that can be derived from separating the host galaxies into early and late types. On average, early-type galaxies have their stars formed earlier than late-type galaxies, and this difference, together with the time delay between progenitor formation and short GRB outburst, leads to different burst rates in the two types of hosts. Presently available data suggest, but do not yet prove, that the local short GRB rate in early-type galaxies may be comparable to that in late-type galaxies. This suggests that, unlike Type Ia supernovae, at least half of the short GRB progenitors that can outburst within a Hubble time have lifetimes greater than about 7 Gyr. Models of the probability distribution of time delays, here parameterized as P( tau ) proportional to tau super(n), with R unk -1 are favored. This apparent long time delay and the fact that early-type galaxies in clusters make a substantial contribution to the local stellar mass inventory can explain the observed preponderance of short GRBs in galaxy clusters.
We present the results of a systematic numerical study of the onset of mass transfer in double degenerate binary systems and its impact on the subsequent evolution. All investigated systems belong to ...the regime of direct impact, unstable mass transfer. In all of the investigated cases, even those considered unstable by conventional stability analysis, we find a long-lived mass transfer phase continuing for as many as several dozen orbital periods. This settles a recent debate sparked by a discrepancy between earlier smoothed particle hydrodynamics (SPH) calculations that showed disruptions after a few orbital periods and newer grid-based studies in which mass transfer continued for tens of orbits. The number of orbits a binary survives sensitively depends on the exact initial conditions. We find that the approximate initial conditions that have been used in most previous SPH calculations have a serious impact on all stages of the evolution from the onset of mass transfer up to the final structure of the remnant. We compare 'approximate' initial conditions where spherical stars are placed at an initial separation obtained from an estimate of the Roche lobe size with 'accurate' initial conditions that were constructed by carefully driving the binary system to equilibrium by a relaxation scheme. Simulations that use the approximate initial conditions underestimate the initial separation when mass transfer sets in, which yields a binary that only survives for only a few orbits and thus a rapidly fading gravitational wave signal. Conversely, the accurate initial conditions produce a binary system in which the mass transfer phase is extended by almost two orders of magnitude in time, resulting in a gravitational wave signal with amplitude and frequency that remain essentially constant up until merger. As we show that these binaries can survive at small separation for hundreds of orbital periods, their associated gravitational wave signal should be included when calculating the gravitational wave foreground (although expected to be below Laser Interferometer Space Antenna's sensitivity at these high frequencies). We also show that the inclusion of the entropy increase associated with shock heating of the accreted material reduces the number of orbits a binary survives given the same initial conditions, although the effect is not as pronounced when using the appropriate initial conditions. The use of accurate initial conditions and a correct treatment of shock heating allows for a reliable time evolution of the temperature, density, and angular momentum, which are important when considering thermonuclear events that may occur during the mass transfer phase and/or after merger. Our treatment allows us to accurately identify when surface detonations may occur in the lead-up to the merger, as well as the properties of final merger products.
ABSTRACT
The discovery of GRB 170817A, the first unambiguous off-axis short gamma-ray burst (sGRB) arising from a neutron star merger, has challenged our understanding of the angular structure of ...relativistic jets. Studies of the jet propagation usually assume that the jet is ejected from the central engine with a top-hat structure and its final structure, which determines the observed light curve and spectra, is primarily regulated by the interaction with the nearby environment. However, jets are expected to be produced with a structure that is more complex than a simple top-hat, as shown by global accretion simulations. We present numerical simulations of sGRBs launched with a wide range of initial structures, durations, and luminosities. We follow the jet interaction with the merger remnant wind and compute its final structure at distances ≳1011 cm from the central engine. We show that the final jet structure, as well as the resulting afterglow emission, depends strongly on the initial structure of the jet, its luminosity, and duration. While the initial structure of the jet is preserved for long-lasting sGRBs, it is strongly modified for jets barely making their way through the wind. This illustrates the importance of combining the results of global simulations with propagation studies in order to better predict the expected afterglow signatures from neutron star mergers. Structured jets provide a reasonable description of the GRB 170817A afterglow emission with an off-axis angle θobs ≈ 22.5°.