Relativistic simulations of the phase-transition-induced collapse of neutron stars

An increase in the central density of a neutron star may trigger a phase transition from hadronic matter to deconfined quark matter in the core, causing it to collapse to a more compact hybrid-star configuration. We present a study of this, building on previous work by Lin et al. (2006). We follow them in considering a supersonic phase transition and using a simplified equation of state, but our calculations are general relativistic (using 2D simulations in the conformally flat approximation) as compared with their 3D Newtonian treatment. We also improved the treatment of the initial phase transformation, avoiding the introduction of artificial convection. As before, we find that the emitted gravitational-wave spectrum is dominated by the fundamental quasi-radial and quadrupolar pulsation modes but the strain amplitudes are much smaller than suggested previously, which is disappointing for the detection prospects. However, we see significantly smaller damping and observe a nonlinear mode resonance which substantially enhances the emission in some cases. We explain the damping mechanisms operating, giving a different view from the previous work. Finally, we discuss the detectability of the gravitational waves, showing that the signal-to-noise ratio for current or second generation interferometers could be high enough to detect such events in our Galaxy, although third generation detectors would be needed to observe them out to the Virgo cluster, which would be necessary for having a reasonable event rate.Comment: 28 pages, 27 figures. Minor changes to be consistent with published versio

The runaway instability in general relativistic accretion disks

When an accretion disk falls prey to the runaway instability, a large portion of its mass is devoured by the black hole within a few dynamical times. Despite decades of effort, it is still unclear under what conditions such an instability can occur. The technically most advanced relativistic simulations to date were unable to find a clear sign for the onset of the instability. In this work, we present three-dimensional relativistic hydrodynamics simulations of accretion disks around black holes in dynamical space-time. We focus on the configurations that are expected to be particularly prone to the development of this instability. We demonstrate, for the first time, that the fully self-consistent general relativistic evolution does indeed produce a runaway instability.Comment: 5 pages, 3 figures, minor corrections to match published version in MNRAS, +link to animatio

Evolution of open clusters with or without black holes

Binary black holes (BHs) can be formed dynamically in the centers of star clusters. The high natal kicks for stellar-mass BHs used in previous works made it hard to retain BHs in star clusters. Recent studies of massive star evolution and supernovae (SN) propose kick velocities that are lower due to the fallback of the SN ejecta. We study the impact of these updates by performing $N$-body simulations following instantaneous gas expulsion. For comparison, we simulate two additional model sets with the previous treatment of stars: one with high kicks and another with artificial removal of the kicks. Our model clusters initially consist of about one hundred thousand stars, formed with centrally-peaked efficiency. We find that the updated treatment of stars, due to the fallback-scaled lower natal kicks, allows clusters to retain SN remnants after violent relaxation. The mass contribution of the retained remnants does not exceed a few percent of the total bound cluster mass during the early evolution. For this reason, the first giga year of evolution is not affected significantly by this effect. Nevertheless, during the subsequent long-term evolution, the retained BHs accelerate mass segregation, leading to the faster dissolution of the clusters.Comment: 8 pages, 7 figures, submitted to MNRAS, comments are welcom

Correlated Gravitational Wave and Neutrino Signals from General-Relativistic Rapidly Rotating Iron Core Collapse

We present results from a new set of 3D general-relativistic hydrodynamic simulations of rotating iron core collapse. We assume octant symmetry and focus on axisymmetric collapse, bounce, the early postbounce evolution, and the associated gravitational wave (GW) and neutrino signals. We employ a finite-temperature nuclear equation of state, parameterized electron capture in the collapse phase, and a multi-species neutrino leakage scheme after bounce. The latter captures the important effects of deleptonization, neutrino cooling and heating and enables approximate predictions for the neutrino luminosities in the early evolution after core bounce. We consider 12-solar-mass and 40-solar-mass presupernova models and systematically study the effects of (i) rotation, (ii) progenitor structure, and (iii) postbounce neutrino leakage on dynamics, GW, and, neutrino signals. We demonstrate, that the GW signal of rapidly rotating core collapse is practically independent of progenitor mass and precollapse structure. Moreover, we show that the effects of neutrino leakage on the GW signal are strong only in nonrotating or slowly rotating models in which GW emission is not dominated by inner core dynamics. In rapidly rotating cores, core bounce of the centrifugally-deformed inner core excites the fundamental quadrupole pulsation mode of the nascent protoneutron star. The ensuing global oscillations (f~700-800 Hz) lead to pronounced oscillations in the GW signal and correlated strong variations in the rising luminosities of antineutrino and heavy-lepton neutrinos. We find these features in cores that collapse to protoneutron stars with spin periods <~ 2.5 ms and rotational energies sufficient to drive hyper-energetic core-collapse supernova explosions. Hence, joint GW + neutrino observations of a core collapse event could deliver strong evidence for or against rapid core rotation. [abridged]Comment: 29 pages, 14 figures. Replaced with version matching published versio

On the Shear Instability in Relativistic Neutron Stars

We present new results on instabilities in rapidly and differentially rotating neutron stars. We model the stars in full general relativity and describe the stellar matter adopting a cold realistic equation of state based on the unified SLy prescription. We provide evidence that rapidly and differentially rotating stars that are below the expected threshold for the dynamical bar-mode instability, beta_c = T/|W| ~ 0.25, do nevertheless develop a shear instability on a dynamical timescale and for a wide range of values of beta. This class of instability, which has so far been found only for small values of beta and with very small growth rates, is therefore more generic than previously found and potentially more effective in producing strong sources of gravitational waves. Overall, our findings support the phenomenological predictions made by Watts, Andersson and Jones on the nature of the low-T/|W|.Comment: 20 pages; accepted to the Classical and Quantum Gravity special issue for MICRA200

Searching for prompt signatures of nearby core-collapse supernovae by a joint analysis of neutrino and gravitational-wave data

We discuss the science motivations and prospects for a joint analysis of gravitational-wave (GW) and low-energy neutrino data to search for prompt signals from nearby supernovae (SNe). Both gravitational-wave and low-energy neutrinos are expected to be produced in the innermost region of a core-collapse supernova, and a search for coincident signals would probe the processes which power a supernova explosion. It is estimated that the current generation of neutrino and gravitational-wave detectors would be sensitive to Galactic core-collapse supernovae, and would also be able to detect electromagnetically dark SNe. A joint GW-neutrino search would enable improvements to searches by way of lower detection thresholds, larger distance range, better live-time coverage by a network of GW and neutrino detectors, and increased significance of candidate detections. A close collaboration between the GW and neutrino communities for such a search will thus go far toward realizing a much sought-after astrophysics goal of detecting the next nearby supernova.Comment: 10 pages, 3 figures. To appear in Class. Quantum Gra

Stability of general-relativistic accretion disks

Self-gravitating relativistic disks around black holes can form as transient structures in a number of astrophysical scenarios such as binary neutron star and black hole-neutron star coalescences, as well as the core-collapse of massive stars. We explore the stability of such disks against runaway and non-axisymmetric instabilities using three-dimensional hydrodynamics simulations in full general relativity using the THOR code. We model the disk matter using the ideal fluid approximation with a $\Gamma$-law equation of state with $\Gamma=4/3$. We explore three disk models around non-rotating black holes with disk-to-black hole mass ratios of 0.24, 0.17 and 0.11. Due to metric blending in our initial data, all of our initial models contain an initial axisymmetric perturbation which induces radial disk oscillations. Despite these oscillations, our models do not develop the runaway instability during the first several orbital periods. Instead, all of the models develop unstable non-axisymmetric modes on a dynamical timescale. We observe two distinct types of instabilities: the Papaloizou-Pringle and the so-called intermediate type instabilities. The development of the non-axisymmetric mode with azimuthal number m = 1 is accompanied by an outspiraling motion of the black hole, which significantly amplifies the growth rate of the m = 1 mode in some cases. Overall, our simulations show that the properties of the unstable non-axisymmetric modes in our disk models are qualitatively similar to those in Newtonian theory.Comment: 30 pages, 21 figure

Nickel-Rich Outflows Produced by the Accretion-Induced Collapse of White Dwarfs: Lightcurves and Spectra

The accretion-induced collapse (AIC) of a white dwarf to form a neutron star can leave behind a rotationally supported disk with mass of up to ~ 0.1 M_sun. The disk is initially composed of free nucleons but as it accretes and spreads to larger radii, the free nucleons recombine to form helium, releasing sufficient energy to unbind the remaining disk. Most of the ejected mass fuses to form Ni56 and other iron group elements. We present spherically symmetric radiative transfer calculations of the transient powered by the radioactive heating of this ejecta. For an ejecta mass of 1e-2 M_sun (3e-3 M_sun), the lightcurve peaks after <~ 1 day with a peak bolometric luminosity ~ 2e41 erg/s (~ 5e40 erg/s), i.e., a "kilonova"; the decay time is ~ 4 (2) days. Overall, the spectra redden with time reaching U-V ~ 4 after ~ 1 day; the optical colors (B-V) are, however, somewhat blue. Near the peak in the lightcurve, the spectra are dominated by Doppler broadened Nickel features, with no distinct spectral lines present. At ~ 3-5 days, strong Calcium lines are present in the infrared, although the Calcium mass fraction is only ~ 1e-4.5. If rotationally supported disks are a common byproduct of AIC, current and upcoming transient surveys such as the Palomar Transient Factory should detect a few AIC per year for an AIC rate of ~ 1e-2 of the Type Ia rate. We discuss ways of distinguishing AIC from other rapid, faint transients, including .Ia's and the ejecta from binary neutron star mergers.Comment: 9 pages, 4 figures, 1 table, now accepted to MNRA

Crucial Physical Dependencies of the Core-Collapse Supernova Mechanism

We explore with self-consistent 2D F{\sc{ornax}} simulations the dependence of the outcome of collapse on many-body corrections to neutrino-nucleon cross sections, the nucleon-nucleon bremsstrahlung rate, electron capture on heavy nuclei, pre-collapse seed perturbations, and inelastic neutrino-electron and neutrino-nucleon scattering. Importantly, proximity to criticality amplifies the role of even small changes in the neutrino-matter couplings, and such changes can together add to produce outsized effects. When close to the critical condition the cumulative result of a few small effects (including seeds) that individually have only modest consequence can convert an anemic into a robust explosion, or even a dud into a blast. Such sensitivity is not seen in one dimension and may explain the apparent heterogeneity in the outcomes of detailed simulations performed internationally. A natural conclusion is that the different groups collectively are closer to a realistic understanding of the mechanism of core-collapse supernovae than might have seemed apparent.Comment: 25 pages; 10 figure