Showing 1–50 of 118 results for author: Müller, B

The minimum neutron star mass in neutrino-driven supernova explosions

Authors: Bernhard Müller, Alexander Heger, Jade Powell

Abstract: Supernova theory has struggled to explain the lightest known neutron star candidate with an accurate mass determination, the $1.174\mathrm{M}_\odot$ companion in the eccentric compact binary system J0453+1559. To improve the theoretical lower limit for neutron star birth masses, we perform 3D supernova simulations for five stellar models close to the minimum mass for iron core collapse. We obtain… ▽ More Supernova theory has struggled to explain the lightest known neutron star candidate with an accurate mass determination, the $1.174\mathrm{M}_\odot$ companion in the eccentric compact binary system J0453+1559. To improve the theoretical lower limit for neutron star birth masses, we perform 3D supernova simulations for five stellar models close to the minimum mass for iron core collapse. We obtain a record-low neutron star mass of $1.192\mathrm{M}_\odot$ and a substantial kick of $\mathord{\sim} 100\,\mathrm{km}\,\mathrm{s}^{-1}$. Given residual uncertainties in stellar evolution, a neutron star origin for the $1.174\mathrm{M}_\odot$ object remains plausible. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: 6 pages, 3 figures. Submitted to Physical Review Letters

The gravitational-wave emission from the explosion of a 15 solar mass star with rotation and magnetic fields

Authors: Jade Powell, Bernhard Müller

Abstract: Gravitational waveform predictions from 3D simulations of explosions of non-rotating massive stars with no magnetic fields have been extensively studied. However, the impact of magnetic fields and rotation on the core-collapse supernova gravitational-wave signal is not well understood beyond the core-bounce phase. Therefore, we perform four magnetohydrodynamical simulations of the explosion of a… ▽ More Gravitational waveform predictions from 3D simulations of explosions of non-rotating massive stars with no magnetic fields have been extensively studied. However, the impact of magnetic fields and rotation on the core-collapse supernova gravitational-wave signal is not well understood beyond the core-bounce phase. Therefore, we perform four magnetohydrodynamical simulations of the explosion of a $15\,M_{\odot}$ star with the SFHx and SFHo equations of state. All of the models start with a weak magnetic field strength of $10^{8}$\,G, and two of the models are rapidly rotating. We discuss the impact of the rotation and magnetic fields on the gravitational-wave signals. We find that the weak pre-collapse fields do not have a significant impact on the gravitational-wave signal amplitude. With rapid rotation, the f/g-mode trajectory can change in shape, and the dominant emission band becomes broader. We include the low-frequency memory component of the gravitational-wave signal from both matter motions and neutrino emission anisotropy. We show that including the gravitational waves from anisotropic neutrino emission increases the supernova detection distances for the Einstein Telescope, and would also be detectable out to Mpc distances by a moon-based gravitational-wave detector. △ Less

Submitted 22 July, 2024; v1 submitted 13 June, 2024; originally announced June 2024.

Universality in supernova gravitational waves with proto-neutron star properties

Authors: Hajime Sotani, Bernhard Müller, Tomoya Takiwaki

Abstract: Gravitational wave signals from core-collapse supernovae are one of the important observables for extracting the information of dense matter. To extract the properties of proto-neutron stars produced via core-collapse supernovae by asteroseismology, we perform a linear perturbation analysis using data obtained from two-dimensional numerical simulations. We employ 12 and 20 solar-mass progenitors a… ▽ More Gravitational wave signals from core-collapse supernovae are one of the important observables for extracting the information of dense matter. To extract the properties of proto-neutron stars produced via core-collapse supernovae by asteroseismology, we perform a linear perturbation analysis using data obtained from two-dimensional numerical simulations. We employ 12 and 20 solar-mass progenitors and compare two different treatments of gravity. One is a general relativistic one with a conformal flatness condition and the other is an effective gravitational potential mimicking the Tolman-Oppenheimer-Volkoff solution. We discuss how the frequencies of the proto-neutron star oscillations corresponding to the gravitational wave signals in the simulations depend on the proto-neutron star properties. In our models, we find that the gravitational wave frequencies of the proto-neutron stars determined with the Cowling approximation can be expressed to very good approximation as a function of the proto-neutron star average density almost independently of the progenitor mass, treatment of gravity in the simulations, and the interpolations in the simulations. On the other hand, if one considers the gravitational wave frequencies as a function of the surface gravity of proto-neutron stars, such a relation appears sensitive to the treatment of gravity and other numerical details in the simulations. Thus, the average density of proto-neutron stars seems more suitable for universally expressing the supernova gravitational wave frequencies, instead of the surface gravity. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: Accepted for publication in RPD

Report number: RIKEN-iTHEMS-Report-24

Convection and the Core $g$-mode in Proto-Compact Stars -- A detailed analysis

Authors: Pia Jakobus, Bernhard Mueller, Alexander Heger

Abstract: We present a detailed analysis of the dynamics of proto-compact star (PCS) convection and the core ${}^2\!g_1$-mode in core-collapse supernovae based on general relativistic 2D and 3D neutrino hydrodynamics simulations. Based on 2D simulations, we derive a mode relation for the core $g$-mode frequency in terms of PCS and equation of state parameters, and discuss its limits of accuracy. This relati… ▽ More We present a detailed analysis of the dynamics of proto-compact star (PCS) convection and the core ${}^2\!g_1$-mode in core-collapse supernovae based on general relativistic 2D and 3D neutrino hydrodynamics simulations. Based on 2D simulations, we derive a mode relation for the core $g$-mode frequency in terms of PCS and equation of state parameters, and discuss its limits of accuracy. This relation may prove useful for parameter inference from future supernova gravitational wave (GW) signals if the core $g$-mode or an emission gap at the avoided crossing with the fundamental mode can be detected. The current 3D simulation does not show GW emission from the core $g$-mode due to less power in high-frequency convective motions to excite the mode, however. Analysing the dynamics of PCS convection in 3D, we find that simple scaling laws for convective velocity from mixing-length theory (MLT) do not apply. Energy and lepton number transport is instead governed by a more complex balance between neutrino fluxes and turbulent fluxes that results in roughly uniform rates of change of entropy and lepton number in the PCS convection zone. Electron fraction and enthalpy contrasts in PCS convection are not well captured by the MLT gradient approximation. We find distinctly different spectra for the turbulent kinetic energy and turbulent fluctuations in the electron fraction, which scale approximately as $l^{-1}$ without a downturn at low $l$. We suggest that the different turbulence spectrum of the electron fraction is naturally expected for a passive scalar quantity. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Supernova Simulations

Authors: Bernhard Müller

Abstract: Magnetohydrodynamic simulations of core-collapse supernovae have become increasingly mature and important in recent years. Magnetic fields take center stage in scenarios for explaining hypernova explosions, but are now also considered in supernova theory more broadly as an important factor even in neutrino-driven explosions, especially in the context of neutron star birth properties. Here we prese… ▽ More Magnetohydrodynamic simulations of core-collapse supernovae have become increasingly mature and important in recent years. Magnetic fields take center stage in scenarios for explaining hypernova explosions, but are now also considered in supernova theory more broadly as an important factor even in neutrino-driven explosions, especially in the context of neutron star birth properties. Here we present an overview of simulation approaches currently used for magnetohydrodynamic supernova simulations and sketch essential physical concepts for understanding the role of magnetic fields in supernovae of slowly or rapidly rotating massive stars. We review progress on simulations of neutrino-driven supernovae, magnetorotational supernovae, and the relevant field amplification processes. Recent results on the nucleosynthesis and gravitational wave emission from magnetorotational supernovae are also discussed. We highlight efforts to provide better initial conditions for magnetohydrodynamic supernova models by simulating short phases of the progenitor evolution in 3D to address uncertainties in the treatment of rotation and magnetic fields in current stellar evolution models. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: 37 pages, 4 figures. Invited chapter for the edited book "New Frontiers in GRMHD Simulations" (Eds. C. Bambi, Y. Mizuno, S. Shashank and F. Yuan, Springer Singapore, expected in 2024)

Nucleosynthesis in the Innermost Ejecta of Magnetorotational Supernova Explosions in 3-dimensions

Authors: Shuai Zha, Bernhard Müller, Jade Powell

Abstract: Core-collapse supernova (CCSN) explosions powered by rotation and magnetic fields present an interesting astrophysical site for nucleosynthesis that potentially contributes to the production of $r$-process elements. Here we present yields of the innermost ejecta in 3D magnetorotational CCSN models simulated using the CoCoNuT-FMT code. Strong magnetic fields tap the rotational energy of the proto-n… ▽ More Core-collapse supernova (CCSN) explosions powered by rotation and magnetic fields present an interesting astrophysical site for nucleosynthesis that potentially contributes to the production of $r$-process elements. Here we present yields of the innermost ejecta in 3D magnetorotational CCSN models simulated using the CoCoNuT-FMT code. Strong magnetic fields tap the rotational energy of the proto-neutron star and lead to earlier and more energetic ($\sim 3\times 10^{51}$ erg) explosions than typical neutrino-driven CCSNe. Compared to a reference non-magnetic model, the ejecta in the magnetorotational models have much more neutron-rich components with Ye down to $\sim$0.25. Our post-processing calculations with the reaction network SkyNet show significant production of weak $r$-process elements up to mass number $\sim$130. We find negligible differences in the synthesis of heavy elements between two magnetorotational models with different initial field strength of 10$^{10}$ and 10$^{12}$ G, in accord with their similar explosion dynamics. The magnetorotational models produce about $\sim$0.19 and 0.14 Msun of radioactive $^{56}$Ni, on the low end of inferred hypernova nickel masses. The yields are publicly available at Zenodo: doi:10.5281/zenodo.10578981 for comparison with stellar abundance patterns, inclusion in modelling galactic chemical evolution, and comparison with other yield calculations. Our results add to the yet restricted corpus of nucleosynthesis yields from 3D magnetorotational supernova simulations and will help quantify yield uncertainties. △ Less

Submitted 11 May, 2024; v1 submitted 4 March, 2024; originally announced March 2024.

Comments: 19 pages, 9 figures, accepted by ApJ,yields available at doi.org/10.5281/zenodo.10578981

Determining the core-collapse supernova explosion mechanism with current and future gravitational-wave observatories

Authors: Jade Powell, Alberto Iess, Miquel Llorens-Monteagudo, Martin Obergaulinger, Bernhard Müller, Alejandro Torres-Forné, Elena Cuoco, José A. Font

Abstract: Gravitational waves are emitted from deep within a core-collapse supernova, which may enable us to determine the mechanism of the explosion from a gravitational-wave detection. Previous studies suggested that it is possible to determine if the explosion mechanism is neutrino-driven or magneto-rotationally powered from the gravitational-wave signal. However, long duration magneto-rotational wavefor… ▽ More Gravitational waves are emitted from deep within a core-collapse supernova, which may enable us to determine the mechanism of the explosion from a gravitational-wave detection. Previous studies suggested that it is possible to determine if the explosion mechanism is neutrino-driven or magneto-rotationally powered from the gravitational-wave signal. However, long duration magneto-rotational waveforms, that cover the full explosion phase, were not available during the time of previous studies, and explosions were just assumed to be magneto-rotationally driven if the model was rapidly rotating. Therefore, we perform an updated study using new 3D long-duration magneto-rotational core-collapse supernova waveforms that cover the full explosion phase, injected into noise for the Advanced LIGO, Einstein Telescope and NEMO gravitational-wave detectors. We also include a category for failed explosions in our signal classification results. We then determine the explosion mechanism of the signals using three different methods: Bayesian model selection, dictionary learning, and convolutional neural networks. The three different methods are able to distinguish between neutrino-driven explosions and magneto-rotational explosions, even if the neutrino-driven explosion model is rapidly rotating. However they can only distinguish between the non-exploding and neutrino-driven explosions for signals with a high signal to noise ratio. △ Less

Submitted 28 February, 2024; v1 submitted 29 November, 2023; originally announced November 2023.

Comments: Accepted in Physical Review D

Fallback onto Kicked Neutron Stars and its Effect on Spin-Kick Alignment

Authors: B. Müller

Abstract: Fallback in core-collapse supernova explosions is potentially of significant importance for the birth spins of neutron stars and black holes. It has recently been pointed out that the angular momentum imparted onto a compact remnant by fallback material is subtly intertwined with its kick because fallback onto a moving neutron star or black hole will preferentially come for a conical region around… ▽ More Fallback in core-collapse supernova explosions is potentially of significant importance for the birth spins of neutron stars and black holes. It has recently been pointed out that the angular momentum imparted onto a compact remnant by fallback material is subtly intertwined with its kick because fallback onto a moving neutron star or black hole will preferentially come for a conical region around its direction of travel. We show that contrary to earlier expectations such one-sided fallback accretion onto a neutron star will tend to produce spin-kick misalignment. Since the baroclinic driving term in the vorticity equation is perpendicular to the nearly radial pressure gradient, convective eddies in the progenitor as well as Rayleigh-Taylor plumes growing during the explosion primarily carry angular momentum perpendicular to the radial direction. Fallback material from the accretion volume of a moving neutron star therefore carries substantial angular momentum perpendicular to the kick velocity. We estimate the seed angular momentum fluctuations from convective motions in core-collapse supernova progenitors and argue that accreted fallback material will almost invariably be accreted with the maximum permissible specific angular momentum for reaching the Alfvén radius. This imposes a limit of $\mathord{\sim}10^{-2}M_\odot$ of fallback accretion for fast-spinning young neutron stars with periods of $\mathord{\sim}20\,\mathrm{ms}$ and less for longer birth spin periods. △ Less

Submitted 26 September, 2023; v1 submitted 16 August, 2023; originally announced August 2023.

Comments: 9 pages, 4 figures, accepted for publication in MNRAS. Fixed typo and added two references

3D Simulations of Magnetoconvection in a Rapidly Rotating Supernova Progenitor

Authors: Vishnu Varma, Bernhard Mueller

Abstract: We present a first 3D magnetohydrodynamic (MHD) simulation of oxygen, neon and carbon shell burning in a rapidly rotating 16 M_sun core-collapse supernova progenitor. We also run a purely hydrodynamic simulation for comparison. After 180s (15 and 7 convective turnovers respectively), the magnetic fields in the oxygen and neon shells achieve saturation at 10^{11}G and 5 x 10^{10}G. The strong Maxwe… ▽ More We present a first 3D magnetohydrodynamic (MHD) simulation of oxygen, neon and carbon shell burning in a rapidly rotating 16 M_sun core-collapse supernova progenitor. We also run a purely hydrodynamic simulation for comparison. After 180s (15 and 7 convective turnovers respectively), the magnetic fields in the oxygen and neon shells achieve saturation at 10^{11}G and 5 x 10^{10}G. The strong Maxwell stresses become comparable to the radial Reynolds stresses and eventually suppress convection. The suppression of mixing by convection and shear instabilities results in the depletion of fuel at the base of the burning regions, so that the burning shell eventually move outward to cooler regions, thus reducing the energy generation rate. The strong magnetic fields efficiently transport angular momentum outwards, quickly spinning down the rapidly rotating convective oxygen and neon shells and forcing them into rigid rotation. The hydrodynamic model shows complicated redistribution of angular momentum and develops regions of retrograde rotation at the base of the convective shells. We discuss implications of our results for stellar evolution and for the subsequent core-collapse supernova. The rapid redistribution of angular momentum in the MHD model casts some doubt on the possibility of retaining significant core angular momentum for explosions driven by millisecond magnetars. However, findings from multi-D models remain tentative until stellar evolution calculations can provide more consistent rotation profiles and estimates of magnetic field strengths to initialise multi-D simulations without substantial numerical transients. We also stress the need for longer simulations, resolution studies, and an investigation of non-ideal effects. △ Less

Submitted 23 October, 2023; v1 submitted 10 July, 2023; originally announced July 2023.

Comments: Submitted to MNRAS (14 pages, 11 Figures)

Synthetic Light Curves and Spectra from a Self-Consistent 2D Simulation of an Ultra-strippped Supernova

Authors: Thomas Maunder, Bernhard Müller, Fionntan Callan, Stuart Sim, Alexander Heger

Abstract: Spectroscopy is an important tool for providing insights into the structure of core-collapse supernova explosions. We use the Monte Carlo radiative transfer code ARTIS to compute synthetic spectra and light curves based on a two-dimensional explosion model of an ultra-stripped supernova. These calculations are designed both to identify observable fingerprints of ultra-stripped supernovae and as a… ▽ More Spectroscopy is an important tool for providing insights into the structure of core-collapse supernova explosions. We use the Monte Carlo radiative transfer code ARTIS to compute synthetic spectra and light curves based on a two-dimensional explosion model of an ultra-stripped supernova. These calculations are designed both to identify observable fingerprints of ultra-stripped supernovae and as a proof-of-principle for using synthetic spectroscopy to constrain the nature of stripped-envelope supernovae more broadly. We predict very characteristic spectral and photometric features for our ultra-stripped explosion model, but find that these do not match observed ultra-stripped supernova candidates like SN 2005ek. With a peak bolometric luminosity of $6.8\times10^{41}\,\mathrm{erg}\,\mathrm{s}^{-1}$, a peak magnitude of $-15.9\,\mathrm{mag}$ in R-band, and $Δm_{15,\mathrm{R}}=3.50$, the model is even fainter and evolves even faster than SN 2005ek as the closest possible analogue in photometric properties. The predicted spectra are extremely unusual. The most prominent features are Mg II lines at 2,800 Angstrom and 4,500 Angstrom and the infrared Ca triplet at late times. The Mg lines are sensitive to the multi-dimensional structure of the model and are viewing-angle dependent. They disappear due to line blanketing by Fe group elements in a spherically averaged model with additional microscopic mixing. In future studies, multi-D radiative transfer calculations need to be applied to a broader range of models to elucidate the nature of observed Type Ib/c supernovae. △ Less

Submitted 27 May, 2023; originally announced May 2023.

Comments: 13 pages, 14 figures, submitted to MNRAS

Black holes as the end state of stellar evolution: Theory and simulations

Authors: Alexander Heger, Bernhard Müller, Ilya Mandel

Abstract: The collapse of massive stars is one of the most-studied paths to black hole formation. In this chapter, we review black hole formation during the collapse of massive stars in the broader context of single and binary stellar evolution and the theory of supernova explosions. We provide a concise overview of the evolutionary channels that may lead to black hole formation -- the classical route of ir… ▽ More The collapse of massive stars is one of the most-studied paths to black hole formation. In this chapter, we review black hole formation during the collapse of massive stars in the broader context of single and binary stellar evolution and the theory of supernova explosions. We provide a concise overview of the evolutionary channels that may lead to black hole formation -- the classical route of iron core collapse, collapse due to pair instability in very massive stars, and the hypothetical scenario of supermassive star collapse. We then review the current understanding of the parameter space for black hole formation and black hole birth properties that has emerged from theoretical and computational modelling of supernova explosions and transient observations. Finally, we discuss what the intricate interplay between stellar evolution, stellar explosions, and binary interactions implies for the formation of stellar-mass black holes. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: 53 pages, 9 figures. This chapter is the pre-print of the version currently in production. Please cite this chapter as the following: A.Heger, B. Müller, and I. Mandel. "Black holes as the end state of stellar evolution: Theory and simulations," in The Encyclopedia of Cosmology (Set 2): Black Holes, edited by Z. Haiman (World Scientific, New Jersey, 2023)

Simulations of the progenitors of black hole-neutron star gravitational wave sources

Authors: Long Jiang, Wen-Cong Chen, Thomas M. Tauris, Bernhard Muller, Xiang-Dong Li

Abstract: Recent discoveries of gravitational wave (GW) events most likely originating from black hole (BH) + neutron star (NS) mergers reveal the existence of BH+NS binaries. The formation of BH+NS binaries and their merger rates through isolated binary evolution have been investigated extensively with population synthesis simulations. A detailed stellar evolution modelings of the formation of this populat… ▽ More Recent discoveries of gravitational wave (GW) events most likely originating from black hole (BH) + neutron star (NS) mergers reveal the existence of BH+NS binaries. The formation of BH+NS binaries and their merger rates through isolated binary evolution have been investigated extensively with population synthesis simulations. A detailed stellar evolution modelings of the formation of this population, however, is missing in the literature. In this work, we perform the first complete 1D model of more than 30 BH+NS progenitor systems which are calculated self-consistently until the iron core collapse with infall velocity exceeds 1000 km s^-1. Focusing on the progenitors of BH- NS GW sources, we apply the MESA code starting from a post-common envelope binary with short orbital period (< 1 day) consisting of a BH and a zero-age main-sequence helium star that experiences stable mass transfer. These NS masses could be significantly larger depending on the exact mass cut during the supernova explosion. These BH+NS systems are likely to merge and produce GW events within a Hubble time. System C is a potential progenitor of a GW200115-like event, while Systems A and B are possible candidates for a GW200105-like event and may represent the final destiny of the X-ray binary SS433. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 14pages, 9 figures, Accepted to be published on ApJ

Gravitational Waves from a Core g-Mode in Supernovae as Probes of the High-Density Equation of State

Authors: Pia Jakobus, Bernhard Müller, Alexander Heger, Shuai Zha, Jade Powell, Anton Motornenko, Jan Steinheimer, Horst Stoecker

Abstract: Using relativistic supernova simulations of massive progenitor stars with a quark-hadron equation of state (EoS) and a purely hadronic EoS, we identify a distinctive feature in the gravitational-wave signal that originates from a buoyancy-driven mode (g-mode) below the proto-neutron star convection zone. The mode frequency lies in the range $200\lesssim f\lesssim 800\,\text{Hz}$ and decreases with… ▽ More Using relativistic supernova simulations of massive progenitor stars with a quark-hadron equation of state (EoS) and a purely hadronic EoS, we identify a distinctive feature in the gravitational-wave signal that originates from a buoyancy-driven mode (g-mode) below the proto-neutron star convection zone. The mode frequency lies in the range $200\lesssim f\lesssim 800\,\text{Hz}$ and decreases with time. As the mode lives in the core of the proto-neutron star, its frequency and power are highly sensitive to the EoS, in particular the sound speed around twice saturation density. △ Less

Submitted 30 September, 2023; v1 submitted 16 January, 2023; originally announced January 2023.

Comments: Accepted Phys. Rev. Lett. (29 September 2023)

Light Curves of Type IIP Supernovae from Neutrino-driven Explosions of Red Supergiants Obtained by a Semi-analytic Approach

Authors: Shuai Zha, Bernhard Müller, Amy Weir, Alexander Heger

Abstract: Type IIP supernovae (SNe IIP) mark the explosive death of red supergiants (RSGs), evolved massive stars with an extended hydrogen envelope. They are the most common supernova type and allow for benchmarking of supernova explosion models by statistical comparison to observed population properties rather than comparing individual models and events. We construct a large synthetic set of SNe IIP light… ▽ More Type IIP supernovae (SNe IIP) mark the explosive death of red supergiants (RSGs), evolved massive stars with an extended hydrogen envelope. They are the most common supernova type and allow for benchmarking of supernova explosion models by statistical comparison to observed population properties rather than comparing individual models and events. We construct a large synthetic set of SNe IIP light curves (LCs) using the radiation hydrodynamics code \texttt{SNEC} and explosion energies and nickel masses obtained from an efficient semi-analytic model for two different sets of stellar progenitor models. By direct comparison we demonstrate that the semi-analytic model yields very similar predictions as alternative phenomenological explosion models based on one-dimensional simulations. We find systematic differences of a factor of $\mathord{\sim}2$ in plateau luminosities between the two progenitor sets due to different stellar radii, which highlights the importance of the RSG envelope structure as a major uncertainty in interpreting LCs of SNe IIP. A comparison to a volume-limited sample of observed SNe IIP shows decent agreement in plateau luminosity, plateau duration and nickel mass for at least one of the synthetic LC sets. The models, however, do not produce sufficient events with very small nickel mass $M_\mathrm{Ni}<0.01\,M_\odot$ and predict an anticorrelation between plateau luminosity and plateau duration that is not present in the observed sample, a result that warrants further study. Our results suggest that a better understanding of RSG stellar structure is no less important for reliably explaining the light curves of SNe IIP than the explosion physics. △ Less

Submitted 25 May, 2023; v1 submitted 1 January, 2023; originally announced January 2023.

Comments: 21 pages, 10 figures, 3 tables, accepted by ApJ, light curve parameter available at Zenodo:10.5281/zenodo.7354733

Three dimensional magnetorotational core-collapse supernova explosions of a 39 solar mass progenitor star

Authors: Jade Powell, Bernhard Mueller, David R. Aguilera-Dena, Norbert Langer

Abstract: We perform three-dimensional simulations of magnetorotational supernovae using a $39\,M_{\odot}$ progenitor star with two different initial magnetic field strengths of $10^{10}$ G and $10^{12}$ G in the core. Both models rapidly undergo shock revival and their explosion energies asymptote within a few hundred milliseconds to values of $\gtrsim 2\times10^{51}$ erg after conservatively correcting fo… ▽ More We perform three-dimensional simulations of magnetorotational supernovae using a $39\,M_{\odot}$ progenitor star with two different initial magnetic field strengths of $10^{10}$ G and $10^{12}$ G in the core. Both models rapidly undergo shock revival and their explosion energies asymptote within a few hundred milliseconds to values of $\gtrsim 2\times10^{51}$ erg after conservatively correcting for the binding energy of the envelope. Magnetically collimated, non-relativistic jets form in both models, though the jets are subject to non-axisymmetric instabilities. The jets do not appear crucial for driving the explosion, as they only emerge once the shock has already expanded considerably. Our simulations predict moderate neutron star kicks of about $150\, \mathrm{km}\,\mathrm{s}^{-1}$, no spin-kick alignment, and rapid early spin-down that would result in birth periods of about $20\, \mathrm{ms}$, too slow to power an energetic gamma-ray burst jet. More than $0.2\,M_\odot$ of iron-group material are ejected, but we estimate that the mass of ejected $^{56}\mathrm{Ni}$ will be considerably smaller as the bulk of this material is neutron-rich. Explosive burning does not contribute appreciable amounts of $^{56}\mathrm{Ni}$ because the burned material originates from the slightly neutron-rich silicon shell. The iron-group ejecta also show no pronounced bipolar geometry by the end of the simulations. The models thus do not immediately fit the characteristics of observed hypernovae, but may be representative of other transients with moderately high explosion energies. The gravitational-wave emission reaches high frequencies of up to 2000 Hz and amplitudes of over 100 cm. The gravitational-wave emission is detectable out to distances of $\sim4$ Mpc in the planned Cosmic Explorer detector. △ Less

Submitted 4 May, 2023; v1 submitted 30 November, 2022; originally announced December 2022.

A Black-Hole Excision Scheme for General Relativistic Core-Collapse Supernova Simulations

Authors: B. Sykes, B. Mueller, I. Cordero-Carrión, P. Cerdá-Durán, J. Novak

Abstract: Fallback supernovae and the collapsar scenario for long-gamma ray burst and hypernovae have received considerable interest as pathways to black-hole formation and extreme transient events. Consistent simulations of these scenarios require a general relativistic treatment and need to deal appropriately with the formation of a singularity. Free evolution schemes for the Einstein equations can handle… ▽ More Fallback supernovae and the collapsar scenario for long-gamma ray burst and hypernovae have received considerable interest as pathways to black-hole formation and extreme transient events. Consistent simulations of these scenarios require a general relativistic treatment and need to deal appropriately with the formation of a singularity. Free evolution schemes for the Einstein equations can handle the formation of black holes by means of excision or puncture schemes. However, in constrained schemes, which offer distinct advantages in long-term numerical stability in stellar collapse simulations over well above $10^{4}$ light-crossing time scales, the dynamical treatment of black-hole spacetimes is more challenging. Building on previous work on excision in conformally flat spacetimes, we here present the implementation of a black-hole excision scheme for supernova simulations with the CoCoNuT-FMT neutrino transport code. We describe in detail a choice of boundary conditions that ensures long-time numerical stability, and also address upgrades to the hydrodynamics solver that are required to stably evolve the relativistic accretion flow onto the black hole. The scheme is currently limited to a spherically symmetric metric, but the hydrodynamics can be treated multi-dimensionally. For demonstration, we present a spherically symmetric simulation of black-hole formation in an $85 M_\odot$ star, as well as a two-dimensional simulation of the fallback explosion of the same progenitor. These extend past 9s and 0.3s after black-hole formation, respectively. △ Less

Submitted 22 May, 2024; v1 submitted 23 October, 2022; originally announced October 2022.

Comments: 19 pages, 12 figures. Accepted by PRD. Updated to correct typographical error in funding acknowledgements

Calibration of neutron star natal kick velocities to isolated pulsar observations

Authors: Veome Kapil, Ilya Mandel, Emanuele Berti, Bernhard Müller

Abstract: Current prescriptions for supernova natal kicks in rapid binary population synthesis simulations are based on fits of simple functions to single pulsar velocity data. We explore a new parameterization of natal kicks received by neutron stars in isolated and binary systems developed by Mandel & Müller, which is based on 1D and 3D supernova simulations and accounts for the physical correlations betw… ▽ More Current prescriptions for supernova natal kicks in rapid binary population synthesis simulations are based on fits of simple functions to single pulsar velocity data. We explore a new parameterization of natal kicks received by neutron stars in isolated and binary systems developed by Mandel & Müller, which is based on 1D and 3D supernova simulations and accounts for the physical correlations between progenitor properties, remnant mass, and the kick velocity. We constrain two free parameters in this model using very long baseline interferometry velocity measurements of Galactic single pulsars. We find that the inferred values of natal kick parameters do not differ significantly between single and binary evolution scenarios. The best-fit values of these parameters are $v_{\rm ns} = 520$ km s$^{-1}$ for the scaling pre-factor for neutron star kicks, and $σ_{\rm ns}=0.3$ for the fractional stochastic scatter in the kick velocities. △ Less

Submitted 28 February, 2023; v1 submitted 19 September, 2022; originally announced September 2022.

Comments: 8 pages, 8 figures, 1 table. Updated to match MNRAS version

Journal ref: MNRAS Volume 519, Issue 4, March 2023, Pages 5893-5901

Laboratory spectroscopy of theoretical ices: Predictions for JWST and test for astrochemical models

Authors: B. Müller, B. M. Giuliano, A. Vasyunin, G. Fedoseev, P. Caselli

Abstract: Context. The gas and ice-grain chemistry of the pre-stellar core L1544 has been the subject of several observations and modelling studies conducted in the past years. The chemical composition of the ice mantles reflects the environmental physical changes along the temporal evolution. The investigation outcome hints at a layered structure of interstellar ices with mainly H$_2$O in the inner layers… ▽ More Context. The gas and ice-grain chemistry of the pre-stellar core L1544 has been the subject of several observations and modelling studies conducted in the past years. The chemical composition of the ice mantles reflects the environmental physical changes along the temporal evolution. The investigation outcome hints at a layered structure of interstellar ices with mainly H$_2$O in the inner layers and an increasing amount of CO near the surface. The morphology of interstellar ice analogues can be investigated experimentally. Aims. This research presents a new approach of a three-dimensional fit where observational results are first fitted with a gas-grain chemical model. Then, based on the numerical results the laboratory IR spectra are recorded for interstellar ice analogues in a layered and in a mixed morphology. These results can then be compared with future James Webb Space Telescope (JWST) observations. Special attention is paid to the inclusion of the IR inactive species N$_2$ and O$_2$. Methods. Ice analogue spectra containing the most abundant predicted molecules were recorded at a temperature of 10 K using a Fourier transform infrared spectrometer. In the case of layered ice we deposited a H$_2$O-CO-N$_2$-O$_2$ mixture on top of a H2O-CH$_3$OH-N$_2$ ice, while in the case of mixed ice we examined a H$_2$O-CH$_3$OH-N$_2$-CO composition. Results. Following the changing composition and structure of the ice, we find differences in the absorption bands for most of the examined vibrational modes. The extent of observed changes in the IR band profiles will allow us to analyse the structure of ice mantles in L1544 from future observations by the JWST. Conclusions. The comparison of our spectroscopic measurements with upcoming JWST observations is crucial in order to put stringent constraints on the chemical and physical structure of dust icy mantles, and to explain surface chemistry. △ Less

Submitted 5 September, 2022; v1 submitted 16 August, 2022; originally announced August 2022.

Comments: Accepted for publication in A&A, 13 pages, 5 figures

Journal ref: A&A 668, A46 (2022)

Broadband spectroscopy of astrophysical ice analogues: II. Optical constants of CO and CO$_2$ ices in the terahertz and infrared ranges

Authors: A. A. Gavdush, F. Kruczkiewicz, B. M. Giuliano, B. Müller, G. A. Komandin, T. Grassi, P. Theulé, K. I. Zaytsev, A. V. Ivlev, P. Caselli

Abstract: Context: Broadband optical constants of astrophysical ice analogues in the infrared (IR) and terahertz (THz) ranges are required for modeling the dust continuum emission and radiative transfer in dense and cold regions, where thick icy mantles are formed on the surface of dust grains. Aims: In this paper, the THz time-domain spectroscopy (TDS) and the Fourier-transform IR spectroscopy (FTIR) are c… ▽ More Context: Broadband optical constants of astrophysical ice analogues in the infrared (IR) and terahertz (THz) ranges are required for modeling the dust continuum emission and radiative transfer in dense and cold regions, where thick icy mantles are formed on the surface of dust grains. Aims: In this paper, the THz time-domain spectroscopy (TDS) and the Fourier-transform IR spectroscopy (FTIR) are combined to study optical constants of CO and CO$_2$ ices in the broad THz-IR spectral range. Methods: The measured ices are grown at cryogenic temperatures by gas deposition on a cold Si window. A method to quantify the broadband THz-IR optical constants of ices is developed based on the direct reconstruction of the complex refractive index of ices in the THz range from the TDS data, and the use of the Kramers-Kronig relation in the IR range for the reconstruction from the FTIR data. Uncertainties of the Kramers-Kronig relation are eliminated by merging the THz and IR spectra. The reconstructed THz-IR response is then analyzed using classical models of complex dielectric permittivity. Results: The complex refractive index of CO and CO$_2$ ices deposited at the temperature of $28$ K is obtained in the range of 0.3--12.0 THz. Based on the measured dielectric constants, opacities of the astrophysical dust with CO and CO$_2$ icy mantles are computed. Conclusions: The developed method can be used for a model-independent reconstruction of optical constants of various astrophysical ice analogs in a broad THz-IR range. Such data can provide important benchmarks to interpret the broadband observations from the existing and future ground-based facilities and space telescopes. The reported results will be useful to model sources that show a drastic molecular freeze-out, such as central regions of prestellar cores and mid-planes of protoplanetary disks, as well as CO and CO$_2$ snow lines in disks. △ Less

Submitted 22 July, 2022; originally announced July 2022.

Comments: Accepted for publication in A&A, 9 pages, 7 figures

Journal ref: A&A 667, A49 (2022)

3D Simulations of Strongly Magnetised Non-Rotating Supernovae: Explosion Dynamics and Remnant Properties

Authors: Vishnu Varma, Bernhard Mueller, Fabian R. N. Schneider

Abstract: We investigate the impact of strong initial magnetic fields in core-collapse supernovae of non-rotating progenitors by simulating the collapse and explosion of a 16.9 Msun star for a strong- and weak-field case assuming a twisted-torus field with initial central field strengths of ~1012 G and ~106 G. The strong-field model has been set up with a view to the fossil-field scenario for magnetar forma… ▽ More We investigate the impact of strong initial magnetic fields in core-collapse supernovae of non-rotating progenitors by simulating the collapse and explosion of a 16.9 Msun star for a strong- and weak-field case assuming a twisted-torus field with initial central field strengths of ~1012 G and ~106 G. The strong-field model has been set up with a view to the fossil-field scenario for magnetar formation and emulates a pre-collapse field configuration that may occur in massive stars formed by a merger. This model undergoes shock revival already 100 s after bounce and reaches an explosion energy of 9.3 x 10^50 erg at 310 ms, in contrast to a more delayed and less energetic explosion in the weak-field model. The strong magnetic fields help trigger a neutrino-driven explosion early on, which results in a rapid rise and saturation of the explosion energy. Dynamically, the strong initial field leads to a fast build-up of magnetic fields in the gain region to 40% of kinetic equipartition and also creates sizable pre-shock ram pressure perturbations that are known to be conducive to asymmetric shock expansion. For the strong-field model, we find an extrapolated neutron star kick of ~350 km/s, a spin period of ~70 ms, and no spin-kick alignment. The dipole field strength of the proto-neutron star is 2 x 10^14 G by the end of the simulation with a declining trend. Surprisingly, the surface dipole field in the weak-field model is stronger, which argues against a straightforward connection between pre-collapse fields and the birth magnetic fields of neutron stars. △ Less

Submitted 23 November, 2022; v1 submitted 23 April, 2022; originally announced April 2022.

Comments: 16 pages, 20 figures Accepted for publication in MNRAS

The Role of the Hadron-Quark Phase Transition in Core-Collapse Supernovae

Authors: Pia Jakobus, Bernhard Mueller, Alexander Heger, Anton Motornenko, Jan Steinheimer, Horst Stoecker

Abstract: The hadron-quark phase transition in quantum chromodyanmics has been suggested as an alternative explosion mechanism for core-collapse supernovae. We study the impact of three different hadron-quark equations of state (EoS) with first-order (DD2F\_SF, STOS-B145) and second-order (CMF) phase transitions on supernova dynamics by performing 97 simulations for solar- and zero-metallicity progenitors i… ▽ More The hadron-quark phase transition in quantum chromodyanmics has been suggested as an alternative explosion mechanism for core-collapse supernovae. We study the impact of three different hadron-quark equations of state (EoS) with first-order (DD2F\_SF, STOS-B145) and second-order (CMF) phase transitions on supernova dynamics by performing 97 simulations for solar- and zero-metallicity progenitors in the range of $14\texttt{-}100\,\text{M}_\odot$. We find explosions only for two low-compactness models ($14 \text{M}_\odot$ and $16\,\text{M}_\odot$) with the DD2F\_SF EoS, both with low explosion energies of $\mathord{\sim}10^{50}\,\mathrm{erg}$. These weak explosions are characterised by a neutrino signal with several mini-bursts in the explosion phase due to complex reverse shock dynamics, in addition to the typical second neutrino burst for phase-transition driven explosions. The nucleosynthesis shows significant overproduction of nuclei such as $^{90}\mathrm{Zr}$ for the $14\,\text{M}_\odot$ zero-metallicity model and $^{94}\mathrm{Zr}$ for the $16\,\text{M}_\odot$ solar-metallicity model, but the overproduction factors are not large enough to place constraints on the occurrence of such explosions. Several other low-compactness models using the DD2F\_SF EoS and two high-compactness models using the STOS EoS end up as failed explosions and emit a second neutrino burst. For the CMF EoS, the phase transition never leads to a second bounce and explosion. For all three EoS, inverted convection occurs deep in the core of the proto-compact star due to anomalous behaviour of thermodynamic derivatives in the mixed phase, which heats the core to entropies up to $4k_\text{B}/\text{baryon}$ and may have a distinctive gravitational wave signature, also for a second-order phase transition. △ Less

Submitted 22 August, 2022; v1 submitted 21 April, 2022; originally announced April 2022.

Stripped-envelope stars in different metallicity environments. II. Type I supernovae and compact remnants

Authors: David R. Aguilera-Dena, Bernhard Müller, John Antoniadis, Norbert Langer, Luc Dessart, Alejandro Vigna-Gómez, Sung-Chul Yoon

Abstract: Stripped-envelope stars can be observed as Wolf-Rayet (WR) stars, or as less luminous hydrogen-poor stars with low mass loss rates and transparent winds. Both types are potential progenitors of Type I core-collapse supernovae (SNe). We use grids of core-collapse models obtained from helium stars at different metallicities to study the effects of metallicity on the transients and remnants these sta… ▽ More Stripped-envelope stars can be observed as Wolf-Rayet (WR) stars, or as less luminous hydrogen-poor stars with low mass loss rates and transparent winds. Both types are potential progenitors of Type I core-collapse supernovae (SNe). We use grids of core-collapse models obtained from helium stars at different metallicities to study the effects of metallicity on the transients and remnants these stars produce. We characterise the surface and core properties of our core collapse models, and investigate their explodability employing three criteria. In cases where explosions are predicted, we estimate the ejecta mass, explosion energy, nickel mass and neutron star (NS) mass. Otherwise, we predict the mass of the resulting black hole (BH). We construct a simplified population model, and find that the properties SNe and compact objects depend strongly on metallicity. Ejecta masses and explosion energies for Type Ic SNe are best reproduced by models with Z=0.04 which exhibit strong winds during core helium burning. This implies that either their mass loss rates are underestimated, or that Type Ic SN progenitors experience mass loss through other mechanisms before exploding. The distributions of ejecta masses, explosion energies and nickel mass for Type Ib SNe are not well reproduced by progenitor models with WR mass loss, but are better reproduced if we assume no mass loss in progenitors with luminosities below the minimum WR star luminosity. We find that Type Ic SNe become more common as metallicity increases, and that the vast majority of progenitors of Type Ib SNe must be transparent-wind stripped-envelope stars. We find several models with pre-collapse CO-masses of up to $\sim 30 M_{\odot}$ may form $\sim 3 M_{\odot}$ BHs in fallback SNe. This may carry important consequences for our understanding of SNe, binary BH and NS systems, X-ray binary systems and gravitational wave transients. △ Less

Submitted 31 March, 2022; originally announced April 2022.

Comments: 30 pages, 20 figures, submitted to Astronomy & Astrophysics

Journal ref: A&A 671, A134 (2023)

Inferring Astrophysical Parameters of Core-Collapse Supernovae from their Gravitational-Wave Emission

Authors: Jade Powell, Bernhard Müller

Abstract: Nearby core-collapse supernovae (CCSNe) are powerful multi-messenger sources for gravitational-wave, neutrino and electromagnetic telescopes as they emit gravitational waves in the ideal frequency band for ground based detectors. Once a CCSN gravitational-wave signal is detected, we will need to determine the parameters of the signal, and understand how those parameters relate to the source's expl… ▽ More Nearby core-collapse supernovae (CCSNe) are powerful multi-messenger sources for gravitational-wave, neutrino and electromagnetic telescopes as they emit gravitational waves in the ideal frequency band for ground based detectors. Once a CCSN gravitational-wave signal is detected, we will need to determine the parameters of the signal, and understand how those parameters relate to the source's explosion, progenitor and remnant properties. This is a challenge due to the stochastic nature of CCSN explosions, which is imprinted on their time series gravitational waveforms. In this paper, we perform Bayesian parameter estimation of CCSN signals using an asymmetric chirplet signal model to represent the dominant high-frequency mode observed in spectrograms of CCSN gravitational-wave signals. We use design sensitivity Advanced LIGO noise and CCSN waveforms from four different hydrodynamical supernova simulations with a range of different progenitor stars. We determine how well our model can reconstruct time-frequency images of the emission modes, and show how well we can determine parameters of the signal such as the frequency, amplitude, and duration. We show how the parameters of our signal model may allow us to place constraints on the proto-neutron star mass and radius, the turbulent kinetic energy onto the proto-neutron star, and the time of shock revival. △ Less

Submitted 2 March, 2022; v1 submitted 4 January, 2022; originally announced January 2022.

Radial Distribution of the Dust Comae of Comets 45P/Honda-Mrkos-Pajduusáková and 46P/Wirtanen

Authors: C. Lejoly, W. Harris, N. Samarasinha, B. E. A. Mueller, E. Howell, J. Bodnarik, A. Springmann, T. Kareta, B. Sharkey, J. Noonan, L. R. Bedin, J. -G. Bosch, A. Brosio, E. Bryssinck, J. -B. de Vanssay, F. -J. Hambsch, O. Ivanova, V. Krushinsky, Z. -Y. Lin, F. Manzini, A. Maury, N. Moriya, P. Ochner, V. Oldani

Abstract: There was an unprecedented opportunity to study the inner dust coma environment, where the dust and gas are not entirely decoupled, of comets 45P/Honda-Mrkos-Pajduusáková (45P/HMP) from Dec. 26, 2016 - Mar. 15, 2017, and 46P/Wirtanen from Nov. 10, 2018 - Feb. 13, 2019, both in visible wavelengths. The radial profile slopes of these comets were measured in the R and HB-BC filters most representativ… ▽ More There was an unprecedented opportunity to study the inner dust coma environment, where the dust and gas are not entirely decoupled, of comets 45P/Honda-Mrkos-Pajduusáková (45P/HMP) from Dec. 26, 2016 - Mar. 15, 2017, and 46P/Wirtanen from Nov. 10, 2018 - Feb. 13, 2019, both in visible wavelengths. The radial profile slopes of these comets were measured in the R and HB-BC filters most representative of dust, and deviations from a radially expanding coma were identified as significant. The azimuthally averaged radial profile slope of comet 45P/HMP gradually changes from -1.81 $\pm$ 0.20 at 5.24 days pre-perihelion to -0.35 $\pm$ 0.16 at 74.41 days post perihelion. Contrastingly, the radial profile slope of 46P/Wirtanen stays fairly constant over the observed time period at -1.05 $\pm$ 0.05. Additionally, we find that the radial profile of 46P/Wirtanen is azimuthally dependent on the skyplane-projected solar position angle, while that of 45P/HMP is not. These results suggest that comet 45P/HMP and 46P/Wirtanen have vastly different coma dust environments and that their dust properties are distinct. As evident from these two comets, well-resolved inner comae are vital for detailed characterization of dust environments. △ Less

Submitted 20 December, 2021; originally announced December 2021.

Comments: 21 pages, 13 figures, to be published in the Planetary Science Journal

The Great Slump: Mrk926 reveals discrete and varying Balmer line satellite components during a drastic phase of decline

Authors: Wolfram Kollatschny, Martin W. Ochmann, Shai Kaspi, Claas Schumacher, Ehud Behar, Doron Chelouche, Keith Horne, Bjoern Mueller, Stephen E. Rafter, Rolf Chini, Martin Haas, Malte A. Probst

Abstract: Mrk926 is known to be a highly variable AGN with very broad line profiles. We studied the continuum and line profile variations of this object with high temporal resolution in order to determine its broad-line region structure and to derive its BH mass. We carried out a high-cadence spectroscopic variability campaign of Mrk926 with the 10m HET telescope, aided by photometric V-band data taken at t… ▽ More Mrk926 is known to be a highly variable AGN with very broad line profiles. We studied the continuum and line profile variations of this object with high temporal resolution in order to determine its broad-line region structure and to derive its BH mass. We carried out a high-cadence spectroscopic variability campaign of Mrk926 with the 10m HET telescope, aided by photometric V-band data taken at the Wise Observatory, over a period of about five months. We extracted spectroscopic continuum and line light curves, and computed CCFs as well as velocity-resolved CCFs with respect to the combined continuum V-band light curve. The continuum luminosity of Mrk926 showed a drastic decrease during our campaign: it dropped to less than 50% of its original value within only 2.5 months. The spectra show complex and very broad Balmer line profiles, including outer Balmer satellites ranging from +-5000 to +-13,000 kms-1. The Ha, Hb, and HeI5876 line light curves are delayed relative to the continuum light curve. The Ha and Hb lines show two velocity-delay structures in the central part of their line profile (within +- 5000 kms-1), at ~10 and ~57 ld and at ~5 and ~48 ld, respectively. These structures might be interpreted as the signature of a line-emitting ring, inclined by ~50 deg to the line of sight and orbiting the BH at radii, R, of 33.5 and 26.5 ld. We determined continuum luminosities, log(lambda\,L_(lambda)/erg s-1), of 43.68 to 44.13, which are in good agreement with the established R_(BLR)-L_(AGN) relation. We derive a black hole mass of 1.1 +- 0.2 * 10^8 M_solar; this indicates a low Eddington ratio, which decreased from 8 to 3 percent during our campaign. The Balmer satellite components show a response to the continuum variations on the order of only 3-5 days. We attribute this to the central line segment and the Balmer satellites having different, spatially distinct regions of origin. △ Less

Submitted 31 March, 2022; v1 submitted 19 October, 2021; originally announced October 2021.

Comments: 23 pages, 14 figures, published in Astronomy & Astrophysics

Journal ref: A&A 657, A122 (2022)

Explodability fluctuations of massive stellar cores enable asymmetric compact object mergers such as GW190814

Authors: John Antoniadis, David R. Aguilera-Dena, Alejandro Vigna-Gómez, Michael Kramer, Norbert Langer, Bernhard Müller, Thomas M. Tauris, Chen Wang, Xiao-Tian Xu

Abstract: The first three observing runs with Advanced LIGO and Virgo have resulted in the detection of binary black hole mergers (BBH) with highly unequal mass components, which are difficult to reconcile with standard formation paradigms. The most representative of these is GW190814, a highly asymmetric merger between a 23 M$_{\odot}$ black hole and a 2.6 M$_{\odot}$ compact object. Here, we explore recen… ▽ More The first three observing runs with Advanced LIGO and Virgo have resulted in the detection of binary black hole mergers (BBH) with highly unequal mass components, which are difficult to reconcile with standard formation paradigms. The most representative of these is GW190814, a highly asymmetric merger between a 23 M$_{\odot}$ black hole and a 2.6 M$_{\odot}$ compact object. Here, we explore recent results suggesting that a sizeable fraction of stars with pre-collapse carbon-oxygen core masses above 10 M$_{\odot}$, and extending up to at least 30 M$_{\odot}$, may produce objects inside the so-called lower mass gap that bridges the division between massive pulsars and BHs in Galactic X-ray binaries. We demonstrate that such an explosion landscape would naturally cause a fraction of massive binaries to produce GW190814-like systems instead of symmetric-mass BBHs. We present examples of specific evolutionary channels leading to the formation of GW190814 and GW200210, a 24+2.8 M$_{\odot}$ merger discovered during the O3b observing run. We estimate the merger-rate density of these events in our scenario to be $\mathcal{O}$(5%) of the total BBH merger rate. Finally, we discuss the broader implications of this formation channel for compact object populations, and its possible relevance to less asymmetric merger events such as GW200105 and GW200115 △ Less

Submitted 21 December, 2021; v1 submitted 4 October, 2021; originally announced October 2021.

Comments: A&A Letters, in press

Journal ref: A&A 657, L6 (2022)

A Comparison of 2D Magnetohydrodynamic Supernova Simulations with the CoCoNuT-FMT and Aenus-Alcar Codes

Authors: Vishnu Varma, Bernhard Mueller, Martin Obergaulinger

Abstract: Code comparisons are a valuable tool for the verification of supernova simulation codes and the quantification of model uncertainties. Here we present a first comparison of axisymmetric magnetohydrodynamic (MHD) supernova simulations with the CoCoNuT-FMT and Aenus-Alcar codes, which use distinct methods for treating the MHD induction equation and the neutrino transport. We run two sets of simulati… ▽ More Code comparisons are a valuable tool for the verification of supernova simulation codes and the quantification of model uncertainties. Here we present a first comparison of axisymmetric magnetohydrodynamic (MHD) supernova simulations with the CoCoNuT-FMT and Aenus-Alcar codes, which use distinct methods for treating the MHD induction equation and the neutrino transport. We run two sets of simulations of a rapidly rotating 35M gamma-ray burst progenitor model with different choices for the initial field strength, namely 10^12 G for the maximum poloidal and toroidal field in the strong-field case and 10^10 G in the weak-field case. We also investigate the influence of the Riemann solver and the resolution in CoCoNuT-FMT. The dynamics is qualitatively similar for both codes and robust with respect to these numerical details, with a rapid magnetorotational explosion in the strong-field case and a delayed neutrino-driven explosion in the weak-field case. Despite relatively similar shock trajectories, we find sizeable differences in many other global metrics of the dynamics, like the explosion energy and the magnetic energy of the proto-neutron star. Further differences emerge upon closer inspection, for example, the disk-like surface structure of the proto-neutron star proves highly sensitives to numerical details. The electron fraction distribution in the ejecta as a crucial determinant for the nucleosynthesis is qualitatively robust, but the extent of neutron- or proton-rich tails is sensitive to numerical details. Due to the complexity of the dynamics, the ultimate cause of model differences can rarely be uniquely identified, but our comparison helps gauge uncertainties inherent in current MHD supernova simulations. △ Less

Submitted 15 October, 2021; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: 16 pages, 20 figures, accepted for publication in MNRAS

Repeating gas ejection events from comet 45P/Honda-Mrkos-Pajdušáková

Authors: Alessondra Springmann, Walter M. Harris, Erin L. Ryan, Cassandra Lejoly, Ellen S. Howell, Beatrice E. A. Mueller, Nalin H. Samarasinha, Laura M. Woodney, Jordan K. Steckloff

Abstract: Studying materials released from Jupiter-family comets (JFCs) -- as seen in their inner comæ, the envelope of gas and dust that forms as the comet approaches the Sun -- improves the understanding of their origin and evolutionary history. As part of a coordinated, multi-wavelength observing campaign, we observed comet 45P/Honda-Mrkos-Pajdušáková during its close approach to Earth in February 2017.… ▽ More Studying materials released from Jupiter-family comets (JFCs) -- as seen in their inner comæ, the envelope of gas and dust that forms as the comet approaches the Sun -- improves the understanding of their origin and evolutionary history. As part of a coordinated, multi-wavelength observing campaign, we observed comet 45P/Honda-Mrkos-Pajdušáková during its close approach to Earth in February 2017. Narrowband observations were taken using the Bok 90" telescope at Kitt Peak National Observatory on February 16 and 17 UT, revealing gas and dust structures. We observed different jet directions for different volatile species, implying source region heterogeneity, consistent with other ground-based and \textit{in situ} observations of other comet nuclei. A repeating feature visible in CN and C$_2$ images on February 16 was also observed on February 17 with an interval of $7.6\pm0.1$ hours, consistent with the rotation period of the comet derived from Arecibo Observatory radar observations. The repeating feature's projected gas velocity away from the nucleus is 0.8 km s$^{-1}$, with an expansion velocity of 0.5 km s$^{-1}$. A bright compact spot adjacent to the nucleus provides a lower limit of the amount of material released in one cycle of $\sim$9.2 kg, depending on composition -- a quantity small enough to be produced by repeated exposure of nucleus ices to sunlight. This repeating CN jet, forming within 400 km of the nucleus, may be typical of inner coma behavior in JFCs; however, similar features could be obscured by other processes and daughter product species when viewed from distances further than the scale length of CN molecules. △ Less

Submitted 28 November, 2021; v1 submitted 21 July, 2021; originally announced July 2021.

Comments: 11 pages, 1 table, 7 figures. Accepted at The Planetary Science Journal on November 28, 2021

Spectroscopic measurements of CH$_3$OH in layered and mixed interstellar ice analogues

Authors: B. Müller, B. M. Giuliano, M. Goto, P. Caselli

Abstract: Context. The molecular composition of interstellar ice mantles is defined by gas-grain processes in molecular clouds, with the main components being $H_2O$, $CO$, and $CO_2$. $CH_3OH$ ice is detected towards the denser regions, where large amounts of $CO$ freeze out and get hydrogenated. Heating from nearby protostars can further change the ice structure and composition. Despite the several observ… ▽ More Context. The molecular composition of interstellar ice mantles is defined by gas-grain processes in molecular clouds, with the main components being $H_2O$, $CO$, and $CO_2$. $CH_3OH$ ice is detected towards the denser regions, where large amounts of $CO$ freeze out and get hydrogenated. Heating from nearby protostars can further change the ice structure and composition. Despite the several observations of icy features towards molecular clouds and along the line of site of protostars, it is not yet clear if interstellar ices are mixed or if they have a layered structure. Aims. We aim to examine the effect of mixed and layered ice growth in ice mantle analogues, with focus on the position and shape of methanol infrared bands, so future observations could shed light on the structure of interstellar ices in different environments. Methods. Mixed and layered ice samples were deposited on a cold substrate kept at T = 10 K using a closed-cycle cryostat placed in a vacuum chamber. The spectroscopic features were analysed by FTIR spectroscopy. Different proportions of the most abundant four molecules in ice mantles, namely $H_2O$, $CO$, $CO_2$, and $CH_3OH$, were investigated, with special attention on the analysis of the $CH_3OH$ bands. Results. We measure changes in the position and shape of the CH and CO stretching bands of $CH_3OH$ depending on the mixed or layered nature of the ice sample. Spectroscopic features of methanol are also found to change due to heating. Conclusions. A layered ice structure best reproduces the $CH_3OH$ band position recently observed towards a pre-stellar core and in star-forming regions. Based on our experimental results, we conclude that observations of $CH_3OH$ ices can provide information about the structure of interstellar ices, and we expect JWST to put stringent constraints on the layered or mixed nature of ices in different interstellar environments. △ Less

Submitted 20 July, 2021; originally announced July 2021.

Comments: 11 pages, 8 figures, 3 tables

Journal ref: A&A 652, A126 (2021)

Differential Rotation in a 3D Simulation of Oxygen Shell Burning

Authors: Lucy O. McNeill, Bernhard Müller

Abstract: We study differential rotation in late-stage shell convection in a 3D hydrodynamic simulation of a rapidly rotating $16M_\odot$ helium star with a particular focus on the convective oxygen shell. We find that the oxygen shell develops a quasi-stationary pattern of differential rotation that is neither described by uniform angular velocity as assumed in current stellar evolution models of supernova… ▽ More We study differential rotation in late-stage shell convection in a 3D hydrodynamic simulation of a rapidly rotating $16M_\odot$ helium star with a particular focus on the convective oxygen shell. We find that the oxygen shell develops a quasi-stationary pattern of differential rotation that is neither described by uniform angular velocity as assumed in current stellar evolution models of supernova progenitors, nor by uniform specific angular momentum. Instead, the oxygen shell develops a positive angular velocity gradient with faster rotation at the equator than at the pole by tens of percent. We show that the angular momentum transport inside the convection zone is not adequately captured by a diffusive mixing-length flux proportional to the angular velocity or angular momentum gradient. Zonal flow averages reveal stable large-scale meridional flow and an entropy deficit near the equator that mirrors the patterns in the angular velocity. The structure of the flow is reminiscent of simulations of stellar surface convection zones and the differential rotation of the Sun, suggesting that similar effects are involved; future simulations will need to address in more detail how the interplay of buoyancy, inertial forces, and turbulent stresses shapes differential rotation during late-stage convection in massive stars. If convective regions develop positive angular velocity gradients, angular momentum could be shuffled out of the core region more efficiently, potentially making the formation of millisecond magnetars more difficult. Our findings have implications for neutron star birth spin periods and supernova explosion scenarios that involve rapid core rotation. △ Less

Submitted 9 October, 2022; v1 submitted 30 June, 2021; originally announced July 2021.

Comments: 13 pages, 6 figures. Published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 509, Issue 1, January 2022, pp.818-830

Current-Conserving Relativistic Linear Response for Collisional Plasmas

Authors: Martin Formanek, Christopher Grayson, Johann Rafelski, Berndt Müller

Abstract: We investigate the response of a relativistic plasma to electromagnetic fields in the framework of the Boltzmann equation incorporating a collision term in the relaxation rate approximation selected in a form assuring current conservation. We obtain an explicit solution for the linearized perturbation of the Fermi-Dirac equilibrium distribution in terms of the average relaxation rate $κ$. We study… ▽ More We investigate the response of a relativistic plasma to electromagnetic fields in the framework of the Boltzmann equation incorporating a collision term in the relaxation rate approximation selected in a form assuring current conservation. We obtain an explicit solution for the linearized perturbation of the Fermi-Dirac equilibrium distribution in terms of the average relaxation rate $κ$. We study the resulting covariant, gauge invariant, and current conserving form of the polarization tensor in the ultrarelativistic and non-relativistic limits. We evaluate the susceptibility in the ultrarelativistic limit and explore their dependence on $κ$. Finally, we study the dispersion relations for the longitudinal and transverse poles of the propagator. We show that for $κ> 2ω_p$, where $ω_p$ is the plasma frequency, the plasma wave modes are overdamped. In the opposite case, $κ\ll ω_p$, the propagating plasma modes are weakly damped. △ Less

Submitted 9 September, 2021; v1 submitted 14 May, 2021; originally announced May 2021.

Comments: 14 pages, 9 figures

Gravitational Wave Signals from Two-Dimensional Core-Collapse Supernova Models with Rotation and Magnetic Fields

Authors: Rylan Jardine, Jade Powell, Bernhard Müller

Abstract: We investigate the impact of rotation and magnetic fields on the dynamics and gravitational wave emission in 2D core-collapse supernova simulations with neutrino transport. We simulate 17 different models of $15\,M_\odot$ and $39\,M_\odot$ progenitor stars with various initial rotation profiles and initial magnetic fields strengths up to $10^{12}\, \mathrm{G}$, assuming a dipolar field geometry in… ▽ More We investigate the impact of rotation and magnetic fields on the dynamics and gravitational wave emission in 2D core-collapse supernova simulations with neutrino transport. We simulate 17 different models of $15\,M_\odot$ and $39\,M_\odot$ progenitor stars with various initial rotation profiles and initial magnetic fields strengths up to $10^{12}\, \mathrm{G}$, assuming a dipolar field geometry in the progenitor. Strong magnetic fields generally prove conducive to shock revival, though this trend is not without exceptions. The impact of rotation on the post-bounce dynamics is more variegated, in line with previous studies. A significant impact on the time-frequency structure of the gravitational wave signal is found only for rapid rotation or strong initial fields. For rapid rotation, the angular momentum gradient at the proto-neutron star surface can appreciably affect the frequency of the dominant mode, so that known analytic relations for the high-frequency emission band no longer hold. In case of two magnetorotational explosion models, the deviation from these analytic relations is even more pronounced. One of the magnetorotational explosions has been evolved to more than half a second after the onset of the explosion and shows a subsidence of high-frequency emission at late times. Its most conspicuous gravitational wave signature is a high-amplitude tail signal. We also estimate the maximum detection distances for our waveforms. The magnetorotational models do not stick out for higher detectability during the post-bounce and explosion phase. △ Less

Submitted 22 December, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 19 pages, 12 figures, substantially revised version. Model m15afb12 has been recalculated to remove numerical contamination of high-frequency tail of GW spectrum. One new model m15nrb11 and additional analysis added. Accepted for publiation in MNRAS

The final core collapse of pulsational pair instability supernovae

Authors: Jade Powell, Bernhard Müller, Alexander Heger

Abstract: We present 3D core-collapse supernova simulations of massive Pop-III progenitor stars at the transition to the pulsational pair instability regime. We simulate two progenitor models with initial masses of $85\,\mathrm{M}_{\odot}$ and $100\,\mathrm{M}_\odot$ with the LS220, SFHo, and SFHx equations of state. The $85\,\mathrm{M}_{\odot}$ progenitor experiences a pair instability pulse coincident wit… ▽ More We present 3D core-collapse supernova simulations of massive Pop-III progenitor stars at the transition to the pulsational pair instability regime. We simulate two progenitor models with initial masses of $85\,\mathrm{M}_{\odot}$ and $100\,\mathrm{M}_\odot$ with the LS220, SFHo, and SFHx equations of state. The $85\,\mathrm{M}_{\odot}$ progenitor experiences a pair instability pulse coincident with core collapse, whereas the $100\,\mathrm{M}_{\odot}$ progenitor has already gone through a sequence of four pulses $1\mathord,500$ years before collapse in which it ejected its H and He envelope. The $85\,\mathrm{M}_{\odot}$ models experience shock revival and then delayed collapse to a black hole (BH) due to ongoing accretion within hundreds of milliseconds. The diagnostic energy of the incipient explosion reaches up to $2.7\times10^{51}\,\mathrm{erg}$ in the SFHx model. Due to the high binding energy of the metal core, BH collapse by fallback is eventually unavoidable, but partial mass ejection may be possible. The $100\,\mathrm{M}_\odot$ models have not achieved shock revival or undergone BH collapse by the end of the simulation. All models exhibit relatively strong gravitational-wave emission both in the high-frequency g-mode emission band and at low frequencies. The SFHx and SFHo models show clear emission from the standing accretion shock instability. For our models, we estimate maximum detection distances of up to $\mathord{\sim}46\,\mathrm{kpc}$ with LIGO and $\mathord{\sim} 850\,\mathrm{kpc}$ with Cosmic Explorer. △ Less

Submitted 26 February, 2021; v1 submitted 18 January, 2021; originally announced January 2021.

3D Simulations of Oxygen Shell Burning with and without Magnetic Fields

Authors: Vishnu Varma, Bernhard Müller

Abstract: We present a first 3D magnetohydrodynamic (MHD) simulation of convective oxygen and neon shell burning in a non-rotating $18\, M_\odot$ star shortly before core collapse to study the generation of magnetic fields in supernova progenitors. We also run a purely hydrodynamic control simulation to gauge the impact of the magnetic fields on the convective flow and on convective boundary mixing. After a… ▽ More We present a first 3D magnetohydrodynamic (MHD) simulation of convective oxygen and neon shell burning in a non-rotating $18\, M_\odot$ star shortly before core collapse to study the generation of magnetic fields in supernova progenitors. We also run a purely hydrodynamic control simulation to gauge the impact of the magnetic fields on the convective flow and on convective boundary mixing. After about 17 convective turnover times, the magnetic field is approaching saturation levels in the oxygen shell with an average field strength of $\mathord{\sim}10^{10}\, \mathrm{G}$, and does not reach kinetic equipartition. The field remains dominated by small to medium scales, and the dipole field strength at the base of the oxygen shell is only $10^{9}\, \mathrm{G}$. The angle-averaged diagonal components of the Maxwell stress tensor mirror those of the Reynolds stress tensor, but are about one order of magnitude smaller. The shear flow at the oxygen-neon shell interface creates relatively strong fields parallel to the convective boundary, which noticeably inhibit the turbulent entrainment of neon into the oxygen shell. The reduced ingestion of neon lowers the nuclear energy generation rate in the oxygen shell and thereby slightly slows down the convective flow. Aside from this indirect effect, we find that magnetic fields do not appreciably alter the flow inside the oxygen shell. We discuss the implications of our results for the subsequent core-collapse supernova and stress the need for longer simulations, resolution studies, and an investigation of non-ideal effects for a better understanding of magnetic fields in supernova progenitors. △ Less

Submitted 25 March, 2021; v1 submitted 1 January, 2021; originally announced January 2021.

Comments: 12 pages, 9 figures, revised version accepted for publication in MNRAS

Self-consistent 3D Supernova Models From -7 Minutes to +7 Seconds: a 1-bethe Explosion of a ~19 Solar-mass Progenitor

Authors: R. Bollig, N. Yadav, D. Kresse, H. -Th. Janka, B. Mueller, A. Heger

Abstract: To date, modern three-dimensional (3D) supernova (SN) simulations have not demonstrated that explosion energies of 10^{51} erg (=1 bethe = 1B) or more are possible for neutrino-driven SNe of non/slow-rotating M < 20 solar-mass progenitors. We present the first such model, considering a non-rotating, solar-metallicity 18.88 solar-mass progenitor, whose final 7 minutes of convective oxygen-shell bur… ▽ More To date, modern three-dimensional (3D) supernova (SN) simulations have not demonstrated that explosion energies of 10^{51} erg (=1 bethe = 1B) or more are possible for neutrino-driven SNe of non/slow-rotating M < 20 solar-mass progenitors. We present the first such model, considering a non-rotating, solar-metallicity 18.88 solar-mass progenitor, whose final 7 minutes of convective oxygen-shell burning were simulated in 3D and showed a violent oxygen-neon shell merger prior to collapse. A large set of 3D SN-models was computed with the Prometheus-Vertex code, whose improved convergence of the two-moment equations with Boltzmann closure allows now to fully exploit the implicit neutrino-transport treatment. Nuclear burning is treated with a 23-species network. We vary the angular grid resolution and consider different nuclear equations of state and muon formation in the proto-neutron star (PNS), which requires six-species transport with coupling of all neutrino flavors across all energy-momentum groups. Elaborate neutrino transport was applied until ~2 seconds after bounce. In one case the simulation was continued to >7 seconds with an approximate treatment of neutrino effects that allows for seamless continuation without transients. A spherically symmetric neutrino-driven wind does not develop. Instead, accretion downflows to the PNS and outflows of neutrino-heated matter establish a monotonic rise of the explosion energy until ~7 seconds post bounce, when the outgoing shock reaches about 50,000 km and enters the He-layer. The converged value of the explosion energy at infinity (with overburden subtracted) is roughly 1B and the ejected 56Ni mass up to 0.087 solar masses, both within a few 10 percent of the SN 1987A values. The final NS mass and kick are about 1.65 solar masses and over 450 km/s, respectively. △ Less

Submitted 15 April, 2021; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 29 pages, 17 figures, 1 table, 3 appendices; considerably extended due to comments by referee and readers; accepted by ApJ

The chemical signature of jet-driven hypernovae

Authors: J. J. Grimmett, Bernhard Müller, Alexander Heger, Projjwal Banerjee, Martin Obergaulinger

Abstract: Hypernovae powered by magnetic jets launched from the surface of rapidly rotating millisecond magnetars are one of the leading models to explain broad-lined Type Ic supernovae (SNe Ic-BL), and have been implicated as an important source of metal enrichment in the early Universe. We investigate the nucleosynthesis in such jet-driven hypernovae using a parameterised, but physically motivated, approa… ▽ More Hypernovae powered by magnetic jets launched from the surface of rapidly rotating millisecond magnetars are one of the leading models to explain broad-lined Type Ic supernovae (SNe Ic-BL), and have been implicated as an important source of metal enrichment in the early Universe. We investigate the nucleosynthesis in such jet-driven hypernovae using a parameterised, but physically motivated, approach that analytically relates an artificially injected jet energy flux to the power available from the energy in differential rotation in the proto-neutron star. We find ejected $^{56}\mathrm{Ni}$ masses of $0.05\,\mathrm{M}_\odot - 0.45\,\mathrm{M}_\odot$ in our most energetic models with explosion energy $>10^{52}\,\mathrm{erg}$. This is in good agreement with the range of observationally inferred values for SNe Ic-BL. The $^{56}\mathrm{Ni}$ is mostly synthesised in the shocked stellar envelope, and is therefore only moderately sensitive to the jet composition. Jets with a high electron fraction $Y_\mathrm{e}=0.5$ eject more $^{56}\mathrm{Ni}$ by a factor of 2 than neutron-rich jets. We can obtain chemical abundance profiles in good agreement with the average chemical signature observed in extremely metal-poor (EMP) stars presumably polluted by hypernova ejecta. Notably, $\mathrm{[Zn/Fe]} \gtrsim 0.5$ is consistently produced in our models. For neutron-rich jets, there is a significant r-process component, and agreement with EMP star abundances in fact requires either a limited contribution from neutron-rich jets or a stronger dilution of r-process material in the interstellar medium than for the slow SN ejecta outside the jet. The high $\mathrm{[C/Fe]}\gtrsim 0.7$ observed in many EMP stars cannot be consistently achieved due to the large mass of iron in the ejecta, however, and remains a challenge for jet-driven hypernovae based on the magneto-rotational mechanism. △ Less

Submitted 6 December, 2020; v1 submitted 13 October, 2020; originally announced October 2020.

Nucleosynthesis of an $11.8\,M_\odot$ Supernova with 3D Simulation of the Inner Ejecta: Overall Yields and Implications for Short-Lived Radionuclides in the Early Solar System

Authors: Andre Sieverding, Bernhard Mueller, Yong-Zhong Qian

Abstract: Based on a 3D supernova simulation of an $11.8\,M_\odot$ progenitor model with initial solar composition, we study the nucleosynthesis using tracers covering the innermost $0.1\,M_\odot$ of the ejecta. These ejecta are mostly proton-rich and contribute significant amounts of $^{45}$Sc and $^{64}$Zn. The production of heavier isotopes is sensitive to the electron fraction and hence the neutrino emi… ▽ More Based on a 3D supernova simulation of an $11.8\,M_\odot$ progenitor model with initial solar composition, we study the nucleosynthesis using tracers covering the innermost $0.1\,M_\odot$ of the ejecta. These ejecta are mostly proton-rich and contribute significant amounts of $^{45}$Sc and $^{64}$Zn. The production of heavier isotopes is sensitive to the electron fraction and hence the neutrino emission from the proto-neutron star. The yields of these isotopes are rather uncertain due to the approximate neutrino transport used in the simulation. In order to obtain the total yields for the whole supernova, we combine the results from the tracers with those for the outer layers from a suitable 1D model. Using the yields of short-lived radionuclides (SLRs), we explore the possibility that an $11.8\,M_\odot$ supernova might have triggered the formation of the solar system and provided some of the SLRs measured in meteorites. In particular, we discuss two new scenarios that can account for at least the data on $^{41}$Ca, $^{53}$Mn, and $^{60}$Fe without exceeding those on the other SLRs. △ Less

Submitted 17 December, 2020; v1 submitted 28 August, 2020; originally announced August 2020.

Comments: 16 pages, 9 Figures, 4 Tables

Journal ref: The Astrophysical Journal, Volume 904, Issue 2, id.163, 14 pp. (2020)

Pre-collapse Properties of Superluminous Supernovae and Long Gamma-Ray Burst Progenitor Models

Authors: David R. Aguilera-Dena, Norbert Langer, John Antoniadis, Bernhard Müller

Abstract: We analyze the properties of 42 rapidly rotating, low metallicity, quasi-chemically homogeneously evolving stellar models in the mass range between 4 and 45 $\,\mathrm{M}_\odot$ at the time of core collapse. Such models were proposed as progenitors for both superluminous supernovae (SLSNe) and long duration gamma-ray bursts (lGRBs), and the Type Ic-BL supernovae (SNe) that are associated with them… ▽ More We analyze the properties of 42 rapidly rotating, low metallicity, quasi-chemically homogeneously evolving stellar models in the mass range between 4 and 45 $\,\mathrm{M}_\odot$ at the time of core collapse. Such models were proposed as progenitors for both superluminous supernovae (SLSNe) and long duration gamma-ray bursts (lGRBs), and the Type Ic-BL supernovae (SNe) that are associated with them. Our findings suggest that whether these models produce a magnetar driven SLSN explosion or a near-critically rotating black hole (BH) is not a monotonic function of the initial mass. Rather, their explodability varies non-monotonically depending on the late core evolution, once chemical homogeneity is broken. Using different explodability criteria we find that our models have a clear preference to produce SLSNe at lower masses, and lGRBs at higher masses; but find several exceptions, expecting lGRBs to form from stars as low as 10 $\,\mathrm{M}_\odot$, and SLSNe with progenitors as massive as 30 $\,\mathrm{M}_\odot$. In general, our models reproduce the predicted angular momenta, ejecta masses and magnetic field strengths at core collapse inferred for SLSNe and lGRBs, and suggest significant interaction with their circumstellar medium, particularly for explosions with low ejecta mass. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Comments: Accepted for publication in ApJ. 17 pages, 11 figures, 2 tables

Pre-supernova evolution, compact object masses and explosion properties of stripped binary stars

Authors: F. R. N. Schneider, Ph. Podsiadlowski, B. Müller

Abstract: Most massive stars are born in binary or higher-order multiple systems and exchange mass with a companion during their lives. In particular, the progenitors of a large fraction of compact object mergers, and Galactic neutron stars (NSs) and black holes (BHs) have been stripped off their envelopes by a binary companion. Here, we study the evolution of single and stripped binary stars up to core col… ▽ More Most massive stars are born in binary or higher-order multiple systems and exchange mass with a companion during their lives. In particular, the progenitors of a large fraction of compact object mergers, and Galactic neutron stars (NSs) and black holes (BHs) have been stripped off their envelopes by a binary companion. Here, we study the evolution of single and stripped binary stars up to core collapse with the stellar evolution code MESA and their final fates with a parametric supernova (SN) model. We find that stripped binary stars can have systematically different pre-SN structures compared to genuine single stars and thus also different SN outcomes. The bases of these differences are already established by the end of core helium burning and are preserved up to core collapse. We find a non-monotonic pattern of NS and BH formation as a function of CO core mass that is different in single and stripped binary stars. In terms of initial masses, single stars of >35 Msun all form BHs, while this transition is only at 70 Msun in stripped stars. On average, stripped stars give rise to lower NS and BH masses, higher explosion energies, higher kick velocities and higher nickel yields. Within a simplified population synthesis model, we show that our results lead to a significant reduction of the rates of BH-NS and BH-BH mergers with respect to typical assumptions made on NS and BH formation. Therefore, we predict lower detection rates of such merger events by, e.g., advanced LIGO than is often considered. We further show how features in the NS-BH mass distribution of single and stripped stars relate to the chirp-mass distribution of compact object mergers. Further implications of our findings are discussed with respect to the missing red-supergiant problem, a possible mass gap between NSs and BHs, X-ray binaries and observationally inferred nickel masses from Type Ib/c and IIP Sne. [abridged] △ Less

Submitted 22 October, 2020; v1 submitted 19 August, 2020; originally announced August 2020.

Comments: 25 pages (incl. appendix), 17 figures, 2 tables; final version accepted for publication in A&A

Journal ref: A&A 645, A5 (2021)

A 3D Simulation of a Neutrino-Driven Supernova Explosion Aided By Convection and Magnetic Fields

Authors: Bernhard Müller, Vishnu Varma

Abstract: We study the impact of a small-scale dynamo in core-collapse supernovae using a 3D neutrino magnetohydrodynamics simulation of a $15 M_\odot$ progenitor. The weak seed field is amplified exponentially in the gain region once neutrino-driven convection develops, and remains dominated by small-scale structures. About $250\, \mathrm{ms}$ after bounce, the field energy in the gain region reaches… ▽ More We study the impact of a small-scale dynamo in core-collapse supernovae using a 3D neutrino magnetohydrodynamics simulation of a $15 M_\odot$ progenitor. The weak seed field is amplified exponentially in the gain region once neutrino-driven convection develops, and remains dominated by small-scale structures. About $250\, \mathrm{ms}$ after bounce, the field energy in the gain region reaches $\mathord{\sim} 50\%$ of kinetic equipartition. This supports the development of a neutrino-driven explosion with modest global anisotropy, which does not occur in a corresponding model without magnetic fields. Our results suggest that magnetic fields may play a beneficial subsidiary role in neutrino-driven supernovae even without rapid progenitor rotation. Further investigation into the nature of magnetohydrodynamic turbulence in the supernova core is required. △ Less

Submitted 9 August, 2020; v1 submitted 9 July, 2020; originally announced July 2020.

Comments: 6 pages, 3 figures, accepted for publication in MNRAS Letters

Neutron Star Extreme Matter Observatory: A kilohertz-band gravitational-wave detector in the global network

Authors: K. Ackley, V. B. Adya, P. Agrawal, P. Altin, G. Ashton, M. Bailes, E. Baltinas, A. Barbuio, D. Beniwal, C. Blair, D. Blair, G. N. Bolingbroke, V. Bossilkov, S. Shachar Boublil, D. D. Brown, B. J. Burridge, J. Calderon Bustillo, J. Cameron, H. Tuong Cao, J. B. Carlin, S. Chang, P. Charlton, C. Chatterjee, D. Chattopadhyay, X. Chen , et al. (139 additional authors not shown)

Abstract: Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provid… ▽ More Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica. △ Less

Submitted 5 November, 2020; v1 submitted 6 July, 2020; originally announced July 2020.

Comments: Accepted for publication in PASA

Journal ref: PASA (2020) 37, e047

Simple recipes for compact remnant masses and natal kicks

Authors: Ilya Mandel, Bernhard Müller

Abstract: Based on recent results from three-dimensional supernova simulations and semi-analytical parametrised models, we develop analytical prescriptions for the dependence of the mass of neutron stars and black holes and the natal kicks, if any, on the pre-supernova carbon-oxygen core and helium shell masses. Our recipes are probabilistic rather than deterministic in order to account for the intrinsic st… ▽ More Based on recent results from three-dimensional supernova simulations and semi-analytical parametrised models, we develop analytical prescriptions for the dependence of the mass of neutron stars and black holes and the natal kicks, if any, on the pre-supernova carbon-oxygen core and helium shell masses. Our recipes are probabilistic rather than deterministic in order to account for the intrinsic stochasticity of stellar evolution and supernovae. We anticipate that these recipes will be particularly useful for rapid population synthesis, and we illustrate their application to distributions of remnant masses and kicks for a population of single stars. △ Less

Submitted 29 September, 2020; v1 submitted 15 June, 2020; originally announced June 2020.

Comments: Version accepted to MNRAS

Hydrodynamics of core-collapse supernovae and their progenitors

Authors: B. Müller

Abstract: Multi-dimensional fluid flow plays a paramount role in the explosions of massive stars as core-collapse supernovae. In recent years, three-dimensional (3D) simulations of these phenomena have matured significantly. Considerable progress has been made towards identifying the ingredients for shock revival by the neutrino-driven mechanism, and successful explosions have already been obtained in a num… ▽ More Multi-dimensional fluid flow plays a paramount role in the explosions of massive stars as core-collapse supernovae. In recent years, three-dimensional (3D) simulations of these phenomena have matured significantly. Considerable progress has been made towards identifying the ingredients for shock revival by the neutrino-driven mechanism, and successful explosions have already been obtained in a number of self-consistent 3D models. These advances also bring new challenges, however. Prompted by a need for increased physical realism and meaningful model validation, supernova theory is now moving towards a more integrated view that connects multi-dimensional phenomena in the late convective burning stages prior to collapse, the explosion engine, and mixing instabilities in the supernova envelope. Here we review our current understanding of multi-D fluid flow in core-collapse supernovae and their progenitors. We start by outlining specific challenges faced by hydrodynamic simulations of core-collapse supernovae and of the late convective burning stages. We then discuss recent advances and open questions in theory and simulations. △ Less

Submitted 9 June, 2020; originally announced June 2020.

Comments: Invited review article for Living Reviews in Computational Astrophysics. 100 pages, 15 figures

Journal ref: Living Rev Comput Astrophys 6, 3 (2020)

Stochastic Core Spin-Up in Massive Stars -- Implications of 3D Simulations of Oxygen Shell Burning

Authors: L. O. McNeill, B. Müller

Abstract: It has been suggested based on analytic theory that even in non-rotating supernova progenitors stochastic spin-up by internal gravity waves (IGWs) during the late burning stages can impart enough angular momentum to the core to result in neutron star birth spin periods below 100ms, and a relatively firm upper limit of 500ms for the spin period. We here investigate this process using a 3D simulatio… ▽ More It has been suggested based on analytic theory that even in non-rotating supernova progenitors stochastic spin-up by internal gravity waves (IGWs) during the late burning stages can impart enough angular momentum to the core to result in neutron star birth spin periods below 100ms, and a relatively firm upper limit of 500ms for the spin period. We here investigate this process using a 3D simulation of oxygen shell burning in a $3M_\odot$ He star. Our model indicates that stochastic spin-up by IGWs is less efficient than previously thought. We find that the stochastic angular momentum flux carried by waves excited at the shell boundary is significantly smaller for a given convective luminosity and turnover time than would be expected from simple dimensional analysis. This can be explained by noting that the waves launched by overshooting convective plumes contain modes of opposite angular wave number with similar amplitudes, so that the net angular momentum of excited wave packets almost cancels. We find that the wave-mediated angular momentum flux from the oxygen shell follows a random walk, but again dimensional analysis overestimates the random walk amplitudes since the correlation time is only a fraction of the convective turnover time. Extrapolating our findings over the entire life time of the last burning stages prior to collapse, we predict that the core angular momentum from stochastic spin-up would translate into long birth spin periods of several seconds for low-mass progenitors and no less than 100ms even for high-mass progenitors. △ Less

Submitted 30 July, 2020; v1 submitted 3 June, 2020; originally announced June 2020.

Comments: 12 pages, 6 figures. MNRAS accepted

Searching for Electromagnetic Counterparts to Gravitational-wave Merger Events with the Prototype Gravitational-wave Optical Transient Observer (GOTO-4)

Authors: B. P. Gompertz, R. Cutter, D. Steeghs, D. K. Galloway, J. Lyman, K. Ulaczyk, M. J. Dyer, K. Ackley, V. S. Dhillon, P. T. O'Brien, G. Ramsay, S. Poshyachinda, R. Kotak, L. Nuttall, R. P. Breton, E. Pallé, D. Pollacco, E. Thrane, S. Aukkaravittayapun, S. Awiphan, M. J. I. Brown, U. Burhanudin, P. Chote, A. A. Chrimes, E. Daw , et al. (27 additional authors not shown)

Abstract: We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light… ▽ More We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light curves of on- and off-axis gamma-ray bursts and kilonovae. In cases where the source region was observable immediately, GOTO-4 was able to respond to a GW alert in less than a minute. The average time of first observation was $8.79$ hours after receiving an alert ($9.90$ hours after trigger). A mean of $732.3$ square degrees were tiled per event, representing on average $45.3$ per cent of the LVC probability map, or $70.3$ per cent of the observable probability. This coverage will further improve as the facility scales up alongside the localisation performance of the evolving gravitational-wave detector network. Even in its 4-telescope prototype configuration, GOTO is capable of detecting AT2017gfo-like kilonovae beyond 200~Mpc in favourable observing conditions. We cannot currently place meaningful electromagnetic limits on the population of distant ($\hat{D}_L = 1.3$~Gpc) binary black hole mergers because our test models are too faint to recover at this distance. However, as GOTO is upgraded towards its full 32-telescope, 2 node (La Palma \& Australia) configuration, it is expected to be sufficiently sensitive to cover the predicted O4 binary neutron star merger volume, and will be able to respond to both northern and southern triggers. △ Less

Submitted 28 July, 2020; v1 submitted 31 March, 2020; originally announced April 2020.

Comments: 15 pages, 7 figures, 3 tables. Accepted for publication in MNRAS. Author's final submitted version

The impact of fallback on the compact remnants and chemical yields of core-collapse supernovae

Authors: Conrad Chan, Bernhard Mueller, Alexander Heger

Abstract: Fallback in core-collapse supernovae plays a crucial role in determining the properties of the compact remnants and of the ejecta composition. We perform three-dimensional simulations of mixing and fallback for selected non-rotating supernova models to study how explosion energy and asymmetries correlate with the remnant mass, remnant kick, and remnant spin. We find that the strongest kick and spi… ▽ More Fallback in core-collapse supernovae plays a crucial role in determining the properties of the compact remnants and of the ejecta composition. We perform three-dimensional simulations of mixing and fallback for selected non-rotating supernova models to study how explosion energy and asymmetries correlate with the remnant mass, remnant kick, and remnant spin. We find that the strongest kick and spin is imparted by partial fallback in an asymmetric explosion. Black hole (BH) kicks of several hundred $\mathrm{km}\,\mathrm{s}^{-1}$ and spin parameters of $\mathord{\sim}0.25$ can be obtained in this scenario. If the initial explosion energy barely exceeds the envelope binding energy, stronger fallback results, and the remnant kick and spin remain small. If the explosion energy is high with respect to the envelope binding energy, there is little fallback with a small effect on the remnant kick, but the spin-up by fallback can be substantial. For a non-rotating $12\,\mathrm{M}_\odot$ progenitor, we find that the neutron star (NS) is spun up to millisecond periods. The high specific angular momentum of the fallback material can also lead to disk formation around black holes. Fallback may thus be a pathway towards millisecond-magnetar or collapsar-type engines for hypernovae and gamma-ray bursts that does not require rapid progenitor rotation. Within our small set of simulations, none reproduced the peculiar layered fallback necessary to explain the metal-rich iron-poor composition of many carbon-enhanced metal-poor (CEMP) stars. Models with different explosion energy and different realisations of asymmetries may, however, be compatible with CEMP abundance patterns. △ Less

Submitted 24 May, 2020; v1 submitted 9 March, 2020; originally announced March 2020.

Comments: 12 pages, 13 figures, accepted for publication in MNRAS

A Novel Multi-Dimensional Boltzmann Neutrino Transport Scheme for Core-Collapse Supernovae

Authors: Conrad Chan, Bernhard Mueller

Abstract: We introduce a new discrete-ordinate scheme for solving the general relativistic Boltzmann transport equation in the context of core-collapse supernovae. Our algorithm avoids the need to spell out the complicated advection terms in energy and angle that arise when the transport equation is formulated in spherical polar coordinates, in the comoving frame, or in a general relativistic spacetime. We… ▽ More We introduce a new discrete-ordinate scheme for solving the general relativistic Boltzmann transport equation in the context of core-collapse supernovae. Our algorithm avoids the need to spell out the complicated advection terms in energy and angle that arise when the transport equation is formulated in spherical polar coordinates, in the comoving frame, or in a general relativistic spacetime. We instead approach the problem by calculating the advection of neutrinos across momentum space using an intuitive particle-like approach that has excellent conservation properties and fully accounts for Lorentz boosts, general relativistic effects, and grid geometry terms. In order to avoid the need for a global implicit solution, time integration is performed using a locally implicit Lax-Wendroff scheme that correctly reproduces the diffusion limit. This will facilitate the use of our method on massively parallel distributed-memory architectures. We have verified the accuracy and stability of our scheme with a suite of test problems in spherical symmetry and axisymmetry. To demonstrate that the new algorithm works stably in core-collapse supernova simulations, we have coupled it to the general relativistic hydrodynamics code CoCoNuT and present a first demonstration run of a $20M_\odot$ progenitor with a reduced set of neutrino opacities. △ Less

Submitted 26 June, 2020; v1 submitted 5 March, 2020; originally announced March 2020.

Comments: 22 pages, 15 figures, accepted for publication in MNRAS

Three-dimensional core-collapse supernova simulations of massive and rotating progenitors

Authors: Jade Powell, Bernhard Müller

Abstract: We present three-dimensional simulations of the core-collapse of massive rotating and non-rotating progenitors performed with the general relativistic neutrino hydrodynamics code CoCoNuT-FMT and analyse their explosion properties and gravitational-wave signals. The progenitor models include Wolf-Rayet stars with initial helium star masses of $39\,M_{\odot}$ and $20\,M_{\odot}$, and an… ▽ More We present three-dimensional simulations of the core-collapse of massive rotating and non-rotating progenitors performed with the general relativistic neutrino hydrodynamics code CoCoNuT-FMT and analyse their explosion properties and gravitational-wave signals. The progenitor models include Wolf-Rayet stars with initial helium star masses of $39\,M_{\odot}$ and $20\,M_{\odot}$, and an $18\,M_{\odot}$ red supergiant. The $39\,M_{\odot}$ model is a rapid rotator, whereas the two other progenitors are non-rotating. Both Wolf-Rayet models produce healthy neutrino-driven explosions, whereas the red supergiant model fails to explode. By the end of the simulations, the explosion energies have already reached $1.1\times 10^{51}\,\mathrm{erg}$ and $0.6\times 10^{51}\,\mathrm{erg}$ for the $39\,M_{\odot}$ and $20\,M_{\odot}$ model, respectively. The explosions produce neutron stars of relatively high mass, but with modest kicks. Due to the alignment of the bipolar explosion geometry with the rotation axis, there is a relatively small misalignment of $30^\circ$ between the spin and the kick in the $39\,M_{\odot}$ model. In terms of gravitational-wave signals, the massive and rapidly rotating $39\,M_{\odot}$ progenitor stands out by large gravitational-wave amplitudes that would make it detectable out to almost 2 Mpc by the Einstein Telescope. For this model, we find that rotation significantly changes the dependence of the characteristic gravitational-wave frequency of the f-mode on the proto-neutron star parameters compared to the non-rotating case. The other two progenitors have considerably smaller detection distances, despite significant low-frequency emission in the most sensitive frequency band of current gravitational-wave detectors due to the standing accretion shock instability in the $18\,M_{\odot}$ model. △ Less

Submitted 24 February, 2020; originally announced February 2020.

The Chemical Evolution of Iron-Peak Elements with Hypernovae

Authors: J. J. Grimmett, Amanda I. Karakas, Alexander Heger, Bernhard Müller, Christopher West

Abstract: We calculate the mean evolution of the iron-peak abundance ratios [(Cr,Mn,Co,Zn)/Fe] in the Galaxy, using modern supernova and hypernova chemical yields and a Galactic Chemical Evolution code that assumes homogeneous chemical evolution. We investigate a range of hypernova occurrence rates and are able to produce a chemical composition that is a reasonable fit to the observed values in metal-poor s… ▽ More We calculate the mean evolution of the iron-peak abundance ratios [(Cr,Mn,Co,Zn)/Fe] in the Galaxy, using modern supernova and hypernova chemical yields and a Galactic Chemical Evolution code that assumes homogeneous chemical evolution. We investigate a range of hypernova occurrence rates and are able to produce a chemical composition that is a reasonable fit to the observed values in metal-poor stars. This requires a hypernova occurence rate that is large (50%) in the early Universe, decreasing throughout evolution to a value that is within present day observational constraints (>~ 1%). A large hypernova occurence rate is beneficial to matching the high [Zn/Fe] observed in the most metal-poor stars, although including hypernovae with progenitor mass >= 60 solar masses is detrimental to matching the observed [(Mn,Co)/Fe] evolution at low [Fe/H]. A significant contribution from HNe seems to be critical for producing supersolar [(Co,Zn)/Fe] at low metallicity, though more work will need to be done in order to match the most extreme values. We also emphasise the need to update models for the enrichment sources at higher metallicity, as the satisfactory recovery of the solar values of [(Cr,Mn,Co,Zn)/Fe] still presents a challenge. △ Less

Submitted 3 June, 2020; v1 submitted 13 November, 2019; originally announced November 2019.

Comments: 17 pages, 14 figures

Broadband spectroscopy of astrophysical ice analogs. I. Direct measurement of complex refractive index of CO ice using terahertz time-domain spectroscopy

Authors: B. M. Giuliano, A. A. Gavdush, B. Müller, K. I. Zaytsev, T. Grassi, A. V. Ivlev, M. E. Palumbo, G. A. Baratta, C. Scirè, G. A. Komandin, S. O. Yurchenko, P. Caselli

Abstract: Context: Reliable, directly measured optical properties of astrophysical ice analogs in the infrared (IR) and terahertz (THz) range are missing. These parameters are of great importance to model the dust continuum radiative transfer in dense and cold regions, here thick ice mantles are present, and are necessary for the interpretation of future observations planned in the far-IR region. Aims: Cohe… ▽ More Context: Reliable, directly measured optical properties of astrophysical ice analogs in the infrared (IR) and terahertz (THz) range are missing. These parameters are of great importance to model the dust continuum radiative transfer in dense and cold regions, here thick ice mantles are present, and are necessary for the interpretation of future observations planned in the far-IR region. Aims: Coherent THz radiation allows direct measurement of the complex dielectric function (refractive index) of astrophysically relevant ice species in the THz range. Methods: The time-domain waveforms and the frequency-domain spectra of reference samples of CO ice, deposited at a temperature of 28.5 K and annealed to 33 K at different thicknesses, have been recorded. A new algorithm is developed to reconstruct the real and imaginary parts of the refractive index from the time-domain THz data. Results: The complex refractive index in the wavelength range of 1 mm - 150 $μ$m (0.3 - 2.0 THz) has been determined for the studied ice samples, and compared with available data found in the literature. Conclusions: The developed algorithm of reconstructing the real and imaginary parts of the refractive index from the time-domain THz data enables, for the first time, the determination of optical properties of astrophysical ice analogs without using the Kramers-Kronig relations. The obtained data provide a benchmark to interpret the observational data from current ground based facilities as well as future space telescope missions, and have been used to estimate the opacities of the dust grains in presence of CO ice mantles. △ Less

Submitted 26 July, 2019; originally announced July 2019.

Journal ref: A&A 629, A112 (2019)