Showing 1–50 of 168 results for author: Röpke, F

Large-scale ordered magnetic fields generated in mergers of helium white dwarfs

Authors: Rüdiger Pakmor, Ingrid Pelisoli, Stephen Justham, Abinaya S. Rajamuthukumar, Friedrich K. Röpke, Fabian R. N. Schneider, Selma E. de Mink, Sebastian T. Ohlmann, Philipp Podsiadlowski, Javier Moran Fraile, Marco Vetter, Robert Andrassy

Abstract: Stellar mergers are one important path to highly magnetised stars. Mergers of two low-mass white dwarfs may create up to every third hot subdwarf star. The merging process is usually assumed to dramatically amplify magnetic fields. However, so far only four highly magnetised hot subdwarf stars have been found, suggesting a fraction of less than $1\%$. We present two high-resolution magnetohydrod… ▽ More Stellar mergers are one important path to highly magnetised stars. Mergers of two low-mass white dwarfs may create up to every third hot subdwarf star. The merging process is usually assumed to dramatically amplify magnetic fields. However, so far only four highly magnetised hot subdwarf stars have been found, suggesting a fraction of less than $1\%$. We present two high-resolution magnetohydrodynamical (MHD) simulations of the merger of two helium white dwarfs in a binary system with the same total mass of $0.6\,M_\odot$. We analyse one equal-mass merger with two $0.3\,M_\odot$ white dwarfs, and one unequal-mass merger with a $0.25\,M_\odot$ white dwarf and a $0.35\,M_\odot$ white dwarf. We simulate the inspiral, merger, and further evolution of the merger remnant for about $50$ rotations. We find efficient magnetic field amplification in both mergers via a small-scale dynamo, reproducing previous results of stellar merger simulations. The magnetic field saturates at similar strength for both simulations. We then identify a second phase of magnetic field amplification in both merger remnants that happens on a timescale of several tens of rotational periods of the merger remnant. This phase generates a large-scale ordered azimuthal field. We identify it as a large-scale dynamo driven by the magneto-rotational instability (MRI). Finally, we suggest that in the unequal-mass merger remnant, helium burning will eventually start in a shell around a cold core. The convection zone this generates will coincide with the region that contains most of the magnetic energy, probably erasing the strong, ordered field. The equal-mass merger remnant instead will probably ignite burning in the center, retaining its ordered field. Therefore, the mass ratio of the initial merger could be the selecting factor that decides if a merger remnant will stay highly magnetised long after the merger. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 12 pages, 9 figures, submitted to A&A, comments welcome

Gravitational-wave model for neutron star merger remnants with supervised learning

Authors: Theodoros Soultanis, Kiril Maltsev, Andreas Bauswein, Katerina Chatziioannou, Friedrich K. Roepke, Nikolaos Stergioulas

Abstract: We present a time-domain model for the gravitational waves emitted by equal-mass binary neutron star merger remnants for a fixed equation of state. We construct a large set of numerical relativity simulations for a single equation of state consistent with current constraints, totaling 157 equal-mass binary neutron star merger configurations. The gravitational-wave model is constructed using the su… ▽ More We present a time-domain model for the gravitational waves emitted by equal-mass binary neutron star merger remnants for a fixed equation of state. We construct a large set of numerical relativity simulations for a single equation of state consistent with current constraints, totaling 157 equal-mass binary neutron star merger configurations. The gravitational-wave model is constructed using the supervised learning method of K-nearest neighbor regression. As a first step toward developing a general model with supervised learning methods that accounts for the dependencies on equation of state and the binary masses of the system, we explore the impact of the size of the dataset on the model. We assess the accuracy of the model for a varied dataset size and number density in total binary mass. Specifically, we consider five training sets of $\{ 20,40, 60, 80, 100\}$ simulations uniformly distributed in total binary mass. We evaluate the resulting models in terms of faithfulness using a test set of 30 additional simulations that are not used during training and which are equidistantly spaced in total binary mass. The models achieve faithfulness with maximum values in the range of $0.980$ to $0.995$. We assess our models simulating signals observed by the three-detector network of Advanced LIGO-Virgo. We find that all models with training sets of size equal to or larger than $40$ achieve an unbiased measurement of the main gravitational-wave frequency. We confirm that our results do not depend qualitatively on the choice of the (fixed) equation of state. We conclude that training sets, with a minimum size of $40$ simulations, or a number density of approximately $11$ simulations per $0.1\,M_\odot$ of total binary mass, suffice for the construction of faithful templates for the post-merger signal for a single equation of state and equal-mass binaries (abbreviated). △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 20 pages, 12 figures, submitted to Phys. Rev. D

Including a Luminous Central Remnant in Radiative Transfer Simulations for Type Iax Supernovae

Authors: F. P. Callan, S. A. Sim, C. E. Collins, L. J. Shingles, F. Lach, F. K. Roepke, R. Pakmor, M. Kromer, S. Srivastav

Abstract: Type Iax supernovae (SNe Iax) are proposed to arise from deflagrations of Chandrasekhar mass white dwarfs (WDs). Previous deflagration simulations have achieved good agreement with the light curves and spectra of intermediate-luminosity and bright SNe Iax. However, the model light curves decline too quickly after peak, particularly in red optical and near-infrared (NIR) bands. Deflagration models… ▽ More Type Iax supernovae (SNe Iax) are proposed to arise from deflagrations of Chandrasekhar mass white dwarfs (WDs). Previous deflagration simulations have achieved good agreement with the light curves and spectra of intermediate-luminosity and bright SNe Iax. However, the model light curves decline too quickly after peak, particularly in red optical and near-infrared (NIR) bands. Deflagration models with a variety of ignition configurations do not fully unbind the WD, leaving a remnant polluted with $^{56}\mathrm{Ni}$. Emission from such a remnant may contribute to the luminosity of SNe Iax. Here we investigate the impact of adding a central energy source, assuming instantaneous powering by $^{56}\mathrm{Ni}$ decay in the remnant, in radiative transfer calculations of deflagration models. Including the remnant contribution improves agreement with the light curves of SNe Iax, particularly due to the slower post-maximum decline of the models. Spectroscopic agreement is also improved, with intermediate-luminosity and faint models showing greatest improvement. We adopt the full remnant $^{56}\mathrm{Ni}$ mass predicted for bright models, but good agreement with intermediate-luminosity and faint SNe Iax is only possible for remnant $^{56}\mathrm{Ni}$ masses significantly lower than those predicted. This may indicate that some of the $^{56}\mathrm{Ni}$ decay energy in the remnant does not contribute to the radiative luminosity but instead drives mass ejection, or that escape of energy from the remnant is significantly delayed. Future work should investigate the structure of remnants predicted by deflagration models and the potential roles of winds and delayed energy escape, as well as extend radiative transfer simulations to late times. △ Less

Submitted 19 April, 2024; v1 submitted 22 March, 2024; originally announced March 2024.

Comments: 17 pages, 6 figures. Lightcurves and spectra available at https://hesma.h-its.org

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 530, Issue 2, May 2024, Pages 1457 to 1473

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

Authors: G. Leidi, R. Andrassy, W. Barsukow, J. Higl, P. V. F. Edelmann, F. K. Röpke

Abstract: High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so t… ▽ More High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. In this work, we use our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We consider Mach numbers in the range from $10^{-3}$ to $10^{-1}$ in a well-posed, 2D, Kelvin--Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. We find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a "low-dissipation" Riemann solver and a sextic reconstruction scheme, (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same, (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate, (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow, (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep. △ Less

Submitted 26 February, 2024; originally announced February 2024.

Journal ref: A&A 686, A34 (2024)

Type Ia supernova explosion models are inherently multidimensional

Authors: R. Pakmor, I. R. Seitenzahl, A. J. Ruiter, S. A. Sim, F. K. Roepke, S. Taubenberger, R. Bieri, S. Blondin

Abstract: Theoretical and observational approaches to settling the important questions surrounding the progenitor systems and the explosion mechanism of normal Type Ia supernovae have thus far failed. With its unique capability to obtain continuous spectra through the near- and mid-infrared, JWST now offers completely new insights into Type Ia supernovae. In particular, observing them in the nebular phase a… ▽ More Theoretical and observational approaches to settling the important questions surrounding the progenitor systems and the explosion mechanism of normal Type Ia supernovae have thus far failed. With its unique capability to obtain continuous spectra through the near- and mid-infrared, JWST now offers completely new insights into Type Ia supernovae. In particular, observing them in the nebular phase allows us to directly see the central ejecta and thereby constrain the explosion mechanism. We aim to understand and quantify differences in the structure and composition of the central ejecta of various Type Ia supernova explosion models. We examined the currently most popular explosion scenarios using self-consistent multidimensional explosion simulations of delayed-detonation and pulsationally assisted, gravitationally confined delayed detonation Chandrasekhar-mass models and double-detonation sub-Chandrasekhar-mass and violent merger models. We find that the distribution of radioactive and stable nickel in the final ejecta, both observable in nebular spectra, are significantly different between different explosion scenarios. Therefore, comparing synthetic nebular spectra with JWST observations should allow us to distinguish between explosion models. We show that the explosion ejecta are inherently multidimensional for all models, and the Chandrasekhar-mass explosions simulated in spherical symmetry in particular lead to a fundamentally unphysical ejecta structure. Moreover, we show that radioactive and stable nickel cover a significant range of densities at a fixed velocity of the homologously expanding ejecta. Any radiation transfer postprocessing has to take these variations into account to obtain faithful synthetic observables; this will likely require multidimensional radiation transport simulations. △ Less

Submitted 26 April, 2024; v1 submitted 16 February, 2024; originally announced February 2024.

Comments: 7 pages, 2 figures, accepted by A&A, comments welcome

Journal ref: A&A 686, A227 (2024)

Morphology and Mach Number Dependence of Subsonic Bondi-Hoyle Accretion

Authors: Logan J. Prust, Hila Glanz, Lars Bildsten, Hagai B. Perets, Friedrich K. Roepke

Abstract: We carry out three-dimensional computations of the accretion rate onto an object (of size $R_{\rm sink}$ and mass $m$) as it moves through a uniform medium at a subsonic speed $v_{\infty}$. The object is treated as a fully-absorbing boundary (e.g. a black hole). In contrast to early conjectures, we show that when $R_{\rm sink}\ll R_{A}=2Gm/v^2$ the accretion rate is independent of $v_{\infty}$ and… ▽ More We carry out three-dimensional computations of the accretion rate onto an object (of size $R_{\rm sink}$ and mass $m$) as it moves through a uniform medium at a subsonic speed $v_{\infty}$. The object is treated as a fully-absorbing boundary (e.g. a black hole). In contrast to early conjectures, we show that when $R_{\rm sink}\ll R_{A}=2Gm/v^2$ the accretion rate is independent of $v_{\infty}$ and only depends on the entropy of the ambient medium, its adiabatic index, and $m$. Our numerical simulations are conducted using two different numerical schemes via the Athena++ and Arepo hydrodynamics solvers, which reach nearly identical steady-state solutions. We find that pressure gradients generated by the isentropic compression of the flow near the accretor are sufficient to suspend much of the surrounding gas in a near-hydrostatic equilibrium, just as predicted from the spherical Bondi-Hoyle calculation. Indeed, the accretion rates for steady flow match the Bondi-Hoyle rate, and are indicative of isentropic flow for subsonic motion where no shocks occur. We also find that the accretion drag may be predicted using the Safronov number, $Θ=R_{A}/R_{\rm sink}$, and is much less than the dynamical friction for sufficiently small accretors ($R_{\rm sink}\ll R_{A}$). △ Less

Submitted 15 February, 2024; originally announced February 2024.

Comments: 14 pages, 7 figures. Submitted to ApJ. Comments welcome!

Thermonuclear explosions as Type II supernovae

Authors: Alexandra Kozyreva, Javier Moran-Fraile, Alexander Holas, Vincent A. Bronner, Friedrich K. Roepke, Nikolay Pavlyuk, Alexey Mironov, Dmitriy Tsvetkov

Abstract: We consider a binary stellar system, in which a low-mass, of 0.6 Msun, carbon-oxygen white dwarf (WD) mergers with a degenerate helium core of 0.4 Msun of a red giant. We analyse the outcome of a merger within a common envelope (CE). We predict the observational properties of the resulting transient. We find that the double detonation of the WD, being a pure thermonuclear explosion and embedded in… ▽ More We consider a binary stellar system, in which a low-mass, of 0.6 Msun, carbon-oxygen white dwarf (WD) mergers with a degenerate helium core of 0.4 Msun of a red giant. We analyse the outcome of a merger within a common envelope (CE). We predict the observational properties of the resulting transient. We find that the double detonation of the WD, being a pure thermonuclear explosion and embedded into the hydrogen-rich CE, has a light curve with the distinct plateau shape, i.e. looks like a supernova (SN) Type IIP, with a duration of about 40 days. We find five observed SNe IIP: SN 2004dy, SN 2005af, SN 2005hd, SN 2007aa, and SN 2008bu, that match the V-band light curve of our models. Hence, we show that a thermonuclear explosion within a CE might be mistakenly identified as a SN IIP, which are believed to be an outcome of a core-collapse neutrino-driven explosion of a massive star. We discuss a number of diagnostics, that may help to distinguish this kind of a thermonuclear explosion from a core-collapse SN. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: 12 pp., accepted for publication in A & A

Evolution and final fate of massive post-common-envelope binaries

Authors: Dandan Wei, Fabian R. N. Schneider, Philipp Podsiadlowski, Eva Laplace, Friedrich K. Roepke, Marco Vetter

Abstract: Mergers of neutron stars (NSs) and black holes (BHs) are nowadays observed routinely thanks to gravitational-wave (GW) astronomy. In the isolated binary-evolution channel, a common-envelope (CE) phase of a red supergiant (RSG) and a compact object is crucial to sufficiently shrink the orbit and thereby enable a merger via GW emission. Here, we use the outcomes of two three-dimensional (3D) magneto… ▽ More Mergers of neutron stars (NSs) and black holes (BHs) are nowadays observed routinely thanks to gravitational-wave (GW) astronomy. In the isolated binary-evolution channel, a common-envelope (CE) phase of a red supergiant (RSG) and a compact object is crucial to sufficiently shrink the orbit and thereby enable a merger via GW emission. Here, we use the outcomes of two three-dimensional (3D) magneto-hydrodynamic CE simulations of an initially 10.0 solar-mass RSG with a 5.0 solar-mass BH and a 1.4 solar-mass NS, respectively, to explore the further evolution and final fate of the post-CE binaries. Notably, the 3D simulations reveal that the post-CE binaries are likely surrounded by circumbinary disks (CBDs), which contain substantial mass and angular momentum to influence the subsequent evolution. The binary systems in MESA modelling undergo another phase of mass transfer (MT) and we find that most donor stars do not explode in ultra-stripped supernovae (SNe), but rather in Type Ib/c SNe. The final orbits of our models with the BH companion are too wide, and NS kicks are actually required to sufficiently perturb the orbit and thus facilitate a merger via GW emission. Moreover, by exploring the influence of CBDs, we find that mass accretion from the disk widens the binary orbit, while CBD-binary resonant interactions can shrink the separation and increase the eccentricity depending on the disk mass and lifetime. Efficient resonant contractions may even enable a BH or NS to merge with the remnant He stars before a second SN explosion, which may be observed as gamma-ray burst-like transients, luminous fast blue optical transients and Thorne-Żytkow objects. For the surviving post-CE binaries, the CBD-binary interactions may significantly increase the GW-induced double compact merger fraction. We conclude that accounting for CBD may be crucial to better understand observed GW mergers. △ Less

Submitted 28 May, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Comments: 21 pages, 19 figures, 1 table. Accepted for publication on A&A

Going from 3D to 1D: A one-dimensional approach to common-envelope evolution

Authors: V. A. Bronner, F. R. N. Schneider, Ph. Podsiadlowski, F. K. Roepke

Abstract: The common-envelope (CE) phase is a crucial stage in binary star evolution because the orbital separation can shrink drastically while ejecting the envelope of a giant star. Three-dimensional (3D) hydrodynamic simulations of CE evolution are indispensable to learning about the mechanisms that play a role during the CE phase. While these simulations offer great insight, they are computationally exp… ▽ More The common-envelope (CE) phase is a crucial stage in binary star evolution because the orbital separation can shrink drastically while ejecting the envelope of a giant star. Three-dimensional (3D) hydrodynamic simulations of CE evolution are indispensable to learning about the mechanisms that play a role during the CE phase. While these simulations offer great insight, they are computationally expensive. We propose a one-dimensional (1D) model to simulate the CE phase within the stellar evolution code $\texttt{MESA}$ by using a parametric drag force prescription for dynamical drag and adding the released orbital energy as heat into the envelope. We compute CE events of a $0.97\,\mathrm{M}_\odot$ asymptotic giant-branch star and a point mass companion with mass ratios of 0.25, 0.50, and 0.75, and compare them to 3D simulations of the same setup. The 1D CE model contains two free parameters, which we demonstrate are both needed to fit the spiral-in behavior and the fraction of ejected envelope mass of the 1D method to the 3D simulations. For mass ratios of 0.25 and 0.50, we find good-fitting 1D simulations, while for a mass ratio of 0.75, we do not find a satisfactory fit to the 3D simulation as some of the assumptions in the 1D method are no longer valid. In all our simulations, we find that the released recombination energy is important to accelerate the envelope and drive the ejection. △ Less

Submitted 10 November, 2023; originally announced November 2023.

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

Faint calcium-rich transient from the double-detonation of a $0.6\,M_\odot$ carbon-oxygen white dwarf star

Authors: J. Moran-Fraile, A. Holas, F. K. Roepke, R. Pakmor, F. R. N. Schneider

Abstract: We have computed a three-dimensional hydrodynamic simulation of the merger between a massive ($0.4\,M_\odot$) helium white dwarf (He WD) and a low-mass ($0.6\,M_\odot$) carbon-oxygen white dwarf (CO WD). Despite the low mass of the primary, the merger triggers a thermonuclear explosion as a result of a double detonation, producing a faint transient and leaving no remnant behind. This type of event… ▽ More We have computed a three-dimensional hydrodynamic simulation of the merger between a massive ($0.4\,M_\odot$) helium white dwarf (He WD) and a low-mass ($0.6\,M_\odot$) carbon-oxygen white dwarf (CO WD). Despite the low mass of the primary, the merger triggers a thermonuclear explosion as a result of a double detonation, producing a faint transient and leaving no remnant behind. This type of event could also take place during common-envelope mergers whenever the companion is a CO WD and the core of the giant star has a sufficiently large He mass. The spectra show strong Ca lines throughout the first few weeks after the explosion. The explosion only yields $<0.01\,M_\odot$ of $^{56}$Ni, resulting in a low-luminosity SN Ia-like lightcurve that resembles the Ca-rich transients within this broad class of objects, with a peak magnitude of $M_\mathrm{bol} \approx -15.7\,$mag and a rather slow decline rate of $Δm_{15}^\mathrm{bol}\approx 1.5\,$mag. Both, its lightcurve-shape and spectral appearance, resemble the appearance of Ca-rich transients, suggesting such mergers as a possible progenitor scenario for this class of events. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: Submitted to A&A letters, posted on ArXiv after positive referee report. 7 pages, 4 figures

Self-consistent MHD simulation of jet launching in a neutron star - white dwarf merger

Authors: J. Moran-Fraile, F. K. Roepke, R. Pakmor, M. A. Aloy, S. T. Ohlmann, F. R. N. Schneider, G. Leidi

Abstract: The merger of a white dwarf (WD) and a neutron star (NS) is a relatively common event that will produce an observable electromagnetic signal. Furthermore, the compactness of these stellar objects makes them an interesting candidate for gravitational wave (GW) astronomy, potentially being in the frequency range of LISA and other missions. To date, three-dimensional simulations of these mergers have… ▽ More The merger of a white dwarf (WD) and a neutron star (NS) is a relatively common event that will produce an observable electromagnetic signal. Furthermore, the compactness of these stellar objects makes them an interesting candidate for gravitational wave (GW) astronomy, potentially being in the frequency range of LISA and other missions. To date, three-dimensional simulations of these mergers have not fully modelled the WD disruption, or have used lower resolutions and have not included magnetic fields even though they potentially shape the evolution of the merger remnant. In this work, we simulate the merger of a 1.4$M_\odot$ NS with a 1$M_\odot$ carbon oxygen WD in the magnetohydrodynamic moving mesh code \AREPO. We find that the disruption of the WD forms an accretion disk around the NS, and the subsequent accretion by the NS powers the launch of strongly magnetized, mildly relativistic jets perpendicular to the orbital plane. Although the exact properties of the jets could be altered by unresolved physics around the NS, the event could result in a transient with a larger luminosity than kilonovae. We discuss possible connections to fast blue optical transients (FBOTs) and long-duration gamma-ray bursts. We find that the frequency of GWs released during the merger is too high to be detectable by the LISA mission, but suitable for deci-hertz observatories such as LGWA, BBO or DECIGO. △ Less

Submitted 12 October, 2023; originally announced October 2023.

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

Journal ref: A&A 681, A41 (2024)

Potential of the Cherenkov Telescope Array for studying the young nearby supernova remnants

Authors: Dmitry Prokhorov, Shiu-Hang Lee, Shigehiro Nagataki, Jacco Vink, Donald Ellison, Gilles Ferrand, Daniel Patnaude, Friedrich Roepke, Ivo Seitenzahlh, Patrick Slane

Abstract: Modern imaging atmospheric Cherenkov telescopes have extensively observed young nearby supernova remnants (SNRs), with ages of about 1000 years or less, in the very-high-energy (VHE) gamma-ray band. These efforts resulted in the detection of VHE emission from three young SNRs - Cassiopeia A, Tycho, and SN 1006 - and provided significant evidence for emission from the more distant Kepler's SNR. How… ▽ More Modern imaging atmospheric Cherenkov telescopes have extensively observed young nearby supernova remnants (SNRs), with ages of about 1000 years or less, in the very-high-energy (VHE) gamma-ray band. These efforts resulted in the detection of VHE emission from three young SNRs - Cassiopeia A, Tycho, and SN 1006 - and provided significant evidence for emission from the more distant Kepler's SNR. However, many questions on the production of VHE gamma rays in these remnants remain unanswered. Using detailed physical models for Tycho's SNR based on the CR-hydro-NEI code and physically motivated models for the other young nearby remnants, we simulated observations with the Cherenkov Telescope Array (CTA) of these gamma-ray sources. We highlight properties of these remnants accessible for investigation with future CTA observations and discuss which questions are expected to be answered. △ Less

Submitted 4 October, 2023; originally announced October 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Report number: CTA-ICRC/2023/23

Turbulent dynamo action and its effects on the mixing at the convective boundary of an idealized oxygen-burning shell

Authors: G. Leidi, R. Andrassy, J. Higl, P. V. F. Edelmann, F. K. Röpke

Abstract: Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dyn… ▽ More Convection is one of the most important mixing processes in stellar interiors. Hydrodynamic mass entrainment can bring fresh fuel from neighboring stable layers into a convection zone, modifying the structure and evolution of the star. Under some conditions, strong magnetic fields can be sustained by the action of a turbulent dynamo, adding another layer of complexity and possibly altering the dynamics in the convection zone and at its boundaries. In this study, we used our fully compressible Seven-League Hydro code to run detailed and highly resolved three-dimensional magnetohydrodynamic simulations of turbulent convection, dynamo amplification, and convective boundary mixing in a simplified setup whose stratification is similar to that of an oxygen-burning shell in a star with an initial mass of $25\ M_\odot$. We find that the random stretching of magnetic field lines by fluid motions in the inertial range of the turbulent spectrum (i.e., a small-scale dynamo) naturally amplifies the seed field by several orders of magnitude in a few convective turnover timescales. During the subsequent saturated regime, the magnetic-to-kinetic energy ratio inside the convective shell reaches values as high as $0.33$, and the average magnetic field strength is ${\sim}10^{10}\,\mathrm{G}$. Such strong fields efficiently suppress shear instabilities, which feed the turbulent cascade of kinetic energy, on a wide range of spatial scales. The resulting convective flows are characterized by thread-like structures that extend over a large fraction of the convective shell. The reduced flow speeds and the presence of magnetic fields with strengths up to $60\%$ of the equipartition value at the upper convective boundary diminish the rate of mass entrainment from the stable layer by ${\approx}\,20\%$ as compared to the purely hydrodynamic case. △ Less

Submitted 29 September, 2023; originally announced September 2023.

Journal ref: A&A 679, A132 (2023)

Scalable stellar evolution forecasting: Deep learning emulation vs. hierarchical nearest neighbor interpolation

Authors: K. Maltsev, F. R. N. Schneider, F. K. Roepke, A. I. Jordan, G. A. Qadir, W. E. Kerzendorf, K. Riedmiller, P. van der Smagt

Abstract: Many astrophysical applications require efficient yet reliable forecasts of stellar evolution tracks. One example is population synthesis, which generates forward predictions of models for comparison with observations. The majority of state-of-the-art rapid population synthesis methods are based on analytic fitting formulae to stellar evolution tracks that are computationally cheap to sample stati… ▽ More Many astrophysical applications require efficient yet reliable forecasts of stellar evolution tracks. One example is population synthesis, which generates forward predictions of models for comparison with observations. The majority of state-of-the-art rapid population synthesis methods are based on analytic fitting formulae to stellar evolution tracks that are computationally cheap to sample statistically over a continuous parameter range. The computational costs of running detailed stellar evolution codes, such as MESA, over wide and densely sampled parameter grids are prohibitive, while stellar-age based interpolation in-between sparsely sampled grid points leads to intolerably large systematic prediction errors. In this work, we provide two solutions for automated interpolation methods that offer satisfactory trade-off points between cost-efficiency and accuracy. We construct a timescale-adapted evolutionary coordinate and use it in a two-step interpolation scheme that traces the evolution of stars from ZAMS all the way to the end of core helium burning while covering a mass range from ${0.65}$ to $300 \, \mathrm{M_\odot}$. The feedforward neural network regression model (first solution) that we train to predict stellar surface variables can make millions of predictions, sufficiently accurate over the entire parameter space, within tens of seconds on a 4-core CPU. The hierarchical nearest-neighbor interpolation algorithm (second solution) that we hard-code to the same end achieves even higher predictive accuracy, the same algorithm remains applicable to all stellar variables evolved over time, but it is two orders of magnitude slower. Our methodological framework is demonstrated to work on the MIST (Choi et al. 2016) data set. Finally, we discuss the prospective applications of these methods and provide guidelines for generalizing them to higher dimensional parameter spaces. △ Less

Submitted 27 October, 2023; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: Accepted at A&A

Journal ref: A&A 681, A86 (2024)

Helium as a signature of the double detonation in Type Ia supernovae

Authors: Christine E. Collins, Stuart A. Sim, Luke. J. Shingles, Sabrina Gronow, Friedrich K. Roepke, Ruediger Pakmor, Ivo R. Seitenzahl, Markus Kromer

Abstract: The double detonation is a widely discussed mechanism to explain Type Ia supernovae from explosions of sub-Chandrasekhar mass white dwarfs. In this scenario, a helium detonation is ignited in a surface helium shell on a carbon/oxygen white dwarf, which leads to a secondary carbon detonation. Explosion simulations predict high abundances of unburnt helium in the ejecta, however, radiative transfer… ▽ More The double detonation is a widely discussed mechanism to explain Type Ia supernovae from explosions of sub-Chandrasekhar mass white dwarfs. In this scenario, a helium detonation is ignited in a surface helium shell on a carbon/oxygen white dwarf, which leads to a secondary carbon detonation. Explosion simulations predict high abundances of unburnt helium in the ejecta, however, radiative transfer simulations have not been able to fully address whether helium spectral features would form. This is because helium can not be sufficiently excited to form spectral features by thermal processes, but can be excited by collisions with non-thermal electrons, which most studies have neglected. We carry out a full non-local thermodynamic equilibrium (non-LTE) radiative transfer simulation for an instance of a double detonation explosion model, and include a non-thermal treatment of fast electrons. We find a clear He I λ 10830 feature which is strongest in the first few days after explosion and becomes weaker with time. Initially this feature is blended with the Mg II λ 10927 feature but over time separates to form a secondary feature to the blue wing of the Mg II λ 10927 feature. We compare our simulation to observations of iPTF13ebh, which showed a similar feature to the blue wing of the Mg II λ 10927 feature, previously identified as C I. Our simulation shows a good match to the evolution of this feature and we identify it as high velocity He I λ 10830. This suggests that He I λ 10830 could be a signature of the double detonation scenario. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Comments: 7 pages, accepted by MNRAS

Towards a self-consistent model of the convective core boundary in upper-main-sequence stars

Authors: R. Andrassy, G. Leidi, J. Higl, P. V. F. Edelmann, F. R. N. Schneider, F. K. Roepke

Abstract: There is strong observational evidence that convective cores of intermediate-mass and massive main-sequence stars are substantially larger than standard stellar-evolution models predict. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-time-scale process th… ▽ More There is strong observational evidence that convective cores of intermediate-mass and massive main-sequence stars are substantially larger than standard stellar-evolution models predict. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-time-scale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number Seven-League Hydro (SLH) code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a $15$ M$_\odot$ zero-age main-sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor ranging from $10^3$ to $10^6$. After reaching thermal equilibrium, all of our models show a clear penetration layer. Its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. As the luminosity is decreased, the penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification. This structure corresponds to the adiabatic ,,step overshoot'' model often employed in stellar-evolution calculations. The thickness of the penetration layer slowly decreases with decreasing luminosity. Depending on how we extrapolate our 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from $0.09$ to $0.44$ pressure scale heights, which are broadly compatible with observations. △ Less

Submitted 8 July, 2023; originally announced July 2023.

Comments: 10 pages, 12 figures, submitted to A&A

Type Ia Supernova Explosions in Binary Systems: A Review

Authors: Zheng-Wei Liu, Friedrich K. Roepke, Zhanwen Han

Abstract: SNe Ia play a key role in the fields of astrophysics and cosmology. It is widely accepted that SNe Ia arise from thermonuclear explosions of WDs in binaries. However, there is no consensus on the fundamental aspects of the nature of SN Ia progenitors and their explosion mechanism. This fundamentally flaws our understanding of these important astrophysical objects. We outline the diversity of SNe I… ▽ More SNe Ia play a key role in the fields of astrophysics and cosmology. It is widely accepted that SNe Ia arise from thermonuclear explosions of WDs in binaries. However, there is no consensus on the fundamental aspects of the nature of SN Ia progenitors and their explosion mechanism. This fundamentally flaws our understanding of these important astrophysical objects. We outline the diversity of SNe Ia and the proposed progenitor models and explosion mechanisms. We discuss the recent theoretical and observational progress in addressing the SN Ia progenitor and explosion mechanism in terms of the observables at various stages of the explosion, including rates and delay times, pre-explosion companion stars, ejecta-companion interaction, early excess emission, early radio/X-ray emission from CSM interaction, surviving companions, late-time spectra and photometry, polarization signals, and SNR properties, etc. Despite the efforts from both the theoretical and observational side, the questions of how the WDs reach an explosive state and what progenitor systems are more likely to produce SNe Ia remain open. No single published model is able to consistently explain all observational features and the full diversity of SNe Ia. This may indicate that either a new progenitor paradigm or the improvement of current models is needed if all SNe Ia arise from the same origin. An alternative scenario is that different progenitor channels and explosion mechanisms contribute to SNe Ia. In the next decade, the ongoing campaigns with the JWST, Gaia and the ZTF, and upcoming extensive projects with the LSST and the SKA will allow us to conduct not only studies of individual SNe Ia in unprecedented detail but also systematic investigations for different subclasses of SNe Ia. This will advance theory and observations of SNe Ia sufficiently far to gain a deeper understanding of their origin and explosion mechanism. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: An invited review (accepted for publication in RAA)

Gravitational wave emission from dynamical stellar interactions

Authors: Javier Moran-Fraile, Fabian R. N. Schneider, Friedrich K. Roepke, Sebastian T. Ohlmann, Ruediger Pakmor, Theodoros Soultanis, Andreas Bauswein

Abstract: We are witnessing the dawn of gravitational wave (GW) astronomy. With currently available detectors, observations are restricted to GW frequencies in the range between ${\sim} 10\,\mathrm{Hz}$ and $10\,\mathrm{kHz}$, which covers the signals from mergers of compact objects. The launch of the space observatory LISA will open up a new frequency band for the detection of stellar interactions at lower… ▽ More We are witnessing the dawn of gravitational wave (GW) astronomy. With currently available detectors, observations are restricted to GW frequencies in the range between ${\sim} 10\,\mathrm{Hz}$ and $10\,\mathrm{kHz}$, which covers the signals from mergers of compact objects. The launch of the space observatory LISA will open up a new frequency band for the detection of stellar interactions at lower frequencies. In this work, we predict the shape and strength of the GW signals associated with common-envelope interaction and merger events in binary stars, and we discuss their detectability. Previous studies estimated these characteristics based on semi-analytical models. In contrast, we used detailed three-dimensional magnetohydrodynamic simulations to compute the GW signals. We show that for the studied models, the dynamical phase of common-envelope events and mergers between main-sequence stars lies outside of the detectability band of the LISA mission. We find, however, that the final stages of common-envelope interactions leading to mergers of the stellar cores fall into the frequency band in which the sensitivity of LISA peaks, making them promising candidates for detection. These detections can constrain the enigmatic common-envelope dynamics. Furthermore, future decihertz observatories such as DECIGO or BBO would also be able to observe this final stage and the post-merger signal, through which we might be able to detect the formation of Thorne-Żytkow objects. △ Less

Submitted 9 March, 2023; originally announced March 2023.

Comments: Accepted for publication in A&A, 12 pages, 8 figures

Journal ref: A&A 672, A9 (2023)

Simulations of common-envelope evolution in binary stellar systems: physical models and numerical techniques

Authors: Friedrich K. Roepke, Orsola De Marco

Abstract: When the primary star in a close binary system evolves into a giant and engulfs its companion, its core and the companion temporarily orbit each other inside a common envelope. Drag forces transfer orbital energy and angular momentum to the envelope material. Depending on the efficiency of this process, the envelope may be ejected leaving behind a tight remnant binary system of two stellar cores,… ▽ More When the primary star in a close binary system evolves into a giant and engulfs its companion, its core and the companion temporarily orbit each other inside a common envelope. Drag forces transfer orbital energy and angular momentum to the envelope material. Depending on the efficiency of this process, the envelope may be ejected leaving behind a tight remnant binary system of two stellar cores, or the cores merge retaining part of the envelope material. The exact outcome of common-envelope evolution is critical for in the formation of X-ray binaries, supernova progenitors, the progenitors of compact-object mergers that emit detectable gravitational waves, and many other objects of fundamental astrophysical relevance. The wide ranges of spatial and temporal timescales that characterize common-envelope interactions and the lack of spatial symmetries present a substantial challenge to generating consistent models. Therefore, these critical phases are one of the largest sources for uncertainty in classical treatments of binary stellar evolution. Three-dimensional hydrodynamic simulations of at least part of the common envelope interaction are the key to gain predictive power in modeling common-envelope evolution. We review the development of theoretical concepts and numerical approaches for such three-dimensional hydrodynamic simulations. The inherent multi-physics, multi-scale challenges have resulted in a wide variety of approximations and numerical techniques to be exercised on the problem. We summarize the simulations published to date and their main results. Given the recent rapid progress, a sound understanding of the physics of common-envelope interactions is within reach and thus there is hope that one of the remaining fundamental problems of stellar astrophysics may be solved before long. △ Less

Submitted 14 December, 2022; originally announced December 2022.

Comments: submitted to Living Reviews in Computational Astrophysics, comments welcome (can be taken into account in the resubmission if received within 3 days after this posting)

A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors

Authors: G. Leidi, C. Birke, R. Andrassy, J. Higl, P. V. F. Edelmann, G. Wiest, C. Klingenberg, F. K. Röpke

Abstract: Fully compressible magnetohydrodynamic (MHD) simulations are a fundamental tool for investigating the role of dynamo amplification in the generation of magnetic fields in deep convective layers of stars. The flows that arise in such environments are characterized by low (sonic) Mach numbers (M_son < 0.01 ). In these regimes, conventional MHD codes typically show excessive dissipation and tend to b… ▽ More Fully compressible magnetohydrodynamic (MHD) simulations are a fundamental tool for investigating the role of dynamo amplification in the generation of magnetic fields in deep convective layers of stars. The flows that arise in such environments are characterized by low (sonic) Mach numbers (M_son < 0.01 ). In these regimes, conventional MHD codes typically show excessive dissipation and tend to be inefficient as the Courant-Friedrichs-Lewy (CFL) constraint on the time step becomes too strict. In this work we present a new method for efficiently simulating MHD flows at low Mach numbers in a space-dependent gravitational potential while still retaining all effects of compressibility. The proposed scheme is implemented in the finite-volume Seven-League Hydro (SLH) code, and it makes use of a low-Mach version of the five-wave Harten-Lax-van Leer discontinuities (HLLD) solver to reduce numerical dissipation, an implicit-explicit time discretization technique based on Strang splitting to overcome the overly strict CFL constraint, and a well-balancing method that dramatically reduces the magnitude of spatial discretization errors in strongly stratified setups. The solenoidal constraint on the magnetic field is enforced by using a constrained transport method on a staggered grid. We carry out five verification tests, including the simulation of a small-scale dynamo in a star-like environment at M_son ~ 0.001 . We demonstrate that the proposed scheme can be used to accurately simulate compressible MHD flows in regimes of low Mach numbers and strongly stratified setups even with moderately coarse grids. △ Less

Submitted 4 October, 2022; originally announced October 2022.

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

Double detonations: variations in Type Ia supernovae due to different core and He shell masses -- II: synthetic observables

Authors: Christine E. Collins, Sabrina Gronow, Stuart A. Sim, Friedrich K. Roepke

Abstract: Double detonations of sub-Chandrasekhar mass white dwarfs are a promising explosion scenario for Type Ia supernovae, whereby a detonation in a surface helium shell triggers a secondary detonation in a carbon-oxygen core. Recent work has shown that low mass helium shell models reproduce observations of normal SNe Ia. We present 3D radiative transfer simulations for a suite of 3D simulations of the… ▽ More Double detonations of sub-Chandrasekhar mass white dwarfs are a promising explosion scenario for Type Ia supernovae, whereby a detonation in a surface helium shell triggers a secondary detonation in a carbon-oxygen core. Recent work has shown that low mass helium shell models reproduce observations of normal SNe Ia. We present 3D radiative transfer simulations for a suite of 3D simulations of the double detonation explosion scenario for a range of shell and core masses. We find light curves broadly able to reproduce the faint end of the width-luminosity relation shown by SNe Ia, however, we find that all of our models show extremely red colours, not observed in normal SNe Ia. This includes our lowest mass helium shell model. We find clear Ti II absorption features in the model spectra, which would lead to classification as peculiar SNe Ia, as well as line blanketing in some lines of sight by singly ionised Cr and Fe-peak elements. Our radiative transfer simulations show that these explosion models remain promising to explain peculiar SNe Ia. Future full non-LTE simulations may improve the agreement of these explosion models with observations of normal SNe Ia. △ Less

Submitted 5 July, 2023; v1 submitted 9 September, 2022; originally announced September 2022.

Comments: Spectra available at https://doi.org/10.5281/zenodo.7997388

General relativistic moving-mesh hydrodynamics simulations with AREPO and applications to neutron star mergers

Authors: Georgios Lioutas, Andreas Bauswein, Theodoros Soultanis, Rüdiger Pakmor, Volker Springel, Friedrich K. Röpke

Abstract: We implement general relativistic hydrodynamics in the moving-mesh code AREPO. We also couple a solver for the Einstein field equations employing the conformal flatness approximation. The implementation is validated by evolving isolated static neutron stars using a fixed metric or a dynamical spacetime. In both tests the frequencies of the radial oscillation mode match those of independent calcula… ▽ More We implement general relativistic hydrodynamics in the moving-mesh code AREPO. We also couple a solver for the Einstein field equations employing the conformal flatness approximation. The implementation is validated by evolving isolated static neutron stars using a fixed metric or a dynamical spacetime. In both tests the frequencies of the radial oscillation mode match those of independent calculations. We run the first moving-mesh simulation of a neutron star merger. The simulation includes a scheme to adaptively refine or derefine cells and thereby adjusting the local resolution dynamically. The general dynamics are in agreement with independent smoothed particle hydrodynamics and static-mesh simulations of neutron star mergers. Coarsely comparing, we find that dynamical features like the post-merger double-core structure or the quasi-radial oscillation mode persist on longer time scales, possibly reflecting a low numerical diffusivity of our method. Similarly, the post-merger gravitational wave emission shows the same features as observed in simulations with other codes. In particular, the main frequency of the post-merger phase is found to be in good agreement with independent results for the same binary system, while, in comparison, the amplitude of the post-merger gravitational wave signal falls off slower, i.e. the post-merger oscillations are less damped. The successful implementation of general relativistic hydrodynamics in the moving-mesh AREPO code, including a dynamical spacetime evolution, provides a fundamentally new tool to simulate general relativistic problems in astrophysics. △ Less

Submitted 9 January, 2024; v1 submitted 8 August, 2022; originally announced August 2022.

Comments: 24 pages, 18 figures, accepted for publication in MNRAS

Long-term evolution of post-explosion Helium-star Companions of Type Iax Supernovae

Authors: Yaotian Zeng, Zheng-Wei Liu, Alexander Heger, Curtis McCully, Friedrich K. Röpke, Zhanwen Han

Abstract: Supernovae of Type Iax (SNe Iax) are an accepted faint subclass of hydrogen-free supernovae. Their origin, the nature of the progenitor systems, however, is an open question. Recent studies suggest that the weak deflagration explosion of a near-Chandrasekhar-mass white dwarf in a binary system with a helium star donor could be the origin of SNe Iax. In this scenario, the helium star donor is expec… ▽ More Supernovae of Type Iax (SNe Iax) are an accepted faint subclass of hydrogen-free supernovae. Their origin, the nature of the progenitor systems, however, is an open question. Recent studies suggest that the weak deflagration explosion of a near-Chandrasekhar-mass white dwarf in a binary system with a helium star donor could be the origin of SNe Iax. In this scenario, the helium star donor is expected to survive the explosion. We use the one-dimensional stellar evolution codes \textsc{MESA} and \textsc{Kepler} to follow the post-impact evolution of the surviving helium companion stars. The stellar models are based on our previous hydrodynamical simulations of ejecta-donor interaction, and we explore the observational characteristics of these surviving helium companions. We find that the luminosities of the surviving helium companions increase significantly after the impact: They could vary from $2\mathord,500\,\mathrm{L_{\odot}}$ to $16\mathord,000\,\mathrm{L_{\odot}}$ for a Kelvin-Helmholtz timescale of about $10^{4}\,\mathrm{yr}$. After the star reaches thermal equilibrium, it evolves as an O-type hot subdwarf (sdO) star and continues its evolution along the evolutionary track of a normal sdO star with the same mass. Our results will help to identify the surviving helium companions of SNe Iax in future observations and to place new constraints on their progenitor models. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 13 pages, 7 figures, accepted for publication in ApJ

Modelling the ionisation state of Type Ia supernovae in the nebular-phase

Authors: Luke J. Shingles, Andreas Flörs, Stuart A. Sim, Christine E. Collins, Friedrich K. Roepke, Ivo R. Seitenzahl, Ken J. Shen

Abstract: The nebular spectra of Type Ia supernovae ($\gtrapprox$ 100 days after explosion) consist mainly of emission lines from singly- and doubly-ionised Fe-group nuclei. However, theoretical models for many scenarios predict that non-thermal ionisation leads to multiply-ionised species whose recombination photons ionise and deplete Fe$^{+}$ , resulting in negligible [Fe II] emission. We investigate a me… ▽ More The nebular spectra of Type Ia supernovae ($\gtrapprox$ 100 days after explosion) consist mainly of emission lines from singly- and doubly-ionised Fe-group nuclei. However, theoretical models for many scenarios predict that non-thermal ionisation leads to multiply-ionised species whose recombination photons ionise and deplete Fe$^{+}$ , resulting in negligible [Fe II] emission. We investigate a method to determine the collisional excitation conditions from [Fe II] line ratios independently from the ionisation state and find that it cannot be applied to highly-ionised models due to the influence of recombination cascades on Fe$^{+}$ level populations. When the ionisation state is artificially lowered, the line ratios (and excitation conditions) are too similar to distinguish between explosion scenarios. We investigate changes to the treatment of non-thermal energy deposition as a way to reconcile over-ionised theoretical models with observations and find that a simple work function approximation provides closer agreement with the data for sub-Mch models than a detailed Spencer-Fano treatment with widely-used cross section data. To quantify the magnitude of additional heating processes that would be required to sufficiently reduce ionisation from fast leptons, we artificially boost the rate of energy loss to free electrons. We find that the equivalent of as much as an eight times increase to the plasma loss rate would be needed to reconcile the sub-Mch model with observed spectra. Future studies could distinguish between reductions in the non-thermal ionisation rates and increased recombination rates, such as by clumping. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 13 pages. Accepted by MNRAS

On the fate of the secondary white dwarf in double-degenerate double-detonation Type Ia supernovae

Authors: R. Pakmor, F. P. Callan, C. E. Collins, S. E. de Mink, A. Holas, W. E. Kerzendorf, M. Kromer, P. G. Neunteufel, John T. O'Brien, F. K. Roepke, A. J. Ruiter, I. R. Seitenzahl, Luke J. Shingles, S. A. Sim, S. Taubenberger

Abstract: The progenitor systems and explosion mechanism of Type Ia supernovae are still unknown. Currently favoured progenitors include double-degenerate systems consisting of two carbon-oxygen white dwarfs with thin helium shells. In the double-detonation scenario, violent accretion leads to a helium detonation on the more massive primary white dwarf that turns into a carbon detonation in its core and exp… ▽ More The progenitor systems and explosion mechanism of Type Ia supernovae are still unknown. Currently favoured progenitors include double-degenerate systems consisting of two carbon-oxygen white dwarfs with thin helium shells. In the double-detonation scenario, violent accretion leads to a helium detonation on the more massive primary white dwarf that turns into a carbon detonation in its core and explodes it. We investigate the fate of the secondary white dwarf, focusing on changes of the ejecta and observables of the explosion if the secondary explodes as well rather than survives. We simulate a binary system of a $1.05\,\mathrm{M_\odot}$ and a $0.7\,\mathrm{M_\odot}$ carbon-oxygen white dwarf with $0.03\,\mathrm{M_\odot}$ helium shells each. We follow the system self-consistently from inspiral to ignition, through the explosion, to synthetic observables. We confirm that the primary white dwarf explodes self-consistently. The helium detonation around the secondary white dwarf, however, fails to ignite a carbon detonation. We restart the simulation igniting the carbon detonation in the secondary white dwarf by hand and compare the ejecta and observables of both explosions. We find that the outer ejecta at $v~>~15\,000$\,km\,s$^{-1}$ are indistinguishable. Light curves and spectra are very similar until $\sim~40\,\mathrm{d}$ after explosion and the ejecta are much more spherical than violent merger models. The inner ejecta differ significantly slowing down the decline rate of the bolometric light curve after maximum of the model with a secondary explosion by $\sim20$ per cent. We expect future synthetic 3D nebular spectra to confirm or rule out either model. △ Less

Submitted 25 October, 2022; v1 submitted 28 March, 2022; originally announced March 2022.

Comments: 12 pages, 7 figures, accepted to MNRAS, comments welcome

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

Signatures of a surviving helium-star companion in Type Ia supernovae and constraints on the progenitor companion of SN 2011fe

Authors: Zheng-Wei Liu, Friedrich K. Roepke, Yaotian Zeng

Abstract: Single-degenerate (SD) binary systems composed of a white dwarf and a non-degenerate helium (He)-star companion have been proposed as the potential progenitors of Type Ia supernovae (SNe Ia). The He-star companions are expected to survive the SN Ia explosion in this SD progenitor model. In the present work, we map the surviving He-star companion models computed from our previous three-dimensional… ▽ More Single-degenerate (SD) binary systems composed of a white dwarf and a non-degenerate helium (He)-star companion have been proposed as the potential progenitors of Type Ia supernovae (SNe Ia). The He-star companions are expected to survive the SN Ia explosion in this SD progenitor model. In the present work, we map the surviving He-star companion models computed from our previous three-dimensional hydrodynamical simulations of ejecta-companion interaction into the one-dimensional stellar evolution code MESA to follow their long-term evolution to make predictions on their post-impact observational properties, which can be helpful for searches of such surviving He-star companions in future observations. By comparing with the very late-epoch light curve of the best observed SN Ia, SN 2011fe, we find that our surviving He-star companions become significantly more luminous than SN 2011fe about 1000d after the maximum light. This suggests that a He star is very unlikely to be a companion to the progenitor of SN 2011fe. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: 14 pages, 8 figures, 1 table, accepted for publication in ApJ

Models of pulsationally assisted gravitationally confined detonations with different ignition conditions

Authors: F. Lach, F. P. Callan, S. A. Sim, F. K. Roepke

Abstract: Over the past decades, many explosion scenarios for Type Ia supernovae have been proposed and investigated including various combinations of deflagrations and detonations in white dwarfs of different masses up to the Chandrasekhar mass. One of these is the gravitationally confined detonation model. In this case a weak deflagration burns to the surface, wraps around the bound core, and collides at… ▽ More Over the past decades, many explosion scenarios for Type Ia supernovae have been proposed and investigated including various combinations of deflagrations and detonations in white dwarfs of different masses up to the Chandrasekhar mass. One of these is the gravitationally confined detonation model. In this case a weak deflagration burns to the surface, wraps around the bound core, and collides at the antipode. A subsequent detonation is then initiated in the collision area. Since the parameter space for this scenario, that is, varying central densities and ignition geometries, has not been studied in detail, we used pure deflagration models of a previous parameter study dedicated to Type Iax supernovae as initial models to investigate the gravitationally confined detonation scenario. We aim to judge whether this channel can account for one of the many subgroups of Type Ia supernovae, or even normal events. To this end, we employed a comprehensive pipeline for three-dimensional Type Ia supernova modeling that consists of hydrodynamic explosion simulations, nuclear network calculations, and radiative transfer. The observables extracted from the radiative transfer are then compared to observed light curves and spectra. The study produces a wide range in masses of synthesized 56 Ni ranging from 0.257 to 1.057 $M_\odot$ , and, thus, can potentially account for subluminous as well as overluminous Type Ia supernovae in terms of brightness. However, a rough agreement with observed light curves and spectra can only be found for 91T-like objects. Although several discrepancies remain, we conclude that the gravitationally confined detonation model cannot be ruled out as a mechanism to produce 91T-like objects. However, the models do not provide a good explanation for either normal Type Ia supernovae or Type Iax supernovae. △ Less

Submitted 19 March, 2022; v1 submitted 29 November, 2021; originally announced November 2021.

Comments: 15 pages, 12 figures, published in A&A, corrected reference

Journal ref: A&A 659, A27 (2022)

From 3D hydrodynamic simulations of common-envelope interaction to gravitational-wave mergers

Authors: Melvin M. Moreno, Fabian R. N. Schneider, Friedrich K. Roepke, Sebastian T. Ohlmann, Ruediger Pakmor, Philipp Podsiadlowski, Christian Sand

Abstract: Modeling the evolution of progenitors of gravitational-wave merger events in binary stars faces two major uncertainties: the common-envelope phase and supernova kicks. These two processes are critical for the final orbital configuration of double compact-object systems with neutron stars and black holes. Predictive one-dimensional models of common-envelope interaction are lacking and multidimensio… ▽ More Modeling the evolution of progenitors of gravitational-wave merger events in binary stars faces two major uncertainties: the common-envelope phase and supernova kicks. These two processes are critical for the final orbital configuration of double compact-object systems with neutron stars and black holes. Predictive one-dimensional models of common-envelope interaction are lacking and multidimensional simulations are challenged by the vast range of relevant spatial and temporal scales. Here, we present three-dimensional hydrodynamic simulations of the common-envelope interaction of an initially $10\,M_{\odot}$ red supergiant primary star with a black-hole and a neutron-star companion. We show that the high-mass regime is accessible to full ab-initio simulations. Nearly complete envelope ejection is reached assuming that all recombination energy still available at the end of our simulation continues to help unbinding the envelope. In contrast to previous assumptions, we find that the dynamical plunge-in of both companions terminates at orbital separations too wide for gravitational waves to merge the systems in a Hubble time. We discuss the further evolution of the system based on analytical estimates. A subsequent mass-transfer episode from the remaining $3\,M_{\odot}$ core of the supergiant to the compact companion does not shrink the orbit sufficiently either. A neutron-star--neutron-star and neutron-star--black-hole merger is still expected for a fraction of the systems if the supernova kick aligns favorably with the orbital motion. For double neutron star (neutron-star--black-hole) systems we estimate mergers in about $9 \%$ ($1 \%$) of cases while about $77 \%$ ($94 \%$) of binaries are disrupted, i.e., supernova kicks actually enable gravitational-wave mergers in our cases; however, we expect a reduction in predicted gravitational-wave merger events. (abbr.) △ Less

Submitted 13 July, 2022; v1 submitted 23 November, 2021; originally announced November 2021.

Comments: 16 pages, 11 figures, accepted by A&A

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

Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes

Authors: R. Andrassy, J. Higl, H. Mao, M. Mocák, D. G. Vlaykov, W. D. Arnett, I. Baraffe, S. W. Campbell, T. Constantino, P. V. F. Edelmann, T. Goffrey, T. Guillet, F. Herwig, R. Hirschi, L. Horst, G. Leidi, C. Meakin, J. Pratt, F. Rizzuti, F. K. Roepke, P. Woodward

Abstract: Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility o… ▽ More Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within $\lesssim 3σ$ of the mean of all simulations on a given grid ($128^3$ and $256^3$ grid cells), where $σ$ describes the statistical variation due to the flow's time dependence. They also agree well with a $512^3$ reference run. The $128^3$ and $256^3$ simulations agree within $9\%$ and $4\%$, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study. △ Less

Submitted 26 January, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 20 pages, 19 figures, new Fig. B.2, new clarifications, electronic materials updated, accepted for publication in A&A

Journal ref: A&A 659, A193 (2022)

Bipolar planetary nebulae from common envelope evolution of binary stars

Authors: P. A. Ondratschek, F. K. Roepke, F. R. N. Schneider, C. Fendt, C. Sand, S. T. Ohlmann, R. Pakmor, V. Springel

Abstract: Asymmetric shapes and evidence for binary central stars suggest a common-envelope origin for many bipolar planetary nebulae. The bipolar components of the nebulae are observed to expand faster than the rest and the more slowly expanding material has been associated with the bulk of the envelope ejected during the common-envelope phase of a stellar binary system. Common-envelope evolution in genera… ▽ More Asymmetric shapes and evidence for binary central stars suggest a common-envelope origin for many bipolar planetary nebulae. The bipolar components of the nebulae are observed to expand faster than the rest and the more slowly expanding material has been associated with the bulk of the envelope ejected during the common-envelope phase of a stellar binary system. Common-envelope evolution in general remains one of the biggest uncertainties in binary star evolution and the origin of the fast outflow has not been explained satisfactorily. We perform three-dimensional magnetohydrodynamic simulations of common-envelope interaction with the moving-mesh code AREPO. Starting from the plunge-in of the companion into the envelope of an asymptotic giant branch star and covering hundreds of orbits of the binary star system, we are able to follow the evolution to complete envelope ejection. We find that magnetic fields are strongly amplified in two consecutive episodes. First, when the companion spirals in the envelope and, second, when it forms a contact binary with the core of the former giant star. In the second episode, a magnetically-driven, high-velocity outflow of gas is launched self-consistently in our simulations. The outflow is bipolar and the gas is additionally collimated by the ejected common envelope. The resulting structure reproduces typical morphologies and velocities observed in young planetary nebulae. We propose that the magnetic driving mechanism is a universal consequence of common envelope interaction responsible for a substantial fraction of observed planetary nebulae. Such a mechanism likely also exists in the common-envelope phase of other binary stars that lead to the formation of Type Ia supernovae, X-ray binaries and gravitational-wave merger events. △ Less

Submitted 17 March, 2022; v1 submitted 25 October, 2021; originally announced October 2021.

Comments: 9 pages, 11 figures; accepted for publication by A&A

Journal ref: A&A 660, L8 (2022)

Observational signatures of the surviving donor star in the double detonation model of Type Ia supernovae

Authors: Zheng-Wei Liu, Friedrich K. Roepke, Yaotian Zeng, Alexander Heger

Abstract: The sub-Chandrasekhar mass double-detonation (DDet) scenario is a contemporary model for SNe Ia. The donor star in the DDet scenario is expected to survive the explosion and to be ejected at the high orbital velocity of a compact binary system. For the first time, we consistently perform 3D hydrodynamical simulations of the interaction of SN ejecta with a helium (He) star companion within the DDet… ▽ More The sub-Chandrasekhar mass double-detonation (DDet) scenario is a contemporary model for SNe Ia. The donor star in the DDet scenario is expected to survive the explosion and to be ejected at the high orbital velocity of a compact binary system. For the first time, we consistently perform 3D hydrodynamical simulations of the interaction of SN ejecta with a helium (He) star companion within the DDet scenario. We map the outcomes of 3D impact simulations into 1D stellar evolution codes and follow the long-term evolution of the surviving He-star companions. Our main goal is to provide the post-impact observable signatures of surviving He-star companions of DDet SNe Ia, which will support the search for such companions in future observations. We find that our surviving He-star companions become significantly overluminous for about 1e6 yr during the thermal re-equilibration phase. After the star re-establishes thermal equilibrium, its observational properties are not sensitive to the details of the ejecta-donor interaction. We apply our results to hypervelocity star US 708, which is the fastest unbound star in our Galaxy, travelling with a velocity of about 1200 km/s, making it natural candidate for an ejected donor remnant of a DDet SN Ia. We find that a He-star donor with an initial mass of >0.5 Msun is needed to explain the observed properties of US 708. Based on our detailed binary evolution calculations, however, the progenitor system with such a massive He-star donor cannot get close enough at the moment of SN explosion to explain the high velocity of US 708. Instead, if US 708 is indeed the surviving He-star donor of a DDet SN~Ia, it would require the entire pre-SN progenitor binary to travel at a velocity of about 400 km/s. It could, for example, have been ejected from a globular cluster in the direction of the current motion of the surviving donor star. △ Less

Submitted 21 September, 2021; originally announced September 2021.

Comments: 11 pages, 9 figures, 1 table, accepted for publication in A&A

Journal ref: A&A 654, A103 (2021)

Type Iax supernovae from deflagrations in Chandrasekhar mass white dwarfs

Authors: F. Lach, F. P. Callan, D. Bubeck, F. K. Roepke, S. A. Sim, M. Schrauth, S. T. Ohlmann, M. Kromer

Abstract: Due to the increasing number of observations Type Ia supernovae are nowadays regarded as a heterogeneous class of objects consisting of several subclasses. One of the largest of these is the class of Type Iax supernovae (SNe Iax) which have been suggested to originate from pure deflagrations in CO Chandrasekhar-mass white dwarfs (WDs). Although a few deflagration studies have been carried out, the… ▽ More Due to the increasing number of observations Type Ia supernovae are nowadays regarded as a heterogeneous class of objects consisting of several subclasses. One of the largest of these is the class of Type Iax supernovae (SNe Iax) which have been suggested to originate from pure deflagrations in CO Chandrasekhar-mass white dwarfs (WDs). Although a few deflagration studies have been carried out, the full diversity of the class is not captured yet. We therefore present a parameter study of single-spot ignited deflagrations with varying ignition locations, central densities, metallicities and compositions. We also explore a rigidly rotating progenitor and carry out 3D hydrodynamic simulations, nuclear network calculations and radiative transfer. The new models extend the range in brightness covered by previous studies to the lower end. Our explosions produce $^{56}$Ni masses from $5.8 \times 10^{-3}$ to $9.2 \times 10^{-2}\,M_\odot$. In spite of the wide exploration of the parameter space the main characteristics of the models are primarily driven by the mass of $^{56}$Ni. Secondary parameters have too little impact to explain the observed trend among faint SNe~Iax. We report kick velocities of the bound explosion remnants from $6.9$ to $369.8\,$km$\,s^{-1}$. The wide exploration of the parameter space and viewing-angle effects in the radiative transfer lead to a significant spread in the synthetic observables. The trends towards the faint end of the class are, however, not reproduced. This motivates a quantification of the systematic uncertainties in the modeling procedure and the influence of the $^{56}$Ni-rich bound remnant. While the pure deflagration scenario remains a favorable explanation for bright and intermediate luminosity SNe~Iax, the possibility that SNe~Iax do not consist of a single explosion scenario needs to be considered. △ Less

Submitted 3 December, 2021; v1 submitted 7 September, 2021; originally announced September 2021.

Comments: 27 pages, 17 figures, 3 tables, accepted for publication in A&A

Journal ref: A&A 658, A179 (2022)

Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star

Authors: L. Horst, R. Hirschi, P. V. F. Edelmann, R. Andrassy, F. K. Roepke

Abstract: Context. Multidimensional hydrodynamic simulations of convection in stellar interiors are numerically challenging, especially for flows at low Mach numbers. Methods. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional si… ▽ More Context. Multidimensional hydrodynamic simulations of convection in stellar interiors are numerically challenging, especially for flows at low Mach numbers. Methods. We explore the benefits of using a low-Mach hydrodynamic flux solver and demonstrate its usability for simulations in the astrophysical context. The time-implicit Seven-League Hydro (SLH) code was used to perform multidimensional simulations of convective helium shell burning based on a 25 M$_\odot$ star model. The results obtained with the low-Mach AUSM$^{+}$-up solver were compared to results when using its non low-Mach variant AUSM$_\mathrm{B}^{+}$-up. We applied well-balancing of the gravitational source term to maintain the initial hydrostatic background stratification. The computational grids have resolutions ranging from $180 \times 90^2$ to $810 \times 540^2$ cells and the nuclear energy release was boosted by factors of $3 \times 10^3$, $1 \times 10^4$, and $3 \times 10^4$ to study the dependence of the results on these parameters. Results. The boosted energy input results in convection at Mach numbers in the range of $10^{-2}$ to $10^{-3}$. Standard mixing-length theory (MLT) predicts convective velocities of about $1.6 \times 10^{-4}$ if no boosting is applied. Simulations with AUSM$^{+}$-up show a Kolmogorov-like inertial range in the kinetic energy spectrum that extends further toward smaller scales compared with its non low-Mach variant. The kinetic energy dissipation of the AUSM$^{+}$-up solver already converges at a lower resolution compared to AUSM$^{+}_{\mathrm{B}}$ -up. The extracted entrainment rates at the boundaries of the convection zone are well represented by the bulk Richardson entrainment law and the corresponding fitting parameters are in agreement with published results for carbon shell burning. △ Less

Submitted 5 July, 2021; originally announced July 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 653, A55 (2021)

Metallicity-dependent nucleosynthetic yields of Type Ia supernovae originating from double detonations of sub-M$_{\text{Ch}}$ white dwarfs

Authors: Sabrina Gronow, Benoit Cote, Florian Lach, Ivo R. Seitenzahl, Christine E. Collins, Stuart A. Sim, Friedrich K. Roepke

Abstract: Double detonations in sub-Chandrasekhar mass carbon-oxygen white dwarfs with helium shell are a potential explosion mechanism for a Type Ia supernova (SNe Ia). It comprises a shell detonation and subsequent core detonation. The focus of our study is on the effect of the progenitor metallicity on the nucleosynthetic yields. For this, we compute and analyse a set of eleven different models with vary… ▽ More Double detonations in sub-Chandrasekhar mass carbon-oxygen white dwarfs with helium shell are a potential explosion mechanism for a Type Ia supernova (SNe Ia). It comprises a shell detonation and subsequent core detonation. The focus of our study is on the effect of the progenitor metallicity on the nucleosynthetic yields. For this, we compute and analyse a set of eleven different models with varying core and shell masses at four different metallicities each. This results in a total of 44 models at metallicities between 0.01$Z_\odot$ and 3$Z_\odot$. Our models show a strong impact of the metallicity in the high density regime. The presence of $^{22}$Ne causes a neutron-excess which shifts the production from $^{56}$Ni to stable isotopes such as $^{54}$Fe and $^{58}$Ni in the $α$-rich freeze-out regime. The isotopes of the metallicity implementation further serve as seed nuclei for additional reactions in the shell detonation. Most significantly, the production of $^{55}$Mn increases with metallicity confirming the results of previous work. A comparison of elemental ratios relative to iron shows a relatively good match to solar values for some models. Super-solar values are reached for Mn at 3$Z_\odot$ and solar values in some models at $Z_\odot$. This indicates that the required contribution of SNe Ia originating from Chandrasekhar mass WDs can be lower than estimated in orevious work to reach solar values of [Mn/Fe] at [Fe/H]$=0$. Our galactic chemical evolution models suggest that SNe Ia from sub-Chandrasekhar mass white dwarfs, along with core-collapse supernovae, could account for more than 80% of the solar Mn abundance. Using metallicity-dependent SN Ia yields helps to reproduce the upward trend of [Mn/Fe] as a function of metallicity for the solar neighborhood. These chemical evolution predictions, however, depend on the massive star yields adopted in the calculations. △ Less

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

Comments: accepted for publication by A&A

Journal ref: A&A 656, A94 (2021)

Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers

Authors: P. V. F. Edelmann, L. Horst, J. P. Berberich, R. Andrassy, J. Higl, G. Leidi, C. Klingenberg, F. K. Roepke

Abstract: Accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and the pressure gradient, such simulations place heavy demands on numerical hydrodynamics schemes. We demonstrate how discretization errors on grids of reasonable s… ▽ More Accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and the pressure gradient, such simulations place heavy demands on numerical hydrodynamics schemes. We demonstrate how discretization errors on grids of reasonable size can lead to spurious flows orders of magnitude faster than the physical flow. Well-balanced numerical schemes can deal with this problem. Three such schemes were applied in the implicit, finite-volume Seven-League Hydro (SLH) code in combination with a low-Mach-number numerical flux function. We compare how the schemes perform in four numerical experiments addressing some of the challenges imposed by typical problems in stellar hydrodynamics. We find that the $α$-$β$ and deviation well-balancing methods can accurately maintain hydrostatic solutions provided that gravitational potential energy is included in the total energy balance. They accurately conserve minuscule entropy fluctuations advected in an isentropic stratification, which enables the methods to reproduce the expected scaling of convective flow speed with the heating rate. The deviation method also substantially increases accuracy of maintaining stationary orbital motions in a Keplerian disk on long timescales. The Cargo-LeRoux method fares substantially worse in our tests, although its simplicity may still offer some merits in certain situations. Overall, we find the well-balanced treatment of gravity in combination with low Mach number flux functions essential to reproducing correct physical solutions to challenging stellar slow-flow problems on affordable collocated grids. △ Less

Submitted 8 July, 2021; v1 submitted 25 February, 2021; originally announced February 2021.

Comments: Accepted for publication in A&A

Report number: LA-UR-21-21056

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

Double detonations of sub-Chandrasekhar mass CO white dwarfs: Can different core and He shell masses explain variations of Type Ia supernovae?

Authors: Sabrina Gronow, Christine E. Collins, Stuart A. Sim, Friedrich K. Roepke

Abstract: Sub-Chandrasekhar mass carbon-oxygen white dwarfs (CO WDs) with a surface helium (He) shell have been proposed as progenitors of Type Ia supernovae (SNe Ia). If true, the resulting thermonuclear explosions should be able to account for at least some of the range of SNe Ia observables. To study this, we conduct a parameter study based on 3D simulations of double detonations in CO WDs with a He shel… ▽ More Sub-Chandrasekhar mass carbon-oxygen white dwarfs (CO WDs) with a surface helium (He) shell have been proposed as progenitors of Type Ia supernovae (SNe Ia). If true, the resulting thermonuclear explosions should be able to account for at least some of the range of SNe Ia observables. To study this, we conduct a parameter study based on 3D simulations of double detonations in CO WDs with a He shell, assuming different core and shell masses. An admixture of C to the shell and solar metallicity are included in the models. The hydrodynamic simulations are carried out using the AREPO code. This allows us to follow the He shell detonation with high numerical resolution, and improves the reliability of predicted nucleosynthetic shell detonation yields. The addition of C to the shell leads to a lower production of 56Ni while including solar metallicity increases the production of IMEs. The production of higher mass elements is further shifted to stable isotopes at solar metallicity. Moreover, we find different core detonation ignition mechanisms depending on the core and shell mass configuration. This has an influence on the ejecta structure. We present the bolometric light curves predicted from our explosion simulations using the radiative transfer code ARTIS, and make comparisons with SNe Ia data. The bolometric light curves of our models show a range of brightnesses, able to account for sub-luminous to normal SNe Ia. We show the model bolometric width-luminosity relation compared to data for a range of viewing angles. We find that, on average, our brighter models lie within the observed data. The ejecta asymmetries produce a wide distribution of observables, which might account for outliers in the data. However, the models overestimate the extent of this compared to data. We also find the bolometric decline rate over 40 days appears systematically faster than data. (abridged) △ Less

Submitted 12 February, 2021; originally announced February 2021.

Comments: submitted to A&A, second revision after referee's comments

Journal ref: A&A 649, A155 (2021)

From supernova to supernova remnant: comparison of thermonuclear explosion models

Authors: Gilles Ferrand, Donald C. Warren, Masaomi Ono, Shigehiro Nagataki, Friedrich K. Roepke, Ivo R. Seitenzahl, Florian Lach, Hiroyoshi Iwasaki, Toshiki Sato

Abstract: Progress in the three-dimensional modeling of supernovae (SN) prompts us to revisit the supernova remnant (SNR) phase. We continue our study of the imprint of a thermonuclear explosion on the SNR it produces, that we started with a delayed-detonation model of a Chandrasekhar-mass white dwarf. Here we compare two different types of explosion models, each with two variants: two delayed detonation mo… ▽ More Progress in the three-dimensional modeling of supernovae (SN) prompts us to revisit the supernova remnant (SNR) phase. We continue our study of the imprint of a thermonuclear explosion on the SNR it produces, that we started with a delayed-detonation model of a Chandrasekhar-mass white dwarf. Here we compare two different types of explosion models, each with two variants: two delayed detonation models (N100ddt, N5ddt) and two pure deflagration models (N100def, N5def), where the N number parametrizes the ignition. The output of each SN simulation is used as input of a SNR simulation carried on until 500 yr after the explosion. While all SNR models become more spherical over time and overall display the theoretical structure expected for a young SNR, clear differences are visible amongst the models, depending on the geometry of the ignition and on the presence or not of detonation fronts. Compared to N100 models, N5 models have a strong dipole component, and produce asymmetric remnants. N5def produces a regular-looking, but offset remnant, while N5ddt produces a two-sided remnant. Pure deflagration models exhibit specific traits: a central over-density, because of the incomplete explosion, and a network of seam lines across the surface, boundaries between burning cells. Signatures from the SN dominate the morphology of the SNR up to 100 yr to 300 yr after the explosion, depending on the model, and are still measurable at 500 yr, which may provide a way of testing explosion models. △ Less

Submitted 9 November, 2020; originally announced November 2020.

Comments: Accepted by ApJ. Animated figures will be available on the online version

Report number: RIKEN-iTHEMS-Report-20

Nucleosynthesis imprints from different Type Ia Supernova explosion scenarios and implications for galactic chemical evolution

Authors: F. Lach, F. K. Roepke, I. R. Seitenzahl, B. Coté, S. Gronow, A. J. Ruiter

Abstract: We analyze the nucleosynthesis yields of various Type Ia supernova explosion simulations including pure detonations in sub- Chandrasekhar mass white dwarfs, double detonations and pure helium detonations of sub-Chandrasekhar mass white dwarfs with an accreted helium envelope, a violent merger model of two white dwarfs and deflagrations as well as delayed detonations in Chandrasekhar mass white dwa… ▽ More We analyze the nucleosynthesis yields of various Type Ia supernova explosion simulations including pure detonations in sub- Chandrasekhar mass white dwarfs, double detonations and pure helium detonations of sub-Chandrasekhar mass white dwarfs with an accreted helium envelope, a violent merger model of two white dwarfs and deflagrations as well as delayed detonations in Chandrasekhar mass white dwarfs. We focus on the iron peak elements Mn, Zn and Cu. To this end, we also briefly review the different burning regimes and production sites of these elements as well as the results of abundance measurements and several galactic chemical evolution studies. We find that super-solar values of [Mn/Fe] are not restricted to Chandrasekhar mass explosion models. Scenarios including a helium detonation can significantly contribute to the production of Mn, in particular the models proposed for calcium-rich transients. Although Type Ia supernovae are often not accounted for as production sites of Zn and Cu, our models involving helium shell detonations can produce these elements in super-solar ratios relative to Fe. Our results suggest a re-consideration of Type Ia supernova yields in galactic chemical evolution models. A detailed comparison with observations can provide new insight into the progenitor and explosion channels of these events. △ Less

Submitted 27 October, 2020; originally announced October 2020.

Comments: 15 paged, 8 figures, accepted for publication in A&A

Journal ref: A&A 644, A118 (2020)

Common-envelope evolution with an asymptotic giant branch star

Authors: C. Sand, S. T. Ohlmann, F. R. N. Schneider, R. Pakmor, F. K. Roepke

Abstract: Common-envelope phases are decisive for the evolution of many binary systems. Of particular interest are cases with asymptotic giant branch (AGB) primary stars, because they are thought to be progenitors of various astrophysical transients. In three-dimensional hydrodynamic simulations with the moving-mesh code AREPO, we study the common-envelope evolution of a $1.0\,M_{\odot}$ early-AGB star with… ▽ More Common-envelope phases are decisive for the evolution of many binary systems. Of particular interest are cases with asymptotic giant branch (AGB) primary stars, because they are thought to be progenitors of various astrophysical transients. In three-dimensional hydrodynamic simulations with the moving-mesh code AREPO, we study the common-envelope evolution of a $1.0\,M_{\odot}$ early-AGB star with companions of different masses. Although the stellar envelope of the AGB star is less tightly bound than that of a red giant, we find that the release of orbital energy of the core binary is insufficient to eject more than about twenty percent of the envelope mass. Ionization energy released in the expanding envelope, however, can lead to complete envelope ejection. Because recombination proceeds largely at high optical depths in our simulations, it is likely that this effect indeed plays a significant role in the considered systems. The efficiency of mass loss and the final orbital separation of the core binary system depend on the mass ratio between the companion and the primary star. Our results suggest a linear relation between the ratio of final to initial orbital separation and this parameter. △ Less

Submitted 2 November, 2020; v1 submitted 21 July, 2020; originally announced July 2020.

Comments: 12 pages, 9 figures, 5 tables; accepted for publication by A&A

Journal ref: A&A 644, A60 (2020)

Formation of sdB-stars via common envelope ejection by substellar companions

Authors: M. Kramer, F. R. N. Schneider, S. T. Ohlmann, S. Geier, V. Schaffenroth, R. Pakmor, F. K. Roepke

Abstract: Common envelope (CE) phases in binary systems where the primary star reaches the tip of the red giant branch are discussed as a formation scenario for hot subluminous B-type (sdB) stars. For some of these objects, observations point to very low-mass companions. In hydrodynamical CE simulations with the moving-mesh code AREPO, we test whether low-mass objects can successfully unbind the envelope. T… ▽ More Common envelope (CE) phases in binary systems where the primary star reaches the tip of the red giant branch are discussed as a formation scenario for hot subluminous B-type (sdB) stars. For some of these objects, observations point to very low-mass companions. In hydrodynamical CE simulations with the moving-mesh code AREPO, we test whether low-mass objects can successfully unbind the envelope. The success of envelope removal in our simulations critically depends on whether or not the ionization energy released by recombination processes in the expanding material is taken into account. If this energy is thermalized locally, envelope ejection eventually leading to the formation of an sdB star is possible with companion masses down to the brown dwarf range. For even lower companion masses approaching the regime of giant planets, however, envelope removal becomes increasingly difficult or impossible to achieve. Our results are consistent with current observational constraints on companion masses of sdB stars. Based on a semianalytic model, we suggest a new criterion for the lowest companion mass that is capable of triggering a dynamical response of the primary star thus potentially facilitating the ejection of a common envelope. This gives an estimate consistent with the findings of our hydrodynamical simulations. △ Less

Submitted 30 June, 2020; originally announced July 2020.

Comments: 12 pages, 8 figures, 3 tables; submitted to A&A

Journal ref: A&A 642, A97 (2020)

Photometric detection of internal gravity waves in upper main-sequence stars. II. Combined TESS photometry and high-resolution spectroscopy

Authors: D. M. Bowman, S. Burssens, S. Simón-Díaz, P. V. F. Edelmann, T. M. Rogers, L. Horst, F. K. Roepke, C. Aerts

Abstract: Context. Massive stars are predicted to excite internal gravity waves (IGWs) by turbulent core convection and from turbulent pressure fluctuations in their near-surface layers. These IGWs are extremely efficient at transporting angular momentum and chemical species within stellar interiors, but they remain largely unconstrained observationally. Aims. We aim to characterise the photometric detectio… ▽ More Context. Massive stars are predicted to excite internal gravity waves (IGWs) by turbulent core convection and from turbulent pressure fluctuations in their near-surface layers. These IGWs are extremely efficient at transporting angular momentum and chemical species within stellar interiors, but they remain largely unconstrained observationally. Aims. We aim to characterise the photometric detection of IGWs across a large number of O and early-B stars in the Hertzsprung-Russell diagram, and explain the ubiquitous detection of stochastic variability in the photospheres of massive stars. Methods. We combined high-precision time-series photometry from the NASA Transiting Exoplanet Survey Satellite with high-resolution ground-based spectroscopy of 70 stars with spectral types O and B to probe the relationship between the photometric signatures of IGWs and parameters such as spectroscopic mass, luminosity, and macroturbulence. Results. A relationship is found between the location of a star in the spectroscopic Hertzsprung-Russell diagram and the amplitudes and frequencies of stochastic photometric variability in the light curves of massive stars. Furthermore, the properties of the stochastic variability are statistically correlated with macroturbulent velocity broadening in the spectral lines of massive stars. Conclusions. The common ensemble morphology for the stochastic low-frequency variability detected in space photometry and its relationship to macroturbulence is strong evidence for IGWs in massive stars, since these types of waves are unique in providing the dominant tangential velocity field required to explain the observed spectroscopy. △ Less

Submitted 8 June, 2020; v1 submitted 4 June, 2020; originally announced June 2020.

Comments: Accepted for publication in A&A. 10 pages and 5 figures (and an additional 7 pages of appendix tables and figures). Resolution of appendix figures have been downgraded to meet arXiv's maximum file size of 15Mb. This version (2ver) is post A&A language editing

Journal ref: A&A 640, A36 (2020)

Fully compressible simulations of waves and core convection in main-sequence stars

Authors: L. Horst, P. V. F. Edelmann, R. Andrassy, F. K. Roepke, D. M. Bowman, C. Aerts, R. P. Ratnasingam

Abstract: Context. Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and rather high artificial viscosity and thermal diffusivity are needed for numerical stability. Direct comparison with obse… ▽ More Context. Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and rather high artificial viscosity and thermal diffusivity are needed for numerical stability. Direct comparison with observations is thus limited. Aims. We explore the capabilities of our compressible multidimensional hydrodynamics code SLH to simulate stellar oscillations. Methods. We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic $3\mathrm{M}_\odot$ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results. Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with observed spectra of massive stars. Conclusion. The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows the coupling of gravity and pressure waves to be studied too. △ Less

Submitted 14 September, 2020; v1 submitted 4 June, 2020; originally announced June 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 641, A18 (2020)

Long-term evolution of a magnetic massive merger product

Authors: F. R. N. Schneider, S. T. Ohlmann, Ph. Podsiadlowski, F. K. Roepke, S. A. Balbus, R. Pakmor

Abstract: About 10% of stars more massive than ${\approx}\,1.5\,\mathrm{M}_\odot$ have strong, large-scale surface magnetic fields and are being discussed as progenitors of highly-magnetic white dwarfs and magnetars. The origin of these fields remains uncertain. Recent 3D magnetohydrodynamical simulations have shown that strong magnetic fields can be generated in the merger of two massive stars. Here, we fo… ▽ More About 10% of stars more massive than ${\approx}\,1.5\,\mathrm{M}_\odot$ have strong, large-scale surface magnetic fields and are being discussed as progenitors of highly-magnetic white dwarfs and magnetars. The origin of these fields remains uncertain. Recent 3D magnetohydrodynamical simulations have shown that strong magnetic fields can be generated in the merger of two massive stars. Here, we follow the long-term evolution of such a 3D merger product in a 1D stellar evolution code. During a thermal relaxation phase after the coalescence, the merger product reaches critical surface rotation, sheds mass and then spins down primarily because of internal mass readjustments. The spin of the merger product after thermal relaxation is mainly set by the co-evolution of the star-torus structure left after coalescence. This evolution is still uncertain, so we also consider magnetic braking and other angular-momentum-gain and -loss mechanisms that may influence the final spin of the merged star. Because of core compression and mixing of carbon and nitrogen in the merger, enhanced nuclear burning drives a transient convective core that greatly contributes to the rejuvenation of the star. Once the merger product relaxed back to the main sequence, it continues its evolution similar to that of a genuine single star of comparable mass. It is a slow rotator that matches the magnetic blue straggler $τ$ Sco. Our results show that merging is a promising mechanism to explain some magnetic massive stars and it may also be key to understand the origin of the strong magnetic fields of highly-magnetic white dwarfs and magnetars. △ Less

Submitted 11 May, 2020; originally announced May 2020.

Comments: 17 pages (incl. appendix), 14 figures, 2 tables; accepted for publication in MNRAS

White dwarf deflagrations for Type Iax supernovae: Polarisation signatures from the explosion and companion interaction

Authors: M. Bulla, Z. -W. Liu, F. K. Roepke, S. A. Sim, M. Fink, M. Kromer, R. Pakmor, I. R. Seitenzahl

Abstract: Growing evidence suggests that Type Iax supernovae might be the result of thermonuclear deflagrations of Chandrasekhar-mass white dwarfs in binary systems. We carry out Monte Carlo radiative transfer simulations and predict spectropolarimetric features originating from the supernova explosion and subsequent ejecta interaction with the companion star. Specifically, we calculate viewing-angle depend… ▽ More Growing evidence suggests that Type Iax supernovae might be the result of thermonuclear deflagrations of Chandrasekhar-mass white dwarfs in binary systems. We carry out Monte Carlo radiative transfer simulations and predict spectropolarimetric features originating from the supernova explosion and subsequent ejecta interaction with the companion star. Specifically, we calculate viewing-angle dependent flux and polarisation spectra for a 3D model simulating the deflagration of a Chandrasekhar-mass white dwarf and, for a second model, simulating the ejecta interaction with a main-sequence star. We find that the intrinsic signal is weakly polarised and only mildly viewing-angle dependent, owing to the overall spherical symmetry of the explosion and the depolarising contribution of iron-group elements dominating the ejecta composition. The interaction with the companion star carves out a cavity in the ejecta and produces a detectable, but modest signal that is significant only at relatively blue wavelengths ($\lesssim$ 5000 $\unicode{x212B}$). In particular, increasingly fainter and redder spectra are predicted for observer orientations further from the cavity, while a modest polarisation signal $P\sim0.2$ per cent is found at blue wavelengths for orientations 30$^\circ$ and 45$^\circ$ away from the cavity. We find a reasonable agreement between the interaction model viewed from these orientations and spectropolarimetric data of SN 2005hk and interpret the maximum-light polarisation signal seen at blue wavelengths for this event as a possible signature of the ejecta-companion interaction. We encourage further polarimetric observations of SNe Iax to test whether our results can be extended and generalised to the whole SN Iax class. △ Less

Submitted 27 February, 2020; v1 submitted 25 February, 2020; originally announced February 2020.

Comments: 7 pages, 3 figures; accepted for publication in A&A; metadata update

Report number: NORDITA 2020-019

Journal ref: A&A 635, A179 (2020)

SNe Ia from double detonations: Impact of core-shell mixing on the carbon ignition mechanism

Authors: Sabrina Gronow, Christine Collins, Sebastian T. Ohlmann, Ruediger Pakmor, Markus Kromer, Ivo R. Seitenzahl, Stuart A. Sim, Friedrich K. Roepke

Abstract: Sub-Chandrasekhar mass white dwarfs accreting a helium shell on a carbon-oxygen core are potential progenitors of normal Type Ia supernovae. This work focuses on the details of the onset of the carbon detonation in the double detonation sub-Chandrasekhar model. In order to simulate the influence of core-shell mixing on the carbon ignition mechanism, the helium shell and its detonation are followed… ▽ More Sub-Chandrasekhar mass white dwarfs accreting a helium shell on a carbon-oxygen core are potential progenitors of normal Type Ia supernovae. This work focuses on the details of the onset of the carbon detonation in the double detonation sub-Chandrasekhar model. In order to simulate the influence of core-shell mixing on the carbon ignition mechanism, the helium shell and its detonation are followed with an increased resolution compared to the rest of the star treating the propagation of the detonation wave more accurately. This significantly improves the predictions of the nucleosynthetic yields from the helium burning. The simulations were carried out with the AREPO code. A carbon-oxygen core with a helium shell was set up in one dimension and mapped to three dimensions. We ensured the stability of the white dwarf with a relaxation step before the hydrodynamic detonation simulation started. Synthetic observables were calculated with the radiative transfer code ARTIS. An ignition mechanism of the carbon detonation was observed, which received little attention before. In this "scissors mechanism", the impact the helium detonation wave has on unburnt material when converging opposite to its ignition spot is strong enough to ignite a carbon detonation. This is possible in a carbon enriched transition region between the core and shell. The detonation mechanism is found to be sensitive to details of the core-shell transition and our models illustrate the need to consider core-shell mixing taking place during the accretion process. Even though the detonation ignition mechanism differs form the converging shock mechanism, the differences in the synthetic observables are not significant. Though they do not fit observations better than previous simulations, they illustrate the need for multi-dimensional simulations. △ Less

Submitted 3 February, 2020; originally announced February 2020.

Comments: Accepted for publication in A&A

Journal ref: A&A 635, A169 (2020)

Monte Carlo radiative transfer for the nebular phase of Type Ia supernovae

Authors: L. J. Shingles, S. A. Sim, M. Kromer, K. Maguire, M. Bulla, C. Collins, C. P. Ballance, A. S. Michel, C. A. Ramsbottom, F. K. Roepke, I. R. Seitenzahl, N. B. Tyndall

Abstract: We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium (non-LTE) population and ionisation solver, a new multi-frequency radiation field model, and a new atomic dataset with forbidden transitions. We treat collisions with non-the… ▽ More We extend the range of validity of the ARTIS 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium (non-LTE) population and ionisation solver, a new multi-frequency radiation field model, and a new atomic dataset with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionisation, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf in which the possible effects of gravitational settling of Ne22 prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M$_\odot$ white dwarf model obtained when 5.5 Gyr of very-efficient settling is assumed to a similar model without settling. We find that this degree of Ne22 settling has only a modest effect on the resulting nebular spectra due to increased Ni58 abundance. Due to the high ionisation in sub-Chandrasekhar models, the nebular [Ni II] emission remains negligible, while the [Ni III] line strengths are increased and the overall ionisation balance is slightly lowered in the model with Ne22 settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe III] emission relative to [Fe II] in comparison to observations of normal Type Ia supernovae. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Comments: 15 pages, 12 figures. Accepted for publication by MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 492, Issue 2, February 2020, Pages 2029-2043

How mergers magnetise massive stars

Authors: Fabian R. N. Schneider, Sebastian T. Ohlmann, Philipp Podsiadlowski, Friedrich K. Röpke, Steven A. Balbus, Rüdiger Pakmor, Volker Springel

Abstract: Magnetic fields are ubiquitous in the Universe. The Sun's magnetic field drives the solar wind and causes solar flares and other energetic surface phenomena that profoundly affect space weather here on Earth. The first magnetic field in a star other than the Sun was detected in 1947 in the peculiar A-type star 78 Vir. It is now known that the magnetic fields of the Sun and other low-mass stars (<1… ▽ More Magnetic fields are ubiquitous in the Universe. The Sun's magnetic field drives the solar wind and causes solar flares and other energetic surface phenomena that profoundly affect space weather here on Earth. The first magnetic field in a star other than the Sun was detected in 1947 in the peculiar A-type star 78 Vir. It is now known that the magnetic fields of the Sun and other low-mass stars (<1.5 solar masses) are generated in-situ by a dynamo process in their turbulent, convective envelopes. Unlike such stars, intermediate-mass and high-mass stars (>1.5 solar masses; referred to as "massive" stars here) have relatively quiet, radiative envelopes where a solar-like dynamo cannot operate. However, about 10% of them, including 78 Vir, have strong, large-scale surface magnetic fields whose origin has remained a major mystery. The massive star $τ$ Sco is a prominent member of this group and appears to be surprisingly young compared to other presumably coeval members of the Upper Scorpius association. Here, we present the first 3D magneto-hydrodynamical simulations of the coalescence of two massive main-sequence stars and 1D stellar evolution computations of the subsequent evolution of the merger product that can explain $τ$ Sco's magnetic field, apparent youth and other observed characteristics. We argue that field amplification in stellar mergers is a general mechanism to form strongly-magnetised massive stars. These stars are promising progenitors of those neutron stars that host the strongest magnetic fields in the Universe, so-called magnetars, and that may give rise to some of the enigmatic fast radio bursts. Strong magnetic fields affect the explosions of core-collapse supernovae and, moreover, those magnetic stars that have rapidly-rotating cores at the end of their lives might provide the right conditions to power long-duration gamma-ray bursts and super-luminous supernovae. △ Less

Submitted 30 October, 2019; originally announced October 2019.

Comments: This is the authors' initially submitted version of a paper that is published in the October 10th 2019 issue of Nature at https://www.nature.com/articles/s41586-019-1621-5

Journal ref: Nature 574, 211-214 (2019)

A New Model for Electron-Capture Supernovae in Galactic Chemical Evolution

Authors: S. Jones, B. Côté, F. K. Roepke, S. Wanajo

Abstract: We examine the contribution of electron-capture supernovae (ECSNe), low-mass SNe from collapsing Fe cores (FeCCSNe), and rotating massive stars to the chemical composition of the Galaxy. Our model includes contributions to chemical evolution from both thermonuclear ECSNe (tECSNe) and gravitational collapse ECSNe (cECSNe). We show that if ECSNe are predominantly gravitational collapse SNe but about… ▽ More We examine the contribution of electron-capture supernovae (ECSNe), low-mass SNe from collapsing Fe cores (FeCCSNe), and rotating massive stars to the chemical composition of the Galaxy. Our model includes contributions to chemical evolution from both thermonuclear ECSNe (tECSNe) and gravitational collapse ECSNe (cECSNe). We show that if ECSNe are predominantly gravitational collapse SNe but about 15% are partial thermonuclear explosions, the model is able to reproduce the solar abundances of several important and problematic isotopes including $^{48}$Ca, $^{50}$Ti and $^{54}$Cr together with $^{58}$Fe, $^{64}$Ni, $^{82}$Se and $^{86}$Kr and several of the Zn--Zr isotopes. A model in which no cECSNe occur, only tECSNe with low-mass FeCCSNe or rotating massive stars, proves also very successful at reproducing the solar abundances for these isotopes. Despite the small mass range for the progenitors of ECSNe and low-mass FeCCSNe, the large production factors suffice for the solar inventory of the above isotopes. Our model is compelling because it introduces no new tensions with the solar abundance distribution for a Milky Way model -- only tending to improve the model predictions for several isotopes. The proposed astrophysical production model thus provides a natural and elegant way to explain one of the last uncharted territories on the periodic table of astrophysical element production. △ Less

Submitted 6 August, 2019; originally announced August 2019.

Comments: 10 pages, 4 figures; accepted for publication in the Astrophysical Journal

Discovery of an Exceptionally Strong $β$-Decay Transition of $^{20}$F and Implications for the Fate of Intermediate-Mass Stars

Authors: O. S. Kirsebom, S. Jones, D. F. Strömberg, G. Martínez-Pinedo, K. Langanke, F. K. Roepke, B. A. Brown, T. Eronen, H. O. U. Fynbo, M. Hukkanen, A. Idini, A. Jokinen, A. Kankainen, J. Kostensalo, I. Moore, H. Möller, S. T. Ohlmann, H. Penttilä, K. Riisager, S. Rinta-Antila, P. C. Srivastava, J. Suhonen, W. H. Trzaska, J. Äystö

Abstract: A significant fraction of stars between 7-11 solar masses are thought to become supernovae, but the explosion mechanism is unclear. The answer depends critically on the rate of electron capture on $^{20}$Ne in the degenerate oxygen-neon stellar core. However, due to the unknown strength of the transition between the ground states of $^{20}$Ne and $^{20}$F, it has not previously been possible to fu… ▽ More A significant fraction of stars between 7-11 solar masses are thought to become supernovae, but the explosion mechanism is unclear. The answer depends critically on the rate of electron capture on $^{20}$Ne in the degenerate oxygen-neon stellar core. However, due to the unknown strength of the transition between the ground states of $^{20}$Ne and $^{20}$F, it has not previously been possible to fully constrain the rate. By measuring the transition, we have established that its strength is exceptionally large and enhances the capture rate by several orders of magnitude. This has a decisive impact on the evolution of the core, increasing the likelihood that the star is (partially) disrupted by a thermonuclear explosion rather than collapsing to form a neutron star. Importantly, our measurement resolves the last remaining nuclear physics uncertainty in the final evolution of degenerate oxygen-neon stellar cores, allowing future studies to address the critical role of convection, which at present is poorly understood. △ Less

Submitted 2 November, 2019; v1 submitted 22 May, 2019; originally announced May 2019.

Comments: 6 pages, 5 figures, supplemental material; submitted to Science on January 26, 2019; revised version submitted to Physical Review Letters on May 22, 2019; re-submitted to Physical Review Letters on August 15, 2019, after minor revisions; revised title

Journal ref: Phys. Rev. Lett. 123, 262701 (2019)