No Arabic abstract
We continue our investigations of the magnetorotational collapse of stellar cores discussing simulations performed with a modified Newtonian gravitational potential that mimics general relativistic effects. The approximate TOV potential used in our simulations catches several features of fully relativistic simulations quite well. It is able to correctly reproduce the behavior of models which show a qualitative change both of the dynamics and the gravitational wave signal when switching from Newtonian to fully relativistic simulations. If this is not the case, the Newtonian and the approximate TOV models differ quantitatively. The collapse proceeds to higher densities with the approximate TOV potential allowing for a more efficient amplification of the magnetic field by differential rotation. Sufficiently strong magnetic fields brake down the cores rotation and trigger a contraction phase to higher densities. Several models exhibit two different kinds of shock generation. Due to magnetic braking, a first shock wave created during the initial centrifugal bounce does not suffice to eject any mass, and the core continues to collapse to supranuclear densities. Another stronger shock wave is generated during the second bounce as the core exceeds nuclear matter density. The gravitational wave signal of these models does not fit into the standard classification. Instead it belongs to the signal type IV introduced by us in the first paper of this series. This signal type is more frequent for the approximate relativistic potential than for the Newtonian one. Strongly magnetized models emit a substantial fraction of their GW power at very low frequencies. A flat spectrum between 10 Hz and > 100 kHz denotes the generation of a jet-like outflow. [Abstract abbreviated]
We study the three-dimensional (3D) hydrodynamics of the post-core-bounce phase of the collapse of a 27-solar-mass star and pay special attention to the development of the standing accretion shock instability (SASI) and neutrino-driven convection. To this end, we perform 3D general-relativistic simulations with a 3-species neutrino leakage scheme. The leakage scheme captures the essential aspects of neutrino cooling, heating, and lepton number exchange as predicted by radiation-hydrodynamics simulations. The 27-solar-mass progenitor was studied in 2D by B. Mueller et al. (ApJ 761:72, 2012), who observed strong growth of the SASI while neutrino-driven convection was suppressed. In our 3D simulations, neutrino-driven convection grows from numerical perturbations imposed by our Cartesian grid. It becomes the dominant instability and leads to large-scale non-oscillatory deformations of the shock front. These will result in strongly aspherical explosions without the need for large-scale SASI shock oscillations. Low-l-mode SASI oscillations are present in our models, but saturate at small amplitudes that decrease with increasing neutrino heating and vigor of convection. Our results, in agreement with simpler 3D Newtonian simulations, suggest that once neutrino-driven convection is started, it is likely to become the dominant instability in 3D. Whether it is the primary instability after bounce will ultimately depend on the physical seed perturbations present in the cores of massive stars. The gravitational wave signal, which we extract and analyze for the first time from 3D general-relativistic models, will serve as an observational probe of the postbounce dynamics and, in combination with neutrinos, may allow us to determine the primary hydrodynamic instability.
(Abridged.) The accretion-induced collapse (AIC) of a white dwarf (WD) may lead to the formation of a protoneutron star and a collapse-driven supernova explosion. This process represents a path alternative to thermonuclear disruption of accreting white dwarfs in Type Ia supernovae. Neutrino and gravitational-wave (GW) observations may provide crucial information necessary to reveal a potential AIC. Motivated by the need for systematic predictions of the GW signature of AIC, we present results from an extensive set of general-relativistic AIC simulations using a microphysical finite-temperature equation of state and an approximate treatment of deleptonization during collapse. Investigating a set of 114 progenitor models in rotational equilibrium, with a wide range of rotational configurations, temperatures and central densities, we extend previous Newtonian studies and find that the GW signal has a generic shape akin to what is known as a Type III signal in the literature. We discuss the detectability of the emitted GWs, showing that the signal-to-noise ratio for current or next-generation interferometer detectors could be high enough to detect such events in our Galaxy. Some of our AIC models form massive quasi-Keplerian accretion disks after bounce. In rapidly differentially rotating models, the disk mass can be as large as ~0.8-Msun. Slowly and/or uniformly rotating models produce much smaller disks. Finally, we find that the postbounce cores of rapidly spinning white dwarfs can reach sufficiently rapid rotation to develop a nonaxisymmetric rotational instability.
This paper presents results from axisymmetric simulations of magneto-rotational stellar core collapse to neutron stars in general relativity using the passive field approximation for the magnetic field. These simulations are performed using a new general relativistic numerical code specifically designed to study this astrophysical scenario. The code is based on the conformally-flat approximation of Einsteins field equations and conservative formulations of the magneto-hydrodynamics equations. The code has been recently upgraded to incorporate a tabulated, microphysical equation of state and an approximate deleptonization scheme. This allows us to perform the most realistic simulations of magneto-rotational core collapse to date, which are compared with simulations employing a simplified (hybrid) equation of state, widely used in the relativistic core collapse community. Furthermore, state-of-the-art (unmagnetized) initial models from stellar evolution are used. In general, stellar evolution models predict weak magnetic fields in the progenitors, which justifies our simplification of performing the computations under the approach that we call the passive field approximation for the magnetic field. Our results show that for the core collapse models with microphysics the saturation of the magnetic field cannot be reached within dynamical time scales by winding up the poloidal magnetic field into a toroidal one. We estimate the effect of other amplification mechanisms including the magneto-rotational instability (MRI) and several types of dynamos.
We apply our recently developed code for spherically symmetric, fully general relativistic (GR) Lagrangian hydrodynamics and multigroup flux-limited diffusion neutrino transport to examine the effects of GR on the hydrodynamics and transport during collapse, bounce, and the critical shock reheating phase of core collapse supernovae. Comparisons of models computed with GR versus Newtonian hydrodynamics show that collapse to bounce takes slightly less time in the GR limit, and that the shock propagates slightly farther out in radius before receding. After a secondary quasistatic rise in the shock radius, the shock radius declines considerably more rapidly in the GR simulations than in the corresponding Newtonian simulations. During the shock reheating phase, core collapse computed with GR hydrodynamics results in a substantially more compact structure from the center out to the stagnated shock. The inflow speed of material behind the shock is also increased. Comparisons also show that the luminosity and rms energy of any neutrino flavor during the shock reheating phase increases when switching from Newtonian to GR hydrodynamics, and decreases when switching from Newtonian to GR transport. This latter decrease in neutrino luminosities and rms energies is less in magnitude than the increase that arise when switching from Newtonian to GR hydrodynamics, with the result that a fully GR simulation gives higher neutrino luminosities and harder neutrino spectra than a fully Newtonian simulation of the same precollapse model.
We report on a set of long-term general-relativistic three-dimensional (3D) multi-group (energy-dependent) neutrino-radiation hydrodynamics simulations of core-collapse supernovae. We employ a full 3D two-moment scheme with the local M1 closure, three neutrino species, and 12 energy groups per species. With this, we follow the post-core-bounce evolution of the core of a nonrotating $27$-$M_odot$ progenitor in full unconstrained 3D and in octant symmetry for $gtrsim$$ 380,mathrm{ms}$. We find the development of an asymmetric runaway explosion in our unconstrained simulation. We test the resolution dependence of our results and, in agreement with previous work, find that low resolution artificially aids explosion and leads to an earlier runaway expansion of the shock. At low resolution, the octant and full 3D dynamics are qualitatively very similar, but at high resolution, only the full 3D simulation exhibits the onset of explosion.