Search

The double detonation is a widely discussed explosion mechanism for Type Ia supernovae, whereby a helium shell detonation ignites a secondary detonation in the carbon/oxygen core of a white dwarf. Even for modern models that invoke relatively small He shell masses, many previous studies have found that the products of the helium shell detonation lead to discrepancies with normal Type Ia supernovae, such as strong Ti II absorption features, extremely red light curves and too large a variation with viewing direction. It has been suggested that non local thermodynamic equilibrium (non-LTE) effects may help to reduce these discrepancies with observations. Here we carry out full non-LTE radiative transfer simulations for a recent double detonation model with a relatively small helium shell mass of 0.05 M$_\odot$. We construct 1D models representative of directions in a 3D explosion model to give an indication of viewing angle dependence. The full non-LTE treatment leads to improved agreement between the models and observations. The light curves become less red, due to reduced absorption by the helium shell detonation products, since these species are more highly ionised. Additionally, the expected variation with observer direction is reduced. The full non-LTE treatment shows promising improvements, and reduces the discrepancies between the double detonation models and observations of normal Type Ia supernovae., Comment: Submitted to MNRAS

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 analysed an equal-mass merger with two $0.3\,M_\odot$ white dwarfs, and an unequal-mass merger with white dwarfs of $0.25\,M_\odot$ and $0.35\,M_\odot$. We simulated the inspiral, merger, and further evolution of the merger remnant for about $50$ rotations. We found efficient magnetic field amplification in both mergers via a small-scale dynamo, reproducing previous results of stellar merger simulations. The magnetic field saturates at a similar strength for both simulations. We then identified 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 via a large-scale dynamo driven by the magneto-rotational instability. Finally, we speculate that in the unequal-mass merger remnant, helium burning will initially start in a shell around a cold core, rather than in the centre. This forms a convection zone that coincides with the region that contains most of the magnetic energy, and likely destroys the strong, ordered field. Ohmic resistivity might then quickly erase the remaining small-scale 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., Comment: 13 pages, 10 figures, accepted by A&A, comments welcome

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)., Comment: 20 pages, 12 figures, submitted to Phys. Rev. D

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, $\Theta=R_{A}/R_{\rm sink}$, and is much less than the dynamical friction for sufficiently small accretors ($R_{\rm sink}\ll R_{A}$)., Comment: 14 pages, 7 figures. Submitted to ApJ. Comments welcome!

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., Comment: 12 pp., accepted for publication in A & A

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-\.Zytkow 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., Comment: 21 pages, 19 figures, 1 table. Accepted for publication on A&A

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 {\lambda} 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 {\lambda} 10927 feature but over time separates to form a secondary feature to the blue wing of the Mg II {\lambda} 10927 feature. We compare our simulation to observations of iPTF13ebh, which showed a similar feature to the blue wing of the Mg II {\lambda} 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 {\lambda} 10830. This suggests that He I {\lambda} 10830 could be a signature of the double detonation scenario., Comment: 7 pages, accepted by MNRAS

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., Comment: An invited review (accepted for publication in RAA)

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-\.Zytkow objects., Comment: Accepted for publication in A&A, 12 pages, 8 figures

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., Comment: 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)

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., Comment: Spectra available at https://doi.org/10.5281/zenodo.7997388

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., Comment: 24 pages, 18 figures, accepted for publication in MNRAS

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., Comment: 13 pages, 7 figures, accepted for publication in ApJ

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., Comment: 13 pages. Accepted by MNRAS

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.

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., Comment: 14 pages, 8 figures, 1 table, accepted for publication in ApJ

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.), Comment: 16 pages, 11 figures, accepted by A&A

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., Comment: 11 pages, 9 figures, 1 table, accepted for publication in A&A

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 $\alpha$-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., Comment: accepted for publication by A&A

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), Comment: submitted to A&A, second revision after referee's comments

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., Comment: Accepted by ApJ. Animated figures will be available on the online version

Instrumental variable methods are popular choices in combating unmeasured confounding to obtain less biased effect estimates. However, we demonstrate that alternative methods may give less biased estimates depending on the nature of unmeasured confounding. Treatment preferences of clusters (e.g., physician practices) are the most f6requently used instruments in instrumental variable analyses (IVA). These preference-based IVAs are usually conducted on data clustered by region, hospital/facility, or physician, where unmeasured confounding often occurs within or between clusters. We aim to quantify the impact of unmeasured confounding on the bias of effect estimators in IVA, as well as alternative methods including ordinary least squares regression, linear mixed models (LMM) and fixed effect models (FE) to study the effect of a continuous exposure (e.g., treatment dose). We derive bias formulae of estimators from these four methods in the presence of unmeasured within- and/or between-cluster confounders. We show that IVAs can provide consistent estimates when unmeasured within-cluster confounding exists, but not when between-cluster confounding exists. On the other hand, FEs and LMMs can provide consistent estimates when unmeasured between-cluster confounding exits, but not for within-cluster confounding. Whether IVAs are advantageous in reducing bias over FEs and LMMs depends on the extent of unmeasured within-cluster confounding relative to between-cluster confounding. Furthermore, the impact of unmeasured between-cluster confounding on IVA estimates is larger than the impact of unmeasured within-cluster confounding on FE and LMM estimates. We illustrate these methods through data applications. Our findings provide guidance for choosing appropriate methods to combat the dominant types of unmeasured confounders and help interpret statistical results in the context of unmeasured confounding

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., Comment: Accepted for publication in A&A

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 $\tau$ 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 $\tau$ 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., Comment: 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

Recent progress in the three-dimensional modeling of supernovae (SN) has shown the importance of asymmetries for the explosion. This calls for a reconsideration of the modeling of the subsequent phase, the supernova remnant (SNR), which has commonly relied on simplified ejecta models. In this paper we bridge SN and SNR studies by using the output of a SN simulation as the input of a SNR simulation carried on until 500~yr. We consider the case of a thermonuclear explosion of a carbon-oxygen white dwarf star as a model for a Type Ia SN; specifically we use the N100 delayed detonation model of Seitenzahl et al 2013. In order to analyze the morphology of the SNR, we locate the three discontinuities that delineate the shell of shocked matter: the forward shock, the contact discontinuity, and the reverse shock, and we decompose their radial variations as a function of angular scale and time. Assuming a uniform ambient medium, we find that the impact of the SN on the SNR may still be visible after hundreds of years. Previous 3D simulations aiming at reproducing Tycho's SNR, that started out from spherically symmetric initial conditions, failed to reproduce structures at the largest angular scales observed in X-rays. Our new simulations strongly suggest that the missing ingredient was the initial asymmetries from the SN itself. With this work we establish a way of assessing the viability of SN models based on the resulting morphology of the SNR., Comment: Accepted by ApJ. Animated figures will be available on the online version

We develop a second order well-balanced finite volume scheme for compressible Euler equations with a gravitational source term. The well-balanced property holds for arbitrary hydrostatic solutions of the corresponding Euler equations without any restriction on the equation of state. The hydrostatic solution must be known a priori either as an analytical formula or as a discrete solution at the grid points. The scheme can be applied on curvilinear meshes and in combination with any consistent numerical flux function and time stepping routines. These properties are demonstrated on a range of numerical tests., Comment: 32 pages