Dynamical instabilities in protoneutron stars may produce gravitational waves whose observation could shed light on the physics of core-collapse supernovae. When born with sufficient differential rotation, these stars are susceptible to a shear insta
bility (the low-T/|W| instability), but such rotation can also amplify magnetic fields to strengths where they have a considerable impact on the dynamics of the stellar matter. Using a new magnetohydrodynamics module for the Spectral Einstein Code, we have simulated a differentially-rotating neutron star in full 3D to study the effects of magnetic fields on this instability. Though strong toroidal fields were predicted to suppress the low-T/|W| instability, we find that they do so only in a small range of field strengths. Below 4e13 G, poloidal seed fields do not wind up fast enough to have an effect before the instability saturates, while above 5e14 G, magnetic instabilities can actually amplify a global quadrupole mode (this threshold may be even lower in reality, as small-scale magnetic instabilities remain difficult to resolve numerically). Thus, the prospects for observing gravitational waves from such systems are not in fact diminished over most of the magnetic parameter space. Additionally, we report that the detailed development of the low-T/|W| instability, including its growth rate, depends strongly on the particular numerical methods used. The high-order methods we employ suggest that growth might be considerably slower than found in some previous simulations.
We present a first exploration of the results of neutron star-black hole mergers using black hole masses in the most likely range of $7M_odot-10M_odot$, a neutrino leakage scheme, and a modeling of the neutron star material through a finite-temperatu
re nuclear-theory based equation of state. In the range of black hole spins in which the neutron star is tidally disrupted ($chi_{rm BH}gtrsim 0.7$), we show that the merger consistently produces large amounts of cool ($Tlesssim 1,{rm MeV}$), unbound, neutron-rich material ($M_{rm ej}sim 0.05M_odot-0.20M_odot$). A comparable amount of bound matter is initially divided between a hot disk ($T_{rm max}sim 15,{rm MeV}$) with typical neutrino luminosity $L_ usim 10^{53},{rm erg/s}$, and a cooler tidal tail. After a short period of rapid protonization of the disk lasting $sim 10,{rm ms}$, the accretion disk cools down under the combined effects of the fall-back of cool material from the tail, continued accretion of the hottest material onto the black hole, and neutrino emission. As the temperature decreases, the disk progressively becomes more neutron-rich, with dimmer neutrino emission. This cooling process should stop once the viscous heating in the disk (not included in our simulations) balances the cooling. These mergers of neutron star-black hole binaries with black hole masses $M_{rm BH}sim 7M_odot-10M_odot$ and black hole spins high enough for the neutron star to disrupt provide promising candidates for the production of short gamma-ray bursts, of bright infrared post-merger signals due to the radioactive decay of unbound material, and of large amounts of r-process nuclei.
Neutrino emission significantly affects the evolution of the accretion tori formed in black hole-neutron star mergers. It removes energy from the disk, alters its composition, and provides a potential power source for a gamma-ray burst. To study thes
e effects, simulations in general relativity with a hot microphysical equation of state and neutrino feedback are needed. We present the first such simulation, using a neutrino leakage scheme for cooling to capture the most essential effects and considering a moderate mass (1.4 M_{odot} neutron star, 5.6 M_{odot} black hole), high spin (black hole J/M^2=0.9) system with the K_0=220 MeV Lattimer-Swesty equation of state. We find that about 0.08 M_{odot} of nuclear matter is ejected from the system, while another 0.3 M_{odot} forms a hot, compact accretion disk. The primary effects of the escaping neutrinos are (i) to make the disk much denser and more compact, (ii) to cause the average electron fraction Y_e of the disk to rise to about 0.2 and then gradually decrease again, and (iii) to gradually cool the disk. The disk is initially hot (T~6 MeV) and luminous in neutrinos (L_{ u}~10^{54} erg s^{-1}), but the neutrino luminosity decreases by an order of magnitude over 50 ms of post-merger evolution.
Black hole-neutron star mergers resulting in the disruption of the neutron star and the formation of an accretion disk and/or the ejection of unbound material are prime candidates for the joint detection of gravitational-wave and electromagnetic sign
als when the next generation of gravitational-wave detectors comes online. However, the disruption of the neutron star and the properties of the post-merger remnant are very sensitive to the parameters of the binary. In this paper, we study the impact of the radius of the neutron star and the alignment of the black hole spin for systems within the range of mass ratio currently deemed most likely for field binaries (M_BH ~ 7 M_NS) and for black hole spins large enough for the neutron star to disrupt (J/M^2=0.9). We find that: (i) In this regime, the merger is particularly sensitive to the radius of the neutron star, with remnant masses varying from 0.3M_NS to 0.1M_NS for changes of only 2 km in the NS radius; (ii) 0.01-0.05M_sun of unbound material can be ejected with kinetic energy >10^51 ergs, a significant increase compared to low mass ratio, low spin binaries. This ejecta could power detectable optical and radio afterglows. (iii) Only a small fraction (<3%) of the Advanced LIGO events in this parameter range have gravitational-wave signals which could offer constraints on the equation of state of the neutron star. (iv) A misaligned black hole spin works against disk formation, with less neutron star material remaining outside of the black hole after merger, and a larger fraction of that material remaining in the tidal tail instead of the forming accretion disk. (v) Large kicks (v>300 km/s) can be given to the final black hole as a result of a precessing BHNS merger, when the disruption of the neutron star occurs just outside or within the innermost stable spherical orbit.
We present a code for solving the coupled Einstein-hydrodynamics equations to evolve relativistic, self-gravitating fluids. The Einstein field equations are solved in generalized harmonic coordinates on one grid using pseudospectral methods, while th
e fluids are evolved on another grid using shock-capturing finite difference or finite volume techniques. We show that the code accurately evolves equilibrium stars and accretion flows. Then we simulate an equal-mass nonspinning black hole-neutron star binary, evolving through the final four orbits of inspiral, through the merger, to the final stationary black hole. The gravitational waveform can be reliably extracted from the simulation.
